Fully computable a posteriori error bounds for stabilised FEM approximations of convection–reaction–diffusion problems in three dimensions

Fully computable upper bounds are developed for the discretisation error measured in the natural (energy) norm for convection–reaction–diffusion problems in three dimensions. The upper bounds are genuine upper bounds in the sense that the numerical value of the estimated error exceeds the actual numerical value of the true error regardless of the coarseness of the mesh or the nature of the data for the problem. All constants appearing in the bounds are fully specified. Examples show the estimator to be reliable and accurate even in the case of complicated three‐dimensional problems. Copyright © 2013 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.     

Development of a class of multiple time‐stepping schemes for convection–diffusion equations in two dimensions

In this paper we present a class of semi‐discretization finite difference schemes for solving the transient convection–diffusion equation in two dimensions. The distinct feature of these scheme developments is to transform the unsteady convection–diffusion (CD) equation to the inhomogeneous steady convection–diffusion‐reaction (CDR) equation after using different time‐stepping schemes for the time derivative term. For the sake of saving memory, the alternating direction implicit scheme of Peaceman and Rachford is employed so that all calculations can be carried out within the one‐dimensional framework. For the sake of increasing accuracy, the exact solution for the one‐dimensional CDR equation is employed in the development of each scheme. Therefore, the numerical error is attributed primarily to the temporal approximation for the one‐dimensional problem. Development of the proposed time‐stepping schemes is rooted in the Taylor series expansion. All higher‐order time derivatives are replaced with spatial derivatives through use of the model differential equation under investigation. Spatial derivatives with orders higher than two are not taken into account for retaining the linear production term in the convection–diffusion‐reaction differential system. The proposed schemes with second, third and fourth temporal accuracy orders have been theoretically explored by conducting Fourier and dispersion analyses and numerically validated by solving three test problems with analytic solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

A multi‐dimensional monotonic finite element model for solving the convection–diffusion‐reaction equation

In this paper, we develop a finite element model for solving the convection–diffusion‐reaction equation in two dimensions with an aim to enhance the scheme stability without compromising consistency. Reducing errors of false diffusion type is achieved by adding an artificial term to get rid of three leading mixed derivative terms in the Petrov–Galerkin formulation. The finite element model of the Petrov–Galerkin type, while maintaining convective stability, is modified to suppress oscillations about the sharp layer by employing the M‐matrix theory. To validate this monotonic model, we consider test problems which are amenable to analytic solutions. Good agreement is obtained with both one‐ and two‐dimensional problems, thus validating the method. Other problems suitable for benchmarking the proposed model are also investigated. Copyright © 2002 John Wiley & Sons, Ltd.     

Fourier analysis of several finite difference schemes for the one‐dimensional unsteady convection–diffusion equation

This paper reports a comparative study on the stability limits of nine finite difference schemes to discretize the one‐dimensional unsteady convection–diffusion equation. The tested schemes are: (i) fourth‐order compact; (ii) fifth‐order upwind; (iii) fourth‐order central differences; (iv) third‐order upwind; (v) second‐order central differences; and (vi) first‐order upwind. These schemes were used together with Runge–Kutta temporal discretizations up to order six. The remaining schemes are the (vii) Adams–Bashforth central differences, (viii) the Quickest and (ix) the Leapfrog central differences. In addition, the dispersive and dissipative characteristics of the schemes were compared with the exact solution for the pure advection equation, or simple first or second derivatives, and numerical experiments confirm the Fourier analysis. The results show that fourth‐order Runge–Kutta, together with central schemes, show good conditional stability limits and good dispersive and dissipative spectral resolution. Overall the fourth‐order compact is the recommended scheme. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

Parallel finite element computation of unsteady incompressible flows

A parallel semi-explicit iterative finite element computational procedure for modelling unsteady incompressible fluid flows is presented. During the procedure, element flux vectors are calculated in parallel and then assembled into global flux vectors. Equilibrium iterations which introduce some ‘local implicitness’ are performed at each time step. The number of equilibrium iterations is governed by an implicitness parameter. The present technique retains the advantages of purely explicit schemes, namely (i) the parallel speed-up is equal to the number of parallel processors if the small communication overhead associated with purely explicit schemes is ignored and (ii) the computation time as well as the core memory required is linearly proportional to the number of elements. The incompressibility condition is imposed by using the artificial compressibility technique. A pressure-averaging technique which allows the use of equal-order interpolations for both velocity and pressure, this simplifying the formulation, is employed. Using a standard Galerkin approximation, three benchmark steady and unsteady problems are solved to demonstrate the accuracy of the procedure. In all calculations the Reynolds number is less than 500. At these Reynolds numbers it was found that the physical dissipation is sufficient to stabilize the convective term with no need for additional upwind-type dissipation. © 1998 John Wiley & Sons, Ltd.     

An unsteady incompressible Navier–Stokes solver for large eddy simulation of turbulent flows

An unsteady incompressible Navier–Stokes solver that uses a dual time stepping method combined with spatially high‐order‐accurate finite differences, is developed for large eddy simulation (LES) of turbulent flows. The present solver uses a primitive variable formulation that is based on the artificial compressibility method and various convergence–acceleration techniques are incorporated to efficiently simulate unsteady flows. A localized dynamic subgrid model, which is formulated using the subgrid kinetic energy, is employed for subgrid turbulence modeling. To evaluate the accuracy and the efficiency of the new solver, a posteriori tests for various turbulent flows are carried out and the resulting turbulence statistics are compared with existing experimental and direct numerical simulation (DNS) data. Copyright © 1999 John Wiley & Sons, Ltd.     

A simplified circular function–based gas kinetic scheme for simulation of incompressible flows

In this paper, the circular function–based gas kinetic scheme (GKS), which is often applied for simulation of compressible flows, is simplified to improve computational efficiency for simulation of incompressible flows. In the original circular function–based GKS, the integral domain along the circle for computing conservative variables and numerical fluxes is usually not symmetric at the cell interface. This leads to relatively complicated formulations for computing the numerical flux at the cell interface. As shown in this work, for incompressible flows, the circle at the cell interface can be approximately considered to be symmetric. As a consequence, the simple expressions for calculation of conservative variables and numerical fluxes at the cell interface can be obtained, and computational efficiency is greatly improved. In the meanwhile, like the original circular function–based GKS, the discontinuity of conservative variables and their derivatives at the cell interface is still kept in the present scheme to keep good numerical stability at high Reynolds numbers. Several numerical examples, including decaying vortex flow, lid‐driven cavity flow, and flow past a stationary and rotating circular cylinder, are tested to validate the accuracy, efficiency, and stability of the present scheme.     

A new weighted upwind finite volume element method based on non‐standard covolume for time‐dependent convection–diffusion problems

In modern numerical simulation of problems in energy resources and environmental science, it is important to develop efficient numerical methods for time‐dependent convection–diffusion problems. On the basis of nonstandard covolume grids, we propose a new kind of high‐order upwind finite volume element method for the problems. We first prove the stability and mass conservation in the discrete forms of the scheme. Optimal second‐order error estimate in L₂‐norm in spatial step is then proved strictly. The scheme is effective for avoiding numerical diffusion and nonphysical oscillations and has second‐order accuracy. Numerical experiments are given to verify the performance of the scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

A note on two upwind strategies for RBF‐based grid‐free schemes to solve steady convection–diffusion equations

In this paper, two radial basis function (RBF)‐based local grid‐free upwind schemes have been discussed for convection–diffusion equations. The schemes have been validated over some convection–diffusion problems with sharp boundary layers. It is found that one of the upwind schemes realizes the boundary layers more accurately than the rest. Comparisons with the analytical solutions demonstrate that the local RBF grid‐free upwind schemes based on the exact velocity direction are stable and produce accurate results on domains discretized even with scattered distribution of nodal points. Copyright © 2009 John Wiley & Sons, Ltd.     

Bubble‐based stabilized finite element methods for time‐dependent convection–diffusion–reaction problems

In this paper, we propose a numerical algorithm for time‐dependent convection–diffusion–reaction problems and compare its performance with the well‐known numerical methods in the literature. Time discretization is performed by using fractional‐step θ‐scheme, while an economical form of the residual‐free bubble method is used for the space discretization. We compare the proposed algorithm with the classical stabilized finite element methods over several benchmark problems for a wide range of problem configurations. The effect of the order in the sequence of discretization (in time and in space) to the quality of the approximation is also investigated. Numerical experiments show the improvement through the proposed algorithm over the classical methods in either cases. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of a convection–diffusion‐reaction magnetohydrodynamic solver on non‐staggered grids

This paper presents a convection–diffusion‐reaction (CDR) model for solving magnetic induction equations and incompressible Navier–Stokes equations. For purposes of increasing the prediction accuracy, the general solution to the one‐dimensional constant‐coefficient CDR equation is employed. For purposes of extending this discrete formulation to two‐dimensional analysis, the alternating direction implicit solution algorithm is applied. Numerical tests that are amenable to analytic solutions were performed in order to validate the proposed scheme. Results show good agreement with the analytic solutions and high rate of convergence. Like many magnetohydrodynamic studies, the Hartmann–Poiseuille problem is considered as a benchmark test to validate the code. Copyright © 2004 John Wiley & Sons, Ltd.     

An upwind finite‐volume element scheme and its maximum‐principle‐preserving property for nonlinear convection–diffusion problem

For a class of nonlinear convection–diffusion equation in multiple space dimensions, a kind of upwind finite‐volume element (UFVE) scheme is put forward. Some techniques, such as calculus of variations, commutating operators and prior estimates, are adopted. It is proved that the UFVE scheme is unconditionally stable and satisfies maximum principle. Optimal‐order estimates in H¹‐norm are derived to determine the error in the approximate solution. Numerical results are presented to observe the performance of the scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

A transformation‐free HOC scheme for steady convection–diffusion on non‐uniform grids

A higher order compact (HOC) finite difference solution procedure has been proposed for the steady two‐dimensional (2D) convection–diffusion equation on non‐uniform orthogonal Cartesian grids involving no transformation from the physical space to the computational space. Effectiveness of the method is seen from the fact that for the first time, an HOC algorithm on non‐uniform grid has been extended to the Navier–Stokes (N–S) equations. Apart from avoiding usual computational complexities associated with conventional transformation techniques, the method produces very accurate solutions for difficult test cases. Besides including the good features of ordinary HOC schemes, the method has the advantage of better scale resolution with smaller number of grid points, with resultant saving of memory and CPU time. Gain in time however may not be proportional to the decrease in the number of grid points as grid non‐uniformity imparts asymmetry to some of the associated matrices which otherwise would have been symmetric. The solution procedure is also highly robust as it computes complex flows such as that in the lid‐driven square cavity at high Reynolds numbers (Re), for which no HOC results have so far been seen. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of Galerkin and control volume finite element for advection–diffusion problems

The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass‐lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov–Galerkin streamline upwind schemes in the context of advection–diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection‐dominated problems. Published in 2005 by John Wiley & Sons, Ltd.     

A multiscale control volume finite element method for advection–diffusion equations

We present a new stabilized method for advection–diffusion equations, which combines a control volume FEM formulation of the governing equations with a novel multiscale approximation of the total flux. The latter incorporates information about the exact solution that cannot be represented on the mesh. To define this flux, we solve the governing equations along suitable mesh segments under the assumption that the flux varies linearly along these segments. This procedure yields second‐order accurate fluxes on the edges of the mesh. Then, we use curl‐conforming elements of the same order to lift these edge fluxes into the mesh elements. In so doing, we obtain a stabilized control volume FEM formulation that is second‐order accurate and does not require mesh‐dependent stabilization parameters. Numerical convergence studies on uniform and nonuniform grids along with several standard advection tests illustrate the computational properties of the new method. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.