CN112182985B - 一种控制细长回转体边界层保持层流不分离流动的方法 - Google Patents
一种控制细长回转体边界层保持层流不分离流动的方法 Download PDFInfo
- Publication number
- CN112182985B CN112182985B CN202010845434.7A CN202010845434A CN112182985B CN 112182985 B CN112182985 B CN 112182985B CN 202010845434 A CN202010845434 A CN 202010845434A CN 112182985 B CN112182985 B CN 112182985B
- Authority
- CN
- China
- Prior art keywords
- revolving body
- flow
- slender
- elongated
- revolution
- 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.)
- Active
Links
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02T70/00—Maritime or waterways transport
- Y02T70/10—Measures concerning design or construction of watercraft hulls
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种控制细长回转体边界层保持层流不分离流动的方法,涉及流体控制技术领域,该方法通过优化细长回转体的形状和表面特性,尽量增大顺压梯度范围,且在逆压梯度区域内降低流体流动速度的变化率,有利于使细长回转体在设定速度范围内保持层流流动,避免细长回转体表面发生流动分离或转捩,可以大幅减小细长回转体的摩擦阻力和压差阻力,有很好的减阻效果;而且,实现减阻的同时并没有消耗额外的能量,达到了节能和提高推进效率的双重效果;此外,没有发生边界层分离的细长回转体的水动力噪声也会降低。
Description
技术领域
本发明涉及流体控制技术领域,具体涉及一种控制细长回转体边界层保持层流不分离流动的方法。
背景技术
航行体在水中运动时,由于水的粘性的存在,会在其表面附近形成边界层,如图1所示,为航行体(为细长回转体)的传统标准形状,即常见的三段式设计,由进流段、平行中体和去流段组成,其边界层很容易发生转捩或流动分离,造成航行体的阻力增加;边界层作用于航行体的阻力主要分为摩擦阻力和压差阻力,阻力大小受流体运动状态的影响:边界层在航行体表面的逆压梯度区域内可能会发生流动分离,这时的压差阻力会大幅增加;湍流边界层虽然不易产生流动分离,却增加了航行体所受的摩擦阻力。因此,使航行体在层流边界层中保持不分离运动是航行体在摩擦阻力不增加的前提下、所受形状阻力最小的方法。
目前已有很多关于避免和延缓边界层转捩和流动分离的研究,其中不乏主动流动控制方法(如边界层抽吸方法),但主动流动控制方法意味着额外的能量损失,所以如果单以航行体的表面特性实现边界层在层流状态下的流动不分离会有很大的经济优势。
发明内容
有鉴于此,本发明提供了一种控制细长回转体边界层保持层流不分离流动的方法,通过优化细长回转体的形状和表面特性,使细长回转体在设定速度范围内保持层流流动,避免细长回转体表面发生流动分离或转捩。
本发明的技术方案为:一种控制细长回转体边界层保持层流不分离流动的方法,包括以下步骤:
第一步:预设细长回转体的长度、长细比和航速指标,根据细长体理论和势流理论得到细长回转体表面流场的势流方程;
第二步:求解第一步中的势流方程;
第三步:引入高阶连续物面曲线,计算其中的中间变量E的值,并计算细长回转体表面的压强系数Cp(x);
第四步:确定高阶连续物面曲线中两段分段函数的交点x*和待定系数的初始值;
第五步:求解细长回转体在x*处的表面流函数及压强系数的表达式;
第六步:根据细长回转体的最大半径位置和压强系数最小位置,以及细长回转体表面的流函数Ψ(x,R(x))=0,对第五步中的流函数和压强系数进行反复迭代,动态确定x*和待定系数的值,直至达到优化条件的设定精度要求,计算得到细长回转体所需的表面流函数及压强系数;
其中,优化条件为:
其中,x为柱坐标系下细长回转体的轴向坐标,r为柱坐标系下细长回转体的径向坐标,ε为细长回转体的厚度系数,Rmax为细长回转体最大直径;
第七步:通过积分估算细长回转体的体积,判断细长回转体边界层的实际雷诺数与临界雷诺数之间的关系,当实际雷诺数小于临界雷诺数时,细长回转体边界层能够保持为层流且不发生流动分离;否则,返回第一步,重新选择细长回转体的长度、长细比和航速指标。
优选地,第一步中,所述势流方程为:
其中,M为马赫数,M=U∞/Ua<0.9;U∞为细长回转体的航行速度,Ua为声音在水中的传播速度,Φ为无量纲化速度势。
优选地,第二步中,所述求解包括:
Φ(x,r,ε)=x+ε2lnεA(x)+ε2{A(x)ln(ωr*)+B(x)}+O(ε4ln2ε) (2)
优选地,第三步包括:
对第二步中的(2)式进行伯努利积分,可得细长回转体表面的压强系数为:
忽略(3)式中的无穷小量O(ε2),引入高阶连续物面曲线的函数关系式为:
得到压强系数为:
其中,a,c,a1为待定系数。
优选地,第四步中,所述确定的过程包括以下步骤:
1)根据细长回转体外形的单调性,确定a,c,a1,x*的取值范围为:
2)假设细长回转体的外形的最大半径位置和压强系数最小位置相同,由此得到c,x*的初始值;
在交点x*处,利用半径及压强系数连续的条件得到a,c,a1,x*之间的关系为:
优选地,第五步中,所述求解过程包括:
在x*处,对于不可压缩流通过使用分布在对称轴上的源和汇获得拉普拉斯方程的精确解:
即亚音速下的流函数表示如下:
其中,源和汇的强度为:
联立(4)、(8)和(9)式,得到流函数为:
Ψ(x,r)=0.5r2+β1u(x*)-0.75E{a[F1(x*)-F1(0)]+a1[F1(1)-F1(x*)] (10)
+2(ax+2c)[F2(x*)-F2(0)]+2a1(x-1)[F2(1)-F2(x*)]}
可得,细长回转体表面的压强系数为:
优选地,第七步中,实际雷诺数的计算公式如下:
临界雷诺数的计算公式如下:
其中,V为细长回转体的体积,ν为流体的运动粘度。
有益效果:
本发明的方法通过优化细长回转体的形状和表面特性,尽量增大顺压梯度范围,且在逆压梯度区域内降低流体流动速度的变化率,有利于使细长回转体在设定速度范围内保持层流流动,避免细长回转体表面发生流动分离或转捩,可以大幅减小细长回转体的摩擦阻力和压差阻力,有很好的减阻效果;而且,实现减阻的同时并没有消耗额外的能量,达到了节能和提高推进效率的双重效果;此外,没有发生边界层分离的细长回转体的水动力噪声也会降低。
附图说明
图1为现有技术中细长回转体的标准形状(一般为三段式结构,长细比为6)。
图2为本发明中控制边界层保持层流不分离流动对应细长回转体的一种优化形状(长细比为6)。
图3为本发明细长回转体形状的压力系数随长度变化的示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本实施例提供了一种控制细长回转体边界层保持层流不分离流动的方法,通过优化细长回转体的形状和表面特性,使细长回转体在设定速度范围内保持层流流动,避免细长回转体表面发生流动分离或转捩。
该方法包括以下步骤:
第一步:预设细长回转体的长度、长细比(取较大值)和航速指标,根据细长体理论和势流理论得到细长回转体表面流场的势流方程;
具体为:因为流体发生流动分离的两个条件为:(1)流体具有粘性;(2)物体表面存在逆压梯度,由达朗贝尔悖论可以知道不存在整个表面都是顺压梯度的物体,所以确定物体(细长回转体)外形的压力分布特点是流动控制的关键;
当长细比较大时,假设边界层厚度足够小,甚至可以忽略,且假定在边界层之外的流体是无粘的,对于水下细长回转体则假设都处于亚音速范围内,此时柱坐标系下的势流方程为:
其中,M为马赫数,M=U∞Ua<0.9;U∞为细长回转体的航行速度,Ua为声音在水中的传播速度,Φ为无量纲化速度势,x为柱坐标系下细长回转体的轴向坐标,r为柱坐标系下细长回转体的径向坐标;
第二步:求解势流方程:
Φ(x,r,ε)=x+ε2lnεA(x)+ε2{A(x)ln(ωr*)+B(x)}+O(ε4ln2ε) (2)
第三步:引入高阶连续物面曲线,计算其中的中间变量E的值,并计算细长回转体表面的压强系数;
对(2)式进行伯努利积分,可得细长回转体表面的压强系数为:
忽略(3)式中的无穷小量O(ε2),采用高阶连续物面曲线以得到光滑的压力系数曲线,取高阶连续物面曲线的函数关系式为:
得到压强系数为:
其中,a,c,a1为待定系数;x*为(4)式中两段分段函数的交点;
第四步:确定高阶连续物面曲线中两段分段函数的交点x*和待定系数的初始值;
1)根据细长回转体外形的单调性,确定a,c,a1,x*的取值范围为:
2)假设优化后的细长回转体的外形的最大半径位置和压强系数最小位置相同,由此得到c,x*的初始值;
在交点x*处,利用半径及压强系数连续的条件得到a,c,a1,x*之间的关系为:
第五步:求解细长回转体在x*处的表面流函数及压强系数的表达式;
因为上述细长回转体表面函数的选取,在x*处细长回转体的外形函数不连续,会导致边界层流动分离,所以在此处需要得到势流方程的精确解,对于不可压缩流(M=0)通过使用分布在对称轴上的源和汇获得拉普拉斯方程的精确解:
即亚音速下的流函数表示如下:
其中,源和汇的强度为:
联立(4)、(8)和(9)式,得到流函数为:
Ψ(x,r)=0.5r2+β1u(x*)-0.75E{a[F1(x*)-F1(0)]+a1[F1(1)-F1(x*)] (10)
+2(ax+2c)[F2(x*)-F2(0)]+2a1(x-1)[F2(1)-F2(x*)]}
可得,所设计的细长回转体表面的压强系数为:
第六步:根据细长回转体的最大半径位置和压强系数最小位置,以及细长回转体表面的流函数Ψ(x,R(x))=0,对(10)式中的流函数方程和(11)式中的压强系数进行反复迭代优化,直至达到优化条件的设定精度要求(无量纲化下,精度控制在10e-4以下),计算得到所需的细长回转体外形流函数及其表面压力系数;
其中,优化条件为:
如图2所示,为一种优化后的细长回转体外形,它不再存在平行中体,由进流段和去流段两部分构成,这两部分分别由两个三次函数拟合而成;如图3所示,为细长回转体表面压强系数随长度的变化关系,由图中可以看出,优化后的细长回转体外形去流段占比大、尾椎角非常小,尽可能使细长回转体保证尖尾、减缓速度梯度,从而尽量避免流动分离;
第七步:通过积分估算细长回转体的体积,判断细长回转体边界层的实际雷诺数与临界雷诺数之间的关系,当实际雷诺数小于临界雷诺数时,细长回转体边界层能够保持为层流且不发生流动分离;否则,返回第一步,重新选择细长回转体的长度、长细比和航速指标;
其中,实际雷诺数的计算公式如下:
临界雷诺数的计算公式如下:
其中,V为细长回转体的体积,ν为流体的运动粘度;
现有技术中标准型细长回转体的体积阻力系数为其中,S为表面积;优化后的细长回转体的体积阻力系数为以航速15m/s,总长1.2m,长细比为6为例,现有技术中标准型细长回转体的体积阻力系数(0.01892)约为优化后的细长回转体的体积阻力系数(0.002477)的八倍,减阻效果明显。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种控制细长回转体边界层保持层流不分离流动的方法,其特征在于,包括以下步骤:
第一步:预设细长回转体的长度、长细比和航速指标,根据细长体理论和势流理论得到细长回转体表面流场的势流方程为:
其中,M为马赫数,M=U∞/Ua<0.9;U∞为细长回转体的航行速度,Ua为声音在水中的传播速度,Φ为无量纲化速度势;
第二步:求解第一步中的势流方程,所述求解包括:
Φ(x,r,ε)=x+ε2lnεA(x)+ε2{A(x)ln(ωr*)+B(x)}+O(ε4ln2ε) (2)
第三步:引入高阶连续物面曲线,对第二步中的(2)式进行伯努利积分,可得细长回转体表面的压强系数为:
忽略(3)式中的无穷小量O(ε2),引入高阶连续物面曲线的函数关系式为:
得到压强系数为:
其中,a,c,a1为待定系数;
第四步:确定高阶连续物面曲线中两段分段函数的交点x*和待定系数的初始值,所述确定的过程包括以下步骤:
1)根据细长回转体外形的单调性,确定a,c,a1,x*的取值范围为:
2)假设细长回转体的外形的最大半径位置和压强系数最小位置相同,由此得到c,x*的初始值;
在交点x*处,利用半径及压强系数连续的条件得到a,c,a1,x*之间的关系为:
第五步:求解细长回转体在x*处的表面流函数及压强系数的表达式,所述求解过程包括:
在x*处,对于不可压缩流通过使用分布在对称轴上的源和汇获得拉普拉斯方程的精确解:
即亚音速下的流函数表示如下:
其中,源和汇的强度为:
联立(4)、(8)和(9)式,得到流函数为:
Ψ(x,r)=0.5r2+β1u(x*)-0.75E{a[F1(x*)-F1(0)]+a1[F1(1)-F1(x*)] (10)+2(ax+2c)[F2(x*)-F2(0)]+2a1(x-1)[F2(1)-F2(x*)]}
可得,细长回转体表面的压强系数为:
第六步:根据细长回转体的最大半径位置和压强系数最小位置,以及细长回转体表面的流函数Ψ(x,R(x))=0,对第五步中的流函数和压强系数进行反复迭代,动态确定x*和待定系数的值,直至达到优化条件的设定精度要求,计算得到细长回转体所需的表面流函数及压强系数;
其中,优化条件为:
其中,x为柱坐标系下细长回转体的轴向坐标,r为柱坐标系下细长回转体的径向坐标,ε为细长回转体的厚度系数,Rmax为细长回转体最大直径;
第七步:通过积分估算细长回转体的体积,判断细长回转体边界层的实际雷诺数与临界雷诺数之间的关系,当实际雷诺数小于临界雷诺数时,细长回转体边界层能够保持为层流且不发生流动分离;否则,返回第一步,重新选择细长回转体的长度、长细比和航速指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010845434.7A CN112182985B (zh) | 2020-08-20 | 2020-08-20 | 一种控制细长回转体边界层保持层流不分离流动的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010845434.7A CN112182985B (zh) | 2020-08-20 | 2020-08-20 | 一种控制细长回转体边界层保持层流不分离流动的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112182985A CN112182985A (zh) | 2021-01-05 |
CN112182985B true CN112182985B (zh) | 2022-08-09 |
Family
ID=73924147
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010845434.7A Active CN112182985B (zh) | 2020-08-20 | 2020-08-20 | 一种控制细长回转体边界层保持层流不分离流动的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112182985B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4664345A (en) * | 1983-11-24 | 1987-05-12 | Messerschmitt-Boelkow-Blohm Gesellschaft Mit Beschrankter Haftung | Method for stabilizing laminar separated boundary layers |
US6131529A (en) * | 1997-05-31 | 2000-10-17 | The East Group | Water going vessel hull and method for hull design |
CN101318553A (zh) * | 2008-05-28 | 2008-12-10 | 中国航天空气动力技术研究院 | 一种机翼流动控制装置 |
CN107487440A (zh) * | 2017-09-21 | 2017-12-19 | 贵州大学 | 一种无人侦察机的可变形伸缩储油机翼 |
CN107651165A (zh) * | 2017-09-21 | 2018-02-02 | 贵州大学 | 一种无人侦察机的厚度可调机翼 |
CN107878728A (zh) * | 2016-09-29 | 2018-04-06 | 北京航空航天大学 | 机翼结构及飞行器 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9542506B2 (en) * | 2012-11-13 | 2017-01-10 | Exa Corporation | Computer simulation of physical processes including modeling of laminar-to-turbulent transition |
US20170206291A1 (en) * | 2016-01-20 | 2017-07-20 | Soliton Holdings Corporation, Delaware Corporation | Method for computational fluid dynamics and apparatuses for jet-effect use |
-
2020
- 2020-08-20 CN CN202010845434.7A patent/CN112182985B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4664345A (en) * | 1983-11-24 | 1987-05-12 | Messerschmitt-Boelkow-Blohm Gesellschaft Mit Beschrankter Haftung | Method for stabilizing laminar separated boundary layers |
US6131529A (en) * | 1997-05-31 | 2000-10-17 | The East Group | Water going vessel hull and method for hull design |
CN101318553A (zh) * | 2008-05-28 | 2008-12-10 | 中国航天空气动力技术研究院 | 一种机翼流动控制装置 |
CN107878728A (zh) * | 2016-09-29 | 2018-04-06 | 北京航空航天大学 | 机翼结构及飞行器 |
CN107487440A (zh) * | 2017-09-21 | 2017-12-19 | 贵州大学 | 一种无人侦察机的可变形伸缩储油机翼 |
CN107651165A (zh) * | 2017-09-21 | 2018-02-02 | 贵州大学 | 一种无人侦察机的厚度可调机翼 |
Non-Patent Citations (6)
Title |
---|
"Flow separation from submarine shaped bodies of revolution in steady turning";A.G.L.Holloway等;《Ocean Engineering》;20150905;第108卷;第426-438页 * |
"TVC effects on flow separation on slender cylinders";A.Mehmood等;《Thermophysics and Aeromechanics》;20180109;第25卷(第4期);第507-513页 * |
"基于粘性伴随方法的旋翼先进气动外形优化设计分析";吴琪;《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》;20160315(第03期);第2-6章 * |
"水下高速运动体超空泡及可靠性研究";朱荣荣;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20070415(第04期);第2-4章 * |
"空化器设计及超空泡参数控制";傅慧萍 等;《舰船科学技术》;20031031;第25卷(第5期);第49-51页 * |
"超空泡外形和阻力及空化器形状优化研究";刘磊;《中国优秀硕士学位论文全文数据库 基础科学辑》;20070415(第04期);第2-4章 * |
Also Published As
Publication number | Publication date |
---|---|
CN112182985A (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Batchelor | A proposal concerning laminar wakes behind bluff bodies at large Reynolds number | |
Radespiel et al. | Efficient cell-vertex multigrid scheme for the three-dimensional Navier-Stokes equations | |
Anderson et al. | Airfoil design on unstructured grids for turbulent flows | |
Bontempo et al. | Performance analysis of ducted marine propellers. Part I–Decelerating duct | |
CN112182985B (zh) | 一种控制细长回转体边界层保持层流不分离流动的方法 | |
Liu et al. | A numerical study on the natural transition locations in the flat-plate boundary layers on superhydrophobic surfaces | |
Nesteruk | Rigid bodies without boundary-layer separation | |
Viswanath | Drag reduction of afterbodies by controlled separated flows | |
CN111723447A (zh) | 一种动静腔流动微沟槽减阻设计方法 | |
Issa | Rise of total pressure in frictional flow | |
Kim et al. | Numerical analysis of supercavitating flow around axisymmetric underwater objects using a viscous-potential method | |
Garrett | Linear growth rates of types I and II convective modes within the rotating-cone boundary layer | |
Nesteruk et al. | Turbulent skin-friction drag on a slender body of revolution and Gray's Paradox | |
Liu et al. | Influence of Reynolds number on the natural transition of boundary layers over underwater axisymmetric bodies | |
Kim et al. | Design of propeller geometry using streamline-adapted blade sections | |
Amromin | Vehicles drag reduction with control of critical Reynolds number | |
Zheng et al. | Computational analysis of conical forebody flow at high alpha with transitional model | |
Liu et al. | Numerical investigation of the natural transition in boundary layers on underwater axisymmetric bodies with superhydrophobic surfaces | |
Han et al. | Role of wall roughness on interaction of leakage flow and main flow in a mixed flow pump with tip clearance | |
Tatsukawa et al. | Aerodynamic design exploration for reusable launch vehicle using genetic algorithm with navier stokes solver | |
Pearce et al. | Numerical analysis of base-ventilated intercepted supercavitating hydrofoil sections | |
Celik et al. | Investigation of the optimum duct Geometry for A Passenger ferry | |
CN111737836B (zh) | 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 | |
Hocking et al. | The effect of gravity on flow past a semi-circular cylinder with a constant pressure wake | |
Dong et al. | Numerical Study on Flap Style Rudder Hydrodynamic Characteristics With Different Connections |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |