Método de Lax–Wendroff

O método de Lax–Wendroff, em homenagem a Peter Lax e Burton Wendroff, é um método numérico para a resolução de equações hiperbólicas em derivadas parciais, baseado em diferenças finitas. É um método de segunda ordem no tempo e no espaço. Lax e Wendroff[1] apresentaram um método de discretização de segunda ordem para a solução de equações hiperbólicas, o que substituiu o método de Lax-friedrichs.

Ilustração do método

Para determinar o método de Lax-Wendroff, podemos expandir a variável u {\displaystyle u} em séries de Taylor e truncar os termos até a segunda ordem:

( 1 ) u i n + 1 = u i n + Δ t ( u t ) + Δ t 2 2 ( 2 u t 2 ) {\displaystyle (1)\qquad u_{i}^{n+1}=u_{i}^{n}+\Delta t\left({\frac {\partial u}{\partial t}}\right)+{\frac {\Delta t^{2}}{2}}\left({\frac {\partial ^{2}u}{\partial t^{2}}}\right)}

Relacionando as derivadas do tempo e do espaço:

u t = c u x {\displaystyle {\frac {\partial u}{\partial t}}=-c{\frac {\partial u}{\partial x}}}
2 u t 2 = t ( u t ) = c x ( c u x ) = c 2 2 u x 2 {\displaystyle {\frac {\partial ^{2}u}{\partial t^{2}}}={\frac {\partial }{\partial t}}\left({\frac {\partial u}{\partial t}}\right)=-c{\frac {\partial }{\partial x}}\left(-c{\frac {\partial u}{\partial x}}\right)=c^{2}{\frac {\partial ^{2}u}{\partial x^{2}}}}

Podemos fazer substituições na equação (1), obtendo:

( 2 ) u i n + 1 = u i n c Δ t ( u x ) + c 2 Δ t 2 2 ( 2 u x 2 ) {\displaystyle (2)\qquad u_{i}^{n+1}=u_{i}^{n}-c\Delta t\left({\frac {\partial u}{\partial x}}\right)+c^{2}{\frac {\Delta t^{2}}{2}}\left({\frac {\partial ^{2}u}{\partial x^{2}}}\right)}

Usando diferenças centradas de primeira e segunda ordem em relação ao espaço:

u x = u i + 1 n u i 1 n 2 Δ x {\displaystyle {\frac {\partial u}{\partial x}}={\frac {u_{i+1}^{n}-u_{i-1}^{n}}{2\Delta x}}}
2 u x 2 = u i + 1 n 2 u i n + u i 1 n Δ x 2 {\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}={\frac {u_{i+1}^{n}-2u_{i}^{n}+u_{i-1}^{n}}{\Delta x^{2}}}}

E substituindo em (2), obtemos assim o método de Lax-Wendroff:

u i n + 1 = u i n c Δ t 2 Δ x ( u i + 1 n u i 1 n ) + c 2 Δ t 2 2 Δ x 2 ( u i + 1 n 2 u i n + u i 1 n ) {\displaystyle u_{i}^{n+1}=u_{i}^{n}-{\frac {c\Delta t}{2\Delta x}}\left(u_{i+1}^{n}-u_{i-1}^{n}\right)+{\frac {c^{2}\Delta t^{2}}{2\Delta x^{2}}}\left(u_{i+1}^{n}-2u_{i}^{n}+u_{i-1}^{n}\right)}

O qual também pode ser mostrado em relação ao número de Courant–Friedrichs–Lewy(CFL):

u i n + 1 = u i n σ 2 ( u i + 1 n u i 1 n ) + σ 2 2 ( u i + 1 n 2 u i n + u i 1 n ) {\displaystyle u_{i}^{n+1}=u_{i}^{n}-{\frac {\sigma }{2}}\left(u_{i+1}^{n}-u_{i-1}^{n}\right)+{\frac {\sigma ^{2}}{2}}\left(u_{i+1}^{n}-2u_{i}^{n}+u_{i-1}^{n}\right)}

Referências

  1. Peter D. Lax, Burton Wendroff (1960). Systems of conservation laws, Communications in Pureand Applied Mathematics. [S.l.: s.n.] p. 217-237 
  • P.D Lax; B. Wendroff (1960). «Systems of conservation laws». Commun. Pure Appl Math. 13 (2): 217–237. doi:10.1002/cpa.3160130205  A referência emprega parâmetros obsoletos |coautor= (ajuda)
  • Michael J. Thompson, An Introduction to Astrophysical Fluid Dynamics, Imperial College Press, London, 2006.
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). «Section 20.1. Flux Conservative Initial Value Problems». Numerical Recipes: The Art of Scientific Computing 3rd ed. New York: Cambridge University Press. p. 1040. ISBN 978-0-521-88068-8 
  • v
  • d
  • e
Resolução de equações diferenciais parciais
Método das diferenças finitas
Relativos à equação do calor: esquema FTCS · Método de Crank–Nicolson Hiperbólicos: Método de Lax–Friedrichs · Método de Lax–Wendroff · Método de MacCormack · Esquema Upwind · Outros: Método da direção implícita alternada · Método FDTD
Método dos volumes finitos
Esquema de Godunov · Esquema de Alta-resolução · Esquema MUSCL · AUSM · Riemann solver
Método dos elementos finitos
hp-FEM · Método dos elementos finitos estendido · Método de Galerkin Descontínuo · Método dos elementos espectrais · Métodos Meshfree · Métodos Mortar
Outros métodos
Método espectral · Método pseudo espectral · Método das linhas · Método Multigrid · Método da colocação · Método Level set · Método dos elementos de contorno · Método de Fronteira Imersa · Método de elementos analíticos · Método Particle-in-cell · Análise Isogeométrica
Métodos de decomposição de domínios
Método de Schur · Método dos domínios fictícios · Método alternante de Schwarz · Método aditivo de Schwarz · Método aditivo abstrato de Schwarz · Método de Neumann–Dirichlet · Métodos Neumann–Neumann · Operador Poincaré–Steklov · Balancing domain decomposition · BDDC · FETI · FETI-DP
Ícone de esboço Este artigo sobre matemática é um esboço. Você pode ajudar a Wikipédia expandindo-o.
  • v
  • d
  • e