Numerical Simulation of Flows in Domains with Moving Boundaries
Participants
Dipl.Math. Martin Engel, Prof. Dr. Michael Griebel
Description
In many applications of Computational Fluid Dynamics (CFD) the domain
under consideration is not constant in time, but the boundaries are
moving. There exists a variety of different numerical methods for the
solution of flow problems in timedependent domains. One widespread
approach is to discretize the Eulerian formulation of the
NavierStokes equations on fixed structured grids. The actual fluid
domain is embedded into the grid and described by additional, in
general passively transported, quantities, e.g. massless marker
particles (Marker and Cell), scalar color functions (Volume of Fluid)
or levelsets of higherdimensional functions. Most of these
techniques allow for an easy tracking of the fluid domain even when it
undergoes changes in topology like merging or breaking up of
interfaces. They are well suited for the simulation of free surface or
multiphase flows. However, since the interface is only implicitly given,
prescription of boundary conditions is difficult and depends on the accurate
reconstruction of the interface.
In this project we use another approach, which corresponds to the
discretization of the Arbitrary Lagrange Euler formulation of the
NavierStokes equations. A fixed logical domain built of unit cubes
is mapped onto the fluid domain in physical space by a
timedependent coordinate transformation. The timederivative of this
mapping enters the transformed equations as the grid velocity. To
achieve a conservative discretization, additional socalled Geometric
Conservation Laws (GCL) have to be fulfilled.
We implemented a flow solver for the parallel solution of the
incompressible NavierStokes equations in three space dimensions. The solver
is based on a Chorin projection method for decoupling velocity and
pressure variables. The discretization is carried out on a staggered
mesh using a finite volume method and weighted gridoriented velocity
components as the primary variables. The GCL is reformulated and its
discretization is used to compute the mesh velocity. The whole
solution process is parallelized by decomposing the fluid domain in
several structured zones.
Examples
The above image demonstrates the blockstructured approach which is used to
handle complicated geometries and provides a means for parallelization as well.
Shown is a transported scalar quantity in a flow field around a cylinder.
Click the image to view the corresponding mpeganimation (mpg, 4.4 MB).
The following series of images shows colorcodings of the horizontal
velocity component for a flow through a channel with a moving
indentation. Only the part of the channel downstream from the
indentation is shown. Click on an image to view the mpeg animation (mpg, 350 KB).
References
[1] 
M. Engel, M. Griebel. Flow simulation on moving boundaryfitted grids and application to fluidstructure interaction problems. International Journal for Numerical Methods in Fluids, 2005. To Appear.

[2] 
M. Engel. Numerische Simulation von Strömungen in zeitabhängigen Gebieten und Anwendung auf FluidStrukturWechselwirkungsprobleme. Diplomarbeit, Institut für Angewandte Mathematik, Universität Bonn, Bonn, Germany, 2002.

Related Projects
 Numerical Simulation of FluidStructure Interaction