CN115455864B 一种模拟全流域流动的高效数值方法、系统、装置及介质 （西北工业大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN115455864B(45)授权公告日2025.07.04

(21)申请号202211254874.0

(22)申请日2022.10.13

(65)同一申请的已公布的文献号申请公布号CN115455864A

(43)申请公布日2022.12.09

G06F119/14(2020.01)

(56)对比文件

CN112613243A,2021.04.06CN110348059A,2019.10.18审查员曹青

(73)专利权人西北工业大学

地址710072陕西省西安市碑林区友谊西

路127号

(72)发明人刘沙陈健锋钟诚文卓丛山

(74)专利代理机构西安通大专利代理有限责任公司61200

专利代理师张宇鸽

(51)Int.CI.

GO6F30/28(2020.01)

G06F113/08(2020.01)权利要求书3页说明书12页附图5页

(54)发明名称

一种模拟全流域流动的高效数值方法、系统、装置及介质

(57)摘要

CN115455864B本发明公开了一种模拟全流域流动的高效数值方法、系统、装置及介质，包括：获取飞行器外形几何参数和气体流动信息；构建飞行器的物理网格；基于物理网格中所存储的宏观量，获取速度网格的平衡态分布函数；基于物理网格所存储的宏观量和速度网格的分布函数，获取自适应准则值；将自适应准则值与速度网格自适应的粗化阈值和细化阈值分别进行比较，判断是否进行速度网格自适应；对重构后的速度网格中单元格心分布函数再次进行更新，获取更新后的流场物理量；判断流场是否进行收敛；对收敛后的流场进行处理，获取流动特性参数。本发明避免了在物理网格和对应的速度网格交换信息时，通过插

CN115455864B

第一获取模块第二获取模块更新

第一获取模块

第二获取

模块

更新模块

构建模块

判断模块

CN115455864B权利要求书1/3页

2

1.一种模拟全流域流动的高效数值方法，其特征在于，包括：

步骤1:获取飞行器外形几何参数和气体流动信息；基于所获取的飞行器外形几何参数和气体流动信息，构建飞行器的物理网格；

步骤2:基于物理网格中所存储的宏观量，获取速度网格的平衡态分布函数；

具体为：

(1)

其中，U,R,T和U分别是离散速度、气体常数、温度和速度；所述宏观量包括气体的密度、速度和温度；

步骤3:基于物理网格所存储的宏观量和速度网格的分布函数，获取若干个数值不等的自适应准则值；选取若干个自适应准则值中最大数值的自适应准则值，将选取出的自适应准则值与速度网格自适应的粗化阈值和细化阈值分别进行比较，判断是否进行速度网格自适应，若是，重构速度网格中的单元格心分布函数；若否，不进行任何操作；其中，初始时刻的速度网格分布函数等于速度网格的平衡态分布函数；

所述基于物理网格所存储的宏观量和速度网格的分布函数，获取若干个自适应准则值；

具体为：

M1,m?=Aμfim/P(2)

(3)

其中，Au为速度网格单元uk的面积或体积，Uk为第k个速度网格单元；fi,u为格心分布

函数；Pi、Ui、pE分别为第1个物理网格单元中气体的密度、速度、动能和总能；步骤4:基于物理网格所存储的宏观量、速度网格的平衡态分布函数和速度网格的分布

函数，对重构后的速度网格中单元格心分布函数再次进行更新，获取更新后的流场物理量；步骤5:基于经过若干次的更新，基于更新后的流场物理量，判断流场是否进行收敛；若

是，执行步骤6;若否，跳转至步骤3和步骤4,直至流场收敛；

步骤6:对收敛后的流场进行处理，获取流动特性参数。

2.根据权利要求1所述的模拟全流域流动的高效数值方法，其特征在于，自适应准则值的个数为物理网格单元的个数乘以2;每个物理网格单元包含两个自适应准则值。

3.根据权利要求2所述的模拟全流域流动的高效数值方法，其特征在于，所述选取若干个自适应准则值中最大数值的自适应准则值，具体为：

M?=max(Mu,M?m)