基于一般时间根方尺度的雷诺应力湍流模型的数值方法

1.本技术涉及计算流体力学领域，特别是涉及一种基于一般时间根方尺度的雷诺应力湍流模型的数值方法、装置、计算机设备和存储介质。


背景技术：

2.随着计算技术的发展，计算流体力学(cfd)已经成为研究湍流结构生成演化机理，解决实际工程湍流问题的三大手段(理论、实验、计算)之一。尤其是cfd软件的蓬勃发展，为设计师和科研人员带来了极大的便利。湍流的数值模拟方法传统意义上包括三大类：雷诺平均模拟(rans)方法、大涡模拟(les)方法和直接数值模拟(dns)方法。由于计算机硬件条件的限制，rans方法仍然为解决工程问题的主要方法。
3.rans方法可追溯到一百多年前，boussinesq提出著名的涡粘性假设来模拟雷诺应力。其后，诸多流体力学先驱发展了一系列半经验理论，但由于只考虑了一阶湍流统计量，这类型湍流模型也被称为一阶矩模型。物理上更为完善的封闭建模方法，通常称之为高阶矩模型，如雷诺应力模型(rsm)。1940年，我国著名力学大师周培源建立了一般湍流的雷诺应力所满足的输运微分方程组，并对其中出现的三元速度关联等新未知量提出了封闭性假设。1951年，rotta进一步发展了周培源的工作，并提出了完整的rsm模型。其后，donaldson等提出了模型不变性概念，lumley发展的封闭性近似等，但都尚处于初级研究阶段。
4.到了上世纪90年代，计算机水平的飞速发展和一阶矩模型的不足使得人们又开始重新重视二阶矩模型。其中，speziale认为rsm给出了用于复杂湍流流动计算的les方法和rans方法之间的桥梁。他进一步讨论了雷诺应力模型的可实现性问题，并提出了简化的设计准则。其后，诸多学者基于上述准则，发展了多种rsm模型，如lrr-ε模型、lrr-ω模型、ssg-ε模型、ssg/lrr-ω模型等。这些rsm模型被众多商业cfd软件集成，并开展了广泛的应用。
5.在飞行器湍流流场数值模拟的工程应用中，通常使用上述rsm模型，但上述rsm模型多基于ε尺度或ω尺度，在黏性壁面处不具有自然的边界条件，会在采用高精度离散或复杂结构/非结构网格时导致一些数值不稳定。因此，现有技术存在适应性不佳的问题。


技术实现要素：

6.基于此，有必要针对上述技术问题，提供一种能够解决飞行器流场模拟采用高精度离散或复杂结构/非结构网格时数值不稳定问题的基于一般时间根方尺度的雷诺应力湍流模型的数值方法、装置、计算机设备和存储介质。
7.一种基于一般时间根方尺度的雷诺应力湍流模型的数值方法，所述方法包括：
8.获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程；ω表示比耗散率尺度；
9.获取预先推导的一般时间根方尺度和ω尺度的关系式，根据所述关系式和所述ω尺度方程得到湍流关于尺度的尺度方程；其中，τ表示含能时间尺度，n表示n次根方，为正整数；
10.获取湍流基于ssg/lrr-ω模型的方程；表示雷诺应力方程分量；
11.将所述方程的相关项修改与所述尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型；
12.获取雷诺平均navier-stokes(rans)方程组，与所述基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组；
13.构建待模拟的飞行器流场的网格数据，根据所述网格数据对所述耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
14.在其中一个实施例中，还包括：获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程为：
[0015][0016]
其中，为时间平均密度，t为时间，i,j,k做下标时为坐标索引，为favre平均速度三分量，xj(j＝1,2,3)为三方向坐标分量，α
ω
为比耗散率尺度方程生成项系数，p
kk
/2＝(p
11
p
22
p
33
)/2为湍动能生成，其中p
11
、p
22
、p
33
为三个方向雷诺正应力的生成，β
ω
为比耗散率尺度方程耗散项系数，μ为动力学黏性系数，σ
ω
为比耗散率尺度方程扩散项系数，σd为比耗散率尺度方程交叉导数项系数，为湍动能，为三个方向雷诺正应力。
[0017]
在其中一个实施例中，还包括：所述预先推导的一般时间根方尺度和ω尺度的关系式为：
[0018][0019]
其中，n为时间根方尺度和ω尺度的调节系数。
[0020]
在其中一个实施例中，还包括：获取预先推导的时间根方尺度和ω尺度的关系式，根据所述关系式和所述ω尺度方程得到湍流关于尺度的尺度方程为：
[0021][0022]
其中，表示一般时间根方尺度。
[0023]
在其中一个实施例中，还包括：获取湍流ssg/lrr-ω模型的方程为：
[0024]
[0025]
其中，为生成项，为再分配项，为耗散项，为扩散项。
[0026]
在其中一个实施例中，还包括：
[0027]
所述生成项为：
[0028][0029]
所述再分配项为：
[0030][0031]
其中，为各向异性张量，为应变量，为特殊应变量，c1、c2、c3、c4、c5为预知的ssg/lrr雷诺应力方程中的系数，为涡量。
[0032]
所述耗散项为：
[0033][0034]
其中，各向同性耗散率c
μ
为耗散系数。
[0035]
所述扩散项为：
[0036][0037]
其中，为涡黏性系数。
[0038]
在其中一个实施例中，还包括：将所述方程的相关项修改与所述尺度方程进行耦合，基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型为：
[0039][0040]
在其中一个实施例中，还包括：获取雷诺平均navier-stokes(rans)方程组为：
[0041][0042]
其中，为时间平均压力，为favre平均温度，为favre平均总能，c
p
为定压比热比，pr为层流普朗特常数，pr
t
为湍流普朗特常数；
[0043]
其中为黏性应力张量：
[0044][0045]
τ
ij
为雷诺应力张量，与的关系为：
[0046][0047]
从而完成rans方程组与ssg/lrr-模型耦合。
[0048]
在其中一个实施例中，还包括：根据所述网格数据，利用针对偏微分方程组的数值法对所述耦合方程组中包含的12个独立变量：：和λ进行数值求解；再通过关系式推导得到其他变量的数值解。
[0049]
一种基于一般时间根方尺度的雷诺应力湍流模型的数值装置，所述装置包括：
[0050]
ω尺度方程获取模块，用于获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程；ω表示比耗散率尺度；
[0051]
一般时间根方尺度方程确定模块，用于获取预先推导的一般时间根方尺度和ω尺度的关系式，根据所述关系式和所述ω尺度方程得到湍流关于尺度方程；其中，τ表示含能时间尺度，n表示n次根方，为正整数；
[0052]
基于一般时间根方尺度的雷诺应力湍流模型确定模块，用于获取湍流基于ssg/lrr-ω模型的方程；表示雷诺应力方程分量；将所述方程的相关项修改与所述尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型；获取雷诺平均navier-stokes(rans)方程组，与所述基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组；
[0053]
数值模拟模块，用于构建待模拟的飞行器(或其他空气动力学问题)流场的网格数据，根据所述网格数据对所述耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
[0054]
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
[0055]
获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程；ω表示比耗散率尺度；
[0056]
获取预先推导的时间根方尺度和ω尺度的关系式，根据所述关系式和所述ω尺度方程得到湍流关于尺度的尺度方程；其中，τ表示含能时间尺度，n表示n次根方，为正整数；
[0057]
获取湍流基于ssg/lrr-ω模型的方程；表示雷诺应力方程分量；
[0058]
将所述方程的相关项修改与所述尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型；
[0059]
获取雷诺平均navier-stokes(rans)方程组，与所述基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组；
[0060]
构建待模拟的飞行器(或其他空气动力学问题)流场的网格数据，根据所述网格数据对所述耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
[0061]
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
[0062]
获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程；ω表示比耗散率尺度；
[0063]
获取预先推导的一般时间根方尺度和ω尺度的关系式，根据所述关系式和所述ω尺度方程得到湍流关于尺度的尺度方程；其中，τ表示含能时间尺度，n表示n次根方，为正整数；
[0064]
获取湍流基于ssg/lrr-ω模型的方程；表示雷诺应力方程分量；
[0065]
将所述方程的相关项修改与所述尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型；
[0066]
获取雷诺平均navier-stokes(rans)方程组，与所述基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组；
[0067]
构建待模拟的飞行器流场的网格数据，根据所述网格数据对所述耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
[0068]
上述基于一般时间根方尺度的雷诺应力湍流模型的数值方法、装置、计算机设备和存储介质，基于现有的ssg/lrr-ω雷诺应力模型的ω尺度方程，根据推导得到的一般时间根方尺度和ω尺度的关系式，将两者结合得到湍流关于尺度的尺度方程；将尺度方程和预知的湍流基于ssg/lrr-ω模型的方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型，获取雷诺平均navier-stokes(rans)方程组，与基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组，在进行飞行器流场数值模拟时，构建待模拟的飞行器湍流流场的网格数据，根据网格数据对耦合方程组进
行数值求解，得到飞行器湍流流场的数值模拟结果。本发明所提出的ssg/lrr-模型中，n为正整数时壁面边界条件严格为0，有利于降低方程对壁面信息的依赖程度；进一步n为大于2的正偶数时，λ的计算值将不会影响雷诺应力方程耗散项的符号，这对于rsm的得到满足物理的解十分有利，能够在采用高精度离散或复杂结构/非结构网格时实现数值稳定。
附图说明
[0069]
图1为一个实施例中基于一般时间根方尺度的雷诺应力湍流模型的数值方法的流程示意图；
[0070]
图2为一个具体实施例中尖前缘三角翼绕流：表面网格与站位划分示意图，其中(a)为表面网络示意图，(b)为占位划分示意图；
[0071]
图3为一个具体实施例中尖前缘三角翼绕流中ssg/lrr-模型参数n不同取值的收敛结果程示意图；
[0072]
图4为一个具体实施例中尖前缘三角翼绕流中ssg/lrr-模型在不同站位得到的压力分布示意图，其中，(a)为x/cr＝0.40得到的压力分布示意图，(b)为x/cr＝0.60得到的压力分布示意图，(c)为x/cr＝0.80得到的压力分布示意图，(d)为x/cr＝0.95得到的压力分布示意图；
[0073]
图5为一个实施例中基于一般时间根方尺度的雷诺应力湍流模型的数值装置的结构框图；
[0074]
图6为一个实施例中计算机设备的内部结构图。
具体实施方式
[0075]
为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
[0076]
在一个实施例中，如图1所示，提供了一种基于一般时间根方尺度的雷诺应力湍流模型的数值方法，包括以下步骤：
[0077]
步骤102，获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程。
[0078]
ω表示比耗散率尺度。
[0079]
飞行器湍流流场的数值模拟基于所提出的物理模型，现有的多种rsm模型多是基于ε尺度或ω尺度，在黏性壁面处不具有自然的边界条件，会在采用高精度离散或复杂结构/非结构网格时导致数值不稳定，对飞行器流场的数值模拟产生不利影响。本发明提出一种在黏性壁面具备自然边界的尺度方程，基于一般时间根方尺度并耦合雷诺应力六方程，形成ssg/lrr-雷诺应力模型。
[0080]
具体地，获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程为：
[0081][0082]
其中，为时间平均密度，t为时间，i,j,k做下标时为坐标索引，为
favre平均速度三分量，xj(j＝1,2,3)为三方向坐标分量，α
ω
为比耗散率尺度方程生成项系数，p
kk
/2＝(p
11
p
22
p
33
)/2为湍动能生成，其中p
11
、p
22
、p
33
为三个方向雷诺正应力的生成，β
ω
为比耗散率尺度方程耗散项系数，μ为动力学黏性系数，σ
ω
为比耗散率尺度方程扩散项系数，σd为比耗散率尺度方程交叉导数项系数，为湍动能，为三个方向雷诺正应力。
[0083]
步骤104，获取预先推导的一般时间根方尺度和ω尺度的关系式，根据关系式和ω尺度方程得到湍流关于尺度的尺度方程。
[0084]
其中，τ表示含能时间尺度，n表示n次根方，为正整数。
[0085]
本发明通过推导，得出时间根方尺度和ω尺度的关系式为：
[0086][0087]
结合ω尺度方程得到关于尺度方程：
[0088][0089]
其中，表示一般时间根方尺度。
[0090]
其中n为正整数，用于调节该尺度方程的数值特性。n＝1时是一种特殊情况，即还原回spezia1e在1992年提出的时间尺度方程。该方程的右端第二项为常数，即无生成项，会造成严重的数值刚性问题。
[0091]
步骤106，获取湍流基于ssg/lrr-ω模型的方程。
[0092]
表示雷诺应力方程分量。
[0093]
具体地，方程为：
[0094][0095]
其中，为生成项，为再分配项，为耗散项，为扩散项。
[0096]
生成项为：
[0097][0098]
再分配项为：
[0099][0100]
其中，为各向异性张量，为应变量，为特殊应变量，c1、c2、c3、c4、c5为预知的ssg/lrr雷诺应力方程中的系数，为涡量。
[0101]
耗散项为：
[0102][0103]
其中，c
μ
为耗散系数。
[0104]
扩散项为：
[0105][0106]
其中，为涡黏性系数。
[0107]
上述雷诺应力方程中的系数通过过渡函数f1加权获得：
[0108]
φ＝f1φ
(lrr)
(1-f1)φ
(ssg)
[0109]
同样，尺度方程中的系数也通过过渡函数f1加权获得：
[0110]
φ＝f1φ
(ω)
(1-f1)φ
(ε)
[0111]
其中φ代表任意系，具体值分别参见表1和表2。
[0112]
表1雷诺应力方程中的系数
[0113][0114]
表2一般时间根方尺度方程中的系数
[0115][0116]
其余系数还有：
[0117][0118]cμ
＝0.09
[0119]
过渡函数f1为：
[0120][0121]
其中：
[0122]
arg1＝min[max(term1,term2),term3]
[0123][0124][0125]
最后，需要强调的是雷诺应力六分量与尺度变量在黏性壁面处均取为0。
[0126]
步骤108，将方程的相关项修改与尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型。
[0127]
基于一般时间根方尺度的雷诺应力湍流模型即ssg/lrr-模型。具体为：
[0128][0129]
步骤110，获取雷诺平均navier-stokes(rans)方程组，与基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组。
[0130]
步骤112，构建待模拟的飞行器流场的网格数据，根据网格数据对耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
[0131]
根据网格数据，利用针对偏微分方程组的数值法对耦合方程组中包含的12个独立变量：和λ进行数值求解；
[0132]
再通过关系式推导得到其他变量的数值解。
[0133]
在构建好基于一般时间根方尺度的雷诺应力湍流模型后，根据所要模拟的飞行器流场，建立对应的网格数据，网格数据包括网格的形状以及网格节点的位置。根据网格数据对基于一般时间根方尺度的雷诺应力湍流模型进行数值求解，得到飞行器湍流流场的分析结果。
[0134]
上述基于一般时间根方尺度的雷诺应力湍流模型的数值方法中，基于现有的ssg/lrr-ω雷诺应力模型的ω尺度方程，根据推导得到的一般时间根方尺度和ω尺度的关系式，将两者结合得到湍流关于尺度的尺度方程；将尺度方程和预知的湍流基于ssg/lrr-ω模型的方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型，获取雷诺平均navier-stokes(rans)方程组，与基于一般时间根方尺度
的雷诺应力湍流模型耦合，得到耦合方程组，在进行飞行器流场数值模拟时，构建待模拟的飞行器湍流流场的网格数据，根据网格数据对耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。本发明所提出的ssg/lrr-模型中，n为正整数时壁面边界条件严格为0，有利于降低方程对壁面信息的依赖程度；进一步n为大于2的正偶数时，λ的计算值将不会影响雷诺应力方程耗散项的符号，这对于rsm的得到满足物理的解十分有利，能够在采用高精度离散或复杂结构/非结构网格时实现数值稳定。
[0135]
应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0136]
在一个具体实施例中，nasa langley中心利用国家跨声速设施(national transonic facility，ntf)对65
°
后掠三角翼开展了多因素的风洞实验研究，具体考察了马赫数、雷诺数、攻角和前缘钝度等对于流动的影响。随着自由流马赫数达到跨声速，三角翼前缘分离涡形态逐渐发生变化，特别是激波使得流场变得极其复杂，出现了许多与亚声速流场不同的现象。因此，第二届国际涡流动实验(the second international vortex flow experiment，vfe)就专门对跨声速条件下三角翼上表面出现的激波/涡干扰和涡破裂现象开展了一系列研究。
[0137]
本实施例利用本发明所提方法对尖前缘三角翼外形在雷诺数rec＝6
×
106，马赫数ma
ref
＝0.85和攻角α＝22.6
°
的条件下开展数值模拟研究。图2给出了所使用的多块结构网格示意图，其网格量约为235万。
[0138]
图3给出了尖前缘三角翼绕流中参数n采用不同取值时的计算收敛结果。对于ssg/lrr-模型，只有n＝2时高阶计算才会发散，而采用其它三个值则可对二阶muscl格式、五阶wcns-e5格式、七阶wncs-e7格式和九阶wcns-e9格式等都实现计算收敛。
[0139]
由于本算例的计算攻角小于涡破裂临界攻角，因而流动仍可看作存在复杂涡系结构的定常流动。具体地，当流动经过前缘时，三角翼背风面会形成一个主体分离涡和一个二次分离涡。这两者可通过壁面压力分布中的吸力峰来分辨。图4展示了ssg/lrr-模型计算的壁面压力分布，其中n＝8。通过与实验值对比对比发现，在该网格上使用二阶muscl格式不能清楚地解析二次涡的结构。相反地，高阶wcns格式的结果与实验数据吻合较好。
[0140]
在一个实施例中，如图5所示，提供了一种基于一般时间根方尺度的雷诺应力湍流模型的数值装置，包括：ω尺度方程获取模块502、一般时间根方尺度方程确定模块504、基于一般时间根方尺度的雷诺应力湍流模型确定模块506和数值模拟模块508，其中：
[0141]
ω尺度方程获取模块502，用于获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程；ω表示比耗散率尺度；
[0142]
一般时间根方尺度方程确定模块504，用于获取预先推导的一般时间根方尺度和ω尺度的关系式，根据关系式和ω尺度方程得到湍流关于尺度的尺度方程；其中，
τ表示含能时间尺度，n表示n次根方，为正整数；
[0143]
基于一般时间根方尺度的雷诺应力湍流模型确定模块506，用于获取湍流基于ssg/lrr-ω模型的方程；表示雷诺应力方程分量；将所述方程的相关项修改与所述尺度方程进行耦合，得到基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型；获取雷诺平均navier-stokes(rans)方程组，与所述基于一般时间根方尺度的雷诺应力湍流模型耦合，得到耦合方程组；
[0144]
数值模拟模块508，用于构建待模拟的飞行器(或其他空气动力学问题)流场的网格数据，根据所述网格数据对所述耦合方程组进行数值求解，得到飞行器湍流流场的数值模拟结果。
[0145]
ω尺度方程获取模块502还用于获取湍流基于ssg/lrr-ω雷诺应力模型的ω尺度方程为：
[0146][0147]
其中，为时间平均密度，t为时间，i,j,k做下标时为坐标索引，为favre平均速度三分量，xj(j＝1,2,3)为三方向坐标分量，α
ω
为比耗散率尺度方程生成项系数，p
kk
/2＝(p
11
p
22
p
33
)/2为湍动能生成，其中p
11
、p
22
、p
33
为三个方向雷诺正应力的生成，β
ω
为比耗散率尺度方程耗散项系数，μ为动力学黏性系数，σ
ω
为比耗散率尺度方程扩散项系数，σd为比耗散率尺度方程交叉导数项系数，为湍动能，为三个方向雷诺正应力。
[0148]
一般时间根方尺度方程确定模块504还用于获取预先推导的一般时间根方尺度和ω尺度的关系式，根据关系式和ω尺度方程得到湍流关于尺度的尺度方程为：
[0149][0150]
其中，表示一般时间根方尺度。
[0151]
基于一般时间根方尺度的雷诺应力湍流模型确定模块506还用于获取湍流ssg/lrr-ω模型的方程为：
[0152][0153]
其中，为生成项，为再分配项，为耗散项，为扩散项。
[0154]
基于一般时间根方尺度的雷诺应力湍流模型确定模块506还用于将方程的相关
项修改与尺度方程进行耦合，基于一般时间根方尺度的雷诺应力湍流模型，即ssg/lrr-模型为：
[0155][0156]
基于一般时间根方尺度的雷诺应力湍流模型确定模块506还用于获取雷诺平均navier-stokes(rans)方程组：
[0157][0158]
其中，为时间平均压力，为favre平均温度，为favre平均总能，c
p
为定压比热比，pr为层流普朗特常数，pr
t
为湍流普朗特常数；
[0159]
其中为黏性应力张量：
[0160][0161]
τ
ij
为雷诺应力张量，与的关系为：
[0162][0163]
从而完成rans方程组与ssg/lrr-模型耦合。
[0164]
数值模拟模块508还用于根据所述网格数据，利用针对偏微分方程组的数值法对所述耦合方程组中包含的12个独立变量：：和λ进行数值求解；再通过关系式推导得到其他变量的数值解。
[0165]
关于基于一般时间根方尺度的雷诺应力湍流模型的数值装置的具体限定可以参见上文中对于基于一般时间根方尺度的雷诺应力湍流模型的数值方法的限定，在此不再赘述。上述基于一般时间根方尺度的雷诺应力湍流模型的数值装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执
行以上各个模块对应的操作。
[0166]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图6所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于一般时间根方尺度的雷诺应力湍流模型的数值方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0167]
本领域技术人员可以理解，图6中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0168]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述方法实施例中的步骤。
[0169]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例中的步骤。
[0170]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
[0171]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0172]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。