失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

旋翼/机身气动干扰的数值模拟

通过求解粘性非定常可压缩N-S方程,发展了旋翼/机身气动干扰的数值模拟方法,主要技术措施包括:(1)通过ADT数据结构,建立了高效的多块对接网格动态重叠技术;(2)在可压缩方程求解中引入低速预处理方法,克服低速流场求解时的刚性问题;(3)加入了多重网格方法,提高收敛效率;(4)湍流模型方程求解需要反复计算壁面最近距离,采用Wigton的优化算法很大程度上提高了计算速度。利用建立的计算方法对佐治亚大学(GIT)旋翼/机身干扰模型进行了数值模拟,分别比较了机身表面的平均压力分布、旋翼的下洗流场以及非定常压力脉动,对干扰现象及其流动机理进行了分析。     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

用预处理方法数值模拟街区峡谷内的污染物扩散

通过耦合求解二维可压缩Navier-Stokes(NS)方程和k-ε湍流模型方程及平流扩散方程，数值模拟了街区谷（Canyon)内的低速湍流流场及大气污染物在其中的扩散。为了加快计算的收敛，对控制方程进行了预处理。预处理后的NS方程的求解采用LU－SGS隐式方法和JST中心差分格式，数据模拟结果表明预处理方法求解低速流场可加快敛，同时就街峡欲内的污染物扩散状况得出了有意义的结论。     

跨洞庭湖大桥主梁节段定常气动力研究

采用风洞实验与数值模拟的相结合的手段,对跨洞庭湖大桥主梁节段的定常气动力特性进行了研究;流场数值模拟是采用ｋ－ε湍流模型封闭的雷诺时均Ｎ－Ｓ方程,求解方法是建立在交错网格了的压力－速度场迭代求解的ＳＩＭＰＬＥ算法;实验是在国防科技大学ＫＤ－０３低速风洞进行,研究表明实验与计算结果较为一致。     

多段翼型低速流场与气动特性数值计算研究

在结构化网格求解RANS方程的基础上,将嵌套网格技术应用在低速绕流气动特性计算中,并以NLR-7301两段翼型和某运输机5段翼型在起飞构型下对低速绕流流场进行数值模拟,将气动特性及流场特性的数值模拟数据与对接网格计算以及风洞实验数据进行了比较.结果表明,该气动特性计算方法有着较高的精度和可靠性.     

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。     

空化流动隐式求解方法研究

空化流动问题本质上是可压缩流动,应用可压缩方法开展数值模拟研究更符合物理实际。气体和液体压缩性的显著差别使得低速空化流动数值模拟的刚性问题非常突出,通过引入预处理技术解决低速问题中由于特征值量值不一致导致的收敛刚性问题,提高收敛速度。同时,以预处理后的特征值构造Roe格式耗散项,提高低速流动计算精度。鉴于自然空化流动中气液组分转换现象和物质输运现象并存,且气体和液体的密度在常温状态下存在显著差别,预处理后的源项雅克比矩阵的特征值与无粘通量雅克比矩阵的特征值存在量级的差异,这会使得求解过程不稳定或收敛速度极慢,即出现"源项刚性"问题。为此,本文系统推导了预处理框架下的气液两相流隐式求解方法,采用点隐式方法处理源项,通过直接求逆的方式增强算法的稳定性。研究中分别考察了三种不同的算子分裂方式在LU-SGS(Lower-Upper Symmetric-Gauss-Seidel)隐式迭代中的模拟效果,并提出适用于DDADI(Diagonally Dominant Alternating Direction Implicit)的算子分裂方式,在NACA0015水翼空化算例模拟中综合比较了上述四种隐式迭代方式对低速空化问题收敛性的影响。最后,应用本方法对三维尖锥自然空化算例进行了考核模拟,捕捉到自然空泡流场中的主要特征,计算结果与试验测量结果吻合。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

A new accelerating technique for low speed flow: pseudo high speed method

This paper proposes a new accelerating technique for simulating low speed flows, termed as pSeudo High Speed method (SHS), which uses governing equations and numerical methods of compressible flows. SHS method has advantages of simple formula, easy manipulation, and only need to modify flux of Euler equations. It can directly employ the existing well-developed schemes of hyperbolic conservation laws. To verify the technique, several numerical experiments are performed, such as: flow past airfoils and flow past a cylinder. Analysis of SHS method and comparisons with some precondition methods are made numerically. All tests show that SHS method can greatly improve the efficiency of compressible method simulating low speed flow fields, which exhibits in accelerating the convergence rate and increasing the accuracy of the numerical results.     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.