二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

隐式多重网格法求解叶轮机械三维跨声速湍流流场

本文提出了一种高效率求解叶轮机械三维跨声速湍流流场的隐工多重网格算法。它以矢通量分裂格式为基础，采用ＦＡＳ型多重网格加速技术。     

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

分层Reynolds应力模式

本文根据切变湍流结构的特征提出一种分层Reynolds应力模式。在切变湍流中存在湍能生成项达到狭峰的极大值层次(本文称它为活跃层),根据这一特征,在活跃层中保留微分Reynolds应力模式,但它可以简化为常微分方程。在两层活跃层之间以及活跃层和壁面之间的Reynolds应力分布,用满足衔接相容条件的多项式近似。时均湍流场将由以上得到的Reynolds应力代入Reynolds方程求出。用本方法计算了两个比较复杂的湍流例子,结果表明,本方法较常用的k-ε模式有明显的优越性,它的计算时间并不比k-ε模式多,而计算结果与实测结果的吻合则较k-ε模式好得多。另一方面本方法较之完全的微分Reynolds应力模式可大大节省计算时间。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

模型燃气轮机燃烧室中三维气相反应流的数值计算

 本文给出一种模型燃气轮机燃烧室中三维气相反应流流场的数值计算结果。计算中采用了著名的SIMPLE算法,改进的κ-ε双方程湍流模型及湍流扩散燃烧模型。研究了燃烧过程对流动图形的影响。计算结果经过实验验证,两者符合较好。     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。     

二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

带尾支杆的超声速弹丸湍流底压计算

尽管在超声速湍流情况下，尾支杆对弹丸底部压力的影响没有层流情况下的影响大。但是实验研究表明，带尾支杆的超声速弹丸湍流底部压力随湍流边界层相对厚度增加而增加，随天平尾支杆的相对直径增加而减小。本文介绍一种考虑圆柱体上边界层状态以及尾支杆对绕弹丸超声速湍流底压影响的计算方法，计算结果和实验数据吻合较好，还根据计算和实验结果分析了尾支杆对底压系数的影响。     

任意边界扩压器内紊流流动的时间推进数值计算

研究了隐式ＡＦ有限差分理论。指出薄层抛物化控制方程中的右端粘性项是主元素变量矢的一个简单线性函数。给出了任意边界扩压器内紊流流动的数值计算结果，并和实验数据进行了对照。两者吻合得相当好。研究结果对ＡＦ有限差分理论以及推广至三维复杂管道内的非定常流动计算具有重要的参考意义。     

轴对称超临界粘性喷流的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式中的矢通量分裂技术，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了收敛喷管的内流场和喷射流场，通过对趋临界工况下喷管跨音速流场的数值计算，表明本文的结果与试验结果是吻合的。     

三维控制翼周围粘性分离流动的数值模拟

本文首先分析了中心型隐式格式和迎风型隐式格式的优缺点,在此基础上发展了一种自由参数可确定的反扩散隐式格式。用此格式计算了绕三维控制翼的层流和湍流分离流动。出发方程是三维简化N-S方程,用Baldwin-Lomax湍流模型。结果表明,这种隐式格式效率高,能满意地刻划流场。     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

计算全马赫数下粘性流动的隐式算法

把预处理方法引入到求解三维粘性流动的Navier-Stokes方程隐式格式算法中,并应用这种方法对平板算例在层流和湍流两种不同流态下及四个不同Ma数下进行了数值模拟,计算结果与理论值及相关公式吻合,同时对比分析了有无预处理情况下的收敛性,证明了本文所采用的预处理隐式算法对于任意马赫数都是适合的,而且可以提高残差精度量级,加快收敛。     

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.