A Lagrangian finite element approach for the analysis of fluid–structure interaction problems

A Lagrangian finite element method for the analysis of incompressible Newtonian fluid flows, based on a continuous re‐triangulation of the domain in the spirit of the so‐called Particle Finite Element Method, is here revisited and applied to the analysis of the fluid phase in fluid–structure interaction problems. A new approach for the tracking of the interfaces between fluids and structures is proposed. Special attention is devoted to the mass conservation problem. It is shown that, despite its Lagrangian nature, the proposed combined finite element‐particle method is well suited for large deformation fluid–structure interaction problems with evolving free surfaces and breaking waves. The method is validated against the available analytical and numerical benchmarks. Copyright © 2010 John Wiley & Sons, Ltd.     

Unified fractional step method for Lagrangian analysis of quasi‐incompressible fluid and nonlinear structure interaction using the PFEM

The fractional step method (FSM) is an efficient solution technique for the particle finite element method, a Lagrangian‐based approach to simulate fluid–structure interaction (FSI). Despite various refinements, the applicability of the FSM has been limited to low viscosity flow and FSI simulations with a small number of equations along the fluid–structure interface. To overcome these limitations, while incorporating nonlinear response in the structural domain, an FSM that unifies structural and fluid response in the discrete governing equations is developed using the quasi‐incompressible formulation. With this approach, fluid and structural particles do not need to be treated separately, and both domains are unified in the same system of equations. Thus, the equations along the fluid–structure interface do not need to be segregated from the fluid and structural domains. Numerical examples compare the unified FSM with the non‐unified FSM and show that the computational cost of the proposed method overcomes the slow convergence of the non‐unified FSM for high values of viscosity. As opposed to the non‐unified FSM, the number of iterations required for convergence with the unified FSM becomes independent of viscosity and time step, and the simulation run time does not depend on the size of the FSI interface. Copyright © 2016 John Wiley & Sons, Ltd.     

A coupled BEM and arbitrary Lagrangian–Eulerian FEM model for the solution of two‐dimensional laminar flows in external flow fields

This paper describes a new computational model developed to solve two‐dimensional incompressible viscous flow problems in external flow fields. The model based on the Navier–Stokes equations in primitive variables is able to solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the pressure projection method. The external flow field is simulated using the boundary element method by solving a pressure Poisson equation that assumes the pressure as zero at the infinite boundary. The momentum equation of the flow motion is solved using the three‐step finite element method. The arbitrary Lagrangian–Eulerian method is incorporated into the model, to solve the moving boundary problems. The present model is applied to simulate various external flow problems like flow across circular cylinder, acceleration and deceleration of the circular cylinder moving in a still fluid and vibration of the circular cylinder induced by the vortex shedding. The simulation results are found to be very reasonable and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Updated Lagrangian/Arbitrary Lagrangian–Eulerian framework for interaction between a compressible neo‐Hookean structure and an incompressible fluid

We propose a numerical method for a fluid–structure interaction problem. The material of the structure is homogeneous, isotropic, and it can be described by the compressible neo‐Hookean constitutive equation, while the fluid is governed by the Navier–Stokes equations. Our study does not use turbulence model. Updated Lagrangian method is used for the structure and fluid equations are written in Arbitrary Lagrangian–Eulerian coordinates. One global moving mesh is employed for the fluid–structure domain, where the fluid–structure interface is an ‘interior boundary’ of the global mesh. At each time step, we solve a monolithic system of unknown velocity and pressure defined on the global mesh. The continuity of velocity at the interface is automatically satisfied, while the continuity of stress does not appear explicitly in the monolithic fluid–structure system. This method is very fast because at each time step, we solve only one linear system. This linear system was obtained by the linearization of the structure around the previous position in the updated Lagrangian formulation and by the employment of a linear convection term for the fluid. Numerical results are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

Combined three‐dimensional lattice Boltzmann method and discrete element method for modelling fluid–particle interactions with experimental assessment

A general algorithmic framework is established in this paper for numerical simulations of three‐dimensional fluid–particle interaction problems with a large number of moving particles in turbulent flows using a combined lattice Boltzmann method (LBM) and discrete element method (DEM). In this approach, the fluid field is solved by the extended three‐dimensional LBM with the incorporation of the Smagorinsky turbulence model, while particle interactions are modelled by the DEM. The hydrodynamic interactions between fluid and particles are realized through the extension of an existing two‐dimensional fluid–particle hydrodynamic interaction scheme. The main computational aspects comprise the lattice Boltzmann formulation for the solution of fluid flows, the incorporation of a large eddy simulation‐based turbulence model within the framework of the three‐dimensional LBM for turbulent flows, the moving boundary condition for hydrodynamic interactions between fluid and moving particles, and the discrete element modelling of particle‐particle interactions. To assess the solution accuracy of the proposed approach, a much simplified laboratory model of vacuum dredging systems for mineral recovery is employed. The numerical results are compared with the experimental data available. It shows that the overall correspondence between numerical results and experimental measurements is good and thus indicates, to a certain extent, the solution accuracy of the proposed methodology. Copyright © 2009 John Wiley & Sons, Ltd.     

Velocity‐based formulations for standard and quasi‐incompressible hypoelastic‐plastic solids

We present three velocity‐based updated Lagrangian formulations for standard and quasi‐incompressible hypoelastic‐plastic solids. Three low‐order finite elements are derived and tested for non‐linear solid mechanics problems. The so‐called V‐element is based on a standard velocity approach, while a mixed velocity–pressure formulation is used for the VP and the VPS elements. The two‐field problem is solved via a two‐step Gauss–Seidel partitioned iterative scheme. First, the momentum equations are solved in terms of velocity increments, as for the V‐element. Then, the constitutive relation for the pressure is solved using the updated velocities obtained at the previous step. For the VPS‐element, the formulation is stabilized using the finite calculus method in order to solve problems involving quasi‐incompressible materials. All the solid elements are validated by solving two‐dimensional and three‐dimensional benchmark problems in statics as in dynamics. Copyright © 2016 John Wiley & Sons, Ltd.     

A coupled SDPH–FVM method for gas‐particle multiphase flow: Methodology

This paper describes a novel methodology that combines smoothed discrete particle hydrodynamics (SDPH) and finite volume method (FVM) to enhance the effective performance in solving the problems of gas‐particle multiphase flow. To describe the collision and fluctuation of particles, this method also increases a new parameter, namely, granular temperature, according to the kinetic theory of granular flow. The coupled framework of SDPH–FVM has been established, in which the drag force and pressure gradient act on the SDPH particles and the momentum sources of drag force are added back onto the FVM mesh. The proposed technique is a coupled discrete‐continuum method based on the two‐fluid model. To compute for the discrete phase, its SDPH is developed from smoothed particle hydrodynamics (SPH), in which the properties of SPH are redefined with some new physical quantities added into the traditional SPH parameters, so that it is more beneficial for SDPH in representing the particle characteristics. For the continuum phase, FVM is employed to discretize the continuum flow field on a stationary grid by capturing fluid characteristics. The coupled method exhibits strong efficiency and accuracy in several two‐dimensional numerical simulations. Copyright © 2016 John Wiley & Sons, Ltd.     

A particle method for two‐phase flows with large density difference

A new numerical approach for solving incompressible two‐phase flows is presented in the framework of the recently developed Consistent Particle Method (CPM). In the context of the Lagrangian particle formulation, the CPM computes spatial derivatives based on the generalized finite difference scheme and produces good results for single‐phase flow problems. Nevertheless, for two‐phase flows, the method cannot be directly applied near the fluid interface because of the abrupt discontinuity of fluid density resulting in large change in pressure gradient. This problem is resolved by dealing with the pressure gradient normalized by density, leading to a two‐phase CPM of which the original singlephase CPM is a special case. In addition, a new adaptive particle selection scheme is proposed to overcome the problem of ill‐conditioned coefficient matrix of pressure Poisson equation when particles are sparse and non‐uniformly spaced. Numerical examples of Rayleigh–Taylor instability, gravity current flow, water‐air sloshing and dam break are presented to demonstrate the accuracy of the proposed method in wave profile and pressure solution. Copyright © 2015 John Wiley & Sons, Ltd.     

A saturated discrete particle model and characteristic-based SPH method in granular materials

Based on the discrete particle model for solid-phase deformation of granular materials consisting of dry particulate assemblages, a discrete particle–continuum model for modelling the coupled hydro-mechanical behaviour in saturated granular materials is developed. The motion of the interstitial fluid is described by two parallel continuum schemes governed by the averaged incompressible N–S equations and Darcy's law, respectively, where the latter one can be regarded as a degraded case of the former. Owing to the merits in both Lagrangian and mesh-free characters, the characteristic-based smoothed particle hydrodynamics (SPH) method is proposed in this paper for modelling pore fluid flows relative to the deformed solid phase that is modelled as packed assemblages of interacting discrete particles. It is assumed that the formulation is Lagrangian with the co-ordinate system transferring with the movement of the solid particles. The assumed continuous fluid field is discretized into a finite set of Lagrangian (material) points with their number equal to that of solid particles situated in the computational domain. An explicit meshless scheme for granular materials with interstitial water is formulated. Numerical results illustrate the capability and performance of the present model in modelling the fluid–solid interaction and deformation in granular materials saturated with water. Copyright © 2007 John Wiley & Sons, Ltd.     

A PARTICLE TRACKING TECHNIQUE FOR THE LAGRANGIAN–EULERIAN FINITE ELEMENT METHOD IN MULTI-DIMENSIONS

This paper presents a multi-dimensional particle tracking technique for applying the Lagrangian–Eulerian finite element method to solve transport equations in transient-state simulations. In the Lagrangian– Eulerian approach, the advection term is handled in the Lagrangian step so that the associated numerical errors can be considerably reduced. It is important to have an adequate particle tracking technique for computing advection accurately in the Lagrangian step. The particle tracking technique presented here is designed to trace fictitious particles in the real-world flow field where the flow velocity is either measured or computed at a limited number of locations. The technique, named ‘in-element’ particle tracking, traces fictitious particles on an element-by-element basis. Given a velocity field, a fictitious particle is traced one element by one element until either a boundary is encountered or the available time is completely consumed. For the tracking within an element, the element is divided into a desired number of subelements with the interpolated velocity computed at all nodes of the subelements. A fictitious particle, thus, is traced one subelement by one subelement within the element. The desired number of subelements can be determined based on the complexity of the flow field being considered. The more complicated the flow field is, the more subelements are needed to achieve accurate particle tracking results. A single-velocity approach can be used to efficiently perform particle tracking in a smooth flow field, while an average-velocity approach can be employed to increase the tracking accuracy for more complex flow fields.     

A multimass correction for multicomponent fluid flow simulation using smoothed particle hydrodynamics

In multicomponent fluid flow simulations using smoothed particle hydrodynamics, the Lagrangian particles used are mostly of equal mass. This is preferred over multimass particle setup (particles with different values of mass), as it resolves the fluid interfaces comparatively better. But the flip side of using uniform mass particle setup is that it may not be computationally economical in situations with large‐density ratios. Hence, using multimass particle setup is both economical and perhaps inevitable. An attractive feature of multimass particle setup is that it allows uniform resolution in regions with different values of density. To take advantage of the multimass setup, it is therefore imperative to reduce the error associated with its usage. In this work, we present suitable multimass correction terms and assess its effectiveness using the ?h–smooth particle hydrodynamics scheme. Standard benchmark problems, viz, shock tube test, triple‐point shock test, Rayleigh‐Taylor instability, and Kelvin‐Helmholtz instability were solved with multimass particle setup, where significant improvements could be achieved in resolving the associated contact discontinuities.     

Source terms in Eulerian–Lagrangian contaminant transport simulation

Standard Eulerian treatment of source terms in Eulerian–Lagrangian numerical simulations results in poor performance at higher Courant numbers. To regain the customary high accuracy of Eulerian–Lagrangian methods under these conditions, a Lagrangian treatment of source terms is needed. It is also important to include the effects of fluid sources as well as contaminant sources. A new Lagrangian source formulation is presented, which has been implemented in a finite element simulator for contaminant transport in rivers and estuaries. Test problems demonstrate the high accuracy of the technique under a range of conditions, and its applicability to general multi‐dimensional problems and unstructured grids. Copyright © 2004 John Wiley & Sons, Ltd.     

Non-linear dynamic analysis of axisymmetric shells

Incremental equations of motion are derived from a Lagrangian variational formulation for the large displacement elastic-plastic and elastic-viscoplastic dynamic analysis of deformable bodies. The material constitutive behaviour is described in terms of the symmetric Piola–Kirchhoff stress and Lagrangian strain tensors. Degenerate isoparametric elements, permitting relaxation of the Kirchhoff–Love hypothesis, are used in a finite element formulation specialized for the analysis of shells of revolution subjected to axisymmetric loading. The linearized incremental equations of motion are solved using direct integration procedures, with added accuracy obtained from application of equilibrium correction at each step. The effectiveness of the numerical techniques is illustrated by the dynamic response analyses carried out on a shallow spherical cap subjected to uniform external step pressure loadings.     

A Jacobian‐free‐based IIM for incompressible flows involving moving interfaces with Dirichlet boundary conditions

In this paper, a finite difference marker‐and‐cell (MAC) scheme is presented for the steady Stokes equations with moving interfaces and Dirichlet boundary condition. The moving interfaces are represented by Lagrangian control points and their position is updated implicitly using a Jacobian‐free approach within each time step. The forces at the moving interfaces are calculated from the position of the interfaces and interpolated using cubic splines and then applied to the fluid through the related jump conditions. The proposed Jacobian‐free Newton–generalized minimum residual (GMRES) method avoids the need to form and store the matrix explicitly in the computation of the inverse of the Jacobian and betters numerical stability. The Stokes equations are discretized on a MAC grid via a second‐order finite difference scheme with the incorporation of jump contributions and the resulting saddle point system is solved by the conjugate gradient Uzawa‐type method. Numerical results demonstrate very well the accuracy and effectiveness of the proposed method. The present algorithm has been applied to solve incompressible Navier–Stokes flows with moving interfaces. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigation of fluid–structure interactions using a velocity‐linked P2/P1 finite element method and the generalized‐ α method

A velocity‐linked algorithm for solving unsteady fluid–structure interaction (FSI) problems in a fully coupled manner is developed using the arbitrary Lagrangian–Eulerian method. The P2/P1 finite element is used to spatially discretize the incompressible Navier–Stokes equations and structural equations, and the generalized‐ α method is adopted for temporal discretization. Common velocity variables are employed at the fluid–structure interface for the strong coupling of both equations. Because of the velocity‐linked formulation, kinematic compatibility is automatically satisfied and forcing terms do not need to be calculated explicitly. Both the numerical stability and the convergence characteristics of an iterative solver for the coupled algorithm are investigated by solving the FSI problem of flexible tube flows. It is noteworthy that the generalized‐ α method with small damping is free from unstable velocity fields. However, the convergence characteristics of the coupled system deteriorate greatly for certain Poisson's ratios so that direct solvers are essential for these cases. Furthermore, the proposed method is shown to clearly display the advantage of considering FSI in the simulation of flexible tube flows, while enabling much larger time‐steps than those adopted in some previous studies. This is possible through the strong coupling of the fluid and structural equations by employing common primitive variables. Copyright © 2012 John Wiley & Sons, Ltd.     

Formulation of equations of motion of finite element form for vehicle–track–bridge interaction system with two types of vehicle model

Vehicle, track and bridge are considered as an entire system in this paper. Two types of vertical vehicle model are described. One is a one foot mass–spring–damper system having two‐degree‐of‐freedom, and the other is four‐wheelset mass–spring–damper system with two‐layer suspension systems possessing 10‐degree‐of‐freedom. For the latter vehicle model, the eccentric load of car body is taken into account. The rails and the bridge deck are modelled as an elastic Bernoulli–Euler upper beam with finite length and a simply supported Bernoulli–Euler lower beam, respectively, while the elasticity and damping properties of the rail bed are represented by continuous springs and dampers. The dynamic contact forces between the moving vehicle and rails are considered as internal forces, so it is not necessary to calculate the internal forces for setting up the equations of motion of the vehicle–track–bridge interaction system. The two types of equations of motion of finite element form for the entire system are derived by means of the principle of a stationary value of total potential energy of dynamic system. The proposed method can set up directly the equations of motion for sophisticated system, and these equations can be solved by step‐by‐step integration method, to obtain simultaneously the dynamic responses of vehicle, of track and of bridge. Illustration examples are given. Copyright 2004 © John Wiley & Sons, Ltd.     

A particle‐in‐cell solution to the silo discharging problem

The problem of flow of a granular material during the process of discharging a silo is considered in the present paper. The mechanical behaviour of the material is described by the use of the model of the elastic–plastic solid with the Drucker–Prager yield condition and the non‐associative flow rule. The phenomenon of friction between the stored material and the silo walls is taken into account—the Coulomb model of friction is used in the analysis. The problem is analysed by means of the particle‐in‐cell method—a variant of the finite element method which enables to solve the pertinent equations of motion on an arbitrary computational mesh and trace state variables at points of the body chosen independently of the mesh. The method can be regarded as an arbitrary Lagrangian–Eulerian formulation of the finite element method, and overcomes the main drawback of the updated Lagrangian formulation of FEM related to mesh distortion. The entire process of discharging a silo can be analysed by this approach. The dynamic problem is solved by the use of the explicit time‐integration scheme. Several numerical examples are included. The plane strain and axisymmetric problems are solved for silos with flat bottoms and conical hoppers. Some results are compared with experimental ones. Copyright © 1999 John Wiley & Sons, Ltd.     

Explicit Reproducing Kernel Particle Methods for large deformation problems

The explicit Reproducing Kernel Particle Method (RKPM) is presented and applied to the simulations of large deformation problems. RKPM is a meshless method which does not need a mesh structure in its formulation. Because of this mesh-free property, RKPM is able to simulate large deformation problems without remeshing which is often required for the mesh-based methods such as the finite element method. The RKPM shape function and its derivatives are constructed by imposing the consistency conditions. An efficient treatment of essential boundary conditions is also proposed for explicit time integration. The Lagrangian method based on the reference configuration is employed for the RKPM simulation of large deformation problems. Several examples of non-linear elastic materials are solved to demonstrate the performance of the method. The numerical experiment for the problem of underwater bubble explosion is also performed using the explicit Lagrangian RKPM formulation. © 1998 John Wiley & Sons, Ltd.