Fourth-Order Conservative Transport on Overset Grids Using Multi-Moment Constrained Finite Volume Scheme for Oceanic Modeling

With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.

Stability of Semi-implicit Finite Volume Scheme for Level Set Like Equation

We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation （image selective smoothing model） given by Alvarez et al. （Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29： 845-866）. Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1＇2 （W1＇1） sense in isotropic （anisotropic） diffu- sion domain.

A Compact Finite Volume Scheme for the Multi-Term Time Fractional Sub-Diffusion Equation

In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L∞-norm. The convergence order is O(τ2-α + h4). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.

AMaximum-Principle-Preserving Finite Volume Scheme for Diffusion Problems on Distorted Meshes

In this paper,we propose an approach for constructing conservative and maximum-principle-preserving finite volume schemes by using the method of undetermined coefficients,which depend nonlinearly on the linear non-conservative onesided fluxes.In order to facilitate the derivation of expressions of these undetermined coefficients,we explicitly provide a simple constriction condition with a scaling parameter.Such constriction conditions can ensure the final schemes are exact for linear solution problems and may induce various schemes by choosing different values for the parameter.In particular,when this parameter is taken to be 0,the nonlinear terms in our scheme degenerate to a harmonic average combination of the discrete linear fluxes,which has often been used in a variety of maximum-principle-preserving finite volume schemes.Thus our method of determining the coefficients of the nonlinear terms is more general.In addition,we prove the convergence of the proposed schemes by using a compactness technique.Numerical results demonstrate that our schemes can preserve the conservation property,satisfy the discrete maximum principle,possess a second-order accuracy,be exact for linear solution problems,and be available for anisotropic problems on distorted meshes.

Non-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes

This is the continuation of the paper”Central discontinuous Galerkin methods on overlapping cells with a non-oscillatory hierarchical reconstruction”by the same authors.The hierarchical reconstruction introduced therein is applied to central schemes on overlapping cells and to finite volume schemes on non-staggered grids.This takes a new finite volume approach for approximating non-smooth solutions.A critical step for high-order finite volume schemes is to reconstruct a non-oscillatory high degree polynomial approximation in each cell out of nearby cell averages.In the paper this procedure is accomplished in two steps:first to reconstruct a high degree polynomial in each cell by using e.g.,a central reconstruction,which is easy to do despite the fact that the reconstructed polynomial could be oscillatory;then to apply the hierarchical reconstruction to remove the spurious oscillations while maintaining the high resolution.All numerical computations for systems of conservation laws are performed without characteristic decomposition.In particular,we demonstrate that this new approach can generate essentially non-oscillatory solutions even for 5th-order schemes without characteristic decomposition.

COMPACT FINITE VOLUME SCHEMES AND THEIR APPLICATIONS

A class of Compact Finite Volume Schemes (CFVS) with small support stencils are developed based on a new reconstruction method for cell face variables. The accumulative errors of these schemes for a scalar wave equation are analyzed and compared. The established compact schemes are proved suitable for steady and unsteady flow simulations by several numerical experiments in this paper.

Hyperbolic Conservation Laws on Manifolds.An Error Estimate for Finite Volume Schemes

Following Ben-Artzi and LeFloch, we consider nonlinear hyperbolic conservation laws posed on a Riemannian manifold, and we establish an L1-error estimate for a class of finite volume schemes allowing for the approximation of entropy solutions to the initial value problem. The error in the L1 norm is of order h1/4 at most, where h represents the maximal diameter of elements in the family of geodesic triangulations. The proof relies on a suitable generalization of Cockburn, Coquel, and LeFloch＇s theory which was originally developed in the Euclidian setting. We extend the arguments to curved manifolds, by taking into account the effects to the geometry and overcoming several new technical difficulties.

Monotone Finite Volume Scheme for Three Dimensional Diffusion Equation on Tetrahedral Meshes

We construct a nonlinear monotone finite volume scheme for threedimensional diffusion equation on tetrahedral meshes.Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme,we present a new efficient eliminating method.The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously.The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes,and also show that our scheme appears to be approximate second-order accuracy for solution.

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper,we present a unified finite volume method preserving discrete maximum principle(DMP)for the conjugate heat transfer problems with general interface conditions.We prove the existence of the numerical solution and the DMP-preserving property.Numerical experiments show that the nonlinear iteration numbers of the scheme in[24]increase rapidly when the interfacial coefficients decrease to zero.In contrast,the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero,which reveals that the unified scheme is more robust than the scheme in[24].The accuracy and DMP-preserving property of the scheme are also veri ed in the numerical experiments.

A New Fifth-Order Finite Volume Central WENO Scheme for Hyperbolic Conservation Laws on Staggered Meshes

In this paper,a new fifth-order finite volume central weighted essentially non-oscillatory(CWENO)scheme is proposed for solving hyperbolic conservation laws on staggered meshes.The high-order spatial reconstruction procedure using a convex combination of a fourth degree polynomial with two linear polynomials(in one dimension)or four linear polynomials(in two dimensions)in a traditional WENO fashion and a time discretization method using the natural continuous extension(NCE)of the Runge-Kutta method are applied to design this new fifth-order CWENO scheme.This new finite volume CWENO scheme uses the information defined on the same largest spatial stencil as that of the same order classical CWENO schemes[37,46]with the application of smaller number of unequal-sized spatial stencils.Since the new nonlinear weights are adopted,the new finite volume CWENO scheme could obtain the same order of accuracy and get smaller truncation errors in L1 and L¥norms in smooth regions,and control the spurious oscillations near strong shocks or contact discontinuities.The new CWENO scheme has advantages over the classical CWENO schemes[37,46]on staggered meshes in its simplicity and easy extension to multi-dimensions.Some one-dimensional and two-dimensional benchmark numerical examples are provided to illustrate the good performance of this new fifthorder finite volume CWENO scheme.

A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a Class of Hyperbolic Systems with Source Terms

Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.

High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows

In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.

Computation of the stability derivatives via CFD and the sensitivity equations

The method to calculate the aerodynamic stability derivates of aircrafts by using the sensitivity equations is ex- tended to flows with shock waves in this paper. Using the newly developed second-order cell-centered finite volume scheme on the unstructured-grid, the unsteady Euler equations and sensitivity equations are solved simultaneously in a non-inertial frame of reference, so that the aerodynamic stability derivatives can be calculated for aircrafts with complex geometries. Based on the numerical results, behavior of the aerodynamic sensitivity parameters near the shock wave is discussed. Furthermore, the stability derivatives are analyzed for supersonic and hypersonic flows. The numerical results of the stability derivatives are found in good agree- ment with theoretical results for supersonic flows, and variations of the aerodynamic force and moment predicted by the stability derivatives are very close to those obtained by CFD simulation for both supersonic and hypersonic flows.

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE WENO SCHEMES BASED ON DIFFERENT NUMERICAL FLUXES FOR THE SHALLOW WATER EQUATIONS

In this paper we investigate the performance of the weighted essential non-oscillatory （WENO） methods based on different numerical fluxes, with the objective of obtaining better performance for the shallow water equations by choosing suitable numerical fluxes. We consider six numerical fluxes, i.e., Lax-Friedrichs, local Lax-Friedrichs, Engquist-Osher, Harten-Lax-van Leer, HLLC and the first-order centered fluxes, with the WENO finite volume method and TVD Runge-Kutta time discretization for the shallow water equations. The detailed numerical study is performed for both one-dimensional and two-dimensional shallow water equations by addressing the property, and resolution of discontinuities. issues of CPU cost, accuracy, non-oscillatory

HERMITE WENO SCHEMES WITH STRONG STABILITY PRESERVING MULTI-STEP TEMPORAL DISCRETIZATION METHODS FOR CONSERVATION LAWS

Based on the work of Shu [SIAM J. Sci. Stat. Comput, 9 （1988）, pp.1073-1084], we construct a class of high order multi-step temporal discretization procedure for finite volume Hermite weighted essential non-oscillatory （HWENO） methods to solve hyperbolic conservation laws. The key feature of the multi-step temporal discretization procedure is to use variable time step with strong stability preserving （SSP）. The multi-step tem- poral discretization methods can make full use of computed information with HWENO spatial discretization by holding the former computational values. Extensive numerical experiments are presented to demonstrate that the finite volume HWENO schemes with multi-step diseretization can achieve high order accuracy and maintain non-oscillatory properties near discontinuous region of the solution.

Analysis of the nonlinear scheme preserving the maximum principle for the anisotropic diffusion equation on distorted meshes

In this paper,a nonlinear finite volume scheme preserving the discrete maximum principle for the anisotropic diffusion equation on distorted meshes is described.We prove the coercivity of the scheme under some constraints on the cell deformation and the diffusion coefficient.Numerical results show that the scheme is indeed coercive and satisfies the discrete maximum principle,and the accuracy of this scheme is remarkably better than that of an existing scheme preserving the discrete maximum principle on random triangular meshes.

Mesh Adaptation for Curing the Pathological Behaviors of an Upwind Scheme

In present paper,mesh adaptation is applied for curing the pathological behaviors of the enhanced time-accurate upwind scheme(Loh&Jorgenson,AIAAJ 2016).In the original ETAU(enhanced time-accurate upwind)scheme,a multidimensional dissipation model is required to cure the pathological behaviors.The multi-dimensional dissipation model will increase the global dissipation level reducing numerical resolution.In present work,the metric-based mesh adaptation strategy provides an alternative way to cure the pathological behaviors of the shock capturing.The Hessian matrix of flow variables is applied to construct the metric,which represents the curvature of the physical solution.The adapting operation can well refine the anisotropic meshes at the location with large gradients.The numerical results show that the adaptation of mesh provides a possible way to cure the pathological behaviors of upwind schemes.

Hydrodynamic Regimes,Knudsen Layer,Numerical Schemes:Definition of Boundary Fluxes

We propose a numerical solution to incorporate in the simulation of a system of conservation laws boundary conditions that come from a microscopic modeling in the small mean free path regime.The typical example we discuss is the derivation of the Euler system from the BGK equation.The boundary condition relies on the analysis of boundary layers formation that accounts from the fact that the incoming kinetic flux might be far from the thermodynamic equilibrium.

Conservative Semi-Lagrangian Finite Difference WENO Formulations with Applications to the Vlasov Equation

In this paper,we propose a new conservative semi-Lagrangian(SL)finite difference(FD)WENO scheme for linear advection equations,which can serve as a base scheme for the Vlasov equation by Strang splitting[4].The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume(FV)WENO scheme[3].However,instead of inputting cell averages and approximate the integral form of the equation in a FV scheme,we input point values and approximate the differential form of equation in a FD spirit,yet retaining very high order(fifth order in our experiment)spatial accuracy.The advantage of using point values,rather than cell averages,is to avoid the second order spatial error,due to the shearing in velocity(v)and electrical field(E)over a cell when performing the Strang splitting to the Vlasov equation.As a result,the proposed scheme has very high spatial accuracy,compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson(VP)system.We perform numerical experiments on linear advection,rigid body rotation problem;and on the Landau damping and two-stream instabilities by solving the VP system.For comparison,we also apply(1)the conservative SL FD WENO scheme,proposed in[22]for incompressible advection problem,(2)the conservative SL FD WENO scheme proposed in[21]and(3)the non-conservative version of the SL FD WENO scheme in[3]to the same test problems.The performances of different schemes are compared by the error table,solution resolution of sharp interface,and by tracking the conservation of physical norms,energies and entropies,which should be physically preserved.

On Arbitrary-Lagrangian-Eulerian One-Step WENO Schemes for Stiff Hyperbolic Balance Laws

In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws.High order accuracy in space is obtained with a standard WENO reconstruction algorithm and high order in time is obtained using the local space-time discontinuous Galerkinmethod recently proposed in[20].In the Lagrangian framework considered here,the local space-time DG predictor is based on a weak formulation of the governing PDE on a moving space-time element.For the spacetime basis and test functions we use Lagrange interpolation polynomials defined by tensor-product Gauss-Legendre quadrature points.The moving space-time elements are mapped to a reference element using an isoparametric approach,i.e.the spacetime mapping is defined by the same basis functions as the weak solution of the PDE.We show some computational examples in one space-dimension for non-stiff and for stiff balance laws,in particular for the Euler equations of compressible gas dynamics,for the resistive relativistic MHD equations,and for the relativistic radiation hydrodynamics equations.Numerical convergence results are presented for the stiff case up to sixth order of accuracy in space and time and for the non-stiff case up to eighth order of accuracy in space and time.