非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧,提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步,在粗网格上采用Oseen迭代法求解非线性问题,在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散,时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.     

基于激光雷达的垂直能见度反演算法及其误差评估

针对大气垂直方向上消光系数分布不均匀难以用传统方法直接测量垂直能见度的问题,提出了一种基于激光雷达探测垂直能见度的计算方法。根据大气辐射传输基本原理,借助于辐射传输方程,推导出了垂直能见度的计算公式;然后利用激光雷达原理方程和Klett算法反演出大气垂直方向上的消光系数分布,基于此提出了垂直能见度的迭代算法。最后,利用灰色模型GM(1,1)和批统计算法,对激光雷达反演得到的后向散射系数进行了评估,给出了误差置信区间为（0.760±0.339）×10⁻⁴（srad·km）⁻¹。结果表明,该方法是一种特别有效的计算垂直能见度的方法,符合探测的基本需求,且误差小精度高。     

基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

暗怪波的相互作用

以耦合非线性薛定谔方程为理论模型,数值研究了两个一阶暗怪波在正常色散单模光纤中的相互作用.基于一阶暗怪波精确解,采用分步傅里叶数值模拟法,从间距、相位差和振幅系数比方面讨论相邻两个一阶暗怪波之间的相互作用.基于二阶暗怪波精确解,讨论了两个一阶暗怪波的非线性相互作用.研究结果表明:同相位情况下,间距参数T1为0、5、20时,相邻两个一阶暗怪波相互作用激发产生“扭结型”暗怪波.相比较于单个暗怪波发生能量的弥散,“扭结型”暗怪波分裂形成多个次暗怪波.反相位情况下,间距参数T1为2、7、12时,相邻两个一阶暗怪波相互作用也可以激发产生“扭结型”暗怪波.并且“扭结型”暗怪波初始激发的空间位置偏离原始单个暗怪波的位置5.振幅系数比越大,该空间位置越接近5.二阶暗怪波可以看作是两个一阶暗怪波的非线性叠加,复合型和三组分型二阶暗怪波与相邻两个一阶暗怪波的相互作用略有相似.     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

随机薛定谔级联方程：理论和应用

随机薛定谔级联方程是一种基于波函数的严格量子动力学方法,它可用于研究耦合到玻色子热库的复杂体系中的量子动力学过程. 本综述从开放量子体系费曼路径积分的影响泛函出发,概述了随机薛定谔级联方程的一般理论框架和各种具体表述形式,并通过对复杂体系中超快激发能量转移过程的模拟来展示方法的应用范例和计算效率.     

对薄层柱壳爆炸膨胀断裂过程的研究

 本文提出了一个用于描述动态破坏发展过程的损伤度函数。从这个损伤度函数出发,把材料特征性方程取为强化粘塑性本构方程形式,导出了薄层柱壳爆炸膨胀运动在两种近似下（恒定膨胀速度近似合恒定应变速率近似）断裂判据的解析表达式。结果分析表明,在上述条件下,存在着一个动态断裂“塑性峰”,在这个峰值条件的应变率下,柱壳出现贯穿断裂时刻的应变最大。以软钢为算例,本断裂判据可以比较好地解释Иванов和陈大年等给出的实验结果。这时,动态断裂“塑性峰”对应的应变率为4×10⁴ s^-1,相应的应变约为60%~80%。     

原子分子光子系统的耗散相互作用和退相干

原子分子系统与量子化的电磁场或光子模式耦合的系统是非相对论量子力学理论研究和实验研究的主要对象和模型. 现实系统必然与外界环境耦合,且即便原子隔绝较好、光学腔壁品质因子足够高,原子系统也不等价于少数几个能级构成的简单模型：它仍然有不为零的几率跃迁到不可控的能级空间、与原子相互作用的自由空间真空场的量子效应也必须考虑. 本文将结合开放量子系统理论的基本要素与原子光子的基本模型,对原子分子系统在电磁场中发生的耗散以及量子退相干过程做简单综述,并重点介绍描述量子系统退相干过程的主流理论工具——主方程.     

受匀速转动约束的变质量非完整力学系统的Appell方程

本文建立受匀速转动约束的变质量非完整力学系统l...     

从相关测量的数据恢复超短激光脉冲的波形

本文证明了,只要测量了脉冲光强的二阶相关函数G~(2)(τ)和单延迟的三阶相关函数G~(3)(τ),便可以恢复超短激光脉冲的强度波形。并提出一种恢复脉冲波形的计算方法,编制了计算机程序,对实验数据作了处理,证实这种计算方法是可行的。在此基础上建议一种通过相关测量确定超短激光脉冲波形的新型仪器。