Non-neutral Plasma Shock

1. Non-neutral Plasma Shock 工 HU Xiwei（胡希伟） HE Yong （何勇） Hu Yemin（胡业民） Huazhong University of Science and Techonology 2006.10.25 Zhejiang University, Hangzhou, China

2. Outline • Introduction • Motivation • One dimensional case • Two dimensional case (with magnetic field produced by shock current) • The electrostatic instabilities on the non-neutral shock front

3. Introductionto non-neutral plasma shock

4. Plasma shock • Plasma shock can arise: when a fluid velocity is larger than the ion sound velocity (without magnetic field) or the Alfven velocity (with magnetic field) . • Asteady shock (or shock front, shock profile) is the result of a balance between compressive and dissipative effects.

5. Front, Upstream and Downstream • Shock front (or shock profile): a steady wave-front, propagating at supersonic speed through the undisturbed fluid. • Upstream⑴: the unshocked (undisturbed) fluid before the shock front • Downstream⑵: the shocked fluid (disturbed) behind the front front

6. The neutral plasma shock • In most cases, we describe the plasma shock with the single (neutral) fluid or magnetohydrodynamic (MHD) equations. So the plasma shock is neutral essentially. • In neutral plasma shock there may be plasma current and magnetic field, but there is notcharge separation and electrostatic field.

7. Electric field in the front • The ion fluid elements will run ahead of the electron fluid elements due to its larger inertial. • This causes a charge separation, then a electrostatic field E and a potential differenceΔφ between the up- and down stream.

8. Effects of electric field in front (I) • The electric field E will hold back the ion fluid elements and draw the electron fluid elements. This effect is similar the role of dissipation in the formation of steady shock (shock front or profile). So, the E could play very important role in the formation of plasma shock profile, especially the collisionless (non-dissipative) shock profile.

9. Effects of electric field in front (II) • The potential difference Δφwill modify the jump (discontinuity) relations – the so called Hugoniot relationsand the critical Mach number Mc (the least velocity for shock emergency) . • The electric field will drive a new kind of electrostatic instability in the shock front as the role of density gradient in the Rayleigh-Taylor instability.

10. Motivation • In our simulation for the imploding shock in a single sonoluminescing bubble, we found that there is a extreme strong electrostatic field in the shock front. • The electric field is caused by the charge separation. • The charge separation is due to the difference of electron fluid and various ion fluid velocities.

11. Profiles of the charge number density Δn, electric fieldE, and electric current densityJstriding on the shock front at the moment near the emergence of the maximum electric field.

13. Electric field in strong laser and in imploding shock

14. One dimensional case • Plasma shock carries low current. • The effects of magnetic field produced by the shock current can be neglected.

15. Static coupled equations in a reference frame on the front The double fluid equations and Poisson equation where is the equivalent heat flow

16. Profiles: ━ M1=1.6　 ┉ M1=1.5

17. Jump conditions across the shock Particle fluxes Momentum fluxes Energy fluxes Where

18. Critical upstream Mach number • When , from we can obtain the minimum upstream velocity for the non-neutral shock emergency

19. Two dimensional case • Cylindrical Non-Neutral Plasma Shock with Current, Electric and magnetic Field. • The current and electric field is in axial (z) direction. • The magnetic field produced by the shock current is in poloidal (θ) direction.

20. Basic steady equations for shock front

21. Jump conditionsof particle, momentumin randzdirection, andenergyfluxes

22. Profiles of velocity and temperature

23. Downstream profiles━ 2D case,　 ┉ 1D case

24. Shock current changes the

25. Summary (I) • The jump conditions across the shock will be changed that the particle flux and energy flux are no longer conservative. • The profiles of plasma parameters in downstream are no longer as same as the profiles in upstream.

26. Summary (II) • The critical Mach number ( the minimum upstream velocity for the non-neutral shock emergency) is larger than its value (1.0) in neutral shock. When J0=0, • The positive shock current (from up- to down-stream) will reduce the , and the negative current will increase the .

27. Electrostatic instabilitiesin shock front(1D case)

28. The linear instability analysis x -- the direction of the shock propagating, y – the direction perpendicular to shock propagating.

29. A ExampleDispersion relation aboutkxDisturbance is in the direction of x– the shock propagating direction

30. Profiles of equilibrium parameters labeled x = 0.1, 0.3, 0.5, 0.7, 0.9 as No.1, No.3, No.5, No.7, No.9

31. The numerical results • The high frequency ( ) approximation vs. low frequency ( ) approximation. • With electric field (E≠ 0 ) vs. E = 0. • With dissipation – viscous and friction -- (α≠0) vs. α=0. • Adiabatic process approximation vs. diabatic (energy equation) case.

32. The high frequency approximation

36. Adiabatic vs. Diabatic (E≠ 0)

37. Summary • The Ecan drive the electrostatic instability in the shock front in both parallel (kx) and perpendicular (ky) directions. • The E and density gradient are the destabilizing factors and dissipation (the viscosity etc.) is the stabilizing factor. • The increasing rate • There is evidence of zero frequency instable mode , which is a kind of absolute instability.

38. Prospect • The electrostatic instability in the shock propagating direction is comparable with the Rayleigh-Taylor instability in neutral plasma shock. • The electrostatic instability with zero frequency in the perpendicular shock direction is comparable with the Richtmyer-Meshkov instability in neutral plasma shock. • These works are in process.

39. THANKS