On introduit ici une classe de schémas symplectiques pour l'intégration numérique de systèmes hamiltoniens hautement oscillants. L'approche est basée sur la formulation Hamilton–Jacobi des équations du mouvement. En appliquant un développement à deux échelles à la solution des équations de Hamilton–Jacobi elle-même, on obtient ainsi, via la fonction génératrice, une classe de schémas symplectiques par construction. Un exemple de schéma ainsi construit est présenté ici sur un cas test habituel de système hautement oscillant, démontrant l'efficacité de l'approche. La dérivation d'autres schémas, et leurs tests dans des situations plus générales, feront l'objet d'une autre publication.
We introduce here a class of symplectic schemes for the numerical integration of highly oscillatory Hamiltonian systems. The bottom line for the approach is to exploit the Hamilton–Jacobi form of the equations of motion. Because we perform a two-scale expansion of the solution of the Hamilton–Jacobi equations itself, we readily obtain, after an appropriate discretization, a symplectic integration scheme. An example of such an integration scheme, following the general approach, is presented here on a specific commonly used test case. The efficiency of the approach is demonstrated. Further developments will be presented elsewhere.
Accepté le :
Publié le :
@article{CRMATH_2007__344_4_277_0,
author = {Le Bris, Claude and Legoll, Fr\'ed\'eric},
title = {D\'erivation de sch\'emas num\'eriques symplectiques pour des syst\`emes hamiltoniens hautement oscillants},
journal = {Comptes Rendus. Math\'ematique},
pages = {277--282},
publisher = {Elsevier},
volume = {344},
number = {4},
year = {2007},
doi = {10.1016/j.crma.2006.12.012},
language = {fr},
url = {http://www.numdam.org/articles/10.1016/j.crma.2006.12.012/}
}
TY - JOUR AU - Le Bris, Claude AU - Legoll, Frédéric TI - Dérivation de schémas numériques symplectiques pour des systèmes hamiltoniens hautement oscillants JO - Comptes Rendus. Mathématique PY - 2007 SP - 277 EP - 282 VL - 344 IS - 4 PB - Elsevier UR - http://www.numdam.org/articles/10.1016/j.crma.2006.12.012/ DO - 10.1016/j.crma.2006.12.012 LA - fr ID - CRMATH_2007__344_4_277_0 ER -
%0 Journal Article %A Le Bris, Claude %A Legoll, Frédéric %T Dérivation de schémas numériques symplectiques pour des systèmes hamiltoniens hautement oscillants %J Comptes Rendus. Mathématique %D 2007 %P 277-282 %V 344 %N 4 %I Elsevier %U http://www.numdam.org/articles/10.1016/j.crma.2006.12.012/ %R 10.1016/j.crma.2006.12.012 %G fr %F CRMATH_2007__344_4_277_0
Le Bris, Claude; Legoll, Frédéric. Dérivation de schémas numériques symplectiques pour des systèmes hamiltoniens hautement oscillants. Comptes Rendus. Mathématique, Tome 344 (2007) no. 4, pp. 277-282. doi : 10.1016/j.crma.2006.12.012. http://www.numdam.org/articles/10.1016/j.crma.2006.12.012/
[1] Homogenization of hamiltonian systems with a strong constraining potential, Phys. D, Volume 102 (1997), pp. 57-77
[2] F. Castella, Ph. Chartier, E. Faou, en préparation
[3] F. Castella, Ph. Chartier, E. Faou, C. Le Bris, F. Legoll, A. Murua, travaux en cours
[4] Difference schemes for Hamiltonian formalism and symplectic geometry, J. Comp. Math., Volume 4 (1986), pp. 279-289
[5] Error analysis of exponential integrators for oscillatory second-order differential equations, J. Phys. A, Volume 39 (2006), pp. 5495-5507
[6] Generalized Verlet algorithm for efficient molecular dynamics simulations with long range interaction, Molecular Simulation, Volume 6 (1991) no. 1–3, pp. 121-142
[7] Geometric Numerical Integration, Springer, 2006
[8] Reversible multiple time scale molecular dynamics, J. Chem. Phys., Volume 97 (1992), pp. 1990-2001
Cité par Sources :
⁎ Les auteurs ont bénéficié du support financier du Ministère de la recherche et des nouvelles technologies, dans le cadre de l'action « Nouvelles interfaces des mathématiques » SIMUMOL, et de l'Agence nationale de la recherche, programme non thématique INGEMOL. Ils remercient Philippe Chartier pour ses commentaires dans l'élaboration de ce texte.
