三阶线性常微分方程Sinc方程组的结构预处理方法

三阶线性常微分方程在天文学和流体力学等学科的研究中有着广泛的应用.本文介绍求解三阶线性常微分方程由Sinc方法离散所得到的线性方程组的结构预处理方法.首先, 我们利用Sinc方法对三阶线性常微分方程进行离散,证明了离散解以指数阶收敛到原问题的精确解.针对离散后线性方程组的系数矩阵的特殊结构, 提出了结构化的带状预处理子,并证明了预处理矩阵的特征值位于复平面上的一个矩形区域之内.然后, 我们引入新的变量将三阶线性常微分方程等价地转化为由两个二阶线性常微分方程构成的常微分方程组, 并利用Sinc方法对降阶后的常微分方程组进行离散.离散后线性方程组的系数矩阵是分块2×2的, 且每一块都是Toeplitz矩阵与对角矩阵的组合.为了利用Krylov子空间方法有效地求解离散后的线性方程组,我们给出了块对角预处理子, 并分析了预处理矩阵的性质.最后, 我们对降阶后二阶线性常微分方程组进行了一些比较研究.数值结果证实了Sinc方法能够有效地求解三阶线性常微分方程.     

一种解非线性三阶多点边值问题的数值算法

主要研究了一类带有多点边值条件的非线性三阶微分方程的求解方法.利用迭代技巧和再生核(RKM)理论相结合来求解此类问题,同时给出了一些算例来说明方法的有效性.     

非线性三阶常微分方程的非线性三点边值问题解的存在性

基于上下解方法,在一定条件下,得到了一类带有非线性混合边界条件的三阶常微分方程的非线性三点边值问题解的存在性,作为上述存在性结果的应用,在推论中给出了一类三阶非线性微分方程三点边值问题解的存在性.     

求解二阶刚性微分方程的对角隐式Runge-Kutta-Nystr(o)m方法

在本文中,主要研究二级三阶对角隐式Runge-Kutta-Nystr(o)m(DIRKN)方法关于二阶刚性常微分方程的R-稳定性,P-稳定性以及相延迟性质.我们获得了该方法的R-稳定域,并构造了R-稳定的二级三阶、相延迟阶为四阶的DIRKN方法.P-稳定的二级三阶DIRKN方法被证明是不存在的.我们还构造了相延迟阶为6阶和8阶的二级三阶DIRKN方法,但是这些方法不是R-稳定的.这推广了文献中的单对角隐式Runge-Kutta-Nystr(o)m(SDIRKN)方法的相关结果.     

三阶非线性微分方程周期解的非退化和存在唯一性北大核心CSCD

该文虑了一类三阶线性微分方程x′″(t)+a_(2)x″(t)+a_(1)x’′t)=a_(0)(t)x(t)的非退化性.利用Writinger不等式,给出该方程的非退化条件.再利用三阶线性微分方程的非退化性,证明了三阶非线性微分方程在半线性条件和超线性条件下周期解的存在唯一性.     

一类具三维核共振条件多点边值问题

本文利用Mawhin的重合度理论建立了一类三阶非线性微分方程多点边值问题解的存在性,而后给出一个例子说明该存在性结果的可用性.     

三阶非线性泛函微分方程振动的比较准则

本文利用比较法,与一阶时滞微分方程的振动性相比较,研究一类三阶非线性泛函微分方程的振动性,建立该类方程振动的新比较振动准则.本文结果推广和改进了最近文献中的一些新结果.     

半直线上三阶微分方程的Laguerre-Petrov-Galerkin谱方法

对半无界区域上的三阶方程提出了Laguerre-Petrov-Galerkin谱逼近方法,选取了相同的试探空间和检验空间.通过构造该空间上的基函数,离散问题所对应的线性系统的系数矩阵是半稀疏的.数值算例验证了该方法的有效性和高精度.     

一类三阶偏微分方程的对称约化

主要利用广义条件对称方法研究一类三阶偏微分方程的对称约化问题,首先给出所研究三阶偏微分方程允许的广义条件对称的分类,并根据所允许的对称将偏微分方程的初值问题约化为常微分方程的初值问题并得到其解.     

三阶Airy方程满足两个守恒律的非线性差分格式

近年来,学者们对发展型偏微分方程设计了一种能保持多个守恒律的数值方法,这类方法无论在解的精度还是长时间的数值模拟方面都表现出非常好的性质.将这类思想应用到三阶Airy方程,即三阶散射方程,对其设计了满足两个守恒律的非线性差分格式.该格式不仅计算数值解,同时计算数值能量,并且保证数值解和数值能量同时守恒.从数值结果可以看出,该格式在长时间的数值模拟中具有更好的保结构性质.     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

半直线上三阶方程的Legendre-Laguerre耦合谱元法

针对建立在半直线上的三阶微分方程,提出Legendre-Laguerre耦合谱元法.通过构造满足试探函数空间和检验函数空间的基函数,分解得到的线性系统的系数矩阵是稀疏的,可以有效地进行求解.数值例子验证了方法的有效性和高精度.     

A characteristic nonoverlapping domain decomposition method for multidimensional convection‐diffusion equations

We develop a quasi‐two‐level, coarse‐mesh‐free characteristic nonoverlapping domain decomposition method for unsteady‐state convection‐diffusion partial differential equations in multidimensional spaces. The development of the domain decomposition method is carried out by utilizing an additive Schwarz domain decomposition preconditioner, by using an Eulerian‐Lagrangian method for convection‐diffusion equations and by delicately choosing appropriate interface conditions that fully respect and utilize the hyperbolic nature of the governing equations. Numerical experiments are presented to illustrate the method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A coupled finite element-differential quadrature element method and its accuracy for moving load problem

In this paper a mixed method, which combines the finite element method and the differential quadrature element method (DQEM), is presented for solving the time dependent problems. In this study, the finite element method is first used to discretize the spatial domain. The DQEM is then employed as a step-by-step DQM in time domain to solve the resulting initial value problem. The resulting algebraic equations can be solved by either direct or iterative methods. Two general formulations using the DQM are also presented for solving a system of linear second-order ordinary differential equations in time. The application of the formulation is then shown by solving a sample moving load problem. Numerical results show that the present mixed method is very efficient and reliable.     

A collocation approach for solving linear complex differential equations in rectangular domains

In this paper, a collocation method is presented to find the approximate solution of high‐order linear complex differential equations in rectangular domain. By using collocation points defined in a rectangular domain and the Bessel polynomials, this method transforms the linear complex differential equations into a matrix equation. The matrix equation corresponds to a system of linear equations with the unknown Bessel coefficients. The proposed method gives the analytic solution when the exact solutions are polynomials. Numerical examples are included to demonstrate the validity and applicability of the technique and the comparisons are made with existing results. The results show the efficiency and accuracy of the present work. All of the numerical computations have been performed on a computer using a program written in MATLAB v7.6.0 (R2008a). Copyright © 2012 John Wiley & Sons, Ltd.     

Comparative study on order-reduced methods for linear third-order ordinary differential equations

The linear third-order ordinary differential equation (ODE) can be transformed into a system of two second-order ODEs by introducing a variable replacement, which is different from the common order-reduced approach. We choose the functions p(x) and q(x) in the variable replacement to get different cases of the special order-reduced system for the linear third-order ODE. We analyze the numerical behavior and algebraic properties of the systems of linear equations resulting from the sinc discretizations of these special second-order ODE systems. Then the block-diagonal preconditioner is used to accelerate the convergence of the Krylov subspace iteration methods for solving the discretized system of linear equation. Numerical results show that these order-reduced methods are effective for solving the linear third-order ODEs.     

ON THE HYERS-ULAM STABILITY OF A THIRD-ORDER NONLINEAR DIFFERENTIAL EQUATION

In this paper, the Hyers-Ulam stability of a third-order nonlinear differential equation is investigated. By the integrating method and a Gronwall type inequality, the stability results are obtained in different situations on a bounded domain. Then, the study is extended to nth-order nonlinear differential equations.     

Numerical Solution of Partial Differential Equations in Arbitrary Shaped Domains Using Cartesian Cut-Stencil Finite Difference Method. Part I: Concepts and Fundamentals

A new finite difference (FD) method, referred to as "Cartesian cut-stencil FD", is introduced to obtain the numerical solution of partial differential equations on any arbitrary irregular shaped domain. The 2^nd-order accurate two-dimensional Cartesian cut-stencil FD method utilizes a 5-point stencil and relies on the construction of a unique mapping of each physical stencil, rather than a cell, in any arbitrary domain to a generic uniform computational stencil. The treatment of boundary conditions and quantification of the solution accuracy using the local truncation error are discussed. Numerical solutions of the steady convection-diffusion equation on sample complex domains have been obtained and the results have been compared to exact solutions for manufactured partial differential equations (PDEs) and other numerical solutions.     

On new third-order convergent iterative formulas

In this paper, we present new iteration methods with cubic convergence for solving nonlinear equations. The main advantage of the new methods are free from second derivatives and it permit that the first derivative is zero in some points. Analysis of efficiency shows that the new methods can compete with Newton’s method and the classical third-order methods. Numerical results indicate that the new methods are effective and have definite practical utility.      

Discontinuous Legendre wavelet element method for elliptic partial differential equations

By incorporating the Legendre multiwavelet into the discontinuous Galerkin (DG) method, this paper presents a novel approach for solving Poisson’s equation with Dirichlet boundary, which is known as the discontinuous Legendre multiwavelet element (DLWE) method, derive an adaptive algorithm for the method, and estimate the approximating error of its numerical fluxes. One striking advantage of our method is that the differential operator, boundary conditions and numerical fluxes involved in the elementwise computation can be done with lower time cost. Numerical experiments demonstrate the validity of this method. Furthermore, this paper generalizes the DLWE method to the general elliptic equations defined on a bounded domain and describes the possibilities of constructing optimal adaptive algorithm. The proposed method and its generalizations are also applicable to some other kinds of partial differential equations.