筒状大长度固体药柱在内压载荷下的响应

本文研究不可压缩的粘弹性圆筒，在轴对称的平面应变条件下承受内压的问题，粘弹性材料采用分数阶导数模型。并用分数阶导数模型拟合固体推进剂松弛模量曲线，用Ｌａｐｌａｃｅ变换和其反变换计算了固体火箭发动机筒段在点火内压建立过程中的位移和应力，得到它们随位置和时间变化的曲线。     

二维粘性不可压缩流场的分块耦合求解方法

本文用分块耦合的方法，并采用S.Abdallah^[1,2]和M.L.Mansour^[3]的在非交错网络上隐式求解原始变量的Navier-Stokes方程的提法和离散格式，成功地解决了二维复杂流场的数值求解问题，获得了空腔流、T型叉管流、NACA0020翼型绕流和方块钝体绕流的稳态解，并进行了讨论。此方法还可推广到三维复杂流场的计算问题。     

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

有攻角细长旋成体极端地面效应分析

本文用匹配摄动方法研究了细长旋成体贴近地面运动时诱导的不可压缩势流问题,它是经典细长体理论的推广和发展。文中将流场分为内外两个渐近区域,建立了问题的匹配摄动理论,并利用文献[1]中的双圆诱导不可压缩势流的精确解,给出了该问题的首阶解析解。然后又导出了细长体上的吸力和俯仰力矩的计算公式。在此基础上,文中对细长旋成体极端地面效应的运动学规律和动力学规律作了定性和定量两个方面的分析。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。     

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

二维不可压流程序在PVM环境下的并行计算

本文在PVM环境下进行了二维不可压流的并行数值模拟,应用并行虚拟机模型对机器环境和并行程序的特性进行了分析,研究了并行程序和并行机器的匹配问题.在MPP环境下的测试表明,改进后的并行程序并行效率高.