Parallel‐in‐time simulation of the unsteady Navier–Stokes equations for incompressible flow

In this paper, we present an application of a parallel‐in‐time algorithm for the solution of the unsteady Navier–Stokes model equations that are of parabolic–elliptic type. This method is based on the alternated use of a coarse global sequential solver and a fine local parallel one. A standard finite volume/finite differences first‐order approach is used for discretization of the unsteady two‐dimensional Navier–Stokes equations. The Taylor vortex decay problem and the confined flow around a square cylinder were selected as unsteady flow examples to illustrate and analyse the properties of the parallel‐in‐time method through numerical experiments. The influence of several parameters on the computing time required to perform a parallel‐in‐time calculation on a PC cluster was verified. Among them we have analysed the influence of the number of processors, the number of iterations for convergence, the resolution of the spatial domain and the influence of the time‐step sizes ratio between the coarse and fine grids. Significant computer time saving was achieved when compared with the single processor computing time, particularly when the spatial dimension of the problem is low and the temporal scale is large. Copyright © 2004 John Wiley & Sons, Ltd.     

Two‐dimensional implicit time‐dependent calculations on adaptive unstructured meshes with time evolving boundaries

An implicit multigrid‐driven algorithm for two‐dimensional incompressible laminar viscous flows has been coupled with a solution adaptation method and a mesh movement method for boundary movement. Time‐dependent calculations are performed implicitly by regarding each time step as a steady‐state problem in pseudo‐time. The method of artificial compressibility is used to solve the flow equations. The solution mesh adaptation method performs local mesh refinement using an incremental Delaunay algorithm and mesh coarsening by means of edge collapse. Mesh movement is achieved by modeling the computational domain as an elastic solid and solving the equilibrium equations for the stress field. The solution adaptation method has been validated by comparison with experimental results and other computational results for low Reynolds number flow over a shedding circular cylinder. Preliminary validation of the mesh movement method has been demonstrated by a comparison with experimental results of an oscillating airfoil and with computational results for an oscillating cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

Unstructured pressure‐correction solver based on a consistent discretization of the Poisson equation

A finite volume incompressible flow solver is presented for three‐dimensional unsteady flows based on an unstructured tetrahedral mesh, with collocation of the flow variables at the cell vertices. The solver is based on the pressure‐correction method, with an explicit prediction step of the momentum equations followed by a Poisson equation for the correction step to enforce continuity. A consistent discretization of the Poisson equation was found to be essential in obtaining a solution. The correction step was solved with the biconjugate gradient stabilized (Bi‐CGSTAB) algorithm coupled with incomplete lower–upper (ILU) preconditioning. Artificial dissipation is used to prevent the formation of instabilities. Flow solutions are presented for a stalling airfoil, vortex shedding past a bridge deck and flow in model alveoli. Copyright © 2000 John Wiley & Sons, Ltd.     

An unstructured finite volume technique for the 3D Poisson equation on arbitrary geometry using a σ‐coordinate system

We present a solver for a three‐dimensional Poisson equation issued from the Navier–Stokes equations applied to model rivers, estuaries, and coastal flows. The three‐dimensional physical domain is composed of an arbitrary domain in the horizontal direction and is bounded by an irregular free surface and bottom in the vertical direction. The equations are transformed vertically to the σ‐coordinate system to obtain an accurate representation of top and bottom topographies. The method is based on a second‐order finite volume technique on prisms consisting of triangular grids in the horizontal direction. The algorithm is accompanied by an analysis of different linear system solvers in order to achieve fast solutions. Numerical experiments are conducted to test the numerical accuracy and the computational efficiency of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Coupled lubrication and Stokes flow finite elements

A method is developed for performing a local reduction of the governing physics for fluid problems with domains that contain a combination of narrow and non‐narrow regions, and the computational accuracy and performance of the method are measured. In the narrow regions of the domain, where the fluid is assumed to have no inertia and the domain height and curvature are assumed small, lubrication, or Reynolds, theory is used locally to reduce the two‐dimensional Navier–Stokes equations to the one‐dimensional Reynolds equation while retaining a high degree of accuracy in the overall solution. The Reynolds equation is coupled to the governing momentum and mass equations of the non‐narrow region with boundary conditions on the mass and momentum flux. The localized reduction technique, termed ‘stitching,’ is demonstrated on Stokes flow for various geometries of the hydrodynamic journal bearing—a non‐trivial test problem for which a known analytical solution is available. The computational advantage of the coupled Stokes–Reynolds method is illustrated on an industrially applicable fully‐flooded deformable‐roll coating example. The examples in this paper are limited to two‐dimensional Stokes flow, but extension to three‐dimensional and Navier–Stokes flow is possible. Copyright © 2003 John Wiley & Sons, Ltd.     

A second look at the role of the fast Fourier transform as an elliptic solver

A fast cosine transform (FCT) is coupled with a tridiagonal solver for the purpose of solving the Poisson equation on irregular and non‐uniform rectangular staggered grids. This kind of solution is required for the pressure field during the simulation of the incompressible Navier–Stokes equations when using the projection method. A new technique using the FCT–tridiagonal solver is derived for the cases where the boundaries of the flow regime do not coincide with the boundaries of the computational domain and for non‐uniform grids. The technique is based on an iterative procedure where a defect equation is solved in every iteration, followed by a relaxation procedure. The method is investigated analytically and numerically to show that the solution converges as a geometric series. The method is further investigated for the effects of the relative size of the rigid body, the grid stretching, size and aspect ratio. The new solver is incorporated with the direct numerical simulation (DNS) and large eddy simulation (LES) techniques to simulate the flows around a backward‐facing step and a 3D rectangular obstacle, yielding results that qualitatively compare well with known results. Copyright © 2005 John Wiley & Sons, Ltd.     

Coupled Navier–Stokes—Molecular dynamics simulations using a multi‐physics flow simulation framework

Simulation of nano‐scale channel flows using a coupled Navier–Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi‐physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub‐domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier–Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier–Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro‐scale MD problems are constructed to provide viscous stresses for the Navier–Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite‐volume discretization are evaluated using MD. The second example deals with fully developed three‐dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier–Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite‐volume discretization of the macro‐domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro‐scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd.     

A matrix‐free preconditioned Newton/GMRES method for unsteady Navier–Stokes solutions

The unsteady compressible Reynolds‐averaged Navier–Stokes equations are discretized using the Osher approximate Riemann solver with fully implicit time stepping. The resulting non‐linear system at each time step is solved iteratively using a Newton/GMRES method. In the solution process, the Jacobian matrix–vector products are replaced by directional derivatives so that the evaluation and storage of the Jacobian matrix is removed from the procedure. An effective matrix‐free preconditioner is proposed to fully avoid matrix storage. Convergence rates, computational costs and computer memory requirements of the present method are compared with those of a matrix Newton/GMRES method, a four stage Runge–Kutta explicit method, and an approximate factorization sub‐iteration method. Effects of convergence tolerances for the GMRES linear solver on the convergence and the efficiency of the Newton iteration for the non‐linear system at each time step are analysed for both matrix‐free and matrix methods. Differences in the performance of the matrix‐free method for laminar and turbulent flows are highlighted and analysed. Unsteady turbulent Navier–Stokes solutions of pitching and combined translation–pitching aerofoil oscillations are presented for unsteady shock‐induced separation problems associated with the rotor blade flows of forward flying helicopters. Copyright © 2000 John Wiley & Sons, Ltd.     

3D free‐surface flow computation using a RANSE/Fourier–Kochin coupling

A coupling method for numerical calculations of steady free‐surface flows around a body is presented. The fluid domain in the neighbourhood of the hull is divided into two overlapping zones. Viscous effects are taken in account near the hull using Reynolds‐averaged Navier–Stokes equations (RANSE), whereas potential flow provides the flow away from the hull. In the internal domain, RANSE are solved by a fully coupled velocity, pressure and free‐surface elevation method. In the external domain, potential‐flow theory with linearized free‐surface condition is used to provide boundary conditions to the RANSE solver. The Fourier–Kochin method based on the Fourier–Kochin formulation, which defines the velocity field in a potential‐flow region in terms of the velocity distribution at a boundary surface, is used for that purpose. Moreover, the free‐surface Green function satisfying this linearized free‐surface condition is used. Calculations have been successfully performed for steady ship‐waves past a serie 60 and then have demonstrated abilities of the present coupling algorithm. Copyright © 2003 John Wiley & Sons, Ltd.     

A modular,operator‐splitting scheme for fluid–structure interaction problems with thick structures

We present an operator‐splitting scheme for fluid–structure interaction (FSI) problems in hemodynamics, where the thickness of the structural wall is comparable to the radius of the cylindrical fluid domain. The equations of linear elasticity are used to model the structure, while the Navier–Stokes equations for an incompressible viscous fluid are used to model the fluid. The operator‐splitting scheme, based on the Lie splitting, separates the elastodynamics structure problem from a fluid problem in which structure inertia is included to achieve unconditional stability. We prove energy estimates associated with unconditional stability of this modular scheme for the full nonlinear FSI problem defined on a moving domain, without requiring any sub‐iterations within time steps. Two numerical examples are presented, showing excellent agreement with the results of monolithic schemes. First‐order convergence in time is shown numerically. Modularity, unconditional stability without temporal sub‐iterations, and simple implementation are the features that make this operator‐splitting scheme particularly appealing for multi‐physics problems involving FSI. Copyright © 2013 John Wiley & Sons, Ltd.     

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 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.     

Fast boundary‐domain integral algorithm for the computation of incompressible fluid flow problems

The paper deals with the numerical solution of fluid dynamics using the boundary‐domain integral method (BDIM). A velocity–vorticity formulation of the Navier–Stokes equations is adopted, where the kinematic equation is written in its parabolic form. Computational aspects of the numerical simulation of two‐dimensional flows is described in detail. In order to lower the computational cost, the subdomain technique is applied. A preconditioned Krylov subspace method (PKSM) is used for the solution of systems of linear equations. Level‐based fill‐in incomplete lower upper decomposition (ILU) preconditioners are developed and their performance is examined. Scaling of stopping criteria is applied to minimize the number of iterations for the PKSM. The effectiveness of the proposed method is tested on several benchmark test problems. Copyright © 1999 John Wiley & Sons, Ltd.     

Implicit application of non‐reflective boundary conditions for Navier–Stokes equations in generalized coordinates

The non‐reflective boundary conditions (NRBC) for Navier–Stokes equations originally suggested by Poinsot and Lele (J. Comput. Phys. 1992; 101 :104–129) in Cartesian coordinates are extended to generalized coordinates. The characteristic form Navier–Stokes equations in conservative variables are given. In this characteristic‐based method, the NRBC is implicitly coupled with the Navier–Stokes flow solver and are solved simultaneously with the flow solver. The calculations are conducted for a subsonic vortex propagating flow and the steady and unsteady transonic inlet‐diffuser flows. The results indicate that the present method is accurate and robust, and the NRBC are essential for unsteady flow calculations. Copyright © 2005 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.     

Nonlinear acceleration of coupled fluid–structure transient thermal problems by Anderson mixing

Conjugate heat‐transfer problems are typically solved using partitioned methods where fluid and solid subdomains are evaluated separately by dedicated solvers coupled through a common boundary. Strongly coupled schemes for transient analysis require fluid and solid problems to be solved many times each time step until convergence to a steady state. In many practical situations, a fairly simple and frequently employed fixed‐point iteration process is rather ineffective; it leads to a large number of iterations per time step and consequently to long simulation times. In this article, Anderson mixing is proposed as a fixed‐point convergence acceleration technique to reduce computational cost of thermal coupled fluid–solid problems. A number of other recently published methods with applications to similar fluid–structure interaction problems are also reviewed and analyzed. Numerical experiments are presented to illustrate relative performance of these methods on a test problem of rotating pre‐swirl cavity air flow interacting with a turbine disk. It is observed that performance of Anderson mixing method is superior to that of other algorithms in terms of total iteration counts. Additional computational savings are demonstrated by reusing information from previously solved time steps. Copyright © All rights reserved 2012 Rolls‐Royce plc.     

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of unsteady multidimensional free surface motions by level set method

This paper presents a numerical method that couples the incompressible Navier–Stokes equations with the level set method in a curvilinear co‐ordinate system for study of free surface flows. The finite volume method is used to discretize the governing equations on a non‐staggered grid with a four‐step fractional step method. The free surface flow problem is converted into a two‐phase flow system on a fixed grid in which the free surface is implicitly captured by the zero level set. We compare different numerical schemes for advection of the level set function in a generalized curvilinear format, including the third order quadratic upwind interpolation for convective kinematics (QUICK) scheme, and the second and third order essentially non‐oscillatory (ENO) schemes. The level set equations of evolution and reinitialization are validated with benchmark cases, e.g. a stationary circle, a rotating slotted disk and stretching of a circular fluid element. The coupled system is then applied to a travelling solitary wave, and two‐ and three‐dimensional dam breaking problems. Some interesting free surface phenomena are revealed by the computational results, such as, the large free surface vortices, air entrapment and splashing of the water surge front. The computational results are in excellent agreement with theoretical predictions and experimental data, where they are available. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

Finite element and implicit Runge–Kutta implementation of an acoustics–convection upstream resolution algorithm for the time‐dependent two‐dimensional Euler equations

The second of a two‐paper series, this paper details a solver for the characteristics‐bias system from the acoustics–convection upstream resolution algorithm for the Euler and Navier–Stokes equations. An integral formulation leads to several surface integrals that allow effective enforcement of boundary conditions. Also presented is a new multi‐dimensional procedure to enforce a pressure boundary condition at a subsonic outlet, a procedure that remains accurate and stable. A classical finite element Galerkin discretization of the integral formulation on any prescribed grid directly yields an optimal discretely conservative upstream approximation for the Euler and Navier–Stokes equations, an approximation that remains multi‐dimensional independently of the orientation of the reference axes and computational cells. The time‐dependent discrete equations are then integrated in time via an implicit Runge–Kutta procedure that in this paper is proven to remain absolutely non‐linearly stable for the spatially‐discrete Euler and Navier–Stokes equations and shown to converge rapidly to steady states, with maximum Courant number exceeding 100 for the linearized version. Even on relatively coarse grids, the acoustics–convection upstream resolution algorithm generates essentially non‐oscillatory solutions for subsonic, transonic and supersonic flows, encompassing oblique‐ and interacting‐shock fields that converge within 40 time steps and reflect reference exact solutions. Copyright © 2005 John Wiley & Sons, Ltd.