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.     

Negative wake generation of FENE-CR fluids in uniform and Poiseuille flows past a cylinder

The effect of flow conditions on the negative wake generation (longitudinal velocity overshoot behind a cylinder in the viscoelastic fluid flow along the centerline) has been investigated. FENE-CR model that predicts constant shear viscosity and controlled extensional viscosity was considered as a constitutive equation. The discrete elastic viscous split stress-G/streamline upwind Petrov–Galerkin (DEVSS-G/SUPG) formulation was employed and the high-resolution solutions were obtained with an efficient iterative solver based on the incomplete LU(0)-type preconditioner and BiCGSTAB. We found that the negative wake generation was more obvious in uniform flow conditions than in Poiseuille flow, which suggests that the experimentally unrevealed negative wake generation of Boger fluids could be partially attributed to the geometrical effect of Poiseuille flow. The negative wake generation was more discernable at low extensibility and high value of viscosity ratio, which agrees well with the previous studies. In addition, we could observe an undershoot phenomenon in Poisseuille flow condition, which has never been reported.     

Direct Finite Element Computation of Periodic Adsorption Processes

This article presents a direct method for computing time-periodic solutions of adsorption processes as an alternative to prolonged dynamic simulation of the natural evolution to periodicity. Direct computation of periodicity is established by discretization on a two-dimensional space-time grid that is periodic in time. Petrov-Galerkin (SUPG) finite element approximation is applied for consistent stabilization of convective terms in the governing hyperbolic equations. Newton iteration with Gaussian elimination (frontal method) is used to solve the resulting set of nonlinear algebraic equations. Computations match exact solutions on simple adsorption cycles, and capture shock layers with as few as two elements. In its present form, the direct method is more efficient than dynamic simulation when the natural evolution to periodicity extends over hundreds of cycles, and will likely be even faster with superior discretization and solution techniques.     

Comments on ‘Adiabatic shock capturing in perfect gas hypersonic flows’

Some comments are provided on the shock‐capturing techniques and stabilization parameters used in a recent paper (B.S. Kirk, Int. J. Numer. Meth. Fluids 2009; DOI: 10.1002/fld.2195 ) in conjunction with the SUPG formulation of compressible flows. Copyright © 2010 John Wiley & Sons, Ltd.     

Adiabatic shock capturing in perfect gas hypersonic flows

This paper considers the streamline‐upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock‐capturing operator, which is required to regularize the scheme in the presence of strong, shock‐induced gradients. A standard shock‐capturing operator that has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator that is designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite‐dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock‐capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator. Published in 2009 by John Wiley & Sons, Ltd.     

Investigation of the numerical accuracy of shallow water flow estimation based on the ensemble Kalman filter using the SUPG FEM

ABSTRACT

The shallow water equation is employed as the governing equation to simulate flow behaviour in shallow water flow regions. The SUPG and the backward Euler methods, respectively, were employed for discretisation in space and time. In this paper, we carry out investigations on the numerical accuracy of shallow water flow estimation based on the ensemble Kalman filter using the SUPG FEM, and show the results of numerical experiments. The open channel model was employed as the numerical example in this study. The extended Kalman filter is generally employed to solve parameter identification problems. However, linearisation of the governing equation is carried out to apply parameter identification: it is known to be difficult to carry out computation reliably if problems that include a high degree of non-linearity need to be solved. On the other hand, linearisation of the governing equation is not carried out in the ensemble Kalman filter, so there is potential for the unknown parameter to be simultaneously identified in the computation of the data assimilation. The distribution of the unknown kinematic viscosity coefficient was therefore also investigated.     

一种新的求解对流占优问题的自适应网格细化方法

在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。     

多维SUPG格式不可压流动高精度算法研究

在传统流线迎风Petrov-Galerkin (SUPG)有限元法基础上,通过对稳定因子关键参数进行分析,提出了基于流向投影的最优特征高度参数确定方法,同时针对不可压流动引入变量分裂算法,发展了一种可用于高雷诺数不可压流动计算的高精度稳定化SUPG方法．详细地给出了三角形单元基于流向的最优特征高度确定方法的分析过程,并给出了基于分裂算法的有限元计算步骤和公式．采用该方法对典型的方腔拖曳、圆柱绕流流动问题进行了分析,在网格较稀疏,且雷诺数较大的情况下,依然可以得到稳定的计算结果,从而验证了该方法的稳定性、有效性,对实际工程应用具有积极的意义．     

Implicit Subgrid-Scale Modeling of a Mach 2.5 Spatially Developing Turbulent Boundary Layer

We employ numerically implicit subgrid-scale modeling provided by the well-known streamlined upwind/Petrov–Galerkin stabilization for the finite element discretization of advection–diffusion problems in a Large Eddy Simulation (LES) approach. Whereas its original purpose was to provide sufficient algorithmic dissipation for a stable and convergent numerical method, more recently, it has been utilized as a subgrid-scale (SGS) model to account for the effect of small scales, unresolvable by the discretization. The freestream Mach number is 2.5, and direct comparison with a DNS database from our research group, as well as with experiments from the literature of adiabatic supersonic spatially turbulent boundary layers, is performed. Turbulent inflow conditions are generated via our dynamic rescaling–recycling approach, recently extended to high-speed flows. Focus is given to the assessment of the resolved Reynolds stresses. In addition, flow visualization is performed to obtain a much better insight into the physics of the flow. A weak compressibility effect is observed on thermal turbulent structures based on two-point correlations (IC vs. supersonic). The Reynolds analogy (u′ vs. t′) approximately holds for the supersonic regime, but to a lesser extent than previously observed in incompressible (IC) turbulent boundary layers, where temperature was assumed as a passive scalar. A much longer power law behavior of the mean streamwise velocity is computed in the outer region when compared to the log law at Mach 2.5. Implicit LES has shown very good performance in Mach 2.5 adiabatic flat plates in terms of the mean flow (i.e., Cf and UVD+). iLES significantly overpredicts the peak values of u′, and consequently Reynolds shear stress peaks, in the buffer layer. However, excellent agreement between the turbulence intensities and Reynolds shear stresses is accomplished in the outer region by the present iLES with respect to the external DNS database at similar Reynolds numbers.     

SUPG-Stabilized Virtual Element Method for Optimal Control Problem Governed by a Convection Dominated Diffusion Equation

In this paper, the streamline upwind/Petrov Galerkin (SUPG) stabilized virtual element method (VEM) for optimal control problem governed by a convection dominated diffusion equation is investigated. The virtual element discrete scheme is constructed based on the first-optimize-then-discretize strategy and SUPG stabilized virtual element approximation of the state equation and adjoint state equation. An a priori error estimate is derived for both the state, adjoint state, and the control. Numerical experiments are carried out to illustrate the theoretical findings.