用摄动配置方法求解含时薛定谔方程

将摄动配置方法应用到含时薛定谔方程,在计算实现的基础上结合摄动配置的特征提出了一类新的数值积分方法,并给出了一个2级2阶和一个3级4阶的辛摄动配置方法对含时薛定谔方程的数值算例.为了检验新的数值积分方法,我们还给出了与两个辛摄动配置格式在理论上等价的辛龙格-库塔方法以及同阶的非辛方法的数值模拟.展示了一些数值结果,并给出了一些分析.     

数据采集控制器ADAM—5510和上位机通信的实现

介绍了在VISUAL BASIC环境中研华公司的ADAM-5510数据采集控制器和上位机通信实现的方法。     

单晶Si Nd∶YAG激光连续退火的数值计算

应用固相外延模型来模拟单晶Si的连续Nd∶YAG激光退火过程,在低功率密度连续激光退火下,用准静态模型模拟辐照区向非辐照区的径向传导散热.在数值计算中,应用部分线性法处理非线性非齐次热传导方程,得到相应的隐格式差分方程,再用追赶法求解隐格式差分方程,得出绝热边界条件下的温度的时间和空间分布,从而得出激光退火的再结晶厚度.当激光波长λ=1.06μm、功率密度i0=700W/cm2,预热温度T0=523K时,经过0.7秒,表面温度度升到1 290K左右,再结晶厚度约为0.5μm.     

高精度非线性格式WCNS的分析研究与其应用

首先将Fourier方法推广于高维方向研究了五阶精度WCNS的特性,并与其他高阶格式进行比较。分析结果表明WCNS的高精度特性普遍接近甚至好于迎风偏置五阶显式格式EUW-5与Pade′标准格式。然后开展了WCNS的应用研究,采用高效率的WCNS-E-5数值模拟了含强激波的高维复杂流场。算例包括二维高超声速边界层对自由流扰动波的吸收问题以及三维球头绕流问题。计算结果反映出WCNS-E-5对激波等间断的良好捕捉能力,图像清晰光滑,数据准确可靠。     

随机结构非线性动力响应的概率密度演化分析

提出了随机结构非线性动力响应分析的概率密度演化方法．根据结构动力响应的随机状态方程，利用概率守恒原理，建立了随机结构非线性动力响应的概率密度演化方程．结合Newmark-Beta时程积分方法与Lax-Wendroff差分格式，提出了概率密度演化方程的数值分析方法．通过与Monte Carlo分析方法对比，表明所给出的概率密度演化方法具有良好的计算精度和较小的计算工作量．研究表明：随机结构非线性动力响应概率密度具有典型的演化特征，随着时间增长，概率密度曲线分布趋于复杂．     

平面叶栅跨音流场气动特性的数值研究

通过对模型方程的分析，给出了一种新的隐格式构造思想。将它运用到关通量分裂格式中，可得到无近似因子分解、无矩阵运算的高效二阶精度隐式矢通量分裂差分格式，并用来直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组。数值计算标明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

类升力体外形高超声速粘性绕流拓扑结构分析

采用交替方向隐式分解的隐式NND格式求解全N-S方程模拟“类升力体”外形在高超声速下的大攻角流动,给出了“类升力体”外形表面极限流线随攻角变化的拓扑结构及40°攻角下垂直于体轴的横截面流线拓扑结构。结果表明:类升力体外形三维流场结构十分复杂,攻角从0°～40°变化时,背风面表面极限流线依次由不分离、开式分离向起始于鞍、结点组合的高阶奇点的分离方式转化,翼面横向分离亦随攻角增大而增大;垂直于体轴的横截面流动拓扑结构与文献[2]给出的理论分析一致。     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

基于基元反应模型的H2-O2-N2爆轰数值模拟

 采用12组分、23个化学反应的基元化学反应模型,用5阶加权本质无震荡格式（WENO）、3阶TVD Runge-Kutta格式,对H₂-O₂-N₂混合气体胞格爆轰进行了数值模拟。研究了一维ZND爆轰、自维持爆轰的详细结构以及三波点附近的流动结构。计算结果表明：由横波的压力可以显著促进二维爆轰波波阵面的形成;横波的运动生成三波点,三波点造成了爆轰的自维持传播。     

带三次非线性项的四阶Schr(o)dinger方程的分裂多辛算法

对一类带三次非线性项的四阶Schr(o)dinger方程提出分裂多辛格式.其基本思想是将多辛算法和分裂方法相结合,既具有多辛格式固有的保多辛几何结构的特性,又发挥了分裂方法在计算上灵活高效的特点.数值实验结果表明,分裂多辛格式比其它传统的多辛格式更节约计算时间和计算机的内存,从而更加优越.