CN105631100B - 水场景目标红外尾迹特性的流体模拟方法 - Google Patents
水场景目标红外尾迹特性的流体模拟方法 Download PDFInfo
- Publication number
- CN105631100B CN105631100B CN201510981178.3A CN201510981178A CN105631100B CN 105631100 B CN105631100 B CN 105631100B CN 201510981178 A CN201510981178 A CN 201510981178A CN 105631100 B CN105631100 B CN 105631100B
- Authority
- CN
- China
- Prior art keywords
- particle
- liquid
- infrared
- temperature
- fluid
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水场景目标红外尾迹特性的流体模拟方法。其步骤为:1)对场景中的模型进行体素化产生固体粒子;2)在存在液体的区域添加液体粒子;3)对每一时间帧,利用SPH对流体的力学特性进行数值模拟;4)对每一时间帧,利用瞬态热传导方程对流体的热力学特性进行数值模拟及红外特征图形绘制。本发明解决了流动介质中计算热力学特性的需求。本方法能够真实模拟水场景中的力学特征和红外特征,并在实验室环境下成功模拟了近年新发现的红外现象:船只在航行过程产生的红外尾迹。
Description
技术领域
本发明涉及流体模拟的数值计算领域与热力学的数值计算领域。具体涉及一种水场景目标红外尾迹特性的流体模拟方法。
背景技术
近年来,随着红外热像仪成像分辨率的提高,在遥感及航空拍摄中,发现了船只航行时的红外新特征:船只在水场景中航行时会留下红外尾迹。然而现有的热力学计算方法普遍针对于静态物质进行数值模拟计算,并不能很好的模拟这一特征;而流体模拟的数值计算领域的方法也主要针对于计算模拟流体的力学特征,无法模拟热力学的热传导特征。因此,难以模拟计算并绘制出该红外特征现象。
下面介绍已有的流体模拟的数值计算方法与热力学的数值计算方法:
1)热力学的数值计算方法
目前在国内外的瞬态热传导问题的数值方法,主要基于网格的有限元法、有限差分法和有限体积法。近年来也提出了基于粒子的无网格方法。然而研究过程中,都没有考虑液体运动对其造成的影响,因此也无法模拟水体目标运动造成的红外现象。
2)流体模拟的数值计算方法
这类方法主要可以细分成为三类:欧拉法、拉格朗日法和混合法。欧拉法是一种基于网格的方法,它使用2D或者3D的网格,从流体空间中的各个固定点着手,分析每个固定点上的流体速度、压强、密度等参数随时间和空间的变化。拉格朗日法是一种基于粒子的方法,用一系列遵循物理规则的粒子来表示流体,把单个流体粒子作为研究对象,研究其运动要素(位置、速度等)的变化过程,并通过综合各个流体粒子的运动来获得一定空间内所有流体质点的运动规律。欧拉法与拉格朗日法混合的主要思想是对三维水体用欧拉网格法表示,再用拉格朗日法进行小尺度的泡沫和浪花等效果的模拟。这些方法能够很好的模拟水体的力学性质,并实现良好的可见光水体模拟效果,但无法得到热力学性质,因此无法模拟红外特征现象。
发明内容
本发明的目的在模拟水场景目标在水中的力学特征和红外特征,解决现有流体方法无法模拟红外特征和现有热力学方法无法模拟水体流动特征的问题,提供一种水场景中目标红外尾迹特性的流体模拟方法。
一种水场景目标红外尾迹特性的流体模拟方法包括以下步骤:
1)对场景中的模型进行体素化产生固体粒子;
2)在存在液体的区域添加液体粒子;
3)对每一时间帧,利用SPH对流体的力学特性进行数值模拟;
4)对每一时间帧,利用瞬态热传导方程对流体的热力学特性进行数值模拟及红外特征图形绘制。
所述步骤1)为:
(2.1)根据模型顶点位置数据计算模型的AABB包围盒Cb;
(2.2)将Cb平均划分为X*Y*Z个子正方体Ci,j,k
其中,X、Y、Z为空间在xyz三个方向上的分辨率,i∈[0,X),j∈[0,Y),k∈[0,Z);
(2.3)对每个Ci,j,k与模型进行求交,若相交则该区域存在固体体素,否则该区域固体体素为空,每个求交计算所得的固体体素标记为Si;
(2.4)根据模型的运动状态,设定Si每一时刻的速度值。
所述步骤2)为:
(3.1)设定待模拟的水场景液体表面曲面方程为F(p),其中p∈(R,R,R),R为有理数集合,F(p)<0表示点p在液体内部,F(p)=0表示点p在液体边界,F(p)>0表示点p在液体外部,
(3.2)将空间平均划分为X`*Y`*Z`个子正方体C`i,j,k
其中,X`、Y`、Z`为空间在xyz三个方向上的液体分辨率,i∈[0,X`),j∈[0,Y`),k∈[0,Z`);
(3.3)对每个C`i,j,k,计算F(i,j,k),若小于等于0则该区域存在液体粒子,若大于0则该区域不存在液体粒子,每个液体粒子标记为Li;
(3.4)根据水场景的运动特性与水体温度分布初始化Li的力学参数与热力学参数:
(a)根据初始化条件,初始化每个Li的速度;
(b)根据温度初始化方程T(h)=T0–ta*h,初始化每个Li的初始温度;
其中,h为Li的水深,T0为液体表面温度,ta为温度随深度的变化参数。
所述步骤3)为:
(4.1)对于每一个粒子Li,根据SPH方程定义,利用所有光滑核半径内的粒子Lj或Sj计算液体的密度ρi:
ρi=ρ(ri)=315m/(64πh9)∑j(h2-|ri-rj|2)2
其中ri是Li在空间中的位置,m是Li的质量,h为光滑核半径,rj为粒子Lj或Sj的空间位置;
(4.2)对于每一个粒子Li,根据理想气体状态方程计算液体的压强:
piM=ρiRTi
其中ρi,Ti,pi分别为Li的密度、温度、压强,M为平均摩尔质量,R为理想气体常数;
(4.3)根据SPH公式,利用所有光滑核半径内的粒子Lj或Sj计算每一个粒子Li的加速度ai:
ai=a(ri)=g+m*45/(πh6)∑j((pi+pj)/(2ρiρj)*(h-r)2*(ri-rj)/r)+mμ
*45/(πh6)∑j(uj-ui)/(ρiρj)*(h-r)
其中,r=|ri-rj|,ρj是粒子Lj或Sj的密度,aj为粒子Lj或Sj的加速度,uj为粒子Lj或Sj的速度,ai,ui分别为粒子Li的加速度和速度;
(4.4)利用加速度与速度,根据牛顿第二定律计算模拟每个粒子Li或Si的运动轨迹,重复步骤(4.1)~(4.4)即可计算得到每一时刻每个粒子的力学特性。
所述步骤4)为:
(5.1)对于每个粒子Li,根据笛卡尔坐标系下的热传导方程,利用狄利克雷边界条件计算以下方程得到温度变化:
其中W为光滑粒子核函数W(x),Wij=W(|ri-rj|),q,T为积分点的热流密度与温度;
(5.2)根据计算得到的温度场T的分布,利用Planck公式计算以下方程得到其红外辐射量E:
其中ε0为材料的表面发射率;C1为第一辐射常数,其值为3.742×10-16W·m2;C2为第二辐射常数,其值为1.4388×10-2m·K。λ1和λ2为红外探测仪的探测波段的两段波长;
(5.3)若场景中最高和最低辐射温度分别为Tmax和Tmin,由上式算得的相应的辐射强度分别为Emax和Emin,则对于辐射强度为Ei的表面,其对应的计算灰度值为:
利用计算机图形学中的颜色均衡方法,绘制得到红外特征图形效果。
本发明的有点在于:
传统的流体模拟方法,其基于流体力学分析,从而求得了流体的力学性质,但无法模拟热力学性质产生的红外现象;传统的热力学模拟方法,基于热传导方程,可计算物质的瞬态热力学性质,从而模拟红外现象,但无法模拟流体的运动现象。
本发明的方法提出了水场景目标红外尾迹特性的流体模拟方法。通过体素化将场景划分为粒子,利用SPH估值近似求解偏微分方程,同时求解流体的力学性质与热力学性质。这样我们就产生了可以同时具有力学特征和红外特征的水场景。
本方法提出了一种使用体素化统一划分场景的方法,来使得水场景中的液体和固体能够使用统一的方程进行物理计算,避免了复杂的边界条件计算。
本方法提出了一种使用SPH方法同时计算流体的力学性质与热力学性质的方法。同时模拟的液体的力学性质与热力学性质,能同时模拟水场景的红外特征。
总之,本发明给出了一种水场景目标红外尾迹特性的流体模拟方法,该方法对比目前方法,能够同时计算水场景的力学性质与热力学性质,可以模拟绘制出其尾迹的红外图形特征效果。
附图说明
图1为船只在海面航行时的液体分布图;
图2为船只在海面航行时的红外特征俯视图形。
具体实施方式
水场景目标红外尾迹特性的流体模拟方法包括以下步骤:
1)对场景中的模型进行体素化产生固体粒子;
2)在存在液体的区域添加液体粒子;
3)对每一时间帧,利用SPH对流体的力学特性进行数值模拟;
4)对每一时间帧,利用瞬态热传导方程对流体的热力学特性进行数值模拟及红外特征图形绘制;
所述步骤1)为:
(1.1)根据模型顶点位置数据计算模型的AABB包围盒Cb;
(1.2)将Cb平均划分为X*Y*Z个子正方体Ci,j,k
其中,X、Y、Z为空间在xyz三个方向上的分辨率,i∈[0,X),j∈[0,Y),k∈[0,Z);
(1.3)对每个Ci,j,k求与模型进行求交,若相交则该区域存在固体体素,否则该区域固体体素为空,每个计算所得的固体体素标记为Si。
(1.4)根据模型的运动状态,设定Si每一时刻的速度值;
所述步骤2)为:
(2.1)使用方程F(p)表达液体表面方程
其中p∈(R,R,R),F(p)<0表示点p在液体内部,F(p)=0表示点p在液体边界,F(p)>0表示点p在液体外部。F(p)的方程形式可根据实际实验需要选取;
(2.2)将空间平均划分为X`*Y`*Z`个子正方体C`i,j,k
其中,X`、Y`、Z`为空间在xyz三个方向上的液体分辨率,i∈[0,X`),j∈[0,Y`),k∈[0,Z`);
(2.3)对每个Ci,j,k,计算F(i,j,k),若大于0则该区域存在液体粒子,若小于0则该区域不存在液体粒子,每个液体粒子标记为Li;
(2.4)根据水场景的运动特性与水体温度分布初始化Li的力学参数与热力学参数:
(a)根据初始化条件,初始化每个Li的速度。若无特殊要求,Li的速度初始化为0(水面静止);
(b)根据温度初始化方程T(h)=T0–ta*h,初始化每个Li的初始温度;
其中,h为Li的水深,T0为液体表面温度,ta为温度随深度的变化参数,T0与ta根据实际实验需要选取。
所述步骤3)为:
(3.1)对于每一个粒子Li,根据SPH方程定义,利用所有光滑核半径内的粒子Lj或Sj计算液体的密度:
ρi=ρ(ri)=315m/(64πh9)∑j(h2-|ri-rj|2)2
其中ri是Li在空间中的位置,m是Li的质量,h为光滑核半径,rj为粒子Lj或Sj的空间位置;
(3.2)对于每一个粒子Li,根据理想气体状态方程计算液体的压强:
piM=ρiRTi
其中ρi,Ti,pi分别为Li的密度、温度、压强,M为平均摩尔质量,R为比例常数,约为8.31441±0.00026J/(mol·K);
(3.3)根据SPH公式,利用所有光滑核半径内的粒子Lj或Sj计算每一个粒子Li的加速度:
ai=a(ri)=g+m*45/(πh6)∑j((pi+pj)/(2ρiρj)*(h-r)2*(ri-rj)/r)+mμ
*45/(πh6)∑j(uj-ui)/(ρiρj)*(h-r)
其中,r=|ri-rj|,ai,ui分别为粒子Li在的加速度,速度,aj,uj分别为粒子Lj或Sj的加速度或速度;
(3.4)利用加速度与速度,根据牛顿第二定律计算模拟每个粒子Li或Si的运动轨迹,重复步骤(1)~(4)即可计算得到每一时刻每个粒子的力学特性。海浪形状稳定后的船只在海面航行时的液体分布图如图1所示;
所述步骤4)为:
(4.1)对于每个粒子Li,根据笛卡尔坐标系下的热传导方程,利用狄利克雷边界条件计算以下方程得到温度变化:
其中W为光滑粒子核函数W(x),Wij=W(|ri-rj|),可根据实验需要选用,本利中的计算结果采用lucy核函数进行数值计算。q,T为积分点的热流密度与温度。
(4.2)根据计算得到的温度场T的分布,利用Planck公式计算以下方程得到其红外辐射量E:
其中ε0为材料的表面发射率;C1为第一辐射常数,其值为3.742×10-16W·m2;C2为第二辐射常数,其值为1.4388×10-2m·K。λ1和λ2为红外探测仪的探测波段的两段波长。
(4.3)若场景中最高和最低辐射温度分别为Tmax和Tmin,由上式算得的相应的辐射强度分别为Emax和Emin(ε0的值为1)。则对于辐射强度为Ei的表面,其对应的计算灰度值为:
利用计算机图形学中的颜色均衡方法,本实施例中采用Gouraud Shading明暗处理方法,可绘制得到红外特征图形结果,结果如图2所示。
Claims (3)
1.一种水场景目标红外尾迹特性的流体模拟方法,其特征在于包括以下步骤:
1)对场景中的模型进行体素化产生固体粒子;
2)在存在液体的区域添加液体粒子;
3)对每一时间帧,利用SPH对流体的力学特性进行数值模拟;
所述步骤3)为:
(3.1)对于每一个液体粒子Li,根据SPH方程定义,利用所有光滑核半径内的液体粒子Lj或固体粒子Sj计算液体的密度ρi:
ρi=ρ(ri)=315m/(64πh9)∑j(h2-|ri-rj|2)2
其中ri是Li在空间中的位置,m是Li的质量,h为光滑核半径,rj为粒子Lj或Sj的空间位置;
(3.2)对于每一个粒子Li,根据理想气体状态方程计算液体的压强:
piM=ρiRTi
其中ρi,Ti,pi分别为Li的密度、温度、压强,M为平均摩尔质量,R为理想气体常数;
(3.3)根据SPH公式,利用所有光滑核半径内的粒子Lj或Sj计算每一个粒子Li的加速度ai:
ai=a(ri)=g+m*45/(πh6)∑j((pi+pj)/(2ρiρj)*(h-r)2*(ri-rj)/r)+
mμ*45/(πh6)∑j(uj-ui)/(ρiρj)*(h-r)
其中,r=|ri-rj|,ρj是粒子Lj或Sj的密度,uj为粒子Lj或Sj的速度,ai,ui分别为粒子Li的加速度和速度;
(3.4)利用加速度与速度,根据牛顿第二定律计算模拟每个粒子Li或Si的运动轨迹,重复步骤(3.1)~(3.4)即可计算得到每一时刻每个粒子的力学特性;
4)对每一时间帧,利用瞬态热传导方程对流体的热力学特性进行数值模拟及红外特征图形绘制;
所述步骤4)为:
(4.1)对于每个粒子Li,根据笛卡尔坐标系下的热传导方程,利用狄利克雷边界条件计算以下方程得到温度变化:
其中W为光滑粒子核函数W(x),Wij=W(|ri-rj|),q,T为积分点的热流密度与温度;
(4.2)根据计算得到的温度场T的分布,利用Planck公式计算以下方程得到其红外辐射量E:
其中ε0为材料的表面发射率;C1为第一辐射常数,其值为3.742×10-16W·m2;C2为第二辐射常数,其值为1.4388×10-2m·K,λ1和λ2为红外探测仪的探测波段的两段波长;
(4.3)若场景中最高和最低辐射温度分别为Tmax和Tmin,由上式算得的相应的辐射强度分别为Emax和Emin,则对于辐射强度为Ei的表面,其对应的计算灰度值为:
利用计算机图形学中的颜色均衡方法,绘制得到红外特征图形效果。
2.根据权利要求1所述的一种水场景目标红外尾迹特性的流体模拟方法,其特征在于所述步骤1)为:
(1.1)根据模型顶点位置数据计算模型的AABB包围盒Cb;
(1.2)将Cb平均划分为X*Y*Z个子正方体Ci,j,k;
其中,X、Y、Z为空间在xyz三个方向上的分辨率,i∈[0,X),j∈[0,Y),k∈[0,Z);
(1.3)对每个Ci,j,k与模型进行求交,若相交则该区域存在固体体素,否则该区域固体体素为空,每个求交计算所得的固体体素标记为Si;
(1.4)根据模型的运动状态,设定Si每一时刻的速度值。
3.根据权利要求1所述的一种水场景目标红外尾迹特性的流体模拟方法,其特征在于所述步骤2)为:
(2.1)设定待模拟的水场景液体表面曲面方程为F(p),其中p∈(R,R,R),R为有理数集合,F(p)<0表示点p在液体内部,F(p)=0表示点p在液体边界,F(p)>0表示点p在液体外部,
(2.2)将空间平均划分为X`*Y`*Z`个子正方体C`i,j,k;
其中,X`、Y`、Z`为空间在xyz三个方向上的液体分辨率,i∈[0,X`),j∈[0,Y`),k∈[0,Z`);
(2.3)对每个C`i,j,k,计算F(i,j,k),若小于等于0则该区域存在液体粒子,若大于0则该区域不存在液体粒子,每个液体粒子标记为Li;
(2.4)根据水场景的运动特性与水体温度分布初始化Li的力学参数与热力学参数:
(a)根据初始化条件,初始化每个Li的速度;
(b)根据温度初始化方程T(h)=T0–ta*h,初始化每个Li的初始温度;
其中,h为Li的水深,T0为液体表面温度,ta为温度随深度的变化参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510981178.3A CN105631100B (zh) | 2015-12-23 | 2015-12-23 | 水场景目标红外尾迹特性的流体模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510981178.3A CN105631100B (zh) | 2015-12-23 | 2015-12-23 | 水场景目标红外尾迹特性的流体模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631100A CN105631100A (zh) | 2016-06-01 |
CN105631100B true CN105631100B (zh) | 2018-10-26 |
Family
ID=56046029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510981178.3A Active CN105631100B (zh) | 2015-12-23 | 2015-12-23 | 水场景目标红外尾迹特性的流体模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631100B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096215B (zh) * | 2016-07-28 | 2019-04-30 | 华东师范大学 | 一种涉及热传导与动态黏度的真实感流体仿真方法 |
CN106951599A (zh) * | 2017-02-27 | 2017-07-14 | 南京航空航天大学 | 一种飞机凝结尾生成对全球地表温度影响的评估方法 |
CN109237999B (zh) * | 2018-09-19 | 2020-09-08 | 中国电子科技集团公司第二十八研究所 | 一种批量三维态势目标尾迹实时绘制方法和系统 |
CN115719012B (zh) * | 2023-01-06 | 2023-04-14 | 山东科技大学 | 基于无人机遥感与多相sph算法的尾矿库放矿布置方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN104574518A (zh) * | 2014-12-26 | 2015-04-29 | 浙江大学 | 具有丰富细节的大规模海洋场景半物理绘制方法 |
-
2015
- 2015-12-23 CN CN201510981178.3A patent/CN105631100B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN104574518A (zh) * | 2014-12-26 | 2015-04-29 | 浙江大学 | 具有丰富细节的大规模海洋场景半物理绘制方法 |
Non-Patent Citations (2)
Title |
---|
GPU中的流体场景实时模拟算法;陈曦等;《计算机辅助设计与图形学学报》;20100331;第396-405页 * |
Real-time fluid simulation with adaptive SPH;He Yan等;《Computer Animation and Virtual Worlds》;20090630;第417-426页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105631100A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772682B (zh) | 一种动目标的红外辐射光谱特性仿真分析方法 | |
Windsor et al. | The flow fields involved in hydrodynamic imaging by blind Mexican cave fish (Astyanax fasciatus). Part I: open water and heading towards a wall | |
CN105631100B (zh) | 水场景目标红外尾迹特性的流体模拟方法 | |
Takizawa et al. | Computer modeling techniques for flapping-wing aerodynamics of a locust | |
Harris | Real-time cloud simulation and rendering | |
Reid et al. | Fluid dynamics of moving fish in a two-dimensional multiparticle collision dynamics model | |
CN114168796B (zh) | 一种建立飞行器高空气动力数据库的方法 | |
CN109658496B (zh) | 一种飞机红外纹理图像生成方法 | |
Ferragut et al. | Comparison between 2.5-D and 3-D realistic models for wind field adjustment | |
CN115994496B (zh) | 城市公园高分辨率大气co2浓度三维场的数值模拟方法 | |
US10282885B2 (en) | Computer graphic system and method of multi-scale simulation of smoke | |
Duarte et al. | Real-time simulation of cumulus clouds through skewt/logp diagrams | |
Goswami et al. | Real-time landscape-size convective clouds simulation and rendering | |
CN112347708B (zh) | 一种基于物理过程的实时光学烟幕的渲染方法 | |
Qiu et al. | Dispersion simulation and visualization for urban security | |
Serón et al. | The evolution of a wildland forest fire front | |
CN110992399B (zh) | 一种高精度目标大气扰动检出方法 | |
Tagliavini et al. | Wake characteristics of complex-shaped snow particles: Comparison of numerical simulations with fixed snowflakes to time-resolved particle tracking velocimetry experiments with free-falling analogs | |
Jang et al. | Simulating sph fluid with multi-level vorticity | |
Levis et al. | An Efficient Approach for Optical Radiative Transfer Tomography using the Spherical Harmonics Discrete Ordinates Method | |
Jiang et al. | Measurement of the movement parameters of saltating sand over a flat sand bed using a high-speed digital camera | |
Bao et al. | A Fast Calculation and Rendering Method for Infrared Characteristics of Hypersonic Vehicle | |
Tan et al. | Real-time snowing simulation based on particle systems | |
Valentin et al. | Validation of a new solver based on the vortex particle method for wings, propellers and rotors | |
CN115438393A (zh) | 一种混合域红外辐射特性的仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |