97．5．8华南地区暴雨过程的数值模拟

用非静力平衡的中尺度数值模式MM5V3版，模拟研究了1997年5月8日OOUTC至5月9日OOUTC，发生在华南沿海地区的暴雨过程。结果表明：该模式比较成功地模拟预报了产生该次暴雨的环流形势、急流位置及其强度、暴雨的范围、雨量中心。通过分析数值试验的结果表明：凝结潜热释放的加热作用对暴雨产生和华南低空气流的维持是十分重要的。     

中心分级燃烧室贫油熄火试验研究与数值模拟

为研究三级旋流器下中心分级燃烧室贫油熄火性能，常压下针对在不同进气速度以及旋流器几何参数下中心分级燃烧室贫油熄火特性进行了试验研究与数值模拟，获得了在一定条件下的贫油熄火规律。试验结果表明：随着进气速度增大到35m/s时,贫油熄火油气比减小至0.0070,说明进气速度对贫油熄火起重要作用;值班级旋向相反时贫油熄火油气比(A方案为0.0066，B方案为0.0071)均小于同向时(C方案为0.0078，D方案为0.0074)，且A方案最小，说明值班级旋向对贫油熄火性能有明显的影响，值班级反向时熄火特性最好，且主燃级旋向对其作用不大，最佳旋向组合为A方案，即值班级1，2级叶片旋向反向，值班级2级叶片与主燃级叶片旋向相同；值班级1，2级叶片角度增大贫油熄火油气比分别下降了将近4%和9%，其熄火特性变好。数值模拟计算结果与试验结果熄火规律变化一致。     

高强化活塞振荡冷却的场协同分析

数值模拟了高强化活塞内冷油腔中机油的振荡冷却过程,定义了有效速度的概念,通过数据处理得到了内冷油腔壁面的各场参数在时间和空间上的分布规律,分析了壁面表面传热系数与场协同性能间的联系.结果表明:在面相位为60°~120°和240°~300°范围内,场协同匹配性能较好;有效速度的分布规律与协同角余弦值的分布规律一致;内冷油腔壁面表面传热系数分布规律与协同角余弦值和有效速度的分布规律一致,表明可以使用场协同机理阐述机油振荡冷却的对流换热规律.      

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

青岛冷涡雷暴分析

通过对2003年6月13日青岛地区冷涡雷暴天气过程从大尺度环流活动、区域地面流场、物理量场、卫星云图等方面进行了诊断分析，结果表明：贝加尔阻塞高压的稳定维持，是大尺度环流的主要特征，截断低压形成的东北冷涡是触发雷暴的天气尺度系统，Ekman边界层中正涡度的动力强迫对垂直运动和水汽的垂直输送起着主导作用。而雷暴的空间分布又具有明显的中尺度特征，中尺度云团和中尺度雨团，它们与物理量的有利配置是产生此次雷暴的根本原因。     

倾斜单孔射流轨迹与速度分布的研究

本文采用实验和数值模拟的方法,在吹风比分别为0．5、1．0、1．5、2．0的情况下,研究了圆柱单孔二次射流的贴壁性及速度分布。通过对比不同截面位置的流动轨迹和速度分布,发现圆柱孔射流流场的速度分布与圆柱扰流流场的类似。射流出口的最大速度并不位于射流孔中心位置,u值也不是在射流孔中心位置最高。     

地面粗糙度对下击暴流风剖面特征影响

地面粗糙度是大气边界层中反映下垫面形态的重要指标,也是影响近地风场特征的重要因素。为研究地面粗糙度对下击暴流风剖面特性的影响,基于计算流体力学方法建立了下击暴流三维足尺模型,通过实验对数值模型进行了验证。通过调整粗糙元高度及分布密度来模拟自然界地面不同的粗糙类别,数值模拟了具有不同地面粗糙长度的下击暴流近地风场。结果表明:在距离风暴中心较近的位置(r≤1.0Djet),地面粗糙度对下击暴流风场的影响并不明显,各径向位置的最大风速值和最大风速所在高度都基本不受地面粗糙度影响;在下击暴流冲击地面后沿径向发展的过程中,经过粗糙的地面,产生能量耗散效应,地面粗糙度对于下击暴流风剖面特征的影响逐渐显著,不同地貌下的竖直风剖面产生较大差异;在近地面高度,地面粗糙度对下击暴流径向风剖面影响显著且影响范围大,沿径向发展的方向下击暴流的风速随着地面粗糙长度的增加而下降更迅速;随着距离地面高度增加,地面粗糙度对径向风剖面的影响主要体现在远离风暴中心的区域。尽管下击暴流形成后沿径向扩散过程存在强度自然衰减过程,但在距离风暴中心较远的径向位置,当遭遇强下击暴流时,下击暴流引起的近地面强风仍然具有很大的威胁和破坏性,因此地面粗糙度对下击暴流风剖面特性的影响不能忽略,需要在风剖面模型中考虑地面粗糙度的修正。     

民机风洞试验半模垫板高度对气动特性的影响

半模作为提高大型商用飞机风洞试验雷诺数的一种模拟手段而被广泛应用。首先回顾了半模试验的模拟方式及其优劣,进而选取当前发展趋势的附面层垫板作为研究对象,采用数值模拟研究了垫板高度变化对气动特性影响的内在机理。数值模拟结果和试验吻合较好,数值计算采用速度分布入口可以较好模拟风洞核心段边界层厚度,计算值和试验值更加接近;垫板高度的增加使得升力系数增加、阻力系数减小及俯仰力矩系数增加;垫板在机翼上游区引起的上洗使得机翼沿展向各剖面当地迎角增加5%、动压增加1%,从而使得机翼上翼面压力分布朝负值方向移动。有别于以往认为垫板的洗流只影响内侧机翼,结果表明垫板影响范围扩展至全翼展,当地迎角的增加是主要影响因素,垫板对机翼展向各剖面影响量值不一致,对内侧机翼影响较大。所得结论可更好用于民机半模风洞试验的开展,具有一定的工程实用性。     

AECSC-J ASMIN湍流燃烧仿真软件研发和检验

航空发动机燃烧室几何结构复杂,湍流和化学反应存在强烈非线性相互作用,需要对流动和燃烧及其相互作用进行高精度高时空分辨率的刻画,目前燃烧室湍流燃烧数值模拟仍然是高难度的瓶颈问题之一。介绍了由北京航空航天大学航空发动机数值仿真研究中心、北京应用物理与计算数学研究所和中国工程物理研究院高性能数值模拟软件中心联合研发的AECSC-JASMIN软件主要框架、算法以及针对该软件的算例检验。在Sandia射流火焰、支板火焰和单头部燃烧室检验算例中,对比实验数据,射流和支板火焰预测结果与实验值一致;支板算例的平均相对误差在15%之内;单头部燃烧室模拟结果符合物理实际,总压损失与实验值基本吻合。说明AECSC-JASMIN软件可用于复杂结构高分辨率高精度湍流燃烧数值模拟。     

限带白噪声过程雨流幅值分布

为了在频域内解决限带白噪声载荷作用下结构的疲劳寿命估算问题,提出了限带白噪声过程的雨流幅值分布模型,并给出了参数估算方法。假设雨流幅值可以用双参数Weibull分布来描述,形状参数与带宽系数呈线性关系;尺度参数与均方根呈比例关系;比例因子与不规则因子呈线性关系。数值模拟结果表明：Weibull分布能够很好地描述限带白噪声过程的雨流幅值分布,分布参数能够根据随机过程的谱参数计算得出。     

随机场TPDF在非结构网格平台的开发及检验

为实现对燃烧室内部燃烧过程的高精度数值模拟,增强航空发动机燃烧室研究手段,采用开源的非结构网格平台Saturne,开发可以耦合详细化学反应机理,高精度解析湍流燃烧过程的概率密度函数输运方程（TPDF）湍流燃烧模型。采用随机场方法求解TPDF方程,在原有程序基础上研发了TPDF程序模块及配合的加速算法等模拟单元,发展了针对燃烧室燃烧过程模拟的软件功能。采用射流火焰Flame D和旋流火焰TECFLAM对新软件进行了测试,结果表明：模拟结果与试验数据趋势一致,精度尚可。对某型燃烧室性能模拟验证中,计算的出口温度分布与试验结果一致,平均温度误差小于3.8%,新的软件可进一步应用于燃烧室性能的研究。     

C-H-O-N体系液体火箭发动机喷管流场数值模拟

采用弱耦合点隐式方法的MacCormack格式对C－H－O－N体系的液体火箭发动机喷管粘性流场进行了数值模拟，采用了考虑化学反应和不考虑化学反应的两种模型，其中化学动力学模型为C－H－O－N体系的12组分、14反应有限速率的化学反应模型。数值模拟得到了流场参数在喷管中的分布，结果表明，数值模拟结果与理论分析一致，为发动机的设计及试验提供了参考。     

密封气流流量对轴承腔外壁滑油运动的影响

为了研究密封气流质量流量对轴承腔外壁面上滑油分布及温度分布的影响,基于拉格朗日离散相模型（DPM）和液膜模型,利用STAR CCM+商业软件对轴承腔内油滴运动、油滴向油膜的转化以及油膜在轴承腔内的运动开展非稳态数值模拟计算,并与德国Karlsruhe大学轴承腔油膜厚度试验结果进行对比。结果表明：计算值与试验值在高转速下一致性较好,平均相对误差13.6%;低转速下由于空气分布均匀性较差,个别工况点计算值与试验值存在一定误差,但最小相对误差在4%,总体上具有较好一致性;随着轴承腔密封气流的质量流量的增加,轴承腔内空气的平均流动速度提高,对轴承腔外壁面上油膜的剪切作用和扰动能力增强,导致轴承腔的外壁面油膜厚度和稳定性降低;轴承腔外壁面上的温度分布与油膜的厚度分布一致,最小温度分布在壁面上被油滴冲击位置;随着轴承腔密封气流的质量流量的增加,轴承腔外壁面上低温分布范围增大,但温度分布均匀性变差。     

常规装药近地爆炸冲击波效应数值模拟

针对现有经验公式在计算爆炸地面冲击波超压值时存在的局限,采用NND2M差分格式对TNT柱形装药在小爆高状态下的爆炸冲击波效应进行了数值模拟研究,并与试验实测值及经验公式计算值进行了比较,分析了存在差异的原因,证明了本文采用的计算是合理可靠的.     

混合相/冰晶条件下的结冰研究综述

针对混合相/冰晶结冰问题，阐述了与过冷水结冰的差异及其危害。从混合相/冰晶结冰地面模拟试验设备、试验现象以及数值模拟3个方面详细论述了混合相/冰晶结冰的研究进展。其中，地面模拟试验设备以Cox冰风洞、NRC高空试验设备、NASA推进系统试验室、德国布伦瑞克冰风洞为例，总结了各自的冰晶生成方法和试验能力；分别介绍了以上4个地面模拟设备中开展的混合相/冰晶试验，分析了其试验现象的异同和产生的原因；针对混合相/冰晶结冰数值模拟中面临的4个关键问题：冰晶运动融化相变、冰晶黏附、侵蚀、结冰过程模拟，总结了现有模型和模拟方法，提出了不足之处。最后，对混合相/冰晶结冰有待进一步重点关注的研究方向进行了展望。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景,采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型,采用高精度WENO格式进行空间离散,并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致,来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构,并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析,揭示了流场结构特征及其时空演化规律,并获得了流场密度脉动特性。      

超高温升中心分级燃烧室设计及计算分析

针对航空发动机高推重比、高温升的需求,提出1种中心分级旋流燃烧室的设计方案。在保证与现有单环腔燃烧室(SAC)进出口尺寸、机匣尺寸限制不变的情况下,对设计模型进行了3维数值模拟,并与现有的单环腔燃烧室数值模拟结果及试验结果进行了对比分析。研究结果表明:设计油气比为0.045时,设计中心分级燃烧室温升可达1356 K,出口温度分布可达0.137,出口径向温度分布可达0.096;此外,与SAC相比,中心分级燃烧室可获得更低的总压损失,更低的出口温度分布系数以及高工况下可获得更高的燃烧效率;污染排放性能表明,中心分级燃烧室在慢车点CO排放比SAC的稍高,在设计点NOx排放按g/kg燃油计比SAC的低。     

抽拉速率对定向凝固DZ4125合金温度场及晶粒竞争生长的影响

为探究定向凝固DZ4125柱晶高温合金中温度场分布情况，并研究不同抽拉速率对温度梯度及糊状区变化的影响规律，利用ProCAST软件对抽拉速率分别为3、5、7 mm/min的定向凝固过程进行模拟。结果表明：当试棒位于保温区时，等温线呈外侧低内侧高的倾斜分布，当试棒位于冷却区时，等温线呈外侧高内侧低的倾斜分布；随着抽拉速率的增加，温度梯度逐渐减小；糊状区形状取决于抽拉速率和距水冷盘的距离，当抽拉速率为5 mm/min时，凝固前沿水平状保持的时间最长；在此基础上，对定向凝固试棒进行精铸实验，并研究3mm/min和7 mm/min下试棒液相线形状与晶粒结构之间的关系，发现液相线呈凹形和凸形时，晶粒分别向试棒中心和表面方向生长，水平的液相线能促进晶粒沿试棒轴线方向生长。     

航空发动机振动突增问题分析

针对某型航空发动机在试车过程中的稳态振动突增及高压转子轴心轨迹变化,基于动力学普遍方程及实际发动机的具 体情况,开展了故障因素排查;采用数值仿真计算方法建立了双转子发动机动力学分析模型,分析了K 5 支撑刚度对整机振动的影 响。结果表明：K 5 支撑刚度由各向同性到各向异性状态且刚度值发生变化,发动机支点的振动值将增大且轴心轨迹发生改变,仿 真结果与试车结果一致。分析认为稳态振动突增及转子轴心轨迹变化的主要原因是支撑刚度变化造成的,其主要由于发动机涡 轮后机匣在工作过程中受热变形,促使拉杆由自由向拉紧状态变化,导致K 5 支撑刚度发生变化。研究方法及分析结果对发动机整 机振动异常排故具有指导意义。     

燃气源开槽导管出口烧蚀特性分析

针对某火箭弹控制系统燃气源在点火试验过程中出现的导管出口烧蚀现象,利用数值计算模拟气流流动参数分布, 考虑到导管管壁的非稳态传热过程, 计算出不同时刻在导管开槽出口附近管壁的温度分布, 得出与试验一致的结论