457 research outputs found
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\'{e} 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
Model-reduced variational fluid simulation
We present a model-reduced variational Eulerian integrator for incompressible fluids, which combines the efficiency gains of dimension reduction, the qualitative robustness of coarse spatial and temporal resolutions of geometric integrators, and the simplicity of sub-grid accurate boundary conditions on regular grids to deal with arbitrarily-shaped domains. At the core of our contributions is a functional map approach to fluid simulation for which scalar- and vector-valued eigenfunctions of the Laplacian operator can be easily used as reduced bases. Using a variational integrator in time to preserve liveliness and a simple, yet accurate embedding of the fluid domain onto a Cartesian grid, our model-reduced fluid simulator can achieve realistic animations in significantly less computational time than full-scale non-dissipative methods but without the numerical viscosity from which current reduced methods suffer. We also demonstrate the versatility of our approach by showing how it easily extends to magnetohydrodynamics and turbulence modeling in 2D, 3D and curved domains
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
A Discrete Geometric Optimal Control Framework for Systems with Symmetries
This paper studies the optimal motion control of
mechanical systems through a discrete geometric approach. At
the core of our formulation is a discrete Lagrange-d’Alembert-
Pontryagin variational principle, from which are derived discrete
equations of motion that serve as constraints in our optimization
framework. We apply this discrete mechanical approach to
holonomic systems with symmetries and, as a result, geometric
structure and motion invariants are preserved. We illustrate our
method by computing optimal trajectories for a simple model of
an air vehicle flying through a digital terrain elevation map, and
point out some of the numerical benefits that ensue
- …