三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

流体网络法在发动机空气冷却系统设计中的应用

将一维稳态流体网络方法应用到发动机空气冷却系统沿程压力、温度和流量预测中.该方法将空气系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,在压力修正、流热耦合方法基础上引入系数迭代修正求解.程序应用到有经验关联式的模型中验证,计算结果与经验关系式吻合得较好.验证后的程序被应用到某发动机空气系统初步设计方案的校核和改进中,在模拟结果指导下的参数调整很好地满足了设计要求.      

非稳态流体网络方法在发动机空气冷却系统中的应用

将流体网络法与非稳态模拟方法相结合应用到发动机空气冷却系统沿程温度、压力、流量的非稳态特性预测中.本方法将空气系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,针对传统解法的计算稳定性问题引入了压力修正和积分方法相结合的求解方法.程序应用到可参考对照的模型中得到验证.验证后的程序被应用到发动机轴部冷却通道空气系统非稳态特性预测中.      

离心叶栅湍流流场的数值研究

本文用二维不可压平均定常的 N- S方程组结合两方程 k-ε湍流模型 ,在任意适体曲线坐标系上计算了旋转离心叶轮叶片通道内的相对流场 ,结果与实验数据吻合良好。分析表明 ,由于离心力和哥氏力构成动量方程中的一个很大源项 ,因而较之非旋转流动需要采用较低的低松弛因子 ,并相应地提高了对初始流场和初始压力场的要求     

斜网格修正的低速计算方法及应用

采用无坐标变换的一种有限体积法对RANS方程进行离散求解,离散方法考虑了对网格非正交情况的二阶修正,通过构造一个动量插值实现了同位网格下的压力修正算法,并发展了相应的分区算法.利用以上方法对轿车绕流和三峡坝区主流摆动问题进行了数值计算,计算结果与实验吻合较好.     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

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

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

非稳态流体网络模拟新方法及其应用

推导了一种新的非稳态流体网络模拟方法应用到流体系统沿程温度、压力、流量的非稳态特性预测中.本方法将流体系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,应用压力修正和微积分方法相结合的综合方法进行模拟.程序应用到NASA的推进剂压缩模型中,并与经验关联式和NASA的程序模拟结果进行了对照,结果吻合地较好,证明了程序的可行性,验证后的程序被应用到航空发动机冷却流路模拟中.      

三种跨声速洞壁干扰修正方法及其在小展弦比飞翼标模试验中的应用

在 FL-24风洞进行了带壁压信息测量的小展弦比飞翼标模测力试验,并在 FL-26风洞进行了洞壁干扰验证试验。本文利用小扰动位流壁压信息法、全速势位流方法、基于 RANS 的壁压信息法三种方法开展飞翼标模的洞壁干扰修正研究,并与试验结果进行了比较。结果表明,飞翼标模洞壁干扰呈现与翼身尾常规布局不同的规律,三种方法在飞翼标模洞壁干扰修正中有各自的适用性。     

圆管内不可压粘流的分布式并行算法

用并行算法求解了圆柱通道内不可压层流流动。采用有重叠的条形区域分裂法时计算域进行分区，分别求解动量方程和压力校正方程，对压力校正方程采用块修正法来加速线代数方程的收敛速度。计算结果表明有较好的并行加速比和并行效率。     

钝体前激波的数值模拟

本文通过直接求解N-S方程,研究激波问题,对是否可以和如何用直接求解N-S方程来处理超跨声速流进行了探索,在对粘性项的处理建立可压缩流体的压力校正方程和计算方法等方面均提出并实施了新的方法,最终得到了基本合理的数值计算结果。     

人－椅系统的动压计算问题

通过大量的分析、计算,提出了作用在人－椅系统上的动压修正系数F的半经验公式。给出选取修正因子ε、f（ε）的方法,仅需４个不同马赫数的阻力系数实验值。原则上,若实验结果精度较高,对阻力系数实验值的选取并无限制。但实验数据精度不高,则上述计算的阻力系数的均方根误差就会大些。此方法不仅适用于计算像人－椅系统这样的大钝体,而且也适用于计算其它大钝体在任意马赫数范围的阻力系数,其均方根误差小于01 05。     

复杂区域内粘性流动的数值计算

本文应用原始变量法(SIMPLEC)计算了复杂区域内的不可压和可压缩流动问题。计算时,仍使用交错网格,但在求解动量方程时采用协变分量做为原始变量,有效地消除了在任意非正交坐标系下压力和笛卡尔速度分量不耦合而产生的压力波动。在压力修正方程中考虑了密度的影响并构造出一种密度的迎风格式,保证了计算的稳定与收敛,使之预示马赫数变化范围较大的流动成为可能。对由NACA65(12A10)10构成的平面叶栅、由NGTE10C4/30C50构成的平面叶栅的不可压缩流动和平面超声速缩放喷管进行了计算。结果表明,本方法有较好的收敛性和准确性,能够适用于从亚音速到超音速的流动场合,有一定实用价值。      

飞机空调系统中饱和水蒸气压的计算

提出了方便实用的饱和水蒸气压方程, 与空调系统计算中原来使用的方程相比, 新方程在计算工作量相当的情况下, 显着提高了计算准确度, 能在飞机空调系统可能的饱和温度范围内, 获得满意的计算结果。      

预估-校正LU-SGS的隐式算法

为了进一步提高Euler方程求解效率,在传统LU-SGS格式基础上提出了预估-校正LU-SGS。结合高分辨率迎风通量格式,发展了一套高效的Euler方程求解器。以NACA0012翼型和RAE2822翼型跨音速无粘流动作为算例,对比分析了预估-校正LU-SGS、传统LU-SGS和GMRES格式的计算效率。结果表明:预估-校正LU-SGS算法能显著提高求解效率,其效率接近GMRES算法,比传统LU-SGS方法提高了近三倍。     

考虑低温增压真实气体效应的运输机气动特性数值模拟研究

为了给低温风洞试验数据修正提供参考,本文利用数值模拟手段研究了低温真实气体效应相比较于完全气体对飞行器气动特性的影响,以及该影响与雷诺数影响相比所占比例的大小等问题。文章应用 Aungier-Redlich-Kwongz 方程,发展了适用于模拟氮气低温高压真实气体效应的 RANS 求解软件。与 NIST 数据的对比表明,该状态方程在5倍大气压下,比热比等参数误差在0.3%以下。同时,标模测试结果表明本文软件计算精度与国外软件相当。应用本文方法研究低温增压风洞中氮气真实气体效应：以典型运输机构型 DLRF6为模型,分别计算了高速、低速状态下各种不同温度和压力工况下的流动。计算结果表明,在低温增压情况下,真实气体效应引起的气动力差异很小,升力、阻力、力矩最大相对误差均在0.3%以下,与雷诺数效应引起的偏差相比可以忽略不计。因此,以氮气为介质的低温风洞试验研究可以采用完全气体假设。