畸变进气下压气机三维黏性数值计算模型

通过对压气机环形无叶片通道求解带源项的N-S(Navier-Stokes)方程,实现对压气机全通道流场的计算,建立了压气机三维黏性非定常数值计算模型.源项是由三维计算流体动力学(CFD)求解器定常计算的气流参数计算得到,求得不同工作点源项大小,然后根据进口参数对源项插值得到当前计算点源项.采用该模型对某跨声速单级压气机在均匀进气下性能和流场的计算结果与实验结果进行了对比,两者能够较好的一致,验证了模型的可靠性;之后应用该模型对周向总压畸变进气下的流场和性能进行了计算,分析了畸变进气对压气机三维流场和性能的影响,并引入压气机失稳判据分析畸变对压气机稳定性的影响.      

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验     

风扇进气畸变三维非定常数值模拟技术研究

采用三维欧拉方程加体积力模型的方法, 发展了一套三维非定常进气畸变数值模拟程序, 以实现风扇/压气机在进气畸变条件下的全环面三维非定常数值模拟.利用该程序模拟计算了某航空发动机风扇进口稳态总压畸变条件下特性和流场结构, 并与试验结果进行了对比, 结果表明该程序用较小的计算机资源即可以实现对风扇全环面三维非定场畸变流场的模拟, 能够以较高的精度预估气流畸变对压气机/风扇特性的影响.      

单级跨声速压气机进气畸变数值模拟

基于彻体力模型的基本思想,发展了一套能够反映进气畸变对风扇/压气机气动性能影响的三维数值计算程序CSAC.首先详细研究了彻体力模型的具体建模方法,然后发展了相应的三维数值计算程序.利用该程序针对某单级跨声速压气机在均匀进气下的流场进行了计算验证,计算结果与雷诺平均Navier-Stokes(RANS)计算吻合得很好,流场主要参数分布与相关实验数据也有很好的一致性.利用该程序分别计算分析了进口存在稳态周向总压畸变及周向总温畸变的全环三维畸变流场,结果表明:压气机进口畸变流场存在大尺度、强三维特征,将导致压气机总压比及稳定裕度下降.所有计算结果均显示该模型能够在降低对计算资源和工程经验依赖的同时,有效分析进气畸变对压气机特性的影响.      

风扇/压气机稳定性三维模型的研究

发展了一个基于体积力的可用于分析风扇/压气机性能和稳定性的三维可压缩模型。该模型先是利用三维CFD求解器求得在同一转速线下不同工作状态下的每一叶片排的叶轮机源项,然后通过求解带源项的三维非定常Eu ler方程,获得对风扇/压气机内部三维流场和性能的模拟。利用该模型对某一跨声速风扇级的三维流场和性能进行了数值模拟分析,特别是分析对比了在进口无畸变和进口有畸变情况下的风扇内部三维流场和气动性能。研究结果表明,所发展的计算模型的数值结果与实验结果吻合得很好。与实验结果相比,预测总温比的误差为1.2%,而预测总压比的误差小于4%。     

固体火箭发动机燃烧室三维流场数值模拟

利用ＳＩＭＰＬＥ方法数值求解层流Ｎ－Ｓ方程，模拟了具有移动边界的固体火箭发动机燃烧室内的三维流场。网格划分采用交错网格技术，流场计算时不进行网格划分，而是直接读取网格点的值，将网格划分与流场计算分开，解决了现有计算机内存不够的问题，计算获得了满意结果。     

压缩性修正对γ-Re_θ转捩模型的影响研究

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

波瓣混合流场流向涡的数值分析

借助N-S方程对波瓣喷管混合装置的流场进行了数值分析。采用三维贴体曲线坐标,非交错网格法和K－ε两方程湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段流场统一进行了数值计算。得出了速度场和流向涡及正交涡的分布,并对流向涡的发展规律及其对正交涡分布规律的影响作了分析。并对混合管壁面的压力分布和出口速度型与实验值进行了比较。     

飞船高超声速粘性绕流的数值模拟

本文采用稳式ＮＮＤ格式，通过求解三维薄层近似Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了飞船高超声速粘性流场。在计算中充分考虑了网格雷诺数和网格分布对计算结果的影响，不但得到了满意的表面压力、热流分布和流场的物理量分布，而且得到与实验观察一致的表面油流图像和分离结构。     

基于分区的隐式求解二维不可压NS方程的并行实现

本文通过对目前常用并行算法的分析,基于区域分解的思想,给出了适用于隐式求解原始变量二维不可压NS方程压力修正方法的并行算法,并在网络异构编程PVM环境下进行了数值验证,讨论了计算结果,结果显示对任务级的并行,并行效率也同样随粒度的增加而提高.另外此方法还可推广到三维复杂流场的计算问题.     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

超燃冲压发动机三维喷流干扰流场的数值研究

用隐式有限体积法求解三维N－S方程数值研究超燃冲压发动机三维平行、垂直喷流干扰流场。计算中采用了具有当地极值递减（LED）性质的JST和SLIP中心差分格式。为了减缓计算网格数很多时对计算机内存的压力以及提高计算效率，运用了通量守恒的分区计算方法。     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

三维进气道粘性流场数值模拟

对三维进气道流场进行了数值模拟，将三维进口激波系的数值模拟结果与二维模拟及激波理论估算值进行了对比。说明可以用现有的计算方法对三维进口波系进行评估。对三维进气道内管道流场的模拟由于受诸多因素的影响，与试验结果相比会有一定差别，本文对此进行了分析。同时对计算的边界条件及收敛特性也进行了分析。     

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟.数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象.本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响.     

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

超声速H_2/Air湍流扩散燃烧RANS数值模拟

为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。     

烧蚀外形三维粘性流场数值模拟

本文通过求解层流Navier-Stokes方程,数值模拟了三维奶头型烧蚀外形粘性流场。差分格式为Yee和Harten的迎风隐式TVD格式,在空间方向具有二阶精度。计算得到了高质量的激波和合理的流场结构。物面压力分布与实验结果进行比较,符合较好。