High order positivity-preserving finite volume WENO schemes for a hierarchical size-structured population model

In this paper we develop high order positivity-preserving finite volume weighted essentially non-oscillatory (WENO) schemes for solving a hierarchical size-structured population model with nonlinear growth, mortality and reproduction rates. We carefully treat the technical complications in boundary conditions and global integration terms to ensure high order accuracy and the positivity-preserving property. Comparing with the previous high order difference WENO scheme for this model, the positivity-preserving finite volume WENO scheme has a comparable computational cost and accuracy, with the added advantages of being positivity-preserving and having L¹ stability. Numerical examples, including that of the evolution of the population of Gambusia affinis, are presented to illustrate the good performance of the scheme.     

High order positivity-preserving finite volume WENO schemes for a hierarchical size-structured population model

In this paper we develop high order positivity-preserving finite volume weighted essentially non-oscillatory (WENO) schemes for solving a hierarchical size-structured population model with nonlinear growth, mortality and reproduction rates. We carefully treat the technical complications in boundary conditions and global integration terms to ensure high order accuracy and the positivity-preserving property. Comparing with the previous high order difference WENO scheme for this model, the positivity-preserving finite volume WENO scheme has a comparable computational cost and accuracy, with the added advantages of being positivity-preserving and having $L^1$ stability. Numerical examples, including that of the evolution of the population of Gambusia affinis, are presented to illustrate the good performance of the scheme.     

High order Hybrid central—WENO finite difference scheme for conservation laws

In this article we present a high resolution hybrid central finite difference—WENO scheme for the solution of conservation laws, in particular, those related to shock–turbulence interaction problems. A sixth order central finite difference scheme is conjugated with a fifth order weighted essentially non-oscillatory WENO scheme in a grid-based adaptive way. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme while the smooth regions are computed with the more efficient and accurate central finite difference scheme. The application of high order filtering to mitigate the dispersion error of central finite difference schemes is also discussed. Numerical experiments with the 1D compressible Euler equations are shown.     

An adaptive mesh method for 1D hyperbolic conservation laws

An adaptive method is developed for solving one-dimensional systems of hyperbolic conservation laws, which combines the rezoning approach with the finite volume weighted essentially non-oscillatory (WENO) scheme. An a posteriori error estimate, used to equidistribute the mesh, is obtained from the differences between respective numerical solutions of 5th-order WENO (WENO5) and 3rd-order ENO (ENO3) schemes. The number of grids can be adaptively readjusted based on the solution structure. For higher efficiency, mesh readjustment is performed every few time steps rather than every time step. In addition, a high order conservative interpolation is used to compute the physical solutions on the new mesh from old mesh based on the finite volume ENO reconstruction. Extensive examples suggest that this adaptive method exhibits more accurate resolution of discontinuities for a similar level of computational time comparing with that on a uniform mesh.     

Hybrid WENO schemes with different indicators on curvilinear grids

In {J. Comput. Phys. 229 (2010) 8105-8129}, we studied hybrid weighted essentially non-oscillatory (WENO) schemes with different indicators for hyperbolic conservation laws on uniform grids for Cartesian domains. In this paper, we extend the schemes to solve two-dimensional systems of hyperbolic conservation laws on curvilinear grids for non-Cartesian domains. Our goal is to obtain similar advantageous properties as those of the hybrid WENO schemes on uniform grids for Cartesian domains. Extensive numerical results strongly support that the hybrid WENO schemes with discontinuity indicators on curvilinear grids can also save considerably on computational cost in contrast to the pure WENO schemes. They also maintain the essentially non-oscillatory property for general solutions with discontinuities and keep the sharp shock transition.     

Binary weighted essentially non-oscillatory (BWENO) approximation

Weighted essentially non-oscillatory (WENO) schemes have been mainly used for solving hyperbolic partial differential equations (PDEs). Such schemes are capable of high order approximation in smooth regions and non-oscillatory sharp resolution of discontinuities. The base of the WENO schemes is a non-oscillatory WENO approximation procedure, which is not necessarily related to PDEs. The typical WENO procedures are WENO interpolation and WENO reconstruction. The WENO algorithm has gained much popularity but the basic idea of approximation did not change much over the years. In this paper, we first briefly review the idea of WENO interpolation and propose a modification of the basic algorithm. New approximation should improve basic characteristics of the approximation and provide a more flexible framework for future applications. New WENO procedure involves a binary tree weighted construction that is based on key ideas of WENO algorithm and we refer to it as the binary weighted essentially non-oscillatory (BWENO) approximation. New algorithm comes in a rational and a polynomial version. Furthermore, we describe the WENO reconstruction procedure, which is usually involved in the numerical schemes for hyperbolic PDEs, and propose the new reconstruction procedure based on the described BWENO interpolation. The obtained numerical results show that the newly proposed procedures perform very well on the considered test examples.     

A simple weighted essentially non-oscillatory limiter for the correction procedure via reconstruction (CPR) framework on unstructured meshes

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes on two-dimensional unstructured triangular meshes [39], to the correction procedure via reconstruction (CPR) framework for solving nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes with straight or curved edges. This is an extension of our earlier work [4] in which the WENO limiter was designed for the CPR framework on regular meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper uses information only from the target cell and its immediate neighbors. Hence, it is particularly simple to implement and will not harm the conservativeness and compactness of the CPR framework. Since the CPR framework with this WENO limiter does not in general satisfy the positivity preserving property, we also extend the positivity-preserving limiters [36], [33] to the CPR framework. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of this procedure.     

A Vertex-Centered and Positivity-Preserving Finite Volume Scheme for Two-Dimensional Three-Temperature Radiation Diffusion Equations on General Polygonal Meshes

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multimaterial system and explain the exchange of energy among electrons, ions and photons. In this paper, we suggest a new positivity-preserving finite volume scheme for 2-D 3-T radiation diffusion equations on general polygonal meshes. The vertex unknowns are treated as primary ones for which the finite volume equations are constructed. The edge-midpoint and cell-centered unknowns are used as auxiliary ones and interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. By comparison, most existing positivity-preserving finite volume schemes are cell-centered and based on the convex decomposition of the co-normal. Here, the co-normal decomposition is not convex in general, leading to a fixed stencil of the flux approximation and avoiding a certain search algorithm on complex grids. Moreover, the new scheme effectively alleviates the numerical heat-barrier issue suffered by most existing cell-centered or hybrid schemes in solving strongly nonlinear radiation diffusion equations. Numerical experiments demonstrate the second-order accuracy and the positivity of the solution on various distorted grids. For the problem without analytic solution, the contours of the numerical solutions obtained by our scheme on distorted meshes accord with those on smooth quadrilateral meshes.     

A class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes

In this paper, we developed a class of the fourth order accurate finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes based on the work (Computers & Fluids, 34: 642–663 (2005)) by Qiu and Shu, with Total Variation Diminishing Runge-Kutta time discretization method for the two-dimensional hyperbolic conservation laws. The key idea of HWENO is to evolve both with the solution and its derivative, which allows for using Hermite interpolation in the reconstruction phase, resulting in a more compact stencil at the expense of the additional work. The main difference between this work and the formal one is the procedure to reconstruct the derivative terms. Comparing with the original HWENO schemes of Qiu and Shu, one major advantage of new HWENO schemes is its robust in computation of problem with strong shocks. Extensive numerical experiments are performed to illustrate the capability of the method. Corresponding author This work was partially supported by the National Natural Science Foundation of China (Grant No. 10671097), the European project ADIGMA on the development of innovative solution algorithms for aerodynamic simulations, Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry and the Natural Science Foundation of Jiangsu Province (Grant No. BK2006511)     

A minimum entropy principle of high order schemes for gas dynamics equations

The entropy solutions of the compressible Euler equations satisfy a minimum principle for the specific entropy (Tadmor in Appl Numer Math 2:211–219, 1986). First order schemes such as Godunov-type and Lax-Friedrichs schemes and the second order kinetic schemes (Khobalatte and Perthame in Math Comput 62:119–131, 1994) also satisfy a discrete minimum entropy principle. In this paper, we show an extension of the positivity-preserving high order schemes for the compressible Euler equations in Zhang and Shu (J Comput Phys 229:8918–8934, 2010) and Zhang et?al. (J Scientific Comput, in press), to enforce the minimum entropy principle for high order finite volume and discontinuous Galerkin (DG) schemes.     

A high order central-upwind scheme for hyperbolic conservation laws

A high order central-upwind scheme for approximating hyperbolic conservation laws is proposed. This construction is based on the evaluation of the local propagation speeds of the discontinuities and Peer's fourth order non-oscillatory reconstruction. The presented scheme shares the simplicity of central schemes, namely no Riemann solvers are involved. Furthermore, it avoids alternating between two staggered grids, which is particularly a challenge for problems which involve complex geometries and boundary conditions. Numerical experiments demonstrate the high resolution and non-oscillatory properties of our scheme.     

可压缩流体力学高精度拉格朗日格式及其保正性质

本文对可压缩流体力学高精度拉格朗日格式及其保正性质近年来的发展给出回顾与综述.文中分别介绍了一维、二维可压缩流体力学方程中心型拉格朗日格式的设计步骤,回顾了高精度拉格朗日格式以及高精度保正拉格朗日格式的研究进展.     

Volume-preserving algorithms for source-free dynamical systems

Summary. In this paper, we first expound why the volume-preserving algorithms are proper for numerically solving source-free systems and then prove all the conventional methods are not volume-preserving. Secondly, we give a general method of constructing volume-preserving difference schemes for source-free systems on the basis of decomposing a source-free vector field as a finite sum of essentially 2-dimensional Hamiltonian fields and of composing the corresponding essentially symplectic schemes into a volume-preserving one. Lastly, we make some special discussions for so-called separable source-free systems for which arbitrarily high order explicit revertible volume-preserving schemes can be constructed. Received March 22, 1994 / Revised version received January 25, 1995     

基于ENO格式的三阶修正系数格式

不增加基点,仅摄动二阶ENO格式的系数(简记为MCENO),得到一类求解双曲型守恒律方程的三阶MCENO格式．由MCENO格式的构造过程可以看出,MCENO格式保留了ENO格式的许多性质,例如本质无振荡性、TVB性质等,且能提高一阶精度．进一步,利用MCENO格式模拟二维Rayleigh-Taylor(RT)不稳定性和Lax激波管的数值求解问题．数值结果表明,t=2.0时,MCENO格式的密度曲线处于三阶WENO格式和五阶WENO格式之间,是一个高效高精度格式．值得注意的是,三阶MCENO格式,三阶WENO格式和五阶WENO格式的CPU时间之比为0.62:1:2.19．表明相对于原始ENO格式,MCENO格式在光滑区域有较高精度,能提高格式精度．     

Low dissipative entropy stable schemes using third order WENO and TVD reconstructions

A low dissipative framework is given to construct high order entropy stable flux by addition of suitable numerical diffusion operator into entropy conservative flux. The framework is robust in the sense that it allows the use of high order reconstructions which satisfy the sign property only across the discontinuities. The third order weighted essentially non-oscillatory (WENO) interpolations and high order total variation diminishing (TVD) reconstructions are shown to satisfy the sign property across discontinuities. Third order accurate entropy stable schemes are constructed by using third order WENO and high order TVD reconstructions procedures in the diffusion operator. These schemes are efficient and less diffusive since the diffusion is actuated only in the sign stability region of the used reconstruction which includes discontinuities. Numerical results with constructed schemes for various test problems are given which show the third order accuracy and less dissipative nature of the schemes.     

多边形网格上扩散方程新的单调格式

 本文在星形多边形网格上, 构造了扩散方程新的单调有限体积格式.该格式与现有的基于非线性两点流的单调格式的主要区别是, 在网格边的法向流离散模板中包含当前边上的点, 在推导离散法向流的表达式时采用了定义于当前边上的辅助未知量, 这样既可适应网格几何大变形, 同时又兼顾了当前网格边上物理量的变化. 在光滑解情形证明了离散法向流的相容性.对于具有强各向异性、非均匀张量扩散系数的扩散方程, 证明了新格式是单调的, 即格式可以保持解析解的正性. 数值结果表明在扭曲网格上, 所构造的格式是局部守恒和保正的, 对光滑解有高于一阶的精度, 并且, 针对非平衡辐射限流扩散问题, 数值结果验证了新格式在计算效率和守恒精度上优于九点格式.     

On the positivity of linear weights in WENO approximations

High order accurate weighted essentially non-oscillatory （WENO） schemes have been used extensively in numerical solutions of hyperbolic partial differential equations and other convection dominated problems. However the WENO procedure can not be applied directly to obtain a stable scheme when negative linear weights are present. In this paper, we first briefly review the WENO framework and the role of linear weights, and then present a detailed study on the positivity of linear weights in a few typical WENO procedures, including WENO interpolation, WENO reconstruction and WENO approximation to first and second derivatives, and WENO integration. Explicit formulae for the linear weights are also given for these WENO procedures. The results of this paper should be useful for future design of WENO schemes involving interpolation, reconstruction, approximation to first and second derivatives, and integration procedures.     

High accurate NRK and MWENO scheme for nonlinear degenerate parabolic PDEs

In this paper, we propose a stable high accurate hybrid scheme based on nonstandard Runge–Kutta (NRK) and modified weighted essentially non-oscillatory (MWENO) techniques for nonlinear degenerate parabolic partial differential equations. The necessary stability condition for the combination of a Runge–Kutta and MWENO scheme is given. The stability condition provides a renormalization function such that mixture of explicit NRK and MWENO scheme is unconditionally stable. Novel scheme recovers the sixth order convergent at points of inflection and prevents the appearance of spurious solutions close to discontinuities. The good performance of this scheme is illustrated through five examples. Numerical results are presented.     

A virtual element method-based positivity-preserving conservative scheme for convection–diffusion problems on polygonal meshes

In this paper, we propose a positivity-preserving conservative scheme based on the virtual element method (VEM) to solve convection–diffusion problems on general meshes. As an extension of finite element methods to general polygonal elements, the VEM has many advantages such as substantial mathematical foundations, simplicity in implementation. However, it is neither positivity-preserving nor locally conservative. The purpose of this article is to develop a new scheme, which has the same accuracy as the VEM and preserves the positivity of the numerical solution and local conservation on primary grids. The first step is to calculate the cell-vertex values by the lowest-order VEM. Then, the nonlinear two-point flux approximations are utilized to obtain the nonnegativity of cell-centered values and the local conservation property. The new scheme inherits both advantages of the VEM and the nonlinear two-point flux approximations. Numerical results show that the new scheme can reach the optimal convergence order of the virtual element theory, that is, the second-order accuracy for the solution and the first-order accuracy for its gradient. Moreover, the obtained cell-centered values are nonnegative, which demonstrates the positivity-preserving property of our new scheme.     

Convergence of higher order finite volume schemes on irregular grids

We prove convergence to the entropy solution of a general class of higher order finite volume schemes on unstructured, irregular grids for multidimensional scalar conservation laws. Such grids allow for cells to become flat in the limit. We derive a new entropy inequality for higher order schemes built on Godunov's numerical flux. Our result implies convergence of suitably modified versions of MUSCL-type finite volume schemes, ENO schemes and the discontinuous Galerkin finite element method. Supported by Deutsche Forschungsgemeinschaft, SFB256.