CN107273593A - 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 - Google Patents
一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 Download PDFInfo
- Publication number
- CN107273593A CN107273593A CN201710402397.0A CN201710402397A CN107273593A CN 107273593 A CN107273593 A CN 107273593A CN 201710402397 A CN201710402397 A CN 201710402397A CN 107273593 A CN107273593 A CN 107273593A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mover
- msup
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法,属于高超声速飞行器热防护系统设计领域。本发明首先采用无量纲压强计算网格节点(I,J,K)的光滑因子;并计算该网格节点的衰减函数值,确定流场中强间断的区域;衰减函数耦合k‑ωSST湍流模型,建立了一种用于高马赫数强激波流场气动热预测的湍流模型。与现有的激波间断检测方法和湍流模型相比,本发明所提出的基于光滑因子概念的强激波间断检测方法对于复杂外形钝头体飞行器仍能够实现自动检测;所提出的衰减函数结合k‑ωSST湍流模型实现高精度气动热预测模拟,计算精度显著提高,误差能够降低至10%以内;所提出的模型和方法实用性强;易于融入现代并行化CFD计算程序。
Description
技术领域
本发明属于高超声速飞行器热防护系统设计领域,具体涉及一种可用于高马赫数强激波流场气动热预测的湍流模型。
背景技术
临近空间高超声速飞行器已成为世界强国航空航天技术发展计划的主要目标。此类飞行器由于需要进行长时间的高超声速飞行将产生严重的气动加热问题,因此,飞行器热防护系统的设计成为高超声速飞行器发展的一个关键技术难点。通过准确预测飞行器的气动热环境不仅可以显著提高热防护系统的有效性,还能够在一定程度上增加飞行器的有效载荷。
当高超声速飞行器在较稠密大气层内飞行时,绕飞行器的高超声速流动一般不再为层流状态,此时飞行器的气动热环境预测必须考虑湍流效应。当前计算流体力学(CFD)对高超声速流动中湍流的模拟大多是基于RANS(Reynolds-averaged Navier-Stoke)方法,而当流场中存在强激波间断时,现有湍流模型的模拟将会受到严重影响。例如,将当前应用最为广泛的k-ωSST两方程湍流模型应用于美国NASA火星实验室所设计的椭圆钝化再入飞行器(图1)气动热环境预测模拟时,气动热模拟结果与试验数据相比产生了极大的误差,头部驻点区域误差甚至超过了100%(图2)。这是由于湍流模型中相关变量输运方程中的生成项为非守恒形式,强间断附近的高速度梯度值会引起生成项的过度增加,从而引起湍流变量如湍动能和耗散率等经过激波时会发生非物理的变化,进而导致气动热环境的预测结果出现严重的误差。Sinha等人(见参考文献[1]:A.A.Pasha,K.Sinha,Simulation ofhypersonic shock/turbulent boundary-layer interactions using shock-unsteadiness model,Journal of Propulsion and Power 28(1)(2012)46-60)通过研究发现,若能对激波后湍流变量非物理增长的现象进行修正,可有效改善气动热预测模拟的精度。然而,Sinha等人提出的模型由于引入了需要进行积分的边界层厚度等非当地变量,使得其所提出的修正湍流模型仅适用于简单外形流场的预测模拟,难以应用于现代高超声速复杂外形的飞行器。
发明内容
本发明提出一种用于高马赫数强激波流场气动热预测的湍流模型,首先确定流场中存在强间断的区域,随后构造了一个湍流变量生成项衰减函数以消除激波诱导湍流变量非物理的增长现象,并与k-ωSST湍流模型结合获得了一种适用于高马赫数强激波流场气动热预测的湍流模型,实现高马赫数存在强激波间断流场下提高气动热预测精度的目的。
具体的,本发明提供的适用于高马赫数强激波流场气动热预测的湍流模型通过如下步骤得到:
第一步,计算网格节点(I,J,K)的光滑因子。
第二步,计算该网格节点(I,J,K)的衰减函数值,确定流场中强间断的区域。
第三步,衰减函数耦合k-ωSST湍流模型,建立用于高马赫数强激波流场气动热预测的湍流模型。
本发明的优点在于:
(1)计算精度高。现有的湍流模型(如经典的k-ωSST湍流模型)在应用于存在强激波流场的气动热环境预测时,驻点区域计算误差甚至能达到100%以上,而在相同条件下本发明所提出的湍流模型计算精度显著提高,误差能够降低至10%以内。
(2)实用性强。相比于已有的湍流修正模型(如Sinha的修正模型)仅能够适用于平板等简单外形流动的模拟,本发明所提出的基于光滑因子概念的激波检测方法对于复杂外形钝头体飞行器仍能够自动检测激波,并结合衰减函数实现高精度气动热预测模拟,实用性强。
(3)易于融入现代并行化CFD计算程序。本发明所提出的用于高马赫数强激波流场气动热预测的湍流模型并未引入边界层厚度等需要进行积分的非当地变量,模型中变量均为当地变量,融入现代并行化CFD计算程序难度较低。
附图说明
图1为现有技术中的椭圆钝化飞行器几何外形尺寸;
图2为现有技术中的椭圆钝化飞行器对称面下壁面热流分布与实验数据对比;
图3为对于任意一维网格节点J光滑因子构造中模板的选取方式;
图4为本发明所提出的衰减函数在对称面与马赫数分布云图对比;
图5为本发明所提出的湍流模型和原始模型对称面湍流涡粘性分布云图对比;
图6为本发明所提出的湍流模型和原始模型表面热流预测结果对比。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
本发明提出一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法,所述的建立方法中提出一种适用于复杂外形的强激波间断检测方法,用于确定流场中存在强间断的区域,随后基于此构造了一个湍流变量生成项衰减函数以消除激波诱导湍流变量非物理增长的现象,并与k-ωSST湍流模型结合获得了一种适用于高马赫数强激波流场气动热预测的湍流模型。该湍流模型可显著提高高马赫数强激波流场气动热预测模拟的精度,对复杂钝头体外形实用性强并且易于融入现代并行化CFD计算程序。
本发明提供的用于高马赫数强激波流场气动热预测的湍流模型的建立方法,包括以下步骤:
步骤1:计算网格节点(I,J,K)光滑因子。
为了解决现有湍流模型在强激波附近产生的湍流变量非物理增长现象,本发明首先引入了WENO格式构造中的光滑因子概念以实现流场中强间断区域的自动检测。与WENO格式不同的是,本发明采用了无量纲的压强其中p为流场当地压强,U∞为自由来流速度,ρ∞为自由来流密度,上标“-”代表无量纲)作为检测变量,并在数值计算中针对一维空间任意网格节点J选择了三个具有二阶精度的插值模板获得相应的光滑因子,经理论推导所获得的三个插值模板对应的光滑因子的表达式为:
其中,其中pJ为网格节点J所在流场当地压强,IS1J、IS2J、IS3J分别为网格节点J的三个插值模板1、插值模板2、插值模板3对应的光滑因子。对于三维问题,每一维上的光滑因子以同样的方式构造,将J替换为I和K,即得到相应维度上的光滑因子。
在数值模拟过程中对于任意三维网格节点(I,J,K),采用公式(1)中的表达式对每个维度计算三个插值模板各自对应的光滑因子值,然后对每个插值模板内的三个维度光滑因子按公式(2)的方法计算得到采用该插值模板计算得到的光滑因子。
每个插值模板m(m=1,2,3)的光滑因子值可将三个维度上的对应光滑因子值分别平方并求和后开方得到:
然后,在三个插值模板的光滑因子中选取最大值IS:
IS=max(IS1,IS2,IS3) (3)
步骤2:计算该网格节点的衰减函数值。
构造一个湍流变量输运方程生成项衰减函数fD:
fD=1-tanh(camIS) (4)
其中,cam为衰减函数系数,通过理论分析激波前后压强变化规律确定该衰减函数系数为95。
将公式(3)中三个插值模板的光滑因子中的最大值IS带入公式(4),得到该网格节点的衰减函数值。
若流场中存在强间断,由本发明所建立的上述光滑因子三个维度中最大值IS的取值在1左右,而流场越连续,光滑因子越接近于0。因此,可通过网格节点的三个插值模板中光滑因子的最大值判断流场内强间断的区域,进而实现对强间断区域的自动检测,确定流场中需要对湍流变量非物理增长现象抑制的区域。由于避免了边界层厚度等非当地变量的引入,使得这种基于光滑因子概念的激波检测方法适用于任意一般形状的激波。
步骤3:将所建立的衰减函数fD耦合原始k-ωSST湍流模型的生成项,以自动在强激波间断区域抑制湍流变量的非物理增长,而在其他区域保持原始湍流模型计算方式,建立了一种用于高马赫数强激波流场气动热预测的湍流模型,输运方程具体表达式为:
其中,t代表时间,k为湍动能,ω为定耗散率,ρ为当地流场密度,u为当地流场速度,为流场的剪应力张量,下标i,j为维度哑标并遵循爱因斯坦求和约定,μ为分子粘性,μt为湍流涡粘性,β*、σk、β、σω、σω2、γ和a1为k-ωSST湍流模型的系数,F1为混合函数。其模型系数及混合函数具体表达式与k-ωSST湍流模型相同(见参考文献[2]:F.R.Menter,Two-equation eddy-viscosity turbulence models for engineering applications,AIAAJournal,32(8)(1994)1598-1605)。
实施例:选取了NASA火星实验室所设计的椭圆钝化再入飞行器测试外形作为实施例,其几何外形如图1所示。其中空气来流条件为:马赫数6.03,温度58.6K,压强为2091Pa,壁面温度300K,攻角40°;图4为对称面马赫数与本发明所构造的衰减函数分布云图,图5给出了模型修正前后对称面湍流涡粘性分布云图对比。为了方便对比分析,图4和图5中虚线上半图为钝头体下半部分流场绕x轴旋转了180°。由图4可以看出本发明构造的衰减函数fD在三维模拟中实现了对任意形状激波的准确追踪并在激波区域保持函数值为0而在其它区域为1。在图5中,虚线上半图和虚线下半图分别为原始湍流模型和本发明所提出湍流模型的计算结果,可以看到采用衰减函数后湍流模型模拟的湍流涡粘性在激波处不再产生大幅度的增长,而是集中在边界层内增长,表明本发明所提出的衰减函数对激波处的湍流变量的非物理增长现象起到了很好的抑制作用并且没有影响到其在边界层内的正常发展过程。图6给出了分别采用原始湍流模型、本发明提出的湍流模型以及层流计算所得对称面下壁面的热流分布与实验数据对比。相比于原始的湍流模型,本发明提出的衰减函数湍流模型由于消除了激波诱导的湍流变量的非物理增长,从而大幅改善了湍流模型在头部驻点区热流的预测结果。另外,从层流计算结果与实验数据的对比可以发现,实际流动在x/L=0.13(x为流向坐标,L为该飞行器长度位置开始发生了转捩。由于湍流模型本身无法准确预测高超声边界层转捩的过程,因此修正后的湍流模型所预测的转捩位置提前,转捩区域较短,导致下游湍流区域热流预测结果相比于实验数据略微偏高。然而,与原始的湍流模型预测结果相比,本发明提出的湍流模型热流预测结果基本位于实验数据误差带内或接近误差带的上限,热流预测结果显著优于原始湍流模型,显著提高了气动热环境预测模拟的精度。
Claims (3)
1.一种用于高马赫数强激波流场气动热预测的湍流模型的建立方法,其特征在于:包括如下步骤,
第一步,计算网格节点(I,J,K)的光滑因子;
采用无量纲的压强作为检测变量,并在数值计算中针对一维空间任意网格节点J选择了三个具有二阶精度的插值模板获得相应的光滑因子,经理论推导所获得的三个插值模板对应的光滑因子的表达式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>IS</mi>
<mrow>
<mn>1</mn>
<mi>J</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>4</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>3</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>J</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>13</mn>
<mn>12</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>J</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>IS</mi>
<mrow>
<mn>2</mn>
<mi>J</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>13</mn>
<mn>12</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>J</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>IS</mi>
<mrow>
<mn>3</mn>
<mi>J</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>J</mi>
</msub>
<mo>-</mo>
<mn>4</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>13</mn>
<mn>12</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>J</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>J</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,上标“-”代表无量纲压强,IS1J、IS2J、IS3J分别为网格节点J的三个插值模板1、插值模板2、插值模板3对应的光滑因子;对于三维问题,每一维上的光滑因子以同样的方式构造,将J替换为I和K,即得到相应维度上的光滑因子;
在数值模拟过程中对于任意三维网格节点(I,J,K),采用公式(1)中的表达式对每个维度计算三个插值模板各自对应的光滑因子值,然后对每个插值模板内的三个维度光滑因子按公式(2)的方法计算得到采用该插值模板计算得到的光滑因子:
<mrow>
<msub>
<mi>IS</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>IS</mi>
<mrow>
<mi>m</mi>
<mi>I</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>IS</mi>
<mrow>
<mi>m</mi>
<mi>J</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>IS</mi>
<mrow>
<mi>m</mi>
<mi>K</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
然后,在三个插值模板的光滑因子中选取最大值IS:
IS=max(IS1,IS2,IS3) (3)
步骤2:计算该网格节点的衰减函数值;
构造一个湍流变量输运方程生成项衰减函数fD:
fD=1-tanh(camIS) (4)
其中,cam为衰减函数系数;
将公式(3)中三个插值模板的光滑因子中的最大值IS带入公式(4),得到该网格节点的衰减函数值;
步骤3:将所建立的衰减函数fD耦合原始k-ωSST湍流模型的生成项,以自动在强激波间断区域抑制湍流变量的非物理增长,而在其他区域保持原始湍流模型计算方式,建立了一种用于高马赫数强激波流场气动热预测的湍流模型,输运方程具体表达式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;u</mi>
<mi>j</mi>
</msub>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>f</mi>
<mi>D</mi>
</msub>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<msup>
<mi>&beta;</mi>
<mo>*</mo>
</msup>
<mi>&rho;</mi>
<mi>k</mi>
<mi>&omega;</mi>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>k</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;u</mi>
<mi>j</mi>
</msub>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mi>&gamma;</mi>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
</mfrac>
<msub>
<mi>f</mi>
<mi>D</mi>
</msub>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<msup>
<mi>&beta;&rho;&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>&omega;</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&omega;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&rho;&sigma;</mi>
<mrow>
<mi>&omega;</mi>
<mn>2</mn>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<mi>&omega;</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&omega;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,t代表时间,k为湍动能,ω为定耗散率,ρ为当地流场密度,u为当地流场速度,为流场的剪应力张量,下标i,j为维度哑标并遵循爱因斯坦求和约定,μ为分子粘性,μt为湍流涡粘性,β*、σk、β、σω、σω2、γ和a1为k-ωSST湍流模型的系数,F1为混合函数。
2.根据权利要求1所述的一种用于高马赫数强激波流场气动热预测的湍流模型的建立方法,其特征在于:所述的衰减函数系数为95。
3.一种用于高马赫数强激波流场气动热预测的湍流模型,其特征在于:所述的用于高马赫数强激波流场气动热预测的湍流模型,输运方程具体表达式为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;u</mi>
<mi>j</mi>
</msub>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>f</mi>
<mi>D</mi>
</msub>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<msup>
<mi>&beta;</mi>
<mo>*</mo>
</msup>
<mi>&rho;</mi>
<mi>k</mi>
<mi>&omega;</mi>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>k</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;u</mi>
<mi>j</mi>
</msub>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mi>&gamma;</mi>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
</mfrac>
<msub>
<mi>f</mi>
<mi>D</mi>
</msub>
<msubsup>
<mi>&tau;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<msup>
<mi>&beta;&rho;&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mi>&omega;</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&omega;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&rho;&sigma;</mi>
<mrow>
<mi>&omega;</mi>
<mn>2</mn>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<mi>&omega;</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&omega;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,fD为衰减函数:
fD=1-tanh(camIS)
cam为衰减函数系数;IS为网格节点选取的三个插值模板的光滑因子中的最大值;t代表时间,k为湍动能,ω为定耗散率,ρ为当地流场密度,u为当地流场速度,为流场的剪应力张量,下标i,j为维度哑标并遵循爱因斯坦求和约定,μ为分子粘性,μt为湍流涡粘性,β*、σk、β、σω、σω2、γ和a1为k-ωSST湍流模型的系数,F1为混合函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710402397.0A CN107273593A (zh) | 2017-06-01 | 2017-06-01 | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710402397.0A CN107273593A (zh) | 2017-06-01 | 2017-06-01 | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107273593A true CN107273593A (zh) | 2017-10-20 |
Family
ID=60065836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710402397.0A Withdrawn CN107273593A (zh) | 2017-06-01 | 2017-06-01 | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107273593A (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108255781A (zh) * | 2018-01-04 | 2018-07-06 | 北京环境特性研究所 | 一种高超声速目标表面动态温度建模方法 |
CN110083977A (zh) * | 2019-05-14 | 2019-08-02 | 南京大学 | 一种基于深度学习的大气湍流监测方法 |
CN110309552A (zh) * | 2019-06-10 | 2019-10-08 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN111046570A (zh) * | 2019-12-18 | 2020-04-21 | 北京电子工程总体研究所 | 一种刚性圆柱体涡激振动数值模拟的湍流模型 |
CN111859817A (zh) * | 2020-05-22 | 2020-10-30 | 中国航空工业集团公司西安航空计算技术研究所 | 一种基于cfd软件激波模拟能力的飞行器气动选型方法 |
CN112067239A (zh) * | 2020-07-27 | 2020-12-11 | 中国航天空气动力技术研究院 | 一种基于纹影视频自动判断超声速风洞流场建立的方法 |
CN112270037A (zh) * | 2020-10-15 | 2021-01-26 | 北京空天技术研究所 | 高马赫数下内埋弹舱扰流板高度建模及流动控制方法 |
CN112464583A (zh) * | 2020-11-04 | 2021-03-09 | 空气动力学国家重点实验室 | 一种考虑激波和边界层的网格生成方法 |
CN113065275A (zh) * | 2021-03-05 | 2021-07-02 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN113158339A (zh) * | 2021-04-16 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种针对sst湍流模型的湍流长度尺度修正方法 |
CN113420507A (zh) * | 2021-07-08 | 2021-09-21 | 刘维烁 | 流体流动模拟方法、装置以及介质 |
CN113536461A (zh) * | 2021-07-15 | 2021-10-22 | 北京航空航天大学 | 用于高超声速强激波流场气动热预测的湍流模型修正方法 |
CN113850008A (zh) * | 2021-12-02 | 2021-12-28 | 北京航空航天大学 | 飞行器气动特性预测的自适应网格扰动域更新加速方法 |
CN116070554A (zh) * | 2023-04-06 | 2023-05-05 | 中国人民解放军国防科技大学 | 高超声速飞行器气动热载荷计算方法、装置及设备 |
CN116611369A (zh) * | 2023-07-17 | 2023-08-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN117574029A (zh) * | 2024-01-19 | 2024-02-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105893693A (zh) * | 2016-04-21 | 2016-08-24 | 北京航空航天大学 | 一种基于紧致型高分辨率混合格式的可压缩湍流直接数值模拟方法 |
-
2017
- 2017-06-01 CN CN201710402397.0A patent/CN107273593A/zh not_active Withdrawn
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105893693A (zh) * | 2016-04-21 | 2016-08-24 | 北京航空航天大学 | 一种基于紧致型高分辨率混合格式的可压缩湍流直接数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
ZHICHAO ZHANG等: ""A RANS model correction on unphysical over-prediction of turbulent quantities across shock wave"", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 * |
高振勋 等: ""激波诱导湍流模型变量非物理增长的修正模型"", 《第九届全国流体力学学术会议》 * |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108255781A (zh) * | 2018-01-04 | 2018-07-06 | 北京环境特性研究所 | 一种高超声速目标表面动态温度建模方法 |
CN108255781B (zh) * | 2018-01-04 | 2021-05-11 | 北京环境特性研究所 | 一种高超声速目标表面动态温度建模方法 |
CN110083977A (zh) * | 2019-05-14 | 2019-08-02 | 南京大学 | 一种基于深度学习的大气湍流监测方法 |
CN110083977B (zh) * | 2019-05-14 | 2023-06-06 | 南京大学 | 一种基于深度学习的大气湍流监测方法 |
CN110309552B (zh) * | 2019-06-10 | 2023-04-14 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN110309552A (zh) * | 2019-06-10 | 2019-10-08 | 中国航天空气动力技术研究院 | 一种考虑质量引射效应的飞行器湍流预测方法及系统 |
CN111046570A (zh) * | 2019-12-18 | 2020-04-21 | 北京电子工程总体研究所 | 一种刚性圆柱体涡激振动数值模拟的湍流模型 |
CN111046570B (zh) * | 2019-12-18 | 2023-09-12 | 北京电子工程总体研究所 | 一种刚性圆柱体涡激振动数值模拟的湍流模型 |
CN111859817A (zh) * | 2020-05-22 | 2020-10-30 | 中国航空工业集团公司西安航空计算技术研究所 | 一种基于cfd软件激波模拟能力的飞行器气动选型方法 |
CN111859817B (zh) * | 2020-05-22 | 2024-01-16 | 中国航空工业集团公司西安航空计算技术研究所 | 一种基于cfd软件激波模拟能力的飞行器气动选型方法 |
CN112067239A (zh) * | 2020-07-27 | 2020-12-11 | 中国航天空气动力技术研究院 | 一种基于纹影视频自动判断超声速风洞流场建立的方法 |
CN112270037A (zh) * | 2020-10-15 | 2021-01-26 | 北京空天技术研究所 | 高马赫数下内埋弹舱扰流板高度建模及流动控制方法 |
CN112270037B (zh) * | 2020-10-15 | 2023-11-03 | 北京空天技术研究所 | 高马赫数下内埋弹舱扰流板高度建模及流动控制方法 |
CN112464583B (zh) * | 2020-11-04 | 2023-03-14 | 空气动力学国家重点实验室 | 一种考虑激波和边界层的网格生成方法 |
CN112464583A (zh) * | 2020-11-04 | 2021-03-09 | 空气动力学国家重点实验室 | 一种考虑激波和边界层的网格生成方法 |
CN113065275B (zh) * | 2021-03-05 | 2022-05-24 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN113065275A (zh) * | 2021-03-05 | 2021-07-02 | 浙江大学 | 一种飞行器飞行过程中驻点热流的在线预示方法 |
CN113158339A (zh) * | 2021-04-16 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种针对sst湍流模型的湍流长度尺度修正方法 |
CN113420507A (zh) * | 2021-07-08 | 2021-09-21 | 刘维烁 | 流体流动模拟方法、装置以及介质 |
CN113420507B (zh) * | 2021-07-08 | 2023-11-21 | 刘维烁 | 流体流动模拟方法、装置以及介质 |
CN113536461A (zh) * | 2021-07-15 | 2021-10-22 | 北京航空航天大学 | 用于高超声速强激波流场气动热预测的湍流模型修正方法 |
CN113536461B (zh) * | 2021-07-15 | 2023-11-03 | 北京航空航天大学 | 用于高超声速强激波流场气动热预测的湍流模型修正方法 |
CN113850008B (zh) * | 2021-12-02 | 2022-03-11 | 北京航空航天大学 | 飞行器气动特性预测的自适应网格扰动域更新加速方法 |
CN113850008A (zh) * | 2021-12-02 | 2021-12-28 | 北京航空航天大学 | 飞行器气动特性预测的自适应网格扰动域更新加速方法 |
CN116070554B (zh) * | 2023-04-06 | 2023-06-09 | 中国人民解放军国防科技大学 | 高超声速飞行器气动热载荷计算方法、装置及设备 |
CN116070554A (zh) * | 2023-04-06 | 2023-05-05 | 中国人民解放军国防科技大学 | 高超声速飞行器气动热载荷计算方法、装置及设备 |
CN116611369B (zh) * | 2023-07-17 | 2023-09-29 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN116611369A (zh) * | 2023-07-17 | 2023-08-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN117574029A (zh) * | 2024-01-19 | 2024-02-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
CN117574029B (zh) * | 2024-01-19 | 2024-04-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107273593A (zh) | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 | |
Slater | Improvements in modeling 90-degree bleed holes for supersonic inlets | |
Galindo et al. | Development and validation of a radial variable geometry turbine model for transient pulsating flow applications | |
Kornilov et al. | Flat-plate drag reduction with streamwise noncontinuous microblowing | |
CN104794293A (zh) | 风力机尾流计算方法 | |
Casalino et al. | Lattice-Boltzmann simulation of coaxial jet noise generation | |
Unver Kaynak et al. | Transition modeling for low to high speed boundary layer flows with CFD applications | |
Li et al. | Evaluation method of riblets effects and application on a missile surface | |
Aranake et al. | Automated design optimization of a three-dimensional subsonic diffuser | |
Gopalakrishnan et al. | Computation of high-subsonic and transonic flows by a lattice Boltzmann method | |
Brunet et al. | Experimental and numerical investigations of vortex generators effects | |
Bensayah et al. | Heat transfer in turbulent boundary layers of conical and bell shaped rocket nozzles with complex wall temperature | |
Ouahabi et al. | Application of the turbulence models at low reynolds number for horizontal-axis wind turbine design by using computational fluid dynamics (CFD) | |
Krause et al. | Numerical computations for designing a scramjet intake | |
Baydar et al. | Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet | |
Choe et al. | Supersonic Inlet Design using Bleed Boundary Condition with Porosity Variation and Expansion Wave | |
Xiang et al. | Computational grid dependency in CFD simulation for heat transfer | |
Tao et al. | Catastrophe model for supersonic inlet start/unstart | |
Marcos et al. | Computational and experimental comparison of a powered lift, upper surface blowing configuration | |
Eghlima et al. | Effect of nose shape on the shock standoff distance at nearsonic flows | |
Halila et al. | A numerical study on transitional flows by means of a correlation-based transition model | |
Choe et al. | Numerical investigation of bleed effects on supersonic inlet under various bleed and inlet conditions | |
Yi et al. | Design Optimization of Vortex Generators for a Junction Vortex of Wing-Body Configuration by Discrete Adjoint Approach | |
Schuman | A Computational Evaluation of Transonic Wind Tunnel Wall Interference on High Aspect Ratio Models in the Arnold Engineering Development Complex 16 Foot Transonic Tunnel | |
Kächele et al. | Numerical simulations of an intake-compressor system |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20171020 |