mms_neutrals_circular_m

MMS solutions and sources for braginskii model in circular geometry

Uses

- circular_equilibrium_m
- constants_m
- params_neut_switches_m
- error_handling_grillix_m
- screen_io_m
- rate_coeff_neutrals_m
- precision_grillix_m
- equilibrium_m
- params_neut_model_m
- params_brag_model_m
- mms_braginskii_circular_m
- status_codes_grillix_m

Functions

mms_sol_circular_neutrals_dens mms_sol_circular_neutrals_parmom mms_sol_circular_neutrals_pressure mms_source_circular_neutrals_dens mms_source_circular_neutrals_parmom mms_source_circular_neutrals_pressure

Functions

public function mms_sol_circular_neutrals_dens(equi, x, y, z, t)

MMS solution for neutrals density

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time

Return Value real(kind=gp)

public function mms_sol_circular_neutrals_parmom(equi, x, y, z, t)

MMS solution for neutrals parallel momentum

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time

Return Value real(kind=gp)

public function mms_sol_circular_neutrals_pressure(equi, x, y, z, t)

MMS solution for neutrals temperature

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time

Return Value real(kind=gp)

public function mms_source_circular_neutrals_dens(equi, x, y, z, t, chi, epsinv, neutrals_dens_pen, k_iz_o, k_rec_o)

MMS source for neutral density equation

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time
real(kind=GP),	intent(in)	::	chi	Value of characteristc function of penalisation
real(kind=GP),	intent(in)	::	epsinv	Inverse of penalisation epsilon
real(kind=GP),	intent(in)	::	neutrals_dens_pen	Target value to be penalised towards
type(rate_coeff_neutrals_t),	intent(in)	::	k_iz_o	Ionization rate coefficient object
type(rate_coeff_neutrals_t),	intent(in)	::	k_rec_o	Recombination rate coefficient object

Return Value real(kind=gp)

public function mms_source_circular_neutrals_parmom(equi, x, y, z, t, chi, epsinv, neutrals_parmom_pen, k_iz_o, k_rec_o)

MMS source for neutral parallel momentum equation

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time
real(kind=GP),	intent(in)	::	chi	Value of characteristc function of penalisation
real(kind=GP),	intent(in)	::	epsinv	Inverse of penalisation epsilon
real(kind=GP),	intent(in)	::	neutrals_parmom_pen	Target value to be penalised towards
type(rate_coeff_neutrals_t),	intent(in)	::	k_iz_o	Ionization rate coefficient object
type(rate_coeff_neutrals_t),	intent(in)	::	k_rec_o	Recombination rate coefficient object

Return Value real(kind=gp)

public function mms_source_circular_neutrals_pressure(equi, x, y, z, t, chi, epsinv, neutrals_pressure_pen, k_iz_o, k_rec_o)

MMS source for neutral temperature equation

Arguments

Type	Intent		Name
class(equilibrium_t),	intent(in)	::	equi	Equilibrium
real(kind=GP),	intent(in)	::	x	x-coordinate
real(kind=GP),	intent(in)	::	y	y-coordinate
real(kind=GP),	intent(in)	::	z	Axial coordinate z (represented by toroidal coordinate)
real(kind=GP),	intent(in)	::	t	Time
real(kind=GP),	intent(in)	::	chi	Value of characteristc function of penalisation
real(kind=GP),	intent(in)	::	epsinv	Inverse of penalisation epsilon
real(kind=GP),	intent(in)	::	neutrals_pressure_pen	Target value to be penalised towards
type(rate_coeff_neutrals_t),	intent(in)	::	k_iz_o	Ionization rate coefficient object
type(rate_coeff_neutrals_t),	intent(in)	::	k_rec_o	Recombination rate coefficient object

mms_neutrals_circular_m Module

Contents

Functions

Uses

Contents

Functions

Functions

public function mms_sol_circular_neutrals_dens(equi, x, y, z, t)

Arguments

Return Value real(kind=gp)

public function mms_sol_circular_neutrals_parmom(equi, x, y, z, t)

Arguments

Return Value real(kind=gp)

public function mms_sol_circular_neutrals_pressure(equi, x, y, z, t)

Arguments

Return Value real(kind=gp)

public function mms_source_circular_neutrals_dens(equi, x, y, z, t, chi, epsinv, neutrals_dens_pen, k_iz_o, k_rec_o)

Arguments

Return Value real(kind=gp)

public function mms_source_circular_neutrals_parmom(equi, x, y, z, t, chi, epsinv, neutrals_parmom_pen, k_iz_o, k_rec_o)

Arguments

Return Value real(kind=gp)

public function mms_source_circular_neutrals_pressure(equi, x, y, z, t, chi, epsinv, neutrals_pressure_pen, k_iz_o, k_rec_o)

Arguments

Return Value real(kind=gp)