Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical solution of three‐dimensional velocity–vorticity Navier–Stokes equations by finite difference method

This paper describes the finite difference numerical procedure for solving velocity–vorticity form of the Navier–Stokes equations in three dimensions. The velocity Poisson equations are made parabolic using the false‐transient technique and are solved along with the vorticity transport equations. The parabolic velocity Poisson equations are advanced in time using the alternating direction implicit (ADI) procedure and are solved along with the continuity equation for velocities, thus ensuring a divergence‐free velocity field. The vorticity transport equations in conservative form are solved using the second‐order accurate Adams–Bashforth central difference scheme in order to assure divergence‐free vorticity field in three dimensions. The velocity and vorticity Cartesian components are discretized using a central difference scheme on a staggered grid for accuracy reasons. The application of the ADI procedure for the parabolic velocity Poisson equations along with the continuity equation results in diagonally dominant tri‐diagonal matrix equations. Thus the explicit method for the vorticity equations and the tri‐diagonal matrix algorithm for the Poisson equations combine to give a simplified numerical scheme for solving three‐dimensional problems, which otherwise requires enormous computational effort. For three‐dimensional‐driven cavity flow predictions, the present method is found to be efficient and accurate for the Reynolds number range 100?Re?2000. Copyright © 2004 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

3D two phase flows numerical simulations by SL‐VOF method

This paper describes the development of a semi‐Lagrangian computational method for simulating complex 3D two phase flows. The Navier–Stokes equations are solved separately in both fluids using a robust pseudo‐compressibility method able to deal with high density ratio. The interface tracking is achieved by the segment Lagrangian volume of fluid (SL‐VOF) method. The 2D SL‐VOF method using the concepts of VOF, piecewise linear interface calculation (PLIC) and Lagrangian advection of the interface is herein extended to 3D flows. Three different test cases of SL‐VOF 3D are presented for validation and comparison either with 2D flows or with other numerical methods. A good agreement is observed in each case. Copyright © 2004 John Wiley & Sons, Ltd.     

A fast method for solving fluid–structure interaction problems numerically

The paper presents a semi‐implicit algorithm for solving an unsteady fluid–structure interaction problem. The algorithm for solving numerically the fluid–structure interaction problems was obtained by combining the backward Euler scheme with a semi‐implicit treatment of the convection term for the Navier–Stokes equations and an implicit centered scheme for the structure equations. The structure is governed either by the linear elasticity or by the non‐linear St Venant–Kirchhoff elasticity models. At each time step, the position of the interface is predicted in an explicit way. Then, an optimization problem must be solved, such that the continuity of the velocity as well as the continuity of the stress hold at the interface. During the Broyden, Fletcher, Goldforb, Shano (BFGS) iterations for solving the optimization problem, the fluid mesh does not move, which reduces the computational effort. The term ‘semi‐implicit’ used for the fully algorithm means that the interface position is computed explicitly, while the displacement of the structure, velocity and the pressure of the fluid are computed implicitly. Numerical results are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.     

Vorticity–velocity formulation of the 3D Navier–Stokes equations in cylindrical co‐ordinates

A finite difference method is presented for solving the 3D Navier–Stokes equations in vorticity–velocity form. The method involves solving the vorticity transport equations in ‘curl‐form’ along with a set of Cauchy–Riemann type equations for the velocity. The equations are formulated in cylindrical co‐ordinates and discretized using a staggered grid arrangement. The discretized Cauchy–Riemann type equations are overdetermined and their solution is accomplished by employing a conjugate gradient method on the normal equations. The vorticity transport equations are solved in time using a semi‐implicit Crank–Nicolson/Adams–Bashforth scheme combined with a second‐order accurate spatial discretization scheme. Special emphasis is put on the treatment of the polar singularity. Numerical results of axisymmetric as well as non‐axisymmetric flows in a pipe and in a closed cylinder are presented. Comparison with measurements are carried out for the axisymmetric flow cases. Copyright © 2003 John Wiley & Sons, Ltd.     

Gas‐assisted fluid displacement in a circular tube and a rectangular channel

In this paper the amount of liquid left inside of a circular tube and a rectangular channel when displaced by another immiscible fluid are determined by solving the full creeping‐motion equations. The exact continuity of stress on the free surface is employed with a finite difference method. In order to solve the equations, the steady‐state shape of the interface is guessed and the normal stress boundary condition is dropped. The equations based on a stream function‐vorticity formulation are solved with the aid of elliptic grid generation. The computed results are compared with the experimental results of Taylor (J. Fluid Mech. 1961; 10: 161), the theoretical results of Reinelt and Saffman (SIAM J. Sci. Stat. Comput. 1985; 6: 542) and our experimental data. The computed results are in close agreement with our experimental data and those of previous investigators. Copyright © 2002 John Wiley & Sons, Ltd.     

Coupled fluid–solid interaction under shock wave loading

This paper considers the treatment of fluid–solid interaction problems under shock wave loading, where the solid experiences large bulk Lagrangian displacements. This work addresses the issues associated with using a level set as a generalized interface for fluid–solid coupling where the fluid–solid interface is embedded in an unstructured fluid grid. We outline the formulation used for the edge‐based unstructured‐grid Euler solver. The identification of the fluid–solid interface on the unstructured fluid mesh uses a super‐sampled ??² projection technique, which in conjunction with a Lagrangian interface position, permits fast identification of the interface and the concomitant imposition of boundary conditions. The use of a narrow‐band approach for the identification of the wetted interface is presented with the details of the construction of interface conditions. A series of two and three‐dimensional shock‐body computations are presented to demonstrate the validity of the current approach on problems with static and dynamic interfaces, including projectile/shock interaction simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐resolution method for compressible two‐fluid flows and simulation of three‐dimensional shock–bubble interactions

A high‐resolution method is developed to capture the material interfaces of compressible two‐fluid flows in multiple dimensions. A fluid mixture model system with single velocity and pressure is used, and viscous effect can also be taken into account. A consistent thermodynamic law based on the assumption of pressure equilibrium is employed to describe the thermodynamic behaviors of the pure fluids and mixture of two components. The splitting and unsplit Eulerian formulations of piecewise parabolic method are extended to numerically integrate the hyperbolic part of the model system, whereas the system of diffusion equations is solved using an explicit, central difference scheme. The block‐structured adaptive mesh refinement (AMR) capability is built in the hydrodynamic code to locally improve grid resolution. The resulting method is verified to be at least second‐order accurate in space. Numerical results show that the discontinuities, particularly contact discontinuities, can be resolved sharply. The use of AMR allows flow features at disparate scales to be resolved sufficiently. In addition, three‐dimensional shock–bubble interactions are simulated to investigate effects of Mach number on bubble evolution. The flow structures including those peculiar to three‐dimensional bubble are resolved correctly, and some physical phenomena with increasing Mach number are reported. Copyright © 2012 John Wiley & Sons, Ltd.     

Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods

A robust, accurate, and computationally efficient interface tracking algorithm is a key component of an embedded computational framework for the solution of fluid–structure interaction problems with complex and deformable geometries. To a large extent, the design of such an algorithm has focused on the case of a closed embedded interface and a Cartesian computational fluid dynamics grid. Here, two robust and efficient interface tracking computational algorithms capable of operating on structured as well as unstructured three‐dimensional computational fluid dynamics grids are presented. The first one is based on a projection approach, whereas the second one is based on a collision approach. The first algorithm is faster. However, it is restricted to closed interfaces and resolved enclosed volumes. The second algorithm is therefore slower. However, it can handle open shell surfaces and underresolved enclosed volumes. Both computational algorithms exploit the bounding box hierarchy technique and its parallel distributed implementation to efficiently store and retrieve the elements of the discretized embedded interface. They are illustrated, and their respective performances are assessed and contrasted, with the solution of three‐dimensional, nonlinear, dynamic fluid–structure interaction problems pertaining to aeroelastic and underwater implosion applications. Copyright © 2012 John Wiley & Sons, Ltd.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

Multigrid solution of the incompressible Navier–Stokes equations for three‐dimensional recirculating flow: coupled and decoupled smoothers compared

A comparison of multigrid methods for solving the incompressible Navier–Stokes equations in three dimensions is presented. The continuous equations are discretised on staggered grids using a second‐order monotonic scheme for the convective terms and implemented in defect correction form. The convergence characteristics of a decoupled method (SIMPLE) are compared with those of the cellwise coupled method (SCGS). The convergence rates obtained for computations of the three‐dimensional lid‐driven cavity problem are found to be very similar to those obtained for computations of the corresponding two‐dimensional problem with comparable grid density. Although the convergence rate of SCGS is thus superior to that of SIMPLE, the decoupled method is found to be more efficient computationally and requires less computing time for a given level of convergence. The linewise implementation of the coupled method (CLGS) is also investigated and shown to be more efficient than SCGS, although the convergence rate and computing time required per cycle are both found to depend on the direction of sweep. The optimal implementation of CLGS is found to be only marginally more effective than SIMPLE, but a change to the structure of the data storage would increase the advantage. Copyright © 1999 John Wiley & Sons, Ltd.     

A 3‐D non‐hydrostatic pressure model for small amplitude free surface flows

A three‐dimensional, non‐hydrostatic pressure, numerical model with k–ε equations for small amplitude free surface flows is presented. By decomposing the pressure into hydrostatic and non‐hydrostatic parts, the numerical model uses an integrated time step with two fractional steps. In the first fractional step the momentum equations are solved without the non‐hydrostatic pressure term, using Newton's method in conjunction with the generalized minimal residual (GMRES) method so that most terms can be solved implicitly. This method only needs the product of a Jacobian matrix and a vector rather than the Jacobian matrix itself, limiting the amount of storage and significantly decreasing the overall computational time required. In the second step the pressure–Poisson equation is solved iteratively with a preconditioned linear GMRES method. It is shown that preconditioning reduces the central processing unit (CPU) time dramatically. In order to prevent pressure oscillations which may arise in collocated grid arrangements, transformed velocities are defined at cell faces by interpolating velocities at grid nodes. After the new pressure field is obtained, the intermediate velocities, which are calculated from the previous fractional step, are updated. The newly developed model is verified against analytical solutions, published results, and experimental data, with excellent agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Two‐dimensional non‐Newtonian injection molding with a new control volume FEM/volume of fluid method

A new method based on volume of fluid for interface tracking in the simulation of injection molding is presented. The proposed method is comprised of two main stages: accumulation and distribution of the volume fraction. In the first stage the equation for the volume fraction with a noninterfacial flux condition is solved. In the second stage the accumulated volume of fluid that arises as a consequence of the application of the first one is dispersed. This procedure guarantees that the fluid fills the available space without dispersion of the interface. The mathematical model is based on two‐phase transport equations that are numerically integrated through the control volume finite element method. The numerical results for the interface position are successfully verified with analytical results and numerical data available in the literature for one‐dimensional and two‐dimensional domains. The transient position of the advance fronts showed an effective and consistent simulation of an injection molding process. The nondispersive volume of fluid method here proposed is implemented for the simulation of nonisothermal injection molding in two‐dimensional cavities. The obtained results are represented as transient interface positions, isotherms and pressure distributions during the injection molding of low density polyethylene. Copyright © 2012 John Wiley & Sons, Ltd.     

A note on two upwind strategies for RBF‐based grid‐free schemes to solve steady convection–diffusion equations

In this paper, two radial basis function (RBF)‐based local grid‐free upwind schemes have been discussed for convection–diffusion equations. The schemes have been validated over some convection–diffusion problems with sharp boundary layers. It is found that one of the upwind schemes realizes the boundary layers more accurately than the rest. Comparisons with the analytical solutions demonstrate that the local RBF grid‐free upwind schemes based on the exact velocity direction are stable and produce accurate results on domains discretized even with scattered distribution of nodal points. Copyright © 2009 John Wiley & Sons, Ltd.     

Two‐dimensional compact finite difference immersed boundary method

We present a compact finite differences method for the calculation of two‐dimensional viscous flows in biological fluid dynamics applications. This is achieved by using body‐forces that allow for the imposition of boundary conditions in an immersed moving boundary that does not coincide with the computational grid. The unsteady, incompressible Navier–Stokes equations are solved in a Cartesian staggered grid with fourth‐order Runge–Kutta temporal discretization and fourth‐order compact schemes for spatial discretization, used to achieve highly accurate calculations. Special attention is given to the interpolation schemes on the boundary of the immersed body. The accuracy of the immersed boundary solver is verified through grid convergence studies. Validation of the method is done by comparison with reference experimental results. In order to demonstrate the application of the method, 2D small insect hovering flight is calculated and compared with available experimental and computational results. Copyright © 2009 John Wiley & Sons, Ltd.