A class of explicit forward time-difference square conservative schemes

In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with adjustable time intervals, which is proved to be effective and remarkably time-saving in numerical tests and applications. Partly supported by the State Major Key Project for Basic Researches of China Worked as a post doctor in Computing Center, Chinese Academy of Sciences, when this paper was submitted     

Four-point explicit difference schemes for the dispersive equation

A class of three-level explicit difference schemes for the dispersive equationu_1=au_(xxx)are established These schemes have higher stability and involve four meshpoints at the middle level.Their local truncation errors are O(τ+h)and stabilityconditions are from|R|≤0.25 to|R|≤10,where|R|=|a|τ/h~3,which is muchbetter than|R|≤0.25.     

Explicit characteristic‐based finite volume element method for level set equation

Accurate modeling of interfacial flows requires a realistic representation of interface topology. To reduce the computational effort from the complexity of the interface topological changes, the level set method is widely used for solving two‐phase flow problems. This paper presents an explicit characteristic‐based finite volume element method for solving the two‐dimensional level set equation. The method is applicable for the case of non‐divergence‐free velocity field. Accuracy and performance of the proposed method are evaluated via test cases with prescribed velocity fields on structured grids. By given a velocity field, the motion of interface in the normal direction and the mean curvature, examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time. Copyright © 2010 John Wiley & Sons, Ltd.     

色散方程的高稳定性两层四点显格式的单点精细积分法

基于单点精细积分的思想,对色散方程Ｕｔ＝ａＵｘｘｘ构造了一类高稳定性的两层四点显式差分格式,其局部截断误差为Ｏ（τ＋ｈ）稳定性条件为│Ｒ│＝│ａτ／ｈ＾３│≤ｆ（β）,对任意正实数β为单调递增函数,它们不仅显著地改善了同类格式的稳定性条件│Ｒ│≤０．２５而且也优于众多三层多点（５点或５点以上）显格式的稳定性条件。     

A class of two-level explicit difference schemes for solving three dimensional heat conduction equation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｉｓｐａｐｅｒｄｅａｌｓｗｉｔｈｔｈｅｉｎｉｔｉａｌ＿ｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｏｆｔｈｒｅｅ＿ｄｉｍｅｎｓｉｏｎａｌｈｅａｔｃｏｎｄｕｃｔｉｏｎｅｑｕａｔｉｏｎｉｎｔｈｅｒｅｇｉｏｎＤ :0≤ｘ,ｙ ,ｚ≤Ｌ ,0 ≤ｔ≤Ｔ ｕ ｔ= 2 ｕ ｘ2 2 ｕ ｙ2 2 ｕ ｚ2 ,ｕ｜ｘ=0 =ｆ1(ｙ,ｚ,ｔ) ,　ｕ｜ｘ=Ｌ =ｆ2 (ｙ ,ｚ,ｔ) ,ｕ｜ｙ=0 =ｇ1(ｚ,ｘ,ｔ) ,　ｕ｜ｙ=Ｌ =ｇ2 (ｚ,ｘ,ｔ) ,ｕ｜ｚ=0 =ｈ1(ｘ ,ｙ ,ｔ) ,　ｕ｜ｚ=Ｌ =ｈ2 (ｘ ,ｙ ,ｔ) ,ｕ｜ｔ=0 =φ(ｘ ,ｙ,ｚ) .(1 )(2 )…     

色散方程的高稳定性三层五点显格式的广义单点精细积分法

基于文（１）中的单点精细积分方法,对色散方程Ｕｔ＝ａＵｘｘｘ提出了一种构造高稳定性三层五点（蛙跳）显格式的广义单点精细积分法,文中格式的局部截断误差为Ｏ（ｘ＾２＋ｈ＾２）,而稳定性条件为｜Ｒ｜≤ｇ（β）（其中ｇ对任意正实数是单调递增函数）,同时类格式中最好的。     

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

IntroductionWith the development of modern industry, various pollutants discharge into the air,rivers, lakes and oceans, which makes the environmental qualities worse and has bad effectson the mankind’s health and the sustained development of industry an…     

The stability of difference schemes of a higher dimensional parabolic equation

This paper proposes a new method to improve the stability condition of differencescheme of a parabolic equation.Necessary and sufficient conditions of the stability of thisnew method are given and proved.Some numerical examples show that this method hassome calculation advantages.     

The Lie-group shooting method for multiple-solutions of Falkner-Skan equation under suction-injection conditions

For the Falkner-Skan equation, including the Blasius equation as a special case, we develop a new numerical technique, transforming the governing equation into a non-linear second-order boundary value problem by a new transformation technique, and then solve it by the Lie-group shooting method. The second-order ordinary differential equation is singular, which is, however, much saving computational cost than the original third-order equation defined in a semi-infinite range. In order to overcome the singularity we consider a perturbed equation. The newly developed Lie-group shooting method allows us to search a missing initial slope at the left-end in a compact space of t∈[0,1], and moreover, the initial slope can be expressed as a closed-form function of r∈(0,1), where the best r is determined by matching the right-end boundary condition. All that makes the new method much superior than the conventional shooting method used in the boundary layer equation under imposed boundary conditions. When the initial slope is available we can apply the fourth-order Runge-Kutta method to calculate the solution, which is highly accurate. The present method is very effective for searching the multiple-solutions under very complex boundary conditions of suction or injection, and also allowing the motion of plate.     

A－HIGH－ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THE EQUATION OF TWO－DIMENSIONAL PARABOLIC TYPE

Ａ－ＨＩＧＨ－ＯＲＤＥＲＡＣＣＵＲＡＣＹＥＸＰＬＩＣＩＴＤＩＦＦＥＲＥＮＣＥＳＣＨＥＭＥＦＯＲＳＯＬＶＩＮＧＴＨＥＥＱＵＡＴＩＯＮＯＦＴＷＯ－ＤＩＭＥＮＳＩＯＮＡＬＰＡＲＡＢＯＬＩＣＴＹＰＥＭａＭｉｎｇｓｈｕ（马明书）（ＲｅｃｅｉｖｅｄＪｕｎｅ２,１...     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

A family of high-order accuracy explicit difference schemes with branching stability for solving 3-D parabolic partial differential equation

ＩｎｔｒｏｄｕｃｔｉｏｎＷｅｏｆｔｅｎｍｅｅｔｔｈｅｐｒｏｂｌｅｍｏｆｓｏｌｖｉｎｇｅｑｕａｔｉｏｎｏｆｐａｒａｂｏｌｉｃｔｙｐｅｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｓｅｅｐａｇｅ ,ｄｉｆｆｕｓｉｏｎ ,ｈｅａｔｃｏｎｄｕｃｔｉｏｎａｎｄｓｏｏｎ .Ｉｎｔｈｅｃａｓｅｏｆ3＿ｄｉｍｅｎｓｉｏｎ ,ｔｈｅｍｏｄｅｌｉｓａｎｉｎｉｔｉａｌａｎｄｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍａｓｆｏｌｌｏｗｓ: ｕ ｔ = 2 ｕ ｘ2 2 ｕ ｙ2 2 ｕ ｚ2 　　　　 (0 <ｘ,ｙ,ｚ<1 ;ｔ>0 ) ,ｕ(ｘ ,ｙ,ｚ,0 ) =φ(ｘ ,ｙ ,ｚ)…     

The alternating segment crank-nicolson method for solving convection-diffusion equation with variable coefficient

IntroductionThenumericalmethodsforsolvingtheconvection_diffusionequationwithvariablecoefficientisappliedwidelyintheproblemsofconductionofheatininhomogeneousmediumandionicdiffusion .Forexample ,thetransmissionofthepaleoheatoftheearth’sinterioronthenumericalsimulationforevolutionaryhistoryofbasinisoneoftheseproblems,hencethestudyoftheparallelnumericalmethodoftheproblemanditsapplicationisveryimportantforthedevelopmentofmodernscienceandtechnology .D .J.EvansandZhangBao_linetal.haveobtainedmanyr…     

The necessary and sufficient condition of uniformly convergent difference schemes for the elliptic-parabolic partial differential equation with a small parameter

This paper studies the necessary and sufficient condition of uniformly convergent difference scheme for the elliptic-parabolic partial differential equation with a small parameter.Communicated by Lin Zong-chi.     

Numerical solution of Poisson equation with wavelet bases of Hermite cubic splines on the interval

A new wavelet-based finite element method is proposed for solving the Poisson equation. The wavelet bases of Hermite cubic splines on the interval are employed as the multi-scale interpolation basis in the finite element analysis. The lifting scheme of the wavelet-based finite element method is discussed in detail. For the orthogonal characteristics of the wavelet bases with respect to the given inner product, the corresponding multi-scale finite element equation can be decoupled across scales, totally or partially, and suited for nesting approximation. Numerical examples indicate that the proposed method has the higher efficiency and precision in solving the Poisson equation.     

Improvement on stability and convergence of A. D. I. schemes

IntroductionAlternatingdirectionimplicit(A.D.I.)schemeswhichwasdiscoveredin1950',hasbecomeoneofthemostimportantmethodsintheapproximationofthesolutionsofparabolicpartialdifferentialequationsinmulti-dimensionalspace.Someofresultsaboutstabilityandconvergencearetooweakandincomplete,we'lltrytoimprovetheminthispaper.Considerinitial-boundaryvalueproblemintwospacevariablesLetΩhbeauniformrectangularmeshofO.h>0isthespacestepinxandydireehon*ProjectsupPOI'tedbytheNahonalNaturalScienceFoundahonofChi…     

A least square extrapolation method for the a posteriori error estimate of the incompressible Navier Stokes problem

A posteriori error estimators are fundamental tools for providing confidence in the numerical computation of PDEs. To date, the main theories of a posteriori estimators have been developed largely in the finite element framework, for either linear elliptic operators or non‐linear PDEs in the absence of disparate length scales. On the other hand, there is a strong interest in using grid refinement combined with Richardson extrapolation to produce CFD solutions with improved accuracy and, therefore, a posteriori error estimates. But in practice, the effective order of a numerical method often depends on space location and is not uniform, rendering the Richardson extrapolation method unreliable. We have recently introduced (Garbey, 13th International Conference on Domain Decomposition, Barcelona, 2002; 379–386; Garbey and Shyy, J. Comput. Phys. 2003; 186 :1–23) a new method which estimates the order of convergence of a computation as the solution of a least square minimization problem on the residual. This method, called least square extrapolation, introduces a framework facilitating multi‐level extrapolation, improves accuracy and provides a posteriori error estimate. This method can accommodate different grid arrangements. The goal of this paper is to investigate the power and limits of this method via incompressible Navier Stokes flow computations. Copyright © 2005 John Wiley & Sons, Ltd.     

ON THE NUMERICAL SOLUTION OF QUASILINEAR WAVE EQUATION WITH STRONG DISSIPATIVE TERM

The numerical solution for a type of quasilinear wave equation is studied. The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved. The error of the difference solution is estimated. The theoretical results are controlled on a numerical example.     

三维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解三维非定常变系数对流扩散方程的一种高精度全隐紧致差分格式,该格式在空间上具有四阶精度,时间具有二阶精度。为了克服传统迭代法在每一个时间步上迭代求解隐格式时收敛速度慢的缺点,采用多重网格加速技术,建立了适用于本文高精度全隐紧致格式的多重网格算法,从而大大加快了迭代收敛速度。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

Numerical schemes with high order of accuracy for the computation of shock waves

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｍｐｕｔａｔｉｏｎｏｆｆｌｏｗｆｉｅｌｄｗｉｔｈｓｈｏｃｋｗａｖｅｓｈａｓｂｅｅｎｔｈｅｓｕｂｊｅｃｔｏｆｒｅｓｅａｒｃｈｆｏｒｍａｎｙｙｅａｒｓ.Ｔｈｅｒｅａｒｅｂａｓｉｃａｌｌｙｔｗｏｃａｔｅｇｏｒｉｅｓｏｆｍｅｔｈｏｄｓ,ｎａｍｅｌｙ,ｓｈｏｃｋｆｉｔｔｉｎｇｍｅｔｈｏｄａｎｄｓｈｏｃｋｃａｐｔｕｒｉｎｇｍｅｔｈｏｄ.Ｔｈｅｆｏｒｍｅｒｄｉｖｉｄｅｓｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｉｎｔｏｓｕｂ＿ｄｏｍａｉｎｓｂｙｔｈｅｓｈｏｃｋｗａｖｅｓ.Ｉｎｅ…