TY - JOUR
T1 - Benchmarking the Immersed Finite Element Method for Fluid-Structure Interaction Problems
JF - Computers and Mathematics with Applications 69 (2015) 1167–1188
Y1 - 2015
A1 - Roy Saswati
A1 - Luca Heltai
A1 - Francesco Costanzo
AB - We present an implementation of a fully variational formulation of an immersed methods for fluid-structure interaction problems based on the finite element method. While typical implementation of immersed methods are characterized by the use of approximate Dirac delta distributions, fully variational formulations of the method do not require the use of said distributions. In our implementation the immersed solid is general in the sense that it is not required to have the same mass density and the same viscous response as the surrounding fluid. We assume that the immersed solid can be either viscoelastic of differential type or hyperelastic. Here we focus on the validation of the method via various benchmarks for fluid-structure interaction numerical schemes. This is the first time that the interaction of purely elastic compressible solids and an incompressible fluid is approached via an immersed method allowing a direct comparison with established benchmarks.
U1 - 34633
U2 - Mathematics
U4 - 1
U5 - MAT/08
ER -
TY - RPRT
T1 - A Fully Coupled Immersed Finite Element Method for Fluid Structure Interaction via the Deal.II Library
Y1 - 2012
A1 - Luca Heltai
A1 - Saswati Roy
A1 - Francesco Costanzo
KW - Finite Element Method
KW - Immersed Boundary Method
KW - Immersed Finite Element Method
AB - We present the implementation of a solution scheme for fluid-structure\\r\\ninteraction problems via the finite element software library deal.II. The\\r\\nsolution scheme is an immersed finite element method in which two independent discretizations are used for the fluid and immersed deformable body. In this type of formulation the support of the equations of motion of the fluid is extended to cover the union of the solid and fluid domains. The equations of motion over the extended solution domain govern the flow of a fluid under the action of a body force field. This body force field informs the fluid of the presence of the immersed solid. The velocity field of the immersed solid is the restriction over the immersed domain of the velocity field in the extended equations of motion. The focus of this paper is to show how the determination of the motion of the immersed domain is carried out in practice. We show that our implementation is general, that is, it is not dependent on a specific choice of the finite element spaces over the immersed solid and the extended fluid domains. We present some preliminary results concerning the accuracy of the proposed method.
PB - SISSA
UR - http://hdl.handle.net/1963/6255
N1 - 28 pages, 9 figures
U1 - 6172
U2 - Mathematics
U3 - Functional Analysis and Applications
U4 - 1
U5 - MAT/08 ANALISI NUMERICA
ER -
TY - JOUR
T1 - Variational implementation of immersed finite element methods
JF - Computer Methods in Applied Mechanics and Engineering. Volume 229-232, 1 July 2012, Pages 110-127
Y1 - 2012
A1 - Luca Heltai
A1 - Francesco Costanzo
KW - Turbulent flow
AB - Dirac-delta distributions are often crucial components of the solid-fluid coupling operators in immersed solution methods for fluid-structure interaction (FSI) problems. This is certainly so for methods like the Immersed Boundary Method (IBM) or the Immersed Finite Element Method (IFEM), where Dirac-delta distributions are approximated via smooth functions. By contrast, a truly variational formulation of immersed methods does not require the use of Dirac-delta distributions, either formally or practically. This has been shown in the Finite Element Immersed Boundary Method (FEIBM), where the variational structure of the problem is exploited to avoid Dirac-delta distributions at both the continuous and the discrete level. In this paper, we generalize the FEIBM to the case where an incompressible Newtonian fluid interacts with a general hyperelastic solid. Specifically, we allow (i) the mass density to be different in the solid and the fluid, (ii) the solid to be either viscoelastic of differential type or purely elastic, and (iii) the solid to be and either compressible or incompressible. At the continuous level, our variational formulation combines the natural stability estimates of the fluid and elasticity problems. In immersed methods, such stability estimates do not transfer to the discrete level automatically due to the non- matching nature of the finite dimensional spaces involved in the discretization. After presenting our general mathematical framework for the solution of FSI problems, we focus in detail on the construction of natural interpolation operators between the fluid and the solid discrete spaces, which guarantee semi-discrete stability estimates and strong consistency of our spatial discretization.

PB - Elsevier
UR - http://hdl.handle.net/1963/6462
N1 - 42 pages, 5 figures, Revision 1
U1 - 6389
U2 - Mathematics
U4 - 1
U5 - MAT/05 ANALISI MATEMATICA
ER -