6,476 research outputs found
Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement
We present a new model and a novel loosely coupled partitioned numerical
scheme modeling fluid-structure interaction (FSI) in blood flow allowing
non-zero longitudinal displacement. Arterial walls are modeled by a {linearly
viscoelastic, cylindrical Koiter shell model capturing both radial and
longitudinal displacement}. Fluid flow is modeled by the Navier-Stokes
equations for an incompressible, viscous fluid. The two are fully coupled via
kinematic and dynamic coupling conditions. Our numerical scheme is based on a
new modified Lie operator splitting that decouples the fluid and structure
sub-problems in a way that leads to a loosely coupled scheme which is
{unconditionally} stable. This was achieved by a clever use of the kinematic
coupling condition at the fluid and structure sub-problems, leading to an
implicit coupling between the fluid and structure velocities. The proposed
scheme is a modification of the recently introduced "kinematically coupled
scheme" for which the newly proposed modified Lie splitting significantly
increases the accuracy. The performance and accuracy of the scheme were studied
on a couple of instructive examples including a comparison with a monolithic
scheme. It was shown that the accuracy of our scheme was comparable to that of
the monolithic scheme, while our scheme retains all the main advantages of
partitioned schemes, such as modularity, simple implementation, and low
computational costs
Stationary shapes of deformable particles moving at low Reynolds numbers
Lecture Notes of the Summer School ``Microswimmers -- From Single Particle
Motion to Collective Behaviour'', organised by the DFG Priority Programme SPP
1726 (Forschungszentrum J{\"{u}}lich, 2015).Comment: Pages C7.1-16 of G. Gompper et al. (ed.), Microswimmers - From Single
Particle Motion to Collective Behaviour, Lecture Notes of the DFG SPP 1726
Summer School 2015, Forschungszentrum J\"ulich GmbH, Schriften des
Forschungszentrums J\"ulich, Reihe Key Technologies, Vol 110, ISBN
978-3-95806-083-
A computational framework for the morpho-elastic development of molluskan shells by surface and volume growth
Mollusk shells are an ideal model system for understanding the morpho-elastic
basis of morphological evolution of invertebrates' exoskeletons. During the
formation of the shell, the mantle tissue secretes proteins and minerals that
calcify to form a new incremental layer of the exoskeleton. Most of the
existing literature on the morphology of mollusks is descriptive. The
mathematical understanding of the underlying coupling between pre-existing
shell morphology, de novo surface deposition and morpho-elastic volume growth
is at a nascent stage, primarily limited to reduced geometric representations.
Here, we propose a general, three-dimensional computational framework coupling
pre-existing morphology, incremental surface growth by accretion, and
morpho-elastic volume growth. We exercise this framework by applying it to
explain the stepwise morphogenesis of seashells during growth: new material
surfaces are laid down by accretive growth on the mantle whose form is
determined by its morpho-elastic growth. Calcification of the newest surfaces
extends the shell as well as creates a new scaffold that constrains the next
growth step. We study the effects of surface and volumetric growth rates, and
of previously deposited shell geometries on the resulting modes of mantle
deformation, and therefore of the developing shell's morphology. Connections
are made to a range of complex shells ornamentations.Comment: Main article is 20 pages long with 15 figures. Supplementary material
is 4 pages long with 6 figures and 6 attached movies. To be published in PLOS
Computational Biolog
A monolithic fluid-structure interaction formulation for solid and liquid membranes including free-surface contact
A unified fluid-structure interaction (FSI) formulation is presented for
solid, liquid and mixed membranes. Nonlinear finite elements (FE) and the
generalized-alpha scheme are used for the spatial and temporal discretization.
The membrane discretization is based on curvilinear surface elements that can
describe large deformations and rotations, and also provide a straightforward
description for contact. The fluid is described by the incompressible
Navier-Stokes equations, and its discretization is based on stabilized
Petrov-Galerkin FE. The coupling between fluid and structure uses a conforming
sharp interface discretization, and the resulting non-linear FE equations are
solved monolithically within the Newton-Raphson scheme. An arbitrary
Lagrangian-Eulerian formulation is used for the fluid in order to account for
the mesh motion around the structure. The formulation is very general and
admits diverse applications that include contact at free surfaces. This is
demonstrated by two analytical and three numerical examples exhibiting strong
coupling between fluid and structure. The examples include balloon inflation,
droplet rolling and flapping flags. They span a Reynolds-number range from
0.001 to 2000. One of the examples considers the extension to rotation-free
shells using isogeometric FE.Comment: 38 pages, 17 figure
Hybrid finite difference/finite element immersed boundary method
The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian
description of the structural deformations, stresses, and forces along with an Eulerian description of the
momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary
methods described immersed elastic structures using systems of flexible fibers, and even now, most
immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This
work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian
variables that facilitates independent spatial discretizations for the structure and background grid. This
approach employs a finite element discretization of the structure while retaining a finite difference scheme
for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively
contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases
in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors
that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes.
The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse
structural meshes with the immersed boundary method. This work also contrasts two different weak forms
of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations
facilitated by our coupling approach
- …