Coupling of high order multiplication perturbation method and reduction method for variable coefficient singular perturbation problems

Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.     

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.     

Precise integration method for solving singular perturbation problems

This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.     

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.     

INTERPOLATION PERTURBATION METHOD FOR SOLVING THE BOUNDARY LAYER TYPE PROBLEMS

ＹｕａｎＹｉｗｕ（袁镒吾）（ＲｅｃｅｉｖｅｄＯｃｔ．２，１９９４；ＣｏｍｍｕｎｉｃａｔｅｄｂｙＣｈｉｅｎＷｅｉｚａｎｇ）ＩＮＴＥＲＰＯＬＡＴＩＯＮＰＥＲＴＵＲＢＡＴＩＯＮＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＴＨＥＢＯＵＮＤＡＲＹＬＡＹＥＲＴＹＰＥＰＲＯＢ...     

Singular perturbation for the weakly nonlinear reaction diffusion equation with boundary perturbation

In this paper, a class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation are considered under suitable conditions. Firstly, by dint of the regular perturbation method, the outer solution of the original problem is obtained. Secondly, by using the stretched variable and the expansion theory of power series the initial layer of the solution is constructed. And then, by using the theory of differential inequalities, the asymptotic behavior of the solution for the initial boundary value problems is studied. Finally, using some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation are discussed.     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

ARC-length method for differential equations

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

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

A perturbation boundary integral equation method for analysis of interaction between foundation and nonlinear rock and soil medium

A scheme is developed for analysing the interaction between a foundation and a nonlinear rock and soil medium, in which the foundation is considered as a linear elastic body and a typical boundary integral equation method (BIEM) is employed. On the basis of taking the nonlinear properties of the medium into account, a perturbation BIEM is developed. The fundamental equations for the nonlinear coupling analysis are formulated, and typical problems are solved and discussed by the present method.     

Differential transformation method for studying flow and heat transfer due to stretching sheet embedded in porous medium with variable thickness,variable thermal conductivity,and thermal radiation

This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.     

A uniform high-order method for a singular perturbation problem in conservative form

A uniform high order method is presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems (1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O(h~m+1)accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O(h~m+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.     

Large deflection analysis of rectangular plates by combined perturbation and finite strip method

The perturbation method and finite strip method are combined to solve the largedeflection bending problems of rectangular plates.Perturbation method is used to reducethe nonlinear differential equations into a series of linear differential equations.The finitestrip method is then employed to tackle these linear equations.Some calculation examplesare compared with those got by other methods.     

线性常微分方程灵敏度分析的精细积分法

利用指数矩阵的导数计算来求解一类一阶线性常系数微分方程组对某一设计变量的灵敏度计算问题。对于初值问题,利用指数矩阵的导数,递推得到状态向量的灵敏度;对于线性两点边值问题,通过两点之间的状态向量的导数关系,得到全部初始条件,进而转化为初值问题计算。指数矩阵及其导数阵的高精度计算基于2N类算法,在此基础上可实施灵敏度分析的计算。本文给出了初值和两点边值常微分方程的高精度灵敏度计算方法,计算结果可认为是计算机上的精确解,算例验证了算法的有效性。     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

Singular perturbation of second-order nonlinear system with boundary perturbation

Ｗｅｃｏｎｓｉｄｅｒｉｎｔｈｉｓｐａｐｅｒｔｈｅｓｉｎｇｕｌａｒｐｅｒｔｕｒｂａｔｉｏｎｏｆｓｅｃｏｎｄ＿ｏｒｄｅｒｎｏｎｌｉｎｅａｒｓｙｓｔｅｍｉｎｖｏｌｖｉｎｇｉｎｔｅｒｇｒａｌｏｐｅｒａｔｏｒεｙ″=ｆ(ｔ,ｙ,Ｔｙ,ε)ｙ′ ｇ(ｔ,ｙ,Ｔｙ,ε),(1)ｗｉｔｈｂｏｕｎｄａｒｙｐｅｒｔｕｒｂａｔｉｏｎｙ(ｔ,ε)｜ｔ=φ(ε)=α(ε),ｙ(ｔ,ε)｜ｔ=1 ψ(ε)=β(ε),(2)ｗｈｅｒｅε>0ｉｓａｓｍａｌｌｐａｒａｍｅｔｅｒ,ａｎｄφ(ε),ψ(ε)ａｒｅｂｏｔｈ,ｗｉｔｈｒｅｓｐｅｃｔｔｏε,ｓｕｆｆｉｃｉｅｎｔｌｙｓｍｏ…     

Application of the method of lines to unsteady compressible Euler equations

In the sense of method of lines, numerical solution of the unsteady compressible Euler equations in 1D, 2D and 3D is split into three steps: First, space discretization is performed by the first‐order finite volume method using several approximate Riemann solvers. Second, smoothness and Lipschitz continuity of RHS of the arising system of ordinary dimensional equations (ODEs) is analysed and its solvability is discussed. Finally, the system of ODEs is integrated in time by means of implicit and explicit higher‐order adaptive schemes offered by ODE packages ODEPACK and DDASPK, by a backward Euler scheme based on the linearization of the RHS and by higher‐order explicit Runge–Kutta methods. Time integrators are compared from several points of view, their applicability to various types of problems is discussed, and 1D, 2D and 3D numerical examples are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

比例边界等几何分析方法Ⅰ:波导本征问题

提出比例边界等几何方法 (scaled boundary isogeometric analysis, SBIGA), 并用以求解波导本征值问题. 在比例边界等几何坐标变换的基础上, 利用加权余量法将控制偏微分方程进行离散处理, 半弱化为关于边界控制点变量的二阶常微分方程, 即 TE 波或 TM 波波导的比例边界等几何分析的频域方程以及波导动刚度方程, 同时利用连分式求解波导动刚度矩阵. 通过引入辅助变量进一步得出波导本征方程. 该方法只需在求解域的边界上进行等几何离散, 使问题降低一维, 计算工作量大为节约, 并且由于边界的等几何离散, 使得解的精度更高, 进一步节省求解自由度. 以矩形和 L 形波导的本征问题分析为例, 通过与解析解和其他数值方法比较, 结果表明该方法具有精度高、计算工作量小的优点.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.