GEM WAVE MODELLING CHALLENGE PRELIMINARY NOTES

1. GEM WAVE MODELLING CHALLENGEPRELIMINARY NOTES AGU AUTUMN MEETING SAN FRANCISCO 9 DECEMBER, 2013 David Nunn1 1. School of Electronics and Computer Science Southampton University, U.K.

2. GEM WAVE CHALLENGE • PROVISIONAL SIMULATION RESULTS • Data: L=5; Ne=5.4/cc; Tperp= 62keV; T||=20keV • Maximum linear growth rate at z=0 1146db/s • The code: 1D k||Bo ; VlasovHybrid Simulation • Narrow band code has bandwidth ~100Hz (variable). Centre of code • internal bandpass filter matched to cg of actual spectrum. • It was shown in Katoh et al,2006 that frequency sweep due entirely • to a spatial gradient of frequency in the generation region . • i.e. an advective term. Omura and Nunn(2012),Nunn and Omura,2013 • showed a mechanism for setting up this gradient thru the Jb(z,t) field • Narrow band code can allow linear spatial frequency variation • by using matched filters • Broadband version of code has bandwidth ~1000Hz and can • accommodate spatial frequency variation without matched filtering. • Vlasov code being low noise requires trigger signal. This may be a CW • pulse (narrow band code), or broadband noise across simulation • bandwidth (broadband code)

3. GEM WAVE CHALLENGE • General results • Growth rates are high, emissions readily triggered. Usually risers with • sweep rates ~3-8kHz/s, depending on triggering pulse, simulation • bandwidth, simulation precision etc. • Many runs self saturate in the region of 200pT, but some cases e.g. • broadband code with hiss excitation require artificial saturation to prevent • wave amplitudes reaching unreasonable values. • The narrow band code may show a tendency to trigger fallers or downward • hooks when the trigger pulse is strong and lengthy. • The broadband code with CW wave excitation produces steep risers with a • fairly broad spectrum ~200 Hz. • Excitation by hiss in the broadband code gives a very swishy riser of • bandwidth ~500Hz

4. f/t spectrogram of chorus simulation with broadband vlasovvhs code; steep riser ~ 8kHz/sCW pulse excitation

5. Broadband vlasov simulation-plot of wave amplitude pT in z,t plane for the entire simulation; note spectral structure andgeneration region extends upstream from equator (Helliwell,67)

6. History in z.t plane of current –Je in phase with e field

7. History in z.t plane of current –Jbin phase with b field

8. Exit amplitude in pT at z=+6000kms from equator

9. Quasi-Broadband vlfvlasov code – hiss excitationf/t spectrogram

10. NARROW BAND VLASOV CODE; SPECTROGRAM OF RISING CHORUS ELEMENT AT Z=+6000Km FROM EQUATOR

11. PLOT OF WAVE AMPLITUDE (pT) IN Z,T PLANE FOR SINGLE CHORUSELEMENT EVENT

12. Plot of resonant particle current component parallel to wave E field (Je) in z,t plane for chorus element event. This component gives non linear growth

13. Plot of resonant particle current component parallel to wave B field (Jb) in z,t plane for chorus element event. This component gives phase and frequency shifts in the wavefield

14. Spatial dft of simulation box wavefield of slow fallerproduced by strong key down excitation pulse 5pT; narrow band vlasov code

15. GEM WAVE CHALLENGE • Ongoing research……. • Parallelise code with MPI relative to Vperp • Long runs with higher precision, particularly in Vperp • Investigate triggering at lower frequencies • Investigate saturation problems • Acknowledgements • These computations done on Kyoto University KUDPC Cray and • Southampton University IRIDIS3/4 systems.

17. Spatial dft of simulation box wavefield of slow fallerproduced by strong key down excitation pulse 5pT; narrow band vlasov code