带装药裂纹发动机内流场数值研究

在N－S方程基础上，考虑网格移动，建立了适用于固体火箭发动机内流场的湍流控制方程组，并对带装药裂纹的固体的火箭发动机内流场进行了数值模拟，分析了推进剂中裂纹深度，宽度、位置、角度等多种因素对发动机内流场的影响，计算结果表明：（1）裂纹出口处流速高，大于主通道流速，在裂纹出口附近存在回流区；（2）当裂纹紧靠发动机前封头时，裂纹出口附近回流强度减弱，裂纹对发动机内的流动影响较小；（3）当裂纹深度与裂纹宽度比大于240时，裂纹内压强急剧升高，对发动机装药结构完整性具有重要影响。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采用了自适应网格技术。     

固体火箭发动机两相内流场的连续介质-离散颗粒耦合混合模型数值模拟

为解决含Al推进剂固体火箭发动机内流场两相流动的数值模拟问题,采用将颗粒相视为连续介质和离散颗粒相结合的综合方法,建立起气相与颗粒相的双向耦合混合模型。针对颗粒类型与粒径尺寸的不同,采用拉格朗日法描述大粒径颗粒,平衡欧拉法描述小粒径颗粒,克服了现有模型难以全面考虑颗粒尺寸效应而使模拟精度下降的困难,通过算例验证了该混合模型的有效性和准确性。针对含Al推进剂固体火箭发动机内流场湍流气粒两相流动进行了数值模拟,分析了不同粒径尺寸颗粒的分布及其对发动机内流场和结构性能的影响情况。结果表明,大颗粒粒径在跨喷管段变化明显,平均减小30%。粒径40μm以上的颗粒易破碎,且燃烧效率进入平台区,较10μm颗粒下降50%以上,其Al含量均大于60%。     

模型固体火箭发动机内涡脱落过程的大涡模拟

为研究后向台阶结构引起的涡脱落对固体火箭发动机内流动稳定性的影响，采用大涡模拟技术，对模型发动机内的流动进行了数值模拟，获取了发动机燃烧室内压强随时间的变化曲线及其振荡频率，并和已有的计算结果及实验结果进行了对比分析。结果表明采用LES方法的计算结果明显优于采用k-ε湍流模型的计算结果；不考虑燃烧时模型发动机内由涡脱落引起的压强振荡的振幅较小，但频率较高，且随来流速度的增加而增大；在后向台阶这种结构中，涡一般会在台阶的边沿形成，并在其下游周期性地脱落；在保证y^＋较小的情况下，网格大小对计算结果的影响不明显。     

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

导弹助推器分离过程数值模拟研究

应用了结构网格中的Chimera重叠网格技术和Favre平均三维N-S方程以获得流场解.采用k-ε湍流模型模拟气体的湍流粘性影响,计算中分别考虑空气和两种火箭发动机喷流等三种不同流动介质,采用时间相关边界模拟发动机拖尾段的非定常流动,最终求解带约束的六自由度弹道方程模拟了导弹助推器的分离脱落过程.并对发动机喷流对助推器分离的影响开展研究.所做工作可对于精确确定火箭助推器分离轨迹及姿态提供方法参考.      

子午扩张对涡轮内旋涡结构影响的数值模拟

求解了某型航空发动机的低压涡轮导向器内流场。采用的数值方法为具有 TVD性质的三阶精度GODU NOV格式,湍流模型为 B-L代数模型及 MML模型的混合修正模型。该涡轮导向器的一个显著特点是由于结构上的要求,静叶前子午流道在上端壁具有较大的扩张角。通过数值模拟及流场内涡系结构和损失的分析表明,由于子午扩张造成端壁附面层增厚并在静叶前形成较大的分离区。该分离区的存在不但使静叶前损失急剧增加而且使上端壁马蹄涡和通道涡尺度加大,强度增强,引起流道内损失增加。      

基于大涡模拟分析气膜冷却的湍流场

用大涡模拟考察了单孔平板气膜冷却的湍流场,气膜孔沿流向倾斜30°,气膜出流的雷诺数为2600,吹风比为0.5.计算结果表明:①气膜孔内的流动分离增加了湍流动能,导致气膜冷却效率的降低,在设计中要尽量避免或减少流动分离;②与射流侧向扩展有关的涡黏性系数在气膜孔两侧存在峰值,而采用各向同性的湍流模型预测气膜冷却时,涡黏性系数的峰值出现在射流与主流的剪切区,因此会高估射流的垂直穿透而低估射流的侧向扩展;③可以用大涡模拟辅助建立各向异性的湍流模型,以便提高湍流模型的预测精度.      

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

固体火箭发动机不稳定燃烧时平均压力下降的原因探索

本文提出声湍流概念,指出声湍流引起的能量输运会降低经受适当振幅的推进剂的平均燃速,从而解释了不稳定燃烧推进剂平均压力的下降.本文提出的理论与Morita的实验结果在定性上是一致的,而实际发动机计算与实验结果的比较也证实了本文提出理论的正确性.     

固体火箭发动机射流中存在障碍物时的流场分析

基于N-S控制方程和RNGk-ε湍流模型,采用FLUENT软件对某固体火箭发动机存在球形障碍物的射流流场进行了数值模拟,得到了球形障碍物尺寸及距离喷管出口位置对射流流场的影响规律。计算结果表明,当障碍物处于第一个膨胀压缩过程之后与第二个膨胀压缩波过程之前的特殊位置时,它对整个射流火焰宽度的影响最为显著,该计算结果与某次发动机试车失败时所观测的现象相符。     

固体火箭冲压发动机燃气流量调节特性

为了研究固体火箭冲压发动机燃气流量调节特性,对固体火箭冲压发动机燃气流量控制阀的5种不同开度下的二维流场进行了数值模拟,分析得出了补燃室压力、阀门开度对流量控制阀流场分布、阀头轴向受力和通过控制阀的燃气质量流量的影响,基于模拟得到的数据,分析得到了固体火箭冲压发动机燃气流量调节特性的参数表达式,该表达式可以用来作为固体火箭冲压发动机燃气流量调节的参考。     

固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

旋转射流对含硼固体火箭冲压发动机二次燃烧的影响

为提高固体火箭冲压发动机二次燃烧效率,将旋转射流技术引入固体火箭冲压发动机设计,采用Re-alizable k-ε湍流模型、单步涡团耗散燃烧模型以及KING硼粒子点火和燃烧模型,利用Fluent软件开展了旋转进气和一次燃气旋转含硼固体火箭冲压发动机补燃室三维反应流场流数值分析。研究结果表明,当空气射流切向进入补燃室时,气流产生的旋转均使燃料与空气的混合更充分,燃烧效率更高。当气流切入角度增大时,补燃效率先升后降,对于具体发动机结构,存在一个使燃烧效率最大的切入角,针对研究的模型发动机结构,此值在20°附近;当一次燃气旋流数的增加,二次燃烧效率呈逐渐增高的趋势。     

k-ε（RNG）、LES模拟建筑物周围气流特征及湍流扩散的风洞试验验证

采用k-ε（RNG）与LES湍流模型在来流与建筑物迎风侧呈不同角度的情形下,模拟了位于立方体建筑物顶部污染源所排放污染物的流动和扩散规律,并与相应的风洞试验结果进行了比较。流场分析结果表明：数值模拟能够较好地模拟建筑物顶部回流、背风侧空腔区以及再附着点等。浓度场分析结果表明：来流与建筑物成45°时,建筑物顶部回流区与背风侧空腔区的数值模拟结果略低于风洞试验结果;来流与建筑物成90°时,建筑物顶部回流区数值模拟结果略高于风洞试验结果,而背风侧空腔区的数值模拟结果与风洞试验结果基本一致。综合分析表明：建筑物周围的流场影响浓度场的分布,LES、k-ε（RNG）模型都能够较好地模拟建筑物周围的流动和扩散规律,两种模型相比,LES模型与风洞试验吻合得更好。总之,风洞试验和数值模拟相结合能较好地研究建筑物对流动和扩散的影响。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

基于人工神经网络的可压缩湍流大涡模拟模型

在国家数值风洞（NNW）工程项目的指导下，空间人工神经网络（SANN）模型被用于强可压缩湍流大涡模拟（LES）研究，其中流场的湍流马赫数分别为0.6、0.8、1.0。基于湍流的多尺度空间结构特性和人工神经网络方法发展的高精度空间神经网络（SANN）模型适用于不可压缩湍流和弱可压缩湍流。对于强可压缩湍流，流场中会出现激波结构，给大涡模拟带来了挑战。本文的研究结果表明：SANN模型适用于强可压缩湍流的大涡模拟。在先验分析中，SANN模型预测的亚格子应力和亚格子热流的相关系数超过0.995，远远高于梯度模型和近似反卷积模型等传统模型；传统模型的相对误差大于30%，而SANN模型在这方面有很大的改进，相对误差低于11%。在后验分析中，与隐式大涡模拟（ILES）、动态Smagorinsky模型（DSM）、动态混合模型（DMM）相比，SANN模型能更精确地预测能谱、各类湍流统计特性以及瞬态流动结构。因此，基于湍流多尺度空间结构特性的人工神经网络模型加深了人们对强可压缩湍流亚格子建模的认识，同时可以服务于NNW工程的流体力学模型构造。     

固液发动机实验燃烧器的气固耦合传热计算

针对固液火箭发动机中复杂的燃烧、传热过程,在FLUENT软件的平台上,采用气固耦合传热计算的方法模拟了二维扩散燃烧实验器中的温度、组分分布,并利用动网格方法模拟了固体燃面的不规则热解退移。稳态流场计算的结果和燃面动态退移情况都验证了扩散火焰距燃面的距离和燃面上方主流流速是影响退移速率的主要因素。此方法可以较好地模拟固液火箭发动机中复杂的非定常传热和燃面退移过程,较准确地预示固液发动机的内弹道性能。