Méthode de Crank-Nicolson

Résumé

En mathématiques, en analyse numérique, la méthode de Crank-Nicolson est un algorithme simple permettant de résoudre des systèmes d'équations aux dérivées partielles. Cette méthode utilise les différences finies pour approcher une solution du problème : elle est numériquement stable et quadratique pour le temps. On peut facilement la généraliser à des problèmes à deux ou trois dimensions. Cette méthode, publiée en 1947, est le résultat des travaux de la mathématicienne britannique Phyllis Nicolson (1917 — 1968) et du physicien John Crank (1916 — 2006). Ils l'utilisèrent dans la résolution de l'équation de la chaleur. Son efficacité et sa simplicité en font un outil courant dans les simulations numériques, pour résoudre des problèmes de mécanique quantique, de thermodynamique hors-équilibre, de mécanique des fluides et d'électromagnétisme. Par ailleurs, un certain nombre de phénomènes pouvant être ramenés à l'étude de l'équation de la chaleur, son champ d'application est relativement

Source officielle

https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Crank-Nicolson

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Publications associées (10)

Publications associées

Chargement

Personnes associées (2)

Personnes associées

Chargement

Unités associées (1)

Unités associées

Chargement

Concepts associés (9)

Concepts associés

Chargement

Cours associés (16)

Cours associés

Chargement

Séances de cours associées (38)

Afficher plus

Séances de cours associées

Chargement

Space-time adaptive algorithms for parabolic problems

Virabouth Prachittham

We developed a space and time adaptive method to simulate electroosmosis and mass transport of a sample concentration within a network of microchannels. The space adaptive criteria is based on an error estimator derived using anisotropic interpolation estimates and a post-processing procedure. In order to improve the accuracy of the numerical solution and to reduce even further the computational cost of the numerical simulation, a time adaptive procedure is combined with the one in space. To do so, a time error estimator is derived for a first model problem, the linear heat equation discretized in time with the Crank-Nicolson scheme. The main difficulty is then to obtain an optimal second order error estimator. Applying standard energy techniques with a continuous, piecewise linear approximation in time fail in recovering the optimal order. To restore the appropriate rate of convergence, a continuous piecewise quadratic polynomial function in time is needed. For this purpose, two different quadratic functions are introduced and two different time error estimators are then derived. It turns out that the second error estimator is more efficient than the first one when considering our adaptive algorithm. Thus, using the second quadratic polynomial, an upper bound for the error is derived for a second model problem, the time-dependent convection-diffusion problem discretized in time with the Crank-Nicolson scheme. The corresponding space and time error estimators are finally used for the numerical simulation of mass transport of a sample concentration within a complex network of microchannels driven by an electroosmotic flow and/or by a pressure-driven flow. Numerical results presented show the efficiency and the robustness of this approach.

EPFL2009

Chargement

An a posteriori upper bound is derived for the nonstationary convection-diffusion problem using the Crank-Nicolson scheme and continuous, piecewise linear stabilized finite elements with large aspect ratio. Following Lozinski et al. (2009) [13], a quadratic time reconstruction is used.

Elsevier2009

Chargement

Numerical analysis of optimization-constrained differential equations

Chantal Landry

The modeling of a system composed by a gas phase and organic aerosol particles, and its numerical resolution are studied. The gas-aerosol system is modeled by ordinary differential equations coupled with a mixed-constrained optimization problem. This coupling induces discontinuities when inequality constraints are activated or deactivated. Two approaches for the solution of the optimization-constrained differential equations are presented. The first approach is a time splitting scheme together with a fixed-point method that alternates between the differential and optimization parts. The ordinary differential equations are approximated by the Crank-Nicolson scheme and a primal-dual interior-point method combined with a warm-start strategy is used to solve the minimization problem. The second approach considers the set of equations as a system of differential algebraic equations after replacing the minimization problem by its first order optimality conditions. An implicit 5th-order Runge-Kutta method (RADAU5) is then used. Both approaches are completed by numerical techniques for the detection and computation of the events (activation and deactivation of inequality constraints) when the system evolves in time. The computation of the events is based on continuation techniques and geometric arguments. Moreover the first approach completes the computation with extrapolation polynomials and sensitivity analysis, whereas the second approach uses dense output formulas. Numerical results for gas-aerosol system made of several chemical species are proposed for both approaches. These examples show the efficiency and accuracy of each method. They also indicate that the second approach is more efficient than the first one. Furthermore theoretical examples show that the method for the computation of the activation is of second order for the first approach and exact for the second one.

EPFL2009

Chargement