基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

基于非结构重叠网格的二维N-S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

基于网格变形的伴随方法在翼型优化中的应用

进行了基于网格变形的伴随方法在翼型气动优化设计中应用的研究。随着集成了伴随方法的流场求解器的普及和效率的不断提高,采用梯度法进行翼型气动优化设计的主要计算花费逐渐由计算梯度对流场变量的敏感度（或流场敏感度）转变为计算梯度对网格变形的敏感度（或网格敏感度）,后者的求解通常采用有限差分方法,计算花费较高。在传统伴随方法（或流场伴随方法）的基础上,引入了基于网格变形的伴随方法（或网格伴随方法）,采用网格伴随方法计算梯度,可以大幅度减少梯度计算花费,提高翼型气动优化设计的效率。     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

非结构化多重网格流场数值计算

本文对二维非结构网格无粘流场计算的多重网格方法作了探讨。研究了非结构化多重网格的生成方法,插值算子的构造以及流场计算．以绕流ＮＡＣＡ００１２翼型无粘流场为算例,研究了多重网格法在非结构网格上的实现方法及其加速效果表明本文的方法有显著的加速收敛效果。     

基于MADS算法与嵌套网格技术的多段翼型优化设计方法研究

基于MADS算法、嵌套网格技术以及刚性动网格技术,以求解RANS方程为气动特性分析方法,对某大型客机带增升装置的多段翼型着陆构型下的缝道、重叠量和偏转角等参数进行优化设计,结果显示,该方法有较高的优化效率及计算精度。     

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

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

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

基于重叠网格技术和多重网格算法的悬停旋翼粘性绕流数值模拟

发展了一种基于重叠网格的多重网格算法,并应用于悬停旋翼粘性绕流的Navier-Stokes方程数值模拟.多重网格算法在重叠网格上实施主要面临如下两个问题:粗网格间几何关系的建立,即洞单元、边界单元的确定;洞边界存在的情况下如何进行粗细网格间流场信息的传递.本文针对以上问题给出了相应的解决方案.采用两层网格的V循环,应用有限体积法对一跨声速悬停旋翼粘性流场进行了数值计算.计算结果表明该多重网格方法可以显著加快悬停旋翼粘性流场数值解收敛速度,加速比约为3.     

High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows

In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The high-order variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4est based on the forest of octrees is adopted. Several two- and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.     

适用于混合网格的约束最小二乘重构方法

为强化边界条件对计算流体力学(CFD)计算中流场重构过程的影响,基于三维任意多面体混合网格,对一种约束最小二乘法进行了研究。通过对不同边界类型的透明化处理,统一了约束方程组的形式,从而简化了该方法的应用。为了充分吸收消去法和加权法各自求解约束方程组的优点,提出了混合采用两种方法的思路。针对层流平板研究了加权法中加权系数如何取值问题,数值实验表明加权系数取到5基本已经足够大。最后通过亚声速层流平板和跨声速湍流ONERA-M6算例对比了混合法、加权法与原始最小二乘法,计算结果表明:约束最小二乘法相比原始最小二乘法的流场重构更加准确,特别是对边界质量较差网格的计算优势更为突出;混合消去加权法相比单独采用加权法的计算结果也有所改善。     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

拓扑结构改变流场数值方法研究及应用

物体从相互接触的状态分离时会导致流场拓扑结构的改变,针对拓扑结构改变过程提出了一种非结构重叠网格算法.该方法主要包括"洞边界"自动判定、重叠区形成和流场信息交换三个方面的关键问题.针对三维问题提出了判定洞边界的"曲面判定方法";对于重叠区待插值的网格节点建立了"动态义树"贡献单元搜索方法.洞边界上待插值点的流场参数通过插值方法从其贡献单元上得到.弹丸从炮口发射的过程中,由于弹丸与炮管的分离导致了流场拓扑结构变化,而当炮口带有制退器时将使流场的边界条件更加复杂.利用非结构重叠网格方法对弹丸从带有炮口制退器的火炮中发射过程进行了数值模拟,分析了初始冲击波形成和发展过程,详细描述了运动弹丸对膛口流场的影响.本文提出的算法将能够有效地应用到航空、航天以及兵器技术领域的研究中,特别是物体分离过程导致的流场拓扑结构改变问题.     

稀疏非结构网格上的亚声速流高精度数值模拟

采用高阶间断有限元法在非结构网格上数值求解二维亚声速Euler方程。数值结果表明,尽管采用的非结构网格非常稀疏,但通过采用真实曲线物面边界和高阶的基函数仍可得到高精度的数值解。另外,对于超低速情况,方程无需经过任何特殊处理就可以得到收敛的数值解。由于采用牛顿一般最小余量法(Newton-GMRES)时通常需要较好的初始值,本文设计了一种阶谱循环过程来提高数值求解的鲁棒性。     

飞机最短时间爬升轨迹的格点搜寻算法

以飞机在垂直平面的运动方程为基础,采用线性时间负荷控制的飞行模式,提出了求解飞机最短时间爬升轨迹的格点搜寻算法.该方法视飞机为一质点,以此为基础计算出每个相邻格点的时间间隔,最后用h×Ma格点数搜寻法求出飞机最短时间爬升的最优轨迹.数字仿真表明,在规定的爬升范围内用格点搜寻算法可以实现快速的时间响应,得到的最优轨迹能够准确地跟踪设定的性能指标.     

结构/直角切割网格下时域有限体积法在计算电磁中的应用研究

 采用有限体积法直接求解时域麦克斯韦方程组，研究了二维圆柱、翼型和三维球体、立方体、双球体等典型完全导电目标的电磁散射问题。空间离散使用近似黎曼解构建网格通量，并针对二维结构网格、三维自适应直角切割网格分别给出重构方法，时间方向采用二步龙格库塔法。计算与理论、实验或文献结果符合良好。     

用有限元显式算法解具有烧蚀动边界的传热问题

本文借助于运动坐标系用显式有限元算法求解烧蚀动边界的传热问题。通过坐标变换,在变换平面内确定固定的有限元结点,物理平面内对应的结点随边界的烧蚀而运动。在线性插值函数模式下温度的时间导数在三角形单元中处处相等的假设下导出有限元方程的显式计算格式。对横观各向同性材料的端头帽在再入飞行中的烧蚀和瞬时温度场的计算结果与用有限差分法计算的结果相互吻合。     

冲压发动机进气道压力振荡过程的数值研究

针对来流马赫数3.5,0°迎角飞行条件下冲压发动机的简化模型,数值模拟了进气道整流罩开启过程的非定常流动。计算采用AUSMPW有限体积格式,基于结构与非结构的混合网格技术,求解了非定常N-S方程。计算结果表明,在进气道整流罩打开以后,超声速气流会在燃烧室通道内形成振幅衰减较慢的低频剧烈振荡。在同样计算条件下采用Euler方程进行对比性研究,数值实验发现对这类大雷诺数、强非定常性的问题,Euler方程和N-S方程得到的结果相差不大。     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。