-
公开(公告)号:CN119312720A
公开(公告)日:2025-01-14
申请号:CN202411356403.X
申请日:2024-09-27
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/28 , G06F111/10 , G06F113/08 , G06F119/14
Abstract: 本发明公开了一种六阶精度自适应光滑尺度分离的多阶子模板非线性加权插值方法,涉及高阶精度有限差分格式的非线性加权领域,包括:获得四个多阶子模板,四个多阶子模板包括三个三阶精度子模板以及一个六阶精度子模板;分别计算获得四个多阶子模板的光滑度量因子;基于光滑度量因子构造获得光滑尺度分离器;基于六个模板点离散计算获得半节点处物理量的一阶导数、二阶导数、三阶导数和四阶导数;基于半节点处物理量的一阶导数、二阶导数、三阶导数和四阶导数构造自适应尺度识别器;使用所述自适应尺度识别器作为所述光滑尺度分离器的判定条件,分别计算获得四个多阶子模板的加权系数;基于四个多阶子模板的加权系数将四个多阶子模板加权计算获得求解界面的物理量值;本方法能够解决高阶格式构造中采用多阶子模板在非线性加权插值中遇到的在光滑区未能实现高阶精度插值计算,导致的计算结果可能不准确的问题,进而保障湍流数值模拟的精度。
-
公开(公告)号: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通量通过相邻单元的计算格式对应的插值多项式计算左右值。本发明不仅节省了计算量,而权利要求能够显著提高格式对激波的捕捉能力,并且在降低算法对侦测准确性的要求的同时,还能保证最终算法在激波附近仍然具有较高分辨率。
-
公开(公告)号:CN119026523B
公开(公告)日:2025-01-24
申请号:CN202411497810.2
申请日:2024-10-25
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
Abstract: 本发明公开了适用于多量阶非线性加权方法的光滑因子量阶调节方法,涉及高阶格式非线性加权领域,包括:基于精度K选取非等距模板集;分别计算获得非等距模板集中的3个子模板的光滑因子;确定最长子模板上的光滑因子量阶调节器;根据光滑因子量阶调节器和光滑因子,计算得到修正后的最长子模板上光滑因子;确定最短子模板上的光滑因子量阶调节器;根据光滑因子量阶调节器和光滑因子,计算得到修正后的最短子模板上光滑因子;获得3个子模板最终的光滑因子;本发明能够确保非线性加权高阶格式精度,且修正光滑因子计算的非线性权与控制方程的量纲无关。
-
公开(公告)号:CN118981891A
公开(公告)日:2024-11-19
申请号:CN202411072887.5
申请日:2024-08-06
Applicant: 中国空气动力研究与发展中心计算空气动力研究所
IPC: G06F30/20 , G06F30/15 , G06F30/28 , G06F113/08 , G06F119/14
Abstract: 本发明公开了一种飞行器跨流域气动特性高效预测方法,涉及气动特性预测领域,包括:步骤1:构建飞行器表面网格;确定飞行器物理空间网格的外边界范围;生成飞行器物理空间网格;步骤2:生成结构化速度空间网格,确定速度空间数值积分的权系数;步骤3:基于来流条件获得平衡态分布函数go;步骤4:基于飞行器物理空间网格、结构化速度空间网格和所述速度空间数值积分的权系数,以确定的平衡态分布函数go为行器物理空间网格中所有单元初始值,开始进行迭代求解,获得飞行器跨流域气动特性预测结果,本发明能够高效的获得飞行器跨流域气动特性。
-
公开(公告)号:CN114091376A
公开(公告)日:2022-02-25
申请号: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: 本申请公开了一种基于多重光滑度量因子的数值模拟方法及装置,通过获取待模拟的流场区域,并对待模拟的流场区域进行拆分,得到多个结构网格;根据结构网格,确定每一个结构网格对应的求解点位置和通量点位置;确定每一个求解点位置的节点的物理量数值;采用非线性加权高阶精度插值方法,根据节点的物理量数量,确定半节点的物理量数值;根据节点的物理量数值和半节点的物理量数值,确定节点处的数值通量;根据节点处的数值通量对待模拟的流场区域进行流场数值的模拟,通过本申请实施例,不需要借助多次加权组合;不需要借助逻辑判断实时调整理想权重的数值;能处理两个间断共存的情况,保持无振荡捕捉间断的特性。
-
-
-
-
-
-
-
-
-