Cospar
1 like204 views
The document describes a new code and graphical user interface (GUI) for the European Space Agency's (ESA) meteoroid model. The code includes a FORTRAN calculation engine and a Java-based GUI. It allows users to estimate dust flux, number density, and velocity from meteor showers by specifying a target trajectory and selecting model parameters and output formats through the GUI. Future plans include optimizing the code for parallel processing and adding new modeling capabilities.
Cospar
- 1. The new code for ESA meteoroid model. Alexey Mints1 Valery Dikarev1 Gerhard Drolshagen2 1 University of Bielefeld, Germany 2 European Space Agency, Noordwijk, The Netherlands 18.07.2010
- 2. Problem Outline 1 Problem 2 Code 3 GUI A. Mints, 18.07.2010 The new code for ESA meteoroid model. 2 / 22
- 3. Problem Requirements Input: Target trajectory and geometry, output content specification; Output: Estimated dust flux, number density and average velocity. Dust distributions in mass, incidence direction and velocity. GUI: A tool to set input parameters and view output. A. Mints, 18.07.2010 The new code for ESA meteoroid model. 3 / 22
- 4. Problem Orbital elements A. Mints, 18.07.2010 The new code for ESA meteoroid model. 4 / 22
- 5. Problem Model (as of now) 5 populations, for each: • Rectangular 3D grid in orbital space. Dimensions: pericenter distance (0.05-6 a.u., 50 log-scale bins), eccentricity (0-1, 100 bins) and inclination (0-180, 180 bins); • Mass spectra (200 log-scale bins from 10−18 to 105 grams), independent of orbital elements; Current model file size ∼28 Mbytes. A. Mints, 18.07.2010 The new code for ESA meteoroid model. 5 / 22
- 6. Problem Grids Figure: Possible grids: regular (left), irregular (right) Old IMEM New IMEM Regular (orthogonal) grid in Regular (spherical) grid in orbital elements. incidence velocity. A. Mints, 18.07.2010 The new code for ESA meteoroid model. 6 / 22
- 7. Problem Old and new approaches Old IMEM New IMEM Jacobian diverges, tricks are Symmetries cannot be used; needed; No Jacobian needed; Incidence direction has to be Incidence direction emerges calculated: have to run over naturally from the grid; the whole model range; Reduced calculations for Symmetries can be used; sensitivity; Various scans can be easily implemented; A. Mints, 18.07.2010 The new code for ESA meteoroid model. 7 / 22
- 8. Problem Thresholds C1 < mα V β < C2 • α = 1; β = 0 — mass threshold; • α = 1; β = 1 — momentum threshold; • α = 1; β = 2 — kinetic energy threshold; A. Mints, 18.07.2010 The new code for ESA meteoroid model. 8 / 22
- 9. Problem Calculation grid ∗ grid ∗ Vdust = Vdust + Vtarget Geometrical sensitivity can be calculated from ϕ; grid Knowing Vdust we can calculate mass densities as C2 /V β C1 /V β mα dm; A. Mints, 18.07.2010 The new code for ESA meteoroid model. 9 / 22
- 10. Code Outline 1 Problem 2 Code 3 GUI A. Mints, 18.07.2010 The new code for ESA meteoroid model. 10 / 22
- 11. Code Application layout Java Fortran Interface Model Data module Results Output data datafiles I/O routines Core User Input Task data data datafiles Global settings A. Mints, 18.07.2010 The new code for ESA meteoroid model. 11 / 22
- 12. Code Application composition • Calculation engine — FORTRAN-95 program (∼4000 lines); • GUI — Java graphical interface (∼11000 lines) developed with NetBeans ; • User documentation; • Model file; • Sample task and trajectory files; A. Mints, 18.07.2010 The new code for ESA meteoroid model. 12 / 22
- 13. Code Engine input files • Task file — over 30 parameters, defining task properties and output content and format; • Trajectory file — contains orbital parameters or point-by-point trajectory; • Sensitivity file — optional file containing sensitivity function; • Model file — dust orbital distribution model (binary file); A. Mints, 18.07.2010 The new code for ESA meteoroid model. 13 / 22
- 14. Code Task file D e s c r i p t i o n= P o p u l a t i o n s=∗ M e t e o r o i d m o d e l= T r a j e c t o r y f i l e =x . t r j S e n s i t i v i t y f i l e =t . t s k s e n s P l o t s e t t i n g s=t . t s k p l o t R e s u l t=t . t s k r e s S e n s i t i v i t y p r e s e t =0 P o p u l a t i o n s s t y l e =1 P o i n t s s t y l e =1 C o o r d i n a t e s s t y l e =1 E a r t h c o o r d i n a t e s =0 Scan mode=0 P h y s u n i t s =0 F l u x =1 N u m b e r d e n s i t y=0 A v e r a g e v e l o c i t y =1 T h r e s h o l d s=m, 1 . 0 0 e −18, 1 . 0 0 e+02 I n c i d e n c e r a n g e =0.00:180.000000 V e l o c i t y r a n g e =0.000000:100000.000000 T a r g e t t y p e =0 T a r g e t o r i e n t a t i o n =0 O r b i t k i n d =1 O r b i t s o l a r =0 O r b i t e c l i p t i c a l =0 A. Mints, 18.07.2010 The new code for ESA meteoroid model. 14 / 22
- 15. Code Engine output #IMEM2 o u t p u t f i l e #C r e a t e d =20/ 5/2010 2 0 : 2 7 : 3 5 . 1 6 0 #P h y s u n i t s=F #Scan mode=0 #S c a n r e s o l u t i o n= 20 #P o p u l a t i o n s s t y l e =1 #P o i n t s s t y l e =1 #T h r e s h o l d s i n c o l u m n s=F #C o o r d i n a t e s s t y l e =1 #E a r t h c o o r d i n a t e s=F #N u m b e r d e n s i t y=F #F l u x=T #A v e r a g e v e l o c i t y=T #T a r g e t t y p e =0 #T a r g e t o r i e n t a t i o n =0 #S e n s i t i v i t y p r e s e t =0 #T r a j e c t o r y t y p e =0 #T r a j e c t o r y t y p e= 0.: 180. #M i s s i o n d u r a t i o n= 3.0 #M i s s i o n l a u n c h= 0.0 #P o i n t Thr | Time | asteroids collisions |... # | | Flux | AvgV |... # 1 | 2 | 3 | 4 | 5 |... 1 1 0 . 0 0 0 0 0 0 0 0 E+00 0 . 0 0 0 0 0 0 0 E+000 0 . 0 0 0 0 0 0 0 E+000 . . . 2 1 0 . 1 0 0 0 0 0 0 0 E+01 0 . 0 0 0 0 0 0 0 E+000 0 . 0 0 0 0 0 0 0 E+000 . . . 3 1 0 . 2 0 0 0 0 0 0 0 E+01 0 . 0 0 0 0 0 0 0 E+000 0 . 0 0 0 0 0 0 0 E+000 . . . 4 1 0 . 3 0 0 0 0 0 0 0 E+01 0 . 0 0 0 0 0 0 0 E+000 0 . 0 0 0 0 0 0 0 E+000 . . . . A. .Mints, 18.07.2010 .. The new code for ESA meteoroid model. 15 / 22
- 16. GUI Outline 1 Problem 2 Code 3 GUI A. Mints, 18.07.2010 The new code for ESA meteoroid model. 16 / 22
- 17. GUI Main window A. Mints, 18.07.2010 The new code for ESA meteoroid model. 17 / 22
- 18. GUI Output properties window A. Mints, 18.07.2010 The new code for ESA meteoroid model. 18 / 22
- 19. GUI Progress window A. Mints, 18.07.2010 The new code for ESA meteoroid model. 19 / 22
- 20. GUI Plots A. Mints, 18.07.2010 The new code for ESA meteoroid model. 20 / 22
- 21. GUI Maps A. Mints, 18.07.2010 The new code for ESA meteoroid model. 21 / 22
- 22. GUI Future plans • Final release – September 2010; • OpenMP and MPI extensions; • Web-interface; • More Engine features (for example, meteor flux for a given location on Earth); A. Mints, 18.07.2010 The new code for ESA meteoroid model. 22 / 22