CN105808827B - The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method - Google Patents

The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method Download PDF

Info

Publication number
CN105808827B
CN105808827B CN201610112416.1A CN201610112416A CN105808827B CN 105808827 B CN105808827 B CN 105808827B CN 201610112416 A CN201610112416 A CN 201610112416A CN 105808827 B CN105808827 B CN 105808827B
Authority
CN
China
Prior art keywords
opposing end
end surface
cyclone
dust
efficiency
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
Application number
CN201610112416.1A
Other languages
Chinese (zh)
Other versions
CN105808827A (en
Inventor
张建平
车鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201610112416.1A priority Critical patent/CN105808827B/en
Publication of CN105808827A publication Critical patent/CN105808827A/en
Application granted granted Critical
Publication of CN105808827B publication Critical patent/CN105808827B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Cyclones (AREA)

Abstract

The present invention relates to a kind of determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method, comprising the following steps: 1) construct three-dimensional cyclone dust collectors physical model;2) theoretical model of gas flowfield and particle power field interactions is established;3) different opposing end surfaces are simulated than lower dust removal efficiency in cyclone cleaner, determines that optimal opposing end surface when dust removal efficiency in cyclone cleaner maximum compares region;4) dust removal efficiency in cyclone cleaner is fitted using Fourier's approximating method three times with the change curve of opposing end surface ratio, obtains analytic formula, and verify to it.Compared with prior art, the present invention has many advantages, such as the predictability for greatly improving efficiency of dust collection.

Description

The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method
Technical field
The present invention relates to cyclone dust collectors analogue techniques, more particularly, to a kind of optimal opposing end surface ratio of cyclone dust collectors Determination and three times Fourier's approximating method.
Background technique
Compared with thicker Atmospheric particulates, fine particle partial size is small, often rich in poisonous and harmful substances and in an atmosphere Residence time is long, conveying distance is remote, thus the influence to human health and atmosphere quality is bigger.Cyclone dust collectors for Trapping separation fine particle dust purification efficiency is higher, however fine particle still occupies very big ratio in its discharged particle Weight, therefore the research trapped to fine particle in cyclone dust collectors is particularly significant.Domestic and foreign scholars are directed to cyclone dust collectors at present A series of work has been carried out in fine particle trapping, traps using a variety of different methods, from different sides to fine particle Characteristic explored.
Summary of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide a kind of cyclone dust collectors most The determination of excellent opposing end surface ratio and three times Fourier's approximating method.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method, which is characterized in that packet Include following steps:
1) three-dimensional cyclone dust collectors physical model is constructed;
2) theoretical model of gas flowfield and particle power field interactions is established;
3) different opposing end surfaces are simulated than lower dust removal efficiency in cyclone cleaner, determines dust removal efficiency in cyclone cleaner Optimal opposing end surface when maximum compares region;
4) using Fourier's approximating method three times to dust removal efficiency in cyclone cleaner with opposing end surface ratio change curve into Row fitting, obtains analytic formula, and verify to it.
The theoretical model of the gas flowfield and particle power field interactions specifically:
RSM model is selected to simulate the gas flowfield of cyclone dust collectors, is indicated are as follows:
In formula, ρ is the density of air, kg/m3;U ' bet i, j representation space coordinate;With u 'iIt is particle in the direction i Speed average and pulsating quantity, m/s;CijAnd DijRespectively convective term, DIFFUSION IN TURBULENCE item;PijItem is generated for shear stress;φijFor Pressure-strain item;εijFor viscous dissipation item;
The equation of motion of individual particle is as follows:
In formula, mpAnd upRespectively granular mass and movement velocity;FDjThe drag force of particle is acted on for fluid viscosity;FCj The power of generation is collided between particle, between particle and wall surface.
The solution of the dust removal efficiency in cyclone cleaner is as follows:
The efficiency of dust collection η of deduster refer to dust-contained airflow trapped when passing through cyclone dust collectors the Dust Capacity to get off account for into Enter the percentage of the Dust Capacity of deduster, it is assumed that squeeze into the particle of same particle size equal densities, then efficiency of dust collection η is expressed as trapping The granule number N arrivedtWith incident granule number NiThe ratio between, that is, have:
By change curve of the dust removal efficiency in cyclone cleaner with opposing end surface ratio: gradually with opposing end surface ratio Increase, efficiency of dust collection first increases and then decreases, and opposing end surface ratio efficiency of dust collection in x=5 or so reaches maximum value
Fourier's approximating method is bent with the variation of opposing end surface ratio to dust removal efficiency in cyclone cleaner three times for the utilization Line is fitted specifically:
It is distributed if the efficiency of dust collection η of cyclone dust collectors meets following Fourier three times with opposing end surface ratio x variation:
η (x)=a0+a1cos(ωx)+b1sin(ωx)+a2cos(2ωx)+b2sin(2ωx)+a3cos(3ωx)+ b3sin(3ωx)
In formula, a0、a1、a2、a3、b1、b2、b3, ω be fourier formula three times parameter to be estimated;
It is fitted in conjunction with limited metadata point, and according to Fourier model, it is as follows to obtain fitting formula:
η (x)=0.4061-0.1301cos (x) -0.0383sin (0.5862x)+0.1161cos (1.1724x)
+0.0493sin(1.1724x)+0.0334cos(1.7586x)+0.0452sin(1.7586x)
Compared with prior art, the present invention combines theoretical research and the engineering background of domestic and international cyclone dust collectors, establishes The mathematical model of cyclone dust collectors fine particle efficiency of dust collection, calculates different opposing end surfaces than lower fine particle efficiency of dust collection, And Fourier fitting of distribution, verified approximating method essence are carried out to fine particle efficiency of dust collection change curve based on this Degree is high, can be used as the reference of cyclone dust collectors fine particle efficiency of dust collection calculation method, to greatly improve the pre- of efficiency of dust collection The property surveyed.
Detailed description of the invention
Fig. 1 is cyclone dust collectors physical model schematic diagram;
Fig. 2 is change curve of the efficiency of dust collection with opposing end surface ratio;
Fig. 3 is matched curve and contrast verification curve graph of the efficiency of dust collection with opposing end surface ratio.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment
Using single tube vertical cyclone deduster simplified model as research object, three-dimensional cyclone dust collectors physical model is constructed, is built The theoretical model for having found gas flowfield and particle power field interactions (enters different opposing end surface ratios using SIMPLEC algorithm Mouthful the ratio between sectional area and cylinder sectional area) under dust removal efficiency in cyclone cleaner carried out numerical simulation, it is determined that cyclone dust collectors Optimal opposing end surface when PM2.5 efficiency of dust collection maximum compares region.In addition, using Fourier's approximating method three times to cyclone dust removal Device PM2.5 efficiency of dust collection is fitted with the change curve of opposing end surface ratio, has obtained analytic formula, and tested it Card, the results showed that the approximating method accuracy is high.As a result it can be set for engineering staff about cyclone duster structure and parameter optimization Meter provides Technical Reference.
The solution of 1 cyclone dust collectors fine particle efficiency of dust collection
1.1 gas flowfield
Due to having the stronger strong eddy flow of three-dimensional in cyclone dust collectors, and have apparent anisotropic turbulent flow special Point selects RSM model to simulate the gas flowfield of cyclone dust collectors here, and governing equation is mainly transport equation, writeable are as follows:
In formula, ρ is the density of air, kg/m3;U ' bet i, j, k representation space coordinate;With u 'iIt is particle in the direction i Speed average and pulsating quantity, m/s;CijAnd DijRespectively convective term, DIFFUSION IN TURBULENCE item;PijItem is generated for shear stress;φij For pressure-strain item;εijFor viscous dissipation item.
1.2 particle aerodynamic fields
The individual particle power suffered in deduster flow mainly has: drag force, barometric gradient power, additional mass Power, Saffman lift and gravity etc..Although the power acted on particle is considerably complicated, and not all power is all of equal importance. In Dual-Phrase Distribution of Gas olid, since the density of gas is far smaller than the density of particle, compared with the inertia force of particle itself, buoyancy, The equal very little by comparison in magnitude such as barometric gradient power, additional tension stress and gravity, so that negligible.Therefore, The equation of motion of individual particle can be immediately arrived at from Newton's second law:
In formula, mpAnd upRespectively granular mass and movement velocity;FDjThe drag force of particle is acted on for fluid viscosity;FCj The power of generation is collided between particle, between particle and wall surface.
The solution of 1.3 efficiency of dust collections
The efficiency of dust collection η of deduster refers to that dust-contained airflow traps the Dust Capacity to get off when passing through deduster and accounts for into removing The percentage of the Dust Capacity of dirt device.Assuming that squeezing into the particle of same particle size equal densities, then efficiency of dust collection η, which is represented by, traps Granule number NtWith incident granule number NiThe ratio between, that is, have:
2 apply solution procedure
2.1 physical model
The threedimensional model for the cyclone dust collectors that the present invention selects is as shown in Figure 1, deduster height and diameter are respectively H= 760mm and Do=190mm;Wherein rectangle part is admission line, and admission line height and width are respectively a=95mm and b= 38mm;Big column part is main cylinder, and height h=285mm, the trapezoidal rotary table part in lower part is taper ash bucket, diameter Dc= 72.5mm;The small column on top is outlet pipe, and diameter De=64mm is deep into main cylinder inner part height R=95mm, one Part leakage is in outside, height L=55mm.
2.2 mesh generation
In order to guarantee the precision of numerical result, used maximum mesh unit number is selected when calculating herein here to carry out Simulation calculates.
The setting of 2.3 boundary conditions
Due to inlet velocity it is known that the continuous phase boundray of event import is set as speed import, particle is with air-flow using same Velocity amplitude.In addition, boundary face includes speed entrance, flowing outlet, wall to the numerical simulation that fluid in cyclone dust collectors flows Face, ash discharging hole etc..Specific boundary condition facilities are given in Table 1.
Table 1
3, opposing end surface ratio-efficiency of dust collection curve
In fact, the influence of entrance height and barrel diameter to the inside and outside flow field of cyclone dust collectors is very significant, into And influence fine particle efficiency of dust collection.Opposing end surface ratio is the ratio between cyclone dust collectors entrance section product and cylinder sectional area, is denoted as x.Therefore, efficiency of dust collection is improved and then studying its influence to cyclone dust collectors performance, has become and is badly in need of solving at present The problem of.
Fig. 2 gives dust removal efficiency in cyclone cleaner with the change curve of opposing end surface ratio, it can be found that: with opposite end Face ratio is gradually increased, efficiency of dust collection first increases and then decreases, and opposing end surface ratio efficiency of dust collection in x=5 or so reaches maximum Value.This result provides an important parameter for the size design of cyclone dust collectors, in terms of structure optimization and parameter designing With important directive significance.
4, it is distributed using Fourier three times and carries out opposing end surface ratio-efficiency of dust collection curve matching
Cyclone dust collectors fine particle efficiency of dust collection is convenient for engineering staff with the analytic formula of opposing end surface ratio in order to obtain Design, provides opposing end surface ratio-efficiency of dust collection curve approximating method below with reference to the data point in Fig. 2, utilizes example and right The precision of matched curve is verified.
It is distributed if the efficiency of dust collection η of cyclone dust collectors meets following Fourier three times with opposing end surface ratio x variation:
η (x)=a0+a1cos(ωx)+b1sin(ωx)+a2cos(2ωx)+b2sin(2ωx)+a3cos(3ωx)+ b3sin(3ωx) (4)
In formula, a0、a1、a2、a3、b1、b2、b3, ω be fourier formula three times parameter to be estimated.In conjunction with having in Fig. 2 Metadata point is limited, and is fitted according to Fourier model, it is as follows to obtain fitting formula:
η (x)=0.4061-0.1301cos (x) -0.0383sin (0.5862x)+0.1161cos (1.1724x)
+0.0493sin(1.1724x)+0.0334cos(1.7586x)+0.0452sin(1.7586x) (5)
Fig. 3 gives the data point and its Fourier's matched curve three times of FEM calculation, it can be seen that data point is basic It is distributed in matched curve, coincide very well, this illustrates that fitting formula most preferably reflects efficiency of dust collection with opposing end surface ratio Variation relation.
Example:
As opposing end surface ratio x=4.6, the FEM calculation value of efficiency of dust collection can be obtained by the calculating step of this paper, The value corresponds to the solid dot in Fig. 3, and efficiency of dust collection is η=0.577.X=4.6 is substituted into formula (5), efficiency of dust collection is obtained Match value is η*=0.575.On the basis of FEM calculation value, the efficiency of dust collection that Fourier's fitting of distribution three times can be obtained is opposite Error delta is
It can be seen that efficiency of dust collection calculation method of the invention and three times the determination method of Fourier's fitting of distribution formula It is practical, for efficiency of dust collection η with opposing end surface ratio x change curve, either monotone increasing or monotonic decreasing, Fu three times In leaf fitting of distribution all have very high fitting precision.

Claims (5)

1. a kind of determination method of the optimal opposing end surface ratio of cyclone dust collectors, which comprises the following steps:
1) three-dimensional cyclone dust collectors physical model is constructed;
2) theoretical model of gas flowfield and particle power field interactions is established;
3) different opposing end surfaces are simulated than lower dust removal efficiency in cyclone cleaner, determines dust removal efficiency in cyclone cleaner maximum When optimal opposing end surface compare region;
4) dust removal efficiency in cyclone cleaner is intended with the change curve of opposing end surface ratio using Fourier's approximating method three times It closes, obtains analytic formula, and verify to it.
2. the method according to claim 1, wherein the gas flowfield and particle power field interactions Theoretical model specifically:
RSM model is selected to simulate the gas flowfield of cyclone dust collectors, is indicated are as follows:
In formula, ρ is the density of air, kg/m3;U ' bet i, j representation space coordinate;With u 'iSpeed for particle in the direction i Average value and pulsating quantity, m/s;CijAnd DijRespectively convective term, DIFFUSION IN TURBULENCE item;PijItem is generated for shear stress;φijFor pressure Dependent variable;εijFor viscous dissipation item;
The equation of motion of individual particle is as follows:
In formula, mpAnd upRespectively granular mass and movement velocity;FDjThe drag force of particle is acted on for fluid viscosity;FCjFor The power of generation is collided between grain, between particle and wall surface.
3. the method according to claim 1, wherein the solution of the dust removal efficiency in cyclone cleaner is as follows:
The efficiency of dust collection η of deduster refers to that dust-contained airflow traps the Dust Capacity to get off when passing through cyclone dust collectors and accounts for into removing The percentage of the Dust Capacity of dirt device, it is assumed that squeeze into the particle of same particle size equal densities, then efficiency of dust collection η is expressed as trapping Granule number NtWith incident granule number NiThe ratio between, that is, have:
4. the method according to claim 1, wherein by dust removal efficiency in cyclone cleaner with opposing end surface ratio Known to change curve: with being gradually increased for opposing end surface ratio, efficiency of dust collection first increases and then decreases, and opposing end surface ratio x=5 is left Efficiency of dust collection reaches maximum value when right.
5. the method according to claim 1, wherein Fourier's approximating method removes whirlwind three times for the utilization Dirt device efficiency of dust collection is fitted with the change curve of opposing end surface ratio specifically:
It is distributed if the efficiency of dust collection η of cyclone dust collectors meets following Fourier three times with opposing end surface ratio x variation:
η (x)=a0+a1cos(ωx)+b1sin(ωx)+a2cos(2ωx)+b2sin(2ωx)+a3cos(3ωx)+b3sin(3ω x) (4)
In formula, a0、a1、a2、a3、b1、b2、b3, ω be fourier formula three times parameter to be estimated;
It is fitted in conjunction with limited metadata point, and according to Fourier model, it is as follows to obtain fitting formula:
η (x)=0.4061-0.1301cos (x) -0.0383sin (0.5862x)+0.1161cos (1.1724x)
+0.0493sin(1.1724x)+0.0334cos(1.7586x)+0.0452sin(1.7586x) (5)。
CN201610112416.1A 2016-02-29 2016-02-29 The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method Active CN105808827B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610112416.1A CN105808827B (en) 2016-02-29 2016-02-29 The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610112416.1A CN105808827B (en) 2016-02-29 2016-02-29 The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method

Publications (2)

Publication Number Publication Date
CN105808827A CN105808827A (en) 2016-07-27
CN105808827B true CN105808827B (en) 2019-04-09

Family

ID=56465886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610112416.1A Active CN105808827B (en) 2016-02-29 2016-02-29 The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method

Country Status (1)

Country Link
CN (1) CN105808827B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106980710A (en) * 2017-02-28 2017-07-25 上海电力学院 ESP operating voltages contribution rate and magnetic induction intensity curve-fitting method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5961701A (en) * 1997-03-24 1999-10-05 Vision Almet Limited Moisture separator for digester gases and landfill gases and raw natural gases
CN101720255A (en) * 2007-04-30 2010-06-02 罗穆阿尔多·路易斯·里韦拉·萨尔塞多 The process of electrostatic recirculation and the device that are used for dedusting and gas purification

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5961701A (en) * 1997-03-24 1999-10-05 Vision Almet Limited Moisture separator for digester gases and landfill gases and raw natural gases
CN101720255A (en) * 2007-04-30 2010-06-02 罗穆阿尔多·路易斯·里韦拉·萨尔塞多 The process of electrostatic recirculation and the device that are used for dedusting and gas purification

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Study on Flow Field Intensification of Dynamic Pre-cyclone Separators;ZOU Jiupeng 等;《Advanced Materials Research》;20120726;第550-553卷;第3210-3213页
旋风除尘器结构尺寸优化设计的研究;赵萍;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20030315;第2013年卷(第01期);第B027-48页

Also Published As

Publication number Publication date
CN105808827A (en) 2016-07-27

Similar Documents

Publication Publication Date Title
Kumar et al. Multi-objective shape optimization of vortex finders in cyclone separators using response surface methodology and genetic algorithms
CN108319793B (en) CFD-based analysis method for free movement of low-speed inclined water inlet
Noh et al. Performance improvement of a cyclone separator using multiple subsidiary cyclones
CN107832494A (en) Hypersonic aircraft leading edge stream thermosetting integration computational methods
CN104200102B (en) A kind of photovoltaic array wind load Forecasting Methodology based on CFD
CN112883667A (en) Membrane structure wind resistance optimization method based on large vortex simulation under vertical pulsating wind load effect
Zhu et al. Improvement of the power extraction performance of a semi-active flapping airfoil by employing two-sided symmetric slot airfoil
Sun et al. Design of vortex finder structure for decreasing the pressure drop of a cyclone separator
CN105808827B (en) The determination of the optimal opposing end surface ratio of cyclone dust collectors and three times Fourier's approximating method
CN102284360B (en) Method for determining trapping efficiency of electrostatic dust collector on inhalable particulate matters (PM10)
Kubilay Numerical simulations and field experiments of wetting of building facades due to wind-driven rain in urban areas
Mishra et al. Prediction of separation induced transition on thick airfoil using non-linear URANS based turbulence model
Westphalen et al. Extreme wave loading on offshore wave energy devices using CFD: a hierarchical team approach
Souza et al. Preliminary results of large eddy simulations of a hydrocyclone
Matvienko Analysis of turbulence models and investigation of the structure of the flow in a hydrocyclone
Brar et al. Performance analysis of multi-inlet cyclone separators considering different shapes and locations of the inlet ducts
Schmidt et al. Simulation of turbulent flow in a cyclonic separator
Liu et al. Review of application of VOF-based NWT on integrated OWC system
Cao et al. Orthogonal optimization for effective classification of different tea leaves by a novel pressure stabilized inclined chamber classifier
Plashikhin Computer simulation of the hydrodynamic processes of cyclone dust collectors
CN104036126A (en) Numerical simulation method and system for particulate matter motion in contracting and expanding process of pulmonary alveoli
Dai et al. Flow and deposition characteristics of nanoparticles in a 90 square bend
Mahmoudi et al. Numerical study of low Reynolds inertial impactors with elliptical impaction plate
Xu et al. Using CFD based virtual sensor data to study the structure of air flow behind a porous fence
Çalişkan et al. Design and performance analysis of a virtual body mini cyclone

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