BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jerrold E. Marsden

Jerrold E. Marsden appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.comp-phmath-phmath.MPmath.DSmath.DGComputationmath.APmath.OCmath.SG

Catalog footprint

What is connected

15works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2011arXiv

From efficient symplectic exponentiation of matrices to symplectic integration of high-dimensional Hamiltonian systems with slowly varying quadratic stiff potentials

We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based on the highly-non-trivial introduction of two efficient symplectic schemes for exponentiations of matrices that only require O(n) matrix multiplications operations at each coarse time step for a preset small number n. The proposed integrator is shown to be (i) uniformly convergent on positions; (ii) symplectic in both slow and fast variables; (iii) well adapted to high dimensional systems. Our framework also provides a general method for iteratively exponentiating a slowly varying sequence of (possibly high dimensional) matrices in an efficient way.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Space-time FLAVORS: finite difference, multisymlectic, and pseudospectral integrators for multiscale PDEs

We present a new class of integrators for stiff PDEs. These integrators are generalizations of FLow AVeraging integratORS (FLAVORS) for stiff ODEs and SDEs introduced in [Tao, Owhadi and Marsden 2010] with the following properties: (i) Multiscale: they are based on flow averaging and have a computational cost determined by mesoscopic steps in space and time instead of microscopic steps in space and time; (ii) Versatile: the method is based on averaging the flows of the given PDEs (which may have hidden slow and fast processes). This bypasses the need for identifying explicitly (or numerically) the slow variables or reduced effective PDEs; (iii) Nonintrusive: A pre-existing numerical scheme resolving the microscopic time scale can be used as a black box and easily turned into one of the integrators in this paper by turning the large coefficients on over a microscopic timescale and off during a mesoscopic timescale; (iv) Convergent over two scales: strongly over slow processes and in the sense of measures over fast ones; (v) Structure-preserving: for stiff Hamiltonian PDEs (possibly on manifolds), they can be made to be multi-symplectic, symmetry-preserving (symmetries are group actions that leave the system invariant) in all variables and variational.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variational integrators for electric circuits

In this contribution, we develop a variational integrator for the simulation of (stochastic and multiscale) electric circuits. When considering the dynamics of an electrical circuit, one is faced with three special situations: 1. The system involves external (control) forcing through external (controlled) voltage sources and resistors. 2. The system is constrained via the Kirchhoff current (KCL) and voltage laws (KVL). 3. The Lagrangian is degenerate. Based on a geometric setting, an appropriate variational formulation is presented to model the circuit from which the equations of motion are derived. A time-discrete variational formulation provides an iteration scheme for the simulation of the electric circuit. Dependent on the discretization, the intrinsic degeneracy of the system can be canceled for the discrete variational scheme. In this way, a variational integrator is constructed that gains several advantages compared to standard integration tools for circuits; in particular, a comparison to BDF methods (which are usually the method of choice for the simulation of electric circuits) shows that even for simple LCR circuits, a better energy behavior and frequency spectrum preservation can be observed using the developed variational integrator.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometric, Variational Discretization of Continuum Theories

This study derives geometric, variational discretizations of continuum theories arising in fluid dynamics, magnetohydrodynamics (MHD), and the dynamics of complex fluids. A central role in these discretizations is played by the geometric formulation of fluid dynamics, which views solutions to the governing equations for perfect fluid flow as geodesics on the group of volume-preserving diffeomorphisms of the fluid domain. Inspired by this framework, we construct a finite-dimensional approximation to the diffeomorphism group and its Lie algebra, thereby permitting a variational temporal discretization of geodesics on the spatially discretized diffeomorphism group. The extension to MHD and complex fluid flow is then made through an appeal to the theory of Euler-Poincaré systems with advection, which provides a generalization of the variational formulation of ideal fluid flow to fluids with one or more advected parameters. Upon deriving a family of structured integrators for these systems, we test their performance via a numerical implementation of the update schemes on a cartesian grid. Among the hallmarks of these new numerical methods are exact preservation of momenta arising from symmetries, automatic satisfaction of solenoidal constraints on vector fields, good long-term energy behavior, robustness with respect to the spatial and temporal resolution of the discretization, and applicability to irregular meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-Dirac Structures and Hamilton-Pontryagin Principles for Lagrange-Dirac Field Theories

The purpose of this paper is to define the concept of multi-Dirac structures and to describe their role in the description of classical field theories. We begin by outlining a variational principle for field theories, referred to as the Hamilton-Pontryagin principle, and we show that the resulting field equations are the Euler-Lagrange equations in implicit form. Secondly, we introduce multi-Dirac structures as a graded analog of standard Dirac structures, and we show that the graph of a multisymplectic form determines a multi-Dirac structure. We then discuss the role of multi-Dirac structures in field theory by showing that the implicit field equations obtained from the Hamilton-Pontryagin principle can be described intrinsically using multi-Dirac structures. Furthermore, we show that any multi-Dirac structure naturally gives rise to a multi-Poisson bracket. We treat the case of field theories with nonholonomic constraints, showing that the integrability of the constraints is equivalent to the integrability of the underlying multi-Dirac structure. We finish with a number of illustrative examples, including time-dependent mechanics, nonlinear scalar fields and the electromagnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-intrusive and structure preserving multiscale integration of stiff ODEs, SDEs and Hamiltonian systems with hidden slow dynamics via flow averaging

We introduce a new class of integrators for stiff ODEs as well as SDEs. These integrators are (i) {\it Multiscale}: they are based on flow averaging and so do not fully resolve the fast variables and have a computational cost determined by slow variables (ii) {\it Versatile}: the method is based on averaging the flows of the given dynamical system (which may have hidden slow and fast processes) instead of averaging the instantaneous drift of assumed separated slow and fast processes. This bypasses the need for identifying explicitly (or numerically) the slow or fast variables (iii) {\it Nonintrusive}: A pre-existing numerical scheme resolving the microscopic time scale can be used as a black box and easily turned into one of the integrators in this paper by turning the large coefficients on over a microscopic timescale and off during a mesoscopic timescale (iv) {\it Convergent over two scales}: strongly over slow processes and in the sense of measures over fast ones. We introduce the related notion of two-scale flow convergence and analyze the convergence of these integrators under the induced topology (v) {\it Structure preserving}: for stiff Hamiltonian systems (possibly on manifolds), they can be made to be symplectic, time-reversible, and symmetry preserving (symmetries are group actions that leave the system invariant) in all variables. They are explicit and applicable to arbitrary stiff potentials (that need not be quadratic). Their application to the Fermi-Pasta-Ulam problems shows accuracy and stability over four orders of magnitude of time scales. For stiff Langevin equations, they are symmetry preserving, time-reversible and Boltzmann-Gibbs reversible, quasi-symplectic on all variables and conformally symplectic with isotropic friction.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stokes-Dirac Structures through Reduction of Infinite-Dimensional Dirac Structures

We consider the concept of Stokes-Dirac structures in boundary control theory proposed by van der Schaft and Maschke. We introduce Poisson reduction in this context and show how Stokes-Dirac structures can be derived through symmetry reduction from a canonical Dirac structure on the unreduced phase space. In this way, we recover not only the standard structure matrix of Stokes-Dirac structures, but also the typical non-canonical advection terms in (for instance) the Euler equation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structure preserving Stochastic Impulse Methods for stiff Langevin systems with a uniform global error of order 1 or 1/2 on position

Impulse methods are generalized to a family of integrators for Langevin systems with quadratic stiff potentials and arbitrary soft potentials. Uniform error bounds (independent from stiff parameters) are obtained on integrated positions allowing for coarse integration steps. The resulting integrators are explicit and structure preserving (quasi-symplectic for Langevin systems).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structure-Preserving Discretization of Incompressible Fluids

The geometric nature of Euler fluids has been clearly identified and extensively studied over the years, culminating with Lagrangian and Hamiltonian descriptions of fluid dynamics where the configuration space is defined as the volume-preserving diffeomorphisms, and Kelvin's circulation theorem is viewed as a consequence of Noether's theorem associated with the particle relabeling symmetry of fluid mechanics. However computational approaches to fluid mechanics have been largely derived from a numerical-analytic point of view, and are rarely designed with structure preservation in mind, and often suffer from spurious numerical artifacts such as energy and circulation drift. In contrast, this paper geometrically derives discrete equations of motion for fluid dynamics from first principles in a purely Eulerian form. Our approach approximates the group of volume-preserving diffeomorphisms using a finite dimensional Lie group, and associated discrete Euler equations are derived from a variational principle with non-holonomic constraints. The resulting discrete equations of motion yield a structure-preserving time integrator with good long-term energy behavior and for which an exact discrete Kelvin's circulation theorem holds.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Temperature and Friction Accelerated Sampling of Boltzmann-Gibbs Distribution

This paper is concerned with tuning friction and temperature in Langevin dynamics for fast sampling from the canonical ensemble. We show that near-optimal acceleration is achieved by choosing friction so that the local quadratic approximation of the Hamiltonian is a critical damped oscillator. The system is also over-heated and cooled down to its final temperature. The performances of different cooling schedules are analyzed as functions of total simulation time.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Dynamics of a Rigid Body in Potential Flow with Circulation

We consider the motion of a two-dimensional body of arbitrary shape in a planar irrotational, incompressible fluid with a given amount of circulation around the body. We derive the equations of motion for this system by performing symplectic reduction with respect to the group of volume-preserving diffeomorphisms and obtain the relevant Poisson structures after a further Poisson reduction with respect to the group of translations and rotations. In this way, we recover the equations of motion given for this system by Chaplygin and Lamb, and we give a geometric interpretation for the Kutta-Zhukowski force as a curvature-related effect. In addition, we show that the motion of a rigid body with circulation can be understood as a geodesic flow on a central extension of the special Euclidian group SE(2), and we relate the cocycle in the description of this central extension to a certain curvature tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Geometric Computational Electrodynamics with Variational Integrators and Discrete Differential Forms

In this paper, we develop a structure-preserving discretization of the Lagrangian framework for electromagnetism, combining techniques from variational integrators and discrete differential forms. This leads to a general family of variational, multisymplectic numerical methods for solving Maxwell's equations that automatically preserve key symmetries and invariants. In doing so, we demonstrate several new results, which apply both to some well-established numerical methods and to new methods introduced here. First, we show that Yee's finite-difference time-domain (FDTD) scheme, along with a number of related methods, are multisymplectic and derive from a discrete Lagrangian variational principle. Second, we generalize the Yee scheme to unstructured meshes, not just in space but in 4-dimensional spacetime. This relaxes the need to take uniform time steps, or even to have a preferred time coordinate at all. Finally, as an example of the type of methods that can be developed within this general framework, we introduce a new asynchronous variational integrator (AVI) for solving Maxwell's equations. These results are illustrated with some prototype simulations that show excellent energy and conservation behavior and lack of spurious modes, even for an irregular mesh with asynchronous time stepping.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Covariant Balance Laws in Continua with Microstructure

The purpose of this paper is to extend the Green-Naghdi-Rivlin balance of energy method to continua with microstructure. The key idea is to replace the group of Galilean transformations with the group of diffeomorphisms of the ambient space. A key advantage is that one obtains in a natural way all the needed balance laws on both the macro and micro levels along with two Doyle-Erickson formulas.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Controlled Lagrangians and Potential Shaping for Stabilization of Discrete Mechanical Systems

The method of controlled Lagrangians for discrete mechanical systems is extended to include potential shaping in order to achieve complete state-space asymptotic stabilization. New terms in the controlled shape equation that are necessary for matching in the discrete context are introduced. The theory is illustrated with the problem of stabilization of the cart-pendulum system on an incline. We also discuss digital and model predictive control.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Abstract mechanical connection and Abelian reconstruction for almost Kahler manifolds

When the phase space P of a Hamiltonian G-system (P, ω, G, J, H) has an almost Kahler structure a preferred connection, called abstract mechanical connection, can be defined by declaring horizontal spaces at each point to be metric orthogonal to the tangent to the group orbit. Explicit formulas for the corresponding connection one-form A are derived in terms of the momentum map, symplectic and complex structures. Such connection can play the role of the reconstruction connection (see A. Blaom, Reconstruction phases via Poisson reduction. Diff. Geom. Appl., 12:231--252, 2000.), thus significantly simplifying computations of the corresponding dynamic and geometric phases for an Abelian group G. These ideas are illustrated using the example of the resonant three-wave interaction. Explicit formulas for the connection one-form and the phases are given together with some new results on the symmetry reduction of the Poisson structure.

0 citations0 reviewsNo code linkedOpen access