1,814 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
Multiphase flow of immiscible fluids on unstructured moving meshes
Figure 1: Multiple fluids with different viscosity coefficients and surface tension densities splashing on the bottom of a cylindrical container. Observe that the simulation has no problem dealing with thin sheets. In this paper, we present a method for animating multiphase flow of immiscible fluids using unstructured moving meshes. Our underlying discretization is an unstructured tetrahedral mesh, the deformable simplicial complex (DSC), that moves with the flow in a Lagrangian manner. Mesh optimization operations improve element quality and avoid element inversion. In the context of multiphase flow, we guarantee that every element is occupied by a single fluid and, consequently, the interface between fluids is represented by a set of faces in the simplicial complex. This approach ensures that the underlying discretization matches the physics and avoids the additional book-keeping required in grid-based methods where multiple fluids may occupy the same cell. Our Lagrangian approach naturally leads us to adopt a finite element approach to simulation, in contrast to the finite volume approaches adopted by a majority of fluid simulation techniques that use tetrahedral meshes. We characterize fluid simulation as an optimization problem allowing for full coupling of the pressure and velocity fields and the incorporation of a second-order surface energy. We introduce a preconditioner based on the diagonal Schur complement and solve our optimization on the GPU. We provide the results of parameter studies as well as
Arbitrary Lagrangian-Eulerian form of flowfield dependent variation (ALE-FDV) method for moving boundary problems
Flowfield Dependent Variation (FDV) method is a mixed explicit-implicit numerical
scheme that was originally developed to solve complex flow problems through the use
of so-called implicitness parameters. These parameters determine the implicitness of
FDV method by evaluating local gradients of physical flow parameters, hence vary
across the computational domain. The method has been used successfully in solving
wide range of flow problems. However it has only been applied to problems where the
objects or obstacles are static relative to the flow. Since FDV method has been proved
to be able to solve many complex flow problems, there is a need to extend FDV
method into the application of moving boundary problems where an object
experiences motion and deformation in the flow. With the main objective to develop a
robust numerical scheme that is applicable for wide range of flow problems involving
moving boundaries, in this study, FDV method was combined with a body
interpolation technique called Arbitrary Lagrangian-Eulerian (ALE) method. The
ALE method is a technique that combines Lagrangian and Eulerian descriptions of a
continuum in one numerical scheme, which then enables a computational mesh to
follow the moving structures in an arbitrary movement while the fluid is still seen in a
Eulerian manner. The new scheme, which is named as ALE-FDV method, is
formulated using finite volume method in order to give flexibility in dealing with
complicated geometries and freedom of choice of either structured or unstructured
mesh. The method is found to be conditionally stable because its stability is dependent
on the FDV parameters. The formulation yields a sparse matrix that can be solved by
using any iterative algorithm. Several benchmark stationary and moving body
problems in one, two and three-dimensional inviscid and viscous flows have been
selected to validate the method. Good agreement with available experimental and
numerical results from the published literature has been obtained. This shows that the
ALE-FDV has great potential for solving a wide range of complex flow problems
involving moving bodies
- …