热流CFD计算中格式和网格效应若干问题研究

通过数值实验、理论分析、试验数据和经验公式的综合研究,探讨了CFD热流计算的格式效应及网格效应:研究了Roe的FDS格式、Van Leer的FVS格式、AUSM 格式及中心差分格式的热流计算分辨率,探讨了不同格式对热流计算精度的影响,认为AUSM 格式在热流计算精确性方面具有优势;在热流绝对值的模拟上,网格的影响很大,能否计算出正确的驻点热流值将是获得精确物面热流解的关键;通过大量的数值计算、理论分析及其同试验数据的对比研究,对可能影响热流计算的网格因素进行了细致的分析,揭示了各因素对热流计算结果的影响规律。     

飞船外形高超音速流场N-S方程数值解

用二阶迎风TVD格式求解N-S访方程,对简化的飞船外形高超音速粘性绕流进行了数值模拟。给出了物面压力分布及流场结构,与实验及其他计算结果符合较好。计算表明,肩部圆弧半径对后体流场影响较小,但严重影响肩部附近热流分布。当肩部圆弧半径是头部半径的20％时,肩部热流低干驻点处热流;半径为5％时,肩部热流超过驻点处热流30％。     

非交错网格求解传热及流体流动

在非交错网格系统上提出了一种新的以控制容积为基础的有限差分方法，来求解传热及流体流动问题。为了避免压力棋盘问题，压力梯度采用偏迎风压力离散格式进行处理，以保持压力－速度的耦合。利用本算法对数值传热的典型问题进行计算，并将计算结果与采用交错网格的PHOENICS通用程序所得结果相比较，二者吻合甚好。     

火星大气非平衡条件下驻点热流计算公式的探讨

针对目前绝大多数火星再入飞行器驻点热流公式未考虑热力学非平衡效应故而适用性及准确性存疑的问题,分析了典型的零攻角驻点热流公式Fay-Riddell公式在火星再入飞行器热力学非平衡驻点热流计算中的局限性;然后采用数值计算方法,对涵盖火星再入飞行器飞行状态的大范围计算状况开展了热力学非平衡模拟,并对模拟得到的驻点热流公式进行拟合,得到了适用于火星再入飞行器热力学非平衡条件下的零攻角驻点热流计算公式;利用Mars Pathfinder飞行数据对公式进行了验证,计算表明,对于验证的计算状态,该公式的计算误差小于10%,符合工程预估的可接受误差范围。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

非均匀笛卡尔坐标系网格优化频散相关保持格式研究

为克服传统频散相关保持格式 (即 DRP格式 )在数值求解非均匀网格问题时由于寄生波的产生而导致的数值不稳定的问题 ,Cheong &Lee提出了一种网格优化频散相关保持格式 (即 GODRP格式 )。根据其优化思想 ,在非均匀笛卡尔坐标系网格条件下 ,推导得到了更为通用的 GODRP格式优化系数公式。针对 Cheong&Lee文中算例比较不是在同一基准上进行的不足之处 ,设计了一个极端恶劣条件下的非均匀笛卡尔网格 ,采用 GODRP格式和 DRP格式分别求解了二维初始扰动波传播的算例。计算结果表明 ,由于 GODRP格式具有与网格相关的自带粘性 ,不用再另外加入人工粘性便可获得精度较高的数值结果 ,避免了传统 DRP格式在给定粘性时人为经验的影响因素。对于非均匀笛卡尔网格 ,GODRP格式确实比传统 DRP格式更加通用 ,并可以得到更为准确的数值模拟效果      

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

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

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

航天器再入大气层热力分析

以OREX（orbital reentry experiment vehicle）的飞行试验数据和相关的CFD数值模拟结果为基础,采用传热理论及相关公式,分析计算了OREX再入大气层过程中的轨迹,驻点处热流密度非平衡假设和平衡假设下的换热问题.所计算的轨迹、热流密度非平衡假设下计算的驻点温度和热流密度值同试验数据及相关的CFD数值结果取得了很好的一致,相应的计算方法可作为航天器驻点热力分析的通式.然而在平衡假设条件下,尽管计算所得到的驻点热流密度与之前的CFD数值结果差别不大,但这种情况下计算得到的温度与试验数据不符,这应归结于计算的热流密度结果对驻点处温度变化的不敏感.比较非平衡假设和平衡假设下的换热计算结果表明,对于航天器再入过程中的热力探讨不能仅仅满足于热流密度分析,对温度的考察或许更重要.      

基于流线跟踪法的气动热工程计算研究

运用结构化网格求解三维Euler方程,计算得到边界层外缘无粘流场气流参数;利用无粘流场气流参数和表面流函数的方法计算了飞行器无粘表面流线分布;在理论和半经验公式的基础上,计算了定比热比和变比热比情况下驻点热流密度,非驻点区域采用参考焓、局部相似性等方法来确定飞行器表面的气动加热,实现了数值算法与工程算法的耦合.上述方法用于求解高超声速钝双锥的表面热流分布,计算结果与经典的热流公式和实验结果进行对比,平均精度为10%左右,满足高超声速飞行器概念研究和初步设计的需要.