三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

中心锥式进气道跨音速内外流场的混合差分计算

一、引言 文献[1]曾用跨音速定常势流的混合差分法成功地计算了皮托管式轴对称进气道内外跨音速流场,本文进一步用上述方法计算了中心锥式进气道的跨音速流场。中心锥式进气道的通道形状大多为收缩-扩张型,内通道壁面扰动较强。同时,由于中心体的存在,使传递进气道出口反压扰动的横向松弛线被隔断,因此,在皮托管式进气道计算中行之     

离子推力器栅极系统电子反流阈值的数值分析

阻止束流等离子体中电子反流到加速栅上游区域是离子推力器加速栅负电压的主要作用之一,能够阻止电子反流的加速栅电压最小值称为电子反流阈值。加速栅电压的选择直接影响到离子推力器的工作性能和运行寿命,电子反流阈值电压是确定加速栅电压的重要参考参数。基于PIC方法计算了20cm氙离子推力器加速栅电子反流阈值,并分析了加速栅孔径、栅间距、单孔引出束流电流大小对加速栅电子反流阈值电压的影响,计算结果与试验测量值符合较好。该数值模型为加速栅参数的选择和降低电子反流失效风险方法提供了参考,为下一步电子反流现象对加速栅寿命的预测分析奠定了基础。     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

叶型反设计计算研究

以某型号发动机压气机叶型生产坐标为设计输入,寻找一种切实可行的叶型反设计方法,以便对该压气机进行反设计计算研究,评估其真实性能,从而为改进设计提供依据.具体反设计过程为:根据等高面叶型生产坐标、设计经验和分析判断给出压气机各级气动参数分布及其子午面布局,得出其S2流面计算结果,然后冉根据S2流面计算结果将生产坐标反设计到流线面上,通过流线上的叶型反设计得出一维、二维及三维计算分析程序所需的计算输入并进行计算分析,对比计算和试验结果,判断反设计方法的有效性.     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

大速度前飞时旋翼性能计算

本文基于压缩性修正,非定常修正,径向流修正,动态失速修正和反流区修正,给出了计算大速度前飞旋翼性能的方法。计算结果表明此方法精度较高,可用于型号设计。     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

动态网格快速生成及组合体跨音速非定常气动力计算

采用代数方法、局部变形技术及快速动态弹性变形技术生成三维翼/体组合的贴体动态网格。用时间精确近似因式分解差分方法求解守恒型非定常全位势方程。考虑附面层影响，通过粘/位迭代得到跨音速翼/体的非定常气动力。所得结果与实验数据相比相当吻合。     

计算三维边界层分离流动的一种双向推进有限元格式

 本文提出一种适用于兰维边界层正、逆向问题的有限元格式。格式采用了对壁面法向进行Galerkin加权,沿壁面进行子域加权并作双向推进的弱解形式。它不仅可节约内存,同时还可通过与等参元素的结合应用给出CFL条件的显式表述,保证了问题的快速求解。进行逆模态计算时,以位移厚度作约束条件,通过Lagrange乘子直接引入边界层方程的有限元类比中进行整体求解。算例表明该格式模拟三维边界层流动是经济有效的。     

边界层分离流边值问题的算法研究

本文给出了边界层方法边值问题的一种快速求解方法,该方法可用于二维与三维问题,它是在正方法边界层算法的基础之上增加了一个矩阵算法,该矩阵很容易逆,总的ＣＰＵ时间仅增加约２０％,使用该算法对反方法问题求解和进行数值稳定性分析都将十分方便。     

三维跨声速流粘性—无粘相互作用问题的分析和计算

本文对跨声速粘性-无粘相互作用问题进行了分析,论证了如何选择湍流边界层反模式的输入量。在FL027多重网格基础上作了局部加密。计算例子表明,本文建议的半逆方法适用于弱激波情形,可用于实际设计之中。     

Two 3-D multi-cascade co-design inverse methods based on streamline curvature approach and Clebsch formulation for single-stage centrifugal compressors

Two strategies extended the single-cascade methods from a compressible three-dimensional inverse method for radial and mixed flow turbomachines to two three-dimensional multi-cascade co-design methods for single-stage centrifugal compressors.These two three-dimensional methods and a typical quasi-three-dimensional streamline curvature through-flow inverse method were employed to design the same subsonic high-speed single-stage centrifugal compressors.The compressor performances were simulated by a commercial Reynolds averaged Navier-Stokes (RANS) equations solver.The studies show that two three-dimensional co-design methods are reasonable and feasible.It was found that:firstly the blade camber angle designed by the three-dimensional methods was larger than that designed by the quasi-three-dimensional method;and secondly with regard to two three-dimensional methods with different boundary conditions,the co-design result differences between the diffusers were small,but those between the deswirlers were relatively large.      

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。