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.
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.
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.
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.
To facilitate the construction of scalable, parallelizeable implementations, we have developed an extension class library for PETSc [8], the Dynamical Systems Toolkit (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.
","abstract_html":"<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>","abstract_has_math":false,"creators":["Karpeev, Dmitry"],"institution":null,"degree_name":"Doctor of Philosophy (PhD)","degree_level":"Dissertation","degree_discipline":"Computer Science","degree_department":null,"school":null,"contributors":["David Keyes","Constance Schober","Chester Grosch","D. Glenn Lasseigne","Alex Pothen"],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2002,"date_issued":"2002-01-01T08:00:00Z","date_published":"2002-01-01T08:00:00Z","updated_at":"2026-07-24T03:34:11Z","subjects":["Geometric integrators","Hamiltonian PDEs","PDEs","Symplectic","Computer Sciences","Mathematics"],"languages":[],"rights":["In Copyright. URI: http://rightsstatements.org/vocab/InC/1.0/ 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).
"],"rights_urls":[],"identifier_entries":[{"key":"dc:identifier","label":"Identifier","values":["9780493883670"],"render_values":[{"text":"9780493883670","href":null,"code":true}]}]},"links":{"outbound_url":"https://digitalcommons.odu.edu/computerscience_etds/76","outbound_label":"Repository record","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["David Keyes","Constance Schober","Chester Grosch","D. Glenn Lasseigne","Alex Pothen"]},{"key":"dc:creator","label":"Author","values":["Karpeev, Dmitry"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.available","label":"Dc Date Available","values":["2019-02-19T08:00:00Z"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Computer Science"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Doctor of Philosophy (PhD)"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Geometric integrators","Hamiltonian PDEs","PDEs","Symplectic","Computer Sciences","Mathematics"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:rights","label":"Dc Rights","values":["In Copyright. URI: http://rightsstatements.org/vocab/InC/1.0/ 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).
"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["9780493883670","https://digitalcommons.odu.edu/computerscience_etds/76"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description.abstract","label":"Abstract","values":["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.
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.
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.
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.
To facilitate the construction of scalable, parallelizeable implementations, we have developed an extension class library for PETSc [8], the Dynamical Systems Toolkit (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.
"]},{"key":"dc:title","label":"Title","values":["Geometric Integrators for Hamiltonian PDEs"]}]}],"canonical_facts":{"dc:contributor":["David Keyes","Constance Schober","Chester Grosch","D. Glenn Lasseigne","Alex Pothen"],"dc:creator":["Karpeev, Dmitry"],"dc:date.available":["2019-02-19T08:00:00Z"],"dc:description.abstract":["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.
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.
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.
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.
To facilitate the construction of scalable, parallelizeable implementations, we have developed an extension class library for PETSc [8], the Dynamical Systems Toolkit (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.
"],"dc:identifier":["9780493883670","https://digitalcommons.odu.edu/computerscience_etds/76"],"dc:rights":["In Copyright. URI: http://rightsstatements.org/vocab/InC/1.0/ 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).
"],"dc:subject":["Geometric integrators","Hamiltonian PDEs","PDEs","Symplectic","Computer Sciences","Mathematics"],"dc:title":["Geometric Integrators for Hamiltonian PDEs"],"thesis:degree_discipline":["Computer Science"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Doctor of Philosophy (PhD)"]},"updated_at":"2026-07-24T03:34:11Z"}