CN103149581A - Computing method of displacement absorbed doses of space three-dimensional structure device - Google Patents

Computing method of displacement absorbed doses of space three-dimensional structure device Download PDF

Info

Publication number
CN103149581A
CN103149581A CN2013100249236A CN201310024923A CN103149581A CN 103149581 A CN103149581 A CN 103149581A CN 2013100249236 A CN2013100249236 A CN 2013100249236A CN 201310024923 A CN201310024923 A CN 201310024923A CN 103149581 A CN103149581 A CN 103149581A
Authority
CN
China
Prior art keywords
space
dimensional structure
displacement
absorbed dose
computing method
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
CN2013100249236A
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN2013100249236A priority Critical patent/CN103149581A/en
Publication of CN103149581A publication Critical patent/CN103149581A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Radiation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a computing method of displacement absorbed doses of a space three-dimensional structure device and relates to the computing method of displacement absorbed doses. The method solves the technical problems, existing in the present displacement absorbed doses, of being high in cost, bad in accuracy. The computing method comprises a step of determining a radiation-receiving energy spectrum of the space three-dimensional structure device; a step of calculating a distribution curve of absorbed doses following with depth of aluminum; a step of ensuring quality of all members of the space three-dimensional structure device; a step of ensuring calculating points of the members; a step of ensuring a plume in a subdivision mode by regarding the calculating points as poles of ball coordination; a step of equating materials on the plume to the aluminum and calculating a displacement absorbed dose valve of the plume; a step of adding up all the displacement absorbed doses values of the plume in order to gain the displacement absorbed doses value of the calculating points; a step of repeating from the fifth step to the seventh step in order to gain the displacement absorbed doses values of all the calculating points of the space three-dimensional structure device. The calculation of the displacement absorbed doses is conducted under an irradiation condition of space rail charge particles.

Description

The computing method of space three-dimensional structure device displacement absorbed dose
Technical field
The present invention relates to the computing method of displacement absorbed dose.
Background technology
Accurately calculate spacecraft and device at rail displacement absorbed dose, it is one of important means of estimating spacecraft and device lifetime/reliability, especially to the device of displacement damage sensitivity (as, solar cell, CCD device, electric light and photoelectric device etc.) all the more so.Comparatively deep in the calculating of rail ionization absorbed dose to spacecraft and device in the world at present.For example, be based on the earliest one dimension Monte Carlo method and develop and be applied to SHIELDOSE on engineering [1]Program.This program can the Computing One-Dimensional housing, the target absorbed dose of unlimited thickness of slab and the protection of limited thickness of slab, and wherein the material of protective layer mainly is thought of as Al, and the material of target can be Si, SiO 2, H 2O, Al etc.This program convenience of calculation, but can not carry out the calculating of 3 D complex structure absorbed dose, can not the displacement calculating absorbed dose.Subsequently, developed NOVICE, MCNP and the Geant4 supervisor of developing based on three-dimensional Monte Carlo method.The computing method of these programs can consider that the ionization absorbed dose of labyrinth calculates, but because the convergence of Monte Carlo computing method is poor especially, if so structure more complicated, quite consuming time, sometimes may all can not get result preferably by some months, this is very unrealistic on scientific research, especially engineering.These years recently, the U.S. and European countries have developed a kind of new ionization absorbed dose computing method---fan-shaped Meshing Method, at first this method is carried out the 3 D complex structure grid and is divided, then divide for every kind of grid, database in conjunction with Monte Carlo method is calculated calculates.The result of calculation of the method is comparatively accurate, and has greatly improved counting yield, is convenient to practical implementation.
Compare with the fast development of ionization absorbed dose computing method, the computing method development of 3 D complex displacement structure absorbed dose is relatively slow.In the calculating of rail displacement absorbed dose, be mainly the displacement damage of estimating spacecraft and device.What be used for the earliest the evaluation displacement damage is equivalent fluence method, and equivalent fluence method is the space service life assessment technique that is used for estimating monocrystalline silicon and unijunction solar cell that proposes in the phase at the beginning of the eighties in last century.This technology can be reinforced the damage of device resist displacement reference frame is provided.Its main shortcoming is: need a large amount of ground simulation tests, and test is to carry out in the mode of ex situ, this just need to provide a large amount of tested device spares, and the test assessment cost is high, and owing to being non-in site measurement, its accuracy is also relatively poor.Development subsequently be that the essence of the NIEL method in NRL laboratory is displacement absorbed dose method, but it does not consider that the displacement absorbed dose is with the distribution of the degree of depth.At present also be not specifically designed in the world the 3 D complex structure devices in the computing method of rail displacement absorbed dose.
Summary of the invention
The method cost that the present invention will solve existing displacement calculating absorbed dose is high, the technical matters of poor accuracy, and the computing method of space three-dimensional structure device displacement absorbed dose are provided.
The computing method of space three-dimensional structure device displacement absorbed dose of the present invention are carried out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
The present invention finds out the calculation level in the space three-dimensional structure instrument one by one, take each calculation level as the utmost point heart, calculate the equivalent thickness of aluminium of all materials that pass of polar curve of this utmost point heart, obtain the displacement absorbed dose of this utmost point heart place calculation level, method is simple, and accurately, the inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.
Description of drawings
Fig. 1 tests the structural representation of 3 d structure model figure in;
Fig. 2 is geostationary orbit in test one, utilizes the AE8 radiation belt electronic model of NASA to calculate acquisition radiation belt electronics differential energy spectrum;
Fig. 3 is geostationary orbit in test one, utilizes the AP8 radiation belt proton model of NASA to calculate acquisition radiation belt proton differential energy spectrum;
Fig. 4 is geostationary orbit in test one, utilizes the JPL91 model to calculate the track solar cosmic ray proton power spectrum that obtains;
Fig. 5 be the test one in the displacement absorbed dose with aluminium depth profile curve D~t.Wherein a is that the displacement absorbed dose of terrestrial radiation having electronic is with the relation curve of aluminium protective thickness; B is that the displacement absorbed dose of radiation belt of the earth proton is with the relation curve of aluminium protective thickness; C is that the displacement absorbed dose of solar cosmic ray proton is with the relation curve of aluminium protective thickness; D is that total displacement absorbed dose is with the relation curve with the aluminium protective thickness;
Embodiment
Embodiment one: present embodiment the computing method of space three-dimensional structure device displacement absorbed dose carry out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
Present embodiment is found out the calculation level in the space three-dimensional structure instrument one by one, take each calculation level as the utmost point heart, calculate the equivalent thickness of aluminium of all materials that pass of polar curve of this utmost point heart, obtain the displacement absorbed dose of this utmost point heart place calculation level, method is simple, the inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.
Embodiment two: what present embodiment was different from embodiment one is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 200~1000 parts, 0~2 π of φ is divided into 200~2000 parts.Other are identical with embodiment one.
Embodiment three: what present embodiment was different from embodiment one is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 250 parts, 0~2 π of φ is divided into 250 parts.Other are identical with embodiment one.
Embodiment four: what present embodiment was different from one of embodiment one to three is: in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision be: interval 0~π of θ is divided into 400 parts, 0~2 π of φ is divided into 300 parts.Other are identical with one of embodiment one to three.
With following verification experimental verification beneficial effect of the present invention:
Test one: the computing method of the space three-dimensional structure device displacement absorbed dose of test one are carried out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, and interval 0~π of θ is divided into 300 parts, 0~2 π of φ is divided into 300 parts, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
In this test the schematic diagram of space three-dimensional structure instrument as shown in Figure 1,
Determine orbit parameter according to the space three-dimensional structure device in rail military service task in this test procedure one, comprise orbit altitude, orbit inclination, excentricity etc., can spectrum model be calculated the power spectrum of the radiation that the track differential energy spectrum accepts by track.Fig. 2 is geostationary orbit, utilizes the AE8 radiation belt electronic model of NASA to calculate acquisition radiation belt electronics differential energy spectrum; Fig. 3 is geostationary orbit, utilizes the AP8 radiation belt proton model of NASA to calculate acquisition radiation belt proton differential energy spectrum; Fig. 4 is geostationary orbit, utilizes the JPL91 model to calculate the track solar cosmic ray proton power spectrum that obtains;
Take the power spectrum of step 1 gained as data source, utilize Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve in step 3, concrete outcome as shown in Figure 5.Wherein a is that the displacement absorbed dose of terrestrial radiation having electronic is with the relation curve of aluminium protective thickness; B is that the displacement absorbed dose of radiation belt of the earth proton is with the relation curve of aluminium protective thickness; C is that the displacement absorbed dose of solar cosmic ray proton is with the relation curve of aluminium protective thickness; D is that total displacement absorbed dose is with the relation curve with the aluminium protective thickness; Curve d is 10 at the thickness of protective layer Al -8~10 -3In time, overlap with curve b, and curve d is 2 * 10 at the thickness of protective layer Al -3~4 * 10 0In time, overlap with curve a.
For space three-dimensional structure instrument shown in Figure 1, define its inner required calculation level, concrete as " A " point in Fig. 1.As seen from the figure, the safeguard structure placement differences of " A " some surrounding is very large, calculates if use Equivalent A l spherical shell, and its error is very large.And its actual labyrinth situation is considered in the calculating of this test one, and result of calculation is comparatively accurate.For " A " calculation level, carry out three-dimensional ball subdivision.Carry out three-dimensional ball subdivision with around calculation level by position angle and polar angle, calculate the corresponding equivalent material of each little grid direction and density.
With the result of calculation combination in step 2 and six, obtain the displacement absorbed dose of the calculation level that defines in the 3 D complex structure.During concrete calculating, need each little grid integration in the hope of last result of calculation.In this example, result of calculation is 6rad (Si).Thus result also as can be known, its equivalent protective thickness is larger, greater than 1g/cm 2
Test in one be with electron device under geostationary orbit, utilize labyrinth displacement absorbed dose computing method, calculate its inner displacement absorbed dose.For other labyrinth, can realize identical target according to the process of the inventive method.The inventive method can be for all spacecraft military service tracks and complicated topology layout, and the displacement absorbed dose that carries out under the radiation environment of space orbit charged particle calculates.

Claims (4)

1. the computing method of space three-dimensional structure device displacement absorbed dose is characterized in that the method carries out according to the following steps:
One, press military service track and the military service period of space three-dimensional structure instrument, the power spectrum of the radiation that mensuration space three-dimensional structure device is accepted, the power spectrum of radiation is radiation belt of the earth proton power spectrum, radiation belt of the earth electronic energy spectrum and solar cosmic ray proton power spectrum;
The power spectrum of the radiation that two, obtains according to step 1 utilizes Monte-Carlo method displacement calculating absorbed dose with aluminium depth profile curve D~t;
Three, determine material, density and the composition of all members in the space three-dimensional structure instrument;
Four, determine the calculation level of member, this calculation level is the arbitrfary point in the space three-dimensional structure instrument;
Five, the limit take calculation level as spherical coordinates, determine spherical coordinate system, spherical coordinate system (r, θ, φ), get θ ∈ 0~π, φ ∈ 0~2 π carries out subdivision for the structure of space three-dimensional structure instrument to θ and φ, determines position and the quantity of polar curve r according to the cell surface of institute's subdivision;
Six, be aluminium with the equivalence of the material on each polar curve, and calculate the thickness t of equivalent aluminium, and according to the D in step 2~t curve, calculate the displacement absorbing agent value on each polar curve;
Seven, the displacement absorbing agent value on all polar curves is added and, obtain the displacement absorbed dose of this calculation level;
Eight, repeating step five to seven, obtain the displacement absorbing agent value of all calculation levels in the space three-dimensional structure instrument.
2. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 200~1000 parts, 0~2 π of φ is divided into 200~2000 parts.
3. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 250 parts, 0~2 π of φ is divided into 250 parts.
4. the computing method of space three-dimensional structure device displacement absorbed dose according to claim 1, it is characterized in that in step 5 for the structure of space three-dimensional structure instrument to the concrete grammar that θ and φ carry out subdivision being: interval 0~π of θ is divided into 400 parts, 0~2 π of φ is divided into 300 parts.
CN2013100249236A 2013-01-23 2013-01-23 Computing method of displacement absorbed doses of space three-dimensional structure device Pending CN103149581A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013100249236A CN103149581A (en) 2013-01-23 2013-01-23 Computing method of displacement absorbed doses of space three-dimensional structure device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013100249236A CN103149581A (en) 2013-01-23 2013-01-23 Computing method of displacement absorbed doses of space three-dimensional structure device

Publications (1)

Publication Number Publication Date
CN103149581A true CN103149581A (en) 2013-06-12

Family

ID=48547763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013100249236A Pending CN103149581A (en) 2013-01-23 2013-01-23 Computing method of displacement absorbed doses of space three-dimensional structure device

Country Status (1)

Country Link
CN (1) CN103149581A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404810A (en) * 2016-10-19 2017-02-15 哈尔滨工业大学 Displacement and radiation damage equivalence evaluation method for aromatic polymer insulating material
CN111581719A (en) * 2020-04-28 2020-08-25 北京天工科仪空间技术有限公司 Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding
CN113543615A (en) * 2021-06-29 2021-10-22 中国科学院长春光学精密机械与物理研究所 Irradiation protection method for space electronic equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2420861A1 (en) * 2010-08-20 2012-02-22 Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO A radiation dose meter for measuring radiation dose in an external magnetic field
CN102564931A (en) * 2010-12-31 2012-07-11 中国航天科技集团公司第五研究院第五一○研究所 Method for simulation test about comprehensive space radiation effect of surface function material for spacecraft

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2420861A1 (en) * 2010-08-20 2012-02-22 Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO A radiation dose meter for measuring radiation dose in an external magnetic field
CN102564931A (en) * 2010-12-31 2012-07-11 中国航天科技集团公司第五研究院第五一○研究所 Method for simulation test about comprehensive space radiation effect of surface function material for spacecraft

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李兴冀: "星用双极型器件带电粒子辐照效应及损伤机理", 《中国博士学位论文全文数据库信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404810A (en) * 2016-10-19 2017-02-15 哈尔滨工业大学 Displacement and radiation damage equivalence evaluation method for aromatic polymer insulating material
CN106404810B (en) * 2016-10-19 2019-04-09 哈尔滨工业大学 The shifted radiation equivalent damage evaluation method of aromatic polymer insulating materials
CN111581719A (en) * 2020-04-28 2020-08-25 北京天工科仪空间技术有限公司 Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding
CN111581719B (en) * 2020-04-28 2024-03-26 北京天工科仪空间技术有限公司 Radiation effect calculation method, device and equipment based on spacecraft three-dimensional shielding
CN113543615A (en) * 2021-06-29 2021-10-22 中国科学院长春光学精密机械与物理研究所 Irradiation protection method for space electronic equipment
CN113543615B (en) * 2021-06-29 2022-11-01 中国科学院长春光学精密机械与物理研究所 Irradiation protection method for space electronic equipment

Similar Documents

Publication Publication Date Title
Shikaze et al. Measurements of 0.2–20 GeV/n cosmic-ray proton and helium spectra from 1997 through 2002 with the BESS spectrometer
Mitrofanov et al. Dynamic albedo of neutrons (DAN) experiment onboard NASA’s Mars Science Laboratory
Cârloganu et al. Towards a muon radiography of the Puy de Dôme
Ackermann et al. Searches for cosmic-ray electron anisotropies with the Fermi Large Area Telescope
CN103116176A (en) Ground equivalent fluence calculating method for electronic component charged particle irradiation effect
Schraube et al. Experimental verification and calculation of aviation route doses
Velinov et al. Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere
CN103149581A (en) Computing method of displacement absorbed doses of space three-dimensional structure device
Mitrica et al. A mobile detector for measurements of the atmospheric muon flux in underground sites
Yáñez et al. A method to measure the integral vertical intensity and angular distribution of atmospheric muons with a stationary plastic scintillator bar detector
Samnøy et al. Microdosimetry with a 3D silicon on insulator (SOI) detector in a low energy proton beamline
CN112051601A (en) Virtual source principle-based source boundary parameter Monte Carlo inversion technology
Gustafsson et al. PHITS simulations of the Matroshka experiment
Tran et al. A novel silicon microdosimeter using 3D sensitive volumes: Modeling the response in neutron fields typical of aviation
CN110057843A (en) Radioactive source shipping container gamma shielding method for testing performance and system based on G (E) function
Kurochkin et al. Study of the radiation environment caused by galactic cosmic rays at flight altitudes, at the summit of the Zugspitze and at PTB Braunschweig
CN107545098B (en) Compound device space continuous spectrum proton displacement damage ground equivalent evaluation method
Oláh et al. Improvement of cosmic-ray muography for Earth sciences and civil engineering
Panuntun et al. Validation of HPGe efficiency modeling by MCNP5 for volume source
Ji et al. Assessment of dose rate of detected gamma emitting nuclides using a carborne survey with a large volume NaI (Tl) detector
Prokopovich et al. Comparison of SOI microdosimeter and tissue equivalent proportional counter measurements at the CERF facility
Chen et al. Monte Carlo simulation for background noise study of the Wudalianchi volcanic density imaging
Mitrica 20 years of cosmic muons research performed in IFIN-HH
Mitrica et al. Studies of the high energy muon flux in Unirea salt mine from Slanic Prahova
CN108710753A (en) Near space Neutron Environment emulation mode

Legal Events

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

Application publication date: 20130612