29,067 research outputs found
A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes
In this paper we show and discuss the use of a versatile interaction
potential approach coupled with an immersed boundary method to simulate a
variety of flows involving deformable bodies. In particular, we focus on two
kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii)
flow in the left ventricle of the heart with either a mechanical or a natural
valve. Both examples have in common the two-way interaction of the flow with a
deformable interface or a membrane. The interaction potential approach (de
Tullio & Pascazio, Jou. Comp. Phys., 2016; Tanaka, Wada and Nakamura,
Computational Biomechanics, 2016) with minor modifications can be used to
capture the deformation dynamics in both classes of problems. We show that the
approach can be used to replicate the deformation dynamics of liquid-liquid
interfaces through the use of ad-hoc elastic constants. The results from our
simulations agree very well with previous studies on the deformation of drops
in standard flow configurations such as deforming drop in a shear flow or a
cross flow. We show that the same potential approach can also be used to study
the flow in the left ventricle of the heart. The flow imposed into the
ventricle interacts dynamically with the mitral valve (mechanical or natural)
and the ventricle which are simulated using the same model. Results from these
simulations are compared with ad- hoc in-house experimental measurements.
Finally, a parallelisation scheme is presented, as parallelisation is
unavoidable when studying large scale problems involving several thousands of
simultaneously deforming bodies on hundreds of distributed memory computing
processors
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
- …