Unified Lagrangian formulation for solid and fluid mechanics and FSI problems

We present a Lagrangian monolithic strategy for solving fluid-structure interaction (FSI) problems. The formulation is called Unified because fluids and solids are solved using the same solution scheme and unknown variables. The method is based on a mixed velocity-pressure formulation. Each time step increment is solved via an iterative partitioned two-step procedure. The Particle Finite Element Method (PFEM) is used for solving the fluid parts of the domain, while for the solid ones the Finite Element Method (FEM) is employed. Both velocity and pressure fields are interpolated using linear shape functions. For quasiincompressible materials, the solution scheme is stabilized via the Finite Calculus (FIC) method. The stabilized elements for quasi-incompressible hypoelastic solids and Newtonian fluids are called VPS/S-element and VPS/F-element, respectively. Other two non-stabilized elements are derived for hypoelastic solids. One is based on a Velocity formulation (V-element) and the other on a mixed Velocity-Pressure scheme (VP-element). The algorithms for coupling the solid elements with the VPS/F fluid element are explained in detail. The Unified formulation is validated by solving benchmark FSI problems and by comparing the numerical solution to the ones published in the literature

3D numerical simulation of free-surface Bingham fluids interacting with structures using the PFEM

This paper presents a purely Lagrangian approach for the 3D simulation of Bingham free-surface uids and their interaction with deformable solid structures. In the proposed numerical strategy, the fluid is handled using the Particle Finite Element Method (PFEM) to tackle the issues resulting from extreme changes of geometry, such as mesh distortion and free-surface evolution. Additionally, the Papanastasiou model is employed as a regularization technique to overcome the computational difficulties associated with the classical Bingham model. The solid structure, on the other hand, is represented by the hypoelastic constitutive model and simulated using the conventional Finite Element Method (FEM). The coupling between the fluid and the structure is achieved via a monolithic approach, called Unified formulation. Several numerical examples are presented to illustrate the correctness and the robustness of the proposed formulation, in 2D and in 3D. Special attention is devoted to the analysis of the convergence behavior of the proposed computational framework, the effect of the regularization on the numerical results and the 3D effects. Moreover, detailed comparisons between the simulated results and experimental data are performed so that the concerned problems and results can serve as benchmarks

Framework for adaptive fluid-structure interaction with industrial applications

We present developments in the Unicorn-HPC framework for unified continuum mechanics, enabling adaptive finite element computation of fluid-structure interaction, and an overview of the larger FEniCS-HPC framework for automated solution of partial diffential equations of which Unicorn-HPC is a part. We formulate the basic model and finite element discretisation method and adaptive algorithms. We test the framework on a 2D model problem consisting of a flexible beam in channel flow, and to illustrate the capabilities of the computational framework, we show two application examples from industry and medicine. We simulate a flexible mixer plate in turbulent flow in an exhaust system where the target output is aeroacoustic quantities. The second example is a self-oscillating vocal fold configuration, where the ultimate goal is to predict how the voice is affected by physiological changes from aerodynamics. Here we give the displacement signal of a point on the folds

Coupled modeling for investigation of blast induced traumatic brain injury

Modeling of human body biomechanics resulting from blast exposure is very challenging because of the complex geometry and the substantially different materials involved. We have developed anatomy based high-fidelity finite element model (FEM) of the human body and finite volume model (FVM) of air around the human. The FEM model was used to accurately simulate the stress wave propagation in the human body under blast loading. The blast loading was generated by simulating C4 explosions, via a combination of 1-D and 3-D computational fluid dynamics (CFD) formulations. By employing the coupled Eulerian-Lagrangian fluid structure interaction (FSI) approach we obtained the parametric response of the human brain by the blast wave impact. We also developed the methodology to solve the strong interaction between cerebrospinal fluids (CSF) and the surrounding tissue for the closed-head impact. We presented both the arbitrary Lagrangian Eulerian (ALE) method and a new unified approach based on the material point method (MPM) to solve fluid dynamics and solid mechanics simultaneously. The accuracy and efficiency of ALE and MPM solvers for the skull-CSF-brain coupling problem was compared. The presented results suggest that the developed coupled models and techniques could be used to predict human biomechanical responses in blast events, and help design the protection against the blast induced TBI

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

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

A Modular Framework for Implicit 3D-0D Coupling in Cardiac Mechanics

In numerical simulations of cardiac mechanics, coupling the heart to a model of the circulatory system is essential for capturing physiological cardiac behavior. A popular and efficient technique is to use an electrical circuit analogy, known as a lumped parameter network or zero-dimensional (0D) fluid model, to represent blood flow throughout the cardiovascular system. Due to the strong physical interaction between the heart and the blood circulation, developing accurate and efficient numerical coupling methods remains an active area of research. In this work, we present a modular framework for implicitly coupling three-dimensional (3D) finite element simulations of cardiac mechanics to 0D models of blood circulation. The framework is modular in that the circulation model can be modified independently of the 3D finite element solver, and vice versa. The numerical scheme builds upon a previous work that combines 3D blood flow models with 0D circulation models (3D fluid - 0D fluid). Here, we extend it to couple 3D cardiac tissue mechanics models with 0D circulation models (3D structure - 0D fluid), showing that both mathematical problems can be solved within a unified coupling scheme. The effectiveness, temporal convergence, and computational cost of the algorithm are assessed through multiple examples relevant to the cardiovascular modeling community. Importantly, in an idealized left ventricle example, we show that the coupled model yields physiological pressure-volume loops and naturally recapitulates the isovolumic contraction and relaxation phases of the cardiac cycle without any additional numerical techniques. Furthermore, we provide a new derivation of the scheme inspired by the Approximate Newton Method of Chan (1985), explaining how the proposed numerical scheme combines the stability of monolithic approaches with the modularity and flexibility of partitioned approaches