基于流场伴随方程的笛卡尔网格自适应加密方法
申请人信息
- 申请人:中国空气动力研究与发展中心计算空气动力研究所
- 申请人地址:621052 四川省绵阳市涪城区二环路南段6号
- 发明人: 中国空气动力研究与发展中心计算空气动力研究所
专利详细信息
| 项目 | 内容 |
|---|---|
| 专利名称 | 基于流场伴随方程的笛卡尔网格自适应加密方法 |
| 专利类型 | 发明授权 |
| 申请号 | CN202311804126.X |
| 申请日 | 2023/12/26 |
| 公告号 | CN117473907B |
| 公开日 | 2024/3/15 |
| IPC主分类号 | G06F30/28 |
| 权利人 | 中国空气动力研究与发展中心计算空气动力研究所 |
| 发明人 | 陈波; 陈浩; 毕林; 黄江涛; 牟斌; 宋超; 罗骁; 刘杨; 庞宇飞; 齐龙 |
| 地址 | 四川省绵阳市涪城区二环路南段6号 |
摘要文本
中国空气动力研究与发展中心计算空气动力研究所取得“一种透气窗帘布”专利技术，本发明提供了基于流场伴随方程的笛卡尔网格自适应加密方法，包括：生成初始笛卡尔网格；设置目标函数的全局误差阈值；在当前笛卡尔网格上求解三维可压缩无粘流动控制方程；在当前笛卡尔网格上求解流场伴随方程；对当前笛卡尔网格进行预加密，采用插值方法计算笛卡尔细网格上的流场插值解和伴随插值解；计算当前笛卡尔网格单元的伴随自适应探测器值，标记需要加密的当前笛卡尔网格单元；对待加密的笛卡尔网格单元进行加密；当所有网格单元的伴随自适应探测器值满足设定阈值要求时终止自适应加密过程。能够准确定位流场中输出函数误差较大的区域，确保网格加密在正确的位置进行, 有效避免了因错误判断加密区域导致的网格过度加密。
专利主权项内容
1.一种基于流场伴随方程的笛卡尔网格自适应加密方法，其特征在于，包括以下步骤：导入计算数模，生成初始笛卡尔网格；设置笛卡尔网格目标函数的全局误差阈值；在当前笛卡尔网格上采用格心格式有限体积法求解三维可压缩无粘流动控制方程；根据三维可压缩无粘流动控制方程的流场解在当前笛卡尔网格上求解流场伴随方程；对当前笛卡尔网格进行预加密得到笛卡尔细网格，采用插值方法计算笛卡尔细网格上的流场插值解和伴随插值解；计算当前笛卡尔网格单元的伴随自适应探测器值，根据自适应判据，标记需要加密的当前笛卡尔网格单元；对标记的待加密的笛卡尔网格单元进行加密；将剖分成细网格后新增加的网格点投影到物面上优化物面网格单元质量；判断当前笛卡尔网格中是否所有网格单元的伴随自适应探测器值满足设定阈值要求，当所有网格单元的伴随自适应探测器值满足设定阈值要求时终止自适应加密过程；所述三维可压缩无粘流动控制方程为：其中，表示方程的解向量，定义为流场守恒变量，即通过三维可压缩无粘流动控制方程求得的流场解；/>为流场分别在x、y、z方向上的对流通量；t为时间；所述在当前笛卡尔网格上采用格心格式有限体积法求解三维可压缩无粘流动控制方程，包括：对流通量的离散采用耦合人工粘性的Jameson中心格式，人工粘性通量包括二阶耗散和四阶耗散，时间推进采用显式四步Runge-Kutta格式；所述流场伴随方程为：其中，为当前笛卡尔网格上的伴随变量，即通过流场伴随方程求得的伴随解；为虚拟时间项，/>为虚拟时间，/>为当前笛卡尔网格单元的体积；/>为当前笛卡尔网格上流场伴随方程的残差；/>为当前笛卡尔网格上三维可压缩无粘流动控制方程的残值，/>为当前笛卡尔网格上的流场守恒变量，/>为当前笛卡尔网格上的目标函数。。数据由马 克 数 据整理