CN108090030B - 一种圆形单桩局部流场的处理方法 - Google Patents
一种圆形单桩局部流场的处理方法 Download PDFInfo
- Publication number
- CN108090030B CN108090030B CN201711414988.6A CN201711414988A CN108090030B CN 108090030 B CN108090030 B CN 108090030B CN 201711414988 A CN201711414988 A CN 201711414988A CN 108090030 B CN108090030 B CN 108090030B
- Authority
- CN
- China
- Prior art keywords
- flow
- pile foundation
- water flow
- local
- mathematical model
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种圆形单桩局部流场的处理方法,包括建立当量尺寸圆柱体平面二维水流数学模型;将雷诺数分级;利用二维水流数学模型计算不同代表雷诺数条件下圆柱体周边的流场;计算确定不同代表雷诺数条件下当量尺寸圆柱体周边任意点绕流场系数;根据需要建立实际工程海域整体平面二维水流数学模型;利用实际二维水流数学模型,不考虑桩基工程,对研究时段内的流场进行模拟计算;提取桩基局部的水流时间过程作为来流条件,结合桩基尺寸,计算并选择代表雷诺数,查找计算的绕流场系数;计算桩基局部的流场;将计算结果代入整体数学模型,替换桩基影响范围内整体模型计算结果。可以节约计算耗时,从而满足河口海岸地区计算时段较长的要求。
Description
技术领域
本发明涉及一种圆形单桩局部流场的处理方法,涉及一种潮流模拟中的桩基概化处理方法,适用于海岸地区的桩基附近的二维水流模拟,属于海岸动力数值模拟领域
背景技术
随着计算机和数值模拟技术的发展,二维水流数值模拟在海岸工程领域得到广泛应用。非结构化网格由于能够准确的拟合岛屿或不规则的岸线边界,而广受研究人员的喜好。采用非结构网格的数值计算方法,目前经常被使用的主要有有限体积法、有限单元法等。在河口海岸地区的水流数值计算中常常遇到海工建筑物的问题,其中桥墩和栈桥、高桩码头的桩基部分处理由于尺度较小,处理相对困难。对于桩基的处理,目前主要存在两种方法:一、加密网格直接模拟。利用非结构网格在桩基局部进行加密,实现刻画出桩基固壁边界,通过有限体积法或有限单元法在模型里进行模拟计算。该方法精度相对较高,其弊端是由于大部分桩基尺度一般在1m左右或更小,直接模拟的网格尺度需要在0.1m左右,远小于海岸地区正常的网格尺度,局部的加密带来了网格数量的急剧增多和时间步长的缩短,从而导致计算耗时大大增加。二、附加阻力概化方法。该方法不直接模拟桩基的边界影响,而采用桩基单元附加阻力,以反映桩基的阻力效应;并对地形进行修正,以反映桩基对过水面积的影响。这种模式不必描述桩基的外形,可以加大计算网格尺度、缩短计算时间。该方法的不足是,将桩基作为完全过水区域,会产生很大的概化误差,体现在流场形态和分流比的严重失真。
严格意义上,桩基局部的水流属于复杂的三维湍流问题,需要通过N-S方程直接模拟、三维紊流模型(k-ε、k-ω等)、大涡模拟方法、离散涡等方法在细密网格或计算点下进行局部的精细模拟,然而这些方法计算过程相对复杂、最主要的是计算量巨大,一般仅在很局部范围内进行模拟,而在河口海岸工程尺度的模型研究范围内难以满足实际应用需要。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种圆形单桩局部流场的处理方法,将局部精细化模拟结果与大范围水流数值模拟方法进行结合,采用一种概化的方法确定桩基局部流场,而其它范围内仍然采用大范围的数值计算结果,从而达到成能够适应较大范围流场计算要求又能体现桩基局部细化的水流分布的目的。同时,由于不需要在模型内对桩基局部网格加密进行直接模拟,可以节约计算耗时,从而满足河口海岸地区计算时段较长的要求。
为了实现上述目标,本发明采用如下的技术方案:一种圆形单桩局部流场的处理方法,其特征在于:包括如下步骤:
步骤1:建立当量尺寸圆柱体平面二维水流数学模型,记录模型范围和网格点坐标;
步骤2:将雷诺数从小到大分级,确定每一分级的代表雷诺数;
步骤3:利用步骤1的模型计算不同代表雷诺数条件下圆柱体周边的流场;
步骤4:计算确定不同代表雷诺数条件下当量尺寸圆柱体周边流场任意点绕流场系数;
步骤5:根据需要建立实际工程海域整体平面二维水流数学模型;
步骤6:利用步骤5所建立数学模型,不考虑桩基工程,模拟研究工程海域不同时间流场过程,对研究时段内的流场进行模拟计算;
步骤7:根据步骤6计算得到的结果,提取桩基局部的水流时间过程作为来流条件;
步骤8:根据步骤7确定的来流条件,结合桩基尺寸,计算不同时间段雷诺数;
步骤9:根据步骤8计算得到的雷诺数,根据步骤2中雷诺数分级选择恰当的代表雷诺数;
步骤10:根据步骤9选择的代表雷诺数,查找步骤4计算的绕流场系数;
步骤11:根据对应的任意点水流与来流关系式计算桩基局部的流场;
步骤12:重复步骤10,计算工程海域整体数学模型中每一步时间过程中的桩基局部流场;
步骤13:将步骤11计算结果代入整体数学模型,替换桩基影响范围内整体模型计算结果。
进一步的,步骤2中,选择河口海岸地区常规的水流雷诺数,结合流体力学绕流特性:当Re<5时,绕流不发生分离;当5<Re<150,涡街层流;当150<Re<300,旋涡由层流向湍流转变;当300<Re<3×105,为亚临界区;当3×105<Re<3×106,为过渡区;当Re>3×106,称为超临界区。可以根据分区特征从小到大进行分级。
进一步的,在每一段分级中,取中间值作为该分段的代表雷诺数。
进一步的,步骤4中,将步骤3计算得到的每个流速点流速值与来流流速进行比较,计算关系式如下:
其中,Cx、Cy为绕流场系数,u、v为流速分量,V0为来流流速;
得到每个流速点的绕流场系数(Cx,Cy)。每个流速点坐标采用与桩基直径的相对值,如采用球坐标系,点坐标表示为(r(d),θ)。储存计算结果。
进一步的,步骤6中,对步骤5所建立的模型给定边界条件和模型参数,利用数值方法,不考虑桩基工程,对研究时段内的流场进行模拟计算。
进一步的,数值方法为有限体积法或有限单元法。
进一步的,步骤7中,选择工程海域整体数学模型中实际桩基的位置处的流速过程,作为桩基绕流的来流流速。
进一步的,步骤11中,根据步骤1的模型的计算范围与当量直径的关系,在大模型里确定同样的范围,将坐标转换为以桩基中心为原点,来流速度方向为x轴,垂直方向为y轴的局部坐标系。
进一步的,步骤13中,忽略原模型中可能存在的非恒定向影响,忽略桩基局部地形变化对水流的影响,在局部模型范围与大模型结合的区域采用线性插值进行平滑过渡。
进一步的,步骤1、步骤2、步骤3、步骤4中计算得到的绕流场系数,可作为参照表格,重复使用。
本发明的有益之处在于:将局部精细化模拟结果与工程海域整体范围水流数值模拟方法进行结合,采用一种概化的方法确定桩基局部流场,其它范围内仍然采用整体范围的数值计算结果,从而达到成能够适应较大范围流场计算要求又能体现桩基局部细化的水流分布的目的。同时,由于不需要在模型内对桩基局部网格加密进行直接模拟,可以节约计算耗时,从而满足河口海岸地区计算时段较长的要求。
附图说明
图1是本发明一实施例圆形单桩局部流场的处理方法的流程图;
图2是本发明一实施例计算流程示意图;
图3是本发明一实施例浅水方程模拟的圆柱体周边近似流场(Re=1.0×107)示意图;
图4是本发明一实施例浅水方程模拟的圆柱体周边流速比分布(Re=1.0×107)示意图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
参照图1圆形单桩局部流场的处理方法的流程图;
参照图2计算流程示意图;
参照图3浅水方程模拟的圆柱体周边近似流场(Re=1.0×107);
参照图4浅水方程模拟的圆柱体周边流速比分布(Re=1.0×107);
本发明一种圆形单桩局部流场的处理方法,包括如下步骤:
S1:建立当量尺寸圆柱体平面二维水流数学模型,记录模型范围和网格点坐标。
优选的,步骤S1中,根据圆柱体绕流的特性,确定圆柱体的当量直径d(可视为1),选择合适的模型范围(x>20d,y>10d)、模型网格尺度(<0.1d),选择合适的数值方法(非结构网格有限体积法、有限元法等),选择常规的模型参数(底部糙率、水流紊动粘性系数),建立圆柱体绕流局部的平面二维水流数学模型。在二维模型中,对模型水深、底部糙率进行概化,选择常规的水深(5m或10m或其它)和常规的糙率(如曼宁系数0.02),忽略这种概化对流场的影响。
优选的,步骤S1中,模型采用局部坐标系,x轴与流速方向平行,y轴与流速方向垂直,原点位于圆柱体中心。
S2:根据流体力学绕流特性对雷诺数进行分级,并确定每一级中的代表雷诺数。
优选的,步骤S2中,选择河口海岸地区常规的水流雷诺数,结合流体力学绕流特性(Re<5时,绕流不发生分离;5<Re<150,涡街层流;150<Re<300,旋涡由层流向湍流转变;300<Re<3×105,为亚临界区;3×105<Re<3×106,为过渡区;Re>3×106,称为超临界区)从小到大进行分级;在每一段分级中,取中间值作为该分段的代表雷诺数
S3:利用步骤S1的模型计算不同代表雷诺数条件下圆柱体周边的流场。
优选的,步骤S3中,根据每一个代表雷诺数,通过雷诺数Re表达式反算来流流速V0;
ρ为水的密度,V0为来流流速,d为桩基直径,μ为水粘滞系数。V0作为入流边界条件,出流边界条件设置成n为出流边界法线方向。采用常规数值方法(如有限体积法、有限元法)计算不同代表雷诺数条件下的流场。采用的平面二维水流数学模型,水平扩散系数采用Smagorinsky公式,近似模拟桩基局部的紊流流场。
其中A水平扩散系数,为cs为常数(一般取0.28),l为特征长度,u、v为x,y方向流速分量。
S4:计算不同代表雷诺数条件下当量尺寸圆柱体周边流场任意点绕流场系数;
优选的,步骤S4中,将步骤3计算得到的每个流速点流速值与来流流速进行比较,计算关系式如下:
其中,Cx、Cy为绕流场系数,u、v为x,y方向流速分量,V0为来流流速;
得到每个流速点的绕流场系数(Cx,Cy)。每个流速点坐标采用与桩基直径的相对值,如采用球坐标系,点坐标表示为(r(d),θ)。储存计算结果。
S5:根据需要建立实际工程海域整体平面二维水流数学模型。
优选的,步骤S5中,数学模型中可对桩基局部进行适当加密。
S6:利用步骤S5所建立的平面二维水流数学模型,不考虑桩基工程,模拟研究海域不同时间流场过程,对研究时段内的流场进行模拟计算。
优选的,步骤S6中,对步骤S5所建立的模型给定边界条件和模型参数,利用数值方法(有限体积法、有限单元法),不考虑桩基工程,对研究时段内的流场进行模拟计算。
S7:根据步骤S6计算得到的流场结果,提取桩基局部的水流时间过程作为来流条件。
优选的,步骤S7中,选择数学模型中实际桩基的位置处的流速过程,作为桩基绕流的来流流速。
S8:根据步骤S7确定的来流条件,结合桩基尺寸,计算不同时间段雷诺数;
S9:根据步骤S8计算得到的不同时间段雷诺数,选择对应的分级,从分级中确定代表雷诺数;
S10:根据步骤S9选择的代表雷诺数,查找步骤S4计算的绕流场系数;
S11:根据对应的任意点水流与来流关系式计算桩基局部的流场。
优选的,步骤S11中,根据步骤S1的模型的计算范围与当量直径的关系,在大模型里确定同样的范围,将坐标转换为以桩基中心为原点,来流速度方向为x轴,垂直方向为y轴的局部坐标系。
S12:重复步骤S10,计算实际工程海域整体数学模型中每一步时间过程中的桩基局部流场;
S13:将步骤S11计算结果代入整体数学模型,替换桩基影响范围内整体数学计算结果。
优选的,步骤S13中,忽略原模型中可能存在的非恒定向影响,忽略桩基局部地形变化对水流的影响。在局部模型范围与大模型结合的区域采用线性插值进行平滑过渡。
优选的,步骤S1、步骤S2、步骤S3、步骤S4中计算得到的绕流场系数,可作为参照表格,重复使用。
本发明将局部精细化模拟结果与大范围水流数值模拟方法进行结合,采用桩基局部,采用一种概化的方法确定流场,而其它范围内仍然采用大范围的数值计算结果,从而达到成能够适应较大范围流场计算要求又能体现桩基局部细化的水流分布的目的,同时,由于不需要在模型内对桩基局部网格加密进行直接模拟,可以节约计算耗时,从而满足河口海岸地区计算时段较长的要求。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (7)
1.一种圆形单桩局部流场的处理方法,其特征在于:包括如下步骤:
步骤1:建立当量尺寸圆柱体平面二维水流数学模型,记录模型范围和网格点坐标;
步骤2:将雷诺数从小到大分级,确定每一分级的代表雷诺数;
步骤3:利用步骤1的模型计算不同代表雷诺数条件下圆柱体周边的流场;
步骤4:计算确定不同代表雷诺数条件下当量尺寸圆柱体周边流场任意点绕流场系数;步骤4中,将步骤3计算得到的每个流速点流速值与来流流速进行比较,计算关系式如下:
其中,Cx、Cy为绕流场系数,u、v为绕流场流速分量,V0为来流流速;
得到每个流速点的绕流场系数(Cx、Cy),每个流速点坐标采用与桩基直径的相对值,储存计算结果;
步骤5:根据需要建立工程海域整体平面二维水流数学模型;
步骤6:利用步骤5所建立工程海域整体平面二维水流数学模型,不考虑桩基工程,对研究时段内的流场进行模拟计算;步骤6中,对步骤5所建立的工程海域整体平面二维水流数学模型给定边界条件和模型参数,利用数值方法,不考虑桩基工程,对研究时段内的流场进行模拟计算;
步骤7:根据步骤6计算得到的结果,提取桩基局部的水流时间过程作为来流条件;
步骤8:根据步骤7确定的来流条件,结合桩基尺寸,计算不同时间段雷诺数;
步骤9:根据步骤8计算得到的雷诺数,根据步骤二中雷诺数分级选择恰当的代表雷诺数;
步骤10:根据步骤9选择的代表雷诺数,查找步骤4计算的绕流场系数;
步骤11:根据对应的任意点水流与来流关系式计算桩基局部的流场;
步骤12:重复步骤10,计算工程海域整体平面二维水流数学模型中每一步时间过程中的桩基局部流场;
步骤13:将步骤11计算结果代入工程海域整体平面二维水流数学模型,替换桩基影响范围内整体数学模型计算结果。
2.根据权利要求1所述的一种圆形单桩局部流场的处理方法,其特征在于:步骤2中,选择河口海岸地区常规的水流雷诺数,结合流体力学绕流特性:当Re<5时,绕流不发生分离;当5<Re<150,涡街层流;当150<Re<300,旋涡由层流向湍流转变;当300<Re<3×105,为亚临界区;当3×105<Re<3×106,为过渡区;当Re>3×106,称为超临界区来从小到大进行分级。
3.根据权利要求2所述的一种圆形单桩局部流场的处理方法,其特征在于:在每一段分级中,取中间值作为该分段的代表雷诺数。
4.根据权利要求1所述的一种圆形单桩局部流场的处理方法,其特征在于:数值方法为有限体积法或有限单元法。
5.根据权利要求1所述的一种圆形单桩局部流场的处理方法,其特征在于:步骤7中,选择工程海域整体平面二维水流数学模型中实际桩基的位置处的流速过程,作为桩基绕流的来流流速。
6.根据权利要求1所述的一种圆形单桩局部流场的处理方法,其特征在于:步骤11中,根据步骤1的模型的计算范围与当量直径的关系,在实际工程海域整体平面二维水流数学模型里确定同样的范围,将坐标转换为以桩基中心为原点,来流速度方向为x轴,垂直方向为y轴的局部坐标系。
7.根据权利要求1所述的一种圆形单桩局部流场的处理方法,其特征在于:步骤13中,忽略步骤5中建立的工程海域整体平面二维水流数学模型可能存在的非恒定向影响,忽略桩基局部地形变化对水流的影响,在局部模型范围与实际工程海域整体平面二维水流数学模型结合的区域采用线性插值进行平滑过渡。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711414988.6A CN108090030B (zh) | 2017-12-22 | 2017-12-22 | 一种圆形单桩局部流场的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711414988.6A CN108090030B (zh) | 2017-12-22 | 2017-12-22 | 一种圆形单桩局部流场的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108090030A CN108090030A (zh) | 2018-05-29 |
CN108090030B true CN108090030B (zh) | 2021-02-19 |
Family
ID=62178917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711414988.6A Active CN108090030B (zh) | 2017-12-22 | 2017-12-22 | 一种圆形单桩局部流场的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108090030B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188020B (zh) * | 2018-09-10 | 2020-09-25 | 中国长江三峡集团有限公司 | 一种基于尾流纹路识别的水面流速测量方法 |
CN115214807A (zh) * | 2022-08-30 | 2022-10-21 | 上海交通大学 | 一种具有气动减阻功能的基于仿生微结构的表面改性技术 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6516292B2 (en) * | 1999-02-10 | 2003-02-04 | Asher Yahalom | Method and system for numerical simulation of fluid flow |
CN105841921A (zh) * | 2016-03-24 | 2016-08-10 | 江苏理工学院 | 一种湍流流场多尺度的测量系统及其测量方法 |
CN107436977A (zh) * | 2017-07-24 | 2017-12-05 | 珠江水利委员会珠江水利科学研究院 | 复杂河网分流的数值模拟方法 |
-
2017
- 2017-12-22 CN CN201711414988.6A patent/CN108090030B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6516292B2 (en) * | 1999-02-10 | 2003-02-04 | Asher Yahalom | Method and system for numerical simulation of fluid flow |
CN105841921A (zh) * | 2016-03-24 | 2016-08-10 | 江苏理工学院 | 一种湍流流场多尺度的测量系统及其测量方法 |
CN107436977A (zh) * | 2017-07-24 | 2017-12-05 | 珠江水利委员会珠江水利科学研究院 | 复杂河网分流的数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
Numerical simulation of the flow around a circular cylinder at high Reynolds numbers;Pietro Catalano 等;《International Journal of Heat and Fluid Flow》;20030831;第24卷(第4期);第463-469页 * |
二维圆柱绕流数值模拟;杨烁 等;《中国造船》;20171130;第48卷;第533-540页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108090030A (zh) | 2018-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Davies | A model of the vertical structure of the wave and current bottom boundary layer | |
Xie et al. | Incorporation of a mass-conserving inundation scheme into a three dimensional storm surge model | |
Schiller et al. | Modeling river plume dynamics with the HYbrid Coordinate Ocean Model | |
CN108090030B (zh) | 一种圆形单桩局部流场的处理方法 | |
CN115221721A (zh) | 一种基于roms的区域海域潮汐潮流数值模拟方法 | |
Markus et al. | A numerical investigation of combined wave–current loads on tidal stream generators | |
Akbari et al. | Tidal constituents in the Persian Gulf, Gulf of Oman and Arabian Sea: a numerical study | |
CN110287504B (zh) | 一种模拟近岸海水交换规律的耦合模型的建立方法 | |
Avdis et al. | Meshing ocean domains for coastal engineering applications | |
CN110990926A (zh) | 一种基于面积修正率的城市地表建筑水动力学仿真方法 | |
Sayehbani et al. | Numerical modeling of wave and current patterns of beris port in east of chabahar-iran | |
Daud et al. | Hydrodynamic modelling of Bidong Island vicinity waters | |
Zhang et al. | Coupling of wave and current numerical model with unstructured quadtree grid for nearshore coastal waters | |
Fattah et al. | Hydrodynamic and sediment transport modelling of suralaya coastal area, Cilegon, Indonesia | |
Skandali | Identification of response amplitude operators for ships based on full scale measurements | |
CN114417752A (zh) | 一种河流入海口处水系的水位与流速模拟方法及系统 | |
Shi et al. | Spatial interpolation of tidal data using a multiple-order harmonic equation for unstructured grids | |
Shi et al. | Wind wave-forced fine sediment erosion during the slack water periods in Hangzhou Bay, China | |
Galperin et al. | Salinity intrusion and residual circulation in Delaware Bay during the drought of 1984 | |
Abdullah et al. | Hydrodynamic performance of cylindrical floating breakwater in waves | |
Aghajanloo et al. | Numerical simulation of tidal currents in Persian gulf | |
Zhang et al. | Effects of wave forces on sediment transport patterns in micro-tidal estuaries | |
Kim et al. | The role of diffraction effects in extreme run-up inundation at Okushiri Island due to 1993 tsunami | |
Roy et al. | MORPHOLOGICAL RESPONSES OF A TIDAL RIVER DUE TO CLIMATE CHANGE: A CASE STUDY FOR KARNAFULI RIVER, BANGLADESH | |
Jones et al. | Storm surge computations in estuarine and near-coastal regions: the Mersey estuary and Irish Sea area |
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 |