Truly monolithic algebraic multigrid for fluid–structure interaction

The coupling of flexible structures to incompressible fluids draws a lot of attention during the last decade. Many different solution schemes have been proposed. In this contribution, we concentrate on the strong coupling fluid–structure interaction by means of monolithic solution schemes. Therein, a Newton–Krylov method is applied to the monolithic set of nonlinear equations. Such schemes require good preconditioning to be efficient. We propose two preconditioners that apply algebraic multigrid techniques to the entire fluid–structure interaction system of equations. The first is based on a standard block Gauss–Seidel approach, where approximate inverses of the individual field blocks are based on a algebraic multigrid hierarchy tailored for the type of the underlying physical problem. The second is based on a monolithic coarsening scheme for the coupled system that makes use of prolongation and restriction projections constructed for the individual fields. The resulting nonsymmetric monolithic algebraic multigrid method therefore involves coupling of the fields on coarse approximations to the problem yielding significantly enhanced performance. Copyright © 2010 John Wiley & Sons, Ltd.     

Robust and efficient monolithic fluid‐structure‐interaction solvers

We construct new robust and efficient preconditioned generalized minimal residual solvers for the monolithic linear systems of algebraic equations arising from the finite element discretization and Newton's linearization of the fully coupled fluid–structure interaction system of partial differential equations in the arbitrary Lagrangian–Eulerian formulation. We admit both linear elastic and nonlinear hyperelastic materials in the solid model and cover a large range of flows, for example, water, blood, and air, with highly varying density. The preconditioner is constructed in form of , where , , and are proper approximations to the matrices L, D, and U in the LDU block factorization of the fully coupled system matrix, respectively. The inverse of the corresponding Schur complement is approximated by applying a few cycles of a special class of algebraic multigrid methods to the perturbed fluid sub‐problem, which is obtained by modifying corresponding entries in the original fluid matrix with an explicitly constructed approximation to the exact perturbation coming from the sparse matrix–matrix multiplications. The numerical studies presented impressively demonstrate the robustness and the efficiency of the preconditioner proposed in the paper. Copyright © 2016 John Wiley & Sons, Ltd.     

A monolithic strategy for fluid–structure interaction problems

In this work, we present a new monolithic strategy for solving fluid–structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time‐stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.     

Topology optimization for stationary fluid–structure interaction problems using a new monolithic formulation

This paper outlines a new procedure for topology optimization in the steady‐state fluid–structure interaction (FSI) problem. A review of current topology optimization methods highlights the difficulties in alternating between the two distinct sets of governing equations for fluid and structure dynamics (hereafter, the fluid and structural equations, respectively) and in imposing coupling boundary conditions between the separated fluid and solid domains. To overcome these difficulties, we propose an alternative monolithic procedure employing a unified domain rather than separated domains, which is not computationally efficient. In the proposed analysis procedure, the spatial differential operator of the fluid and structural equations for a deformed configuration is transformed into that for an undeformed configuration with the help of the deformation gradient tensor. For the coupling boundary conditions, the divergence of the pressure and the Darcy damping force are inserted into the solid and fluid equations, respectively. The proposed method is validated in several benchmark analysis problems. Topology optimization in the FSI problem is then made possible by interpolating Young's modulus, the fluid pressure of the modified solid equation, and the inverse permeability from the damping force with respect to the design variables. Copyright © 2009 John Wiley & Sons, Ltd.     

A hierarchical sequential ALE poromechanics model for tire‐soil‐water interaction on fluid‐infiltrated roads

This paper introduces a hierarchical sequential arbitrary Lagrangian‐Eulerian (ALE) model for predicting the tire‐soil‐water interaction at finite deformations. Using the ALE framework, the interaction between a rolling pneumatic tire and the fluid‐infiltrated soil underneath will be captured numerically. The road is assumed to be a fully saturated two‐phase porous medium. The constitutive response of the tire and the solid skeleton of the porous medium is idealized as hyperelastic. Meanwhile, the interaction between tire, soil, and water will be simulated via a hierarchical operator‐split algorithm. A salient feature of the proposed framework is the steady state rolling framework. While the finite element mesh of the soil is fixed to a reference frame and moves with the tire, the solid and fluid constituents of the soil are flowing through the mesh in the ALE model according to the rolling speed of the tire. This treatment leads to an elegant and computationally efficient formulation to investigate the tire‐soil‐water interaction both close to the contact and in the far field. The presented ALE model for tire‐soil‐water interaction provides the essential basis for future applications, for example, to a path‐dependent frictional‐cohesive response of the consolidating soil and unsaturated soil, respectively. Copyright © 2017 John Wiley & Sons, Ltd.     

Improved fractional step method for simulating fluid‐structure interaction using the PFEM

Numerical difficulties are present in the particle finite element method even though it has been shown to be a powerful and effective approach to simulating fluid‐structure interaction. To overcome problems of mass loss on the free surface and the added‐mass effect, an improved fractional step method (FSM) that handles added‐mass terms in a mathematically exact way is developed. A further benefit is that no assumptions regarding the structural response are made in handling added‐mass terms, thus it is straightforward to incorporate material nonlinearity in fluid‐structure interaction (FSI) under this approach. Patch tests and comparisons with experimental data are presented in order to verify and validate the improved FSM for FSI applications. The computational cost of this approach is shown to be negligible compared with the other aspects of the FSM, particularly when the size of the structure and the fluid‐structure interface is small relative to the volume of fluid. Copyright © 2014 John Wiley & Sons, Ltd.     

Consistent arbitrary Lagrangian Eulerian formulation for large deformation thermo-mechanical analysis

Arbitrary Lagrangian Eulerian (ALE) method is widely used for simulation of large deformation problems, such as metal forming. However, in many such applications, modeling of the heat generation and transfer in conjunction with the stress analysis is necessary. In this work, a fully coupled dynamic ALE formulation is developed. The ALE form of energy balance equation is derived, and is coupled with the dynamic, rate dependent ALE stress analysis. The proposed formulation is used for simulation of a few thermo-mechanical problems. The effectiveness and efficiency of the ALE method is verified by comparing the results of this simulation with available experimental and numerical results.     

The fixed‐mesh ALE approach applied to solid mechanics and fluid–structure interaction problems

In this paper we propose a method to solve Solid Mechanics and fluid–structure interaction problems using always a fixed background mesh for the spatial discretization. The main feature of the method is that it properly accounts for the advection of information as the domain boundary evolves. To achieve this, we use an Arbitrary Lagrangian–Eulerian (ALE) framework, the distinctive characteristic being that at each time step results are projected onto a fixed, background mesh. For solid mechanics problems subject to large strains, the fixed‐mesh (FM)‐ALE method avoids the element stretching found in fully Lagrangian approaches. For FSI problems, FM‐ALE allows for the use of a single background mesh to solve both the fluid and the structure. Copyright © 2009 John Wiley & Sons, Ltd.     

A distributed memory parallel implementation of the multigrid method for solving three‐dimensional implicit solid mechanics problems

We describe the parallel implementation of a multigrid method for unstructured finite element discretizations of solid mechanics problems. We focus on a distributed memory programming model and use the MPI library to perform the required interprocessor communications. We present an algebraic framework for our parallel computations, and describe an object‐based programming methodology using Fortran90. The performance of the implementation is measured by solving both fixed‐ and scaled‐size problems on three different parallel computers (an SGI Origin2000, an IBM SP2 and a Cray T3E). The code performs well in terms of speedup, parallel efficiency and scalability. However, the floating point performance is considerably below the peak values attributed to these machines. Lazy processors are documented on the Origin that produce reduced performance statistics. The solution of two problems on an SGI Origin2000, an IBM PowerPC SMP and a Linux cluster demonstrate that the algorithm performs well when applied to the unstructured meshes required for practical engineering analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

A new monolithic Newton‐multigrid‐based FEM solution scheme for large strain dynamic poroelasticity problems

This paper presents a new efficient monolithic finite element solution scheme to treat the set of PDEs governing a 2D, biphasic, saturated theory of porous media model with intrinsically coupled and incompressible solid and fluid constituents for infinitesimal and large elastic deformation. Our approach, which inherits some of its techniques from CFD, is characterized by the following aspects: (1) it only performs operator evaluation with no additional Gateaux derivatives. In particular, the computations of the time‐consuming material tangent matrix are not involved here; (2) it solves the non‐linear dynamic problem with no restriction on the strength of coupling; (3) it is more efficient than the linear u v p solver discussed in previous works; (4) it requires weaker derivatives, and hence, lower‐order FE can be tested; and (5) the boundary conditions are reduced, solution independent and more convenient to apply than in the old u v p formulation. For the purpose of validation and comparison, prototypical simulations including analytical solutions are carried out, and at the end, an adaptive time stepping procedure is introduced to handle the rapid change in the numbers of nonlinear iterations that may occur. Copyright © 2016 John Wiley & Sons, Ltd.     

A Multiscale/stabilized Formulation of the Incompressible Navier–Stokes Equations for Moving Boundary Flows and Fluid–structure Interaction

This paper presents a multiscale/stabilized finite element formulation for the incompressible Navier–Stokes equations written in an Arbitrary Lagrangian–Eulerian (ALE) frame to model flow problems that involve moving and deforming meshes. The new formulation is derived based on the variational multiscale method proposed by Hughes (Comput Methods Appl Mech Eng 127:387–401, 1995) and employed in Masud and Khurram in (Comput Methods Appl Mech Eng 193:1997–2018, 2006); Masud and Khurram in (Comput Methods Appl Mech Eng 195:1750–1777, 2006) to study advection dominated transport phenomena. A significant feature of the formulation is that the structure of the stabilization terms and the definition of the stabilization tensor appear naturally via the solution of the sub-grid scale problem. A mesh moving technique is integrated in this formulation to accommodate the motion and deformation of the computational grid, and to map the moving boundaries in a rational way. Some benchmark problems are shown, and simulations of an elastic beam undergoing large amplitude periodic oscillations in a viscous fluid domain are presented.     

A general conjugate‐gradient‐based predictor–corrector solver for explicit finite‐element contact

A Lagrange‐multiplier based approach is presented for the general solution of multi‐body contact within an explicit finite element framework. The technique employs an explicit predictor step to permit the detection of interpenetration and then utilizes a corrector step, whose solution is obtained with a pre‐conditioned matrix‐free conjugate gradient projection method, to determine the Lagrange multipliers necessary to eliminate the predicted penetration. The predictor–corrector algorithm is developed for deformable bodies based upon the central difference method, and for rigid bodies from momentum and energy conserving approaches. Both frictionless and Coulomb‐based frictional contact idealizations are addressed. The technique imposes no time‐step constraints and quickly mitigates velocity discontinuities across closed interfaces. Special attention is directed toward contact between rigid bodies. Algorithmic moment arms conserve the translational and angular momentums of the system in the absence of external loads. Elastic collisions are captured with a two‐phase predictor–corrector approach and a geometrically approximate velocity jump criterion. The first step solves the inelastic contact problem and identifies inactive constraints between rigid bodies, while the second step generates the necessary velocity jump condition on the active constraints. The velocity criterion is shown to algorithmically preserve the system kinetic energy for two unconstrained rigid bodies. Copyright © 1999 John Wiley & Sons, Ltd. This paper was produced under the auspices of the U.S. Government and it is therefore not subject to copyright in the U.S.     

A Cartesian ghost‐cell multigrid Poisson solver for incompressible flows

In this paper, a simple Cartesian ghost‐cell multigrid Poisson solver is proposed for simulating incompressible fluid flows. The flow field is discretized efficiently on a rectangular mesh, in which solid bodies are immersed. A small number of ghost mesh cells and their symmetric image cells are distributed in the vicinity of the solid boundary. With the aid of the ghost and image cells, the Dirichlet and Neumann boundary conditions can be implemented effectively. Chorin's fractional‐step projection method is adopted for the coupling of velocity and pressure for the solution of the Navier–Stokes equations. Point‐wise Gauss–Seidel iteration is used to solve the pressure Poisson equation. To speed up the convergence of the solution to the corresponding linear system, sub‐level coarse meshes embedded with ghost and image cells are also introduced and operated in a sequential V‐cycle. Several test cases including the classical ideal incompressible flow around a cylinder, a lid‐driven cavity flow and viscous flow past a fixed/rotating cylinder are presented to demonstrate the accuracy and efficiency of the current approach. Copyright © 2010 John Wiley & Sons, Ltd.     

A particle‐based moving interface method (PMIM) for modeling the large deformation of boundaries in soft matter systems

The mechanics of the interaction between a fluid and a soft interface undergoing large deformations appear in many places, such as in biological systems or industrial processes. We present an Eulerian approach that describes the mechanics of an interface and its interactions with a surrounding fluid via the so‐called Navier boundary condition. The interface is modeled as a curvilinear surface with arbitrary mechanical properties across which discontinuities in pressure and tangential fluid velocity can be accounted for using a modified version of the extended finite element method. The coupling between the interface and the fluid is enforced through the use of Lagrange multipliers. The tracking and evolution of the interface are then handled in a Lagrangian step with the grid‐based particle method. We show that this method is ideal to describe large membrane deformations and Navier boundary conditions on the interface with velocity/pressure discontinuities. The validity of the model is assessed by evaluating the numerical convergence for a axisymmetrical flow past a spherical capsule with various surface properties. We show the effect of slip length on the shear flow past a two‐dimensional capsule and simulate the compression of an elastic membrane lying on a viscous fluid substrate. Copyright © 2015 John Wiley & Sons, Ltd.     

Inexact Schwarz‐algebraic multigrid preconditioners for crack problems modeled by extended finite element methods

Traditional algebraic multigrid (AMG) preconditioners are not well suited for crack problems modeled by extended finite element methods (XFEM). This is mainly because of the unique XFEM formulations, which embed discontinuous fields in the linear system by addition of special degrees of freedom. These degrees of freedom are not properly handled by the AMG coarsening process and lead to slow convergence. In this paper, we proposed a simple domain decomposition approach that retains the AMG advantages on well‐behaved domains by avoiding the coarsening of enriched degrees of freedom. The idea was to employ a multiplicative Schwarz preconditioner where the physical domain was partitioned into “healthy” (or unfractured) and “cracked” subdomains. First, the “healthy” subdomain containing only standard degrees of freedom, was solved approximately by one AMG V‐cycle, followed by concurrent direct solves of “cracked” subdomains. This strategy alleviated the need to redesign special AMG coarsening strategies that can handle XFEM discretizations. Numerical examples on various crack problems clearly illustrated the superior performance of this approach over a brute force AMG preconditioner applied to the linear system. Copyright © 2011 John Wiley & Sons, Ltd.     

A monolithic computational approach to thermo‐structure interaction

In the present work, a monolithic solution approach for thermo‐structure interaction problems motivated by the challenging application of the behaviour of rocket nozzles is proposed. Structural and thermal fields are independently discretised via finite elements. The resulting system of equations is solved via a monolithic thermo‐structure interaction scheme, which is constructed by a block Gauss–Seidel preconditioner in combination with algebraic multigrid methods. The proposed method is tested for four numerical examples, the second Danilovskaya problem, a simplified rocket nozzle configuration, an internally loaded hollow sphere, and a fully three‐dimensional nozzle configuration of a subscale thrust chamber. Good agreement of the numerical results with results from the literature is observed. Furthermore, it is shown that the monolithic solution algorithm can handle the complete range of the parameter spectrum, whereas partitioned algorithms are limited to a certain parameter range only. Moreover, the monolithic algorithm exhibits improved efficiency and robustness compared to partitioned algorithms. Copyright © 2013 John Wiley & Sons, Ltd.     

A radial basis function algorithm for simplified fluid‐structure data transfer

In loosely‐coupled aeroelastic computation, the aerodynamic and elastomechanical models are based on different grids and eventually a simplified structural model, like a shell or stick, is considered. CFD tools are applied over the aerodynamic grid, whereas CSM tools are over the structural one. The meshes are usually non‐conforming; thus, three‐dimensional and non‐intrusive interpolation procedures are necessary to transfer structural deformations to the aerodynamic surface grid and aerodynamic loads to the structure. For that purpose, an interpolator based on radial basis functions has been developed. This tool does not require grid connectivities and can be applied to any three‐dimensional data. We have developed two strategies to deal with two special cases of interest. The first is when the structure is represented by a beam model and there is not a proper structural mesh close to the aerodynamic surface. The second is when managing a full configuration aircraft and the problem has been split into some small blocks in such a way that there are aerodynamic nodes belonging to more than one block. Copyright © 2014 John Wiley & Sons, Ltd.