CN105067621A - Method for judging multi-phase mixing uniformity - Google Patents

Method for judging multi-phase mixing uniformity Download PDF

Info

Publication number
CN105067621A
CN105067621A CN201510503871.XA CN201510503871A CN105067621A CN 105067621 A CN105067621 A CN 105067621A CN 201510503871 A CN201510503871 A CN 201510503871A CN 105067621 A CN105067621 A CN 105067621A
Authority
CN
China
Prior art keywords
pattern
mixing uniformity
moment
calculated
target
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
Application number
CN201510503871.XA
Other languages
Chinese (zh)
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.)
Kunming University of Science and Technology
Yunnan University of Finance and Economics
Original Assignee
Kunming University of Science and Technology
Yunnan University of Finance and Economics
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 Kunming University of Science and Technology, Yunnan University of Finance and Economics filed Critical Kunming University of Science and Technology
Priority to CN201510503871.XA priority Critical patent/CN105067621A/en
Publication of CN105067621A publication Critical patent/CN105067621A/en
Pending legal-status Critical Current

Links

Abstract

The invention discloses a method for judging the multi-phase mixing uniformity, and belongs to the technical field of chemical engineering. The method for judging the multi-phase mixing uniformity based on the uniform experiment design is mainly used for judging a fluid mixing effect and theoretically guiding the design of process parameters in chemical and metallurgical stirring equipment. The method for judging the multi-phase mixing uniformity is simple and convenient and has a very high practical value; the simple and practical method is provided for judging the mixing uniformity and theoretically guiding the design of a stirring reactor in the chemical and metallurgical stirring equipment.

Description

A kind of method judging multi-phase mixing uniformity
Technical field
The present invention relates to a kind of method judging multi-phase mixing uniformity, belong to technical field of chemical engineering.
Background technology
Stirring is one of modal operation in chemical industry, metallurgical process.During 60 ~ eighties, fluid mixing techniques obtains development at full speed greatly, and emphasis of its research mainly carries out experimental study for convention stir oar in the low sticky macroscopic quantity such as agitated flow, incorporation time of gluing in the heterogeneous systems such as non newtonian homogeneous system, solid-liquid suspension and Gas-Liquid Dispersion with height.Although what have a large amount of design experiences and correlation can be used for analysis and prediction mixed system, but stirred reactor is directly amplified to large-scale commercial production from laboratory scale, being still and being not sure, still needing so far to reach mass transfer required by mixing plant, heat transfer and mixing by amplifying step by step.This method not only expends financial resources and a large amount of man power and materials, and the design cycle is very long, and the chemical industry according to the ASSOCIATE STATISTICS display U.S. is about l0 ~ 10,000,000,000 dollar every year because stirred reactor designs unreasonable caused loss.Therefore advanced means of testing and the velocity field set up in rational mathematical model acquisition tank diameter, temperature field and concentration field is adopted, and take the method for effective measuring mixing effect of fluid, not only to the optimal design of mixing apparatus, there is very important economic implications, and to the fundamental research of amplifying and mix, there is the theory significance of reality.Just see at present, the method of measuring mixing effect of fluid is numerous, mainly contains electrical conductivity method, thermocouple method, optical method, decoloring method etc., and wherein electrical conductivity method is applied very wide in being uniformly mixed of low viscosity fluid, but high to the requirement of agitated medium, as required, agitated medium is deionized water; Thermopair rule stream field can produce destruction; General optical method, due to device complexity, is applied also less; It is very effective that decoloring method is used for surveying highly viscous fluid and viscoelastic fluid, but owing to adopting naked eyes to judge, thus with stronger subjectivity.
Summary of the invention
Instant invention overcomes the deficiency of said method, a kind of that have higher using value, simple and feasible method for judging multi-phase mixing uniformity being provided, specifically comprising the following steps:
(1) in the mixed process of transparent tank diameter, obtain continuous print in gas-liquid-solid heterogeneous mixed process and stir real-time pattern and preserve; And utilize digital image processing techniques to obtain the pending pattern of only having powerful connections with target; Using the target in pending pattern as research object;
(2) according to the size of pending original pattern, corresponding from the upper left corner , the lower left corner is corresponding , the lower right corner is corresponding , the upper right corner is corresponding carry out stress and strain model as reference position, divide respectively individual region for natural number, and each zonule is numbered ;
(3) use uniform experiment design theoretical, tolerance tthe method of moment mixing uniformity calculates as follows:
1. the number ratio with the target in overall pattern in a zonule is calculated: ;
2. the area ratio in this zonule and overall pattern is calculated: ;
3. target number is calculated than the absolute deviation with region pattern area ratio: ;
4. the mean value of absolute deviation is calculated: ;
5. tthe uniformity coefficient of moment pattern is: , wherein represent median;
(4) theoretical by uniform experiment design, in incorporation time tfront and back, obvious Changing Pattern can be presented; If experience sufficient incorporation time tafter, numerical values recited fluctuation be bordering on 1, if draw change curve, then exist tafter moment, just be in close proximity to 1 and in time backward fluctuation little, then represent heterogeneous being blended in tthe state mixed is reached after moment.
If the division obtaining real-time pattern grid in mixed process is single from some reference positions, be generally divided into individual region ( for natural number), the division principle in region judges accuracy requirement to determine according to pattern and target sizes, mixing uniformity value size, the larger precision of value higher, the length breadth ratio of zonule is equal with the length breadth ratio of original pattern.
Obtain continuous print in gas-liquid-solid heterogeneous mixed process in step of the present invention (1) stir real-time pattern and preserve, and utilize digital image processing techniques to obtain the pending pattern of only having powerful connections with target; This process is conventional method, such as, utilizes high-speed camera to stir real-time video to obtain continuous print in heterogeneous mixed process in prior art; Video software KMPlayer is utilized to be that real-time pattern preserves by Video Quality Metric.
In statistics, uniform experiment design only considers a kind of test design method of testing site uniformly dispersing in trial stretch; Theoretical according to uniform experiment design, the tolerance of homogeneity adopts number theory method (pseudo-monte carlo method), the most generally adopts deviation; What adopt herein is work as time corresponding deviation; At pilot region (this herein, ) upper cloth nindividual testing site , order , for in by initial point 0 to the rectangle determined.Order for in test count, when during middle distribution even test point, should be with volume be close; In brief, when counting in pilot region than when comparing closely with rectangular area, namely think that testing site is dispersed evenly to pilot region, we are referred to as polyphasic flow and mix.
The invention has the beneficial effects as follows:
(1) the method simple possible, has the uniform experiment design theory of statistics to support;
(2) weak point of electrical conductivity method, thermocouple method, optical method, decoloring method etc. is solved;
(3) economic loss that stirred reactor causes because of irrational design can be reduced.
Accompanying drawing explanation
Fig. 1 is process chart of the present invention.
Fig. 2 is that the present invention is by the upper left corner as the pattern stress and strain model figure of reference position; Filled black part is reference position (S=1).
Fig. 3 is that the present invention is by the lower left corner as the pattern stress and strain model figure of reference position; Filled black part is reference position (S=1).
Fig. 4 is that the present invention is by the lower right corner as the pattern stress and strain model figure of reference position; Filled black part is reference position (S=1).
Fig. 5 is that the present invention is by the upper right corner as the pattern stress and strain model figure of reference position; Filled black part is reference position (S=1).
Fig. 6 is of the present invention variation diagram in time.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail, but protection scope of the present invention is not limited to described content.
Embodiment 1
Use for reference CH 4with the polyphasic flow reaction model of ZnO in fuse salt, built gas top blast and stirred hydrodynamics experiment porch, when being under initial static state, the carrier of oxygen falls to bottom molten bath with solid granular form, when blasting reacting gas CH from reactor top blast 4after, under the effect that is uniformly mixed of gas, the carrier of oxygen can be dispersed in whole molten bath, forms reflection system as fluidized bed as, that is: whole reaction is that a gas-liquid-solid polyphasic flow reacts; For simplifying experimental provision and do not lose representativeness, choose material needed for experiment and reagent is mainly common water, N 2and diameter is the black polystyrene particle of 0.45 millimeter; The different depth immersed by the different flow and spray gun of considering top-blown gas designs multiple stirring condition, then utilizes high-speed camera to take the flowing pattern of solid granule under top-blown gas beating action.
Concrete implementation step is as follows:
(1) N 2sprayed from top to bottom by the glass tube spray gun of a vertical movement, spray gun caliber size immobilizes, but gas flow and spray gun immersion depth can for adjustment; The quality of water and the quality of solid particle are all certain, and wherein polystyrene black particle is the bottom being tiled in glass reactor when testing, and carry out mix and blend, the N in chemical reactor only by injection nitrogen as power 2, granules of polystyrene and water is exactly formed some special flowing patterns by such alr mode.
(2) use the continuous videos of the vertical section of high-speed camera collection shown in Fig. 1, use video software KMPlayer, catch down the pattern of each frame; Then utilize MATLAB to implement digital number image processing techniques, obtain pending image through greyscale transformation, cap transformation, filtering operation, binaryzation and opening operation process.
(2) pending pattern is carried out Region dividing, corresponding from the upper left corner , the lower left corner is corresponding , the lower right corner is corresponding , the upper right corner is corresponding carry out stress and strain model as reference position, divide 16 regions respectively, and each zonule is numbered ;
(3) use uniform experiment design theoretical, tolerance the method of moment mixing uniformity calculates as follows:
1. the number ratio with the target in overall pattern in a zonule is calculated:
.
The upper left corner ( ) number ratio is: 0.0010,0.0010,0.0042,0.0120,0.0345,0.0734,0.1182,0.1747,0.2457,0.3423,0.4432,0.5518,0.6640,0.7793,0.8997,1.0000;
The lower left corner ( ) number ratio is: 0.0002,0.00660.0202,0.0475,0.0883,0.1411,0.2059,0.2789,0.3579,0.4530,0.5648,0.6762,0.7798,0.8718,0.9480,1.0000;
The lower right corner ( ) number ratio is: 0.0018,0.0105,0.0295,0.0574,0.0964,0.1477,0.2193,0.3000,0.3944,0.4924,0.6003,0.7207,0.8208,0.9089,0.9565,1.0000;
The upper right corner ( ) number ratio is: 0.0007,0.0032,0.0142,0.0324,0.0644,0.1139,0.1756,0.2464,0.3262,0.4156,0.5006,0.6068,0.7057,0.8148,0.9063,1.0000;
2. the area ratio in this zonule and overall pattern is calculated:
.
The upper left corner ( ) area ratio is: 0.0039,0.0156,0.0352,0.0625,0.0977,0.1406,0.1914,0.2500,0.3164,0.3906,0.4727,0.5625,0.6602,0.7656,0.8789 and 1.0000;
The lower left corner ( ) area ratio is: 0.0039,0.0156,0.0352,0.0625,0.0977,0.1406,0.1914,0.2500,0.3164,0.3906,0.4727,0.5625,0.6602,0.7656,0.8789 and 1.0000;
The lower right corner ( ) area ratio is: 0.0039,0.0156,0.0352,0.0625,0.0977,0.1406,0.1914,0.2500,0.3164,0.3906,0.4727,0.5625,0.6602,0.7656,0.8789 and 1.0000;
The upper right corner ( ) area ratio is: 0.0039,0.0156,0.0352,0.0625,0.0977,0.1406,0.1914,0.2500,0.3164,0.3906,0.4727,0.5625,0.6602,0.7656,0.8789 and 1.0000;
3. target number is calculated than the absolute deviation with region pattern area ratio:
.
The upper left corner ( ) absolute deviation is: 0.0029,0.0146,0.0309,0.0505,0.0632,0.0672,0.0732,0.0753,0.0707,0.0483,0.0294,0.0107,0.0039,0.0136,0.0208,0;
The lower left corner ( ) area ratio is: 0.0037,0.0090,0.0150,0.0150,0.0093,0.0004,0.0145,0.0289,0.0415,0.0624,0.0921,0.1137,0.1197,0.1062,0.0691,0;
The lower right corner ( ) area ratio is: 0.0021,0.0051,0.0056,0.0051,0.0012,0.0070,0.0279,0.0500,0.0780,0.1017,0.1277,0.1582,0.1607,0.1433,0.0776,0;
The upper right corner ( ) area ratio is: 0.0032,0.0124,0.0210,0.0301,0.0333,0.0267,0.0159,0.0036,0.0098,0.0250,0.0279,0.0443,0.0455,0.0492,0.0274,0;
4. the mean value of absolute deviation is calculated: ;
After the computing carrying out 4 angles, calculating mean absolute deviation is: 0.0030,0.0103,0.0181,0.0252,0.0267,0.0253,0.0329,0.0394,0.0500,0.0594,0.0693,0.0817,0.0824,0.0781,0.0487,0;
5. the uniformity coefficient of moment pattern is: , wherein for median;
Final calculating the uniformity coefficient of moment pattern be 0.9639. 1. similar ~ 5. step, calculate the uniformity coefficient in other moment, partial data (data have 500) is as follows successively: 0.939357235916417, 0.909570362583631, 0.902455628679390, 0.934843059712128, 0.942116295370033, 0.913261325923245, 0.930138086538828, 0.897325073446152, 0.925889388977547, 0.936786099075402, 0.944416821238247, 0.914719842914128, 0.930100811248494, 0.949240266159665, 0.950410111103229, 0.953723600034418, 0.930788402060183, 0.951252482116207, 0.896746433478823, 0.859610547932562, 0.902868476854514, 0.873579129496495, 0.865606883523398, 0.863530648982136, 0.902804123755286, 0.904152042239641, 0.845595571706213, 0.883584062119069, 0.863651742650436, 0.846244221263190, 0.858945595253722, 0.855847575167351, 0.842705478700894, 0.845253026528311, 0.850729807829946, 0.876583288251411, 0.856507593084105, 0.849324879199439, 0.892330341732519, 0.850151657965971, 0.849603273762746, 0.857279523145904, 0.846530705828734,The critical point that uniformity coefficient is occurred in obviously fluctuation change is defined as incorporation time: the 168th frame.As shown in Figure 6.
(4) theoretical by uniform experiment design, in incorporation time tfront and back, obvious Changing Pattern can be presented; If experience sufficient incorporation time tafter, numerical values recited fluctuation be bordering on 1, if draw change curve, then exist tafter moment, just be in close proximity to 1 and in time backward fluctuation little, then represent heterogeneous being blended in tthe state mixed is reached after moment.

Claims (2)

1. judge a method for multi-phase mixing uniformity, it is characterized in that, specifically comprise the following steps:
(1) in the mixed process of transparent tank diameter, obtain continuous print in gas-liquid-solid heterogeneous mixed process and stir real-time pattern and preserve; And utilize digital image processing techniques to obtain the pending pattern of only having powerful connections with target; Using the target in pending pattern as research object;
(2) according to the size of pending original pattern, corresponding from the upper left corner , the lower left corner is corresponding , the lower right corner is corresponding , the upper right corner is corresponding carry out stress and strain model as reference position, divide respectively individual region for natural number, and each zonule is numbered ;
(3) use uniform experiment design theoretical, tolerance tthe method of moment mixing uniformity calculates as follows:
1. the number ratio with the target in overall pattern in a zonule is calculated: ;
2. the area ratio in this zonule and overall pattern is calculated: ;
3. target number is calculated than the absolute deviation with region pattern area ratio: ;
4. the mean value of absolute deviation is calculated: ;
5. tthe uniformity coefficient of moment pattern is: , wherein represent median;
(4) theoretical by uniform experiment design, in incorporation time tfront and back, obvious Changing Pattern can be presented; If experience sufficient incorporation time tafter, numerical values recited fluctuation be bordering on 1, if draw change curve, then exist tafter moment, just be in close proximity to 1 and in time backward fluctuation little, then represent heterogeneous being blended in tthe state mixed is reached after moment.
2. the method for judgement multi-phase mixing uniformity according to claim 1, is characterized in that, the length breadth ratio of step (2) small area is equal with the length breadth ratio of original pattern.
CN201510503871.XA 2015-08-17 2015-08-17 Method for judging multi-phase mixing uniformity Pending CN105067621A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510503871.XA CN105067621A (en) 2015-08-17 2015-08-17 Method for judging multi-phase mixing uniformity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510503871.XA CN105067621A (en) 2015-08-17 2015-08-17 Method for judging multi-phase mixing uniformity

Publications (1)

Publication Number Publication Date
CN105067621A true CN105067621A (en) 2015-11-18

Family

ID=54497045

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510503871.XA Pending CN105067621A (en) 2015-08-17 2015-08-17 Method for judging multi-phase mixing uniformity

Country Status (1)

Country Link
CN (1) CN105067621A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105910996A (en) * 2016-04-07 2016-08-31 昆明理工大学 Method for determining mixing time of multiphase flow by using RGB color model
CN106546591A (en) * 2016-11-07 2017-03-29 首钢总公司 A kind of detection method of sintering deposit uniformity
CN108072740A (en) * 2016-11-10 2018-05-25 内蒙古伊利实业集团股份有限公司 A kind of assay method of the mixture homogeneity of dry mixed powder product
CN109785342A (en) * 2019-01-15 2019-05-21 昆明理工大学 A kind of analysis method of the animal meat quality uniformity
CN109903243A (en) * 2019-02-20 2019-06-18 云南农业大学 A method of multiphase stirring and mixing effect is characterized based on Logistics model
CN110068575A (en) * 2019-04-29 2019-07-30 嘉兴学院 Judgment method, judgment means and the terminal of multi-phase mixing uniformity
CN112669268A (en) * 2020-12-21 2021-04-16 昆明理工大学 Method, system and terminal for evaluating distribution uniformity of multiphase mixed concentration field
CN113144999A (en) * 2021-04-15 2021-07-23 中联重科股份有限公司 Method and system for detecting stirring uniformity of viscous substance and stirring control method
CN114166829A (en) * 2021-12-08 2022-03-11 华中科技大学鄂州工业技术研究院 Slurry uniformity detection system and method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813641A (en) * 2010-05-11 2010-08-25 昆明理工大学 Method for verifying homogeneous state and degree of multiphase stirring and mixing
CN101822957A (en) * 2010-04-07 2010-09-08 昆明理工大学 Method for judging multiphase mixing uniformity based on statistics and ergodic theory
CN102032802A (en) * 2010-07-27 2011-04-27 昆明理工大学 Method for representing combustion non-uniformity characteristic of industrial furnace
CN102156085A (en) * 2010-12-06 2011-08-17 昆明理工大学 Method for detecting uniformity of multi-phase mixed flow field in chemical field

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101822957A (en) * 2010-04-07 2010-09-08 昆明理工大学 Method for judging multiphase mixing uniformity based on statistics and ergodic theory
CN101813641A (en) * 2010-05-11 2010-08-25 昆明理工大学 Method for verifying homogeneous state and degree of multiphase stirring and mixing
CN102032802A (en) * 2010-07-27 2011-04-27 昆明理工大学 Method for representing combustion non-uniformity characteristic of industrial furnace
CN102156085A (en) * 2010-12-06 2011-08-17 昆明理工大学 Method for detecting uniformity of multi-phase mixed flow field in chemical field

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宁建辉: "混料均匀试验设计", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 *
徐建新: "多相体系搅拌混合效果评价方法及其应用研究", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105910996A (en) * 2016-04-07 2016-08-31 昆明理工大学 Method for determining mixing time of multiphase flow by using RGB color model
CN106546591A (en) * 2016-11-07 2017-03-29 首钢总公司 A kind of detection method of sintering deposit uniformity
CN106546591B (en) * 2016-11-07 2019-12-31 首钢集团有限公司 Method for detecting uniformity of sinter
CN108072740A (en) * 2016-11-10 2018-05-25 内蒙古伊利实业集团股份有限公司 A kind of assay method of the mixture homogeneity of dry mixed powder product
CN108072740B (en) * 2016-11-10 2021-08-27 内蒙古伊利实业集团股份有限公司 Method for measuring mixing uniformity of dry-method mixed powder product
CN109785342A (en) * 2019-01-15 2019-05-21 昆明理工大学 A kind of analysis method of the animal meat quality uniformity
CN109903243A (en) * 2019-02-20 2019-06-18 云南农业大学 A method of multiphase stirring and mixing effect is characterized based on Logistics model
CN110068575A (en) * 2019-04-29 2019-07-30 嘉兴学院 Judgment method, judgment means and the terminal of multi-phase mixing uniformity
CN112669268A (en) * 2020-12-21 2021-04-16 昆明理工大学 Method, system and terminal for evaluating distribution uniformity of multiphase mixed concentration field
CN113144999A (en) * 2021-04-15 2021-07-23 中联重科股份有限公司 Method and system for detecting stirring uniformity of viscous substance and stirring control method
CN114166829A (en) * 2021-12-08 2022-03-11 华中科技大学鄂州工业技术研究院 Slurry uniformity detection system and method
CN114166829B (en) * 2021-12-08 2023-09-19 华中科技大学鄂州工业技术研究院 Slurry uniformity detection system and method

Similar Documents

Publication Publication Date Title
CN105067621A (en) Method for judging multi-phase mixing uniformity
Monchaux et al. Preferential concentration of heavy particles: A Voronoï analysis
Madan et al. Modeling of mixing in ladles fitted with dual plugs
CN101822957B (en) Method for judging multiphase mixing uniformity based on statistics and ergodic theory
Ling et al. Investigation on the effect of nozzle number on the recirculation rate and mixing time in the RH process using VOF+ DPM model
Zalc et al. Extensive validation of computed laminar flow in a stirred tank with three Rushton turbines
Hilton et al. Modelling spray coating using a combined CFD–DEM and spherical harmonic formulation
Szalai et al. Mixing analysis in a tank stirred with Ekato Intermig® impellers
CN101813641B (en) Method for verifying homogeneous state and degree of multiphase stirring and mixing
CN101650292B (en) Method for measuring mixing effect of fluid
Witz et al. Local gas holdup simulation and validation of industrial-scale aerated bioreactors
Chibwe et al. CFD modelling of global mixing parameters in a Peirce-Smith converter with comparison to physical modelling
Zhao et al. Gas-liquid mass transfer and flow phenomena in a peirce-smith converter: a numerical model study
Quiyoom et al. Optimization of bottom tuyere configuration for Basic Oxygen Furnace steelmaking through experiments and CFD simulations
CN103279433B (en) A kind of method of characterizing particles aggregate and device
Ansoni et al. Multiobjective optimization of a flat-panel airlift reactor designed by computational fluid dynamics
CN102156085A (en) Method for detecting uniformity of multi-phase mixed flow field in chemical field
Liu et al. Lagrangian particle tracking in mechanically agitated polydisperse suspensions: Multi-component hydrodynamics and spatial distribution
CN102004816A (en) Method for representing multiphase mixing effect by utilizing histogram
Xiao et al. Extraction and evolution of bubbles attributes in a two-phase direct contact evaporator
Hristov et al. A simplified CFD for three-dimensional analysis of fluid mixing, mass transfer and bioreaction in a fermenter equipped with triple novel geometry impellers
Jeng et al. Experimental and numerical study of transient natural convection due to mass transfer in inclined enclosures
CN115980038A (en) Gas-liquid two-phase flow bubble image acquisition and measurement device and measurement method thereof
CN102279091A (en) Test device and test method for simulating molten iron flow at furnace bottom of blast-furnace hearth
Wang et al. Bubble refinement by non-invasive rotating flow

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20151118

RJ01 Rejection of invention patent application after publication