172 research outputs found
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
Discrete exterior calculus (DEC) for the surface Navier-Stokes equation
We consider a numerical approach for the incompressible surface Navier-Stokes
equation. The approach is based on the covariant form and uses discrete
exterior calculus (DEC) in space and a semi-implicit discretization in time.
The discretization is described in detail and related to finite difference
schemes on staggered grids in flat space for which we demonstrate second order
convergence. We compare computational results with a vorticity-stream function
approach for surfaces with genus 0 and demonstrate the interplay between
topology, geometry and flow properties. Our discretization also allows to
handle harmonic vector fields, which we demonstrate on a torus.Comment: 21 pages, 9 figure
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
Discrete Lie Advection of Differential Forms
In this paper, we present a numerical technique for performing Lie advection
of arbitrary differential forms. Leveraging advances in high-resolution finite
volume methods for scalar hyperbolic conservation laws, we first discretize the
interior product (also called contraction) through integrals over Eulerian
approximations of extrusions. This, along with Cartan's homotopy formula and a
discrete exterior derivative, can then be used to derive a discrete Lie
derivative. The usefulness of this operator is demonstrated through the
numerical advection of scalar fields and 1-forms on regular grids.Comment: Accepted version; to be published in J. FoC
Weak Form of Stokes-Dirac Structures and Geometric Discretization of Port-Hamiltonian Systems
We present the mixed Galerkin discretization of distributed parameter
port-Hamiltonian systems. On the prototypical example of hyperbolic systems of
two conservation laws in arbitrary spatial dimension, we derive the main
contributions: (i) A weak formulation of the underlying geometric
(Stokes-Dirac) structure with a segmented boundary according to the causality
of the boundary ports. (ii) The geometric approximation of the Stokes-Dirac
structure by a finite-dimensional Dirac structure is realized using a mixed
Galerkin approach and power-preserving linear maps, which define minimal
discrete power variables. (iii) With a consistent approximation of the
Hamiltonian, we obtain finite-dimensional port-Hamiltonian state space models.
By the degrees of freedom in the power-preserving maps, the resulting family of
structure-preserving schemes allows for trade-offs between centered
approximations and upwinding. We illustrate the method on the example of
Whitney finite elements on a 2D simplicial triangulation and compare the
eigenvalue approximation in 1D with a related approach.Comment: Copyright 2018. This manuscript version is made available under the
CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0
Interactive Blood-Coil Simulation in Real-time during Aneurysm Embolization
International audienceOver the last decade, remarkable progress has been made in the field of endovascular treatment of aneurysms. Technological advances continue to make it possible for a growing number of patients with cerebral aneurysms to be treated with a variety of endovascular strategies, essentially using detachable platinum coils. Yet, coil embolization remains a very complex medical procedure for which careful planning must be combined with advanced technical skills in order to be successful. In this paper, we describe a complete process for patient-specific simulations of coil embolization, from mesh generation with medical datasets to computation of coil-flow bilateral influence. We propose a new method for simulating the complex blood flow patterns that take place within the aneurysm, and for simulating the interaction of coils with this flow. This interaction is twofold, first involving the impact of the flow on the coil during the initial stages of its deployment, and second concerning the decrease of blood velocity within the aneurysm, as a consequence of coil packing. We also propose an approach to achieve real-time computation of coil-flow bilateral influence, necessary for interactive simulation. This in turns allows to dynamically plan coil embolization for two key steps of the procedure: choice and placement of the first coils, and assessment of the number of coils necessary to reduce aneurysmal blood velocity and wall pressure. Finally, we provide the blood flow simulation results on several aneurysms with interesting clinical characteristics both in 2D and 3D, as well as comparisons with a commercial package for validation. The coil embolization procesure is simulated within an aneurysm, and pre- and post-operative status is reported
- …