CN110705185A - 预测管道气锤的方法 - Google Patents
预测管道气锤的方法 Download PDFInfo
- Publication number
- CN110705185A CN110705185A CN201910930682.9A CN201910930682A CN110705185A CN 110705185 A CN110705185 A CN 110705185A CN 201910930682 A CN201910930682 A CN 201910930682A CN 110705185 A CN110705185 A CN 110705185A
- Authority
- CN
- China
- Prior art keywords
- particle
- equation
- information
- gas
- density
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000002245 particle Substances 0.000 claims abstract description 111
- 239000007789 gas Substances 0.000 claims abstract description 48
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 46
- 239000003345 natural gas Substances 0.000 claims abstract description 23
- 239000012530 fluid Substances 0.000 claims abstract description 18
- 230000035939 shock Effects 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000010354 integration Effects 0.000 claims description 13
- 230000001351 cycling effect Effects 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 5
- 238000011144 upstream manufacturing Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及流体力学技术领域,为提出提供一种求解天然气管道传输问题的方法。为此,本发明采取的技术方案是,预测管道气锤的方法,考虑气体的强可压缩性带来的影响,采用光滑粒子流体动力学方法求解移动坐标系统下的激波管物理方程,模拟气体对管道的冲击。本发明主要应用于天然气管道等设计制造场合。
Description
技术领域
本发明涉及流体力学技术领域,具体是涉及一种基于天然气管道传输问题的无网格粒子方法。
背景技术
针对天然气管道气体的输运中,由于气体的可压缩性强,所以针对气体的计算和模拟和其他流体有明显的不同,对于这一问题,常用的方法主要有基于特征线法的界面追踪技术和有限体积法等。为了降低气态垂直界面假设引起的计算误差,通过假定气体状态在单个网格距离内是线性变化的,提出了改进的界面追踪技术,但这影响了计算的效率和稳定性,也可能因插值误差的积累而降低计算精度,针对以上问题,本发明提出了一种基于拉格朗粒子模型的无网格粒子方法。
发明内容
为克服现有技术的不足,本发明旨在提出提供一种求解天然气管道传输问题的方法。为此,本发明采取的技术方案是,预测管道气锤的方法,考虑气体的强可压缩性带来的影响,采用光滑粒子流体动力学方法求解移动坐标系统下的激波管物理方程,模拟气体对管道的冲击。
具体步骤如下:
步骤一,初始化:初始化系统的相关变量和粒子信息;
步骤二,生成粒子信息;
步骤三,列出求解方程并迭代计算:
根据天然气管道问题的原理,得出物理模型的数学模型,即拉格朗日形式下的控制方程为:
其中P是气体压力,ν是气体的速度,ρ是气体的密度,e是气体的能量,对于理想气体,其状态方程为:
p=(γ-1)ρe (4)
将问题简化,相应的状态方程为:
此时,能量项将从控制方程中解耦出来;
为了求解上述方程(5),从而预测天然气管道中气锤问题,采取步骤如下:
光滑粒子流体动力学方法中,函数f(x)的积分表示式为:
f(x)=∫Ωf(x′)δ(x-x′)dx′ (6)
其中δ(x-x′)为狄拉克函数,Ω为包含x的积分体积,若用光滑函数W(x-x′,h)取代狄拉克函数,则f(x)的积分表示式为:
f(x)≈∫Ωf(x′)W(x-x′,h)dx′ (7)
函数导数积分表示为:
又因为
所以
然后进行粒子近似得到:
因为
所以
又有
所以,在粒子i处的函数的粒子近似式写为:
利用光滑粒子流体动力学方法将方程(1)~(3)离散可得(16)~(18)式
p=(γ-1)ρe (19)
其中∏ij为Monaghon型人工粘度,Hij为人工热量项,γ、ξ为常系数;
根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强;
步骤四,输出结果:
1)每一个时间步结束,保存中间结果,并输出中间结果;
2)结束时间循环,并输出最终结果。
根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强,具体计算过程为:
1)循环每个时间步;
2)初始化计算域内所有粒子属性信息后,搜索目标粒子的邻近粒子,获得初始值的速度和质量,然后通过方程(16)计算该粒子的密度的导数,最后通过欧拉时间积分更新粒子的密度信息;
3)通过上一步得到的更新的密度信息、初始设置的压力值和Monaghon型人工粘度方程计算得到速度导数,然后通过欧拉时间积分更新粒子的速度信息,最后通过速度信息更新粒子的位置;
4)根据上面两步得到的密度、速度粒子信息和Monaghon型人工粘度方程及人工热量项方程计算得到能量导数,然后通过欧拉时间积分更新粒子的能量信息;
5)通过上面得到的密度和能量信息结合方程(19),更新粒子的压力;
6)然后通过(20)式更新粒子搜索的光滑长度信息。
本发明的特点及有益效果是:
本发明采用光滑粒子流体动力学方法求解移动坐标系统下的激波管物理方程,充分考虑了气体的强可压缩性带来的影响,在满足数值精度的前提下能够更方便地模拟气体对管道的冲击问题。
附图说明:
图1程序流程图。
图2天然气管道左边冲击右边扩张问题物理模型。
图3 CASE1中天然气管道在0.6s时刻密度、压强、速度和能量随位移的分布图。
图4一维天然气管道气锤问题物理模型。
图5 CASE2中天然气管道中间点压力和密度随时间分布。
具体实施方式
本发明解决的技术问题是提供一种求解天然气管道传输问题的拉格朗日粒子方法,方法采用光滑粒子流体动力学方法求解移动坐标系统下的激波管物理方程,充分考虑了气体的强可压缩性带来的影响,在满足数值精度的前提下能够更方便地模拟气体对管道的冲击问题。
步骤一,初始化:初始化系统的相关变量和粒子信息;
步骤二,生成粒子信息;
步骤三,列出求解方程并迭代计算:
根据天然气管道问题的原理,可以得出物理模型的数学模型,即拉格朗日形式下的控制方程为:
其中P是气体压力,ν是气体的速度,ρ是气体的密度,e是气体的能量。(1)~(3)中的拉格朗日控制方程是未闭合的,需要添加状态方程来保证系统完
整性。对于理想气体,其状态方程为:
p=(γ-1)ρe (4)
对于实际管道中的气锤问题,如果温度变化不大,可以将问题简化,相应的状态方程为:
此时,能量项将从控制方程中解耦出来。
为了求解上述方程,从而预测天然气管道中气锤问题,本发明提出了一种拉格朗日粒子方法。具体内容如下:
光滑粒子流体动力学方法中,函数f(x)的积分表示式为:
f(x)=∫Ωf(x′)δ(x-x′)dx′ (6)
其中δ(x-x′)为狄拉克函数,Ω为包含x的积分体积。若用光滑函数W(x-x′,h)取代狄拉克函数,则f(x)的积分表示式为:
f(x)≈∫Ωf(x′)W(x-x′,h)dx′ (7)
函数导数积分表示为:
又因为
所以
然后进行粒子近似得到:
因为
所以
又有
所以,在粒子i处的函数的粒子近似式可写为:
利用光滑粒子流体动力学方法将方程(1)~(3)离散可得(16)~(18)式
p=(γ-1)ρe (19)
其中Πij为Monaghon型人工粘度,Hij为人工热量项,γ、ξ为常系数。
根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强。具体计算过程如下所示:
具体计算过程为:
6)循环每个时间步;
7)初始化计算域内所有粒子属性信息后,搜索目标粒子的邻近粒子,获得初始值的速度和质量,然后通过方程(16)计算该粒子的密度的导数,最后通过欧拉时间积分更新粒子的密度信息;
8)通过上一步得到的更新的密度信息、初始设置的压力值和Monaghon型人工粘度方程计算得到速度导数,然后通过欧拉时间积分更新粒子的速度信息,最后通过速度信息更新粒子的位置;
9)根据上面两步得到的密度、速度粒子信息和Monaghon型(由Monaghon在1992年提出解决数值震荡的方法)人工粘度方程及人工热量项方程计算得到能量导数,然后通过欧拉时间积分更新粒子的能量信息;
10)通过上面得到的密度和能量信息结合方程(19),更新粒子的压力;
6)然后通过(20)式更新粒子搜索的光滑长度信息。
步骤四,输出结果:
1)每一个时间步结束,保存中间结果,并输出中间结果;
2)结束时间循环,并输出最终结果。
进一步地,在上述方案中,所述初始化变量信息和运行参数具体设置如下:
CASE 1:天然气管道左边冲击右边扩张问题
本实验模拟问题的物理模型见图2,计算域是长度为30m的一维空间,管道的两侧均是打开状态,左端设置为入口,右端设置的是出口,管道的中间位置设置了一层隔板来模拟天然气管道中的阀门。在初始时刻,隔板两侧存在两种不同状态的气体,通过打开中间隔板来模拟实验问题。本实验中我们在出入口都添加了虚粒子,用来弥补边界粒子处的缺失问题。
在CASE1问题中,我们初始条件设置为:PL=7,VL=0,ρL=1,PR=10,ρR=1,
VR=0。在本实验中,我们在管道左端布置1000个质量相同粒子,管道的右端布置1000个质量相同的粒子。实验模拟的时长为1s,时间积分步长为10e-5s。Monaghon型人工粘度的参数α=1,β=1;人工热量项的系数为g1=0.2,g2=0.4。
基于管道方程的SPH离散格式(15)~(17),模拟一维天然气管道的拉格朗日粒子方法为:
步骤一,初始化。初始化系统的相关变量和运行参数,具体包括:
如图1所示,天然气管道长度L为30m,管道左端0~15m处的压力、速度、密度分别为:7Pa、0m/s、1kg/m3,右端15~30m处依次分别为:10Pa、0m/s、1kg/m3。粒子间距Δx0=0.01m,计算时间步长为10e-5s,计算时间为1s。在本实验中,采用可变光滑长度进行计算,初始光滑长度h0=1.5Δx。取三次样条函数作为核函数。
步骤二,生成粒子信息,具体包括:
初始化粒子步骤中,总共生成气体粒子2000个(不包括虚粒子),在计算域的左右两侧初始边界布置虚粒子,保证距左右两测初始处2h范围内的气体粒子能够进行正常运算,虚粒子的密度、速度、压力和能量与气体粒子相同。
步骤三,列出求解方程并迭代计算:
其中P是气体压力,ν是气体的速度,ρ是气体的密度,e是气体的能量。(1)~(3)中的拉格朗日控制方程是未闭合的,需要添加状态方程来保证系统的完整性。对于理想气体,其状态方程为:
p=(γ-1)ρe (4)
对于实际管道中的气锤问题,如果温度变化不大,可以将问题简化,相应的状态方程为:
此时,能量项将从控制方程中解耦出来。
光滑粒子流体动力学方法中,函数f(x)的积分表示式为:
f(x)=∫Ωf(x′)δ(x-x′)dx′ (6)
其中δ(x-x′)为狄拉克函数,Ω为包含x的积分体积。若用光滑函数W(x-x′,h)取代狄拉克函数,则f(x)的积分表示式为:
f(x)≈∫Ωf(x′)W(x-x′,h)dx′ (7)
函数导数积分表示为:
又因为
所以
然后进行粒子近似得到:
因为
所以
又有
所以,在粒子i处的函数的粒子近似式可写为:
利用光滑粒子流体动力学方法将方程(1)~(3)离散可得(16)~(18)式
p=(γ-1)ρe (19)
其中Πij为Monaghon型人工粘度,Hij为人工热量项,γ、ξ为常系数。
根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强。具体计算过程如下所示:
具体计算过程为:
1)循环每个时间步;
2)初始化计算域内所有粒子属性信息后,搜索目标粒子的邻近粒子,获得初始值的速度和质量,然后通过方程(16)计算该粒子的密度的导数,最后通过欧拉时间积分更新粒子的密度信息;
3)通过上一步得到的更新的密度信息、初始设置的压力值和Monaghon型人工粘度方程计算得到速度导数,然后通过欧拉时间积分更新粒子的速度信息,最后通过速度信息更新粒子的位置;
4)根据上面两步得到的密度、速度粒子信息和Monaghon型人工粘度方程及人工热量项方程计算得到能量导数,然后通过欧拉时间积分更新粒子的能量信息;
5)通过上面得到的密度和能量信息结合方程(19),更新粒子的压力;
6)然后通过(20)式更新粒子搜索的光滑长度信息。
步骤四,输出结果:
1)每一个时间步结束,保存中间结果,并输出中间结果;
2)结束时间循环,并输出最终结果。
CASE 2:管道气锤问题
本实验模拟问题的物理模型见图4,天然气管道长度L为37.5m,气罐内压力PR为250Kpa,管道内气体的初始速度为10m/s,管道直径D为0.0221m,气体的密度ρ为1.2kg/m3,摩擦系数f为0.02。天然气传输的初始条件为V(x,0)=V0和0<x≤L。在数值模拟过程中,初始均匀分布了376个粒子(包括上下游的虚粒子),计算时间步长为0.0001s,计算总时间为46.524s。
在一个实例中,具体步骤如下:
步骤一,初始化:初始化相关的变量和粒子(包括添加的虚粒子)信息。具体包括:
1)初始化与问题相关的变量信息:初始化变量信息管道直径D为0.0221m,天然气管道长度L为37.5m,气罐内压力PR为250Kpa,气体的密度ρ为1.2kg/m3,管道内气体的初始速度为10m/s,计算时间步长为0.0001s,计算总时长为46.524s等;
3)初始化虚粒子信息,在流体上下游边界分别布上两层虚粒子并根据边界条件添加初始信息:上下游边界各两个虚粒子,上游虚粒子压力为PR,下游虚粒子压力为0,上下游虚粒子初始速度均为V0。
步骤二,列出求解方程并迭代计算,步骤同CASE1中的步骤三,将状态方程替换成方程(5)。
步骤三同CASE1中步骤四。
尽管上面结合图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以做出很多变形,这些均属于本发明的保护之内。
Claims (3)
1.一种预测管道气锤的方法,其特征是,考虑气体的强可压缩性带来的影响,采用光滑粒子流体动力学方法求解移动坐标系统下的激波管物理方程,模拟气体对管道的冲击。
2.如权利要求1所述的预测管道气锤的方法,其特征是,具体步骤如下:
初始化:初始化系统的相关变量和粒子信息;
生成粒子信息;
列出求解方程并迭代计算:
根据天然气管道问题的原理,得出物理模型的数学模型,即拉格朗日形式下的控制方程为:
其中P是气体压力,ν是气体的速度,ρ是气体的密度,e是气体的能量,对于理想气体,其状态方程为:
p=(γ-1)ρe (4)
将问题简化,相应的状态方程为:
此时,能量项将从控制方程中解耦出来;
为了求解上述方程(5),从而预测天然气管道中气锤问题,采取步骤如下:
光滑粒子流体动力学方法中,函数f(x)的积分表示式为:
f(x)=∫Ωf(x′)δ(x-x′)dx′ (6)
其中δ(x-x′)为狄拉克函数,Ω为包含x的积分体积,若用光滑函数W(x-x′,h)取代狄拉克函数,则f(x)的积分表示式为:
f(x)≈∫Ωf(x′)W(x-x′,h)dx′ (7)
函数导数积分表示为:
又因为
所以
然后进行粒子近似得到:
因为
所以
又有
所以,在粒子i处的函数的粒子近似式写为:
利用光滑粒子流体动力学方法将方程(1)~(3)离散可得(16)~(18)式
p=(γ-1)ρe (19)
其中Πij为Monaghon型人工粘度,Hij为人工热量项,γ、ξ为常系数;
根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强;
输出结果:
1)每一个时间步结束,保存中间结果,并输出中间结果;
2)结束时间循环,并输出最终结果。
3.如权利要求2所述的预测管道气锤的方法,其特征是,根据方程(16)、(17)和(18)计算不同时刻每个粒子的密度、速度和能量信息,进而根据(19)式求解出粒子的压强,具体计算过程为:
1)循环每个时间步;
2)初始化计算域内所有粒子属性信息后,搜索目标粒子的邻近粒子,获得初始值的速度和质量,然后通过方程(16)计算该粒子的密度的导数,最后通过欧拉时间积分更新粒子的密度信息;
3)通过上一步得到的更新的密度信息、初始设置的压力值和Monaghon型人工粘度方程计算得到速度导数,然后通过欧拉时间积分更新粒子的速度信息,最后通过速度信息更新粒子的位置;
4)根据上面两步得到的密度、速度粒子信息和Monaghon型人工粘度方程及人工热量项方程计算得到能量导数,然后通过欧拉时间积分更新粒子的能量信息;
5)通过上面得到的密度和能量信息结合方程(19),更新粒子的压力;
6)然后通过(20)式更新粒子搜索的光滑长度信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910930682.9A CN110705185A (zh) | 2019-09-29 | 2019-09-29 | 预测管道气锤的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910930682.9A CN110705185A (zh) | 2019-09-29 | 2019-09-29 | 预测管道气锤的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110705185A true CN110705185A (zh) | 2020-01-17 |
Family
ID=69196293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910930682.9A Pending CN110705185A (zh) | 2019-09-29 | 2019-09-29 | 预测管道气锤的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110705185A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111339658A (zh) * | 2020-02-25 | 2020-06-26 | 河海大学 | 基于拉格朗日无网格粒子法的水力瞬变模拟方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473418A (zh) * | 2013-09-18 | 2013-12-25 | 梧州学院 | 弯曲边界上的低雷诺数不可压缩流中的压力差的仿真方法 |
CN106570308A (zh) * | 2016-09-21 | 2017-04-19 | 天津大学 | 一种分析含截留气团管道瞬变流的无网格粒子方法 |
CN109948109A (zh) * | 2019-01-31 | 2019-06-28 | 天津大学 | 含截面变化的明渠非恒定流无网格粒子模拟方法 |
CN110287590A (zh) * | 2019-06-24 | 2019-09-27 | 天津大学 | 基于算子分裂及改进半拉格朗日求解污染物传播的方法 |
-
2019
- 2019-09-29 CN CN201910930682.9A patent/CN110705185A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473418A (zh) * | 2013-09-18 | 2013-12-25 | 梧州学院 | 弯曲边界上的低雷诺数不可压缩流中的压力差的仿真方法 |
CN106570308A (zh) * | 2016-09-21 | 2017-04-19 | 天津大学 | 一种分析含截留气团管道瞬变流的无网格粒子方法 |
CN109948109A (zh) * | 2019-01-31 | 2019-06-28 | 天津大学 | 含截面变化的明渠非恒定流无网格粒子模拟方法 |
CN110287590A (zh) * | 2019-06-24 | 2019-09-27 | 天津大学 | 基于算子分裂及改进半拉格朗日求解污染物传播的方法 |
Non-Patent Citations (3)
Title |
---|
PURI K等: "A comparison of SPH schemes for the compressible Euler equations" * |
林昊: "基于光滑粒子流体动力学法的OPenGL可视化" * |
黄春营: "基于拉格朗日粒子模型的气液混流和冲击模拟" * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111339658A (zh) * | 2020-02-25 | 2020-06-26 | 河海大学 | 基于拉格朗日无网格粒子法的水力瞬变模拟方法及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304684B (zh) | 一种火箭发动机尾喷射流仿真方法及系统 | |
Behbahani-Nejad et al. | The accuracy and efficiency of a reduced-order model for transient flow analysis in gas pipelines | |
MX2013008812A (es) | Sistema y metodo para utilizar una red artificial neural para simular tuberias hidraulicas en un simulador de deposito. | |
CN105468844A (zh) | 管道内水-气耦合瞬变流的模拟方法 | |
CN114168796B (zh) | 一种建立飞行器高空气动力数据库的方法 | |
CN113609597B (zh) | 超声速飞行器绕流的时间-空间混合推进扰动域更新方法 | |
CN117408189B (zh) | 高超声速边界层的转捩预测方法、装置、设备及存储介质 | |
CN111209703A (zh) | 一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统 | |
CN110705185A (zh) | 预测管道气锤的方法 | |
Rodriguez et al. | Formulation and implementation of inflow/outflow boundary conditions to simulate propulsive effects | |
CN106777770B (zh) | 基于有限体积法的输水管道中空穴流的模拟方法 | |
WO2021133627A1 (en) | Intelligent completion control in reservoir modeling | |
Zhang et al. | Coupling of level set and volume of fluid methods for simulations of transient internal flow field in solid rocket motors | |
CN115310339A (zh) | 基于物质点法的具有表面张力效应的固液耦合模拟方法 | |
CN114239199A (zh) | 一种考虑凝结水的蒸汽管网动态仿真方法 | |
CN113505544A (zh) | 基于有限体积法自行车运动虚拟数值风洞系统 | |
CN109002419B (zh) | 天然气管网的动态分析方法及装置 | |
Sachdev | Parallel solution-adaptive method for predicting solid propellant rocket motor core flows | |
Vasconcellos et al. | Flow past a circular cylinder: a comparison between commercial finite volume and finite element codes | |
Yu et al. | Large eddy simulation of the cavity shedding characteristics of ventilated cavitation around the underwater vehicle | |
CN110569541A (zh) | 管道水锤分析方法 | |
Boylston | Quasi-One-Dimensional Flow for Use in Real-Time Facility Simulations | |
Yuan et al. | An efficient topology optimization method for steady gas flows in all flow regimes | |
Hendrix et al. | Simulation of slug propagation for by-pass pigging in two-phase stratified pipe flow | |
Alldieck et al. | The effect of rough-wall boundary conditions on rans-based transition prediction |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200117 |