A generalised finite difference scheme based on compact integrated radial basis function for flow in heterogeneous soils

In the present paper, we develop a generalised finite difference approach based on compact integrated radial basis function (CIRBF) stencils for solving highly nonlinear Richards equation governing fluid movement in heterogeneous soils. The proposed CIRBF scheme enjoys a high level of accuracy and a fast convergence rate with grid refinement owing to the combination of the integrated RBF approximation and compact approximation where the spatial derivatives are discretised in terms of the information of neighbouring nodes in a stencil. The CIRBF method is first verified through the solution of ordinary differential equations, 2–D Poisson equations and a Taylor‐Green vortex. Numerical comparisons show that the CIRBF method outperforms some other methods in the literature. The CIRBF method in conjunction with a rational function transformation method and an adaptive time‐stepping scheme is then applied to simulate 1–D and 2–D soil infiltrations effectively. The proposed solutions are more accurate and converge faster than those of the finite different method used with a second‐order central difference scheme. Additionally, the present scheme also takes less time to achieve target accuracy in comparison with the 1D‐IRBF and higher order compact schemes.     

An efficient implicit mesh‐less method for compressible flow calculations

A dual‐time implicit mesh‐less scheme is presented for calculation of compressible inviscid flow equations. The Taylor series least‐square method is used for approximation of spatial derivatives at each node which leads to a central difference discretization. Several convergence acceleration techniques such as local time stepping, enthalpy damping and residual smoothing are adopted in this approach. The capabilities of the method are demonstrated by flow computations around single and multi‐element airfoils at subsonic, transonic and supersonic flow conditions. Results are presented which indicate good agreements with other reliable finite‐volume results. The computational time is considerably reduced when using the proposed mesh‐less method compared with the explicit mesh‐less and finite‐volume schemes using the same point distributions. Copyright © 2010 John Wiley & Sons, Ltd.     

The third‐order polynomial method for two‐dimensional convection and diffusion

Using the upstream polynomial approximation a series of accurate two‐dimensional explicit numerical schemes is developed for the solution of the convection–diffusion equation. A third‐order polynomial approximation (TOP) of the convection term and a consistent second‐order approximation of the diffusion term are combined in a single‐step flux‐difference algorithm. Stability analysis confirms that the TOP‐12 scheme satisfies the CFL condition for two dimensions. Using smaller and narrower flux stencils compared to algorithms of similar accuracy, the TOP‐12 scheme is more efficient in terms of computations per single node. Numerical tests and comparison with other well‐known algorithms show a high performance of the developed schemes. Copyright © 2003 John Wiley & Sons, Ltd.     

Discretization of transport equations on 2D Cartesian unstructured grids using data from remote cells for the convection terms

This paper presents a new finite volume discretization methodology for the solution of transport equations on locally refined or unstructured Cartesian meshes. The implementation of the cell‐face values of the dependent variables enables the employment of data from remote cells and thus the use of higher‐order differencing schemes. It also results in simple and flux‐conservative multiple‐scale stencils for the discretization of the governing equations. The latter are finally cast into a generalized form that does not depend on the local mesh structure. The performance of the numerical model is demonstrated on some classical 2D problems using various gridding techniques and a bounded second‐order upwind scheme. A stable and efficient behaviour of the algorithm is observed in all test cases. The results indicate that the combination in the present model of both local grid refinement and second‐order discretization can produce substantially more accurate solutions than each of the above techniques alone, for the same computational effort. The method is also applicable to turbulent flows and can be easily extended to three‐dimensions. Copyright © 2003 John Wiley & Sons, Ltd.     

An analysis of the stability of the least square finite difference scheme and its shock capturing ability in inviscid flows

A meshless method – The Least Square Finite Difference scheme (LSFD) with diffusion is analyzed and applied to inviscid flows. The scheme is made second-order by using a modified difference in the formulation of LSFD. Several numerical experiments, namely the Sod shock tube and the shallow water problems, are carried out and, in the limelight of the results obtained, the ability of the scheme to resolve shock wave, rarefaction wave, and contact discontinuity is discussed. The conditional stability of the LSFD scheme is established. The LSFD uses weights to diagonalize the least square matrix resulting in the spatial discretization in order to gain computational time. We prove that there exists a unique weight for the resulting optimization problem. The weighted version of LSFD is used to solve the isentropic vortex problem numerically and the results are used to discuss the dissipative nature of the scheme. Five configurations of the two-dimensional Riemann problems are used in our numerical experiments. The capability of the scheme to capture the complex interaction of multiple planar waves is discussed in the limelight of the results on the Riemann problems. The result of the shock reflection problem shows that the scheme is minimally dissipative and leads to sharp and well-resolved shocks.     

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 study of finite element discretization technique in time domain

This paper presents a real time finite element method (TFEM) which unifies discretization techniques both in space domain and time domain and leads to some new possibilities of direct integral schemes. It is proved that some conventional schemes such as central difference and Newmark scheme are only special cases of linear element. The problem of the matching between the division of space mesh and the choice of time step length is also discussed.     

Robust design in aerodynamics using third‐order sensitivity analysis based on discrete adjoint. Application to quasi‐1D flows

In this paper, the second‐order second moment approach, coupled with an adjoint‐based steepest descent algorithm, for the solution of the so‐called robust design problem in aerodynamics is proposed. Because the objective function for the robust design problem comprises first‐order and second‐order sensitivity derivatives with respect to the environmental parameters, the application of a gradient‐based method , which requires the sensitivities of this function with respect to the design variables, calls for the computation of third‐order mixed derivatives. To compute these derivatives with the minimum CPU cost, a combination of the direct differentiation and the discrete adjoint variable method is proposed. This is presented for the first time in the relevant literature and is the most efficient among other possible schemes on condition that the design variables are much more than the environmental ones; this is definitely true in most engineering design problems. The proposed approach was used for the robust design of a duct, assuming a quasi‐1D flow model; the coordinates of the Bézier control points parameterizing the duct shape are used as design variables, whereas the outlet Mach number and the Darcy–Weisbach friction coefficient are used as environmental ones. The extension to 2D and 3D flow problems, after developing the corresponding direct differentiation and adjoint variable methods and software, is straightforward. Copyright © 2011 John Wiley & Sons, Ltd.     

Third‐order Cartesian overset mesh adaptation method for solving steady compressible flows

A third‐order mesh generation and adaptation method is presented for solving the steady compressible Euler equations. For interior points, a third‐order scheme is used on Cartesian and curvilinear meshes. Concerning the mesh adaptation, the method of Meakin is also extended to third order. The accuracy of the new overset mesh adaptation method is demonstrated by a grid convergence study for 2‐D inviscid model problems and results are compared with a second‐order method. Finally, the method is applied to the computation of an inviscid 3‐D flow around a hovering blade of the ONERA 7A helicopter rotor exhibiting an improvement in the wake capture. With a 7 million point mesh, the tip vortex can be followed for more than three rotor revolutions with the third‐order method. The CPU time needed for this calculation is only 3% higher than with a conventional second‐order method. Copyright © 2007 John Wiley & Sons, Ltd.     

A finite volume scheme for solving anisotropic diffusion on ALE‐AMR grids

In the context of High Energy Density Physics and more precisely in the field of laser plasma interaction, Lagrangian schemes are commonly used. The lack of robustness due to strong grid deformations requires the regularization of the mesh through the use of Arbitrary Lagrangian Eulerian methods. Theses methods usually add some diffusion and a loss of precision is observed. We propose to use Adaptive Mesh Refinement (AMR) techniques to reduce this loss of accuracy. This work focuses on the resolution of the anisotropic diffusion operator on Arbitrary Lagrangian Eulerian‐AMR grids. In this paper, we describe a second‐order accurate cell‐centered finite volume method for solving anisotropic diffusion on AMR type grids. The scheme described here is based on local flux approximation which can be derived through the use of a finite difference approximation, leading to the CCLADNS scheme. We present here the 2D and 3D extension of the CCLADNS scheme to AMR meshes. Copyright © 2017 John Wiley & Sons, Ltd.     

Weighted Compact Scheme for Shock Capturing

In this paper a new class of finite difference schemes - the Weighted Compact Schemes are proposed. According to the idea of the WENO schemes, the Weighted Compact Scheme is constructed by a combination of the approximations of derivatives on candidate stencils with properly assigned weights so that the non-oscillatory property is achieve when discontinuities appear. The primitive function reconstruction method of ENO schemes is applied to obtain the conservative form of the Weighted Compact Scheme. This new scheme not only preserves the characteristic of standard compact schemes and achieves high order accuracy and high resolution using a compact stencil, but also can accurately capture shock waves and discontinuities without oscillation. Numerical examples show the new scheme is very promising and successful.     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

A variable relaxation scheme for multiphase,multicomponent flow

We propose a variable relaxation scheme for the simulation of 1D, two-phase, multicomponent flow in porous media. For these strongly nonlinear systems, traditional high order upwind schemes are impractical: Riemann solutions are not directly available when the phase behavior is complex, and the systems are weakly hyperbolic at isolated points. Relaxation schemes avoid the dependency on the eigenstructure and nonlinear Riemann solvers by approximating the original system with a strongly hyperbolic linear system. We exploit the known information about the eigenvalues to construct first order and second order variable relaxation schemes with much reduced numerical diffusion as compared to the standard relaxation formulations. The proposed second order variable relaxation scheme is competitive in accuracy and efficiency with a third order component-wise ENO reconstruction, and performs at least as well as second order component-wise TVD schemes.     

A critical evaluation of seven discretization schemes for convection–diffusion equations

A comparative study of seven discretization schemes for the equations describing convection-diffusion transport phenomena is presented. The (differencing) schemes considered are the conventional central- and upwind-difference schemes, together with the Leonard,¹ Leonard upwind¹ and Leonard super upwind difference¹ schemes. Also tested are the so called locally exact difference scheme² and the quadratic-upstream difference scheme.^3,4 In multidimensional problems errors arise from ‘false-diffusion’ and function approximations. It is asserted that false diffusion is essentially a multidimensional source of error. No mesh constraints are associated with errors in function approximation and discretization. Hence errors associated with discretization only may be investigated via one-dimensional problems. Thus, although the above schemes have been tested for one- and two-dimensional flows with sources, only the former are presented here. For 1D flows, the Leonard super upwind difference scheme and the locally exact scheme are shown to be far superior in accuracy to the others at all Peclet numbers and for most source distributions, for the test cases considered. Furthermore, the latter is shown to be considerably cheaper in computational terms than the former. The stability of the schemes and their CPU time requirements are also discussed.     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids

A novel class of weighted essentially nonoscillatory (WENO) schemes based on Hermite polynomials,termed as HWENO schemes,is developed and applied as limiters for high order discontinuous Galerkin (DG) method on triangular grids.The developed HWENO methodology utilizes high-order derivative information to keep WENO reconstruction stencils in the von Neumann neighborhood.A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils,making higher-order scheme stable and simplifying the reconstruction process at the same time.The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement.Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy,the designed HWENO limiters can simultaneously obtain uniform high order accuracy and sharp,essentially non-oscillatory shock transition.     

一类三维时空二阶精度高分辨率MmB差分格式

直接把Nessyahu和Tadmor^[1,2]的思想推广到三维非线性双曲型守恒律情形，以交错形式Lax—Friedrichs格式为基本模块，使用二阶分片线性逼近代替一阶分片常数逼近，减少了Lax—Friedrichs格式的过多数值粘性，通过对混合导数离散形式的适当处理，构造了一类不须解Riemann问题、具有时空二阶精度高分辨率的MmB差分格式。这些差分格式很容易推广到向量系统中去。最后，一些数值模拟计算结果也证明了这些差分格式的有效性。     

Flux‐difference splitting‐based upwind compact schemes for the incompressible Navier–Stokes equations

Third‐order and fifth‐order upwind compact finite difference schemes based on flux‐difference splitting are proposed for solving the incompressible Navier–Stokes equations in conjunction with the artificial compressibility (AC) method. Since the governing equations in the AC method are hyperbolic, flux‐difference splitting (FDS) originally developed for the compressible Euler equations can be used. In the present upwind compact schemes, the split derivatives for the convective terms at grid points are linked to the differences of split fluxes between neighboring grid points, and these differences are computed by using FDS. The viscous terms are approximated with a sixth‐order central compact scheme. Comparisons with 2D benchmark solutions demonstrate that the present compact schemes are simple, efficient, and high‐order accurate. Copyright © 2008 John Wiley & Sons, Ltd.