微分、差分方程的机械化方法

介绍了微分与差分方程机械化方法研究若干最新进展.主要结果包括: 微分、差分方程的特征列理论与算法,微分、差分方程系统的分解算法以及微分、差分方程解析解求解算法.     

二阶中立型变时滞差分方程解的振动性

研究了二阶中立型变时滞差分方程Δ2(xn+pxn-l)+qnf(xσ(n))=0解的振动性,获得了该类方程全部非平凡解振动的三个定理.所得结果将二阶中立型差分方程已有的振动性的相应结论推广到了二阶中立型变时滞差分方程.     

一般的高阶线性差分方程的求解研究

本文首先研究了k阶线性差分方程的解的结构,其次给出了齐次方程的基础解系的判定方法和存在条件,并得到相应齐次差分方程通解的具体表达式,最后用实例进行了验证.     

关于一类复差分方程的亚纯解

研究了具有允许的亚纯解的复差分方程的形式以及系数的级与解的级两者的关系,得到了两个结果.将复微分方程中一些结果推广至复差分方程.     

具连续变量的变系数偶数阶差分方程的有界振动

研究时滞差分方程解的性质在理论和应用中是非常重要的.本文借助研究离散变量的差分方程振动性的一般方法,研究了一类具有连续变量的变系数偶数阶中立型差分方程的有界解的振动性,给出了有界解振动的几个充分条件.     

关于广义KPP方程的数值解

广义KPP(Kolmogorov-Petrovskii-Piskunov)方程是一个积分微分方程.为了要研究其数值解,我们首先将该方程转化为一个非线性双曲型方程,然后构造了一个线性化的差分格式,得到了差分格式解的存在唯一性,利用能量不等式证明了差分格式二阶收敛性和关于初值的无条件稳定性,数值结果验证了本文提出的方法.     

Kuramoto-Sivashinsky方程的线性化L_∞-差分方法

利用有限差分方法研究Kuramoto-Sivashinsky方程初边值问题的数值解.首先,给出了二阶线性化隐式差分格式,该格式在每一时间层均为线性方程组.其次,给出差分格式的守恒性和数值解的有界性.第三,证明差分格式在最大模意义下的收敛性.最后,通过数值算例验证差分格式的收敛阶,并数值模拟方程的混沌解.     

非定常的Navier-Stokes方程基于特征正交分解的差分格式

用特征正交分解和奇值分解去研究非定常的Navier-Stokes方程的有限差分格式, 并用有限差分格式计算出的非定常的Navier-Stokes方程瞬时解构成数据集合, 再 用特征正交分解和奇值分解求出这数据集合的元素的最优正交基函数. 结合Galerkin投影方法导出了非定常的Navier-Stokes方程具有较高精确度的低维模型. 并给出了特征正交分解格式解与有限差分格式解的误差分析. 数值例子表明特征正交分解格式解和有限差分格式解的误差与理论分析结果是一致的,从而验证特征正交分解的有效性.     

微分差分方程(t)=f(x(t-1))简单周期解的个数

本文对微分差分方程(t)=f(x(t-1))的简单周期解加以分类,通过证明周期解的对称性引理,给出了微分差分方程恰有 n+1个简单周期解的条件.     

二阶线性矩阵差分方程的解及渐近稳定性

讨论了二阶线性矩阵差分方程AXn+2+BXn+1+CXn=0的解及其渐近稳定性.首先,给出了它的特征方程有解的一个充要条件,然后利用特征方程两个相异的解刻划出该矩阵差分方程的通解,并分析其解的渐近稳定性,最后运用一实例验证了相关结果.     

A high resolution finite difference method for a model of structured susceptible‐infected populations coupled with the environment

We develop a general model describing a structured susceptible‐infected (SI) population coupled with the environment. This model applies to problems arising in ecology, epidemiology, and cell biology. The model consists of a system of quasilinear hyperbolic partial differential equations coupled with a system of nonlinear ordinary differential equations that represents the environment. We develop a second‐order high resolution finite difference scheme to numerically solve the model. Convergence of this scheme to a weak solution with bounded total variation is proved. Numerical simulations are provided to demonstrate the high‐resolution property of the scheme and an application to a multi‐host wildlife disease model is explored.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1420–1458, 2017     

A finite difference approach to the equations of age-dependent population dynamics

We establish the convergence of the finite difference scheme for the nonlinear equations of population dynamics proposed by Guertin and MacCamy. The applicability of the discrete equations to establish qualitative properties of the solution to the continuous problem is also illustrated.     

A mapping finite difference model for infinite elastic media

A simple mapping finite difference model is presented for the solution of boundary-value problems in the theory of time-harmonic elastic vibrations. The finite problem domain is condensed by mapping into a smaller finite domain using a suitable coordinate transformation. The field equations and the boundary conditions are also appropriately transformed. The radiation condition at infinity is satisfied through a change of the dependent variable. Finite difference forms of the transformed equations are then solved in the mapped domain, subject to the transformed boundary conditions.     

A compact local one‐dimensional scheme for solving a 3D N‐carrier system with Neumann boundary conditions

We consider a mathematical model for thermal analysis in a 3D N‐carrier system with Neumann boundary conditions, which extends the concept of the well‐known parabolic two‐step model for micro heat transfer. To solve numerically the complex system, we first reduce 3D equations in the model to a succession of 1D equations by using the local one‐dimensional (LOD) method. The obtained 1D equations are then solved using a fourth‐order compact finite difference scheme for the interior points and a second‐order combined compact finite difference scheme for the points next to the boundary, so that the Neumann boundary condition can be applied directly without discretizing. By using matrix analysis, the compact LOD scheme is shown to be unconditionally stable. The accuracy of the solution is tested using two numerical examples. Results show that the solutions obtained by the compact LOD finite difference scheme are more accurate than those obtained by a Crank‐Nicholson LOD scheme, and the convergence rate with respect to spatial variables is about 2.6. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

常系数齐次线性差分方程组的一个新解法

常系数齐次线性差分方程组的求解方法,已有作者讨论过,但都没有给出一个比较简便的计算方法.本文将给出一个十分简明而有效的常系数齐次线性差分方程组的新求解方法.     

A model for peak formation in the two-phase equations

We present a hyperbolic-elliptic model problem related to the equations of two-phase fluid flow. The model problem is solved numerically, and properties of its solution are presented. The model equation is well-posed when linearized around a constant state, but there is a strong focusing effect, and very large solutions exist at certain times. We prove that the model problem has a smooth solution for bounded times.

     

A pair of difference differential equations of Euler-Cauchy type

We study two classes of linear difference differential equations analogous to Euler-Cauchy ordinary differential equations, but in which multiple arguments are shifted forward or backward by fixed amounts. Special cases of these equations have arisen in diverse branches of number theory and combinatorics. They are also of use in linear control theory. Here, we study these equations in a general setting. Building on previous work going back to de Bruijn, we show how adjoint equations arise naturally in the problem of uniqueness of solutions. Exploiting the adjoint relationship in a new way leads to a significant strengthening of previous uniqueness results. Specifically, we prove here that the general Euler-Cauchy difference differential equation with advanced arguments has a unique solution (up to a multiplicative constant) in the class of functions bounded by an exponential function on the positive real line. For the closely related class of equations with retarded arguments, we focus on a corresponding class of solutions, locating and classifying the points of discontinuity. We also provide an explicit asymptotic expansion at infinity.

     

A level set approach for computing discontinuous solutions of Hamilton-Jacobi equations

We introduce two types of finite difference methods to compute the L-solution and the proper viscosity solution recently proposed by the second author for semi-discontinuous solutions to a class of Hamilton-Jacobi equations. By regarding the graph of the solution as the zero level curve of a continuous function in one dimension higher, we can treat the corresponding level set equation using the viscosity theory introduced by Crandall and Lions. However, we need to pay special attention both analytically and numerically to prevent the zero level curve from overturning so that it can be interpreted as the graph of a function. We demonstrate our Lax-Friedrichs type numerical methods for computing the L-solution using its original level set formulation. In addition, we couple our numerical methods with a singular diffusive term which is essential to computing solutions to a more general class of HJ equations that includes conservation laws. With this singular viscosity, our numerical methods do not require the divergence structure of equations and do apply to more general equations developing shocks other than conservation laws. These numerical methods are generalized to higher order accuracy using weighted ENO local Lax-Friedrichs methods as developed recently by Jiang and Peng. We verify that our numerical solutions approximate the proper viscosity solutions obtained by the second author in a recent Hokkaido University preprint. Finally, since the solution of scalar conservation law equations can be constructed using existing numerical techniques, we use it to verify that our numerical solution approximates the entropy solution.

     

On Stability of Solutions to One Class of Nonlinear Difference Systems

We consider some class of systems of nonlinear ordinary differential equations. We adjust the difference schemes corresponding to the equations under study in order to guarantee agreement between differential and difference systems in the sense of stability of the zero solution. We obtain conditions under which perturbations do not violate the asymptotic stability of solutions to difference systems.     

A conservative pressure-correction scheme for transient simulations of reacting flows

An algorithm is presented for numerical simulations of time-dependent low Mach number variable density flows with an arbitrary amount of scalar transport equations and a complex equation of state. The pressure-correction type algorithm is based on a segregated solution formalism. It is conservative and guarantees stable results, regardless of the difference in density between neighboring cells. Furthermore, states are predicted which exactly match the equation of state. In the one-dimensional example, considering non-premixed flames, a simplified flamesheet model is used to describe the combustion of fuel and oxidizer. We demonstrate that the predicted states exactly correspond to the equation of state. We illustrate the accuracy improvement due to higher order formulation and demonstrate grid convergence.