多场耦合计算平台与高超声速热防护结构传热问题研究

从有限元法(FEM)和有限体积法(FVM)的单元特性出发,提出一种具有局部守恒特性的界面载荷插值方法.采用共享内存技术开发适用于通用有限元和计算流体力学(CFD)软件的多场耦合计算平台,并基于分区耦合方式实现流固耦合传热计算.作为验证,分别将ANSYS与Fluent和CFD-FASTRAN软件进行耦合,计算外壁冷却的喷...     

强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

用任意拉格朗日--欧拉方法求解流体力学方程组

从积分形式的二维Lagrange流体力学方程组出发，用有限体积格式进行计算，针对不规则的四边形网格，在用次级网格的方法确定网眼内物理量的梯度的基础上，提出一种体积加权的方法来构造网眼内的线性插值多项式。将新网格进行了细划，先确定新网格小网眼中心点的密度，将小网眼的质量直接计算出来，再计算新网格的质量，然后确定相应的密度，从而实现高精度重映。最后用ALE方法进行数值模拟，得到了预期的效果。     

复合材料层合梁分层问题中应力强度因子与守恒积分的解析变分解法

利用虚功原理导出了涉及两种不同介质的守恒积分表达式以及守恒积分J_Ⅰ与应力强度因子K_Ⅱ的关系式。对于对称正交铺层复合材料层合梁的分层问题,在利用分区广义变分原理确定了应力强度因子K_Ⅱ以及应力与位移本征展开式中的所有待定系数的基础上,进一步求出守恒积分值。最后,作为验证,由J_Ⅰ求出K_Ⅱ,并与由解析变分方法求出的K_Ⅱ相比较。计算表明,本方法前期准备工作少,计算节省机时,结果收敛迅速;沿不同积分路径求出的守恒积分保持恒定,并且用两种方法求出的K_Ⅱ具有很好的一致性。     

计算地球流体力学若干新进展

本文简要介绍了有关计算地球流体力学问题研究的若干新进展,其中包括显式平方守恒格式,经济算法,辛算子法,泥沙冲积和三角洲发育的数值模拟等,从中也可以看到计算地球流体力学的发展是具有重要理论意义和实用价值的。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响

为了研究旋转爆轰发动机燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响,采用守恒元和求解元(CE/SE)方法对带化学反应的汽油/富氧空气两相旋转爆轰理论模型进行求解,获得了气液两相旋转爆轰流场结构,并分析了燃烧室轴向和周向长度对燃烧室流场、爆轰波传播特性以及发动机推力性能的影响。计算结果表明:轴向长度对燃烧室上游流场影响甚小,却对下游流场参数影响较显著。随着轴向长度增加,燃烧室出口压力、温度、密度以及周向速度均降低,轴向速度则逐渐增大,同时发动机平均推力密度和燃料比冲先增大后减小。当周向长度过短,燃烧室内难以形成自持传播的爆轰波,随着周向长度增加,上游爆轰波强度增加,对应的流场参数均有所增大,但发动机推力性能略有所降低。     

旋转蛇形通道内冷气流动和换热的数值模拟

对于具有近真实发动机参数的涡轮叶片蛇形内冷通道,数值模拟了旋转状态下光滑、前后缘带肋及环肋结构的内流场和温度场分布,对比分析了其流动和换热特性。结果表明,带肋通道各面换热均明显好于光滑面通道;旋转对各通道进气段的影响比出气段明显;哥氏力在环肋通道中产生了最为强烈的截面二次流;光滑通道出气段换热受弯管效应的影响最为显著。      

模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.