A co‐ordinate system independent streamwise upwind method for fluid flow computation

A new co‐ordinate invariant streamwise upwind formulation for convection dominated flows is developed. The eddy diffusivity/viscosity is added directly to the equations in order to remove the oscillations in the solution. The equations then can be solved by any high‐order scheme and the solution retains the accuracy of the high‐order scheme. The accuracy and reduced lateral thickness growth rate are demonstrated with several numerical examples, including pure convective flows and lid‐driven cavity flow. The lateral spreading due to the numerical diffusion is controlled by the anisotropic tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.     

对流项占优问题的MLPG/SUPG方法数值模拟

在计算对流项占优问题时易产生假扩散,本文把流线型迎风格式应用于MLPG方法中可以减少对流项的影响,通过两个典型例子（旋转流场问题和Brezzi问题）验证该格式的精度与有效性,并与文献中的迎风格式的计算结果进行比较,计算结果表明,该方法能有效地克服假扩散现象,有较好的稳定性和较高的计算精度。     

Development and validation of a SUPG finite element scheme for the compressible Navier–Stokes equations using a modified inviscid flux discretization

This paper considers the streamline‐upwind Petrov–Galerkin (SUPG) method applied to the unsteady compressible 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, and the second‐order accurate time discretization are presented. The numerical method is discussed in detail. The performance of the algorithm is then investigated by considering inviscid flow past a circular cylinder. Validation of the finite element formulation via comparisons with experimental data for high‐Mach number perfect gas laminar flows is presented, with a specific focus on comparisons with experimentally measured skin friction and convective heat transfer on a 15° compression ramp. Copyright © 2007 John Wiley & Sons, Ltd.     

On streamline diffusion arising in Galerkin FEM with predictor/multi‐corrector time integration

In the present paper, the author shows that the predictor/multi‐corrector (PMC) time integration for the advection–diffusion equations induces numerical diffusivity acting only in the streamline direction, even though the equations are spatially discretized by the conventional Galerkin finite element method (GFEM). The transient 2‐D and 3‐D advection problems are solved with the PMC scheme using both the GFEM and the streamline upwind/Petrov Galerkin (SUPG) as the spatial discretization methods for comparison. The solutions of the SUPG‐PMC turned out to be overly diffusive due to the additional PMC streamline diffusion, while the solutions of the GFEM‐PMC were comparatively accurate without significant damping and phase error. A similar tendency was seen also in the quasi‐steady solutions to the incompressible viscous flow problems: 2‐D driven cavity flow and natural convection in a square cavity. Copyright © 2002 John Wiley & Sons, Ltd.     

An adaptive viscoelastic stress splitting scheme and its applications: AVSS/SI and AVSS/SUPG

We report an adaptive viscoelastic stress splitting (AVSS) scheme, which was found to be extremely robust in the numerical simulation of viscoelastic flow involving steep stress boundary layers. The scheme is different from the elastic viscous split stress (EVSS) in that the local Newtonian component is allowed to depend adaptively on the magnitude of the local elastic stress. Two decoupled versions of the scheme were implemented for the Upper Convected Maxwell (UCM) model: the streamline integration (AVSS/SI), and the streamline upwind Petrov-Galerkin (AVSS/SUPG) methods of integrating the stress. The implementations were benchmarked against the known analytic Poiseuille solution, and no upper limiting Weissenberg number was found (the computation was stopped at Weissenberg number of O(10⁴)). The flow past a sphere in a tube was solved next, giving convergent results up to a Weissenberg number of 3.2 with the AVSS/SI and 1.55 with the AVSS/SUPG (both were decoupled schemes; without the adaptive scheme, the limiting Weissenberg number for the decoupled streamline integration method was about 0.3). These results show that the decoupled AVSS is more stable than the corresponding EVSS, and the SI is more robust than SUPG in solving the constitutive equation of hyperbolic type. In addition, we found a very long stress wake behind the sphere, and a weak vortex in the rear stagnation region at a Weissenberg number above W_i of about 1.6, which does not seem to increase in size or strength with increasing W_i.     

An adaptive wavelet space‐time SUPG method for hyperbolic conservation laws

This article concerns with incorporating wavelet bases into existing streamline upwind Petrov‐Galerkin (SUPG) methods for the numerical solution of nonlinear hyperbolic conservation laws which are known to develop shock solutions. Here, we utilize an SUPG formulation using continuous Galerkin in space and discontinuous Galerkin in time. The main motivation for such a combination is that these methods have good stability properties thanks to adding diffusion in the direction of streamlines. But they are more expensive than explicit semidiscrete methods as they have to use space‐time formulations. Using wavelet bases we maintain the stability properties of SUPG methods while we reduce the cost of these methods significantly through natural adaptivity of wavelet expansions. In addition, wavelet bases have a hierarchical structure. We use this property to numerically investigate the hierarchical addition of an artificial diffusion for further stabilization in spirit of spectral diffusion. Furthermore, we add the hierarchical diffusion only in the vicinity of discontinuities using the feature of wavelet bases in detection of location of discontinuities. Also, we again use the last feature of the wavelet bases to perform a postprocessing using a denosing technique based on a minimization formulation to reduce Gibbs oscillations near discontinuities while keeping other regions intact. Finally, we show the performance of the proposed combination through some numerical examples including Burgers’, transport, and wave equations as well as systems of shallow water equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2062–2089, 2017     

Local projection stabilization for advection–diffusion–reaction problems: One-level vs. two-level approach

Local projection stabilization (LPS) of finite element methods is a new technique for the numerical solution of transport-dominated problems. The main aim of this paper is a critical discussion and comparison of the one- and two-level approaches to LPS for the linear advection–diffusion–reaction problem. Moreover, the paper contains several other novel contributions to the theory of LPS. In particular, we derive an error estimate showing not only the usual error dependence on the mesh width but also on the polynomial degree of the finite element space. Based on this error estimate, we propose a definition of the stabilization parameter depending on the data of the solved problem. Unlike other papers on LPS methods, we observe that the consistency error may deteriorate the convergence order. Finally, we explain the relation between the LPS method and residual-based stabilization techniques for simplicial finite elements.     

对流反应扩散方程的稳定化时间间断时空有限元解的误差估计

在流线迎风Petrov-Galerkin（SUPG）稳定化有限元数值格式的基础上,结合时间方向的变分离散,构造对流反应扩散方程的稳定化时间间断时空有限元格式.该类格式在工程上有一些数值模拟应用,但相关文献没有看到类似数值格式的理论证明.本文以Radau点为节点,构造时间方向的Lagrange插值多项式,证明了稳定化有限元解的稳定性,时间最大模、空间L²（Ω）-模误差估计.文中利用插值多项式和有限元方法相结合的技巧,解耦时空变量,去掉了时空网格的限制条件,提供了时间间断稳定化时空有限元方法的理论证明思路,克服了因时空变量统一导致的实际计算时的复杂性.