Mesh moving techniques in fluid-structure interaction: robustness, accumulated distortion and computational efficiency

An important ingredient of any moving-mesh method for fluid-structure interaction (FSI) problems is the mesh deformation technique (MDT) used to adapt the computational mesh in the moving fluid domain. An ideal technique is computationally inexpensive, can handle large mesh deformations without inverting mesh elements and can sustain an FSI simulation for extensive periods ot time without irreversibly distorting the mesh. Here we compare several commonly used techniques based on the solution of elliptic partial differential equations, including harmonic extension, bi-harmonic extension and techniques based on the equations of linear elasticity. Moreover, we propose a novel technique which utilizes ideas from continuation methods to efficiently solve the equations of nonlinear elasticity and proves to be robust even when the mesh is subject to extreme deformations. In addition to that, we study how each technique performs when combined with the Jacobian-based local stiffening. We evaluate each technique on a popular two-dimensional FSI benchmark reproduced by using an isogeometric partitioned solver with strong coupling.

[1]  Volker John,et al.  Finite Element Methods for Incompressible Flow Problems , 2016 .

[2]  Xiao-Chuan Cai,et al.  A fully implicit domain decomposition based ALE framework for three-dimensional fluid-structure interaction with application in blood flow computation , 2014, J. Comput. Phys..

[3]  Tayfun E. Tezduyar,et al.  SPACE–TIME FLUID–STRUCTURE INTERACTION METHODS , 2012 .

[4]  N. Santos Numerical methods for fluid-structure interaction problems with valves , 2007 .

[5]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[6]  Tayfun E. Tezduyar,et al.  Ventricle-valve-aorta flow analysis with the Space–Time Isogeometric Discretization and Topology Change , 2020, Computational Mechanics.

[7]  Tayfun E. Tezduyar Finite Element Interface-Tracking and Interface-Capturing Techniques for Flows With Moving Boundaries and Interfaces , 2001, Heat Transfer: Volume 3 — Fluid-Physics and Heat Transfer for Macro- and Micro-Scale Gas-Liquid and Phase-Change Flows.

[8]  P. Wriggers Nonlinear Finite Element Methods , 2008 .

[9]  W. Wall,et al.  Fixed-point fluid–structure interaction solvers with dynamic relaxation , 2008 .

[10]  Bert Jüttler,et al.  Geometry + Simulation Modules: Implementing Isogeometric Analysis , 2014 .

[11]  S. Mittal,et al.  Computation of unsteady incompressible flows with the stabilized finite element methods: Space-time formulations, iterative strategies and massively parallel implementations , 1992 .

[12]  Tayfun E. Tezduyar,et al.  Modelling of fluid–structure interactions with the space–time finite elements: Solution techniques , 2007 .

[13]  Jens Gravesen,et al.  Planar Parametrization in Isogeometric Analysis , 2012, MMCS.

[14]  Tayfun E. Tezduyar,et al.  Mesh update strategies in parallel finite element computations of flow problems with moving boundaries and interfaces , 1994 .

[16]  Kenji Takizawa,et al.  Fluid–structure interaction modeling of clusters of spacecraft parachutes with modified geometric porosity , 2013 .

[17]  Tayfun E. Tezduyar,et al.  A low-distortion mesh moving method based on fiber-reinforced hyperelasticity and optimized zero-stress state , 2020, Computational Mechanics.

[18]  John A. Evans,et al.  An Isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD surfaces , 2012 .

[19]  Marek Behr,et al.  Parallel finite-element computation of 3D flows , 1993, Computer.

[20]  Kenji Takizawa,et al.  FSI analysis of the blood flow and geometrical characteristics in the thoracic aorta , 2014 .

[21]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[22]  Olaf Schenk,et al.  Toward the Next Generation of Multiperiod Optimal Power Flow Solvers , 2018, IEEE Transactions on Power Systems.

[23]  Giancarlo Sangalli,et al.  Analysis-suitable G1 multi-patch parametrizations for C1 isogeometric spaces , 2016, Comput. Aided Geom. Des..

[24]  B. Simeon,et al.  A hierarchical approach to adaptive local refinement in isogeometric analysis , 2011 .

[25]  Alfio Quarteroni,et al.  Cardiovascular mathematics : modeling and simulation of the circulatory system , 2009 .

[26]  T. Wick,et al.  Finite elements for fluid–structure interaction in ALE and fully Eulerian coordinates , 2010 .

[27]  Bernd Simeon,et al.  Isogeometric parametrization inspired by large elastic deformation , 2018, Computer Methods in Applied Mechanics and Engineering.

[28]  Victor M. Calo,et al.  Isogeometric Analysis of Hyperelastic Materials Using PetIGA , 2013, ICCS.

[29]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[30]  T. Tezduyar,et al.  Mesh Moving Techniques for Fluid-Structure Interactions With Large Displacements , 2003 .

[31]  S. Turek,et al.  Proposal for Numerical Benchmarking of Fluid-Structure Interaction between an Elastic Object and Laminar Incompressible Flow , 2006 .

[32]  M. Fortin,et al.  Mixed Finite Element Methods and Applications , 2013 .

[33]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[34]  P. Deuflhard Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms , 2011 .

[35]  Yuri Bazilevs,et al.  Computational Fluid-Structure Interaction: Methods and Applications , 2013 .