Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.     

Issue Information

An integrated finite element method (FEM) is proposed to simulate incompressible two‐phase flows with surface tension effects, and three different surface tension models are applied to the FEM to investigate spurious currents and temporal stability. A Q2Q1 element is adopted to solve the continuity and Navier–Stokes equations and a Q2‐iso‐Q1 to solve the level set equation. The integrated FEM solves pressure and velocity simultaneously in a strongly coupled manner; the level set function is reinitialized by adopting a direct approach using interfacial geometry information instead of solving a conventional hyperbolic‐type equation. In addition, a consistent continuum surface force (consistent CSF) model is utilized by employing the same basis function for both surface tension and pressure variables to damp out spurious currents and to estimate the accurate pressure distribution. The model is further represented as a semi‐implicit manner to improve temporal stability with an increased time step. In order to verify the accuracy and robustness of the code, the present method is applied to a few benchmark problems of the static bubble and rising bubble with large density and viscosity ratios. The Q2Q1‐integrated FEM coupled with the semi‐implicit consistent CSF demonstrates the significantly reduced spurious currents and improved temporal stability. The numerical results are in good qualitative and quantitative agreements with those of the existing studies. Copyright © 2015 John Wiley & Sons, Ltd.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Algebraic multigrid and incompressible fluid flow

This paper is concerned with the development of algebraic multigrid (AMG) solution methods for the coupled vector–scalar fields of incompressible fluid flow. It addresses in particular the problems of unstable smoothing and of maintaining good vector–scalar coupling in the AMG coarse‐grid approximations. Two different approaches have been adopted. The first is a direct approach based on a second‐order discrete‐difference formulation in primitive variables. Here smoothing is stabilized using a minimum residual control harness and velocity–pressure coupling is maintained by employing a special interpolation during the construction of the inter‐grid transfer operators. The second is an indirect approach that avoids the coupling problem altogether by using a fourth‐order discrete‐difference formulation in a single scalar‐field variable, primitive variables being recovered in post‐processing steps. In both approaches the discrete‐difference equations are for the steady‐state limit (infinite time step) with a fully implicit treatment of advection based on central differencing using uniform and non‐uniform unstructured meshes. They are solved by Picard iteration, the AMG solvers being used repeatedly for each linear approximation. Both classical AMG (C‐AMG) and smoothed‐aggregation AMG (SA‐AMG) are used. In the direct approach, the SA‐AMG solver (with inter‐grid transfer operators based on mixed‐order interpolation) provides an almost mesh‐independent convergence. In the indirect approach for uniform meshes, the C‐AMG solver (based on a Jacobi‐relaxed interpolation) provides solutions with an optimum scaling of the convergence rates. For non‐uniform meshes this convergence becomes mesh dependent but the overall solution cost increases relatively slowly with increasing mesh bandwidth. Copyright © 2006 John Wiley & Sons, Ltd.     

An attempt to improve accuracy of higher‐order statistics and spectra in direct numerical simulation of incompressible wall turbulence by using the compact schemes for viscous terms

We attempt to improve accuracy in the high‐wavenumber region in DNS of incompressible wall turbulence such as found in fully developed turbulent channel flow. In particular, it is shown that the improvement of accuracy of viscous terms in the Navier–Stokes equations leads to the improvement of accuracy of higher‐order statistics and various spectra. It is emphasized that increase in required computational cost will not be crucial when incompressible flow is simulated, because the introduction of a higher‐order scheme into the viscous terms does not increase computational cost for solving the Poisson equation. We introduced fourth‐order and eighth‐order central compact schemes for discretizing the viscous terms in DNS of a fully developed turbulent channel flow. The results are compared with those using second‐order and fourth‐order central‐difference schemes applied to the viscous terms and those obtained by the spectral method. The results show that accuracy improvement of the viscous terms improve accuracy of higher‐order statistics (i.e., skewness and flatness factors of streamwise velocity fluctuation) and various spectra of velocity and pressure fluctuations in the high‐wavenumber region. Copyright © 2013 John Wiley & Sons, Ltd.     

Incompressible Roe‐like scheme for turbulent flows

In the current study, numerical investigation of incompressible turbulent flow is presented. By the artificial compressibility method, momentum and continuity equations are coupled. Considering Reynolds averaged Navier–Stokes equations, the Spalart–Allmaras turbulence model, which has accurate results in two‐dimensional problems, is used to calculate Reynolds stresses. For convective fluxes a Roe‐like scheme is proposed for the steady Reynolds averaged Navier–Stokes equations. Also, Jameson averaging method was implemented. In comparison, the proposed characteristics‐based upwind incompressible turbulent Roe‐like scheme, demonstrated very accurate results, high stability, and fast convergence. The fifth‐order Runge–Kutta scheme is used for time discretization. The local time stepping and implicit residual smoothing were applied as the convergence acceleration techniques. Suitable boundary conditions have been implemented considering flow behavior. The problem has been studied at high Reynolds numbers for cross flow around the horizontal circular cylinder and NACA0012 hydrofoil. Results were compared with those of others and a good agreement has been observed. Copyright © 2012 John Wiley & Sons, Ltd.     

High‐accuracy upwind method using improved characteristics speeds for incompressible flows

In this study, the Nervier–Stokes equations for incompressible flows, modified by the artificial compressibility method, are investigated numerically. To calculate the convective fluxes, a new high‐accuracy characteristics‐based (HACB) scheme is presented in this paper. Comparing the HACB scheme with the original characteristic‐based method, it is found that the new proposed scheme is more accurate and has faster convergence rate than the older one. The second order averaging scheme is used for estimating the viscose fluxes, and spatially discretized equations are integrated in time by an explicit fourth‐order Runge–Kutta scheme. The lid driven cavity flow and flow in channel with a backward facing step have been used as benchmark problems. It is shown that the obtained results using HACB scheme are in good agreement with the standard solutions. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

An implicit edge‐based ALE method for the incompressible Navier–Stokes equations

A new finite volume method for the incompressible Navier–Stokes equations, expressed in arbitrary Lagrangian–Eulerian (ALE) form, is presented. The method uses a staggered storage arrangement for the pressure and velocity variables and adopts an edge‐based data structure and assembly procedure which is valid for arbitrary n‐sided polygonal meshes. Edge formulas are presented for assembling the ALE form of the momentum and pressure equations. An implicit multi‐stage time integrator is constructed that is geometrically conservative to the precision of the arithmetic used in the computation. The method is shown to be second‐order‐accurate in time and space for general time‐dependent polygonal meshes. The method is first evaluated using several well‐known unsteady incompressible Navier–Stokes problems before being applied to a periodically forced aeroelastic problem and a transient free surface problem. Published in 2003 by John Wiley & Sons, Ltd.     

A second‐order time accurate finite element method for quasi‐incompressible viscous flows

A finite element method for quasi‐incompressible viscous flows is presented. An equation for pressure is derived from a second‐order time accurate Taylor–Galerkin procedure that combines the mass and the momentum conservation laws. At each time step, once the pressure has been determined, the velocity field is computed solving discretized equations obtained from another second‐order time accurate scheme and a least‐squares minimization of spatial momentum residuals. The terms that stabilize the finite element method (controlling wiggles and circumventing the Babuska–Brezzi condition) arise naturally from the process, rather than being introduced a priori in the variational formulation. A comparison between the present second‐order accurate method and our previous first‐order accurate formulation is shown. The method is also demonstrated in the computation of the leaky‐lid driven cavity flow and in the simulation of a crossflow past a circular cylinder. In both cases, good agreement with previously published experimental and computational results has been obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

On the efficiency of linearization schemes and coupled multigrid methods in the simulation of a 3D flow around a cylinder

This paper studies the efficiency of two ways to treat the non‐linear convective term in the time‐dependent incompressible Navier–Stokes equations and of two multigrid approaches for solving the arising linear algebraic saddle point problems. The Navier–Stokes equations are discretized by a second‐order implicit time stepping scheme and by inf–sup stable, higher order finite elements in space. The numerical studies are performed at a 3D flow around a cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

A comparative study of high‐order variable‐property segregated algorithms for unsteady low Mach number flows

In this paper, five different algorithms are presented for the simulation of low Mach flows with large temperature variations, based on second‐order central‐difference or fourth‐order compact spatial discretization and a pressure projection‐type method. A semi‐implicit three‐step Runge–Kutta/Crank–Nicolson or second‐order iterative scheme is used for time integration. The different algorithms solve the coupled set of governing scalar equations in a decoupled segregate manner. In the first algorithm, a temperature equation is solved and density is calculated from the equation of state, while the second algorithm advances the density using the differential form of the equation of state. The third algorithm solves the continuity equation and the fourth algorithm solves both the continuity and enthalpy equation in conservative form. An iterative decoupled algorithm is also proposed, which allows the computation of the fully coupled solution. All five algorithms solve the momentum equation in conservative form and use a constant‐ or variable‐coefficient Poisson equation for the pressure. The efficiency of the fourth‐order compact scheme and the performances of the decoupling algorithms are demonstrated in three flow problems with large temperature variations: non‐Boussinesq natural convection, channel flow instability, flame–vortex interaction. Copyright © 2010 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.     

An accurate pressure–velocity decoupling technique for semi‐implicit rotational projection methods

In this paper, an accurate semi‐implicit rotational projection method is introduced to solve the Navier–Stokes equations for incompressible flow simulations. The accuracy of the fractional step procedure is investigated for the standard finite‐difference method, and the discrete forms are presented with arbitrary orders or accuracy. In contrast to the previous semi‐implicit projection methods, herein, an alternative way is proposed to decouple pressure from the momentum equation by employing the principle form of the pressure Poisson equation. This equation is based on the divergence of the convective terms and incorporates the actual pressure in the simulations. As a result, the accuracy of the method is not affected by the common choice of the pseudo‐pressure in the previous methods. Also, the velocity correction step is redefined, and boundary conditions are introduced accordingly. Several numerical tests are conducted to assess the robustness of the method for second and fourth orders of accuracy. The results are compared with the solutions obtained from a typical high‐resolution fully explicit method and available benchmark reports. Herein, the numerical tests are consisting of simulations for the Taylor–Green vortex, lid‐driven square cavity, and vortex–wall interaction. It is shown that the present method can preserve the order of accuracy for both velocity and pressure fields in second‐order and high‐order simulations. Furthermore, a very good agreement is observed between the results of the present method and benchmark simulations. Copyright © 2016 John Wiley & Sons, Ltd.     

A projection method for the spectral solution of non‐homogeneous and incompressible Navier–Stokes equations

This paper is devoted to the development of a parallel, spectral and second‐order time‐accurate method for solving the incompressible and variable density Navier–Stokes equations. The method is well suited for finite thickness density layers and is very efficient, especially for three‐dimensional computations. It is based on an exact projection technique. To enforce incompressibility, for a non‐homogeneous fluid, the pressure is computed using an iterative algorithm. A complete study of the convergence properties of this algorithm is done for different density variations. Numerical simulations showing, qualitatively, the capabilities of the developed Navier–Stokes solver for many realistic problems are presented. The numerical procedure is also validated quantitatively by reproducing growth rates from the linear instability theory in a three‐dimensional direct numerical simulation of an unstable, non‐homogeneous, flow configuration. It is also shown that, even in a turbulent flow, the spectral accuracy is recovered. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of level set method with good area preservation to predict interface in two‐phase flows

A two‐step conservative level set method is proposed in this study to simulate the gas/water two‐phase flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the coupled compact scheme. For accurately predicting the modified level set function, the dispersion‐relation‐preserving advection scheme is developed to preserve the theoretical dispersion relation for the first‐order derivative terms shown in the pure advection equation cast in conservative form. For the purpose of retaining its long‐time accurate Casimir functionals and Hamiltonian in the transport equation for the level set function, the time derivative term is discretized by the sixth‐order accurate symplectic Runge–Kutta scheme. To resolve contact discontinuity oscillations near interface, nonlinear compression flux term and artificial damping term are properly added to the second‐step equation of the modified level set method. For the verification of the proposed dispersion‐relation‐preserving scheme applied in non‐staggered grids for solving the incompressible flow equations, three benchmark problems have been chosen in this study. The conservative level set method with area‐preserving property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam‐break, Rayleigh–Taylor instability, bubble rising in water, and droplet falling in water problems. Good agreements with the referenced solutions are demonstrated in all the investigated problems. Copyright © 2010 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.     

Three‐dimensional transient Navier–Stokes solvers in cylindrical coordinate system based on a spectral collocation method using explicit treatment of the pressure

A spectral collocation method is developed for solving the three‐dimensional transient Navier–Stokes equations in cylindrical coordinate system. The Chebyshev–Fourier spectral collocation method is used for spatial approximation. A second‐order semi‐implicit scheme with explicit treatment of the pressure and implicit treatment of the viscous term is used for the time discretization. The pressure Poisson equation enforces the incompressibility constraint for the velocity field, and the pressure is solved through the pressure Poisson equation with a Neumann boundary condition. We demonstrate by numerical results that this scheme is stable under the standard Courant–Friedrichs–Lewy (CFL) condition, and is second‐order accurate in time for the velocity, pressure, and divergence. Further, we develop three accurate, stable, and efficient solvers based on this algorithm by selecting different collocation points in r‐, ? ‐, and z‐directions. Additionally, we compare two sets of collocation points used to avoid the axis, and the numerical results indicate that using the Chebyshev Gauss–Radau points in radial direction to avoid the axis is more practical for solving our problem, and its main advantage is to save the CPU time compared with using the Chebyshev Gauss–Lobatto points in radial direction to avoid the axis. Copyright © 2010 John Wiley & Sons, Ltd.