长水波近似方程组的动力学行为和辛几何算法

长水波近似方程组作为一类重要的非线性方程有着许多广泛的应用前景,特别是在浅水非线性色散波的研究中具有重要意义.给出了长水波近似方程组的动力学行为,并基于Hamilton空间体系的多辛理论研究了长水波近似方程组的数值解法,讨论了利用Preissmann方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,格式满足多辛守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

解含有混合导数的多维抛物型方程的两个显格式

迄今为止,以往多见对p=2及p=3的情况用交替方向法求解过。仅作者在[1]中对p=2的情况提出两个显式格式,及文[2]中对p=3的情况提出一个显式格式。本文在此基础上,构造出两个对任意p维空间变量都适用的绝对稳定的显式格式,从而避免了解线代数方程组,大大地减小了计算工作量。     

一些非线性发展方程的显式行波解

该文给出了一种构造非线性发展方程显式行波解的方法并用该方法得到了Hirota-Satsuma方程组,一类非线性常微分方程以及广义耦合标量场方程组的显式行波解.     

带乘性噪声的空间分数阶随机非线性Schrödinger方程的广义多辛算法

带乘性噪声的空间分数阶随机非线性Schrödinger方程是一类重要的方程,可应用于描述开放非局部量子系统的演化过程.该方程为一个无穷维分数阶随机Hamilton系统,且具有广义多辛结构和质量守恒的性质.针对该方程的广义多辛形式,在空间上采用拟谱方法离散分数阶微分算子,在时间上则采用隐式中点格式,构造出一类保持全局质量的广义多辛格式.对行波解和平面波解等进行数值模拟,结果验证了所构造格式的有效性和保结构性质,时间均方收敛阶约在0.5到1之间.     

结构动力学方程的辛RK方法

针对有阻尼和外载荷的线性动力学常微分方程,给出了s级2s阶隐式Gauss-Legendre辛RK(Gauss-Legendre symplectic Runge-Kutta,GLSRK)方法的一种显式高效的执行格式,首次给出了Gauss-Legendre辛RK方法和经典RK方法(classical RK,CRK)的谱半径和单步相位误差的显式表达式,并将两者进行了比较.线性多自由度系统和非线性Rayleigh系统数值算例表明,对结构动力学系统而言,辛RK方法远比经典RK方法优越,在运动学特性和长时间数值模拟方面尤为明显.     

Sine-Gordon方程的多辛Leap-frog格式

非线性发展方程由于具有多种形式的解析解而吸引着众多的研究者,借助多辛保结构理论研究了Sine-Gordon方程的多辛算法．利用Hamilton变分原理,构造出了sine-Gordon方程的多辛格式;采用显辛离散方法得到了Leap-frog多辛离散格式,该格式满足多辛守恒律;数值结果表明leap-frog多辛离散格式能够精确地模拟sine-Gordon方程的孤子解和周期解,模拟结果证实了该离散格式具有良好的数值稳定性．     

浅水长波近似方程组的拟小波解

以浅水长波近似方程组为例,提出了拟小波方法求解(1 1)维非线性偏微分方程组数值解,该方程用拟小波离散格式离散空间导数,得到关于时间的常微分方程组,用四阶Runge-K utta方法离散时间导数,并将其拟小波解与解析解进行比较和验证.     

多维抛物型方程的分支绝对稳定的显式格式

其中及R={0≤x_i≤1,j=1,2,…,p),(?)R只为区域只的边界。 对多维抛物型方程(1)的差分解法,古典显式格式的稳定性条件为r=Δt/(Δx)~2≤1/2p,十分苛刻;古典隐式格式虽是无条件稳定,却需解线性方程组。因此两者的计算量都很大,且它们的精度较低,其局部截断误差仅为O(Δt+(Δx)~2)。因此,对多维抛物型方程而言,构造显式计算、稳定性能良好且精度较高的差分格式便具有十分明显的理论意义和实用价值。本文针对上述古典显式与隐式格式所存在的问题,构造一类对任何p维空间变量的抛物型方程(1)都适用的。分支绝对稳定的显式差分格式,其局部截断误差阶为O((Δt)~2+(Δx)~2),从而避免了解线性代数方程组,大大地减少了计算工作量,且精度较高。 令Δx_k=h_k=Δx=h=1/M(k=1,2,…p)表示空间方向步长,Δt=τ=[T/N]表示时间方向步长,M、N均为正整数。 为简便计,引入下列记号     

二维非线性Schr(o)dinger方程的辛与多辛格式

对满足周期边界条件的二维非线性Schrödinger方程,运用中心差分对该方程进行空间离散, 得到一个有限维Hamilton系统,然后用隐式Euler中点格式进行时间离散得到其辛格式. 针对该方程的多辛形式, 运用有限体积法离散,得到一种直平行六面体上的中点型多辛格式. 用所构造的辛与多辛格式对二维非线性Schrödinger方程的平面波解和奇异解进行数值模拟,结果验证了所构 造格式的有效性. 最后, 根据计算结果,对两种格式进行了分析和比较.       

传统隐式Runge-Kutta方法的转换方法

多级隐式Runge-Kutta(RK)方法簇中,除Gauss类方法是s级2s阶的辛方法以外,Radau类方法和Lobatto类方法既不是s级2s阶的方法也不是辛方法.基于隐式RK方法是一类转换RK方法这一特征,利用V-变换和Pade对角逼近,提出了构造高阶RK方法的转换定理.依据转换定理,导出了s级2s阶的Radau方法和s级2s阶的Lobatto方法.利用V-变换和待定系数法,导出了辛Radau方法和辛Lobatto方法.在此基础上,发现并证明了辛Radau方法是s级2s阶的方法.     

A CLASS OF EXPLICIT RATIONAL SYMPLECTIC INTEGRATORS

In this paper, a class of rational explicit symplectic integrators for one-dimensional oscillatory Hamiltonian problems are presented. These methods are zero-dissipative, and of first algebraic order and high phase-lag order. By means of composition technique, we construct second and fourth order methods with high phase-lag order of this type. Based on our ideas, three applicable explicit symplectic schemes with algebraic order one, two and four are derived, respectively. We report some numerical results to illustrate the good performance of our methods.     

EXISTENCE AND ITERATION OF POSITIVE SYMMETRIC SOLUTIONS TO A MULTI-POINT BOUNDARY VALUE PROBLEM

In this paper,we consider the existence of symmetric solutions to a nonlinear second order multi-point boundary value problem,and establish corresponding iterative schemes based on the monotone iterative method.     

SOME EXACT SOLUTIONS OF 3-DIMENSIONAL ZERO-PRESSURE GAS DYNAMICS SYSTEM

The 3-dimensional zero-pressure gas dynamics system appears in the modeling for the large scale structure formation in the universe. The aim of this paper is to construct spherically symmetric solutions to the system. The radial component of the velocity and density satisfy a simpler one dimensional problem. First we construct explicit solutions of this one dimensional case with initial and boundary conditions. Then we get special radial solutions with different behaviours at the origin.     

Numerical methods for coupled systems of quasilinear elliptic equations with nonlinear boundary conditions

This paper is concerned with numerical solutions of a coupled system of arbitrary number of quasilinear elliptic equations under combined Dirichlet and nonlinear boundary conditions. A finite difference system for a transformed system of the quasilinear equations is formulated, and three monotone iterative schemes for the computation of numerical solutions are given using the method of upper and lower solutions. It is shown that each of the three monotone iterations converges to a minimal solution or a maximal solution depending on whether the initial iteration is a lower solution or an upper solution. A comparison result among the three iterative schemes is given. Also shown is the convergence of the minimal and maximal discrete solutions to the corresponding minimal and maximal solutions of the continuous system as the mesh size tends to zero. These results are applied to a heat transfer problem with temperature dependent thermal conductivity and a Lotka-Volterra cooperation system with degenerate diffusion. This degenerate property leads to some interesting distinct property of the system when compared with the non-degenerate semilinear systems. Numerical results are given to the above problems, and in each problem an explicit continuous solution is constructed and is used to compare with the computed solution     

Symbolic computation of high order compact difference schemes for three dimensional linear elliptic partial differential equations with variable coefficients

We present a symbolic computation procedure for deriving various high order compact difference approximation schemes for certain three dimensional linear elliptic partial differential equations with variable coefficients. Based on the Maple software package, we approximate the leading terms in the truncation error of the Taylor series expansion of the governing equation and obtain a 19 point fourth order compact difference scheme for a general linear elliptic partial differential equation. A test problem is solved numerically to validate the derived fourth order compact difference scheme. This symbolic derivation method is simple and can be easily used to derive high order difference approximation schemes for other similar linear elliptic partial differential equations.     

Absolute stability of explicit Runge-Kutta-Nyström methods for y″=ƒ(x,y, y′)

We examine absolute stability of s-stage explicit Runge-Kutta-Nyström (R-K-N) methods of order s for $\text{s=2, 3, 4}\text{for}\text{y″=?(x, y, y′)}$ by applying these methods to the test equation: y″+2λy′+λ²y=0, λ>0. We show the existence of R-K-N methods of orders two, three and four possessing intervals of absolute stability as large as that of explicit Runge-Kutta (r-K) methods of respective orders.     

一类具非局部边界条件的四阶非线性微分方程的对称正解

本文考虑形如的非线性四阶微分方程非局部边值问题,这里a,b∈L~1[0,1],g:(0,1)→[0,∞)在(0,1)上连续、对称,且可能在t=0和t=1处奇异．f:[0,1]×[0,∞)→[0,∞)连续且对所有x∈[0,∞],f(·,x)在[0,1]上对称．在某些适当的增长性条件下,应用Krasnoselskii不动点定理证明了对称正解的存在性和多重性．     

Regularity Criteria of the Solutions to Axisymmetric Magnetohydrodynamic System

In this paper, we consider Cauchy problem of the axially symmetric Magnetohydrodynamic (MHD) system. By using energy method, we establish some regularity criteria of the solutions for the axisymmetric solutions of the three dimensional incompressible MHD system.     

GLOBAL EXISTENCE OF SOLUTIONS FOR A MULTI-PHASE FLOW: A BUBBLE IN A LIQUID TUBE AND RELATED CASES

In this paper we study the problem of the global existence(in time) of weak,entropic solutions to a system of three hyperbolic conservation laws, in one space dimension,for large initial data. The system models the dynamics of phase transitions in an isothermal fluid; in Lagrangian coordinates, the phase interfaces are represented as stationary contact discontinuities. We focus on the persistence of solutions consisting in three bulk phases separated by two interfaces. Under some stability conditions on the phase configuration and by a suitable front tracking algorithm we show that, if the BV-norm of the initial data is less than an explicit(large) threshold, then the Cauchy problem has global solutions.     

Power law Stokes equations with threshold slip boundary conditions: Numerical methods and implementation

For the power law Stokes equations driven by nonlinear slip boundary conditions of friction type, we propose three iterative schemes based on augmented Lagrangian approach and interior point method to solve the finite element approximation associated to the continuous problem. We formulate the variational problem which in this case is a variational inequality and construct the weak solution of the continuous problem. Next, we formulate two alternating direction methods based on augmented Lagrangian formalism in order to separate the velocity from the symmetric part the velocity gradient and tangential part of the velocity. Thirdly, we present some salient points of a path‐following variant of the interior point method associated to the finite element approximation of the problem. Some numerical experiments are performed to confirm the validity of the schemes and allow us to compare them.