Para equações de difusão (e muitas outras), pode-se provar que o método de Crank–Nicolson é incondicionalmenteestável.[3] Contudo, as soluções aproximadas podem ainda conter oscilações significativas caso a razão entre o passo de tempo e o quadrado do passo de espaço for grande (usualmente maior que 1/2). Por essa razão, sempre que grandes passos de tempo forem tomados, o método menos preciso de euler implícito é frequentemente utilizado, o qual é estável e imune à oscilações.
O método
The Crank–Nicolson stencil for a 1D problem.
O método de Crank-Nicolson é baseado em diferenças centradas no espaço, e na regra trapezoidal no tempo, é de segunda no tempo e no espaço. Por exemplo, para um caso unidimensional, se a equação diferencial parcial for
em seguida, fazendo , a equação para o método de Crank–Nicolson é a combinação do método de euler explícito em e do método de euler implícito em n+1 (deve-se notar, contudo, que o método por si só não é simplesmente a média desses dois métodos, já que a equação tem uma dependência implícita na solução):
Note que este é um método implícito: para conseguir o valor posterior de u no tempo, um sistema de equações algébricas deve ser resolvido. Se a equação diferencial parcial for não-linear, a discretização também deverá ser não-linear para que o avanço no tempo envolva a solução do sistema algébrico não-linear de equações, embora que linearizações sejam possíveis.
Exemplo: Difusão unidimensional
O método de Crank-Nicolson é frequentemente aplicado a problemas de difusão. Como exemplo, para a difusão linear:
a discretização de Crank-Nicolson é dada por:
ou reescrevendo, fazendo :
Referências
↑Tuncer Cebeci (2002). Convective Heat Transfer. [S.l.]: Springer. ISBN 0966846141
↑Crank, J.; Nicolson, P. (1947). «A practical method for numerical evaluation of solutions of partial differential equations of the heat conduction type». Proc. Camb. Phil. Soc. 43 (1): 50–67. doi:10.1007/BF02127704.|given1= e |primeiro1= redundantes (ajuda); |given2= e |primeiro2= redundantes (ajuda).
↑Thomas, J. W. (1995). Numerical Partial Differential Equations: Finite Difference Methods. Col: Texts in Applied Mathematics. 22. Berlin, New York: Springer-Verlag. ISBN 978-0-387-97999-1.. O exemplo 3.3.2 mostra que o método é incondicionalmente estável quando aplicado à .
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
Este artigo sobre matemática é um esboço. Você pode ajudar a Wikipédia expandindo-o.