smod_equil_RTI_KHI Submodule

This submodule defines Rayleigh-Taylor and Kelvin-Helmholtz instabilities in Cartesian geometry. This equilibrium is a specific case of the flow driven instabilities.

This equilibrium is taken from section 13.2, p. 487 (case c) 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 = 1000 : used in the pressure profile, must be large so T0 > 0.
  • delta = -5 : used in the density profile.
  • g = 100 : gravitational constant.
  • alpha = - : magnetic shear.
  • theta = 0 : angle used in the velocity profile.
  • p1 = 1 : , used in the velocity profile.
  • p2 = 2 : , used in the velocity profile.
  • p3 = 1 : , used in the velocity profile.
  • p4 = 0.5 : , used in the magnetic field profile.
  • tau = 4 : 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_rti_khi.html RTI_KHI_eq module subroutine RTI_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