CN117093817A - Radiation transfer factor correction method for non-closed radiation heat exchange system - Google Patents

Radiation transfer factor correction method for non-closed radiation heat exchange system Download PDF

Info

Publication number
CN117093817A
CN117093817A CN202311361158.7A CN202311361158A CN117093817A CN 117093817 A CN117093817 A CN 117093817A CN 202311361158 A CN202311361158 A CN 202311361158A CN 117093817 A CN117093817 A CN 117093817A
Authority
CN
China
Prior art keywords
radiation
transfer factor
closed
heat exchange
exchange system
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
Application number
CN202311361158.7A
Other languages
Chinese (zh)
Other versions
CN117093817B (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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202311361158.7A priority Critical patent/CN117093817B/en
Publication of CN117093817A publication Critical patent/CN117093817A/en
Application granted granted Critical
Publication of CN117093817B publication Critical patent/CN117093817B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

The invention discloses a radiation transfer factor correction method for a non-closed radiation heat exchange system, which relates to the field of radiation heat exchange simulation methods.

Description

Radiation transfer factor correction method for non-closed radiation heat exchange system
Technical Field
The invention relates to the field of radiation heat exchange simulation methods, in particular to a radiation transfer factor correction method for a non-closed radiation heat exchange system.
Background
The Monte Carlo method has the characteristics of simple and clear physical model, flexible application and the like, and is widely applied to the field of heat radiation and heat exchange. However, due to the presence of random errors, the monte carlo method calculates the obtained radiation transfer factor (radiation transfer factor RD i,j Defined as the fraction of the own radiant energy of unit i that is ultimately absorbed by unit j after one or more reflections and scattering by the units of the system) in a radiation transmission system is often not guaranteed to satisfy the reciprocal property, and thus the iterative convergence problem that arises when performing temperature field calculations. The iterative bidirectional statistical Monte Carlo method can effectively reduce the reciprocal error of the radiation transfer factor, thereby effectively improving the solving precision of the heat radiation transmission problem.
The iterative bidirectional statistical Monte Carlo method utilizes the principle of reversible light path during heat radiation transmission to carry out multiple bidirectional statistics on Monte Carlo simulation of heat radiation transmission, fully utilizes the path information of energy beam propagation, improves the reciprocal property of radiation transfer factors, and is applied to numerous engineering problems. However, the iterative bidirectional statistical Monte Carlo method is only suitable for closed radiation heat exchange systems, and for non-closed systems, the reciprocity of the radiation transfer factor cannot be improved, but the calculation result of the radiation transfer factor is deteriorated, and non-physical results are brought.
Disclosure of Invention
In order to solve the problem of radiation transfer factor reciprocal optimization in a non-closed radiation heat exchange system, the invention provides a radiation transfer factor correction method for the non-closed radiation heat exchange system, which comprises the following steps:
step one: calculating to obtain the radiation transfer factor matrix RD of the non-closed radiation transmission system i,j The non-closed radiation transmission system comprises N units, and in the non-closed radiation transmission system, the radiation power coefficients corresponding to the N units are expressed as a set E, E= { E i I=1, 2,3, N, for RD i,j And E, i=1, 2,3,..n, j=1, 2,3,..n, E i The ith radiant power coefficient in set E;
step two: adding a virtual blackbody unit into the non-closed radiation transmission system to enable the non-closed radiation transmission system to be changed into a closed radiation heat exchange system, wherein a radiation transfer factor matrix of the closed radiation heat exchange system is as followsThe closed radiation heat exchange system comprises n+1 units, and in the closed radiation heat exchange system, the radiation power coefficient corresponding to the n+1 units is expressed as a set E ', E' = { }>I=1, 2,3,..n+1 }, for +.>And E', i=1, 2,3,..n+1, j=1, 2,3,..n+1;
step three: RD-based i,j Calculating to obtain the radiation transfer factor matrix corresponding to the ith row and the (n+1) th column
Step four: based onAnd E is i Calculating to obtain radiation power coefficient of virtual unit>
Step five: based on、E i And->Calculation to get->,/>A radiation transfer factor matrix corresponding to the (n+1) th row and the (i) th column;
step six: based on the ith radiant power coefficient in set EThe j-th radiant power coefficient in set EThe kth radiation power coefficient in set E +.>Radiation transfer factor matrix corresponding to ith row and jth column obtained by n-1 th bidirectional statistics>Radiation transfer factor matrix corresponding to jth row and ith column obtained by n-1 th bidirectional statisticsRadiation transfer factor matrix corresponding to kth row and ith column obtained by n-1 th bi-directional statistics +.>Calculated to obtainWherein,/>A radiation transfer factor matrix obtained for the nth bi-directional statistic, k=1, 2,3,..n+1;
step seven: removing the virtual blackbody unit from the non-closed radiation transmission system, and obtaining a corrected radiation transfer factor matrix of the non-closed radiation heat exchange system based on the radiation transfer factor matrix obtained by the nth bidirectional statistics
The invention converts the non-closed system into the closed system by adding the virtual blackbody unit, so as to solve the inapplicability problem of the iterative bidirectional statistical Monte Carlo method under the non-closed system and improve the reciprocal property of the radiation transfer factor of the non-closed radiation heat exchange system.
In some embodiments, the third step calculates the radiation transfer factor matrix corresponding to the ith row and the (n+1) th column by using the following formula
In some embodiments, the fourth step calculates the radiation power coefficient of the virtual cell using the following formula
In some embodiments, the fifth step calculates the radiation transfer factor matrix corresponding to the (n+1) th row and the (i) th column using the following formula
In some embodiments, the iterative formula in the step six is:
wherein,radiation transfer factor matrix obtained for nth bi-directional statistics,>radiation transfer factor matrix corresponding to ith row and jth column obtained for nth-1 bi-directional statistics,/th row and jth column>Radiation transfer factor matrix corresponding to the j-th row and i-th column obtained for the n-1 th bi-directional statistics,/th row>The radiation transfer factor matrix corresponding to the kth row and the ith column obtained for the nth-1 bi-directional statistics, k=1, 2, 3.
In some embodiments, the radiation transfer factor of the non-closed radiation heat exchange system is calculated in the seventh step by using the following formula
The one or more technical schemes provided by the invention have at least the following technical effects or advantages:
according to the characteristics of the radiation heat exchange system, the virtual blackbody unit is added to complement the non-closed radiation heat exchange system, so that the non-closed radiation heat exchange system is formed into a closed system without changing the original radiation heat exchange process, and then the iterative bidirectional statistical Monte Carlo method of the closed system is utilized for correction, so that the problem that the iterative bidirectional statistical Monte Carlo method is not applicable under the non-closed system is solved, and the reciprocal property of the radiation transfer factor of the non-closed radiation heat exchange system is improved.
Drawings
The accompanying drawings, which are included to provide a further understanding of embodiments of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention;
FIG. 1 is a flow chart of a radiation transfer factor correction method for a non-closed radiation heat exchange system;
FIG. 2 is a schematic diagram of a square cavity surface radiant heat exchange system.
Detailed Description
In order that the above-recited objects, features and advantages of the present invention will be more clearly understood, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description. In addition, the embodiments of the present invention and the features in the embodiments may be combined with each other without collision.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, however, the present invention may be practiced in other ways than within the scope of the description, and the scope of the invention is therefore not limited to the specific embodiments disclosed below.
Example 1
According to the embodiment, a virtual unit is added on the basis of the non-closed radiation heat exchange system, the virtual unit is set to be a blackbody unit, the virtual unit can be regarded as a notch set of the non-closed system to form a new closed radiation heat exchange system, then an iterative bidirectional statistical Monte Carlo method is adopted to correct the radiation transfer factor of the closed radiation heat exchange system, and finally the added virtual unit is removed to obtain the radiation transfer factor of the corrected non-closed radiation heat exchange system.
Referring to fig. 1, fig. 1 is a schematic flow chart of a radiation transfer factor correction method for a non-closed radiation heat exchange system, and the specific steps implemented in the present invention are:
step one: for a non-closed radiation delivery system consisting of N cells,the radiation transfer factor matrix calculated by the traditional Monte Carlo method is recorded as RD i,j Where i=1, 2,3,..n, j=1, 2,3,..n. The radiation power coefficients corresponding to the N units are represented as a set { E } i I=1, 2,3, N, for the face unit, the face unit is provided with a plurality of face units,wherein->For emissivity of face unit i, S i For the area of the face unit i, +.>Wherein->For emissivity of volume element i, V i For the volume of the body unit i, wherein the conventional monte carlo method can refer to "iterative bidirectional statistical monte carlo method of heat radiation problem", engineering thermophysics report 2017, 38 (3): 635-639, formulas 1 to 2.
Step two: a virtual blackbody unit is added on the basis of the original non-closed radiation heat exchange system to form a new closed radiation heat exchange system, and the radiation transfer factor matrix of the closed system is recorded asThe corresponding radiation power coefficients of n+1 units are denoted as set { +.>I=1, 2,3,..n+1 }, where i=1, 2,3,..n+1, j=1, 2,3,..n+1. For i=1, 2, 3..n, j=1, 2, 3..n, there is +.>=RD i,j ,/>=E i The method comprises the steps of carrying out a first treatment on the surface of the For i=n+1, j=1, 2,3,..n+1, there are
Among them, a blackbody is an idealized object in the field of radiant heat exchange, capable of absorbing all external electromagnetic radiation (thermal radiation is one of electromagnetic radiation), and without any reflection and transmission. In other words, the absorption coefficient of the black body for electromagnetic waves of any wavelength is 1, and the transmission coefficient is 0.
Step three: calculating the radiation transfer factor of the (n+1) th column:
step four: calculating the radiation power coefficient of the virtual unit:
step five: calculating the radiation transfer factor of the (n+1) th column:
step six: and carrying out n times of bidirectional statistics by using an iterative bidirectional statistics Monte Carlo method, wherein the iterative formula is as follows:
the iterative bidirectional statistical Monte Carlo method can refer to the iterative bidirectional statistical Monte Carlo method of heat radiation problem, and engineering thermophysics report 2017, 38 (3): 635-639.
Step seven: removing the virtual unit to obtain the radiation transfer factor of the modified non-closed radiation heat exchange system
According to the method, the non-closed radiation heat exchange system can be converted into the closed radiation heat exchange system by adding the virtual blackbody unit, so that the problem that the iterative bidirectional statistics Monte Carlo method is not applicable under the non-closed system is solved, and the easiness of radiation transfer factors of the non-closed radiation heat exchange system is improved.
The method can be applied to solving the heat radiation transmission problem so as to improve the solving precision.
Example two
Based on the first embodiment, a square two-dimensional cavity surface radiation heat exchange system is selected in the second embodiment, the side length of the square two-dimensional cavity surface radiation heat exchange system is unit length, each side length is divided into 1 unit, and the surface emissivity of each unit is 1.0, as shown in fig. 2. 1-4 in FIG. 2 represent units 1-4, respectively, assuming that the unit 4 void forms a non-closed radiant heat exchange system comprised of 3 units.
Setting the number of analog energy beams of each unit to 100000, taking radiation transfer factor between the No. 1 unit and the No. 2 unit as an example, wherein the theoretical value of the radiation transfer factor is 0.29289322,1 unit and the radiation transfer factor matrix between the No. 2 unit is RD 1,2, The radiation transfer factor matrix between the No. 2 unit and the No. 1 unit is RD 2,1 RD is calculated by adopting the traditional Monte Carlo method 1,2 0.29362, RD 2,1 A reciprocal error between units 0.29160,1 and 2 ofThe reciprocal error between unit No. 2 and unit No. 1 is +.>The reciprocal error is +.>Wherein the reciprocal error between the i-th unit and the j-th unit is defined as +.>The reciprocal error between the j-th unit and the i-th unit is defined as +.>,/>Is->Is a symmetric matrix of:
the radiation transfer factor matrix of the non-closed radiation transmission system is RD i,j ,RD j,i Is RD i,j Symmetric matrix of E j The j-th radiant power coefficient in set E; after a virtual unit is added to form a closed system, the radiation transfer factor matrix between the modified No. 1 unit and the modified No. 2 unit is as follows after 20 times of iteration bidirectional statisticsThe radiation transfer factor matrix between the modified cell number 2 and cell number 1 is +.>Obtain->0.293189->0.293189, error in easiness of inversion->
The method of the patent can be seen to effectively improve the reciprocal property of the radiation transfer factor of the non-closed radiation heat exchange system, and ensure the calculation precision of the radiation transfer factor.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. It is therefore intended that the following claims be interpreted as including the preferred embodiments and all such alterations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.

Claims (6)

1. A method of modifying a radiation transfer factor for a non-closed radiation heat exchange system, the method comprising:
step one: calculating to obtain the radiation transfer factor matrix RD of the non-closed radiation transmission system i,j The non-closed radiation transmission system comprises N units, and in the non-closed radiation transmission system, the radiation power coefficients corresponding to the N units are expressed as a set E, E= { E i I=1, 2,3, N, for RD i,j And E, i=1, 2,3,..n, j=1, 2,3,..n, E i The ith radiant power coefficient in set E;
step two: adding a virtual blackbody unit into the non-closed radiation transmission system to enable the non-closed radiation transmission system to be changed into a closed radiation heat exchange system, wherein a radiation transfer factor matrix of the closed radiation heat exchange system is as followsThe closed radiation heat exchange system comprises n+1 units, and in the closed radiation heat exchange system, the radiation power coefficient corresponding to the n+1 units is expressed as a set E ', E' = { }>I=1, 2,3,..n+1 }, for +.>And E', i=1, 2,3,..n+1, j=1, 2,3,..n+1;
step three: RD-based i,j Calculating to obtain the radiation transfer factor matrix corresponding to the ith row and the (n+1) th column
Step four: based onAnd E is i Calculating to obtain radiation power coefficient of virtual unit>
Step five: based on、E i And->Calculation to get->,/>A radiation transfer factor matrix corresponding to the (n+1) th row and the (i) th column;
step six: based on the ith radiant power coefficient in set EThe j-th radiant power coefficient in set E->The kth radiation power coefficient in set E +.>Radiation transfer factor matrix corresponding to ith row and jth column obtained by n-1 th bidirectional statistics>Radiation transfer factor matrix corresponding to the jth row and ith column obtained by n-1 th bidirectional statistics>Radiation transfer factor matrix corresponding to kth row and ith column obtained by n-1 th bi-directional statistics +.>Calculation to get->Wherein->A radiation transfer factor matrix obtained for the nth bi-directional statistic, k=1, 2,3,..n+1;
step seven: removing the virtual blackbody unit from the non-closed radiation transmission system, and obtaining a corrected radiation transfer factor matrix of the non-closed radiation heat exchange system based on the radiation transfer factor matrix obtained by the nth bidirectional statistics
2. The method for correcting radiation transfer factor of non-closed radiation heat exchange system according to claim 1, wherein said step three calculates the radiation transfer factor matrix corresponding to the ith row and the n+1th column by using the following formula
3. The method for correcting radiation transfer factor of non-closed radiation heat exchange system according to claim 1, wherein said step four is to calculate the radiation power coefficient of the obtained virtual unit by using the following formula
4. The method for correcting radiation transfer factor of non-closed radiation heat exchange system according to claim 1, wherein the fifth step calculates the radiation transfer factor matrix corresponding to the n+1th row and the i-th column by using the following formula
5. The method for correcting radiation transfer factor of non-closed radiation heat exchange system according to claim 1, wherein the iterative formula in the step six is:
6. the method for correcting radiation transfer factor of non-closed radiation heat exchange system according to claim 1, wherein in the seventh step, the radiation transfer factor matrix of the non-closed radiation heat exchange system is obtained by calculation using the following formula
CN202311361158.7A 2023-10-20 2023-10-20 Radiation transfer factor correction method for non-closed radiation heat exchange system Active CN117093817B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311361158.7A CN117093817B (en) 2023-10-20 2023-10-20 Radiation transfer factor correction method for non-closed radiation heat exchange system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311361158.7A CN117093817B (en) 2023-10-20 2023-10-20 Radiation transfer factor correction method for non-closed radiation heat exchange system

Publications (2)

Publication Number Publication Date
CN117093817A true CN117093817A (en) 2023-11-21
CN117093817B CN117093817B (en) 2023-12-22

Family

ID=88780190

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311361158.7A Active CN117093817B (en) 2023-10-20 2023-10-20 Radiation transfer factor correction method for non-closed radiation heat exchange system

Country Status (1)

Country Link
CN (1) CN117093817B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050088440A1 (en) * 2003-10-22 2005-04-28 Microsoft Corporation Hardware-accelerated computation of radiance transfer coefficients in computer graphics
CN106503425A (en) * 2016-10-10 2017-03-15 北京信息科技大学 A kind of method for calculating dependent radiance transfer factor
CN115964804A (en) * 2022-12-26 2023-04-14 西安航天动力研究所 Numerical simulation method for radiation heat flow of solid wall surface of liquid rocket engine assembly
CN116305892A (en) * 2023-03-07 2023-06-23 西安航空学院 Monte Carlo method-based non-gray gas radiation heat transfer prediction method, system, equipment and medium
CN116883607A (en) * 2023-09-06 2023-10-13 四川物通科技有限公司 Virtual reality scene generation system based on radiation transmission

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050088440A1 (en) * 2003-10-22 2005-04-28 Microsoft Corporation Hardware-accelerated computation of radiance transfer coefficients in computer graphics
CN106503425A (en) * 2016-10-10 2017-03-15 北京信息科技大学 A kind of method for calculating dependent radiance transfer factor
CN115964804A (en) * 2022-12-26 2023-04-14 西安航天动力研究所 Numerical simulation method for radiation heat flow of solid wall surface of liquid rocket engine assembly
CN116305892A (en) * 2023-03-07 2023-06-23 西安航空学院 Monte Carlo method-based non-gray gas radiation heat transfer prediction method, system, equipment and medium
CN116883607A (en) * 2023-09-06 2023-10-13 四川物通科技有限公司 Virtual reality scene generation system based on radiation transmission

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
BAO-HAI GAO 等: "Fast reconstructing two-dimensional temperature distribution in participating media with different surfaces conditions", INFRARED PHYSICS & TECHNOLOGY, vol. 103, pages 103080 *
单喆简: "试论灰体表面的间幅射换热问题――与国内某些《传热学》的编著者商榷", 浙江理工大学学报, no. 03, pages 50 - 54 *
夏新林;任德鹏;郭亮;: "求解介质内热辐射传递的双向统计蒙特卡罗法", 工程热物理学报, no. 2, pages 21 - 24 *
帅永;董士奎;刘林华;谈和平;: "双向蒙特卡罗法模拟参与性介质中的辐射换热", 工程热物理学报, no. 01, pages 109 - 111 *
朱言旦;曾磊;杜雁霞;桂业伟;: "热辐射问题的迭代双向统计蒙特卡罗方法", 工程热物理学报, no. 03, pages 635 - 639 *
董士奎, 帅永, 谈和平, 刘林华: "反向蒙特卡罗法模拟参与性介质中热辐射传递", 哈尔滨工业大学学报, no. 12, pages 1602 - 1604 *
阮立明: "蒙特卡洛法求二维矩形散射性介质内的辐射传递", 燃烧科学与技术, vol. 8, no. 05, pages 390 - 394 *
齐英杰;马洪斌;杨春梅;: "多重蒙特卡罗算法在木片热磨过程中的应用初探", 林产工业, no. 03, pages 29 - 31 *

Also Published As

Publication number Publication date
CN117093817B (en) 2023-12-22

Similar Documents

Publication Publication Date Title
US9292634B2 (en) Compensated molds for manufacturing ion exchange strengthened, 3D glass covers
Choi et al. Modeling line profiles of protostellar collapse in B335 with the Monte Carlo method
CN111709096B (en) Design method of special-shaped fin structure for strengthening natural convection heat transfer
CN104408527B (en) Tower type solar heat and power system Jing Chang focusing strategy optimization method
CN117093817B (en) Radiation transfer factor correction method for non-closed radiation heat exchange system
CN101751697A (en) Three-dimensional scene reconstruction method based on statistical model
CN106839474A (en) A kind of porous media solar heat absorber and its method for designing
CN109117553B (en) Light beam energy distribution determination method suitable for equal-intensity laser impact
CN112231956A (en) Light-heat-force integrated coupling calculation method for tower type photo-thermal power station heat absorber
Johnson et al. Relativistic configuration-interaction calculation of the polarizabilities of heliumlike ions
He et al. A variable-parameter ZNN with predefined-time convergence for dynamic complex-valued Lyapunov equation and its application to AOA positioning
Gao et al. Polarized deep inelastic and elastic scattering from gauge/string duality
CN105303046B (en) A kind of method for reducing reactor by component periphery rod power error in rod calculating
CN113742864B (en) Synergistic control method for curing deformation of composite material based on global compensation quantity
CN114035300A (en) Adaptive correction method based on zoom lens under high-capacity space optical communication link
CN112651170B (en) Efficient characteristic contribution assessment method in longitudinal federal learning scene
CN114626573A (en) Load prediction method for optimizing extreme learning machine based on improved multivariate universe algorithm
Sensui et al. Evolution of clusters of galaxies: mass stripping from galaxies and growth of common halos
CN113252493B (en) Control method of thermal strength test system
CN207198460U (en) Speculum torque tube arrangement under heat flow density analysis method optimal solution
CN109241545A (en) Oval thermal-collecting tube heat flow density analysis method
CN114709286B (en) Photovoltaic module for realizing light transmission and collection of photovoltaic cell gaps
Safavinejad et al. Optimal Boundary Design of Radiant Enclosures Using Micro-Genetic Algorithm (Effects of Refractory Properties and Aspect Ratio of Enclosure on Heaters Setting)
CN113947002A (en) Method for calculating convection heat transfer coefficient of spectrum beam combination semiconductor laser
Chunfeng Boundary corrected cubic smoothing splines

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