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
em séries de Taylor e truncar os termos até a segunda ordem:
![{\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c2bb717a3b5b4a2b6a6b2c7f05d9ba513cb3b90d)
Relacionando as derivadas do tempo e do espaço:
![{\displaystyle {\frac {\partial u}{\partial t}}=-c{\frac {\partial u}{\partial x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10cc27c689797ede7cc676db5e9c859d630c5980)
![{\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}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/446e3b444c7cd69b0ac73e63a67388ec2263a65c)
Podemos fazer substituições na equação (1), obtendo:
![{\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be29fa3953bce1057db9449ab202d7b433eb4763)
Usando diferenças centradas de primeira e segunda ordem em relação ao espaço:
![{\displaystyle {\frac {\partial u}{\partial x}}={\frac {u_{i+1}^{n}-u_{i-1}^{n}}{2\Delta x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf3ba530e001a60e19cb22be6fdf9242625dfc62)
![{\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}={\frac {u_{i+1}^{n}-2u_{i}^{n}+u_{i-1}^{n}}{\Delta x^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd1e1500d6d68aefabb402e7bfe2c59e640abe18)
E substituindo em (2), obtemos assim o método de Lax-Wendroff:
![{\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc6fc264d70cb0067daba6aaec919409641b9c05)
O qual também pode ser mostrado em relação ao número de Courant–Friedrichs–Lewy(CFL):
![{\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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e986c89d9f6597f55de9cad0fe9573417f3c6d5c)
Referências
- ↑ 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
Resolução de equações diferenciais parciais |
---|
Método das diferenças finitas | |
---|
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](//upload.wikimedia.org/wikipedia/commons/thumb/3/35/E-to-the-i-pi.svg/34px-E-to-the-i-pi.svg.png) | Este artigo sobre matemática é um esboço. Você pode ajudar a Wikipédia expandindo-o. |