1,723 research outputs found
Mass preserving distributed langrage multiplier approach to immersed boundary method
This research is devoted to mass conservation and CFL properties of the
Finite Elements Immersed Boundary Method. We first explore an enhanced higher order
scheme applied to the Finite Element Immersed Boundary Method technique introduced
by Boffi and Gastaldi. This technique is based on a Pointwise (PW) formulation of the
kinematic condition, and higher order elements show better conservation properties than
the original scheme. A further improvement with respect to the classical PW formulation
is achieved introducing a fully variational Distributed Lagrange Multiplier (DLM) formulation.
Numerical experiments show that DLM is not affected by any CFL condition.
Furthermore the mass conservation properties of this method are extremely competitive
A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies
We present a moving control volume (CV) approach to computing hydrodynamic
forces and torques on complex geometries. The method requires surface and
volumetric integrals over a simple and regular Cartesian box that moves with an
arbitrary velocity to enclose the body at all times. The moving box is aligned
with Cartesian grid faces, which makes the integral evaluation straightforward
in an immersed boundary (IB) framework. Discontinuous and noisy derivatives of
velocity and pressure at the fluid-structure interface are avoided and
far-field (smooth) velocity and pressure information is used. We re-visit the
approach to compute hydrodynamic forces and torques through force/torque
balance equation in a Lagrangian frame that some of us took in a prior work
(Bhalla et al., J Comp Phys, 2013). We prove the equivalence of the two
approaches for IB methods, thanks to the use of Peskin's delta functions. Both
approaches are able to suppress spurious force oscillations and are in
excellent agreement, as expected theoretically. Test cases ranging from Stokes
to high Reynolds number regimes are considered. We discuss regridding issues
for the moving CV method in an adaptive mesh refinement (AMR) context. The
proposed moving CV method is not limited to a specific IB method and can also
be used, for example, with embedded boundary methods
Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients
The distributed Lagrange multiplier/fictitious domain (DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients. The semi- and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface, where the arbitrary Lagrangian-Eulerian (ALE) technique is employed to deal with the moving and immersed subdomain. Stability and optimal convergence properties are obtained for both schemes. Numerical experiments are carried out for different scenarios of jump coefficients, and all theoretical results are validated
An adaptive fixed-mesh ALE method for free surface flows
In this work we present a Fixed-Mesh ALE method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain. This mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required. Some of the main ingredients of the formulation are the use of an Arbitrary-Lagrangian–Eulerian formulation for computing temporal derivatives, the use of stabilization terms for stabilizing convection, stabilizing the lack of compatibility between velocity and pressure interpolation spaces, and stabilizing the ill-conditioning introduced by the cuts on the background finite element mesh, and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments. Algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved. The method is tested in several numerical examples. The expected convergence rates both in space and time are observed. Smooth solution fields for both velocity and pressure are obtained (as a result of the contribution of the stabilization terms). Finally, a good agreement between the numerical results and the reference experimental data is obtained.Postprint (published version
A partition of unity approach to fluid mechanics and fluid-structure interaction
For problems involving large deformations of thin structures, simulating
fluid-structure interaction (FSI) remains challenging largely due to the need
to balance computational feasibility, efficiency, and solution accuracy.
Overlapping domain techniques have been introduced as a way to combine the
fluid-solid mesh conformity, seen in moving-mesh methods, without the need for
mesh smoothing or re-meshing, which is a core characteristic of fixed mesh
approaches. In this work, we introduce a novel overlapping domain method based
on a partition of unity approach. Unified function spaces are defined as a
weighted sum of fields given on two overlapping meshes. The method is shown to
achieve optimal convergence rates and to be stable for steady-state Stokes,
Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for
FSI in the case of a 2D mock aortic valve simulation. These initial results
point to the potential applicability of the method to a wide range of FSI
applications, enabling boundary layer refinement and large deformations without
the need for re-meshing or user-defined stabilization.Comment: 34 pages, 15 figur
- …