CN100550028C - The method of satellite sun array surface charging numerical simulation prediction - Google Patents
The method of satellite sun array surface charging numerical simulation prediction Download PDFInfo
- Publication number
- CN100550028C CN100550028C CNB2007100384389A CN200710038438A CN100550028C CN 100550028 C CN100550028 C CN 100550028C CN B2007100384389 A CNB2007100384389 A CN B2007100384389A CN 200710038438 A CN200710038438 A CN 200710038438A CN 100550028 C CN100550028 C CN 100550028C
- Authority
- CN
- China
- Prior art keywords
- satellite
- array surface
- sun
- model
- numerical simulation
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of method of satellite sun array surface charging numerical simulation prediction, it comprises the definite residing space environment parameter of solar battery array, sun battle array surfacing characteristic and structural parameters, sets up the satellite sun array surface charging numerical simulation model, model is carried out numerical simulation calculation, sets up five steps of satellite sun array surface charging potential diagram; It is based on dynamic dispatching gas ions model, characteristic in conjunction with satellite structure and material surface, utilize the PIC method, designed the modeling algorithm of particle deposition evaluation algorithm and particle-electric current combination, spacecraft structure and the charged environmental interaction of space plasma to complexity carry out accurate numerical analysis, save satellite Development and design cost greatly, can be the satellite protection design direct reference data is provided, be highly suitable for the charging numerical analysis of the various structured materials in spacecraft surface.
Description
Technical field
The present invention relates to the space application technology, specifically a kind of method of satellite sun array surface charging numerical simulation prediction.
Background technology
Space sun battle array primary power supply system is as one of most important useful load of spacecraft, because it is exposed in the space plasma environment fully, be subjected to the influence of spatial charging environment easily and produce the problem of surface charging, cause spacecraft and the interactional effect that discharges and recharges of space plasma.At present in the world aspect spacecraft and the interactional numerical simulation of space plasma, all adopted advanced unit particle (PIC) method as its numerical simulation calculation method, unit particle (PIC) numerical analysis method is a kind of powerful numerical computation method based on dynamic dispatching gas ions model, it can carry out analytical calculation to the interaction between spacecraft and the space plasma more accurately, particularly the incomparable advantage of other similar numerical methods is arranged handling the target side mask that space sun battle array has labyrinth like this.ESA/TOS-EMA has developed the processing spacecraft plasma interaction and the galvanomagnetic effect software of LEO/PEO (Polar Earth Orbit) track, and wherein numerical evaluation adopts two-dimentional PIC method.By European IRF-K, CNRS-UVSQ/CETP and common PicUp3D/Spis spacecraft and the plasma interaction simulation softward of developing of ESA/TOS-EMA, it is based on 3D PIC (three-dimensional element particle method), in describing, the spacecraft geometric configuration used non-structure graph, the accurate numerical simulation calculation that can be used for spacecraft and plasma interaction, though above simulation method can satisfy computational accuracy, but time, the process of analog computation are oversize, are unfavorable for further carrying out the research work of sun battle array surface charging aspect, space.
Summary of the invention
The method that a kind of satellite sun array surface charging numerical simulation prediction of providing at the deficiencies in the prior art is provided, it utilizes the PIC method, according to the influence of satellite structure and material surface characteristic to surface charging, particle deposition evaluation algorithm and particle-electric current combination algorithm have been designed, and the PIC algorithm optimized, satisfy computational accuracy can be under requiring with the analog computation time decreased to original one tens to several percent.The present invention can directly apply to the research of satellite sun array surface charging rule, has saved satellite Development and design cost.
The technical scheme that realizes the object of the invention is: a kind of method of satellite sun array surface charging numerical simulation prediction, characteristics comprise the following steps:
A. determine electron number density, electron temperature, the ion number density in satellite sun array space of living in, the environmental parameter of ion temperature;
B. determine the maximum secondary electron emission coefficiency δ on vertical incidence satellite sun array surface
m, ENERGY E
mMaterial characteristic parameter and the structural parameters of surfacing present position, thickness and area;
C. based on dynamic plasma model, virtual space plasma and sun battle array structure are set up the satellite sun array surface charging numerical model, and this model resolved into many little grid cells, and determine electric density and electric current thus, obtain the Potential distribution rule on sun battle array surface;
D. the model to above-mentioned foundation carries out numerical simulation calculation, and this calculating comprises:
(i) determine the initial position and the velocity distribution of virtual space plasma particle according to environmental parameter, the time step of computing time and cycle calculations is set;
(ii) Maxwell (Maxwell) system of equations that electromagnetic field develops is described in utilization, obtains the electric field and the magnetic field at each grid cell place, obtains the suffered Lorentz force of each space plasma particle by the Lorentz force formula then;
(iii), obtain this particle's velocity and distribute and space distribution according to the suffered power and the moving situation of space plasma particle of simulation;
(iv) circulation is carried out above-mentioned steps computing (ii)～(iii) till reaching the schedule time, and exports the computational data of each time period inner potential.
E. use MATLAB software that the computational data of above-mentioned output is set up potential diagram, and satellite sun array surface charging numerical value is carried out simulation predicting.
The present invention is based on dynamic dispatching gas ions model, in conjunction with satellite structure and material surface characteristic, utilize the PIC method, designed the modeling algorithm of particle deposition evaluation algorithm and particle-electric current combination, it can carry out accurate numerical analysis to the spacecraft structure and the charged environmental interaction of space plasma of complexity, controlled simulated time effectively, saved satellite Development and design cost greatly, can be the satellite protection design direct reference data is provided, be highly suitable for the charging numerical analysis of the various structured materials in spacecraft surface.
Description of drawings
Fig. 1 is a sun array surface charging potential diagram of the present invention
Embodiment
Consult Fig. 1, the present invention is undertaken by following step:
1, determines the residing space environment parameter of solar battery array
Carry out numerical simulation at the sun array surface charging current potential in the geostationary orbit, adopting abominable its electron temperature of substorm environmental parameter of GEO track is 12keV, and electron density is 1.12 * 10
6/ m
3, ion temperature is 29.5keV, ion concentration is 0.236 * 10
6/ m
3
2, determine solar battery array surfacing characteristic, sun battle array structural parameters
The solar array surfacing comprises cover glass and polyimide (kapton) substrate.Cover glass maximum perpendicular secondary electron yield is set at 5.6, and the corresponding energy of maximum perpendicular secondary electron yield is 800eV; Kapton maximum secondary electron emission coefficiency is 1.7, and the corresponding energy of maximum perpendicular secondary electron yield is 300eV.The cover glass area is 0.019 * 0.019m
2Clearance distance is 0.001m between the cover glass, and clearance material is kapton.
3, set up the satellite sun array surface charging numerical simulation model
According to GEO orbital environment characteristics, sun array surface charging numerical simulation adopts the static model.Set up the satellite sun array surface charging numerical simulation model according to the residing space environment of solar battery array and solar battery array surfacing characteristic, sun battle array structural parameters.Model is based on dynamic plasma model, and virtual space plasma and sun battle array structure also are broken down into many little grid cells.The charged particle that comprises varying number in each grid cell, determine the quantity and the movement locus thereof of motion particle in the grid according to the electric field change rule of each grid cell concrete structure, surfacing characteristic and background plasma environment, and determine electric density and electric current thus, finally obtain the sun battle array surface potential regularity of distribution.
4, the model of setting up is carried out numerical simulation calculation
Come the space plasma particle is carried out motion simulation by following the tracks of the motion in electric field of a large amount of electronics and ion, the numerical simulation calculation process is:
(a) initialization: determine material parameter (comprising secondary electron emission characteristic etc.) and boundary condition, determine simulation particle initial position and velocity distribution, the time step of computing time and cycle calculations is set according to environmental parameter.
(b) solving equation: by the Maxwell system of equations, obtain the electric field at each grid cell place, obtain the suffered Lorentz force of each particle by the Lorentz force formula then;
(c) the stressed and motion of simulation particle:, obtain particle rapidity and distribute and space distribution by simulation particle suffered power and moving situation.
(d) finish computing: judge whether if time is up, then to finish to the schedule time, otherwise, continue loop computation, till reaching the schedule time, the result that output calculates, the result comprises the current potential that calculates in each time period.
5, set up the sun array surface charging potential diagram
Utilization MATLAB software is put in order the computational data of output, sets up sun array surface charging current potential figure, according to potential diagram the satellite sun array surface charging rule is made simulation predicting.
Claims (1)
1, a kind of method of satellite sun array surface charging numerical simulation prediction is characterized in that comprising the following steps:
A. determine electron number density, electron temperature, the ion number density in satellite sun array space of living in, the environmental parameter of ion temperature;
B. determine the maximum secondary electron emission coefficiency δ on vertical incidence satellite sun array surface
mWith the corresponding ENERGY E of vertical incidence maximum secondary electron emission coefficiency
mMaterial characteristic parameter and the structural parameters of satellite sun array surfacing present position, thickness and area;
C. based on dynamic plasma model, virtual space plasma and sun battle array structure are set up the satellite sun array surface charging numerical model, and this model resolved into many little grid cells, and determine electric density and electric current thus, obtain the Potential distribution rule on sun battle array surface;
D. the model to above-mentioned foundation carries out numerical simulation calculation, and this calculating comprises:
(i) determine the initial position and the velocity distribution of virtual space plasma particle according to environmental parameter, the time step of computing time and cycle calculations is set;
(ii) Maxwell (Maxwell) system of equations that electromagnetic field develops is described in utilization, obtains the electric field and the magnetic field at each grid cell place, obtains the suffered Lorentz force of each space plasma particle by the Lorentz force formula then;
(iii), obtain this particle's velocity and distribute and space distribution according to the suffered power and the moving situation of space plasma particle of simulation;
(iv) circulation is carried out above-mentioned steps computing (ii)～(iii) till reaching the schedule time, and exports the computational data of each time period inner potential;
E. use MATLAB software that the computational data of above-mentioned output is set up the charging potential figure on satellite sun array surface, and satellite sun array surface charging numerical value is carried out simulation predicting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100384389A CN100550028C (en) | 2007-03-26 | 2007-03-26 | The method of satellite sun array surface charging numerical simulation prediction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100384389A CN100550028C (en) | 2007-03-26 | 2007-03-26 | The method of satellite sun array surface charging numerical simulation prediction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101276378A CN101276378A (en) | 2008-10-01 |
| CN100550028C true CN100550028C (en) | 2009-10-14 |
Family
ID=39995819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2007100384389A Expired - Fee Related CN100550028C (en) | 2007-03-26 | 2007-03-26 | The method of satellite sun array surface charging numerical simulation prediction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100550028C (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101691140B (en) * | 2009-10-12 | 2012-08-22 | 浙江大学 | Pico-satellite solar cell simulator and simulation method |
| CN102117345B (en) * | 2009-12-31 | 2012-10-31 | 北京卫星环境工程研究所 | Method for surface charging analysis of satellite |
| CN102564931A (en) * | 2010-12-31 | 2012-07-11 | 中国航天科技集团公司第五研究院第五一○研究所 | Method for simulation test about comprehensive space radiation effect of surface function material for spacecraft |
| CN102928714B (en) * | 2012-11-02 | 2014-08-13 | 北京航空航天大学 | Moonlet sun array life forecast method based on I-V curve and energy balance |
| CN102981074A (en) * | 2012-11-20 | 2013-03-20 | 中国航天科技集团公司第五研究院第五一〇研究所 | Interior charging and discharging characteristic stimulation test system and method of high-power part |
| CN103761359A (en) * | 2013-12-24 | 2014-04-30 | 兰州空间技术物理研究所 | Simulated analysis method for suspended potential of low earth orbit satellite |
| CN104615841B (en) * | 2015-03-05 | 2017-08-25 | 哈尔滨工业大学 | Consider the spacecraft solar energy sailboard three dimension dynamic simulation method of occlusion effect |
| CN107657079A (en) * | 2017-08-29 | 2018-02-02 | 兰州空间技术物理研究所 | A kind of satellite charged effect induced malfunction and environments match analysis method |
| CN108460188B (en) * | 2018-02-05 | 2021-06-01 | 电子科技大学 | Charge distribution finite element FEM solving algorithm applied to PIC electrostatic model |
| CN109726458A (en) * | 2018-12-18 | 2019-05-07 | 兰州空间技术物理研究所 | A kind of geostationary orbit material inequality electrification emulation mode |
| CN111913083B (en) * | 2020-08-07 | 2023-07-21 | 许昌学院 | Simulation test method for space charge-discharge effect of multilayer film material |
| CN112071357B (en) * | 2020-08-27 | 2022-08-02 | 南京航天航空大学 | SRAM memory charge-discharge effect test system and method based on FPGA |
| CN112149340B (en) * | 2020-09-11 | 2023-02-24 | 兰州空间技术物理研究所 | Satellite surface charging potential calculation method combining PIC and charging current fitting |
-
2007
- 2007-03-26 CN CNB2007100384389A patent/CN100550028C/en not_active Expired - Fee Related
Non-Patent Citations (10)
| Title |
|---|
| PIC方法在空间太阳阵表面充电数值模拟中的应用. 李凯，王立，买胜利，秦晓刚，柳青.第九届全国抗辐射电子学与电磁脉冲学术年会论文汇编，第2007卷. 2007 |
| PIC方法在空间太阳阵表面充电数值模拟中的应用. 李凯，王立，买胜利，秦晓刚，柳青.第九届全国抗辐射电子学与电磁脉冲学术年会论文汇编，第2007卷. 2007 * |
| 一种计算空间电池阵输出的新方法. 张锦绣，梁新刚，王刚.节能技术，第23卷第130(总)期. 2005 |
| 一种计算空间电池阵输出的新方法. 张锦绣，梁新刚，王刚.节能技术，第23卷第130(总)期. 2005 * |
| 地球同步轨道高压太阳电池阵充放电效应研究. 李凯，王立，秦晓刚，柳青，孙彦铮，崔新宇.航天器环境工程，第25卷第2期. 2008 |
| 地球同步轨道高压太阳电池阵充放电效应研究. 李凯，王立，秦晓刚，柳青，孙彦铮，崔新宇.航天器环境工程，第25卷第2期. 2008 * |
| 材料二次电子发射特性对表面充电过程影响的数值模拟研究. 买胜利，王立，李凯，秦晓刚.真空与低温，第12卷第2期. 2006 |
| 材料二次电子发射特性对表面充电过程影响的数值模拟研究. 买胜利，王立，李凯，秦晓刚.真空与低温，第12卷第2期. 2006 * |
| 空间太阳阵二次放电的模拟实验研究. 李凯，谢二庆，林洪峰，王立，孙彦诤，崔新宇.物理学报，第54卷第5期. 2005 |
| 空间太阳阵二次放电的模拟实验研究. 李凯，谢二庆，林洪峰，王立，孙彦诤，崔新宇.物理学报，第54卷第5期. 2005 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101276378A (en) | 2008-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100550028C (en) | The method of satellite sun array surface charging numerical simulation prediction | |
| Fang et al. | On the effect of the Martian crustal magnetic field on atmospheric erosion | |
| Muranaka et al. | Development of multi-utility spacecraft charging analysis tool (MUSCAT) | |
| Paranicas et al. | Europa’s radiation environment and its effects on the surface | |
| Zimmerman et al. | Grain‐scale supercharging and breakdown on airless regoliths | |
| Boovaragavan et al. | Towards real-time (milliseconds) parameter estimation of lithium-ion batteries using reformulated physics-based models | |
| CN103226638B (en) | Estimation method for numerical simulation of distribution character of plasma generated by electric propulsor | |
| Kallio et al. | Kinetic simulations of finite gyroradius effects in the lunar plasma environment on global, meso, and microscales | |
| Aizawa et al. | Cross-comparison of global simulation models applied to Mercury’s dayside magnetosphere | |
| CN103761359A (en) | Simulated analysis method for suspended potential of low earth orbit satellite | |
| Vichare et al. | Quiet-time low latitude ionospheric electrodynamics in the non-hydrostatic Global Ionosphere–Thermosphere Model | |
| Alken et al. | The ionospheric gravity and diamagnetic current systems | |
| CN104239623A (en) | Method for obtaining satellite surface potential based on multi-time-scale particle push | |
| Jordanova et al. | The RAM-SCB model and its applications to advance space weather forecasting | |
| CN104281740A (en) | Method for acquiring satellite surface potential on basis of non-uniform grid division | |
| Gunell et al. | Planetary ENA imaging: Effects of different interaction models for Mars | |
| Forest et al. | An open-source spacecraft plasma interaction simulation code PicUp3D: tests and validations | |
| Ivan et al. | Three-dimensional MHD on cubed-sphere grids: Parallel solution-adaptive simulation framework | |
| Wang et al. | Modeling the Earth’s magnetosphere under the influence of solar wind with due northward IMF by the AMR-CESE-MHD model | |
| Holmstrom | An energy conserving parallel hybrid plasma solver | |
| Liu et al. | A rapid prediction model of photovoltaic power generation for autonomous long‐duration aerostat | |
| Muranaka et al. | Recent progress of development of multi-utility spacecraft charging analysis tool (MUSCAT) | |
| Song | Evaluation of Coronal and Interplanetary Magnetic Field Extrapolation Using PSP Solar Wind Observation | |
| Barminova et al. | Modeling of intense charged particle bunch dynamics in external magnetic fields | |
| Scharlemann et al. | Influence of the Solar Arrays on the Floating Potential of SMART-1: Numerical Simulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091014 Termination date: 20130326 |