-
公开(公告)号:CN116070071A
公开(公告)日:2023-05-05
申请号:CN202310347874.3
申请日:2023-04-04
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F17/15
Abstract: 本申请公开了一种壁面函数的选择方法、装置、终端设备及存储介质,通过分别使用基于y+的标准壁面函数和基于y*的可压缩修正壁面函数开展模拟及计算物面附近边界层内网格点上物面切向速度的物面法向偏导数,根据基于y+和y*的模拟结果一致性及物面附近边界层内网格点上物面切向速度的物面法向偏导数,确定湍流的流动类型;根据湍流的流动类型,选择与流动类型对应的目标壁面函数,通过本发明实施例,选择与流动类型对应的目标壁面函数,这样可以使得计算结果更加准确,本发明实施例能够在使用壁面函数模拟可压缩流动时自动判断流动类型,选用更适配的壁面函数类型,获得更准确的模拟结果。
-
公开(公告)号:CN114611438B
公开(公告)日:2022-07-22
申请号:CN202210506807.7
申请日:2022-05-11
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明实施例提供了一种湍流流动中目标物的受力状态模拟方法及装置,方法包括:基于目标物的贴体结构网格的网格导数和网格雅可比,迭代求解流动控制方程和湍流模型方程,迭代过程收敛后得到目标流场变量;对目标流场变量进行积分,得到目标物在湍流状态下的受力状态模拟结果;本方案中,同步迭代求解流动控制方程和湍流模型方程,直到收敛后得到湍流状态下的受力状态模拟结果,提高了湍流状态下的模拟精度。
-
公开(公告)号:CN114091376B
公开(公告)日:2022-04-12
申请号:CN202210069705.3
申请日:2022-01-21
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明公开了一种基于子单元加权格式的高精度重构修正激波捕捉方法,属于激波捕捉技术领域,包括步骤:S1,高阶CPR的求解点和通量点的分布以及通量导数的计算;S2,通过激波侦测方法侦测出存在激波的单元;S3,将单元基于通量点剖分为子单元,在子单元上构造具有激波捕捉能力的格式;S4,基于激波侦测方法中的激波侦测因子的值来确定CPR、HNNW和二阶格式的计算单元;S5,单元界面处的Riemann通量通过相邻单元的计算格式对应的插值多项式计算左右值。本发明不仅节省了计算量,而权利要求能够显著提高格式对激波的捕捉能力,并且在降低算法对侦测准确性的要求的同时,还能保证最终算法在激波附近仍然具有较高分辨率。
-
公开(公告)号:CN116882322B
公开(公告)日:2024-02-13
申请号:CN202311141821.2
申请日:2023-09-06
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F113/08 , G06F119/14
Abstract: 本申请公开了一种无粘通量的计算方法、装置、终端设备和存储介质,通过获取待计算的网格单元,并根据网格单元,确定网格单元的求解点和通量点;根据预设初始条件,确定通量点对应的物理信息;根据通量点对应的物理信息和预设的分裂形式,计算通量点对应的无粘通量;根据通量点对应的无粘通量,确定求解点处的守恒形式下的无粘通量的一阶导数和求解点处的非守恒形式下的无粘通量的一阶导数;确定网格单元的交界面处的通量信息,确定守恒形式下的无粘通量的修正信息;并进一步确定待计算的网格单元的无粘通量的离散项;并确定待计算的网格单元的目标无粘通量,不仅提高了计算稳定性,还满足离散守恒律。(56)对比文件Huajun Zhu 等.A viscous-term subcelllimiting approach for high-order FR/CPRmethod in solving compressible Navier-Stokes equations on curvilineargrids.Journal of ComputationalPhysics.2023,第1-38页.Sebastian Hennemann 等.A provablyentropy stable subcell shock capturingapproach for high order split form DG forthe compressible Euler equations.Journalof Computational Physics .2020,第426卷第109935页.Shuai Lou 等.The Flux ReconstructionMethod with Lax–Wendroff Type TemporalDiscretization for HyperbolicConservation Laws.Journal of ScientificComputing.2020,第82卷第42页.刘冉 等.一种基于过积分的能量稳定通量重构方法.计算数学.2023,第45卷(第03期),第368-384页.李津 等.带副翼三维机翼绕流的N-S方程解.航空学报.1999,(第03期),第197-200页.崔树鑫 等.基于CE/SE法的翼型流场激波捕捉计算.北京航空航天大学学报.2012,第38卷(第12期),第1591-1595页.
-
公开(公告)号:CN116611369B
公开(公告)日:2023-09-29
申请号:CN202310872199.6
申请日:2023-07-17
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F17/10 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本申请公开了基于光滑度量量级与候选模板点数的插值方法及装置,通过获取待模拟的流场区域,并对待模拟的流场区域进行拆分,得到多个结构网格;确定目标插值模板,第一数量的候选模板对应的插值精度是第一阶光滑度量因子,全局第二数量的候选模板对应的插值精度是第二阶光滑度量因子;光滑度量因子用于衡量候选模板光滑程度;通过目标插值模板,采用非线性加权高阶精度格式对Navier‑Stokes方程进行求解,新的候选模板集和光滑度量因子计算新方法,实现了新的非线性加权机制,在流场光滑区域有两重机制共同作用来促使加权组合格式逼近线性格式,相比于原有方法分辨率更高而在光滑区耗散水平更低。
-
公开(公告)号:CN116702571A
公开(公告)日:2023-09-05
申请号:CN202310982965.4
申请日:2023-08-07
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/23 , G06F30/28 , G06F113/08 , G06F119/14 , G06F111/10
Abstract: 本申请公开了一种基于多重光滑度量因子的数值模拟方法及装置,通过获取待模拟的流场区域,并对待模拟的流场区域进行拆分,得到多个结构网格;根据结构网格,确定每一个结构网格对应的求解点位置和通量点位置;确定每一个求解点位置的节点的物理量数值;采用非线性加权高阶精度插值方法,根据节点的物理量数量,确定半节点的物理量数值;根据节点的物理量数值和半节点的物理量数值,确定节点处的数值通量;根据节点处的数值通量对待模拟的流场区域进行流场数值的模拟,通过本申请实施例,不需要借助多次加权组合;不需要借助逻辑判断实时调整理想权重的数值;能处理两个间断共存的情况,保持无振荡捕捉间断的特性。
-
公开(公告)号:CN116542184A
公开(公告)日:2023-08-04
申请号:CN202310817635.X
申请日:2023-07-05
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F17/11 , G06F113/08 , G06F119/14
Abstract: 本申请公开了一种粘性通量的计算方法、装置、终端设备和存储介质,通过获取待计算的网格单元,并根据网格单元,确定网格单元的求解点、通量点和辅助点;采用交错直接求导方法,根据辅助点处的物理量和通量点处的物理量,计算各个通量点的修正后一阶导数;根据各个通量点的修正后一阶导数,确定各个通量点处的粘性通量;根据各个通量点的粘性通量,确定求解点处的粘性通量导数;根据各个求解点的粘性通量导数,确定待计算的网格单元求解点处的修正后粘性通量导数,实现带激波问题的NS方程的稳定化粘性通量计算,提高了含激波的粘性问题计算的鲁棒性。
-
公开(公告)号:CN116070071B
公开(公告)日:2023-06-27
申请号:CN202310347874.3
申请日:2023-04-04
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F17/15
Abstract: 本申请公开了一种壁面函数的选择方法、装置、终端设备及存储介质,通过分别使用基于y+的标准壁面函数和基于y*的可压缩修正壁面函数开展模拟及计算物面附近边界层内网格点上物面切向速度的物面法向偏导数,根据基于y+和y*的模拟结果一致性及物面附近边界层内网格点上物面切向速度的物面法向偏导数,确定湍流的流动类型;根据湍流的流动类型,选择与流动类型对应的目标壁面函数,通过本发明实施例,选择与流动类型对应的目标壁面函数,这样可以使得计算结果更加准确,本发明实施例能够在使用壁面函数模拟可压缩流动时自动判断流动类型,选用更适配的壁面函数类型,获得更准确的模拟结果。
-
公开(公告)号:CN116151084A
公开(公告)日:2023-05-23
申请号:CN202310434862.4
申请日:2023-04-21
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/23 , G06F30/28 , G06T17/20 , G06F17/11 , G06F113/08 , G06F111/10 , G06F119/14
Abstract: 本申请公开了一种基于结构网格的模拟方法、装置、终端设备及存储介质,通过获取待模拟的流场区域,并根据待模拟的流场区域确定结构网格单元;根据坐标变换规则,将直角坐标系下的流动控制方程转换成计算坐标系下的流动控制方程,根据计算坐标系下的流动控制方程,确定各个结构网格单元的网格特征信息;根据各个结构网格单元的网格特征信息,采用有限差分方法,以逐维离散求解的方式,确定各个结构网格单元的流场数据,并将待模拟的流场区域的各个结构网格单元的流场数据按照预设规则进行显示输出,在几乎不增加计算量的情况下,能够大幅提升数值方法对固体壁面等边界附近实际应用中重点关注的关键流动区域的模拟能力。
-
公开(公告)号:CN114611438A
公开(公告)日:2022-06-10
申请号:CN202210506807.7
申请日:2022-05-11
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F30/23 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明实施例提供了一种湍流流动中目标物的受力状态模拟方法及装置,方法包括:基于目标物的贴体结构网格的网格导数和网格雅可比,迭代求解流动控制方程和湍流模型方程,迭代过程收敛后得到目标流场变量;对目标流场变量进行积分,得到目标物在湍流状态下的受力状态模拟结果;本方案中,同步迭代求解流动控制方程和湍流模型方程,直到收敛后得到湍流状态下的受力状态模拟结果,提高了湍流状态下的模拟精度。
-
-
-
-
-
-
-
-
-