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 difference method for singular perturbation problem of hyperbolic-parabolic partial differential equation

In this paper we constructed an exponentially fitted difference scheme for singular perturbation problem of hyperbolic-parabolic partial differential equation. Not only do we take a fitting factor in the equation, but also we put one in the approximation of second initial condition. By means of the asymptotic solution of singular perturbation problem we proved the uniform convergence of this scheme with respect to the small parameter.     

A singular perturbation problem for periodic boundary partial differential equation

In this paper,we consider a singular perturbation elliptic-parabolic partial differentialequation for periodic boundary value problem,and construct a difference scheme.Using themethod of decomposing the singular term from its solution and combining an asymptoticexpansion of the equation,we prove that the scheme constructed by this paper convergesuniformly to the solution of its original problem with O(τ h~2).     

The asymptotic expression of the solution of the Cauchy's problem for a higher order linear ordinary differential equation when the limit equation has singularity

In this paper we consider the asymptotic expression of the solution of the Cauchy's problem for a higher order equation when the limit equation has singularity. In order to construct the asymptotic expression of the solution, the region is divided into three sub-areas. In every small region, the solution of the differential equation is different. Project supported by the National Natural Science Foundation of China     

Asymptotic theory of a wing moving near a solid wall

In this study we use the method of matched asymptotic expansions to obtain an approximate solution of the problem of the nonstationary motion of a lifting surface near a solid wall. The region of flow is provisionally subdivided into characteristic zones, in which, using the appropriate coordinates, we construct asymptotic expansions for the velocity potential, which thereafter coalesce in the regions of common validity. In the first approximation (extremely small heights of flight) the problem reduces to the solution of a Poisson equation in a plane region bounded by the contour of the wing in the horizontal plane with boundary conditions established from the coalescence.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 115–124, November–December, 1977.     

轴对称正交异性圆环壳的齐次完全渐近解

承受轴对称载荷的正交异性圆环壳的静力分析,归结为求解一非齐次二阶复变量方程.当所含参数μ较大时,常采用渐近解法.因方程含一阶转点,所以求全域一致有效且达到薄壳理论精度的完全渐近解较为困难.过去,齐次解只求到一级近似.本文采用广义Airy函数方法,求出了高级近似.这样,轴对称正交异性圆环壳的齐次解第一次有了达到薄壳理论精度的完全的渐近展开.     

A study of upstream-weighted high-order differencing for approximation to flow convection

This paper is concerned with a number of upstream-weighted second- and third-order difference schemes. Also considered are the conventional upwind and central difference schemes for comparison. It commences with a general difference equation which unifies all the given first-, second- and third-order schemes. The various schemes are evaluated through the use of the general equation. The unboundedness and accuracy of the solutions by the difference schemes are assessed via various analyses: examination of the coefficients of the difference equation, Taylor series truncation error analysis, study of the upstream connection to numerical diffusion, single-cell analysis. Finally, the difference schemes are tested on one- and two-dimensional model problems. It is shown that the high-order schemes suffer less from the problem of numerical diffusion than the first-order upwind difference scheme. However, unboundedness cannot be avoided in the solutions by these schemes. Among them the linear upwind difference scheme presents the best compromise between numerical diffusion and solution unboundedness.     

Asymptotic Solutions of a First-Order Differential Equation with Deviating Argument and Slowly Varying Coefficients

The object of this paper is to study the problem of constructing an approximate solution of a first-order weakly nonlinear ordinary differential equation with deviating argument and slowly varying coefficients. On the basis of asymptotic techniques in nonlinear mechanics, we construct an algorithm for the asymptotic integration of the differential equation under consideration.__________Published in Neliniini Kolyvannya, Vol. 7, No. 4, pp. 475–486, October–December, 2004.     

Convergence of Anisotropically Decaying Solutions of a Supercritical Semilinear Heat Equation

We consider the Cauchy problem for a semilinear heat equation with a supercritical power nonlinearity. It is known that the asymptotic behavior of solutions in time is determined by the decay rate of their initial values in space. In particular, if an initial value decays like a radial steady state, then the corresponding solution converges to that steady state. In this paper we consider solutions whose initial values decay in an anisotropic way. We show that each such solution converges to a steady state which is explicitly determined by an average formula. For a proof, we first consider the linearized equation around a singular steady state, and find a self-similar solution with a specific asymptotic behavior. Then we construct suitable comparison functions by using the self-similar solution, and apply our previous results on global stability and quasi-convergence of solutions.     

On the Dirichlet problem for quasi-linear elliptic equation with a small parameter

In this paper we deal with the Dirichlet problem for quasilinear elliptic equation with a small parameter at highest derivatives. In case the characteristics of the degenerated equation are curvilinear and the domain, where the problem is defined, is a bounded convex domain, we offer a method to construct the uniformly valid asymptotic solution of this problem, and prove that the solution of this problem really exists, and being uniquely determined as the small parameter is sufficiently small.     

Exponential tracking and approximation of inertial manifolds for dissipative nonlinear equations

In this paper, we study the long-time behavior of a class of nonlinear dissipative partial differential equations. By means of the Lyapunov-Perron method, we show that the equation has an inertial manifold, provided that certain gap condition in the spectrum of the linear part of the equation is satisfied. We verify that the constructed inertial manifold has the property of exponential tracking (i.e., stability with asymptotic phase, or asymptotic completeness), which makes it a faithful representative to the relevant long-time dynamics of the equation. The second feature of this paper is the introduction of a modified Galerkin approximation for analyzing the original PDE. In an illustrative example (which we believe to be typical), we show that this modified Galerkin approximation yields a smaller error than the standard Galerkin approximation.     

Flux difference splitting for 1D open channel flow equations

An upwind finite difference scheme based on flux difference splitting is presented for the solution of the equations governing unsteady open channel hydraulics. An approximate Jacobian needed for splitting the flux differences is defined that satisfies the conditions required to construct a first-order upwind conservative discretization of the equations. Added limited second-order corrections make the resulting scheme robust and accurate for the computation of all regimes of open channel flow. Some numerical results and comparisons with other classical schemes under exacting conditions are presented.     

A uniformly convergent diffrence scheme for the singular perturbation of a self adjoint elliptic partial differential equation

In this paper, we consider a self adjoint elliptic first boundary value problem with a small parameter affecting the highest derivative.In the paper, we set up a new scheme by the asymptotic analysis method, compare asymptotic behavior between the solution of the difference equation and the solution of the differential equation, and show uniform convergence of the new scheme.     

A monotone finite volume scheme for advection–diffusion equations on distorted meshes

A new monotone finite volume method with second‐order accuracy is presented for the steady‐state advection–diffusion equation. The method uses a nonlinear approximation for both diffusive and advective fluxes that guarantee the positivity of the numerical solution. The approximation of the diffusive flux is based on nonlinear two‐point approximation, and the approximation of the advective flux is based on the second‐order upwind method with proper slope limiter. The second‐order convergence rate for concentration and the monotonicity of the nonlinear finite volume method are verified with numerical experiments. Copyright © 2011 John Wiley & Sons, Ltd.     

对流扩散方程的迎风变换及相应有限差分方法

本文提出所谓迎风变换,将对流扩散方程分解为对流迎风函数和扩散方程,并构造相应的有限差分格式。对流迎风函数以简明的指数解析形式反映对流扩散现象的迎风效应,原则上消除了源于不对称对流算子的困难,能够便利对流扩散方程的数值求解。有限差分格式具有二阶精度和无条件稳定性,算例表明其准确性、收敛速度及对边界层效应的适应能力均明显优于中心差分格式和迎风差分格式。     

Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Singular perturbation of general boundary value problem for higher order quasilinear elliptic equation involving many small parameters

In this paper applying M. I. Visik’s and L. R. Lyuster-nik’s^[1] asymptotic method and principle of fixed point of functional analysis, we study the singular perturbation of general boundary value problem for higher order quasilinear elliptic equation in the case of boundary perturbation combined with equation perturbation. We prove the existence and uniqueness of solution for perturbed problem. We give its asymptotic approximation and estimation of related remainder term.     

Numerical solution of a singularly perturbed elliptichyperbolic partial differential equation on a nonuniform discretization mesh

In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter ε, to the solution of problem (1.1). Numerical results are finally provided.     

Asymptotic solution of singular perturbation problems for the fourth-order elliptic differential equations

In this paper we consider the singularly perturbed boundary value problem for the fourth-order elliptic differential equation, establish the energy estimates of the solutionand its derivatives and construct the formal asymptotic solution by Lyuternik- Vishik ’s method.Finally, by means of the energy estimates we obtain the bound of the remainder of the asymptotic solution.     

极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。