分区对接网格在跨流域气体运动论统一算法的应用研究北大核心CSCD

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

外挂物干扰流场特性数值仿真研究

飞行器研制中对于带外挂物的复杂外形气动特性分析研究一直是个难点,对外挂物与挂架不同缝隙的复杂机翼-挂架-外挂物组合外形采用结构重叠网格数值方法,快捷高质地完成了网格生成,通过求解带有k-e湍流模型的Navier-Stokes方程组,得到多个计算工况下的绕流外流场,分析了外挂物在不同工况下的气动特性,结果表明,CFD数值方法可方便快捷地对复杂干扰流场进行数值计算,可为飞行器相关设计提供技术支持。     

带副翼三维机翼绕流的N-S方程解

提出了一种求解带副翼偏转三维机翼绕流的N-S方程计算方法。采用对接分区网格与分区求解算法的结合,有效地求解绕此外形的复杂流动。提出了一种满足通量守恒的内边界耦合条件。数值方法中选用vanLeer分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD性质,采用中心差分格式来离散粘性通量项。数值算例表明该方法是求解带操纵面偏转的机翼绕流的有效方法。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

探月返回器稀薄气体热化学非平衡特性数值模拟

通过发展直角/非结构混合网格下,适于复杂返回器含内能激发的五组元热化学非平衡稀薄气体直接模拟蒙特卡洛(DSMC)方法,数值模拟了钝体返回器稀薄流区7.5 km/s和10.6km/s再入速度下的非平衡流场特性和对物面气动力/热分布的影响。计算结果表明在月球返回速度下稀薄过渡流区存在着强烈的热化学非平衡,近连续滑移区仅在激波层区域存在着较大的平动、转动和振动非平衡度,而在高稀薄流区,热力学非平衡遍布绕流物体四周,包括整个压缩区和尾迹区。高温真实气体效应对表面热环境影响显著,对表面压力和摩擦力的影响相对较弱。计算结果显示探月返回器再入绕流过程100 km飞行高度仍需要考虑热化学非平衡的影响。     

全机复杂流场网格生成及N-S方程数值模拟

全机复杂流场网格生成对计算结果起着举足轻重的作用。本文首先通过合理的拓扑结构,分块对接的方法对复杂流场进行了网格划分,然后采用N-S方程有限体积法及分区对接法对翼—身—尾全机流场进行数值模拟。计算结果验证了网格生成的有效性和计算结果的合理性。     

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

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

分区对接网格算法的应用研究

本文应用分区对接网格技术，由于采用了对接面上网格点完全对应的方法，能自动满足通量守恒和计算的稳定性。采用有限体积算法下的原始变量型NND格式，求解NS方程，数值模拟了复杂外形的高超声速层流流场，计算结果与实验测量符合较好。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究，将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元，采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元，更好地表述了真实物面特征；物面边界相邻体网格单元相应构造为高阶体网格单元，其余体网格单元采用八节点六面体单元，以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解，验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外，由于采用了隐式计算方法，仅需较少的时间步就能迭代收敛。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

DPLR隐格式在多块结构网格的计算实现

利用DPLR隐格式收敛快、适于并行计算的特点,结合DPLR在多块结构网格的对接边界处理,提出了一种隐式对接边界条件数学模型与数值处理方法。以二维和三维球头绕流为算例,研究了网格分块方式对DPLR隐式算法稳定性和收敛性能的影响。通过返回舱复杂绕流问题数值模拟、气动特性检验分析以及与风洞试验对比,考察了算法在复杂外形数值模拟中的能力和高收敛性能特点:在非求解方向上进行网格分块,不会对DPLR算法的收敛性能产生影响;而在求解方向上进行网格分块,特别是分块位置在边界层内,会降低算法的稳定性和收敛性;提出的隐式对接边界条件处理方法,能改善在求解方向上进行网格分块造成的算法稳定性和收敛性能下降的问题。DPLR算法结合隐式对接边界条件能够成功应用于类返回舱外形体复杂流动数值模拟,且收敛速度较快。     

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

开放式气动力数值模拟系统研究

气动力数值模拟系统是ＣＦＤ流场解算技术、网格生成技术、数据可视化技术和网络技术相结合的产物。描述了一种运行于计算机网格环境下多任务、多用户、交互式的数值模拟软件环境－ 气动力数值模拟系统,该环境能计算模拟亚、跨、超音速绕流的气动力,也能方便地集成模拟不同复杂气动布局飞机流场的解算程序。介绍了其网格生成、数据可视化、系统接口等核心部分的结构和功能。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

一种边界非协调数值方法

研究了一种边界非协调的数值方法。该方法的提出是为了能够方便的求解复杂的运动边界绕流问题。空间上采用Galerkin有限体积空间离散格式。数值离散中,网格点分为内部计算点、IB点和外点三种类型。其中临近物面的IB点在嵌入物面边界条件时起到了作用,其流动变量值则是通过流场内解的线性插值获得。时间推进上采用结合人工压缩项的双时间步推进格式。静止和振荡圆柱不可压绕流的数值实验初步验证了本文方法的可靠性和准确性。     

裙后向台阶超声速绕流特性数值研究

本文研究了在中部带有裙体后向台阶的轴对称体的超声速绕流特性。采用空间推进法数值模拟超声速无粘流场，对于轴对称后向台阶的分离区采用Ｃｈａｐｍａｎ－Ｋｏｒｓｔ理论模型进行处理，计算结果给出了与实验值吻合的物面压力分布以及超声速流场的流动图画。