smod_equil_kelvin_helmholtz_cd Submodule

This submodule defines an unperturbed magnetised jet model in cylindrical geometry, giving rise to Kelvin-Helmholtz and current-driven instabilities. The geometry is fixed for this problem; the cylinder wall is dependent on the equilibrium parameters and is given by 2rj.

This equilibrium is taken from Baty, H., & Keppens, R. (2002). Interplay between Kelvin-Helmholtz and current-driven instabilities in jets. The Astrophysical Journal, 580(2), 800.

Note

Default values are given by

  • k2 = -1
  • V = 1.63 : amplitude of the velocity shear
  • cte_p0 = 1 : used to set the pressure.
  • cte_rho0 = 1 : used to set the density.
  • Bz0 = 0.25 : used to set Bz.
  • rc = 0.5 : length for radial variation.
  • rj = 1 : jet radius

and can all be changed in the parfile.

Note

The default setup is HEL2 in the original paper. For HEL1 you can set rc = 2.


Uses


Variables

Type Visibility Attributes Name Initial
real(kind=dp), private :: a

Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function T0(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dT0(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function v03(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dv03(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function B02(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dB02(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function B03()

Arguments

None

Return Value real(kind=dp)


Module Procedures

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

Sets the equilibrium.

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