In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws...In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.展开更多
In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materi...In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.展开更多
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-con...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.展开更多
We modify the construction of the third order finite volume WENO scheme on triangular meshes and present a simplified WENO(SWENO)scheme.The novelty of the SWENO scheme is the less complexity and lower computational co...We modify the construction of the third order finite volume WENO scheme on triangular meshes and present a simplified WENO(SWENO)scheme.The novelty of the SWENO scheme is the less complexity and lower computational cost when deciding the smoothest stencil through a simple mechanism.The LU decomposition with iterative refinement is adopted to implement ill-conditioned interpolation matrices and improves the stability of the SWENOscheme efficiently.Besides,a scaling technique is used to circument the growth of condition numbers as mesh refined.However,weak oscillations still appear when the SWENO scheme deals with complex low density equations.In order to guarantee the maximum-principle-preserving(MPP)property,we apply a scaling limiter to the reconstruction polynomial without the loss of accuracy.A novel procedure is designed to prove this property theoretically.Finally,numerical examples for one-and two-dimensional problems are presented to verify the good performance,maximum principle preserving,essentially non-oscillation and high resolution of the proposed scheme.展开更多
We propose a class of up to fourth-order maximum-principle-preserving and mass-conserving schemes for the conservative Allen-Cahn equation equipped with a non-local Lagrange multiplier.Based on the second-order finite...We propose a class of up to fourth-order maximum-principle-preserving and mass-conserving schemes for the conservative Allen-Cahn equation equipped with a non-local Lagrange multiplier.Based on the second-order finite-difference semidiscretization in the spatial direction,the integrating factor Runge-Kutta schemes are applied in the temporal direction.Theoretical analysis indicates that the proposed schemes conserve mass and preserve the maximum principle under reasonable time step-size restriction,which is independent of the space step size.Finally,the theoretical analysis is verified by several numerical examples.展开更多
基金Project supported by the National Natural Science Foundation of China(No.11571366)the Basic Research Foundation of National Numerical Wind Tunnel Project(No.NNW2018-ZT4A08)
文摘In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.
基金Jie Du is supported by the National Natural Science Foundation of China under Grant Number NSFC 11801302Tsinghua University Initiative Scientific Research Program+1 种基金Eric Chung is supported by Hong Kong RGC General Research Fund(Projects 14304217 and 14302018)The third author is supported by the NSF grant DMS-1818467.
文摘In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.
基金This work of J.Lv was partially supported by the Natural Science Foundation of Jilin province grant(No.20200201259JC)the National Natural Science Foundation of China grant(No.12271209)This work of Z.Sheng was partially supported by the National Natural Science Foundation of China grant(No.12071045).
文摘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.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11501570,91530106 and 11571366)Research Fund of NUDT(Grant Nos.JC15-02-02,ZK16-03-53),and the fund from HPCL.
文摘We modify the construction of the third order finite volume WENO scheme on triangular meshes and present a simplified WENO(SWENO)scheme.The novelty of the SWENO scheme is the less complexity and lower computational cost when deciding the smoothest stencil through a simple mechanism.The LU decomposition with iterative refinement is adopted to implement ill-conditioned interpolation matrices and improves the stability of the SWENOscheme efficiently.Besides,a scaling technique is used to circument the growth of condition numbers as mesh refined.However,weak oscillations still appear when the SWENO scheme deals with complex low density equations.In order to guarantee the maximum-principle-preserving(MPP)property,we apply a scaling limiter to the reconstruction polynomial without the loss of accuracy.A novel procedure is designed to prove this property theoretically.Finally,numerical examples for one-and two-dimensional problems are presented to verify the good performance,maximum principle preserving,essentially non-oscillation and high resolution of the proposed scheme.
基金the National Key R&D Program of China(No.2020YFA0709800)the National Key Project(No.GJXM92579)the National Natural Science Foundation of China(No.12071481)。
文摘We propose a class of up to fourth-order maximum-principle-preserving and mass-conserving schemes for the conservative Allen-Cahn equation equipped with a non-local Lagrange multiplier.Based on the second-order finite-difference semidiscretization in the spatial direction,the integrating factor Runge-Kutta schemes are applied in the temporal direction.Theoretical analysis indicates that the proposed schemes conserve mass and preserve the maximum principle under reasonable time step-size restriction,which is independent of the space step size.Finally,the theoretical analysis is verified by several numerical examples.