100 likes | 217 Views
Kinetic Modeling of the Sheath Scale in the Lunar Plasma Environment. Peter Messmer*, Keegan Amyx, Peter Stoltz, Andrew Poppe, Mihay Horanyi, Scott Robertson, Zoltan Sternovsky messmer@txcorp.com. Tech-X Corporation 5621 Arapahoe Ave., Boulder, CO 80303 http://www.txcorp.com.
E N D
Kinetic Modeling of the Sheath Scale in the Lunar Plasma Environment Peter Messmer*, Keegan Amyx, Peter Stoltz, Andrew Poppe, Mihay Horanyi, Scott Robertson, Zoltan Sternovsky messmer@txcorp.com Tech-X Corporation 5621 Arapahoe Ave., Boulder, CO 80303 http://www.txcorp.com CCLDAS All Hands Meeting, Boulder, CO, July 10, 2009
VORPAL -A Plasma ModelingFramework • Original target applications: • Laser Wakefield Acceleration • PIC, Fluid, Hybird • Electrostatic, EM • Multi-Dimensional (N=1,2,3) • Fully parallel • Scaling for > 32,000 PEs • Flexible domain decomposition • Broad range of physics features: • Complex geometries • Ionization, recombination, CEX physics • Field ionization http://www.txcorp.com/products/VORPAL
Code/setup Validation with 1D Photoelectron Sheath 2D ES simulation, Y periodic , 200 x 10 cells l0/dx = 10 50 particles per cell nominal Monoenergetic Garad&Tunaley Simulation Maxwellian • R. Garad & J Tunaley, JGR 76(10), • 2498, 1971 • A. Poppe & M. Horanyi, WPDP, 2009
2D Monoenergetic Sheath electrons • l0 protons • l0 • 2D ES simulation, left wall = 0V • 200 x 100 cells • Electrons, Protons • Monoenergetic ,V0 = 200 eV (!)
Scenario with Surface Charging Photo- electrons No surface charging Surface-Charging (just for comparison)
2D Thermal Sheath with Surface Charging Electron impact creates “heavy electrons” Electrons get absorbed • 2D ES simulation, left wall = 0V • 200 x 100 cells • Electrons • Heavy Protons, Heavy electrons (m/m0= 5000) • Vsig = 3eV, Vtherm = 3 eV
Electric field mainly due to positive charge of emitting region Charging Non Charging
“Heavy Electrons” follow the electric field lines • 2D ES simulation, left wall = 0V • 200 x 100 cells • Electrons • Heavy Protons, Heavy electrons (m/m0= 5000) • Vsig = 3eV, Vtherm = 3 eV
Initial 3D simulations Charging of surface No charging of surface • 3D ES simulation, bottom wall = 0V • 10 x 100 x 100 cells • Electrons • Protons, Heavy electrons (m/m0= 5000) • Vsig = 3eV, Vtherm = 3 eV
Summary / Conclusions / Future work • Presented VORPAL simulations of plasma sheath • Validated with kinetic theory for 1D sheath • Presented 2D simulation with/without surface charging • “heavy electrons” move in electrostatic field, follow (curved) field lines • Future work: • Convergence studies, more realistic parameters • Inclusion of solar wind • Time dependent problems, angular dependency of photo-currents • Complex geometries (crater, habitat, instrument) • 3D Work supported by CCLDAS and Tech-X Corp.