O que a análise de estabilidade de Von Neumann nos diz sobre equações de diferenças finitas não lineares?

9

Estou lendo um papel [1] onde eles resolver o seguinte não-linear equação usando métodos de diferenças finitas. Eles também analisam a estabilidade dos esquemas usando a análise de estabilidade de Von Neumann. No entanto, como os autores percebem, isso é aplicável apenas aos PDE lineares. Portanto, os autores resolvem isso "congelando" o termo não linear, ou seja, eles substituem o termo por , onde é "considerado como representando localmente valores constantes de

{você}_{t} + {você}_{x} + você {você}_{x} - {você}_{x x t} = 0 0

$\begin{equation} u_t + u_x + uu_x - u_{xxt} = 0 \end{equation}$

u u_{x}

$uu_x$

U u_{x}

$Uu_x$

U

$U$

"

u

$u$

Então, minha pergunta é dupla:

1: como interpretar esse método e por que (não) funciona?

2: poderíamos também substituir o termo termo , onde é "considerado como representando localmente valores constantes de "? $uu_x$ $uU_x$ $U_x$ $u_x$

Referências

Eilbeck, JC e GR McGuire. "Estudo numérico da equação de onda longa regularizada I: métodos numéricos". Journal of Computational Physics 19.1 (1975): 43-57.

— Caçador
fonte

11

Você digitou errado a equação. A equação no trabalho é a equação RLW.

— Ômer 04/03

3

Perguntas relacionadas, sem respostas completas: scicomp.stackexchange.com/q/8717/713 , mathoverflow.net/q/186760 , scicomp.stackexchange.com/q/16142 , scicomp.stackexchange.com/q/6863 . Eu acho que, heuristicamente falando, deve funcionar porque você está interessado na estabilidade dos modos de frequência muito alta (em que os erros ocorrem, comprimento de onda na ordem do espaçamento da malha), enquanto a solução em si varia com frequência muito menor, portanto, não há problema em congelar coeficientes e estudar a estabilidade dos PDE de coeficientes congelados.

— Kirill

2

Dei respostas a algumas das perguntas vinculadas por Kirill. Infelizmente, não estou ciente de nenhum resultado para a equação RLW, mas provavelmente a estabilidade pode ser comprovada desde que a solução seja suave o suficiente.

— David Ketcheson

1

O que você está dizendo é chamado de linearização. É uma técnica comum usada na análise de EDPs não lineares. O que é feito é lançar equações no formato,

$u_t+Au=0$

Aqui A é uma matriz resultante da linearização da equação.

Agora, para suas perguntas,

Como você está pensando, funciona até certo ponto, mas não até certo ponto. A utilidade é que a estabilidade pode ser comprovada para sistemas lineares, mas não prontamente para sistemas não lineares. Portanto, os resultados lineares são estendidos aos sistemas não lineares. Freqüentemente, diferentes métodos são adotados para casos particulares. Por exemplo,

$uu_x = \frac{1}{2}(u^2)_x$

qual é a forma de conservação. Assim,

$u_t+\frac{1}{2}(u^2)_x = 0$

quando representado em um sentido de volume finito, limita a evolução de u.

Qual é a utilidade de fazer a substituição. Você removerá a equação de uma forma de equação de onda. O que significaria que as soluções não se comportariam como uma equação de onda. Portanto, na análise de estabilidade, as soluções de teste também teriam que ser completamente diferentes e pouco físicas.

— Vikram
fonte

2

Para elaborar o argumento de linearização, em uu_x você deseja assumir que u é localmente constante, e não u_x, por duas razões: a) u varia mais lentamente que sua derivada eb) nesse caso em particular, se você assumir que u_x é localmente constante , por definição, você também presume que u é localmente linear, o que significa que derivadas espaciais mais altas são zero e isso não apenas introduz um erro de aproximação adicional, como também pode implicar que você esteja jogando fora o bebê com a água do banho, dependendo da sua equação.

— Domingo Tavella
fonte