粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

高阶Schrdinger方程的高精度辛格式

提出由第三类生成函数法构造高阶Schr dinger方程ut=i(-1)m2mux2m的高精度辛格式.首先,给出它的典则Hamilton方程组;然后,成功地克服了本质上是困难的高阶变分导数的计算,并利用第三类生成函数法得到在时间方向具有任意阶精度的半离散方程,进而得到原始方程相关的修正方程的离散形式,最后得到各种精度的辛格式.数值结果表明该格式是有效的,具有高精度及良好的长时间数值行为等特性.     

求解对流扩散方程的紧致修正方法

提出了求解对流扩散方程的紧致修正方法,该方法是在低阶离散格式的源项中,引入紧致修正项,从而构造高阶紧致修正格式,并进行求解.采用紧致修正方法对典型的对流扩散方程进行计算.结果表明,紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上,但紧致修正方法的精度与紧致方法的精度相同,均具有四阶精度.所以紧致修正方法可以在少网...     

高阶Schr(o)dinger方程的高精度辛格式

提出由第三类生成函数法构造高阶Schroedinger方程δu/δt=i(-1)^nδ^2mu/δx^2m的高精度辛格式．首先，给出它的典则Hamilton方程组；然后，成功地克服了本质上是困难的高阶变分导数的计算，并利用第三类生成函数法得到在时间方向具有任意阶精度的半离散方程，进而得到原始方程相关的修正方程的离散形式，最后得到各种精度的辛格式．数值结果表明该格式是有效的，具有高精度及良好的长时间数值行为等特性．     

扩散方程的有限方向差分方法

研究二维散乱点集上数值求解非线性扩散方程的有限方向差分方法。利用五个邻点信息构造具有最小模板的离散格式,并且离散系数具有显式表达式。另外,利用五点公式获得了间断问题物质界面的离散格式,该格式对界面流的计算具有近似二阶精度。不同计算区域及不同类型的离散点集上的计算结果验证了方法的有效性。     

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

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

守恒双曲方程一类本质非振荡(ENO)格式

用自适应Newton插值,结合自适应模型和重构思想去构造数值流通量,对时间采用Runge-Kutta型离散,得到一类不需"真正"插值和数值微分过程的ENO格式。该格式易于数值实现,数值试验表明,这类格式具有良好的计算结果。     

隐式格式求解拟压缩性非定常不可压Navier-Stokes方程

采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟.     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

Symplectic Euler Method for Nonlinear High Order Schrödinger Equation with a Trapped Term

In this paper, we establish a family of symplectic integrators for a class of high order Schrödinger equations with trapped terms. First, we find its symplectic structure and reduce it to a finite dimensional Hamilton system via spatial discretization. Then we apply the symplectic Euler method to the Hamiltonian system. It is demonstrated that the scheme not only preserves symplectic geometry structure of the original system, but also does not require to resolve coupled nonlinear algebraic equations which is different from the general implicit symplectic schemes. The linear stability of the symplectic Euler scheme and the errors of the numerical solutions are investigated. It shows that the semi-explicit scheme is conditionally stable, first order accurate in time and $2l^{th}$ order accuracy in space. Numerical tests suggest that the symplectic integrators are more effective than non-symplectic ones, such as backward Euler integrators.     

三能级冷原子介质中多个光孤子的相互作用

冷原子介质中的光孤子在电磁感应透明(EIT)的作用下表现出很多奇异的特性,对描述这些特性的理论模型的研究在光信号处理和传输方面具有重要的意义. 描述三能级冷原子EIT介质中空间孤立子演化的二维饱和非线性薛定谔方程被转化成辛结构的Hamilton系统, 利用辛几何算法离散Hamilton系统得到了相应离散的辛格式,并且利用辛格式数值模拟了三能级冷原子EIT介质中在相同振辐不同相位的两个、四个光孤子的相互作用行为. 数值实验结果表明: 冷原子介质中多个光孤子的相互作用行为不但与入射高斯光束的相位有关,还和入射高斯光束的方向有关. 入射的高斯光束能在冷原子介质中形成稳定的孤立子.     

Multi-Symplectic Wavelet Collocation Method for Maxwell's Equations

In this paper, we develop a multi-symplectic wavelet collocation method for three-dimensional (3-D) Maxwell's equations. For the multi-symplectic formulation of the equations, wavelet collocation method based on autocorrelation functions is applied for spatial discretization and appropriate symplectic scheme is employed for time integration. Theoretical analysis shows that the proposed method is multi-symplectic, unconditionally stable and energy-preserving under periodic boundary conditions. The numerical dispersion relation is investigated. Combined with splitting scheme, an explicit splitting symplectic wavelet collocation method is also constructed. Numerical experiments illustrate that the proposed methods are efficient, have high spatial accuracy and can preserve energy conservation laws exactly.     

Symplectic wavelet collocation method for Hamiltonian wave equations

This paper introduces a novel symplectic wavelet collocation method for solving nonlinear Hamiltonian wave equations. Based on the autocorrelation functions of Daubechies compactly supported scaling functions, collocation method is conducted for the spatial discretization, which leads to a finite-dimensional Hamiltonian system. Then, appropriate symplectic scheme is employed for the integration of the Hamiltonian system. Under the hypothesis of periodicity, the properties of the resulted space differentiation matrix are analyzed in detail. Conservation of energy and momentum is also investigated. Various numerical experiments show the effectiveness of the proposed method.     

A Multi-Symplectic Compact Method for the Two-Component Camassa-Holm Equation with Singular Solutions

The two-component Camassa-Holm equation includes many intriguing phenomena. We propose a multi-symplectic compact method to solve the two-component Camassa-Holm equation. Based on its multi-symplectic formulation, the proposed method is derived by the sixth-order compact finite difference method in spatial discretization and the symplectic implicit midpoint scheme in temporal discretization. Numerical experiments finely describe the velocity and density variables in the two-component integrable system and distinctly display the evolvement of the singular solutions. Moreover, the proposed method shows good conservative properties during long-time numerical simulation.     

含时Schrödinger方程的高阶辛FDTD算法研究

提出了一种新的算法——高阶辛时域有限差分法(SFDTD(3, 4): symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散, 空间上采用四阶精度的同位差分格式离散, 建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时, 辛积分可提高算法的稳定度;SFDTD(3, 4)法和FDTD(2, 4)法较传统的FDTD(2, 2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明: SFDTD(3, 4)法较传统的FDTD(2, 2)法及高阶FDTD(2, 4)法有着更好的计算精度和收敛性, 且SFDTD(3, 4)法能够保持量子系统的能量守恒, 适用于长时间仿真.     

Analysis of numerical dispersion properties of SFDTD and MRTD schemes

In this paper, the Maxwell's equations are written as Hamilton canonical equations by using Hamilton functional variation method. Maxwell's equations can be discretized with symplectic propagation technique combined with high-order difference schemes approximations to construct symplectic finite difference time domain (SFDTD) method. The high-order dispersion equations of the scheme for space is deduced. The numerical dispersion analysis is included, and it is compared with the multiresolution time-domain (MRTD) method based on the Daubechies scaling functions. Numerical results show high efficiency and accuracy of the SFDTD method.     

前缘曲率变化对平板边界层感受性问题的影响

边界层感受性问题是层流向湍流转捩的初始阶段,是实现边界层转捩预测和控制的关键环节.目前已有的研究成果显示,在声波扰动或涡波扰动作用下前缘曲率变化对边界层感受性机制有着显著的影响.本文采用直接数值模拟方法,研究了在自由来流湍流作用下具有不同椭圆形前缘平板边界层感受性问题,揭示椭圆形前缘曲率变化对平板边界层内被激发出Tollmien-Schlichting (T-S)波波包的感受性机制以及波包向前传播群速度的影响;通过快速傅里叶分析方法从波包中提取获得了不同频率的T-S波,详细分析了前缘曲率变化对不同频率的T-S波的幅值、色散关系、增长率、相速度以及形状函数的作用;确定了前缘曲率在平板边界层内激发T-S波的感受性过程中所占据的地位.通过上述研究能够进一步认识和理解边界层感受性机制,从而丰富和完善了流动稳定性理论.     

On the quasioptical approximation in dissipative media with spatial dispersion

A new form of the quasioptical equation is proposed to describe the propagation of an electromagnetic wave beam in a stationary smoothly inhomogeneous medium with spatial dispersion and dissipation. The proposed approach guarantees the positive definiteness of the absorbed power in the locally dissipative medium, which is a nontrivial property for inhomogeneous media with spatial dispersion. An efficient numerical scheme is constructed to solve the derived quasioptical equation.     

Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge–Kutta–Nyström (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge–Kutta method, which is equivalent to the well-known Störmer–Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine–Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation.