RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

时间相关法计算喷管跨音速流动

本文用时间相关法完成了定常、无粘跨音速喷管流场的计算,内部点参数用Mac Cormack差分格式计算、壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。 定常、无粘、跨音速喷管流动的控制方程是混合型的(即在亚音速区域是椭圆型,而在超音速区是双曲型的),给数值求解带来很大的困难。为了克服这一困难,广泛采用时间相关法,即认为定常流动方程的解是相应的非定常流方程的解在时间趋于无穷时的渐近解。因为非定常流的控制方程不论在亚音速区,还是在超音速区域,都是双曲型的,可以用统一的方法来求解,而且易于求解;另一方面,它可适用于形状比较复杂的喷管。 本文首先阐述了时间相关法的计算方法,然后列出算例喷管的计算结果。通过与实验数据的比较,证明计算是符合实际的。     

二维湍流边界层雷诺应力方程的数值解

本文给出一种解雷诺应力方程的一种新的显式方法。本方法对坐标和变量作了规格变化,并对方程作了特征变换,与控制方程不作特征变化的方法相比,本方法十分有效。用本方法作了四种算例计算,其计算结果与实验结果十分一致。     

用椭圆型、抛物型两种方法计算平面叶栅湍流尾迹

本文用连续方程、稳定不可压缩椭圆型和抛物型时均N-S方程及K-ε双方程湍流模型计算了平面叶栅湍流尾迹。计算表明,上述两种方法都能得到较好的数值结果,但抛物方法可节省大量计算时间和内存空间。     

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

绕10°、20°锥柱体跨音速流动的混合有限差分计算

对于10°、20°锥柱体,用跨肯速小扰动位势流的混合有限差分格式,用跨音速完全位势流的旋转差分格式与简易差分格式,做了数值计算。完全位势方程旋转差分格式与简易差分格式得到的数值结果是一致的,但前者的计算马赫数范围为宽。20°锥柱体完全位势方程的数值结果与实验结果,基本一致,只在顶部与肩部有些差异。小扰动位势方程的数值计算对于10°锥桂体是适用的,对于20°锥柱体足不适用的。     

激波与壁面粉尘相互作用的数值研究

本文采用厚边界层理论建立了考虑两相耦合效应的气固两相湍流边界层方程用以研究激波与壁面粉尘的相互作用。方程中考虑了压力梯度，Ｓａｆｆｍａｎ升力和湍流扩散等因素对颗粒飞行轨迹的影响，并用Ｍｉｒｅｌｓ方程作为边界条件使方程组封闭可解。数值计算时，经适当变换的气相方程组用ＢＯＸ格式，拉格朗日坐标中的颗粒相方程用Ｒｕｎｇｅ－Ｋｕｔｔａ法，而源项则用颗粒源法和求质量加权的平均值获得。计算结果与某些实验结果较为     

跨音速压气机在非设计工况下动静叶相互作用的数值模拟

本文计算了某跨音速压气机的一级在非设计工况下动静叶干扰产生的非定常流动。计算以准三维的Navier-Stokes方程为模型,湍流模式采用Baldwin-Lomax代数模式。主控方程用高分辨率的MUSCL型Roe格式离散,离散方程用高效的LU-SSOR方法求解。本文提出了一种完全守恒的动、静叶交界面处理方法。周期性边界处理采用Rai的方案。根据计算结果,讨论了非设计工况下动静叶干扰的机理。      

钝体前激波的数值模拟

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

高精度差分格式在机翼分离流计算中的应用

用一个模型方程分析了边界层方程的数值稳定性 ,指出稳定性问题会随计算雷诺数减小而变得严重 ,因而对较低的Re数 ,计算分离边界层流动采用高阶精度差分格式十分必要。本文同时给出用四阶精度差分格式求解在流动分离情况下边界层方程的技术方法 ;对一个大展弦比后掠翼在攻角等于 1 4°时计算了机翼上表面的分离线 ,并与文献上公布的用别的测算方法测出结果做了比较 ,两者总体上符合良好     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

平面叶栅跨声速非等熵流的时间推进解

本文根据Viviand等人提出的求解平面流动问题的拟不定常时间推进法的思路,推导了完整的求解平面叶栅无粘非等熵拟不定常流动的方程系统,应用经过改进了的有限控制面积时间推进格式进行了数值计算。计算结果与精确解或试验结果吻合较好,且能有效地节省计算机时。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。