共查询到19条相似文献,搜索用时 312 毫秒
1.
2.
准确模拟高速流场中的热化学非平衡问题一直都是各种流场数值计算的难点。DSMC作为一种有效的流场模拟方法,从统计学的角度而言,作为统计母体的模拟分子数是影响计算准确性的关键。Bird等[Bird,1990]针对不涉及热化学非平衡的流场提出了一个统一的模拟分子数选取规则。本文在热浴条件下使用DSMC模拟了空气中三种常见氮氧气体的离解反应,对涉及热化学非平衡的流场中模拟分子数对DSMC精度的影响进行了分析。结果表明:与Bird针对普通流场所提出的统一模拟分子数选取规则不同,在不同的流场物理条件(本文主要为温度)下,得到正确结果的模拟分子数阈值与温度成正相关(最小离解度与最大离解度模拟分子数选择下限相差10倍),超出阈值后增加模拟分子数对流场精度并无明显提高,这一结果与Parker的分子能量传递理论相符。从保证计算精度和平衡计算效率的角度,本文的结论对DSMC非平衡流场模拟具有参考价值。 相似文献
3.
通过对基于颗粒输运随机统计模拟的DSMC方法模拟准则的分析研究,定义了该方法中统计量的噪声、相对统计涨落及置信度。引入概率论中的切比雪夫不等式,利用正态分布平衡流统计理论,推导出DSMC方法模拟流场宏观速度和温度的统计噪声理论公式,并以过渡流区平面Couette剪切流和激波结构内流动为算例进行数值模拟验证与理论分析。结果表明,在一定的置信度下,选择适当的模拟参考值,DSMC统计量的噪声除与样本容量(模拟分子数及统计抽样次数)的平方根成反比例外,速度的相对统计涨落还与流场马赫数成反比,而温度的相对统计涨落与流场宏观参数无关。这为控制DSMC模拟结果的相对统计涨落提供了理论指导和计算依据。 相似文献
4.
本文用NND格式求解了化学非平衡流场,利用点隐式处理组分连续方程的刚性,采用全流场激波捕捉办法和特征理论,成功地进行了物面边界条件的处理,模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型,建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。 相似文献
5.
6.
为了研究混合气体分子输运系数的不同近似计算方案对边界层及射流剪切层流场计算的影响,以含双方程k-ωSST湍流模式的质量平均Navier-Stokes方程为控制方程,分别采用两种计算方案:公式计算混合气体导热系数、扩散系数的直接计算方案和由普朗特数和施密特数等关联混合气体导热系数、扩散系数的关联计算方案,对Burrow和Kurkov实验流场进行了数值模拟。结果表明:与采用直接计算方案相比,采用关联计算方案,流场数值模拟在获得相当的计算结果的同时,大大降低了计算量;抛开与控制方程的相容性问题,混合气体扩散系数各直接计算方案对流场参数计算影响不明显;在湍流计算时,边界层/剪切层区域的混合气体湍流输运系数大于分子输运系数,由此弱化了混合气体分子输运系数不同的近似计算方案对流场数值模拟的影响;相对输运系数近似计算方案,湍流可压缩修正模型对组分剖面的影响更为明显。计算表明,在边界层近壁流场、剪切层流场以及层流流场的精细化数值模拟中,混合气体分子输运系数的计算方案有重要的影响。 相似文献
7.
8.
直接模拟蒙特卡罗羽流模拟的两相作用模型 总被引:1,自引:1,他引:1
为了数值模拟包含气体分子和凝相颗粒的两相羽流场,使用自主开发的羽流计算软件(PWS)对两相作用模型进行研究。计算均匀气流对颗粒的加速过程,结果与理论解一致,证明气体分子对颗粒的单向作用模型合理。计算气体分子在颗粒表面的反射时,将全局坐标系转换到颗粒坐标系确定其法线,再转换到法线坐标系进行镜面反射或漫反射。将这种间接法与使用反射向量的分布函数的直接法进行了对比。对同一个气体分子在颗粒表面的反射进行1×106次数值试验,两种方法所得镜面反射的速度向量均符合空间均匀分布,漫反射的速度向量分布相互一致。使用两种方法模拟150N液体火箭发动机的两相羽流场,所得气相羽流场几乎相同,间接法计算效率较高。液滴颗粒对气流场的影响主要集中在轴线附近,与无液滴的气流场相比,有液滴时会使气流温度升高、压强增大、马赫数减小。 相似文献
9.
研究了简化模型超音速燃烧室流场和性能的数值模拟方法,该燃烧室在支板后缘设置一缝隙式喷嘴,平行于超音速空气流喷射氢气。用椭圆型偏微分方程数学模型,MacCormack差分格式成功地计算了两种不同后绿尺寸的流场及其性能,为了对照比较,还计算了另外两种流场。计算结果表明,支板后的流场存在回流区,因喷氢的压力高于超音速空气流的压力和壁面附面层的影响,燃烧室内将出现斜冲波和膨胀波,压力沿横向变化明显,由此证明流场与“边界层流动”有性质上区别。提出了全流场按性质不同分段进行数值模拟的方法,利用计算得到的流场节点状态参数(如温度、压力等),积分获得燃烧室的性能参数,和文献的实验数据对比,计算结果合理、可靠。 相似文献
10.
11.
基于高阶WENO格式的旋翼非定常涡流场数值模拟 总被引:1,自引:0,他引:1
为提高直升机旋翼黏性涡流场计算流体力学(CFD)的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡(WENO)格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T(Caradonna-Tung)旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明:在计算时间增长10%~20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。 相似文献
12.
密切曲面内锥乘波前体进气道设计和试验研究 总被引:3,自引:0,他引:3
介绍了密切曲面内锥乘波前体进气道(Osculating Inward turning Cone Waverider Inlet,OICWI)的一体化设计方法,对该型乘波前体进气道的性能进行了数值分析,针对该型一体化乘波前体进气道完成了风洞试验研究。理论设计结果和设计状态无粘模拟结果一致,设计状态下的计算结果表明,前体进气道具有较高的总压恢复、较好出口流场均匀度及较高的流量捕获率。试验研究结果表明,改型一体化前体进气道在马赫数5~7条件下顺利启动,流场波系及压力分布同数值分析结果吻合。 相似文献
13.
进口型线水平投影可控的变截面内收缩进气道设计 总被引:1,自引:0,他引:1
发展了一种进口型线水平投影可控的三维变截面内收缩进气道设计方法。基于反正切马赫数分布基准流场,在指定进口型线水平投影为椭圆和出口为圆的条件下,结合流线追踪和截面渐变技术设计了光滑过渡的内收缩进气道。在设计点(M_(ai)=5.4)和接力点(M_(ai)=4.0)对其进行数值仿真,计算结果表明,设计点时进气道的主要流动特性与基准流场基本一致,无黏时可以捕获98%的自由来流,喉道性能与基准流场基本相等。相对椭圆进口进气道,截面渐变的椭圆转圆进气道流场结构相似且性能下降较小,有黏条件下设计点和接力点时喉道总压恢复系数分别降低了2.9%和1.2%。此外,该进气道表现出良好的总体性能,接力点的流量系数达0.82。 相似文献
14.
15.
16.
17.
18.
采用隐式近似因子分解法,通过求解三维、非定常、可压缩完全层流N-S方程,对绕大后掠三角翼旋涡破裂的流场进行了数值模拟。结果表明:当三角翼背风区流场发生旋涡破裂时,各流动参变量是非定常的,同时旋涡破裂点的位置沿轴向在一定范围内前后摆动。 相似文献
19.
四段修型弥散反射激波中心体基准流场研究 总被引:3,自引:0,他引:3
针对高超声速内收缩进气道宽马赫数范围工作的要求,设计了一种四段修型弥散反射激波中心体基准流场,可明显提高基准流场来流马赫数高于设计点来流马赫数(6.0)时的压缩效率,巡航点(来流马赫数为7.0)出口总压恢复系数较下凹圆弧中心体基准流场提高了2.3%.此外,基于两种基准流场分别设计了圆形进口内收缩进气道并在来流马赫数为5.0~8.0范围内进行数值计算,结果表明:来流马赫数高于设计点来流马赫数时,四段修型进气道的压缩效率更高.有黏条件下,来流马赫数为8.0时二者的增压比近似相等,四段修型进气道喉道截面出口总压恢复系数相对提高了1.1%. 相似文献