基于常微分方程边值问题的Chebyshev谱方法

常微分方程边值问题的数值解法有多种,其中较常用的是化边值问题为初值问题解法以及边值问题差分解法.常微分方程边值问题数值解的Chebyshev谱方法是近年来出现的一种新解法.作为应用例子,分别采用Chebyshev谱方法、化边值问题为初值问题解法、以及边值问题差分解法对一类二阶常微分方程边值问题进行数值求解,并对数值解的精确性及计算时间定量地比较,从而说明Chebyshev解法是精度很高的一种快捷解法.     

二阶拟线性奇摄动常微分方程的数值解法

本文讨论二阶拟线性奇摄动常微分方程边值问题的数值解法.首先以一个非线性一阶初值问题近似原问题,然后用迭代法求解该近似问题.最后通过迭代法与古典格式得到一个比较满意的结果.     

同伦摄动法在求解四阶边值问题中的应用

将同伦摄动法用于求解常微分方程四阶边值问题.通过将常微分方程边值问题转化为积分方程组,应用同伦摄动法求得近似解.给出同伦摄动法在两个具体的实例中的应用,并将近似解与精确解进行了比较,验证了同伦摄动法对求解线性、非线性常微分方程边值问题是一种非常有效的方法.     

两相渗流方程组混合始边值问题的适定性

经过一个代数转换,将描述水驱油田压力分布的两相渗流方程组转化成拟线性抛物——椭圆型方程组的混合始边值问题。然后利用拟线性抛物型方程组逼近它们。采用Faedo-Галёркцн法证明了弱解的存在性。在给出了一种能量不等式之后,证明了所论问题解的唯一性及对始边值数据的连续依赖性。还给出了所论方程组的极值原理,从而证明了所论方程组的解适合具体物理系统的要求。 自七十年代以来,人们对孔隙介质中的多相渗流问题予以很大的注意,国内外的学者都有很多的工作。但多数研究者往往只重视通过多相渗流方程研究解(压力等)的性质,而对于该方程组的适定性一般不考虑。关于这方面的情况,可以参考[1]及其所引文献。 由于对多相渗流方程组的适定性缺乏研究,它影响着用最优化的方法研究多相渗流问题。 本文在一定的条件下,将描述水驱油田压力分布的两相渗流方程组转化成拟线性抛物——椭圆型方程组的混合始边值问题。然后用一组拟线性抛物型方程组逼近它们,运用Faedo-Галёркцн方法证明了弱解的存在性。 在给出了对这种方程组电适用的一种能量不等式后,证明了所论问题解的唯一性及对始边值数据的连续依赖性。     

线性代数方程组的列摄动解法

解,当指定精度超过计算机容许限度时也不能自动给出无解信息,传统解法实际上往往把‘好解’和‘坏解’混在一起,一个误差很大,甚至不能利用的‘解’也作为解答形式给出,容易使人困惑。 本文对方程组(1)提出的所谓‘列摄动解法’较准确地解决了误差的定量估计问题。列摄动解法的基本思想是用一定的不等式组代替方程组(1),然后求不等式组的一个解     

一类奇异方程两点边值问题的差分—样条校正解

本文考虑一类奇异方程两点边值问题的差分解和样条数值解法,证明了差分解,样条解分别从两侧逼近精确解,从而得到高精度的差分-样条校正解. 考虑如下形式的奇异边值问题:     

一类奇异方程两点边值问题的差分—样条校正解

本文考虑一类奇异方程两点边值问题的差分解和样条数值解法,证明了差分解,样条解分别从两侧逼近精确解,从而得到高精度的差分-样条校正解. 考虑如下形式的奇异边值问题:     

对偶积分方程与对偶积分方程组及其在固体力学与流体力学中的某些应用

本文首先简要地介绍了文献[1、2]关于对偶积分方程的解,在某些实际问题中,出现的是更为复杂的时偶积分方程组。在文献[1、2]的启发下,我们把这种积分方程组化成复数域上的一般函数方程组,并且由此给出形式解。然后介绍我们用上述两种理论计算得到的固体力学与流体力学中某些混合边值问题的实例,其中出现的对偶积分方程组,用本文建议的方法,得到了精确解。     

L~1空间带弱奇异核的第二类Fredholm积分方程解法探究

针对带有弱奇异核的第二类Fredholm积分方程数值解法问题,介绍了两种方法.一种方法是直接用L~1空间中的离散化方法求其数值解;另一种方法是将弱奇异核通过迭代变为连续核,再用L~1空间中的离散化方法求其数值解,且通过对具体算例作图分析,从而得出直接用L~1空间中离散化方法更好.     

多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.     

Numerical treatment of a mathematical model arising from a model of neuronal variability

In this paper, we describe a numerical approach based on finite difference method to solve a mathematical model arising from a model of neuronal variability. The mathematical modelling of the determination of the expected time for generation of action potentials in nerve cells by random synaptic inputs in dendrites includes a general boundary-value problem for singularly perturbed differential-difference equation with small shifts. In the numerical treatment for such type of boundary-value problems, first we use Taylor approximation to tackle the terms containing small shifts which converts it to a boundary-value problem for singularly perturbed differential equation. A rigorous analysis is carried out to obtain priori estimates on the solution of the problem and its derivatives up to third order. Then a parameter uniform difference scheme is constructed to solve the boundary-value problem so obtained. A parameter uniform error estimate for the numerical scheme so constructed is established. Though the convergence of the difference scheme is almost linear but its beauty is that it converges independently of the singular perturbation parameter, i.e., the numerical scheme converges for each value of the singular perturbation parameter (however small it may be but remains positive). Several test examples are solved to demonstrate the efficiency of the numerical scheme presented in the paper and to show the effect of the small shift on the solution behavior.     

摄动有限元法解一般轴对称壳几何非线性问题

本文在处理几何非线性问题时,利用在变分方程中引入振动过程,得到各级变分摄动方程,并通过有限元法求解.由于有限元法能成功地处理各种复杂边界条件、几何形状的力学问题,摄动法又可将非线性问题转化为线性问题求解.若结合这两种方法的优点,将能够解决大量复杂的非线性力学问题.并能够消除单独使用有限元法或摄动法求解复杂非线性问题所出现的困难. 本文应用摄动有限元法求解了一般轴对称壳的几何非线性问题.     

Numerical solution of Hopf bifurcation problems

We present two numerical methods for the solution of Hopf bifurcation problems involving ordinary differential equations. The first one consists in a discretization of the continuous problem by means of shooting or multiple shooting methods. Thus a finite-dimensional bifurcation problem of special structure is obtained. It may be treated by appropriate iterative algorithms. The second approach transforms the Hopf bifurcation problem into a regular nonlinear boundary value problem of higher dimension which depends on a perturbation parameter ?. It has isolated solutions in the ?-domain of interest, so that conventional discretization methods can be applied. We also consider a concrete Hopf bifurcation problem, a biological feedback inhibition control system. Both methods are applied to it successfully.     

Numerical integration of a class of singular perturbation problems

In this paper, we discuss an approximate method for the numerical integration of a class of linear, singularly perturbed two-point boundary-value problems in ordinary differential equations with a boundary layer on the left end of the underlying interval. This method requires a minimum of problem preparation and can be implemented easily on a computer. We replace the original singular perturbation problem by an approximate first-order differential equation with a small deviating argument. Then, we use the trapezoidal formula to obtain the three-term recurrence relationship. Discrete invariant imbedding algorithm is used to solve a tridiagonal algebraic system. The stability of this algorithm is investigated. The proposed method is iterative on the deviating argument. Several numerical experiments have been included to demonstrate the efficiency of the method.The authors wish to express their sincere thanks to Dr. S. M. Roberts for his comments and valuable suggestions.     

Inverse problems of reconstructing parameters of the Navier-Stokes system

An inverse problem of reconstructing parameters not known a priori of the dynamical system described by the boundary-value problem for the Navier-Stokes system is considered. The reconstruction is based on one piece of admissible information or another about the motion of the dynamical system (solution of the corresponding boundary-value problem). In particular, one of the problems considered is the inverse problem consisting of reconstruction of the a priori unknown right-hand side of the Navier-Stokes system. The right-hand side characterizes the density of exterior mass forces acting on the system. This problem, as well as many other similar problems, is ill-posed. Two methods are proposed for its solution: the statistical method and the dynamical method. These methods use different initial information. In solving the problem by using the statistical method, initial information for the solution is the results of approximate measurement (in one sense or another) of the motion of the dynamical system on a given interval of time. Here, the reconstruction is performed after the corresponding interval of time. For solution of the problem by this method, the concepts and constructions of open-loop control theory are used. In solving the problem by using the dynamical method, initial information for its solution is the results of approximate (in one sense or another) measurements of the current states of the system, which are dynamically obtained by the observer. Here, the reconstruction is dynamically performed during the process. For solution of the problem by the dynamical method, the concepts and constructions of the dynamical regularization method based on positional control theory are used. Also, the author considers various modifications and regularizations of the methods for solution of problems proposed that are based on one piece of a priori information or another about the desired solution and solvability conditions of the problem. __________ Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 26, Nonlinear Dynamics, 2005.     

B-Spline Collocation Method for Nonlinear Singularly-Perturbed Two-Point Boundary-Value Problems

A B-spline collocation method is presented for nonlinear singularly-perturbed boundary-value problems with mixed boundary conditions. The quasilinearization technique is used to linearize the original nonlinear singular perturbation problem into a sequence of linear singular perturbation problems. The B-spline collocation method on piecewise uniform mesh is derived for the linear case and is used to solve each linear singular perturbation problem obtained through quasilinearization. The fitted mesh technique is employed to generate a piecewise uniform mesh, condensed in the neighborhood of the boundary layers. The convergence analysis is given and the method is shown to have second-order uniform convergence. The stability of the B-spline collocation system is discussed. Numerical experiments are conducted to demonstrate the efficiency of the method.     

The method of non-linear time transformation in boundary-value problems of potential theory with moving boundaries for the non-linear wave equation

A new method for solving boundary-value problems for the wave equation [1–3] with moving boundaries is used to obtain a solution of a boundary-value problem with boundary conditions of three types [4].     

Initial-Value Technique for Singularly Perturbed Boundary-Value Problems for Second-Order Ordinary Differential Equations Arising in Chemical Reactor Theory

An initial-value technique is presented for solving singularly perturbed two-point boundary-value problems for linear and semilinear second-order ordinary differential equations arising in chemical reactor theory. In this technique, the required approximate solution is obtained by combining solutions of two terminal-value problems and one initial-value problem which are obtained from the original boundary-value problem through asymptotic expansion procedures. Error estimates for approximate solutions are obtained. Numerical examples are presented to illustrate the present technique.     

地下结构和岩土介质弹塑性耦合分析的摄动半解析法*

本文研究了一个用于物理非线性相互作用分析的有效的数值方法。结构和介质耦合分析的弹塑性问题可用摄动法转化为几个线性问题,然后对相应的线性问题分别用有限条和有限层法分析地下结构和岩土介质以达到简化计算的目的。这种方法用了两次半解析技术——摄动和半解析解函数——将三维非线性耦合问题化为一维的数值问题。此外,本法是半解析法结合解析的摄动法应用于非线性问题的新进展,同时也是近年来发展的摄动数值法的一个分支。     

On the numerical solution of loaded systems of ordinary differential equations with nonseparated multipoint and integral conditions

We propose a numerical method of solving systems of loaded linear nonautonomous ordinary differential equations with nonseparated multipoint and integral conditions. This method is based on the convolution of integral conditions to obtain local conditions. This approach allows one to reduce solving the original problem to solving a Cauchy problem for a system of ordinary differential equations and linear algebraic equations. Numerous computational experiments on several test problems with the formulas and schemes proposed for the numerical solution have been carried out. The results of the experiments show that the approach is reasonably efficient.