求解耦合Schrdinger方程组的数值逼近算法(英文)

对耦合Schrdinger方程组提出一个非耦合的线性化差分格式并对其进行分析.证明格式保持原方程组的守恒律,在先验估计的基础上证明格式依L2模的绝对稳定性和无条件二阶收敛性.对孤波碰撞的各种现象进行模拟.     

SRLW方程的多辛Fourier谱格式及其守恒律

通过引进正则动量,将对称正则长波方程（简称SRLW方程）转化成多辛形式的方程组,它具有多辛守恒律;介绍了空间方向满足周期边界条件的函数的Fourier谱方法;对SRLW方程的多辛方程组在空间方向利用Fourer谱方法,时间方向上应用Euler中点格式离散,得到其多辛Fourier拟谱格式;证明此格式的一些离散守恒律．用此格式模拟了SRLW方程的单个孤立波,还模拟了多个孤立波的追赶、碰撞和分离过程．     

加权ENO格式的构造及数值模拟

将加权ENO格式推广到非结构三角形网格上,构造了一类加权ENO有限体积格式,提出的插值多项式的构造方式,可以减少计算时间.对于出现的病态方程组,给出了解决方法.此外还给出了插值点的选取方式及加权因子的构造方法.结合三阶TVD Runge Kutta时间离散,对二维欧拉方程组进行了数值试验.     

光栅作用下双波导间的反向耦合

本文根据非正交模理论导出了描述光栅作用下双波导间反向耦合的耦合模方程组.通过数值计算对这种结构的反向耦合特性进行了初探,并同光栅作用下单波导反向耦合特性作了比较.     

一类广义的高阶非线性Schr.dinger方程组的差分格式

考虑一类广义的高阶非线性Schrodinger方程组,给出蛙跳格式及混合Crank-Nicolson蛙跳格式,并证明其稳定性与收敛性。最后,数值例子表明数值结果与理论分析相一致。     

气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

非线性耦合微分方程组的精确解析解

提出了利用耦合的Riccati方程组的某些特解构造非线性微分方程组精确解析解的一种方法.应用这种方法研究了两个耦合的常微分方程组,系统地获得了它们的一些精确解.给出了非线性浅水波近似方程组和非线性Schr?dinger-KdV方程组若干新的孤波解. 关键词： 非线性耦合方程组 Riccati方程组 符号计算 孤波     

双钨极耦合电弧数值模拟

基于流体力学方程组和麦克斯韦方程组, 在合理的边界条件下, 建立了双钨极耦合电弧三维准静态数学模型. 通过对方程组的迭代求解, 获得了不同钨极间距和电弧长度下耦合电弧的温度场、流场、电弧压力和电流密度分布等重要结果, 与已有的实验研究符合良好. 模拟结果表明: 与相同条件下的钨极惰性气体保护焊电弧相比, 双钨极耦合电弧的最高温度和最大等离子流速较低, 阳极表面电弧压力和电流密度峰值明显减小; 钨极间距和弧长对耦合电弧的温度场、流场、电流密度和电弧压力等都具有显著的影响, 且耦合电弧阳极的电弧压力和电流密度分布不能用高斯近似进行描述. 关键词： 耦合电弧 三维模型 数值模拟     

二维交错网格的GAUSS型格式

利用Gauss型求积公式在交错网格的情况下构造了一类不需解Riemann问题的求解二维双曲守恒律的二阶显式Gauss型差分格式,该格式在CFL条件限制下为MmB格式.并将格式推广到二维方程组,进行了数值试验.     

TVD数值技术在弱电离化学非平衡斜激波反射流场计算中的应用

本文成功地将二阶迎风ＴＶＤ格式运用到弱电离，化学非平衡斜激波反射流场的数值计算中，气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的效方向，文中详细地给出了气体力学方程和组元质量守恒方程组全耦合求解过程，最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论，本文所得结果与文献［４］由Ｉ．Ｉ．Ｇｌａｓｓ提供的实验结果和文献［１］由Ｈ．Ｍ．Ｇｌａｚ提供的数值结果相符合。     

A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis

A high-order accurate hybrid central-WENO scheme is proposed. The fifth order WENO scheme [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] is divided into two parts, a central flux part and a numerical dissipation part, and is coupled with a central flux scheme. Two sub-schemes, the WENO scheme and the central flux scheme, are hybridized by means of a weighting function that indicates the local smoothness of the flowfields. The derived hybrid central-WENO scheme is written as a combination of the central flux scheme and the numerical dissipation of the fifth order WENO scheme, which is controlled adaptively by a weighting function. The structure of the proposed hybrid central-WENO scheme is similar to that of the YSD-type filter scheme [H.C. Yee, N.D. Sandham, M.J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, J. Comput. Phys. 150 (1999) 199–238]. Therefore, the proposed hybrid scheme has also certain merits that the YSD-type filter scheme has. The accuracy and efficiency of the developed hybrid central-WENO scheme are investigated through numerical experiments on inviscid and viscous problems. Numerical results show that the proposed hybrid central-WENO scheme can resolve flow features extremely well.     

基于Voronoi网格的扩散方程差分格式

在Voronoi网格上利用一种基于回路积分法的有限体积法构造扩散方程的的差分格式.在这种特殊的网格上离散扩散方程比通常在四边形网格上离散的格式要简单,不会引进角点未知量,提高了对网格边上的流的离散精度,及差分格式整体精度.这种Voronoi网格上的扩散计算也可以与单元中心流体力学计算耦合.数值算例表明这种格式比四边形网格上的格式精度高,且能更好的应对网格扭曲情形.     

Two expanding forms of a Lie algebra and their application

With the help of a known Lie algebra,two new high order Lie algebras are constructed.It is remarkable that they have different constructing approaches.The first Lie algebra is constructed by the definition of integrable couplings.the second one by an extension of Lie algebra,Then by making use of Tu scheme,a generalized AKNS hierarchy and another new hierarchy are obtained.As a reduction case of the first hierarchy,a kind of coupled KdV equation is presented.As a reduction case of the second one,a new coupled Schroedinger equation is given.     

Conservation Quantities of the Explicit Symplectic Scheme for Time-evolution of Quantum System

ＣｏｎｓｅｒｖａｔｉｏｎＱｕａｎｔｉｔｉｅｓｏｆｔｈｅＥｘｐｌｉｃｉｔＳｙｍｐｌｅｃｔｉｃＳｃｈｅｍｅｆｏｒＴｉｍｅ－ｅｖｏｌｕｔｉｏｎｏｆＱｕａｎｔｕｍＳｙｓｔｅｍＺｈｏｕＺｈｏｎｇｙｕａｎＤｉｎｇＰｅｉｚｈｕＩｎｓｔｉｔｕｔｅｏｆＡｔｏｍｉｃａｎ...     

积分梯度法的两种格式

在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点.     

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(AUSM)思想提出一种通量分裂方法,称为K-CUSP格式.它与传统H-CUSP和E-CUSP格式的最大差异在于总能量的分裂:K-CUSP格式将无粘守恒通量中所有的运动学量分裂到对流项,所有的热力学量分裂到压力项,即总能量被分裂成动能和静焓.对于压力项的数值通量,采用一种新的界面构造方法.数值测试表明:①K-CUSP格式继承了FVS格式的简单性和稳健性.在激波后不易出现压力过冲,在膨胀区域没有振荡,优于AUSM和WPS格式;②K-CUSP格式继承了FDS格式的分辨率.激波间断的分辨率和H-CUSP、E-CUSP格式基本相同,接触间断的分辨率高于FVS格式,低于Roe、AUSM和WPS格式.AUSM和WPS格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.     

Feedback control scheme of traffic jams based on the coupled map car-following model

Based on the pioneering work of Konishi et al. [Phys. Rev. E （1999） 60 4000], a new feedback control scheme is presented to suppress traffic jams based on the coupled map car-following model under the open boundary condition. The effect of the safe headway on the traffic system is considered. According to the control theory, the condition under which traffic jams can be suppressed is analyzed. The results are compared with the previous results concerning congestion control. The simulations show that the suppression performance of our scheme on traffic jams is better than those of the previous schemes, although all the schemes can suppress traffic jams. The simulation results are consistent with theoretical analyses.     

Synchronization of impulsively coupled complex networks

We investigate the synchronization of complex networks,which are impulsively coupled only at discrete instants.Based on the comparison theory of impulsive differential systems,a distributed impulsive control scheme is proposed for complex dynamical networks to achieve synchronization.The proposed scheme not only takes into account the influence of all nodes to network synchronization,which depends on the weight of each node in the network,but also provides us with a flexible method to select the synchronized state of the network.In addition,it is unnecessary for the impulsive coupling matrix to be symmetrical.Finally,the proposed control scheme is applied to a chaotic Lorenz network and Chua’s circuit network.Numerical simulations are used to illustrate the validity of this control scheme.     

A Fifth-Order Low-Dissipative Conservative Upwind Compact Scheme Using Centered Stencil

In this paper, a conservative fifth-order upwind compact scheme using centered stencil is introduced. This scheme uses asymmetric coefficients to achieve the upwind property since the stencil is symmetric. Theoretical analysis shows that the proposed scheme is low-dissipative and has a relatively large stability range. To maintain the convergence rate of the whole spatial discretization, a proper non-periodic boundary scheme is also proposed. A detailed analysis shows that the spatial discretization implemented with the boundary scheme proposed by Pirozzoli [J. Comput. Phys., 178 (2001), pp. 81-117] is approximately fourth-order. Furthermore, a hybrid methodology, coupling the compact scheme with WENO scheme, is adopted for problems with discontinuities. Numerical results demonstrate the effectiveness of the proposed scheme.