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

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

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

A vertex‐centered and positivity‐preserving scheme for anisotropic diffusion equations on general polyhedral meshes

We propose a new nonlinear positivity‐preserving finite volume scheme for anisotropic diffusion problems on general polyhedral meshes with possibly nonplanar faces. The scheme is a vertex‐centered one where the edge‐centered, face‐centered, and cell‐centered unknowns are treated as auxiliary ones that can be computed by simple second‐order and positivity‐preserving interpolation algorithms. Different from most existing positivity‐preserving schemes, the presented scheme is based on a special nonlinear two‐point flux approximation that has a fixed stencil and does not require the convex decomposition of the co‐normal. More interesting is that the flux discretization is actually performed on a fixed tetrahedral subcell of the primary cell, which makes the scheme very easy to be implemented on polyhedral meshes with star‐shaped cells. Moreover, it is suitable for polyhedral meshes with nonplanar faces, and it does not suffer the so‐called numerical heat‐barrier issue. The truncation error is analyzed rigorously, while the Picard method and its Anderson acceleration are used for the solution of the resulting nonlinear system. Numerical experiments are also provided to demonstrate the second‐order accuracy and well positivity of the numerical solution for heterogeneous and anisotropic diffusion problems on severely distorted grids.     

A new difference scheme with high accuracy and absolute stability for solving convection–diffusion equations

In this paper, we use a semi-discrete and a padé approximation method to propose a new difference scheme for solving convection–diffusion problems. The truncation error of the difference scheme is O(h⁴+τ⁵). It is shown through analysis that the scheme is unconditionally stable. Numerical experiments are conducted to test its high accuracy and to compare it with Crank–Nicolson method.     

A parallel finite volume scheme preserving positivity for diffusion equation on distorted meshes

Parallel domain decomposition methods are natural and efficient for solving the implicity schemes of diffusion equations on massive parallel computer systems. A finite volume scheme preserving positivity is essential for getting accurate numerical solutions of diffusion equations and ensuring the numerical solutions with physical meaning. We call their combination as a parallel finite volume scheme preserving positivity, and construct such a scheme for diffusion equation on distorted meshes. The basic procedure of constructing the parallel finite volume scheme is based on the domain decomposition method with the prediction‐correction technique at the interface of subdomains: First, we predict the values on each inner interface of subdomains partitioned by the domain decomposition. Second, we compute the values in each subdomain using a finite volume scheme preserving positivity. Third, we correct the values on each inner interface using the finite volume scheme preserving positivity. The resulting scheme has intrinsic parallelism, and needs only local communication among neighboring processors. Numerical results are presented to show the performance of our schemes, such as accuracy, stability, positivity, and parallel speedup.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2159–2178, 2017     

Calculating the vertex unknowns of nine point scheme on quadrilateral meshes for diffusion equation

In the construction of nine point scheme,both vertex unknowns and cell-centered unknowns are introduced,and the vertex unknowns are usually eliminated by using the interpolation of neighboring cell-centered unknowns,which often leads to lose accuracy.Instead of using interpolation,here we propose a different method of calculating the vertex unknowns of nine point scheme,which are solved independently on a new generated mesh.This new mesh is a Vorono¨i mesh based on the vertexes of primary mesh and some additional points on the interface.The advantage of this method is that it is particularly suitable for solving diffusion problems with discontinuous coeffcients on highly distorted meshes,and it leads to a symmetric positive definite matrix.We prove that the method has first-order convergence on distorted meshes.Numerical experiments show that the method obtains nearly second-order accuracy on distorted meshes.     

A modified Crank-Nicolson scheme with incremental unknowns for convection dominated diffusion equations

A modified Crank-Nicolson scheme based on one-sided difference approximations is proposed for solving time-dependent convection dominated diffusion equations in two-dimensional space. The modified scheme is consistent and unconditionally stable. A priori L² error estimate for the fully discrete modified scheme is derived. With the use of the incremental unknowns preconditioner at each time step, a comparison among several classical numerical schemes has been made and numerical results confirm stability and efficiency of the modified Crank-Nicolson scheme.     

A linear and accurate diffusion scheme respecting the maximum principle on distorted meshes

This Note is devoted to the presentation of a linear diffusion scheme that respects the maximum principle on very distorted meshes. The main idea is to use a classical Finite Volumes scheme and to perform a first integration on the Voronoï mesh based on the centers of the cells. A second integration on the primary mesh is then performed. By construction the scheme preserves the maximum principle. A numerical example is also given. To cite this article: V. Siess, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

A mixed finite volume scheme for anisotropic diffusion problems on any grid

We present a new finite volume scheme for anisotropic heterogeneous diffusion problems on unstructured irregular grids, which simultaneously gives an approximation of the solution and of its gradient. The approximate solution is shown to converge to the continuous one as the size of the mesh tends to 0, and an error estimate is given. An easy implementation method is then proposed, and the efficiency of the scheme is shown on various types of grids and for various diffusion matrices.     

A three-time-level explicit difference scheme of the second order of accuracy for parabolic equations

The difference schemes of Richardson [1] and of Crank-Nicolson [2] are schemes providing second-order approximation. Richardson's three-time-level difference scheme is explicit but unstable and the Crank-Nicolson two-time-level difference scheme is stable but implicit. Explicit numerical methods are preferable for parallel computations. In this paper, an explicit three-time-level difference scheme of the second order of accuracy is constructed for parabolic equations by combining Richardson's scheme with that of Crank-Nicolson. Restrictions on the time step required for the stability of the proposed difference scheme are similar to those that are necessary for the stability of the two-time-level explicit difference scheme, but the former are slightly less onerous.Translated fromMatematicheskie Zametki, Vol. 60, No. 5, pp. 751–759, November, 1996.This research was supported by the Russian Foundation for Basic Research under grant No. 95-01-00489 and by the International Science Foundation under grants No. N8Q300 and No. JBR100.     

A theoretical analysis of a new finite volume scheme for second order hyperbolic equations on general nonconforming multidimensional spatial meshes

We consider the wave equation, on a multidimensional spatial domain. The discretization of the spatial domain is performed using a general class of nonconforming meshes which has been recently studied for stationary anisotropic heterogeneous diffusion problems, see Eymard et al. (IMAJ Numer Anal 30 (2010), 1009–1043). The discretization in time is performed using a uniform mesh. We derive a new implicit finite volume scheme approximating the wave equation and we prove error estimates of the finite volume approximate solution in several norms which allow us to derive error estimates for the approximations of the exact solution and its first derivatives. We prove in particular, when the discrete flux is calculated using a stabilized discrete gradient, the convergence order is begin{align*} h_mathcal{D}end{align*} (resp. k) is the mesh size of the spatial (resp. time) discretization. This estimate is valid under the regularity assumption begin{align*}uin C^3(lbrack 0,Trbrack;C^2(overline{Omega}))end{align*} for the exact solution u. The proof of these error estimates is based essentially on a comparison between the finite volume approximate solution and an auxiliary finite volume approximation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Explicit finite difference schemes with extended stability for advection equations

In this study an explicit central difference approximation of the generalized leap-frog type is applied to the one- and two-dimensional advection equations. The stability of the considered numerical schemes is investigated and the scheme with the largest stable time step is found. For the linear and nonlinear advection equations numerical experiments with different schemes from the considered class are performed in order to evaluate the practical stability of the designed schemes.     

A second-order ADI scheme for three-dimensional parabolic differential equations

A second-order unconditionally stable ADI scheme has been developed for solving three-dimensional parabolic equations. This scheme reduces three-dimensional problems to a succession of one-dimensional problems. Further, the scheme is suitable for simulating fast transient phenomena. Numerical examples show that the scheme gives an accurate solution for the parabolic equation and converges rapidly to the steady state solution. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14:159–168, 1998     

An error estimate for finite volume methods for the Stokes equations on equilateral triangular meshes

We give here an error estimate for a finite volume discretization of the Stokes equations in two space dimensions on equilateral triangular meshes. This work was initiated by an analogous result presented by Alami‐Idrissi and Atounti for general triangular meshes. However, in this latter article, the result is not actually proven. We state here the restricting assumptions (namely equilateral triangles) under which the error estimate holds, using the tools which were introduced by Eymard, Gallouet and Herbin and used by Alami‐Idrissi and Atounti. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

A multisymplectic explicit scheme for the modified regularized long-wave equation

In this paper, we derive a new 10-point multisymplectic scheme for the modified regularized long-wave equation. The new scheme is an explicit scheme in the sense that the third time level does not include nonlinear terms. Numerical results indicate that the new scheme not only provides satisfied numerical solutions, but also preserves three invariants of motion very well.     

四维抛物型方程的高精度分支稳定的显式差分格式

对四维抛物型方程构造了一个高精度显格式,格式的稳定性条件为r=Δt/Δx2=△t/Δy2=△t/△z2=Δt/Δw2<1/2,截断误差阶达到O(Δt2 Δx4),通过数值实验,表明本文理论分析的正确性和文中格式较同类格式的优越性.     

A linearized compact difference scheme for an one-dimensional parabolic inverse problem

The parabolic equation with the control parameter is a class of parabolic inverse problems and is nonlinear. While determining the solution of the problems, we shall determinate some unknown control parameter. These problems play a very important role in many branches of science and engineering. The article is devoted to the following parabolic initial-boundary value problem with the control parameter: ∂u/∂t=∂²u/∂x²+p(t)u+?(x,t),0<x<1,0<t?T$\mathrm{\partial }u/\mathrm{\partial }t={\mathrm{\partial }}^{2}u/\mathrm{\partial }{x}^2+p(t)u+?(x\text{,}t)\text{,}0<x<1\text{,}0<t?T$ satisfying u(x,0)=f(x),0<x<1$u(x\text{,}0)=f(x)\text{,}0<x<1$; u(0,t)=_g0(t)$u(0\text{,}t)=g_0(t)$, u(1,t)=_g1(t)$u(1\text{,}t)=g_1(t)$, u(x^∗,t)=E(t),0?t?T$u(x^{\ast }\text{,}t)=E(t)\text{,}0?t?T$ where ?(x,t),f(x),_g0(t),_g1(t)$?(x\text{,}t)\text{,}f(x)\text{,}{g}_0(t)\text{,}{g}_1(t)$ and E(t)$E(t)$ are known functions, u(x,t)$u(x\text{,}t)$ and p(t)$p(t)$ are unknown functions. A linearized compact difference scheme is constructed. The discretization accuracy of the difference scheme is two order in time and four order in space. The solvability of the difference scheme is proved. Some numerical results and comparisons with the difference scheme given by Dehghan are presented. The numerical results show that the linearized difference scheme of this article improve the accuracy of the space and time direction and shorten computation time largely. The method in this article is also applicable to the two-dimensional inverse problem.     

The explicit solution of a diffusion equation with singularity

We give the explicit solution of the Cauchy problem for the diffusion equation with a singular term:

where . We construct the solution on the basis of a generalization of the Fourier transform. We next show that the solution is expressed by an analytic semigroup, and examine smoothness of and continuity of .

     

Geometric meshes in collocation methods for Volterra integral equations with proportional delays

In this paper we analyse the local superconvergence propertiesof iterated piecewise polynomial collocation solutions for linearsecond-kind Volterra integral equations with (vanishing) proportionaldelays qt (0 < q < 1). It is shown that on suitable geometricmeshes depending on q, collocation at the Gauss points leadsto almost optimal superconvergence at the mesh points. Thiscontrasts with collocation on uniform meshes where the problemregarding the attainable order of local superconvergence remainsopen.