[发明专利]基于归一化物理量间断特征的喷流模拟激波快速判别方法

说明书

本发明涉及计算流体力学领域，公开了一种基于归一化物理量间断特征的喷流模拟激波快速判别方法，该方法采用物理量周围单元的平均物理量对网格单元物理量进行归一化缩放，消除了喷流周围物理量和流场物理量量级差异带来的梯度计算误差，通过间断因子和激波判别指标的引入，建立了物理量间断分布特征的量化表征，实现了含有喷流流场的激波判别。本发明提出的方法避免了传统方法梯度求解，仅依赖局部相邻区域的流场变量；该方法降低了激波判别算法的计算和存储开销，同时防止了喷流出口处大梯度区域的激波误判，为喷流入射高速流动复杂流场提供了一种新的快速激波判别方法。

技术领域

本发明涉及计算流体力学领域，具体涉及一种基于归一化物理量间断特征的喷流模拟激波快速判别方法。

背景技术

激波是超声速飞行器流动的典型物理特征，激波判别是计算流体力学流动显示和算法设计等领域涉及的一项核心技术。当前普遍采用的激波判别方法是流场梯度方法，该方法除了要用到网格单元上的物理量分布之外，还需要基于几何量计算物理量权值并采用较宽的模板对单元的梯度进行重构，由此带来了基于流场梯度方法计算量大的缺点。除此之外，梯度的计算包含了物理量本身的差量和网格相邻几何曲率的影响，在网格相邻几何曲率较大（如驻点或前缘区域）的情况下也会得到较大的梯度而带来激波的误判。更严重的是，基于流场梯度的方法给包含喷流入射的流场计算带来了很大的问题，因为喷流出口附近物理量相对于来流物理量基本上大了数百个量级，按照传统梯度方法计算，在网格尺度基本一致的情况下，喷流出口附近也会计算得到较大的流场梯度，该梯度与激波处梯度数值比较相近，因此采用梯度的判别方法在数学原理上无法对二者进行区分，亟需寻找新的适用于包含喷流入射情况下流场的高效激波判别方法。

发明内容

本发明所要解决的技术问题是：针对上述存在的问题，提供了一种基于归一化物理量间断的喷流激波快速判别方法，该方法考虑到喷流附近区域物理量（压力等）梯度较大是由于其物理量本身相对于来流的放大效应导致的，因此为了消除这种数值放大效应对喷流区域梯度计算的影响，从利用周围物理量特征进行归一化的思想出发，通过间断因子的引入，实现了物理量间断分布特征的量化表征，进而用于激波的判别。

本发明采用的技术方案如下：

步骤1：给定网格单元上的流场变量值为f，针对网格单元i和与其共面的相邻网格单元j，利用网格单元i和其共面的相邻网格单元j的物理量分布平均特征，对流场变量进行归一化缩放：

进行缩放后的物理量分布变化在整个流场区域上不会因为绝对值的大小给间断的计算带来误判。

然后通过归一化后的物理量，计算网格单元的方向间断因子：

通过方向间断因子，来代替传统的梯度求解。因为物理量经过归一化后，物理量本身的变化就表征了流场中物理量的变化剧烈程度，因此可以省去梯度的求解，而直接采用间断因子表征。

步骤2：遍历网格单元i附近所有个和其共面的相邻网格单元，将附近最大的方向间断因子作为网格单元i的间断因子：

。

步骤3：针对具体的计算构型，统计每一个网格单元的间断因子，构成间断因子集合；

在间断因子集合中任意找出一个间断因子，判断该间断因子对应的网格单元是否满足以下要求：

既能标记出所关注的激波位置，又能避免喷流出口区域大梯度区的激波误判；

若满足要求，则将该间断因子标定为临界因子，若不满足要求，则继续重新寻找，直到找出一个满足要求的间断因子，并将其标定为临界因子。