smod_equil_tc_pinch Submodule

This submodule defines a steady Taylor-Couette flow in a cylindrical geometry where a plasma is confined between two (rotating) coaxial cylinders with an azimuthal magnetic field and constant resistivity.

This equilibrium is taken from Shalybkov, Dima. "Rotational stabilization of pinch instabilities in Taylor-Couette flow.", Physical Review E 75, 047302 (2007).



Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=dp), private :: h
real(kind=dp), private :: A
real(kind=dp), private :: B
real(kind=dp), private :: A2
real(kind=dp), private :: B2
real(kind=dp), private :: Bc1
real(kind=dp), private :: Bc2
real(kind=dp), private :: Bc3
real(kind=dp), private :: Tstart
real(kind=dp), private :: Tend
real(kind=dp), private :: Ta
real(kind=dp), private :: Ha
real(kind=dp), private :: Re
real(kind=dp), private :: Pm

Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function v02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dv02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function ddv02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function T0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dT0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function B02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dB02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function ddB02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)


Module Procedures

module procedure tc_pinch_eq module module subroutine tc_pinch_eq(settings, grid, background, physics)

Arguments

TypeIntentOptionalAttributesName
type(settings_t), intent(inout) :: settings
type(grid_t), intent(inout) :: grid
type(background_t), intent(inout) :: background
type(physics_t), intent(inout) :: physics