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

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

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

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

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

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

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

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

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

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

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

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

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

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

圆筒燃烧室紊流流场数值研究

计算方法 本方法不需解压力修正方程,而是通过速度散度余量修正迭代法使连续方程直接得到满足。紊流模型采用k-ε模型,壁面处理采用无滑流固壁条件,避免复杂的壁面函数计算。对燃烧室内火焰稳定器,无论是圆盘或锥形都采用简单的取除格子法,从而大大简化计算程序。同样算例本方法所用计算语句只是SIMPLE算法的1/4～1/3。      

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

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

非均匀来流稳态畸变流场的数值计算

本把TEACH计算机程序推广到三维，用定常，绝热，不可压的时均三维N－S方程和连续方程，并运用紊流动能和紊流耗散率的微分输运方程封闭雷诺方程，模拟计算了非均匀来流稳定畸变流场。数值计算中除了对通用变量进行欠松驰，还对修正压力通过惯量进行了欠松驰，惯量取为压力修正方程中的系数|b|，保证了数值计算的稳定性和收敛性。在进口马赫数小于0．35时，流场畸变指数的计算值和试验值符合很好。     

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

在 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方法提高了近三倍。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。