CN112665820B - 基于变量差及相对位移的r型网格自适应移动方法及设备 - Google Patents
基于变量差及相对位移的r型网格自适应移动方法及设备 Download PDFInfo
- Publication number
- CN112665820B CN112665820B CN202110275460.5A CN202110275460A CN112665820B CN 112665820 B CN112665820 B CN 112665820B CN 202110275460 A CN202110275460 A CN 202110275460A CN 112665820 B CN112665820 B CN 112665820B
- Authority
- CN
- China
- Prior art keywords
- grid
- grid point
- pressure
- points
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Abstract
本发明公开了一种基于变量差及相对位移的r型网格自适应移动方法及设备,以归一化网格点的变量差以及网格点相对位移量这两类参数作为计算网格点移动的权函数,通过加入网格点相对位移量这一权函数,实现网格点的跟随移动。归一化消除所述两类参数的量值及量纲对网格点移动量影响,以所述两类参数形成新的组合权值,以此确定激波附近网格点的位移,并将激波附近网格点的位移信息向网格点的相邻网格点传递,实现激波附近的几排网格点都同步向激波位置移动。本发明通过归一化物理量以及网格点位置移动量,消除两类参数量值及量纲对网格点移动量影响,可有效实现激波处网格点向激波位置靠近,同时能够保证激波附近的几排网格点都同步向激波位置移动。
Description
技术领域
本发明涉及计算流体力学技术领域,尤其涉及一种基于变量差及相对位移的r型网格自适应移动方法及设备。
背景技术
当前计算流体力学中的网格生成工作已成为工程师数值模拟飞行器内外流场时的主要工作之一。网格分布合理性直接影响流场分辨率以及模拟精度,因此CFD工程师在生成网格时需要对当前模拟工况的流场特征有较高认识和理解,包括激波位置、分离区位置和大小等。不同CFD工程师生成的网格分布不一致,导致数值模拟结果有差别,如何减弱网格分布对计算结果的影响是计算流体力学的一个重要研究内容。
网格自适应作为一种依据流场特征自动调整网格分布的方法,通过在数值模拟过程中不断动态调整初始生成的计算网格的分布,实现在流场突变区域的网格自动加密来捕捉流场中的复杂结构,提高流场分辨率和计算精度,能够减弱网格分布对计算结果的影响。
r型网格自适应方法作为网格自适应中的一种方法,通过合理移动网格点位置来实现网格自适应分布,在自适应过程中网格单元的连接关系和网格的空间拓扑结构始终保持不变,流场解算软件针对自适应后的网格不需要调整网格数据结构、不必考虑并行动态负载平衡等问题。r型网格自适应能够自动捕捉流场中的复杂结构,减弱网格分布对计算结果的影响。
r型网格自适应中的关键技术是如何计算每个网格点的移动量,实现强间断例如激波附近网格点逐步向流场中的强间断方向移动和聚集。传统r型网格自适应方法以网格点与周围网格点的变量差如压力、密度、温度为权函数,通过加权周围网格点的位置矢量获得网格点位移。这种移动计算方法能够实现网格点向激波移动,但是位于激波同侧与网格点相连的第二层网格点移动量很小,最终结果是只加密了跨过激波的一层网格,激波附近网格没有跟随加密,流场分辨率有限提升。以附图1为例,0、1、2、3、4表示网格单元编号,a、b、c、d、e、f表示网格点编号,此时激波穿过0,4,3三个网格单元。激波前后的变量如密度和压力间断,网格点a与网格点d的压力值差别很大,而与网格点b的压力值差别小,如果只以压力差作为计算网格点移动的权函数,那么网格点a、d向激波移动,但是与网格点a同侧的网格点b则向激波的移动量很小。
基于以上原因,需要探索r型网格自适应的网格点移动位移计算方法,实现激波附近网格点以及相邻几层网格点同步向激波靠近。
发明内容
为了解决上述问题,本发明提出一种基于变量差及相对位移的r型网格自适应移动方法及设备,技术方案如下:
一种基于变量差及相对位移的r型网格自适应移动方法,以归一化网格点的变量差以及网格点相对位移量这两类参数作为计算网格点移动的权函数,通过加入网格点相对位移量这一权函数,实现网格点的跟随移动;归一化消除所述两类参数的量值及量纲对网格点移动量影响,以所述两类参数形成新的组合权值,以此确定激波附近网格点的位移,并将激波附近网格点的位移信息向网格点的相邻网格点传递,实现激波附近的几排网格点都同步向激波位置移动。
进一步的,以压力这一物理量作为特征变量,计算网格中每个网格点的变量值,包括以下步骤:
步骤一:基于距离权计算所有网格点的压力;
步骤二:统计所有网格点的压力最大值和压力最小值;
步骤三:根据所有网格点的压力最大值和压力最小值归一化网格点上的压力;
步骤四:将网格点与相邻网格点的压力之差作为网格点移动位移的第一个权值;
步骤五:将上一迭代步网格点移动位移与网格点之间的初始位移的比值作为网格点移动位移的另一个权值;
步骤六:通过加权当前网格点到周围所有网格点的矢量获得网格点位移。
进一步的,步骤一中,根据格心型二阶有限体积方法,网格单元的变量存储在网格单元体心,将网格点a周围网格单元体心变量通过距离加权到网格点a:
进一步的,步骤二中,所有网格点的压力最大值p max 和压力最小值p min 为:
其中,j为所有网格点,即j=a,b,c,d,e;b,c,d,e为网格点a的相邻网格点。
进一步的,步骤三中,归一化网格点压力的计算方法为:
进一步的,步骤六中,网格点a向相邻网格点i移动的加权系数为:
通过加权网格点a到相邻网格点i的位置矢量差获得网格点a的移动矢量,计算方法如下:
进一步的,上述一种基于变量差及相对位移的r型网格自适应移动方法适用于以密度或熵这两类物理量作为特征变量。
一种基于变量差及相对位移的r型网格自适应移动设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述一种基于变量差及相对位移的r型网格自适应移动方法。
本发明的有益效果在于:
本发明通过归一化物理量以及网格点位置移动量,消除两类参数量值及量纲对网格点移动量影响,可有效实现激波处网格点向激波位置靠近,同时能够保证激波附近的几排网格点都同步向激波位置移动,为r型网格自适应提供了一种新的网格点移动位移方法。
附图说明
图1是网格单元体心及网格点示意图;
图2是初始计算网格;
图3是r型网格自适应10次后的计算网格;
图4是r型网格自适应后的机翼前缘局部网格。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现说明本发明的具体实施方式。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种基于变量差及相对位移的r型网格自适应移动方法,以归一化网格点的变量差以及网格点相对位移量这两类参数作为计算网格点移动的权函数,通过加入网格点相对位移量这一权函数,实现网格点的跟随移动。具体的,归一化消除这两类参数的量值及量纲对网格点移动量影响,以这两类参数形成新的组合权值,以此确定激波附近网格点的位移,并将激波附近网格点的位移信息向网格点的相邻网格点传递,实现激波附近的几排网格点都同步向激波位置移动。
该方法计算网格中每个网格点的变量值,以压力为例,首先将网格点压力归一化处理,将网格点与相邻网格点的压力之差作为网格点移动位移的权值之一,同时将上一迭代步网格点移动位移与网格点之间的初始位移的比值作为网格点移动位移的另一个权值,最后通过加权当前网格点到周围所有网格点的矢量获得网格点位移。该方法可实现网格点在压力剧烈变化区域的聚集加密,同时将网格点移动信息传播给它的相邻网格点,驱动相邻网格点同步向压力剧烈变化区域移动。该方法能够根据流场特征有效移动网格点,为r型网格自适应提供一种新的网格点移动位移计算方法。
在本发明的一个优选实施例中,以压力这一物理量作为特征变量,网格点移动计算方法如下:
步骤一:基于距离权计算所有网格点的压力。
如图1所示,对于格心型二阶有限体积方法,网格单元的变量存储在网格单元0、1、2、3、4的体心,将网格点a周围网格单元体心变量通过距离加权到网格点a:
步骤二:统计所有网格点的压力最大值和压力最小值。
所有网格点的压力最大值p max 和压力最小值p min 为:
其中,j为所有网格点,即j=a,b,c,d,e;b,c,d,e为网格点a的相邻网格点。
步骤三:根据所有网格点的压力最大值和压力最小值归一化网格点上的压力。
由于本方法采用了两种权重系数的和来确定网格点移动量,而计算两种权重系数的变量分别为网格点压力和网格点移动量,这两类变量具有完全不同的量纲和量值,为消除两类变量量值对网格点移动量影响,将网格点上的两类变量进行归一化处理。归一化网格点上的压力的计算方法为:
步骤四:将网格点与相邻网格点的压力之差作为网格点移动位移的第一个权值。
步骤五:将上一迭代步网格点移动位移与网格点之间的初始位移的比值作为网格点移动位移的另一个权值。
为实现高压力梯度附近的网格点a移动信息传递给其相邻网格点b、c、f,保证网格点a相邻网格点b、c、f同步逐步向高压力梯度附近移动,计算网格点a移动位移的另外一个权值,该权值通过网格点a相邻网格点上一迭代步的位移量来计算,计算方法如下:
步骤六:通过加权当前网格点到周围所有网格点的矢量获得网格点位移
网格点a向相邻网格点i移动的加权系数为:
通过加权网格点a到相邻网格点i的位置矢量差获得网格点a的移动矢量,计算方法如下:
至此完成网格点移动位移的计算。
本发明的网格点移动计算方法也适用于基于密度或者熵这两类物理量作为第一个权值的计算变量。
基于上述优选实施例,选定一个几何构型如图2所示,该构型为并列两段NACA0012超音速绕流。其中一个翼型相对另外一个翼型上移和后移0.5米,计算来流马赫数Ma=2.0,来流攻角0°,来流静压P=26500Pa、来流静温T=223.25K。使用NNW-FlowStar软件计算流场,得到流场变量在网格单元上的分布,图2给出了用户初始生成的计算网格,此时网格从翼型表面逐步向空间逐渐变大。
根据步骤一,首先获得图2所示网格的网格点的压力变量值, 然后按照步骤二获取图2所示网格所有网格点的压力最大最小值,并按照步骤三将网格点上的压力归一化,此时为突出网格点物理量压力在网格点移动位移计算中的主导地位,网格点压力影响系数scale取10。
根据步骤四、步骤五的加权系数计算网格点移动位移,自适应调整网格分布后重新数值模拟流场,然后重复步骤一至步骤六,直到达到设定的自适应。图3给出了r型网格自适应迭代10步后的计算网格,从图3看到,依照本方法,可实现网格点在激波位置聚集,并将激波附近的几排网格同步向激波位置移动。图4给出了r型网格自适应迭代10步后的机翼前缘网格,本方法未在机翼前缘这类高压力梯度的光滑极值区域进行网格自适应。本方法网格自适应后数值模拟得到的流场更加清晰、锐利。
本发明还提供了一种基于变量差及相对位移的r型网格自适应移动设备,该设备包括存储器和处理器,存储器存储有计算机程序,处理器执行计算机程序时实现上述基于变量差及相对位移的r型网格自适应移动方法的步骤。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (8)
1.一种基于变量差及相对位移的r型网格自适应移动方法,其特征在于,以归一化网格点的变量差,以及网格点相对位移量这两类参数作为计算网格点移动的权函数,通过加入网格点相对位移量这一权函数,实现网格点的跟随移动;归一化消除所述两类参数的量值及量纲对网格点移动量影响,以所述两类参数形成新的组合权值,以此确定激波附近网格点的位移,并将激波附近网格点的位移信息向网格点的相邻网格点传递,实现激波附近的几排网格点都同步向激波位置移动;
以压力这一物理量作为特征变量,计算网格中每个网格点的变量值,包括以下步骤:
步骤一:基于距离权计算所有网格点的压力;
步骤二:统计所有网格点的压力最大值和压力最小值;
步骤三:根据所有网格点的压力最大值和压力最小值归一化网格点上的压力;
步骤四:将网格点与相邻网格点的压力之差作为网格点移动位移的第一个权值;
步骤五:将上一迭代步网格点移动位移与网格点之间的初始位移的比值作为网格点移动位移的另一个权值;
步骤六:通过加权当前网格点到周围所有网格点的矢量获得网格点位移;
步骤一中,根据格心型二阶有限体积方法,网格单元的变量存储在网格单元体心,将网格点a周围网格单元体心变量通过距离加权到网格点a:
7.根据权利要求1~6任一项所述的一种基于变量差及相对位移的r型网格自适应移动方法,其特征在于,所述r型网格自适应移动方法适用于以密度或熵这两类物理量作为特征变量。
8.一种基于变量差及相对位移的r型网格自适应移动设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1~6任一项所述的一种基于变量差及相对位移的r型网格自适应移动方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110275460.5A CN112665820B (zh) | 2021-03-15 | 2021-03-15 | 基于变量差及相对位移的r型网格自适应移动方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110275460.5A CN112665820B (zh) | 2021-03-15 | 2021-03-15 | 基于变量差及相对位移的r型网格自适应移动方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112665820A CN112665820A (zh) | 2021-04-16 |
CN112665820B true CN112665820B (zh) | 2021-06-04 |
Family
ID=75399418
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110275460.5A Active CN112665820B (zh) | 2021-03-15 | 2021-03-15 | 基于变量差及相对位移的r型网格自适应移动方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112665820B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113673147B (zh) * | 2021-07-23 | 2023-03-17 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种面向流场数据的激波位置智能预测方法 |
CN114398726B (zh) * | 2022-03-21 | 2022-06-07 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于辅助网格的内流道截面流场提取方法、设备及介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1729234A1 (en) * | 2005-05-31 | 2006-12-06 | Lockheed Martin Corporation | Method for designing streamline traced, mixed compression inlets for aircraft engines |
CN105117585A (zh) * | 2015-07-30 | 2015-12-02 | 江苏大学 | 一种两相流泵的自适应网格动态加权误差评估方法 |
CN106960070A (zh) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | 一种基于有限元‑离散元ct重构煤体的渗流模拟方法 |
CN108197367A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
WO2019006268A1 (en) * | 2017-06-30 | 2019-01-03 | Aerion Corporation | CALCULATION TECHNIQUES FOR THREE-DIMENSIONAL MODELING AND DESIGN ANALYSIS |
CN110688791A (zh) * | 2019-08-30 | 2020-01-14 | 中国航天空气动力技术研究院 | 一种钝头体流场激波自适应结构网格生成方法 |
CN110851929A (zh) * | 2019-11-15 | 2020-02-28 | 中国科学院工程热物理研究所 | 基于自适应网格的二维叶型优化设计方法及装置 |
CN111859763A (zh) * | 2020-07-27 | 2020-10-30 | 上海圣之尧智能科技有限公司 | 有限元模拟方法、系统及介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100145668A1 (en) * | 2008-12-04 | 2010-06-10 | Fisher Mark S | Method for dynamic repartitioning in adaptive mesh processing |
CN102436550B (zh) * | 2011-11-07 | 2014-08-13 | 武汉大学 | 复杂边界及实际地形上溃坝洪水的自适应模拟方法 |
EP3213126A1 (en) * | 2014-10-31 | 2017-09-06 | Exxonmobil Upstream Research Company | Handling domain discontinuity in a subsurface grid model with the help of grid optimization techniques |
US11379640B2 (en) * | 2019-01-16 | 2022-07-05 | Saudi Arabian Oil Company | Reservoir regions management with unstructured grid reservoir simuation model |
CN111783343A (zh) * | 2020-07-08 | 2020-10-16 | 中国人民解放军国防科技大学 | 非结构网格方向模板选择方法与流场求解方法 |
-
2021
- 2021-03-15 CN CN202110275460.5A patent/CN112665820B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1729234A1 (en) * | 2005-05-31 | 2006-12-06 | Lockheed Martin Corporation | Method for designing streamline traced, mixed compression inlets for aircraft engines |
CN105117585A (zh) * | 2015-07-30 | 2015-12-02 | 江苏大学 | 一种两相流泵的自适应网格动态加权误差评估方法 |
CN106960070A (zh) * | 2016-12-28 | 2017-07-18 | 山东科技大学 | 一种基于有限元‑离散元ct重构煤体的渗流模拟方法 |
WO2019006268A1 (en) * | 2017-06-30 | 2019-01-03 | Aerion Corporation | CALCULATION TECHNIQUES FOR THREE-DIMENSIONAL MODELING AND DESIGN ANALYSIS |
CN108197367A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
CN110688791A (zh) * | 2019-08-30 | 2020-01-14 | 中国航天空气动力技术研究院 | 一种钝头体流场激波自适应结构网格生成方法 |
CN110851929A (zh) * | 2019-11-15 | 2020-02-28 | 中国科学院工程热物理研究所 | 基于自适应网格的二维叶型优化设计方法及装置 |
CN111859763A (zh) * | 2020-07-27 | 2020-10-30 | 上海圣之尧智能科技有限公司 | 有限元模拟方法、系统及介质 |
Non-Patent Citations (2)
Title |
---|
"High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows";JianhuaPAN 等;《Chinese Journal of Aeronautics》;20181231;第31卷(第9期);第1829-1841页 * |
"动边界非定常流动网格自适应模拟的研究";彭彪;《中国博士学位论文全文数据库 基础科学辑》;20200115(第01期);第A002-248页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112665820A (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112665820B (zh) | 基于变量差及相对位移的r型网格自适应移动方法及设备 | |
Taylor et al. | An algorithm for computing Fekete points in the triangle | |
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
CN113505443B (zh) | 一种任意外形的三维绕流问题自适应笛卡尔网格生成方法 | |
CN108153984B (zh) | 一种基于流场密度阶跃的高精度间断迦辽金人工粘性激波捕捉方法 | |
CN110851929B (zh) | 基于自适应网格的二维叶型优化设计方法及装置 | |
CN102436550B (zh) | 复杂边界及实际地形上溃坝洪水的自适应模拟方法 | |
EP2525296A1 (en) | Three-dimensional fluid simulation method | |
CN108197072B (zh) | 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 | |
CN116245049B (zh) | 节点式非结构网格的边界修正方法、装置、设备及介质 | |
CN111783343A (zh) | 非结构网格方向模板选择方法与流场求解方法 | |
CN112613206B (zh) | 一种基于各向异性体调和场的附面层网格生成方法 | |
CN103729564A (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
CN113850008A (zh) | 飞行器气动特性预测的自适应网格扰动域更新加速方法 | |
CN109726465A (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
CN116306044A (zh) | 全湍流构型的不确定性分析方法及其梯度优化设计方法 | |
CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
Gao et al. | A B-spline method with AIS optimization for 2-D IoT-based overpressure reconstruction | |
CN112632818B (zh) | 高阶梯度光滑无网格方法、系统、介质、计算机设备 | |
CN111159956A (zh) | 一种基于特征的流场间断捕捉方法 | |
CN104299265B (zh) | 一种流体环境影响下的群体行为控制方法 | |
JP2011065360A (ja) | 非分岐かつ非直交の構造格子の境界条件を設定し、反復計算を用いる流れの数値解析方法 | |
Luo | A finite volume method based on weno reconstruction for compressible flows on hybrid grids | |
CN104937640B (zh) | 定义d维空间E中元件的数据处理方法 | |
Ji | Comparison and analysis of different numerical schemes in sod’s one-dimensional shock tube problems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |