[发明专利]考虑可压缩和热传导效应的壁面第一层网格厚度估算方法

说明书

本发明属于计算流体力学网格前处理技术领域，具体涉及考虑可压缩和热传导效应的壁面第一层网格厚度估算方法。其技术方案为：考虑可压缩和热传导效应的壁面第一层网格厚度估算方法，包括如下步骤：S1：给定计算工况；S2：计算出可压缩摩阻系数C_f；S3：分别计算出壁面无量纲粘性系数壁面无量纲密度和网格雷诺数Reue；S4：根据给定的y⁺值，求出壁面第一层网格厚度y。本发明提供了一种考虑可压缩和热传导效应的壁面第一层网格厚度估算方法。

技术领域

本发明属于计算流体力学网格前处理技术领域，具体涉及考虑可压缩和热传导效应的壁面第一层网格厚度估算方法。

背景技术

CFD数值模拟过程包括网格前处理、流场求解和后处理三部分。计算网格的合理设计和高质量生成是CFD计算可靠性的前提条件，随着复杂流动问题研究的进一步发展，CFD数值模拟对网格生成效率和质量提出了更高的要求，高质量网格是影响CFD计算的关键因素之一。

研究表明，近壁面区域根据流动状态的不同可以分为三层，自壁面向流动核心区分别为：粘性底层(viscous sublayer)、过渡层(buffer layer)和对数率层(log-lawlayer)。对于近壁区域的求解，主要集中在粘性底层的求解上，主要有两种方式，一是加密网格，使网格满足粘性计算需要；另一种是使用壁面函数法(wall function)。对于两种方法的选取，可以通过y⁺来体现。y⁺是与第一层网格厚度y和具体流动特征参数(密度ρ，壁面摩擦速度u_τ和粘性系数μ)有关的无量纲距离，表征了第一层网格高度。其定义为：

不同湍流模型对y⁺的分布有不同的要求。对于低雷诺数湍流模型(如k-ω模型、SA模型等)，不使用壁面函数，需要保证y⁺1为佳；需要用到壁面函数的湍流模型(如k-epsilon模型、雷诺应力模型等)，通常需要确保y⁺为30～200。

为了能在CFD数值模拟中准确计算实际工程所关心的摩阻、热流等物理量，同时控制网格数量以减少计算成本，可以根据y⁺来指导第一层网格的厚度。但所需变量u_τ涉及到壁面剪切应力τ_w的值，壁面剪切应力的计算又涉及到速度梯度的计算，但划分网格时流场还未知，由此陷入死循环。因此，在实际工程应用中，通常引入壁面摩阻系数C_f来计算出壁面剪切应力。

现行方法存在的问题是：计算出的粘性系数μ、流体密度ρ和摩阻系数C_f通常不考虑压缩效应及热传导。对于低马赫数绝热壁问题能保证计算的准确性，但对于高马赫数可压缩或者等温壁的误差很大，估算的网格高度甚至比实际需要的密几个量级，不能正确指导第一层网格的划分，大大浪费了计算资源及计算时间。

发明内容

为了解决现有技术存在的上述问题，本发明目的在于提供一种考虑可压缩和热传导效应的壁面第一层网格厚度估算方法。

本发明所采用的技术方案为：

考虑可压缩和热传导效应的壁面第一层网格厚度估算方法，包括如下步骤：

S1：给定计算工况：流场马赫数M_∞、雷诺数Re_L、边界温度T_∞、壁面温度与绝热壁面温度比值T_w/T_aw；以及湍流复温系数r、气体比热比γ；Sutherland粘性公式常数T₀，μ₀，S；实际尺寸与网格尺寸比值Cr，以及所需y⁺值；

S2：按以下公式计算出可压缩摩阻系数C_f；