smod_equil_KHI Submodule

This submodule defines Kelvin-Helmholtz instabilities in Cartesian geometry. This equilibrium is a specific case of the flow driven instabilities. No magnetic fields are considered in this case (pure HD).

This equilibrium is taken from section 13.2, p. 487 (case b) in Goedbloed, H., Keppens, R., & Poedts, S. (2019). Magnetohydrodynamics of Laboratory and Astrophysical Plasmas. Cambridge University Press. DOI.

Note

Default values are given by

  • k2 = 0
  • k3 = 1
  • cte_rho0 = 1 : used in the density profile.
  • cte_p0 = 10 : used in the pressure profile.
  • delta = 0 : used in the density profile.
  • g = 0 : gravitational constant.
  • alpha = 0 : magnetic shear.
  • theta = 0 : angle used in the velocity profile.
  • p1 = 0 : , used in the velocity profile.
  • p2 = 0 : , used in the velocity profile.
  • p3 = 1 : , used in the velocity profile.
  • p4 = 0 : , used in the magnetic field profile.
  • tau = 11 : used in the velocity profile.

and can all be changed in the parfile.


Uses


Module Procedures

module procedure /home/runner/work/legolas/legolas/docs/ford/module/smod_equil_khi.html KHI_eq module subroutine KHI_eq(settings, grid, background, physics)

Arguments

Type IntentOptional Attributes Name
type(settings_t), intent(inout) :: settings
type(grid_t), intent(inout) :: grid
type(background_t), intent(inout) :: background
type(physics_t), intent(inout) :: physics