CN111881631A - 一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 - Google Patents
一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 Download PDFInfo
- Publication number
- CN111881631A CN111881631A CN202010713015.8A CN202010713015A CN111881631A CN 111881631 A CN111881631 A CN 111881631A CN 202010713015 A CN202010713015 A CN 202010713015A CN 111881631 A CN111881631 A CN 111881631A
- Authority
- CN
- China
- Prior art keywords
- oxygen
- lead
- phase
- oxygen concentration
- oxygen control
- 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.)
- Granted
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法,可用于铅冷快堆气相控氧和固相控氧过程中质量交换器内氧浓度扩散和分布特性数值模拟研究。通过对气相控氧和固相控氧的共同特性进行研究,将控氧过程中的多相流物理现象简化为单相模拟,可以基于一套方法同时完成气相和固相控氧的数值模拟。本发明具有资源占用低、结果可靠性高的特点,可用于铅冷快堆氧控实验预分析,有助于节省实验资源,获得更好的实验效果。
Description
技术领域
本发明涉及铅冷快堆计算流体力学模拟技术领域,具体涉及一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法。
背景技术
铅冷快堆是指以液态金属铅或铅秘合金(统称铅基冷却剂)为冷却剂的一类反应堆,属于第四代反应堆,具有优秀的中子经济性和安全性,在小型化和模块化方面有广阔发展前景。铅基冷却剂对结构材料的腐蚀作用是制约铅冷快堆发展的关键问题,目前针对腐蚀问题最具可行性的解决方法是通过控制铅基冷却剂中的氧浓度实现的。通过将冷却剂中的氧浓度控制在一定范围内,就可以使得结构材料表面生成一层保护性氧化膜,阻止冷却剂对结构材料的进一步腐蚀。同时,氧浓度的范围还必须保证低于氧化铅和氧化铋的溶解度,防止氧化铅和氧化铋析出,堵塞流动通道。因此,可以看出,氧控防腐蚀的实现关键在于对氧浓度的精确控制。
目前控氧技术主要分为两方面,气相控氧和固相控氧。气相控氧通过气体与铅基冷却剂的表面氧气质量交换使氧气溶解于冷却剂中,需要降低氧浓度时,可以通过氢气与铅基冷却剂的反应将多余氧气带走,即可实现氧浓度的控制。气相控氧根据使用气体不同主要分为氢气/氧气控氧和氢气/水控氧两种。固相控氧中,控制铅基冷却剂流经氧化铅质量交换器,通过调节通过的铅基冷却剂温度和流量即可改变质量交换器出口的氧浓度,从而对系统中氧浓度进行控制。
由于铅冷快堆是一个新型的反应堆,各个国家对其研究还尚处于探索阶段。对于氧控技术的研究主要集中在实验方面,而氧控实验需要特制的氧浓度测量仪,造价高且易损坏,因此在实验前需要针对实验工况进行数值模拟预分析。目前尚未报道有关于铅基冷却剂氧浓度扩散数值模拟方法方面的研究。
发明内容
针对上述存在技术问题,本发明提供一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法,可以同时用于气相控氧和固相控氧数值模拟,计算消耗资源低,计算结果可信度高,可以为控氧实验提供预分析,减少实验成本。
为实现上述目的,本发明采取以下技术方案:一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法,步骤如下:
步骤1:氧控装置建模:氧控装置计算模型只包括液体部分,将气体和固体部分简化为给定氧浓度的边界条件;其中气相控氧中认为气液两相交接面维持不变,将气液两相交界面建模为固定壁面;固相控氧中首先采用离散单元法进行氧化铅颗粒随机堆积模拟,利用堆积结果进行流体力学建模,建模中不包括氧化铅颗粒,只针对铅基冷却剂工质流道进行氧控装置建模;
步骤2:网格生成:固相控氧中氧化铅为小球形式,氧化铅小球之间相互接触,因此需要对氧化铅小球直径进行微调便于网格生成,网格生成中需保证氧化铅颗粒表面第一层边界层网格位于层流子层内,便于采用近壁面处理模型进行边界层内氧扩散计算;
步骤3:物性调用和氧控装置计算模型设置:由于氧气在铅基冷却剂中溶解度非常低,因此氧气和铅基冷却剂混合后的物性依然采用铅基冷却剂本身物性;湍流模型选择中,气相控氧模拟选择标准k-ε湍流模型,固相控氧模拟选择SST k-ω湍流模型;边界层处理模型选择中,气相控氧模拟选择Enhanced壁面处理模型,固相控氧模拟使用SST k-ω湍流模型默认的边界层处理模型;对于铅铋工质,湍流施密特数设置为0.9;
步骤4:氧控装置计算模型及边界条件设置:边界条件包括入口流速、入口湍流强度、入口氧浓度、入口温度和边界氧浓度;边界氧浓度根据控氧方式确定,计算公式如下:
若为氢气+氧气型气相控氧或固相控氧,则边界处氧浓度达到饱和,即:
式中:
T——温度,单位为K;
若为氢气+水型气相控氧,则边界处氧浓度为给定氢气和水体积比下的平衡氧气分压对应的平衡氧浓度:
式中:
P0——参考压力,单位为Pa;
K——反应平衡常数;
对应的平衡氧浓度CO,e计算公式:
式中:
CO,e——平衡氧浓度,单位为wt%;
步骤5:采用SIMPLE算法,设置氧气浓度离散格式为二阶迎风格式,开展稳态计算至出口平均氧浓度收敛到一个稳定值。
和现有技术相比较,本发明具备如下优点:
1、本发明的分析方法能够求解铅冷快堆氧控过程中氧浓度扩散分布特性,补全了现有的计算需求空白;
2、本发明针对气相控氧和固相控氧的共性条件进行整合,使得一套计算方法可以同时用于两种控氧方式计算;
3、本发明在充分研究氧控过程的基础上,在维持计算精度的情况下,对氧控物理现象进行了合理简化,从两相计算简化为单相计算,大大节省了计算资源。
附图说明
图1为本发明计算流程框图。
图2为气相控氧装置建模网格划分。
图3为固相控氧装置中氧化铅质量交换器建模。
图4位固相控氧装置中氧化铅小球表面边界层网格设置。
具体实施方式
下面结合附图和具体实施方式对本发明方法作进一步详细说明:
实施案例1:
一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法用于氢气+水气相控氧,步骤如下:
步骤1:氧控装置建模。氧控装置计算模型只包括液体部分,液体部分高度为总高度1/3,将气体部分简化为给定氧浓度的边界条件,认为气液两相交接面维持不变,将气液两相交界面建模为固定壁面;
步骤2:网格生成。如图2所示,图中301为铅铋合金入口,303为铅铋合金出口,302为汽液交界面,在汽液交界面303附近进行了网格的局部加密。
步骤3:物性调用和氧控装置计算模型设置:由于氧气在铅基冷却剂中溶解度非常低,因此氧气和铅基冷却剂混合后的物性依然采用铅基冷却剂本身物性;湍流模型选择中,气相控氧模拟选择标准k-ε湍流模型;边界层处理模型选择中,气相控氧模拟选择Enhanced壁面处理模型;对于铅铋工质,湍流施密特数设置为0.9;
物性参数如下:
ρLBE=11096-1.3236T
Cp,LBE=159-2.72×10-2T+7.12×10-6T2
λLBE=3.61+1.517×10-2T-1.741×10-6T2
μLBE(T)=(4.56-7.03×10-3T+3.61×10-6T2)×10-3
式中:
ρLBE——铅基合金密度/kg·m-3
T——铅基合金温度/K
Cp,LBE——铅基合金热容/J·K-1
λLBE——铅基合金热导率/W·m-1·K-1
μLBE(T)——铅基合金动力粘度/N·s·m-2
DO-LBE(T)——氧气在铅铋合计中扩散系数/m2·s
本实例中,湍流模型选择标准k-ε模型,湍流施密特数为0.9,壁面处理选择Enhanced Wall Treatment。
步骤4:计算模型及边界条件设置。边界条件包括入口流速、入口湍流强度、入口氧浓度、入口温度和边界氧浓度。
本实例中,入口流速0.2m/s,温度643K,湍流强度4.11%,湍流粘度比0.00987,入口氧浓度1×10-8wt%。
本实例采用氢气+水型气相控氧,则边界处氧浓度为给定氢气和水体积比下的平衡氧气分压对应的平衡氧浓度:
式中:
P0——参考压力,单位为Pa;
K——反应平衡常数;
对应的平衡氧浓度CO,e计算公式:
式中:
CO,e——平衡氧浓度,单位为wt%;
步骤5:采用SIMPLE算法,设置氧气浓度离散格式为二阶迎风格式,开展稳态计算至出口平均氧浓度收敛到一个稳定值。
最终计算得到出口出氧浓度为4.01×10-7wt%。
实施案例2:
一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法用于固相控氧,步骤如下:
步骤1:氧控装置建模。氧控装置计算模型只包括液体部分,将气体和固体部分简化为给定氧浓度的边界条件;本实例为固相控氧,因此首先采用离散单元法进行氧化铅颗粒随机堆积模拟,利用堆积结果进行流体力学建模,建模中不包括氧化铅颗粒,只针对铅基冷却剂工质流道进行建模。如图3所示,左侧为铅铋合金入口段,右侧为铅铋合金出口段,黑色阴影球体表示氧化铅颗粒表面。
步骤2:网格生成。固相控氧中氧化铅通常为小球形式,氧化铅小球之间相互接触,因此需要对小球直径进行微调便于网格生成,网格生成中需保证氧化铅颗粒表面第一层边界层网格位于层流子层内,便于采用近壁面处理模型进行边界层内氧扩散计算。图4为氧化铅颗粒网格划分局部放大,图中第一层网格高度为3e-5m。
步骤3:物性调用和计算模型设置。由于氧气在铅基冷却剂中溶解度非常低,因此二者混合后的物性依然采用铅基冷却剂本身物性。湍流模型选择中,固相控氧模拟选择SSTk-ω湍流模型;边界层处理模型选择中,固相控氧模拟使用SST k-ω湍流模型默认的边界层处理模型;对于铅铋工质,湍流施密特数设置为0.9;
步骤4:计算模型及边界条件设置。边界条件包括入口流速、入口湍流强度、入口氧浓度、入口温度和边界氧浓度。
本实例中,入口温度为643K,入口流速0.1m/s,入口氧浓度9.15×10-7wt%。
本实例采用固相控氧,交界面氧浓度为饱和氧浓度,即:
CO,s——边界处氧浓度,单位为wt%;
T——温度,单位为K;
步骤5:采用SIMPLE算法,设置氧气浓度离散格式为二阶迎风格式,开展稳态计算至出口平均氧浓度收敛到一个稳定值。
计算得到出口氧浓度为1.53×10-6wt%。
由以上实例可以看出,本发明可以成功完成气相控氧和固相控氧过程中氧浓度扩散和分布特性,计算结果可以输出质量交换器内氧浓度分布情况和平均出口氧浓度,对质量交换器性能进行评估。
Claims (1)
1.一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法,其特征在于:步骤如下:
步骤1:氧控装置建模:氧控装置计算模型只包括液体部分,将气体和固体部分简化为给定氧浓度的边界条件;其中气相控氧中认为气液两相交接面维持不变,将气液两相交界面建模为固定壁面;固相控氧中首先采用离散单元法进行氧化铅颗粒随机堆积模拟,利用堆积结果进行流体力学建模,建模中不包括氧化铅颗粒,只针对铅基冷却剂工质流道进行氧控装置建模;
步骤2:网格生成:固相控氧中氧化铅为小球形式,氧化铅小球之间相互接触,因此需要对氧化铅小球直径进行微调便于网格生成,网格生成中需保证氧化铅颗粒表面第一层边界层网格位于层流子层内,便于采用近壁面处理模型进行边界层内氧扩散计算;
步骤3:物性调用和氧控装置计算模型设置:由于氧气在铅基冷却剂中溶解度非常低,因此氧气和铅基冷却剂混合后的物性依然采用铅基冷却剂本身物性;湍流模型选择中,气相控氧模拟选择标准k-ε湍流模型,固相控氧模拟选择SST k-ω湍流模型;边界层处理模型选择中,气相控氧模拟选择Enhanced壁面处理模型,固相控氧模拟使用SST k-ω湍流模型默认的边界层处理模型;对于铅铋工质,湍流施密特数设置为0.9;
步骤4:氧控装置计算模型及边界条件设置:边界条件包括入口流速、入口湍流强度、入口氧浓度、入口温度和边界氧浓度;边界氧浓度根据控氧方式确定,计算公式如下:
若为氢气+氧气型气相控氧或固相控氧,则边界处氧浓度达到饱和,即:
式中:
CO,s——边界处氧浓度,单位为wt%;
T——温度,单位为K;
若为氢气+水型气相控氧,则边界处氧浓度为给定氢气和水体积比下的平衡氧气分压对应的平衡氧浓度:
式中:
P0——参考压力,单位为Pa;
K——反应平衡常数;
对应的平衡氧浓度CO,e计算公式:
式中:
CO,e——平衡氧浓度,单位为wt%;
步骤5:采用SIMPLE算法,设置氧气浓度离散格式为二阶迎风格式,开展稳态计算至出口平均氧浓度收敛到一个稳定值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010713015.8A CN111881631B (zh) | 2020-07-22 | 2020-07-22 | 一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010713015.8A CN111881631B (zh) | 2020-07-22 | 2020-07-22 | 一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111881631A true CN111881631A (zh) | 2020-11-03 |
CN111881631B CN111881631B (zh) | 2021-10-19 |
Family
ID=73155341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010713015.8A Active CN111881631B (zh) | 2020-07-22 | 2020-07-22 | 一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111881631B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113486482A (zh) * | 2021-07-12 | 2021-10-08 | 西安交通大学 | 一种液态铅铋外掠螺旋管束速度温度边界层计算方法 |
CN114323991A (zh) * | 2021-11-23 | 2022-04-12 | 华东理工大学 | 一种高温铅铋环境蠕变疲劳试验装置 |
CN114512202A (zh) * | 2021-12-06 | 2022-05-17 | 北京科技大学 | 钛或钛合金粉末表面氧化层形成过程的可视化分析方法 |
CN114722740A (zh) * | 2022-05-05 | 2022-07-08 | 西安交通大学 | 一种液态金属氧化腐蚀数值模拟计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003185788A (ja) * | 2001-12-21 | 2003-07-03 | Mitsui Eng & Shipbuild Co Ltd | 液体金属中の溶解酸素濃度制御方法及び装置 |
JP2003307589A (ja) * | 2002-04-15 | 2003-10-31 | Mitsui Eng & Shipbuild Co Ltd | 液体金属中の溶解酸素濃度制御方法 |
CN105739552A (zh) * | 2016-02-03 | 2016-07-06 | 华北电力大学 | 用于液态铅或液态铅铋合金回路系统的氧控装置及其方法 |
CN106323855A (zh) * | 2016-09-08 | 2017-01-11 | 中国科学院近代物理研究所 | 高温铅铋熔体中自动化控氧/腐蚀装置及方法 |
CN106601314A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院合肥物质科学研究院 | 一种铅基反应堆冷却剂氧控系统和氧控方法 |
EP3163581A1 (en) * | 2014-06-11 | 2017-05-03 | Joint Stock Company "Akme-Engineering" | Method and device for regulating the concentration of oxygen in a reactor facility, and nuclear reactor facility |
-
2020
- 2020-07-22 CN CN202010713015.8A patent/CN111881631B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003185788A (ja) * | 2001-12-21 | 2003-07-03 | Mitsui Eng & Shipbuild Co Ltd | 液体金属中の溶解酸素濃度制御方法及び装置 |
JP2003307589A (ja) * | 2002-04-15 | 2003-10-31 | Mitsui Eng & Shipbuild Co Ltd | 液体金属中の溶解酸素濃度制御方法 |
EP3163581A1 (en) * | 2014-06-11 | 2017-05-03 | Joint Stock Company "Akme-Engineering" | Method and device for regulating the concentration of oxygen in a reactor facility, and nuclear reactor facility |
CN105739552A (zh) * | 2016-02-03 | 2016-07-06 | 华北电力大学 | 用于液态铅或液态铅铋合金回路系统的氧控装置及其方法 |
CN106323855A (zh) * | 2016-09-08 | 2017-01-11 | 中国科学院近代物理研究所 | 高温铅铋熔体中自动化控氧/腐蚀装置及方法 |
CN106601314A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院合肥物质科学研究院 | 一种铅基反应堆冷却剂氧控系统和氧控方法 |
Non-Patent Citations (4)
Title |
---|
A. MARINO等: "Temperature dependence of dissolution rate of a lead oxide mass exchanger in lead–bismuth eutectic", 《JOURNAL OF NUCLEAR MATERIALS》 * |
JINSUO ZHANG: "Oxygen control technology in applications of liquid lead and lead–bismuth systems for mitigating materials corrosion", 《J APPL ELECTROCHEM》 * |
OECD: "Handbook on Lead-bismuth E utectic Alloy and Lead Properties, Materials Compatibility, Thermal hydraulics", 《NUCLEAR ENERGY AGENCY ORGANISATION FOR ECONOMIC CO-OPERATION AND DEVELOPMENT》 * |
魏诗颖等: "铅铋-水直接接触沸水快堆系统子通道耦合热工分析", 《原子能科学技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113486482A (zh) * | 2021-07-12 | 2021-10-08 | 西安交通大学 | 一种液态铅铋外掠螺旋管束速度温度边界层计算方法 |
CN113486482B (zh) * | 2021-07-12 | 2022-10-28 | 西安交通大学 | 一种液态铅铋外掠螺旋管束速度温度边界层计算方法 |
CN114323991A (zh) * | 2021-11-23 | 2022-04-12 | 华东理工大学 | 一种高温铅铋环境蠕变疲劳试验装置 |
CN114512202A (zh) * | 2021-12-06 | 2022-05-17 | 北京科技大学 | 钛或钛合金粉末表面氧化层形成过程的可视化分析方法 |
CN114722740A (zh) * | 2022-05-05 | 2022-07-08 | 西安交通大学 | 一种液态金属氧化腐蚀数值模拟计算方法 |
CN114722740B (zh) * | 2022-05-05 | 2023-07-18 | 西安交通大学 | 一种液态金属氧化腐蚀数值模拟计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111881631B (zh) | 2021-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111881631B (zh) | 一种铅冷快堆控氧过程中氧浓度扩散数值模拟方法 | |
Li et al. | Multi-objective optimization of mini U-channel cold plate with SiO2 nanofluid by RSM and NSGA-II | |
Peng et al. | Improvement of thermal and water management of air-cooled polymer electrolyte membrane fuel cells by adding porous media into the cathode gas channel | |
CN111199122A (zh) | 基于多物理场的质子交换膜燃料电池水淹故障诊断方法 | |
Deymi-Dashtebayaz et al. | Thermo-hydraulic analysis and optimization of CuO/water nanofluid inside helically dimpled heat exchangers | |
US20070065695A1 (en) | Coolant flow estimation for the thermal loop of a fuel cell system using stack loss power | |
Li et al. | Three-dimensional performance simulation of PEMFC of metal foam flow plate reconstructed with improved full morphology | |
Min et al. | Thermal performance analyses of pulsating heat pipe for application in proton exchange member fuel cell | |
CN115358014A (zh) | 一种功率器件风冷散热器模型、优化方法及性能计算方法 | |
Zhang et al. | Experimental and numerical study of an 80-kW zigzag printed circuit heat exchanger for supercritical CO 2 Brayton cycle | |
Abdollahipour et al. | A novel electrochemical refrigeration system based on the combined proton exchange membrane fuel cell-electrolyzer | |
Jiang et al. | Fluid-thermal-mechanical coupled analysis and optimized design of printed circuit heat exchanger with airfoil fins of S-CO2 Brayton cycle | |
Luo et al. | Optimizing a direct flow cooling battery thermal management with bod baffles for electric vehicles: An experimental and simulation study | |
Xie et al. | Influence of perforated gas diffusion layer micropore shape parameters on water removal in a proton exchange membrane fuel cell flow channel | |
Yu et al. | Effects of shark bionic V-groove structures and nanofluids on thermal management of electronic components | |
Ma et al. | Effects of different heat transfer fluids on thermal distribution and electrochemical performance of PEMFC with a non-isothermal multiphase model | |
Liang et al. | Optimization design of proton exchange membrane fuel cell cooling plate based on dual-objective function topology theory | |
Wu et al. | Topological optimization and thermal performance of cold plates for lithium-ion battery with non-uniform heat sources | |
Zhang et al. | Study on heat transfer enhancement performance of cooling channel with elliptical dimples in a proton exchange membrane fuel cell | |
González-Morán et al. | A numerical study on heat transfer for serpentine-type cooling channels in a PEM fuel cell stack | |
CN116881990A (zh) | 一种基于有限元仿真的动力电池液冷板结构优化方法 | |
Yuan et al. | Multi-objective optimization of thermal modules in high heat flux laptops | |
Yang et al. | Effect of cathode operation parameters on the performance of proton exchange membrane fuel cell with curved wing vortex generator by response surface method | |
Sasmito et al. | Numerical investigation of water and temperature distributions for open-cathode polymer electrolyte fuel cell stack with edge cooling | |
CN116631519A (zh) | 一种液态铅铋固态氧控氧传质数值模拟方法 |
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 |