基于一般时间根方尺度的两方程涡黏湍流模型的数值方法

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


背景技术：

2.随着计算技术的发展，计算流体力学(cfd)已经成为研究湍流结构生成演化机理，解决实际工程湍流问题的三大手段(理论、实验、计算)之一。尤其是cfd软件的蓬勃发展，为设计师和科研人员带来了极大的便利。湍流的数值模拟方法传统意义上包括三大类：雷诺平均模拟(rans)方法、大涡模拟(les)方法和直接数值模拟(dns)方法。由于计算机硬件条件的限制，rans方法仍然解决工程问题的主要方法。
3.涡黏性湍流模型是rans方法中最主要的一类湍流模型，基于boussinesq假设：
[0004][0005]
涡黏性系数μ
t
＝ρν
t
是boussinesq近似的比例因子。其中，两方程湍流模型是过去五十年来湍流模式理论研究的基础。这类模型大多求解湍动能k和湍流长度尺度l的等效值。
[0006]
湍流尺度提供变量主要有：动能耗散率类、含能长度尺度类和含能时间尺度类。现在的cfd商业软件中普遍使用第一类。其中各向同性耗散率ε方程起源于对雷诺应力方程或湍动能方程耗散项的直接封闭，是由周培源先生首先得到。
[0007]
经过几十年的发展，ε模型演化出了众多版本，并被广泛使用。
[0008][0009]
但ε作为特征尺度存在两个问题：一是耗散率缺乏自然边界条件；二是在壁面附近需要处理高阶关联。这两个问题对于模型的稳定求解和壁面附近的模拟精度产生了不利的影响。针对这两个问题,研究者们采用了不同的改进方式,比如比耗散率ω：
[0010][0011]
比耗散率的概念由kolmogorov在1942年就提出了，但真正获得成功是通过wilcox、menter等人的工作。ω尺度方程虽然改善了ε方程处理高阶关联引起的数值刚性问题，但由于仍缺乏自然边界条件，在对飞行器或其他空气动力学问题进行湍流流场数值模拟的工程应用中，会在采用高精度离散或复杂结构/非结构网格时导致一些数值不稳定。因此，现有技术存在适应性不佳的问题。


技术实现要素：

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