非线性Schr(o)dinger方程的辛算法

给出了非线性Schrodinger方程的二阶Euler中点格式、四阶Euler中点格式、二阶蛙跳格式和四阶蛙跳格式，并且作了数值实验验证这些格式的可行性并比较其误差．并且对同样截断误差阶的一种辛格式和一种非辛的差分格式进行比较．我们选取二阶蛙跳格式和二阶两层格式作了数值实验并对它们的运行结果作了比较．发现辛算法比同样截断误差的非辛算法误差小，时间越长优势越明显．     

非线性Schrdinger方程的辛算法

给出了非线性Schrdinger方程的二阶Euler中点格式、四阶Euler中点格式、二阶蛙跳格式和四阶蛙跳格式,并且作了数值实验验证这些格式的可行性并比较其误差.并且对同样截断误差阶的一种辛格式和一种非辛的差分格式进行比较.我们选取二阶蛙跳格式和二阶两层格式作了数值实验并对它们的运行结果作了比较.发现辛算法比同样截断误差的非辛算法误差小,时间越长优势越明显.     

一类非线性Sine－Gordon方程的一个“蛙跳”有限差分格式

本文给出了一类非线性Ｓｉｎｅ－Ｇｏｒｄｏｎ方程的一个“蛙跳”有限差分格式．讨论了三对角循环线性系统的微分方程的可能性，进一步的，利用能量不等式的方法研究了差分格式的收敛性和稳定性，最后，给了一个数值解的例子     

数值计算用于瞬态四波混频理论分析

采用差分方法数值求解瞬态四波混频耦合波方程，给出两种差分格式并进行数值计算，所得结果与实验基本相符。     

弹性波方程的紧致差分方法

在对弹性波方程进行数值模拟时 ,低阶差分格式往往产生严重的数值频散 ,高阶显示差分格式需要用较多的网格点 ,不利于边界的处理。而紧致差分格式吸收了它们的优点 ,弥补了它们的不足。为此该文应用紧致差分格式的思想 ,发展了二维情况下弹性波方程初值问题的紧致差分方法 ,研究了它的稳定性 ,并用 Fourier方法分析了显示差分格式和紧致差分格式的相速度误差 ,最后利用紧致差分方法在粗网格条件下对地震波传播进行了数值模拟 ,并同五点四阶中心差分方法的计算结果进行了对比。结果表明 ,求解弹性波方程的紧致差分方法有效 ,且具有比同网格点差分格式更高的计算精度和较小的数值频散。     

Equal-Width波方程的高精度守恒差分格式

对Equal-Width波方程提出一个三层线性高精度守恒差分格式.所建格式满足质量守恒和能量守恒,在时间和空间上分别为二阶和四阶精度.用离散能量法证明了所建差分格式的收敛性和稳定性.数值实验表明该格式是有效的和可靠的.     

Rosenau-Burgers方程的一个新的差分方法

从动力学系统的实际问题出发,针对Rosenau-Burgers方程的初边值问题进行了数值研究,揭示了复杂离散动态系统理论中非线性波耗散问题. 在方程求解的时间和空间区域,采用网格化方法,提出了一个新的三层隐式差分格式,对差分解进行了先验估计,并给出了该格式的稳定性和收敛性的严格理论证明. 数值实验的结果表明,差分格式简单而有效、计算速度快、稳定性好,并且差分格式使用了加权方法,使其具有普遍意义和推广价值.     

对流方程的六阶中心差分格式

有限差分方法是微分方程数值解法中发展最早、理论最完善、应用最广泛的计算方法之一．利用待定系数法构造了对流方程的中心有限差分格式,利用Taylor级数展开推导出了该差分格式的修正偏微分方程（MPDE）,采用数值余项效应分析方法从空间离散方面改进了该格式．利用高阶TVD Runge-Kutta方法从时间离散方面改进了该格式．利用Richardson外推方法在不增加计算复杂度的前提下改革了原格式．数值实验表明本文讨论的3种方法在差分格式改进和优化中的有效性．本文讨论的方法也可以用于其他偏微分方程有限差分方法的构造中．     

求解波动方程的高精度紧致隐式差分方法

基于二阶微商的二阶中心差商和四阶紧致差商逼近公式及其加权平均思想,推导出了数值求解一维波动方程的2种精度分别为O（x^2＋h^4）和O（x^4＋h^4）的三层隐式紧致差分格式,以夏与之相匹配的第一个时间步的同阶离散格式,并采用Fourier方法分析了格式的稳定性．由于每一时间层上最多只用到了3个网格点,所以可采用追赶法直接求解差分方程．数值实验结果验证了所得方法的精确性和可靠性．     

高阶紧致分裂方法在确定性和随机薛定谔方程中的应用

采用高阶紧致分裂方法解决确定性和随机非线性薛定谔方程,其中在时间方向上用Strange-type分裂方法进行离散,而在空间方向上用四阶紧致差分格式离散,并给出相应的电荷守恒量和能量守恒量表达式。最后给出数值实验,观察孤立波在确定情形和随机情形下的表现形式,从数值结果上反映该方法的优越性。     

RK-MRTD算法的稳定性及数值色散性分析

将龙格-库塔(Runge-Kutta)方法引入到时域多分辨分析(MRTD)算法,即在时间上采用Runge-Kutta方法离散,并用此算法解决传统时域有限差分(FDTD)算法较大的色散误差问题.对Runge-Kutta时域多分辨分析算法(RK-MRTD)的稳定性和数值色散性进行系统分析,数值结果表明该方法具有高精度性.     

随机微分方程的Runge-Kutta数值解法

提出了求解一类随机常微分方程(SODEs)的3种Runge-Kutta格式:显式Runge-Kutta格式、半隐式Runge-Kutta格式和隐式Runge-Kutta格式.讨论了这3种Runge-Kutta格式的T稳定条件,并给出了部分数值实验结果.     

三维水动力学模型高精度差分格式和解法研究

采用高精度差分格式对三维水动力学问题的时均化Navier-Stokes方程进行了数值模拟，进而采用广义共轭剩余法（GCR方法）求解压力泊松方程，并采用显式三级二阶Runge-Kutts格式模拟了时间步进过程。傅立叶分析表明，文中所采用的三阶迎风紧致差分格式具有较高的精度，数值实验进一步验证了上述数值模型的准确性和有效性。     

五阶高分辨率熵稳定算法

针对双曲守恒律方程的数值求解问题,构造一种新型的熵稳定算法.新算法空间方向采用五阶中心加权基本无振荡(CWENO)重构格式,时间方向采用四阶强稳定龙格-库塔(Runge-Kutta)方法.将新算法应用于若干一维Burgers方程和Euler方程组问题数值算例的求解.结果表明:新算法精度高,有效抑制了伪振荡的产生,与理论分析的结果一致.     

叶片数约化对叶轮机非定常气动特性影响的研究

通过求解三维非定常N-S方程模拟叶轮机械内部流场。采用有限体积中心差分和四步龙格库塔时间推进格式,并且使用了双时间方法对计算加速。对一级跨音速风扇静叶选取不同的叶片数方案进行非定常流场和性能数值模拟。对于不同的叶片数,叶片几何截面按比例缩放以保证各方案叶片稠度不变。计算结果给出了风扇级的性能、叶片气动载荷和非定常气动力等参数。叶片数目选择影响到非定常效率波动的幅值大小。通过对非定常气动力力的频谱分析得出,不同的叶片数目对非定常气动力的前两阶幅值造成明显的影响。     

多辛Preissman格式及其应用

主要讨论了用于求解多辛哈密尔顿系统的多辛Preissman格式及其简单应用.根据多辛格式必须满足离散的多辛守恒律的基本思想,从Runge-Kutta方法入手,推导出其为多辛格式的充分条件,进而得到了多辛的中点格式,同时举例说明的它在偏微分方程数值求解中的应用.     

线性Boussinesq方程的多辛Runge-Kutta Nystr(o)m算法

考虑线性Boussinesq方程的多辛Hamilton形式, 利用Runge Kutta Nystrom算法离散此多辛结构, 得到了离散多辛守恒律, 并求得一个等价于Runge Kutta Nystrom积分的新格式, 证明了它的稳定性条件. 数值实验结果表明了理论分析的正确性.     

新的多重网格算法在压气机流场计算中的应用

根据在流场数值计算中使用多重网格方法可以加速收敛的原理,在叶栅流场控制方程使用有限体积差分法进行离散的前提下,提出了一种新的4重网格方法来加速求解压气机叶栅可压缩流动,该方法可以与显式时间推进法有效配合使用。通过对NASA67超音速轴流压气机叶栅流场的数值计算表明,该方法明显加快了流场计算的收敛速度,并且计算结果与试验结果吻合较好。     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

求解不可压流体Navier-Stokes方程的四阶精度有限容积紧致格式

采用四阶精度的有限容积紧致格式在交错网格上对二维非定常不可压流体的Navier-Stokes方程中的对流项和扩散项进行离散．压力项则由压力Poisson方程求得，并给出了新的压力Poisson方程的四阶精度有限容积紧致格式的离散表达式．用低存贮的三阶Runge-Kutta方法对Navier-Stokes方程进行时间推进．Fourier分析表明，有限容积紧致格式比一般的有限容积非紧致格式有更高的分辨率．最后以Taylor涡为例，得到了很好的结果．