一种考虑激波和边界层的网格生成方法技术

本发明专利技术公开了一种考虑激波和边界层的网格生成方法。本发明专利技术通过初始流场捕捉激波位置，然后通过插值光滑获得激波面上的网格，保证激波与流向网格平行，同时在垂直表面方向上通过分段函数针对性布置激波和边界层区域的网格，从而生成能够用于开展扰动过激波的感受性过程研究的计算网格。性过程研究的计算网格。性过程研究的计算网格。

【技术实现步骤摘要】
一种考虑激波和边界层的网格生成方法


[0001]本专利技术涉及网格生成
，具体涉及一种考虑激波和边界层的网格生成方法。

技术介绍

[0002]当空气高速掠过物体或者物体在自然界高速移动时，物体前缘附近会受到强烈压缩而形成冲击波，称为激波。激波后的气体继续流动吹向物体表面，因粘性作用气体附着在物体表面而形成附着薄层，称为边界层。理想情况下，边界层流动呈有规律的层流状态。但在外界干扰下，如空气中引入小扰动，层流逐渐转变为“杂乱无章”的湍流状态，这个过程称为转捩。转捩的发生一般由边界层内扰动失稳导致的，而边界层扰动可由自然空气中的扰动激发产生，这个过程称为感受性。空气中扰动激发边界层内扰动的感受性过程的关键阶段是扰动经过激波和激波后扰动进入边界层这两个阶段，因此如何在激波和边界层区域内安排网格分布，是通过CFD仿真计算方法准确模拟感受性的关键所在。
[0003]目前用于CFD计算的网格生成方法主要来自商业软件，如ICEM CFD、Griden、Pointwise、Gambit等，另外还有流体界学者根据自身研究需求编制的典型外形的网格生成技术。但是它们大多是通用技术，并不是专门针对感受性数值模拟研究的网格生成技术。其技术以及生成的网格存在以下不足：
[0004](1)激波与网格不能保证平行。商业软件不能就已有的流场捕捉激波位置，从而不能根据激波位置布置与激波平行的网格，使得激波与网格相互交叉，在数值计算中容易因网格带来的数值扰动，不利于感受性问题研究。
[0005](2)不能有效布置激波区域的网格。激波区域的网格一方面要平行于激波，另一方面在垂直于激波的网格分布上一阶导数和二阶导数要连续且保证网格数量足够密。这样的安排是为了保证扰动过激波时不会因网格不连续而出现扰动间断，或者因网格稀疏而扰动耗散衰减，导致不能准确捕捉激波后的扰动信息。但是现有网格生成技术还做不到这样的效果。
[0006](3)不能保证边界层内在垂直表面方向上的网格数量。边界层内网格数量至少有100个点时才能保证边界层内流场及扰动的信息准确性。一般网格分布在垂直表面方向上从物体表面到上边界呈单一的指数分布，这样的后果是不能保证激波区域的网格数量。而且在贴近物体表面的网格数量要么过稀不能有效捕捉边界层信息，要么过密导致网格尺度过小而影响CFD计算效率。另外，如果分块单独安排边界层区域的网格，则不能保证边界层内外之间的网格连续性，同样不利于感受性CFD计算。

技术实现思路

[0007]针对现有技术中的上述不足，本专利技术提供的一种考虑激波和边界层的网格生成方法解决了激波与网格不能保证平行、不能有效布置激波区域的网格以及不能保证边界层内在垂直表面方向上的网格数量的问题。
[0008]为了达到上述专利技术目的，本专利技术采用的技术方案为：一种考虑激波和边界层的网格生成方法，其特征在于，包括以下步骤：
[0009]S1、利用商业软件画出钝锥的初始三维计算网格，并利用CFD软件在初始三维计算网格上模拟计算得到初始基本流场；
[0010]S2、根据当地马赫数的大小判断激波位置，当当地马赫数小于0.99倍的来流马赫数时，利用初始三维计算网格和初始基本流场识别激波位置，对激波位置在空间进行光滑处理，获得光滑连续的激波曲面；
[0011]S3、预估边界层的厚度，在边界层以及激波附近区域加密网格，将垂直表面方向的网格分为5段，得到五段组成的网格分布函数y(j)；
[0012]所述垂直表面方向的网格包括第一段网格、第二段网格、第三段网格、第四段网格和第五段网格；
[0013]S4、以激波位置为基准，将垂直表面方向的分布函数y(j)映射到具体每个流向和周向的位置，并通过空间插值方法从初始三维计算网格插值获得考虑激波进和边界层的网格。
[0014]进一步地：所述步骤S3中的第一段网格在边界层区域内以指数规律分布，所述第一段网格的分布函数为：
[0015][0016]上式中，y1(j)为第一段网格的分布函数，y
c1
为边界层区域长度/厚度，p1为第一段网格的指数系数，jn1为第一段网格的边界区域网格点数，jn1>100。
[0017]进一步地：所述第三段网格为边界层与激波之间的区域，网格均匀分布，所述第三段网格的分布函数为y3(j)。
[0018]进一步地：所述第一段网格与第三段网格通过第二段网格以过渡方式连接，所述第二段网格的分布函数为：
[0019][0020]上式中，y2(j)为第二段网格分布函数，jn2为第二段网格的边界区域网格点数，p2为第二段网格的指数系数，p2＝(j-1)/(jn
2-1)。
[0021]进一步地：所述第五段网格为激波外区域并以指数规律分布，所述第三段网格与第五段网格通过第四段网格以过渡方式连接。
[0022]进一步地：所述步骤还包括：利用CFD软件在考虑激波进和边界层的网格上模拟计算得到基本流场。
[0023]本专利技术的有益效果为：本专利技术通过初始流场捕捉激波位置，然后通过插值光滑获得激波面上的网格，保证激波与流向网格平行，同时在垂直表面方向上通过分段函数针对性布置激波和边界层区域的网格，从而生成能够用于开展扰动过激波的感受性过程研究的计算网格。
附图说明
[0024]图1为本专利技术流程图；
[0025]图2为本专利技术中初始计算网格的示意图；
[0026]图3为本专利技术中初始基本流场的示意图；
[0027]图4为本专利技术中网格和激波的分布情况示意图；
[0028]图5为本专利技术中生成的垂直表面方向的网格分布规律示意图；
[0029]图6a为本专利技术中生成的计算网格的正面示意图；
[0030]图6b为本专利技术中生成的计算网格的侧面示意图；
[0031]图7a为本专利技术中计算网格下的基本流场示意图；
[0032]图7b为本专利技术中计算网格下的激波与网格分布情况示意图；
[0033]图7c为本专利技术中计算网格中边界层网格分布情况。
具体实施方式
[0034]下面对本专利技术的具体实施方式进行描述，以便于本
的技术人员理解本专利技术，但应该清楚，本专利技术不限于具体实施方式的范围，对本
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本专利技术的精神和范围内，这些变化是显而易见的，一切利用本专利技术构思的专利技术创造均在保护之列。
[0035]如图1所示，一种考虑激波和边界层的网格生成方法，包括以下步骤：
[0036]S1、利用商业软件画出钝锥的初始三维计算网格，如图2所示，并利用CFD软件在初始三维计算网格上模拟计算得到初始基本流场，如图3所示；图4给出初始三维计算网格和基本流场在激波附近的局部放大图，可以看出激波附近网格与激波不平行，且附近的流场出现毛刺、不光滑的现象，激波后的流场出现抖动。
[0037]S2、根据当地马赫数的大小判断激波位置，当当地马赫数小于0.99倍的来流马赫数时，利用初始三维计算网格本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种考虑激波和边界层的网格生成方法，其特征在于，包括以下步骤：S1、利用商业软件画出钝锥的初始三维计算网格，并利用CFD软件在初始三维计算网格上模拟计算得到初始基本流场；S2、根据当地马赫数的大小判断激波位置，当当地马赫数小于0.99倍的来流马赫数时，利用初始三维计算网格和初始基本流场识别激波位置，对激波位置在空间进行光滑处理，获得光滑连续的激波曲面；S3、预估边界层的厚度，在边界层以及激波附近区域加密网格，将垂直表面方向的网格分为5段，得到五段组成的网格分布函数y(j)；所述垂直表面方向的网格包括第一段网格、第二段网格、第三段网格、第四段网格和第五段网格；S4、以激波位置为基准，将垂直表面方向的分布函数y(j)映射到具体每个流向和周向的位置，并通过空间插值方法从初始三维计算网格插值获得考虑激波和边界层的网格。2.根据权利要求1所述的考虑激波和边界层的网格生成方法，其特征在于，所述步骤S3中的第一段网格在边界层区域内以指数规律分布，所述第一段网格的分布函数为：上式中，y1(j)为第一段网格的分布函数，y
...

【专利技术属性】
技术研发人员：万兵兵，李晓虎，陈坚强，袁先旭，涂国华，段茂昌，
申请(专利权)人：空气动力学国家重点实验室，
类型：发明
国别省市：