液体火箭发动机喷管壁面辐射热流的数值计算

为火箭发动机喷管设计及其冷却系统的设计提供有效参考数据作为应用背景,计算喷管壁所受的红外辐射热流。分别推导了有限体积法和反向蒙特卡罗法适用于不含粒子条件下喷管壁面所受红外辐射热流的计算公式,并针对不同的边界条件设定,重点分析了有限体积法的特性,并提出了有限体积法的局限性。最后,以某三维双圆弧喷管模型为例,用有限体积法和反向蒙特卡罗法自主编写程序计算在不同边界设定下壁面受到的红外辐射热流密度。通过计算结果对有限体积法的局限性进行了论证,同时肯定了反向蒙特卡罗法在壁面辐射热流计算上的优势。     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

用于复杂流场数值模拟的重叠／嵌套网格方法

介绍了一类用于复杂流场数值模拟的重叠／嵌套网格方法。它是一类多块结构化网格,可用于有或无多重网格加速收敛技术的有限体积、有限差分等计算。其基本思想是把复杂流场的计算区域化分成相对简单的子区域,每个子区域上的网格可以独立生成;网格界面之间的流场信息交换通过三线性插值来实现,通过修改流场解算方法来排除洞中的点。     

一种适用于进气道/风扇一体化计算的体积力模型

全流道的进气道/风扇一体化计算虽然能捕捉到全面、完整的流场细节,但仍然需要消耗大量的计算资源。为了节约计算成本,利用有限的设备达到快速高精度计算的目的,采用自行发展的块体积力模型替代真实叶片进行一体化计算。该模型将叶片流域沿周向等距分成若干块,每个块的体积力源项关联于叶片前缘参数。研究结果表明:该体积力模型能够准确地模拟出流场特征与细节。在均匀来流下,各参数的相对误差在0.5%以内。同时,在大S弯进气道高畸变来流条件下,与冻结转子方法计算结果相比,总压比、总温比和等熵效率的相对误差分别为4.49%,0.26%和2.38%。     

预测纳米多孔铜弹性模量的随机分布单胞模型

纳米多孔铜的弹性模量实测值远低于分子动力学的模拟结果,且实际韧带尺寸远大于分子动力学模拟结果。通过Python平台在ABAQUS中构建纳米多孔材料的随机分布单胞模型,基于热应力比拟方法,用多尺度均匀化方法求出均匀化弹性参数。首先,利用所提模型预测纳米多孔金的弹性模量并与实测结果进行了对比,两者的高度吻合说明所提模型的有效性。其次,预测并分析了不同体积分数下纳米多孔铜的均匀化弹性模量,发现了纳米多孔铜均匀化弹性模量随着体积分数变化的阈值现象,并从力学角度对其机理进行了解释,分析了预测结果大于实测结果的影响因素。      

三维非结构混合网格生成与N-S方程求解

 提出了一套较为通用的, 完全自动化的非结构混合网格生成方法, 在物面粘性作用区, 运用一种改进的推进层方法生成三棱柱形和金字塔形网格; 在其他流动区域采用阵面推进方法生成四面体网格。采用一种改进精度的格心有限体积法和Spalart-Allmaras 一方程湍流模型对三维N-S 方程进行了求解, 并以M6 机翼、DLR-F4 翼身组合体和某运输机外形的粘性流场作为数值算例, 验证了上述网格生成和流场求解方法的正确性和实用性。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

飞行器复杂外形流场NS方程数值模拟

本文采用分块对接和“O”型网格技术，首先在飞行器表面长出“一层”分块结构网格来，再将计算区域用若干分块结构网格填充，这样既保证了网格的质量和模拟复杂外形的能力，又大大提高了计算效率。在计算中，采用有限体积法，MuScL型迎风差分格式求解Ns方程，根据多块网格的排序，研究了计算方法。最后，对某飞行器外形进行了亚跨超流场数值计算，计算结果与试验符合较好。     

高超声速升力体气动力气动热数值计算

升力体气动布局，由于其良好的高超声速气动特性、有效的内部空间和有益的防热特性，越来越受到气动工作者的重视，本文应用张涵信院士的NND格式思想，用数值模拟的方法对类似于美国X－33的升力体外形进行了气动力、气动热研究。特别是在粘性项的离散和表面热流率的计算中应用了散度和梯度的积分定义，避免了数值计算的奇异性，保证了通量守恒，提高了计算效率。结果表明气动热计算结果与美国同类飞行器结果接近，驻点热流率与经典理论结果一致。