We develop mesh conditions for linear finite volume element approximations of anisotropic diffusionconvectionreaction problems to satisfy the discrete maximum principle.We obtain the sufficient conditions to gurantee ...We develop mesh conditions for linear finite volume element approximations of anisotropic diffusionconvectionreaction problems to satisfy the discrete maximum principle.We obtain the sufficient conditions to gurantee the both upper and lower bounds of the numerical solution when each angle of arbitrary triangle is O(∥q∥_∞h+∥g∥_∞h~2)-acute and h is small enough,where h denotes the mesh size,q and g are coefficients of the convection and reaction terms,respectively.To deal with the convection-dominated problems,we use the upwind triangle technique.For such scheme,the mesh condition can be sharper to O(∥g∥_∞h~2)-acute.Some numerical examples are presented to demonstrate the theoretical results.展开更多
In this paper some new parallel difference schemes with interface extrapolation terms for a quasi-linear parabolic system of equations are constructed. Two types of time extrapolations are proposed to give the interfa...In this paper some new parallel difference schemes with interface extrapolation terms for a quasi-linear parabolic system of equations are constructed. Two types of time extrapolations are proposed to give the interface values on the interface of sub-domains or the values adjacent to the interface points, so that the unconditional stable parallel schemes with the second accuracy are formed. Without assuming heuristically that the original boundary value problem has the unique smooth vector solution, the existence and uniqueness of the discrete vector solutions of the parallel difference schemes constructed are proved. Moreover the unconditional stability of the parallel difference schemes is justified in the sense of the continuous dependence of the discrete vector solution of the schemes on the discrete known data of the original problems in the discrete W2(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the parallel difference schemes with interface extrapolation terms to the unique generalized solution of the original quasi-linear parabolic problem is proved. Numerical results are presented to show the good performance of the parallel schemes, including the unconditional stability, the second accuracy and the high parallelism.展开更多
We study the asymptotic-preserving fully discrete schemes for nonequilibrium radiation diffusion problem in spherical and cylindrical symmetric geometry.The research is based on two-temperature models with Larsen’s f...We study the asymptotic-preserving fully discrete schemes for nonequilibrium radiation diffusion problem in spherical and cylindrical symmetric geometry.The research is based on two-temperature models with Larsen’s flux-limited diffusion operators.Finite volume spatially discrete schemes are developed to circumvent the singularity at the origin and the polar axis and assure local conservation.Asymmetric second order accurate spatial approximation is utilized instead of the traditional first order one for boundary flux-limiters to consummate the schemes with higher order global consistency errors.The harmonic average approach in spherical geometry is analyzed,and its second order accuracy is demonstrated.By formal analysis,we prove these schemes and their corresponding fully discrete schemes with implicitly balanced and linearly implicit time evolutions have first order asymptoticpreserving properties.By designing associated manufactured solutions and reference solutions,we verify the desired performance of the fully discrete schemes with numerical tests,which illustrates quantitatively they are first order asymptotic-preserving and basically second order accurate,hence competent for simulations of both equilibrium and non-equilibrium radiation diffusion problems.展开更多
基金the National Natural Science Foundation of China(No.11301033,11971069 and 12271209)the Natural Science Foundation of Jilin Province(No.20200201259JC)Jilin Province Science and Technology Plan Development Project(No.20210201078GX)。
文摘We develop mesh conditions for linear finite volume element approximations of anisotropic diffusionconvectionreaction problems to satisfy the discrete maximum principle.We obtain the sufficient conditions to gurantee the both upper and lower bounds of the numerical solution when each angle of arbitrary triangle is O(∥q∥_∞h+∥g∥_∞h~2)-acute and h is small enough,where h denotes the mesh size,q and g are coefficients of the convection and reaction terms,respectively.To deal with the convection-dominated problems,we use the upwind triangle technique.For such scheme,the mesh condition can be sharper to O(∥g∥_∞h~2)-acute.Some numerical examples are presented to demonstrate the theoretical results.
基金This work was supported by the Special Funds for Major State Basic Research Projects (Grant No.2005CB321703)the National Natural Science Foundation of China (Grant Nos. 10476002, 60533020)the Science Foundation of CAEP (Grant No. 20060649)
文摘In this paper some new parallel difference schemes with interface extrapolation terms for a quasi-linear parabolic system of equations are constructed. Two types of time extrapolations are proposed to give the interface values on the interface of sub-domains or the values adjacent to the interface points, so that the unconditional stable parallel schemes with the second accuracy are formed. Without assuming heuristically that the original boundary value problem has the unique smooth vector solution, the existence and uniqueness of the discrete vector solutions of the parallel difference schemes constructed are proved. Moreover the unconditional stability of the parallel difference schemes is justified in the sense of the continuous dependence of the discrete vector solution of the schemes on the discrete known data of the original problems in the discrete W2(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the parallel difference schemes with interface extrapolation terms to the unique generalized solution of the original quasi-linear parabolic problem is proved. Numerical results are presented to show the good performance of the parallel schemes, including the unconditional stability, the second accuracy and the high parallelism.
基金The authors are very grateful to the editors and the anonymous referees for helpful suggestions to enhance the paper.This work is supported by the National Natural Science Foundation of China(11271054,11471048,11571048,U1630249)the Science Foundation of CAEP(2014A0202010)the Science Challenge Project(No.JCKY2016212A502)and the Foundation of LCP.
文摘We study the asymptotic-preserving fully discrete schemes for nonequilibrium radiation diffusion problem in spherical and cylindrical symmetric geometry.The research is based on two-temperature models with Larsen’s flux-limited diffusion operators.Finite volume spatially discrete schemes are developed to circumvent the singularity at the origin and the polar axis and assure local conservation.Asymmetric second order accurate spatial approximation is utilized instead of the traditional first order one for boundary flux-limiters to consummate the schemes with higher order global consistency errors.The harmonic average approach in spherical geometry is analyzed,and its second order accuracy is demonstrated.By formal analysis,we prove these schemes and their corresponding fully discrete schemes with implicitly balanced and linearly implicit time evolutions have first order asymptoticpreserving properties.By designing associated manufactured solutions and reference solutions,we verify the desired performance of the fully discrete schemes with numerical tests,which illustrates quantitatively they are first order asymptotic-preserving and basically second order accurate,hence competent for simulations of both equilibrium and non-equilibrium radiation diffusion problems.