CN111859532B - 考虑高超声速化学非平衡效应的改进热壁修正方法 - Google Patents
考虑高超声速化学非平衡效应的改进热壁修正方法 Download PDFInfo
- Publication number
- CN111859532B CN111859532B CN202010547613.2A CN202010547613A CN111859532B CN 111859532 B CN111859532 B CN 111859532B CN 202010547613 A CN202010547613 A CN 202010547613A CN 111859532 B CN111859532 B CN 111859532B
- Authority
- CN
- China
- Prior art keywords
- wall
- heat flow
- temperature
- heat
- field data
- 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
- 239000000126 substance Substances 0.000 title claims abstract description 17
- 230000000694 effects Effects 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 title abstract description 32
- 238000012937 correction Methods 0.000 title abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims abstract description 30
- 238000009413 insulation Methods 0.000 claims abstract description 22
- 238000009792 diffusion process Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000009471 action Effects 0.000 claims description 5
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000003197 catalytic effect Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 238000002715 modification method Methods 0.000 claims description 3
- 238000007254 oxidation reaction Methods 0.000 claims description 3
- 238000006557 surface reaction Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 abstract description 7
- 238000010168 coupling process Methods 0.000 abstract description 7
- 238000005859 coupling reaction Methods 0.000 abstract description 7
- 238000010438 heat treatment Methods 0.000 abstract description 7
- 239000007787 solid Substances 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 230000002925 chemical effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000003779 heat-resistant material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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
- 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)
- Chemical & Material Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Analytical Chemistry (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明涉及热壁修正方法技术领域，公开了一种用于高速飞行器气动加热数值预测的考虑高超声速化学非平衡效应的改进热壁修正方法，本方法将对流热流按物理过程的贡献进行分解，其中温度梯度引起的热流部分满足传统线性的热壁修正方法，而反应扩散引起的热流部分与壁面高温化学作用相关，根据壁面化学反应机制做相应处理。分别求解温度梯度引起的热流和反应扩散引起的热流，再将求解到的热流相加，得到壁面总热流。本发明提出的方法只需要求解冷壁和绝热壁条件下的高超声速流场，不需要流固耦合迭代求解，计算效率高。
Description
技术领域
本发明涉及热壁修正方法技术领域，尤其涉及一种考虑高超声速化学非平衡效应的改进热壁修正方法。
背景技术
在利用计算流体力学技术预测高超声速气动加热特性时，需要解决不同壁面热状态以及与防热结构耦合条件下飞行器所受的气动热载荷，将耗费巨大的计算成本和计算机资源，影响飞行器气动设计周期。针对相同来流和布局条件，考虑到壁面热流和温度在所关注的一定温度范围内近似满足线性关系，即描述对流换热系数的斯坦顿数为定值，当前工程设计中通常通过对冷壁热流进行修正来获得热壁热流或与防热结构耦合条件下的热流，即热壁修正方法。
上述常规的热壁修正方法是针对量热完全气体而发展起来的，未考虑高超声速流动诱发激波后化学非平衡效应带来的影响，需要针对化学非平衡效应的物理机理作改进。
发明内容
本发明所要解决的技术问题是：针对上述存在的问题，提供了一种用于高速飞行器气动加热数值预测的考虑高超声速化学非平衡效应的改进热壁修正方法，该方法通过对壁面热流的来源进行物理过程分解，解决一般壁面条件(热壁或与防热结构耦合)下高超声速化学非平衡流动气动加热快速预测的问题。本方法将对流热流按物理过程的贡献进行分解，其中温度梯度引起的热流部分满足传统线性的热壁修正方法，而反应扩散引起的热流部分与壁面高温化学作用相关，根据壁面化学反应机制做相应处理。
本发明采用的技术方案如下：考虑高超声速化学非平衡效应的改进热壁修正方法，包括以下步骤：
步骤1：根据计算需求，通过数值计算方法，获得壁面热流的流场数据；
步骤2：根据步骤1得到的流场数据，计算获得壁面温度梯度引起的热流；
步骤3：根据步骤1得到的流场数据，计算获得壁面反应扩散引起的热流；
步骤4：将温度梯度引起的热流与反应扩散引起的热流相加，获得壁面的热壁总热流。
通过上述技术手段，将热壁总热流分为两部分进行计算，简化了计算步骤，提高了计算效率，而且在计算壁面反应扩散引起的热流后，得到的热壁总热流更为精准。
进一步的，所述步骤1中，所述流场数据包括冷壁流场数据和绝热壁流场数据。
进一步的，所述步骤2具体包括：
步骤21：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁温度、冷壁热流以及绝热壁温度；
步骤22：通过线性插值，获得给定温度条件下的温度梯度引起的热流。
进一步的，所述步骤22中，温度梯度引起的热流的计算公式如下：
其中，qtemp,h表示所求的温度梯度引起的热流，qtemp,c表示冷壁热流，Tc表示冷壁温度，Tad表示绝热壁温度，Th表示热壁温度。
基于传统的热壁修正方法，求得给定温度条件下的温度梯度引起的热流。
进一步的，所述步骤3具体包括：
步骤31：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁浓度、冷壁温度、绝热壁浓度和绝热壁温度；
步骤32：根据冷壁和绝热壁所处温度条件下的浓度，结合壁面反应质量作用定理，获得热壁条件下的组分浓度和焓值；
步骤33：根据壁面化学反应机理以及求得的热壁条件下的组分和焓值，求解获得壁面反应扩散引起的热流。
进一步的，所述步骤33中，反应扩散引起的热流的计算公式如下：
其中，qchem表示反应扩散引起的热流，所述J表示壁面催化氧化反应引起的质量通量，所述Y表示组分质量分数，所述T表示当地温度，hs表示焓值。
与现有技术相比，采用上述技术方案的有益效果为：本方法只需要求解冷壁和绝热壁条件下的高超声速流场，不需要流固耦合迭代求解，计算效率高。
附图说明
图1是本发明的考虑化学非平衡效应的改进热壁修正方法适用的物理问题示意图。
图2是本发明的考虑化学非平衡效应的改进热壁修正方法的流程图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明提供了一种考虑高超声速化学非平衡效应的改进热壁修正方法，该方法适用的物理问题如图1所示。
本发明的改进热壁修正方法是基于壁面高温化学效应对壁面热流和壁面温度线性度的影响规律而发展起来的。对化学非平衡效应的研究表明，尽管反应扩散引起的对流热流与壁面温度的线性度变差，但温度梯度引起的对流热流与壁面温度依旧保持较高的线性性。因此，可以将对流热流按物理过程的贡献进行分解，其中温度梯度引起的热流部分满足传统线性的热壁修正方法，而反应扩散引起的热流部分与壁面高温化学作用相关，根据壁面化学反应机制做相应处理。
首先，将壁面热流分为温度梯度和反应扩散引起的对流热流，具体公式如下：
qtot,h＝qtemp,h+qchem,h
其中温度梯度引起的热流部分为：
而反应扩散引起的热流部分为：
其中，J为壁面催化氧化反应引起的质量通量，h为气体介质的焓值。
由于qchem是由各个壁面化学反应引起的壁面法向组分扩散造成的，在中等或弱化学反应条件下(典型防热材料的表面催化反应基本处于弱化学反应区段)，组分扩散主要与壁面化学反应速率kf正相关。因此，qchem的热壁效应，与当地温度T和组分质量分数Y有关，可由高焓环境下壁面化学反应独立求得(而非耦合高超声速流场)；而qtemp,h可由改进的热壁修正公式求解。进而根据冷壁的非平衡气动加热结果预测热壁的非平衡气动加热量，获得改进的热壁修正方法。
如图2所示，本发明方法的具体步骤如下：
步骤1：根据计算需求，通过数值计算方法分别获得冷壁流场数据和绝热壁流场数据；
步骤2：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁温度Tc和热流qtemp,c以及绝热壁温度Tad；
步骤3：基于传统的热壁修正方法，通过线性插值，获得给定温度条件下的温度梯度引起的热流qtemp,h；
温度梯度引起的热流qtemp,h的计算公式如下：
步骤4：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁浓度Xc和冷壁温度Tc，以及绝热壁浓度Xad和绝热壁温度Tad；
步骤5：根据冷壁和绝热壁所处温度条件下的浓度Xc和Xad，结合壁面反应质量作用定理Kp，获得热壁条件下的组分浓度Xh和焓值hh；
步骤6：根据壁面化学反应机理kf以及热壁条件下的组分Xh和焓值hh，求解壁面反应扩散引起的热流qchem,h；
反应扩散引起的热流qchem,h的计算公式如下：
步骤7：将步骤3和步骤6获得的热流相加，获得最终需要的热壁总热流qtot,h。具体如下：
qtot,h＝qtemp,h+qchem,h。
本方法通过对壁面热流的来源进行物理过程分解，解决一般壁面条件(热壁或与防热结构耦合)下高超声速化学非平衡流动气动加热快速预测的问题。该方法只需要求解冷壁和绝热壁条件下的高超声速流场，不需要流固耦合迭代求解，计算效率高。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合，以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员，在不脱离本发明的精神所做的非实质性改变或改进，都应该属于本发明权利要求保护的范围。
Claims (1)
1.考虑高超声速化学非平衡效应的改进热壁修正方法，其特征在于，包括以下步骤：
步骤1：根据计算需求，通过数值计算方法，获得壁面热流的流场数据；
所述步骤1中，所述流场数据包括冷壁流场数据和绝热壁流场数据；
步骤2：根据步骤1得到的流场数据，计算获得壁面温度梯度引起的热流；
所述步骤2具体包括：
步骤21：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁温度、冷壁热流以及绝热壁温度；
步骤22：通过线性插值，获得给定温度条件下的温度梯度引起的热流；
所述步骤22中，温度梯度引起的热流的计算公式如下：
其中，表示所求的温度梯度引起的热流，/>表示冷壁热流，/>表示冷壁温度，/>表示绝热壁温度，/>表示热壁温度；
步骤3：根据步骤1得到的流场数据，计算获得壁面反应扩散引起的热流；
步骤4：将温度梯度引起的热流与反应扩散引起的热流相加，获得壁面的热壁总热流；
所述步骤3具体包括：
步骤31：根据冷壁流场数据和绝热壁流场数据，分别获得冷壁浓度、冷壁温度、绝热壁浓度和绝热壁温度；
步骤32：根据冷壁和绝热壁所处温度条件下的浓度，结合壁面反应质量作用定理，获得热壁条件下的组分浓度和焓值；
步骤33：根据壁面化学反应机理以及求得的热壁条件下的组分和焓值，求解获得壁面反应扩散引起的热流；
所述步骤33中，反应扩散引起的热流的计算公式如下：
其中，表示反应扩散引起的热流，所述/>表示壁面催化氧化反应引起的质量通量，所述/>表示组分质量分数，所述/>表示当地温度，/>表示焓值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010547613.2A CN111859532B (zh) | 2020-06-16 | 2020-06-16 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010547613.2A CN111859532B (zh) | 2020-06-16 | 2020-06-16 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111859532A CN111859532A (zh) | 2020-10-30 |
| CN111859532B true CN111859532B (zh) | 2023-11-28 |
Family
ID=72987208
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010547613.2A Active CN111859532B (zh) | 2020-06-16 | 2020-06-16 | 考虑高超声速化学非平衡效应的改进热壁修正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111859532B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112989485B (zh) * | 2021-01-18 | 2022-04-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于热壁修正的沿弹道热流插值方法 |
| CN113065201B (zh) * | 2021-05-08 | 2022-03-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种考虑滑移修正的辐射平衡温度计算方法 |
| CN113782106B (zh) * | 2021-07-29 | 2024-04-16 | 空气动力学国家重点实验室 | 一种用于cfd的高焓离解气体表面催化反应速率的求解方法 |
| CN114818149B (zh) * | 2022-06-28 | 2022-11-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 飞行器表面壁面热流分解方法、分析方法、介质及设备 |
| CN115312139B (zh) * | 2022-09-23 | 2023-01-13 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高超声速流动化学反应模型数据的存取与转换方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106508020B (zh) * | 2010-12-10 | 2014-10-22 | 上海机电工程研究所 | 一种可用于工程设计的复杂飞行器气动加热计算方法 |
| CN105203591A (zh) * | 2015-09-24 | 2015-12-30 | 中国航天空气动力技术研究院 | 航天器试验模型热壁热流密度的测量方法 |
| CN107808065A (zh) * | 2017-11-23 | 2018-03-16 | 南京航空航天大学 | 三维复杂外形高速飞行器流‑固‑热快速计算方法 |
| CN107832494A (zh) * | 2017-10-13 | 2018-03-23 | 南京航空航天大学 | 高超声速飞行器前缘流‑热‑固一体化计算方法 |
| WO2020015375A1 (zh) * | 2018-07-18 | 2020-01-23 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
-
2020
- 2020-06-16 CN CN202010547613.2A patent/CN111859532B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106508020B (zh) * | 2010-12-10 | 2014-10-22 | 上海机电工程研究所 | 一种可用于工程设计的复杂飞行器气动加热计算方法 |
| CN105203591A (zh) * | 2015-09-24 | 2015-12-30 | 中国航天空气动力技术研究院 | 航天器试验模型热壁热流密度的测量方法 |
| CN107832494A (zh) * | 2017-10-13 | 2018-03-23 | 南京航空航天大学 | 高超声速飞行器前缘流‑热‑固一体化计算方法 |
| CN107808065A (zh) * | 2017-11-23 | 2018-03-16 | 南京航空航天大学 | 三维复杂外形高速飞行器流‑固‑热快速计算方法 |
| WO2020015375A1 (zh) * | 2018-07-18 | 2020-01-23 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
Non-Patent Citations (4)
| Title |
|---|
| 持续气动加热环境下的结构热载荷分析与应用;雷桂林 等;科学技术与工程;第13卷(第12期);第3343-3349, 3371页 * |
| 表面催化特性对火星进入气固耦合热效应的影响研究;杨肖峰 等;中国科学:技术科学;第48卷(第09期);第939-949页 * |
| 高超声速气动热的耦合计算方法研究;赵晓利 等;科学技术与工程;第10卷(第22期);第5450-5455, 5461页 * |
| 高超声速流动壁面催化复合气动加热特性;苗文博 等;宇航学报;第34卷(第03期);第442-446页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111859532A (zh) | 2020-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111859532B (zh) | 考虑高超声速化学非平衡效应的改进热壁修正方法 | |
| Schietekat et al. | Computational fluid dynamics‐based design of finned steam cracking reactors | |
| Tang et al. | Experimental investigation of a novel heat pipe thermoelectric generator for waste heat recovery and electricity generation | |
| Liu et al. | Study on safety boundary of flow instability and CHF for parallel channels in motion | |
| CN101961628B (zh) | 一种中小型绝热反应器 | |
| Gauthier et al. | An evaluation of numerical models for temperature-stabilized CH4/air flames in a small channel | |
| Bijjargi et al. | CFD modeling and optimal design of SiC deposition on the fuel combustion nozzle in a commercial CVD reactor | |
| CN102119293B (zh) | 设计极高温反应器用的同轴双管式热气管道的方法 | |
| CN111159875B (zh) | 一种电站凝汽器壳侧和管侧动态耦合数学模型及建模方法 | |
| Arkhazloo et al. | Determination of temperature distribution during heat treatment of forgings: Simulation and experiment | |
| CN116306375A (zh) | 高炉生产运行状态的数字孪生方法 | |
| Kuznetsov et al. | Numerical analysis of basic regularities of heat and mass transfer in a high-temperature heat pipe | |
| Govardhan et al. | Experimental investigations and CFD study of temperature distribution during oscillating combustion in a crucible furnace | |
| Liu et al. | Experimental study on transient heat transfer enhancement from a twisted plate in convection flow of helium gas | |
| Seto et al. | Heat transfer enhancement in sphere-packed pipes under high Reynolds number conditions | |
| Zhang et al. | Heat transfer analysis of the forced air quenching with non-isothermal and non-uniform oxidation | |
| Ahmed et al. | A Truncated transient slab model for a reheating furnace | |
| Dussoubs et al. | Modelling of a moving bed furnace for the production of uranium tetrafluoride Part 1: formulation of the model | |
| Chalaev et al. | Heat transfer enhancement in a corrugated tube heat exchanger | |
| Danov et al. | Exact formulations and nearly exact numerical solutions for convection in turbulent flow between parallel plates | |
| Hauptmann et al. | Axial development of unusual velocity profiles due to heat transfer in variable density fluids | |
| Osawa et al. | Study of Heattransfer During Piercing Process of Oxyfuel Gas Cutting | |
| Unz et al. | Calculation Approach for Ceramic Heat Pipe Heat Exchangers for High‐Temperature Applications | |
| CN113782106B (zh) | 一种用于cfd的高焓离解气体表面催化反应速率的求解方法 | |
| Mai et al. | Transient mixed convection flow instabilities in a vertical pipe |
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 |