A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

Mesh adaptation for simulation of unsteady flow with moving immersed boundaries

In this work, an approach for performing mesh adaptation in the numerical simulation of two‐dimensional unsteady flow with moving immersed boundaries is presented. In each adaptation period, the mesh is refined in the regions where the solution evolves or the moving bodies pass and is unrefined in the regions where the phenomena or the bodies deviate. The flow field and the fluid–solid interface are recomputed on the adapted mesh. The adaptation indicator is defined according to the magnitude of the vorticity in the flow field. There is no lag between the adapted mesh and the computed solution, and the adaptation frequency can be controlled to reduce the errors due to the solution transferring between the old mesh and the new one. The preservation of conservation property is mandatory in long‐time scale simulations, so a P1‐conservative interpolation is used in the solution transferring. A nonboundary‐conforming method is employed to solve the flow equations. Therefore, the moving‐boundary flows can be simulated on a fixed mesh, and there is no need to update the mesh at each time step to follow the motion or the deformation of the solid boundary. To validate the present mesh adaptation method, we have simulated several unsteady flows over a circular cylinder stationary or with forced oscillation, a single self‐propelled swimming fish, and two fish swimming in the same or different directions. Copyright © 2012 John Wiley & Sons, Ltd.     

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 local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical analysis on the propulsive performance and vortex shedding of fish‐like travelling wavy plate

Numerical analysis is carried out to investigate viscous flow over a travelling wavy plate undergoing lateral motion in the form of a streamwise travelling wave, which is similar to the backbone undulation of swimming fish. The two‐dimensional incompressible Navier–Stokes equations are solved using the finite element technique with the deforming‐spatial‐domain/stabilized space–time formulation. The objective of this study is to elucidate hydrodynamic features of flow structure and vortex shedding near the travelling wavy plate and to get into physical insights to the understanding of fish‐like swimming mechanisms in terms of drag reduction and optimal propulsive performance. The effects of some typical parameters, including the phase speed, amplitude, and relative wavelength of travelling wavy plate, on the flow structures, the forces, and the power consumption required for the propulsive motion of the plate are analysed. These results predicted by the present numerical analysis are well consistent with the available data obtained for the wave‐like swimming motion of live fish in nature. 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.     

A hybrid immersed‐boundary and multi‐block lattice Boltzmann method for simulating fluid and moving‐boundaries interactions

A numerical method is developed for modelling the interactions between incompressible viscous fluid and moving boundaries. The principle of this method is introducing the immersed‐boundary concept in the framework of the lattice Boltzmann method, and improving the accuracy and efficiency of the simulation by refining the mesh near moving boundaries. Besides elastic boundary with a constitutive law, the method can also efficiently simulate solid moving‐boundary interacting with fluid by employing the direct forcing technique. The method is validated by the simulations of flow past a circular cylinder, two cylinders moving with respect to each other and flow around a hovering wing. The versatility of the method is demonstrated by the numerical studies including elastic filament flapping in the wake of a cylinder and fish‐like bodies swimming in quiescent fluid. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of flow‐flexible body interactions with large deformation

A modified front‐tracking method was proposed for the simulation of fluid‐flexible body interactions with large deformations. A large deformable body was modeled by restructuring the body using a grid adaptation. Discontinuities in the viscosity at the fluid‐structure interface were incorporated by distributing the viscosity across the interface using an indicator function. A viscosity gradient field was created near the interface, and a smooth transition occurred between the structure and the fluid. The fluid motion was defined on the Eulerian domain and was solved using the fractional step method on a staggered Cartesian grid system. The solid motion was described by Lagrangian variables and was solved by the finite element method on an unstructured triangular mesh. The fluid motion and the structure motion were independently solved, and their interaction force was calculated using a feedback law. The interaction force was the restoring force of a stiff spring with damping, and spread from the Lagrangian coordinates to the Eulerian grid by a smoothed approximation of the Dirac delta function. In the numerical simulations, we validated the effect of the grid adaptation on the solid solver using a vibrating circular ring. The effects of the viscosity gradient field were verified by solving the deformation of a circular disk in a linear shear flow, including an elastic ring moving through a channel with constriction, deformation of a suspended catenary, and a swimming jellyfish. A comparison of the numerical results with the theoretical solutions was presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐order upwind control‐volume method based on integrated RBFs for fluid‐flow problems

This paper is concerned with the development of a high‐order upwind conservative discretization method for the simulation of flows of a Newtonian fluid in two dimensions. The fluid‐flow domain is discretized using a Cartesian grid from which non‐overlapping rectangular control volumes are formed. Line integrals arising from the integration of the diffusion and convection terms over control volumes are evaluated using the middle‐point rule. One‐dimensional integrated radial basis function schemes using the multiquadric basis function are employed to represent the variations of the field variables along the grid lines. The convection term is effectively treated using an upwind scheme with the deferred‐correction strategy. Several highly non‐linear test problems governed by the Burgers and the Navier–Stokes equations are simulated, which show that the proposed technique is stable, accurate and converges well. Copyright © 2010 John Wiley & Sons, Ltd.     

An adaptive Lagrangian boundary element approach for three‐dimensional transient free‐surface Stokes flow as applied to extrusion,thermoforming, and rheometry

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional simulation of confined free‐surface Stokes flow. The method is stable as it includes remeshing capabilities of the deforming free surface and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. Several flow problems are presented to illustrate the utility of the approach, particularly as encountered in polymer processing and rheology. These problems illustrate the transient nature of the flow during the processes of extrusion and thermoforming, the elongation of a fluid sample in an extensional rheometer, and the coating of a sphere. Surface tension effects are also explored. Copyright © 2001 John Wiley & Sons, Ltd.     

Thermo‐hydrodynamic‐mechanical multiphase flow model for CO2 sequestration in fracturing porous media

In this paper, we introduce a fully coupled thermo‐hydrodynamic‐mechanical computational model for multiphase flow in a deformable porous solid, exhibiting crack propagation due to fluid dynamics, with focus on CO₂ geosequestration. The geometry is described by a matrix domain, a fracture domain, and a matrix‐fracture domain. The fluid flow in the matrix domain is governed by Darcy's law and that in the crack is governed by the Navier–Stokes equations. At the matrix‐fracture domain, the fluid flow is governed by a leakage term derived from Darcy's law. Upon crack propagation, the conservation of mass and energy of the crack fluid is constrained by the isentropic process. We utilize the representative elementary volume‐averaging theory to formulate the mathematical model of the porous matrix, and the drift flux model to formulate the fluid dynamics in the fracture. The numerical solution is conducted using a mixed finite element discretization scheme. The standard Galerkin finite element method is utilized to discretize the diffusive dominant field equations, and the extended finite element method is utilized to discretize the crack propagation, and the fluid leakage at the boundaries between layers of different physical properties. A numerical example is given to demonstrate the computational capability of the model. It shows that the model, despite the relatively large number of degrees of freedom of different physical nature per node, is computationally efficient, and geometry and effectively mesh independent. Copyright © 2017 John Wiley & Sons, Ltd.     

An ALE approach to mechano‐chemical processes in fluid–structure interactions

Mathematical modeling and simulation of fluid–structure interaction problems are in the focus of research already for a longer period. However, taking into account also chemical reactions, leading to structural changes, including changes of mechanical properties of the solid phase, is rather new but for many applications is highly important area. This paper formulates a model system for reactive flow and transport in a vessel system, the penetration of chemical substances into the solid wall. Inside the wall, reactions take place that lead to changes of volume and of the mechanical properties of the wall. Numerical algorithms are developed and used to simulate the dynamics of such a mechano‐chemical fluid–structure interaction problem. As a proof of concept scenario, plaque formation in blood vessels is chosen. The arbitrary Lagrangian Eulerian approach (ALE) is chosen to solve the systems numerically. Temporal discretization of the fully coupled monolithic model is accomplished by backward Euler scheme and spatial discretization by stabilized finite elements. The numerical approach is verified by numerical tests, and effective methods to maintain mesh qualities under large deformations are described. For realistic system parameters, the simulations show that the plaque formation in blood vessel is a long‐time effect. The time scale of the formation is in the simulation of comparable order as in reality. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

Improving the CBS‐based partitioned semi‐implicit coupling algorithm for fluid‐structure interaction

We detail in this work 2 simple but effective alternatives to improve the characteristic‐based split–based partitioned semi‐implicit coupling algorithm for fluid‐structure interaction. The basic idea lies in introducing the end‐of‐step velocity into the implicit stages of the 2 algorithms integrating different splits. The algorithm built upon the second‐order pressure split is further stabilized via the pressure gradient projection with particular emphasis on the extremely low mass ratio. The smoothed finite element method is exploited for spatial discretization of fluid and solid equations. Even without any accelerators, both the semi‐implicit solvers incorporating fixed‐point iterations engender visible improvements versus the previously published data for several benchmarks.     

A σ‐coordinate three‐dimensional numerical model for surface wave propagation

A three‐dimensional numerical model based on the full Navier–Stokes equations (NSE) in σ‐coordinate is developed in this study. The σ‐coordinate transformation is first introduced to map the irregular physical domain with the wavy free surface and uneven bottom to the regular computational domain with the shape of a rectangular prism. Using the chain rule of partial differentiation, a new set of governing equations is derived in the σ‐coordinate from the original NSE defined in the Cartesian coordinate. The operator splitting method (Li and Yu, Int. J. Num. Meth. Fluids 1996; 23 : 485–501), which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified NSE. The model is first tested for mass and energy conservation as well as mesh convergence by using an example of water sloshing in a confined tank. Excellent agreements between numerical results and analytical solutions are obtained. The model is then used to simulate two‐ and three‐dimensional solitary waves propagating in constant depth. Very good agreements between numerical results and analytical solutions are obtained for both free surface displacements and velocities. Finally, a more realistic case of periodic wave train passing through a submerged breakwater is simulated. Comparisons between numerical results and experimental data are promising. The model is proven to be an accurate tool for consequent studies of wave‐structure interaction. Copyright © 2002 John Wiley & Sons, Ltd.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

A compressible two‐fluid multiphase model for CO2 leakage through a wellbore

This paper introduces an effectively mesh‐independent and computationally efficient model for CO₂ leakage through wellbores. A one‐dimensional compressible two‐fluid domain, representing a homogeneous air gas and a multiphase CO₂ with a jump at the interface between them, is modeled. The physical domain is modeled using the drift‐flux model, and the governing equations are solved using a mixed finite‐element discretization scheme. The standard Galerkin FEM, the partition of unity method, and the level‐set method are integrated to solve the problem. All important physical phenomena and processes occurring along the wellbore path, including fluid dynamics, buoyancy, phase change, compressibility, thermal interaction, wall friction, and slip between phases, together with the jump in density and enthalpy between air and CO₂, are considered. Two numerical examples illustrating the computational capability and efficiency of the model are presented. Copyright © 2015 John Wiley & Sons, Ltd.