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.

Contents

Variables

a

Functions

rho0 T0 dT0 v03 dv03 B02 dB02 B03

Module Procedures

kh_cd_instability_eq

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	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dT0(r)

Arguments

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

Return Value real(kind=dp)

function v03(r)

Arguments

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

Return Value real(kind=dp)

function dv03(r)

Arguments

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

Return Value real(kind=dp)

function B02(r)

Arguments

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

Return Value real(kind=dp)

function dB02(r)

Arguments

Type	Intent	Optional	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 kh_cd_instability_eq module module subroutine kh_cd_instability_eq(settings, grid, background, physics) Interface →

Sets the equilibrium.

Arguments

Type	Intent	Optional	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