一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法技术

本发明专利技术公开了一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法，包括如下步骤：S1，计算某一位置的数值通量时，对于所给的最高精度WENO格式，根据一个准则判断是否需要降阶；S2，如果需要降阶，则对于降阶后的次高精度WENO格式，根据步骤S1中的同一准则判断是否需要降阶；S3，循环执行步骤S2，直到不需要降阶或者到达最低阶精度的WENO格式；然后采用相应精度的WENO格式计算该位置的数值通量，以进一步计算得到该位置对含激波的可压缩流动问题的数值模拟值。本发明专利技术可用于守恒型差分格式，并且由于预先确定所用格式，因此计算量较小。

【技术实现步骤摘要】
一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法
本专利技术涉及可压缩流动的数值模拟
，尤其是一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法。
技术介绍
在可压缩流动的数值模拟中，可能同时存在间断和其他连续的小尺度结构，高阶精度WENO(Weightedessentiallynon-oscillatory，加权本质无振荡)格式既可以捕捉间断，又能在光滑区域达到高阶精度，常用来模拟此类问题。通常更高阶的格式能在流场连续区域获得更精细的流动结构，并且能更加锐利地捕捉激波。但是，对于强非线性间断问题，采用高阶WENO格式更容易出现非物理震荡，甚至可能计算失败。原因之一是高阶格式的耗散更低，因而可能无法有效地抑制振荡。其二是它们所用的较宽模板可能包含多于一个的间断。那么，所有子模板可能都不光滑，因此无法得到本质无振荡的解。其三是高阶插值或重构中的Runge现象。偏离中心的宽模板上的数值通量可能导致较大的数值误差，进而影响稳定性。对于欧拉方程，由于不同特征方向的相互作用，情况会变得更糟。非常高阶格式的解中可能会出现严重的数值振荡，并且计算更可能会发散。对于非常高阶WENO格式，现有方法采用递归降阶来解决这个问题，该方法需要密度和压强在单元边界处的重构值。然而，高阶WENO差分格式的重构通常是在通量上进行的。因此，该方法不能直接应用到守恒型差分格式中，而需要构造其它方法。
技术实现思路
本专利技术所要解决的技术问题是：针对上述存在的问题，提供一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法。本专利技术采用的技术方案如下：一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法，包括如下步骤：S1，计算某一位置的数值通量时，对于所给的最高精度WENO格式，根据一个准则判断是否需要降阶；S2，如果需要降阶，则对于降阶后的次高精度WENO格式，根据步骤S1中的同一准则判断是否需要降阶；S3，循环执行步骤S2，直到不需要降阶或者到达最低阶精度的WENO格式；然后采用相应精度的WENO格式计算该位置的数值通量，以进一步计算得到该位置对含激波的可压缩流动问题的数值模拟值。进一步的，记所求数值通量的位置为xj+1/2，WENO格式用到的模板为S＝{xj-m+1,xj-m+2,…,xj+m}，则步骤S1与步骤S2中的所述准则为同时满足如下两个要求：(1)当模板点数2m≥4时，压强的光滑因子需满足：β(2m)≤4β(2m-2)；光滑因子β(2m)的计算需用到模板S上所有2m个点；光滑因子β(2m-2)的计算则用到模板S去掉两个端点后的所有2m-2个点；(2)压强p和速度u需满足：i∈(j-m+1,…,j+m-1),l∈(j-m+1,…,j+m-2)；其中，γ为气体比热比；h及Δt分别为当前方向空间网格尺度和时间步长；d为所求问题的空间维数。进一步的，所述光滑因子β(2m)或光滑因子β(2m-2)的计算公式如下：其中，β(n)表示光滑因子β(2m)或光滑因子β(2m-2)；xc为模板的中点，即xj+1/2；σ,δ的定义如下：其中，h为当前方向空间网格尺度。进一步的，所述光滑因子β(2m)或光滑因子β(2m-2)的计算公式如下：β(n)＝(bn)2+|an·cn|其中，β(n)表示光滑因子β(2m)或光滑因子β(2m-2)；xc为模板的中点，即xj+1/2。进一步的，an,bn,cn均可写为如下形式：其中，k为光滑因子所用模板最左端点的下标，z为an,bn或cn；系数αn,l的取值如下：进一步的，n＞3，且n为整数。综上所述，由于采用了上述技术方案，本专利技术的有益效果是：1、本专利技术可以在连续区域采用非常高阶的WENO格式，从而获得更精细的流场结构；2、本专利技术在强激波附近会启用相对低阶的WENO格式，因此具有良好的激波稳定性，即激波附近出现的非物理振荡的可能性更小，以及计算中更不易由于出现负压而导致计算失败；3、本专利技术可用于守恒型差分格式，并且由于预先确定所用格式，因此计算量较小。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本专利技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。图1为本专利技术通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法的流程框图。图2为本专利技术算例一中网格点数为N＝200，计算终止时间T＝1.8时的几种格式计算结果的对比图。图3为本专利技术算例二中网格点数为N＝400，计算终止时间T＝0.038时的几种格式计算结果的对比图。图4a为本专利技术算例三中计算区域是[0,1]×[0,1]，网格点数为N＝400，计算终止时间T＝0.8时，WENO5-M格式的计算结果的密度云图。图4b为本专利技术算例三中计算区域是[0,1]×[0,1]，网格点数为N＝400，计算终止时间T＝0.8时，PORWENO9-S格式的计算结果的密度云图。图4c为本专利技术算例三中计算区域是[0,1]×[0,1]，网格点数为N＝400，计算终止时间T＝0.8时，PORWENO17-S格式的计算结果的密度云图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术，即所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本专利技术实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本专利技术的实施例的详细描述并非旨在限制要求保护的本专利技术的范围，而是仅仅表示本专利技术的选定实施例。基于本专利技术的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。说明：对于欧拉方程的守恒型差分方法求解，通常的数值求解步骤如下：步骤1，网格划分，根据计算要求和计算条件生成网格；步骤2，赋初值，根据初值条件在网格点上给出节点初值；步骤3，设置时间步长，根据问题和所使用方法结合CFL条件得到；步骤4，开始循环如下时间推进过程，直到问题所要求的终止时间：(1)对于内点，采用空间格式进行数值模拟计算；(2)根据边界条件，采用合适的数值边界处理；(3)内点和边界的数值计算空间离散完成后，采用合适的时间格式进行推进。本专利技术主要涉及步骤4中第(1)分步骤采用空间格式进行数值模拟计算的方法，由此本专利技术公开了一种通过高阶WENO格式降阶对可压缩流动问本文档来自技高网...

【技术保护点】
1.一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法，其特征在于，包括如下步骤：/nS1，计算某一位置的数值通量时，对于所给的最高精度WENO格式，根据一个准则判断是否需要降阶；/nS2，如果需要降阶，则对于降阶后的次高精度WENO格式，根据步骤S1中的同一准则判断是否需要降阶；/nS3，循环执行步骤S2，直到不需要降阶或者到达最低阶精度的WENO格式；然后采用相应精度的WENO格式计算该位置的数值通量，以进一步计算得到该位置对含激波的可压缩流动问题的数值模拟值。/n

【技术特征摘要】
1.一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法，其特征在于，包括如下步骤：
S1，计算某一位置的数值通量时，对于所给的最高精度WENO格式，根据一个准则判断是否需要降阶；
S2，如果需要降阶，则对于降阶后的次高精度WENO格式，根据步骤S1中的同一准则判断是否需要降阶；
S3，循环执行步骤S2，直到不需要降阶或者到达最低阶精度的WENO格式；然后采用相应精度的WENO格式计算该位置的数值通量，以进一步计算得到该位置对含激波的可压缩流动问题的数值模拟值。


2.根据权利要求1所述的一种通过高阶WENO格式降阶对可压缩流动问题进行数值模拟的方法，其特征在于，记所求数值通量的位置为xj+1/2，WENO格式用到的模板为S＝{xj-m+1，xj-m+2，…，xj+m}，则步骤S1与步骤S2中的所述准则为同时满足如下两个要求：
(1)当模板点数2m≥4时，压强的光滑因子需满足：
β(2m)≤4β(2m-2)；
光滑因子β(2m)的计算需用到模板S上所有2m个点；
光滑因子β(2m-2)的计算则用到模板S去掉两个端点后的所有2m-2个点；
(2)压强p和速度u需满足：



其中，γ为气体比热比；h及Δt分别为当前方向空间网格尺度和时间步长；d为所求问题的空间维数。


...

【专利技术属性】
技术研发人员：武从海，罗勇，韩帅斌，田浩，张树海，李虎，王益民，马瑞轩，
申请(专利权)人：空气动力学国家重点实验室，
类型：发明
国别省市：四川;51