Geometric Integrators for Hamiltonian PDEs

<p>We consider methods for systematic construction of algorithms for a class of time-dependent PDEs with Hamiltonian structure. These systems possess phase space geometry and constants of the motion that need to be preserved by the integration algorithm to reflect the qualitative features of the system.</p> <p>We exploit the structure of Hamiltonian systems, in particular their variational formulation based on a Lagrangian, and the dual covariant formulation, to expose the geometric features of the system that have natural analogs when discretized. We emphasize the local space-time approach to the constructions, making them amenable to parallelization and preconditioning using domain decomposition methods, and enabling treatment of complex spatial geometries and boundary conditions.</p> <p>These methods are applied to the two representative problems: the Nonlinear Schrödinger equation (NLS) and the Heisenberg magnet model (HM), both of which possess highly nontrivial geometric structures. We treat the “1+1” case of one spatial and one temporal dimension with periodic boundary conditions, but both systems have higher-dimensional “m+1” generalizations and the analyses extend accordingly.</p> <p>In addition, NLS has an integrable semidiscretization due to Ablowitz and Ladik (AL), which, as an ODE, possesses a highly nonlinear symplectic structure, For this system we consider a generating function approach to the nonstandard symplectic form, and derive novel symplectic integrators of arbitrary order exploiting the particular structure of AL where the standard techniques fail.</p> <p>To facilitate the construction of scalable, parallelizeable implementations, we have developed an extension class library for PETSc [8], the <em>Dynamical Systems Toolkit</em> (DST), that implements parallel mesh manipulation and discrete grid function and operator routines, which serve as the basic building blocks for efficient geometric integrators. While a piece of research code, it also facilitates rapid development and testing of algorithms obtained using the methods developed above and can serve as a foundation for further software development in this direction, In addition, the question of the best object framework for expression of natural PDE operations deserves investigation in its own right and we include its discussion.</p>

Name thesis:degree_name

Doctor of Philosophy (PhD)

Level thesis:degree_level

Dissertation

Discipline thesis:degree_discipline

Computer Science

Year dc:date.available

2002

Author dc:creator

• Karpeev, Dmitry

Contributors dc:contributor

• David Keyes
• Constance Schober
• Chester Grosch
• D. Glenn Lasseigne
• Alex Pothen

• Geometric integrators
• Hamiltonian PDEs
• PDEs
• Symplectic
• Computer Sciences
• Mathematics

Statement dc:rights

• <p>In Copyright. URI: <a href="http://rightsstatements.org/vocab/InC/1.0/">http://rightsstatements.org/vocab/InC/1.0/</a> This Item is protected by copyright and/or related rights. You are free to use this Item in any way that is permitted by the copyright and related rights legislation that applies to your use. For other uses you need to obtain permission from the rights-holder(s).</p>

Identifier

9780493883670

OAI identifier oai:identifier

oai:digitalcommons.odu.edu:computerscience_etds-1079