四阶杆振动方程的两类高精度辛格式

In this paper, we present two classes of symplectic schemes with high order accuracy for solving four-order rod vibration equation utt uxxxx=0 via the third type generating function method. First, the equation of four order rod vibration is written into the canonical Hamilton system; second, overcoming successfully the essential difficult on the calculus of high order variations derivative, we get the semi-discretization with arbitrary order of accuracy in time direction for the PDEs by the third type generating function method. Furthermore the discretization of the related modified equation of original equation is obtained. Finally, arbitrary order accuracy symplectic schemes are obtained. Numerical results are also presented to show the effectiveness of the scheme, high order accuracy and properties of excellent long-time numerical behavior.     

Navier-Stokes方程的低阶混合有限元方法

1引言 定常N-s方程是流体力学中一类非常重要的方程,而经典的混合有限元方法要求有限元空间组合满足B-B条件.这一条件限制了工程中常用的低阶有限元空间如:P1/P1,P,/Po等.为了去掉LBB条件限制,产生了一种新的方法--稳定化方法(也成CBB方法).1988年,F.Brezzi和J.Douglas.Jr对线性的Stokes方程建立了一种稳定化格式([2]).对于低阶的有限元应用压力投影稳定项构造了一种稳定化格式,并给出了格式的解存在唯一性,且给出了几种有限元的算例.     

一类反应扩散方程组的保持正性的差分格式

本文研究一类具有正解的反应扩散方程组的有限差分解法.构造了一个保持正性的差分格式.利用离散的最大值原理证明了差分格式解的非负性,有界性及差分格式的无条件稳定性.这些估计的证明不依赖于微分方程的解而仅仅与初边值条件有关.当微分方程的解适当光滑时,证明了差分格式的一致收敛性.最后给出了数值计算结果,并与以往方法进行了比较.计算结果说明了本文给出的方法的有效性.     

一类偏积分微分方程二阶差分空间半离散格式的全局行为

偏积分微分方程产生于许多科学与工程领域,数值求解此类问题具有重要应用.本文给出了数值求解一类长时间偏积分微分方程的二阶差分空间半离散格式.借助于Laplace变换及Parseval等式,给出了全局稳定性的证明、误差估计及全局收敛性的结果.     

对流扩散方程的数值计算

本文研究了对流扩散方程的一种并行格式.利用一组saul'yev型非对称格式进行二次构造,分别得到了一类并行GE格式和GEL、GER格式;进一步推广,得到绝对稳定的交替分组显式AGE格式,并用数值例子检验AGE格式的数值计算效果.     

一类四阶抛物积分微分方程混合元方法的超收敛分析

本文的主要目的是利用双线性元Q_(11)及Q_(01)×Q_(10)元研究一类非线性四阶抛物积分微分方程的混合有限元方法.一方面,利用上述两种元的高精度结果以及对时间t的导数转移技巧,在半离散格式下,导出原始变量u和中间变量w=-?u在H~1-模意义下及流量p(向量)=-?u在(L~2)~2-模意义下具有O(h~2)阶的超逼近性质.进一步地,借助插值后处理技术,得到上述变量的整体超收敛结果.另一方面,建立一个新的向后Euler全离散格式.通过采取新的分裂技术,得到u和w在H~1-模意义下及p在(L~2)~2-模意义下具有O(h~2+?t)阶的超逼近和超收敛结果.这里,h和?t分别表示空间剖分参数和时间步长.最后,给出一个数值算例,计算结果验证了理论分析的正确性.     

五次NLS方程全离散谱格式的大时间性态

Nonlinear Schroedinger equation arises in many physical problems. There are many works in which properties of the solution are studied. In this paper we use fully discrete Fourier spectral method to get an approximation solution of nonlinear weakly dissipative Schroedinger equation with quintic term. We give a large-time error estimate and obtain the existence of the approximate attractor A N^k.     

一种多介质流体计算的三维非守恒型差分格式

Motivated by the need for three-dimensional methods for interface calculations we describe a 3D non-conservative difference scheme in Lagrange coordinates for calculating multi-component fluid dynamics, according to the thoughts of 2D non-conservative difference scheme and making use of difference method of characteristic and Upwind form. Furthermore, we give a three-dimensional numerical simulation of explosion in concrete, and the results of simulation are consistent with examination.     

Accuracy analysis of immersed boundary method using method of manufactured solutions

The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics.This paper analyzes the accuracy of the immersed boundary method.The procedure contains two parts,i.e.,the code verification and the accuracy analysis.The code verification provides the confidence that the code used is free of mistakes,and the accuracy analysis gives the order of accuracy of the immersed boundary method.The method of manufactured solutions is taken as a means for both parts.In the first part,the numerical code employs a second-order discretization scheme,i.e.,it has second-order accuracy in theory.It matches the calculated order of accuracy obtained in the numerical calculation for all variables.This means that the code contains no mistake,which is a premise of the subsequent work.The second part introduces a jump in the manufactured solution for the pressure and adds the corresponding singular forcing terms in the momentum equations.By analyzing the discretization errors,the accuracy of the immersed boundary method is proven to be first order even though the discretization scheme is second order.It has been found that the coarser mesh may not be sensitive enough to capture the influence of the immersed boundary,and the refinement on the Lagrangian markers barely has any effect on the numerical calculation.     

Nonergodic Brownian Motion

The long-time behavior of a system is suggested to confirm nonergodicity of non-Markovian Brownian dynamics, namely, whether the stationary probability density function （PDF） of the system characterized mainly by low moments of variables depends on the initial preparation. Thus we classify nonergodic Brownian motion into two classes： the class-I is that the PDF of a force-free particle depends on the initial velocity and the equilibration can be recovered through a bounded potential; while the PDF in the class-H depends on the initial coordinate and the equilibration can not be approached by introducing any potential. We also compare our result with the conditions of three kinds for ergodicity.