Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

On the role of (weak) compressibility for fluid-structure interaction solvers

In the present study, a weakly compressible formulation of the Navier-Stokes equations is developed and examined for the solution of fluid-structure interaction (FSI) problems. Newtonian viscous fluids under isothermal conditions are considered, and the Murnaghan-Tait equation of state is employed for the evaluation of mass density changes with pressure. A pressure-based approach is adopted to handle the low Mach number regime, ie, the pressure is chosen as primary variable, and the divergence-free condition of the velocity field for incompressible flows is replaced by the continuity equation for compressible flows. The approach is then embedded into a partitioned FSI solver based on a Dirichlet-Neumann coupling scheme. It is analytically demonstrated how this formulation alleviates the constraints of the instability condition of the artificial added mass effect, due to the reduction of the maximal eigenvalue of the so-called added mass operator. The numerical performance is examined on a selection of benchmark problems. In comparison to a fully incompressible solver, a significant reduction of the coupling iterations and the computational time and a notable increase in the relaxation parameter evaluated according to Aitken's Δ² method are observed.     

A collocated finite volume method for predicting flows at all speeds

An existing two-dimensional method for the prediction of steady-state incompressible flows in complex geometry is extended to treat also compressible flows at all speeds. The primary variables are the Cartesian velocity components, pressure and temperature. Density is linked to pressure via an equation of state. The influence of pressure on density in the case of compressible flows is implicitly incorporated into the extended SIMPLE algorithm, which in the limit of incompressible flow reduces to its well-known form. Special attention is paid to the numerical treatment of boundary conditions. The method is verified on a number of test cases (inviscid and viscous flows), and both the results and convergence properties compare favourably with other numerical results available in the literature.     

A multiblock operator‐splitting algorithm for unsteady flows at all speeds in complex geometries

This paper describes a non‐iterative operator‐splitting algorithm for computing all‐speed flows in complex geometries. A pressure‐based algorithm is adopted as the base, in which pressure, instead of density, is a primary variable, thus allowing for a unified formulation for all Mach numbers. The focus is on adapting the method for (a) flows at all speeds, and (b) multiblock, non‐orthogonal, body‐fitted grids for very complex geometries. Key features of the formulation include special treatment of mass fluxes at control volume interfaces to avoid pressure–velocity decoupling for incompressible (low Mach number limit) flows and to provide robust pressure–velocity–density coupling for compressible (high‐speed) flows. The method is shown to be robust for all Mach number regimes for both steady and unsteady flows; it is found to be stable for CFL numbers of order ten, allowing large time steps to be taken for steady flows. Enhancements to the method which allow for stable solutions to be obtained on non‐orthogonal grids are also discussed. The method is found to be very reliable even in complex engineering applications such as unsteady rotor–stator interactions in turbulent, all‐speed turbomachinery flows. Copyright © 2004 John Wiley & Sons, Ltd.     

A Mach‐uniform preconditioner for incompressible and subsonic flows

In this study, a novel Mach‐uniform preconditioning method is developed for the solution of Euler equations at low subsonic and incompressible flow conditions. In contrast to the methods developed earlier in which the conservation of mass equation is preconditioned, in the present method, the conservation of energy equation is preconditioned, which enforces the divergence free constraint on the velocity field even at the limiting case of incompressible, zero Mach number flows. Despite most preconditioners, the proposed Mach‐uniform preconditioning method does not have a singularity point at zero Mach number. The preconditioned system of equations preserves the strong conservation form of Euler equations for compressible flows and recovers the artificial compressibility equations in the case of zero Mach number. A two‐dimensional Euler solver is developed for validation and performance evaluation of the present formulation for a wide range of Mach number flows. The validation cases studied show the convergence acceleration, stability, and accuracy of the present Mach‐uniform preconditioner in comparison to the non‐preconditioned compressible flow solutions. The convergence acceleration obtained with the present formulation is similar to those of the well‐known preconditioned system of equations for low subsonic flows and to those of the artificial compressibility method for incompressible flows. Copyright © 2013 John Wiley & Sons, Ltd.     

A finite element method for compressible viscous fluid flows

This work presents a mixed three‐dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a ‘stable’ numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf–sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady‐state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright © 2010 John Wiley & Sons, Ltd.     

A non‐conforming least‐squares finite element method for incompressible fluid flow problems

In this paper, we develop least‐squares finite element methods (LSFEMs) for incompressible fluid flows with improved mass conservation. Specifically, we formulate a new locally conservative LSFEM for the velocity–vorticity–pressure Stokes system, which uses a piecewise divergence‐free basis for the velocity and standard C⁰ elements for the vorticity and the pressure. The new method, which we term dV‐VP improves upon our previous discontinuous stream‐function formulation in several ways. The use of a velocity basis, instead of a stream function, simplifies the imposition and implementation of the velocity boundary condition, and eliminates second‐order terms from the least‐squares functional. Moreover, the size of the resulting discrete problem is reduced because the piecewise solenoidal velocity element is approximately one‐half of the dimension of a stream‐function element of equal accuracy. In two dimensions, the discontinuous stream‐function LSFEM [1] motivates modification of our functional, which further improves the conservation of mass. We briefly discuss the extension of this modification to three dimensions. Computational studies demonstrate that the new formulation achieves optimal convergence rates and yields high conservation of mass. We also propose a simple diagonal preconditioner for the dV‐VP formulation, which significantly reduces the condition number of the LSFEM problem. Published 2012. This article is a US Government work and is in the public domain in the USA.     

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.     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite strain––isotropic hyperelasticity

This paper presents a strain energy density for isotropic hyperelastic materials. The strain energy density is decomposed into a compressible and incompressible component. The incompressible component is the same as the generalized Mooney expression while the compressible component is shown to be a function of the volume invariant J only. The strain energy density proposed is used to investigate problems involving incompressible isotropic materials such as rubber under homogeneous strain, compressible isotropic materials under high hydrostatic pressure and volume change under uniaxial tension. Comparison with experimental data is good. The formulation is also used to derive a strain energy density expression for compressible isotropic neo-Hookean materials. The constitutive relationship for the second Piola–Kirchhoff stress tensor and its physical counterpart, involves the contravariant Almansi strain tensor. The stress stretch relationship comprises of a component associated with volume constrained distortion and a hydrostatic pressure which results in volumetric dilation. An important property of this constitutive relationship is that the hydrostatic pressure component of the stress vector which is associated with volumetric dilation will have no shear component on any surface in any configuration. This same property is not true for a neo-Hookean Green’s strain–second Piola–Kirchhoff stress tensor formulation.     

A compressible single‐temperature conservative two‐phase model with phase transitions

A model for multidimensional compressible two‐phase flow with pressure and velocity relaxations based on the theory of thermodynamically compatible system is extended to study liquid–gas flows with cavitation. The model assumes for each phase its own pressure and velocity, while a common temperature is considered. The governing equations form a hyperbolic system in conservative form and are derived through the theory of a thermodynamically compatible system. The phase pressure‐equalizing process and the interfacial friction are taken into account in the balance laws for the volume fractions of one phase and for the relative velocity by adding two relaxation source terms, while the phase transition is introduced into the model considering in the balance equation for the mass of one phase the relaxation of the Gibbs free energies of the two phases. A modification of the central finite‐volume Kurganov–Noelle–Petrova method is adopted in this work to solve the homogeneous hyperbolic part, while the relaxation source terms are treated implicitly. In order to investigate the effect of the mass transfer in the solution, a 1D cavitation tube problem is presented. In addition, two 2D numerical simulations regarding cavitation problem are also studied: a cavitating Richtmyer–Meshkov instability and a laser‐induced cavitation problem. Copyright © 2014 John Wiley & Sons, Ltd.     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

GMRES physics-based preconditioner for all Reynolds and Mach numbers: numerical examples

This paper presents several numerical results using a vectorized version of a 3D finite element compressible and nearly incompressible Euler and Navier–Stokes code. The assumptions were set on laminar flows and Newtonian fluids. The goal of this research is to show the capabilities of the present code to treat a wide range of problems appearing in laminar fluid dynamics towards the unification from incompressible to compressible and from inviscid to viscous flow codes. Several authors with different approaches have tried to attain this target in CFD with relative success. At the beginning the methods based on operator splitting and perturbation were preferred, but lately, with the wide usage of time-marching algorithms, the preconditioning mass matrix (PMM) has become very popular. With this kind of relaxation scheme it is possible to accelerate the rate of convergence to steady state solutions with the modification of the mass matrix under certain restrictions. The selection of the mass matrix is not an easy task, but we have certain freedom to define it in order to improve the condition number of the system. In this paper we have used a physics-based preconditioner for the GMRES implicit solver developed previously by us and an SUPG formulation for the semidiscretization of the spatial operator. In sections 2 and 3 we present some theoretical aspects related to the physical problem and the mathematical model, showing the inviscid and viscous flow equations to be solved and the variational formulation involved in the finite element analysis. Section 4 deals with the numerical solution of non-linear systems of equations, with some emphasis on the preconditioned matrix-free GMRES solver. Section 5 shows how boundary conditions were treated for both Euler and Navier–Stokes problems. Section 6 contains some aspects about vectorization on the Cray C90. The performance reached by this implementation is close to 1 Gflop using multitasking. Section 7 presents several numerical examples for both models covering a wide range of interesting problems, such as inviscid low subsonic, transonic and supersonic regimes and viscous problems with interaction between boundary layers and shock waves in either attached or separated flows. © 1997 John Wiley & Sons, Ltd.     

A numerical‐variational procedure for laminar flow in curved square ducts

A new numerical method is presented for the solution of the Navier–Stokes and continuity equations governing the internal incompressible flows. The method denoted as the CVP method consists in the numerical solution of these equations in conjunction with three additional variational equations for the continuity, the vorticity and the pressure field, using a non‐staggered grid. The method is used for the study of the characteristics of the laminar fully developed flows in curved square ducts. Numerical results are presented for the effects of the flow parameters like the curvature, the Dean number and the stream pressure gradient on the velocity distributions, the friction factor and the appearance of a pair of vortices in addition to those of the familiar secondary flow. The accuracy of the method is discussed and the results are compared with those obtained by us, using a variation of the velocity–pressure linked equation methods denoted as the PLEM method and the results obtained by other methods. Copyright © 2004 John Wiley & Sons, Ltd.     

A staggered control volume scheme for unstructured triangular grids

The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2D/axisymmetric incompressible flows. The present study introduces a numerical procedure for solving the Navier–Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a flow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure–velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure field surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume. The method is applied to the simulation of the following test cases: backward-facing step flow, flow over a two-dimensional obstacle and flow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for flow in a driven cavity. Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible flows in complex geometries. © 1997 John Wiley & Sons, Ltd.     

全速解法在湍流跨音速流动中的应用

本文对不可压流动的常用算法SIMPLE算法进行了推广,使其能计算从亚音速到超音速一定马赫数范围内的流动,这里,可压缩流动和不可压缩流动的数值自满实现了统一,称为全速解法本文对全速解法在二维流动计算中的应用性进行了初步的研究,采用了非交错网格的有限体积方法对控制方程进行离散,并用动量插值法来求得连续方程中单元边界上的变量值,本文对全速解法在二维层流的计算效果进行了考核,而后又将此算法在湍流跨音速流动中应用,计算表明本方法是成功的,能够很好地反映各种马赫数下的流场特性。     

A divergence-free semi-implicit finite volume scheme for ideal,viscous, and resistive magnetohydrodynamics

In this paper, we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous, and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum, and total energy, and in multiple space dimensions, it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multidimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes the nonlinear convective terms as well as the time evolution of the magnetic field explicitly, whereas all terms related to the pressure in the momentum equation and the total energy equation are discretized implicitly, making again use of a properly staggered grid for pressure and velocity. Inserting the discrete momentum equation into the discrete energy equation then yields a mildly nonlinear symmetric and positive definite algebraic system for the pressure as the only unknown, which can be efficiently solved with the (nested) Newton method of Casulli et al. The pressure system becomes linear when the specific internal energy is a linear function of the pressure. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Being a semi-implicit pressure-based scheme, our new method is therefore particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate because of the more and more stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense, our new scheme is a true all Mach number flow solver.     

Multiple pressure variables methods for fluid flow at all Mach numbers

In this paper we present a method for the simulation of incompressible as well as compressible unsteady flows. At first we discuss three different forms, i.e. a primitive‐, conservative‐ and a semi‐conservative form of the governing equations. We use a semi‐implicit time integration in such a fashion that the stability is guaranteed independently of the speed of sound and the resulting method is independent of the Mach number range. Moreover, with the application of the so‐called multiple pressure variables (MPV) approach the difficulties with the pressure term can be circumvented as in the incompressible limit the hydrodynamic pressure decouples from the equation of state. Increasing approximation errors in the low Mach number regime are avoided. As a result, the proposed algorithm can also simulate incompressible flows as limit for zero Mach number. Copyright © 2005 John Wiley & Sons, Ltd.     

Non-isothermal polymer melt flow in sudden expansions

This paper presents numerical results for laminar, incompressible and non-isothermal polymer melt flow in sudden expansions. The mathematical model includes the mass, momentum and energy conservation laws within the framework of a generalized Newtonian formulation. Two constitutive relations are adopted to describe the non-Newtonian behavior of the flow, namely Cross and Modified Arrhenius Power-Law models. The governing equations are discretized using the finite difference method based on central, second-order accurate formulas for both convective and diffusive terms. The pressure–velocity coupling is treated by solving a Poisson equation for pressure. The results are presented for two commercial polymers and demonstrate that important flow parameters, such as pressure drop and viscosity distribution, are strongly affected by heat transfer features.