-
公开(公告)号:CN112685978A
公开(公告)日:2021-04-20
申请号:CN202110264089.2
申请日:2021-03-11
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F113/08 , G06F119/14
Abstract: 本发明提供了一种适用于五次样条重构格式的自适应人工粘性控制方法,包括:步骤1、读取流场数据,求解五次样条重构方程组得到单元界面处物理量的状态值及其一阶导数;步骤2、计算出单元界面处的二阶至五阶导数;步骤3、根据波数识别方法,求出单元界面处流场的等效无量纲波数,从而确定该处流场的人工粘性系数;步骤4、根据单元界面处物理量的状态值计算无粘数值通量,并根据人工粘性系数添加六阶人工粘性项,最后采用相对应的时间离散方法在时间上进行推进。将局部瞬时流场特征与格式的解析能力结合起来,通过判断局部瞬时的流场特征是否在格式的可解析范围内,从而确定格式的人工粘性系数的取值,使格式的人工粘性控制更为合理。
-
公开(公告)号:CN119761266A
公开(公告)日:2025-04-04
申请号:CN202510259175.2
申请日:2025-03-06
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/15 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明属于计算流体力学技术领域,公开了一种面向超声速激波湍流干扰的数值模拟方法,包括如下步骤:S100初始化雷诺平均纳维斯托克斯方程中流动物理量,以及#imgabs0#SST湍流模型方程中的湍流特征量;S200求解RANS方程,得到第#imgabs1#时刻的#imgabs2#;根据S200中得到的相关值,计算尺度规约化的无量纲逆压梯度量#imgabs3#;根据S300中得到的#imgabs4#值,计算调节函数#imgabs5#,得到改进的湍动能产生项#imgabs6#;根据S400中得到的#imgabs7#,求解改进的#imgabs8#SST湍流模型方程,得到#imgabs9#;根据S500中得到的#imgabs10#代入S200,迭代求解,直到当RANS方程残差Resn小于一定值ε或达到最大步数N=Nmax时,则停止继续迭代求解。根据S600输出最新流动物理量的平均值#imgabs11#,结束数值模拟。
-
公开(公告)号:CN116384288A
公开(公告)日:2023-07-04
申请号:CN202310654382.9
申请日:2023-06-05
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F17/11 , G06F113/08 , G06F119/14 , G06F111/10
Abstract: 本申请涉及计算流体力学中数值计算方法领域,公开了一种可压缩流动高分辨率数值模拟方法、介质及设备,包括:读取初始流场数据,得到当前时刻计算空间中各网格点的变量和通量;根据得到的变量,采用非线性插值计算半点处的左值和右值,并根据左值和右值,计算半点处的低阶通量;根据得到的通量,计算半点处的修正通量;根据低阶通量和修正通量,获取计算空间中不同方向各节点处的一阶导数,以完成空间离散;采用龙格‑库塔法将当前时刻的变量推进到下一时刻的变量,以完成时空离散;将时间推进至设定时刻结束计算,得到设定时刻的流场数据。这样可以得到合理的可压缩流动数值模拟结果,稳定捕捉相应的可压缩流动,显著提高对流动结构的分辨能力。
-
公开(公告)号:CN114139465B
公开(公告)日:2023-04-07
申请号:CN202111239727.1
申请日:2021-10-25
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明提供了一种脱体涡模拟模型构造方法,包括以下步骤:步骤1、获取RANS模型控制方程中的湍动能破坏项,确定湍动能破坏项中的长度尺度;步骤2、对长度尺度进行修正,完成新的替代破坏项构建;步骤3、对新的替代破坏项进行网格尺度修正,完成脱体涡模拟模型的构造。本发明提出的模型构造方法,在模拟含分离的湍流流动时,能够更合理地自动区分需要使用RANS的近壁区域和需要使用LES的远离壁面区域,使对流动的解析更加精准,能够在较为苛刻的网格分辨率条件下得到最优的数值模拟结果。
-
公开(公告)号:CN117682060A
公开(公告)日:2024-03-12
申请号:CN202410142963.9
申请日:2024-02-01
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: B64C23/00
Abstract: 本发明公开了一种等离子体激励控制激波边界层干扰的减阻装置,包括:第一平板、第二平板与放电电极;第一平板与第二平板固定连接,第一平板在远离其与第二平板连接处的一端,设置有转捩带;且第一平板所在的平面与第二平板所在的平面存在夹角α;第一平板上设置有多个圆柱孔,圆柱孔用于放置放电电极;其中,圆柱孔每两个形成一组,分别用于放置正负放电电极,每组圆柱孔呈交错阵列设置;圆柱孔设置有2*N*(2M+1)个,并在第一平板上按照从转捩带往第二平板的方向依次排布。本发明通过提供一种交错式布局的阵列式等离子体激励器,在消耗同样能量的情况下,可以实现更好的流动控制效果,可以显著降低摩擦阻力,压阻也能进一步降低。
-
公开(公告)号:CN116384288B
公开(公告)日:2023-08-25
申请号:CN202310654382.9
申请日:2023-06-05
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F17/11 , G06F113/08 , G06F119/14 , G06F111/10
Abstract: 本申请涉及计算流体力学中数值计算方法领域,公开了一种可压缩流动高分辨率数值模拟方法、介质及设备,包括:读取初始流场数据,得到当前时刻计算空间中各网格点的变量和通量;根据得到的变量,采用非线性插值计算半点处的左值和右值,并根据左值和右值,计算半点处的低阶通量;根据得到的通量,计算半点处的修正通量;根据低阶通量和修正通量,获取计算空间中不同方向各节点处的一阶导数,以完成空间离散;采用龙格‑库塔法将当前时刻的变量推进到下一时刻的变量,以完成时空离散;将时间推进至设定时刻结束计算,得到设定时刻的流场数据。这样可以得到合理的可压缩流动数值模拟结果,稳定捕捉相应的可压缩流动,显著提高对流动结构的分辨能力。
-
公开(公告)号:CN116306025A
公开(公告)日:2023-06-23
申请号:CN202310534776.0
申请日:2023-05-12
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/20 , G06F30/15 , G06F30/18 , G06F111/10 , G06F119/14
Abstract: 本申请公开了一种湍流生成方法、装置、设备及存储介质,涉及湍流生成领域,该方法包括:对参考网格模拟的参考湍流脉动进行特征模态分解得到参考流场的空间模态和时间系数;基于随机生成的时间窗口的起始时刻与目标流场时刻的相对位置关系对目标流场时刻增加时间变化量得到调整后的目标时刻,遍历时间窗口从参考流场确定第一参考时刻和相邻第二参考时刻;基于重构特征模态数、时间系数和空间模态重构第一、第二参考时刻内的参考流场;对重构后参考流场依次进行时间插值、二维线性插值以及速度和温度脉动调整生成与目标时刻对应的目标流场的湍流脉动。本申请通过特征模态分解提取空间模态和时间系数,缩短过渡段长度,提高湍流生成的效率和准确率。
-
公开(公告)号:CN114139465A
公开(公告)日:2022-03-04
申请号:CN202111239727.1
申请日:2021-10-25
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明提供了一种脱体涡模拟模型构造方法,包括以下步骤:步骤1、获取RANS模型控制方程中的湍动能破坏项,确定湍动能破坏项中的长度尺度;步骤2、对长度尺度进行修正,完成新的替代破坏项构建;步骤3、对新的替代破坏项进行网格尺度修正,完成脱体涡模拟模型的构造。本发明提出的模型构造方法,在模拟含分离的湍流流动时,能够更合理地自动区分需要使用RANS的近壁区域和需要使用LES的远离壁面区域,使对流动的解析更加精准,能够在较为苛刻的网格分辨率条件下得到最优的数值模拟结果。
-
公开(公告)号:CN113468679B
公开(公告)日:2021-11-12
申请号:CN202111036423.5
申请日:2021-09-06
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
Abstract: 本发明提供了一种基于S‑A模型的湍流长度尺度计算方法,包括:步骤1、获取湍动能的理论运输方程中的生成项与耗散项;步骤2、获取S‑A模型控制方程的生成项与破坏项;步骤3、建立名义粘性系数和湍动能之间的近似关系,类比得到S‑A模型中湍动能的生成项和破坏项;步骤4、联立理论运输方程中的生成项与耗散项与S‑A模型中湍动能的生成项和破坏项得到湍流长度尺度代数表达式;步骤5、使用RANS模型对壁湍流进行模拟,得到湍流流动的时均结果,结合湍流长度尺度代数表达式完成湍流长度尺度的计算。本发明能够加快入口边界的人工合成壁湍流向真实壁湍流的发展过程,最终达到提高整个DNS或LES计算精度的目的。
-
公开(公告)号:CN119720860A
公开(公告)日:2025-03-28
申请号:CN202510193236.X
申请日:2025-02-21
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06T17/20 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本申请公开了基于尺度滤波的格式耗散调节方法、装置、设备及介质,涉及流体动力学技术领域,对流场数据中的速度场进行特征分析,并投影至特征空间,得到速度场的特征变量;对特征变量进行移动最小二乘滤波;将第一滤波后的特征变量及第二滤波后的特征变量反投影至物理空间,得到第一滤波后的速度场及第二滤波后的速度场;计算第一高频分量平方差及第二高频分量平方差;对第一高频分量平方差及第二高频分量平方差求和,得到第一求和结果及第二求和结果,计算能量比;建立能量比与耗散调节因子的映射关系,计算耗散调节因子;利用耗散调节因子对格式的耗散进行自适应调节,解决耗散调节过程依赖于具体格式和具体问题,实现格式耗散的自适应调节。
-
-
-
-
-
-
-
-
-