smod_equil_RTI Submodule

This submodule defines Rayleigh-Taylor 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 a) 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 = 15 : gravitational constant.
  • alpha = 0 : magnetic shear.
  • theta = 0.35 : angle used in the velocity profile.
  • p1 = 0.2 : , used in the velocity profile.
  • p2 = 0.6 : , used in the velocity profile.
  • p3 = 0 : , used in the velocity profile.
  • p4 = -0.35 : , used in the magnetic field profile.
  • tau = 0 : 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.html RTI_eq module subroutine RTI_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