examples.dg.mhd.gem_challenge_rmhd Namespace Reference

Variables
str	torder = 'RK4'
str	sorder = 'secondOrder'
bool	useLimiter = False
bool	always_limit = False
float	penalty_beta = -0.5
float	penalty_eta = 1.0
str	gradient_method = "ldg"
float	skin_depth_norm = 1.0
float	nu_p_norm = 0.1
int	gamma = 5./3.
int	z_i = 1
int	A_i = 1
int	theta = 1
int	di_L = 1.
int	Lx = 8*math.pi*di_L
int	Ly = 4*math.pi*di_L
float	sheet_width = 0.5
float	n0 = 1.0
float	n_inf = 0.2
float	B0 = 1.0
float	Ti_Te = 1.0/theta
float	Te = B0**2/2.0/n0/(Ti_Te + 1.0)
float	Ti = Ti_Te * Te
list	T = [Ti, Te]
str	imhd_numerical_flux_type = "HLLD"
warpy	mesh
list	component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']
warpy	fluid = warpy.variable(name='fluid',components=component_names,basis_array_set=sorder)
list	gradient_component_names = [prefix+"_"+suffix for prefix in component_names for suffix in ['x','y','z']]
warpy	fluid_gradient = warpy.variable(name='fluid_gradient', components=gradient_component_names,basis_array_set=sorder)
warpy	fluid_gradient2 = warpy.variable(name='fluid_gradient2', components=gradient_component_names,basis_array_set=sorder)
warpy	gem_ic
list	applicators = [warpy.applicator(spatial_order=sorder,fun=gem_ic,var=fluid,spatial_scheme='Nodal')]
list	apps = []
list	boundary_condition_q_apps = []
list	boundary_condition_grad_q_apps = []
warpy	writer
warpy	spatial_solver
list	variable_adjusters = []
warpy	ti
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=5e-5)
warpy	sim
datetime	now = datetime.datetime.now()
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

A_i

int examples.dg.mhd.gem_challenge_rmhd.A_i = 1

always_limit

bool examples.dg.mhd.gem_challenge_rmhd.always_limit = False

applicators

list examples.dg.mhd.gem_challenge_rmhd.applicators = [warpy.applicator(spatial_order=sorder,fun=gem_ic,var=fluid,spatial_scheme='Nodal')]

apps

list examples.dg.mhd.gem_challenge_rmhd.apps = []

B0

float examples.dg.mhd.gem_challenge_rmhd.B0 = 1.0

boundary_condition_grad_q_apps

list examples.dg.mhd.gem_challenge_rmhd.boundary_condition_grad_q_apps = []

boundary_condition_q_apps

list examples.dg.mhd.gem_challenge_rmhd.boundary_condition_q_apps = []

component_names

list examples.dg.mhd.gem_challenge_rmhd.component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']

detect_nonscalar

examples.dg.mhd.gem_challenge_rmhd.detect_nonscalar

di_L

int examples.dg.mhd.gem_challenge_rmhd.di_L = 1.

dt_controller

warpy examples.dg.mhd.gem_challenge_rmhd.dt_controller = warpy.dt_calc.fixed_dt(init_dt=5e-5)

fluid

warpy examples.dg.mhd.gem_challenge_rmhd.fluid = warpy.variable(name='fluid',components=component_names,basis_array_set=sorder)

fluid_gradient

warpy examples.dg.mhd.gem_challenge_rmhd.fluid_gradient = warpy.variable(name='fluid_gradient', components=gradient_component_names,basis_array_set=sorder)

fluid_gradient2

warpy examples.dg.mhd.gem_challenge_rmhd.fluid_gradient2 = warpy.variable(name='fluid_gradient2', components=gradient_component_names,basis_array_set=sorder)

gamma

int examples.dg.mhd.gem_challenge_rmhd.gamma = 5./3.

gem_ic

warpy examples.dg.mhd.gem_challenge_rmhd.gem_ic

Initial value:

=  warpy.functions.gem_hallmhd(name='gem_ic',
                                     gamma=gamma,
                                     mass=A_i,
                                     temperature=Ti+Te,
                                     # velocity=V,
                                     sheet_width=sheet_width,
                                     number_density=n0,
                                     number_density_infinity=n_inf,
                                     B0=B0,
                                     Lx=Lx,
                                     Ly=Ly)

gen_xdmf

examples.dg.mhd.gem_challenge_rmhd.gen_xdmf

gradient_component_names

list examples.dg.mhd.gem_challenge_rmhd.gradient_component_names = [prefix+"_"+suffix for prefix in component_names for suffix in ['x','y','z']]

gradient_method

str examples.dg.mhd.gem_challenge_rmhd.gradient_method = "ldg"

imhd_numerical_flux_type

str examples.dg.mhd.gem_challenge_rmhd.imhd_numerical_flux_type = "HLLD"

Lx

int examples.dg.mhd.gem_challenge_rmhd.Lx = 8*math.pi*di_L

Ly

int examples.dg.mhd.gem_challenge_rmhd.Ly = 4*math.pi*di_L

mesh

warpy examples.dg.mhd.gem_challenge_rmhd.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-Lx/2.,Lx/2., -Ly/2.,Ly/2.],
                        # NumCells=[64,32],
                        NumCells=[32,16],
                        NodeSets=['Left','Right', 'Lower','Upper'],
                        NumLayers=1,
                        basis_array_set=sorder,
                        PeriodicBoundaries=['Left','Right'])

n0

float examples.dg.mhd.gem_challenge_rmhd.n0 = 1.0

n_inf

float examples.dg.mhd.gem_challenge_rmhd.n_inf = 0.2

now

datetime examples.dg.mhd.gem_challenge_rmhd.now = datetime.datetime.now()

nu_p_norm

float examples.dg.mhd.gem_challenge_rmhd.nu_p_norm = 0.1

penalty_beta

float examples.dg.mhd.gem_challenge_rmhd.penalty_beta = -0.5

penalty_eta

float examples.dg.mhd.gem_challenge_rmhd.penalty_eta = 1.0

sheet_width

float examples.dg.mhd.gem_challenge_rmhd.sheet_width = 0.5

sim

warpy examples.dg.mhd.gem_challenge_rmhd.sim

Initial value:

=  warpy.dg_sim(name='test_0p1_nu_p_tau_negate',
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[ti],
                   writers=[writer],
                   # time=[0,3.25],
                   time=[0,0.04],                   
                   dt_controller=dt_controller,
                   write_steps=100,
                   verbosity='debug')

skin_depth_norm

float examples.dg.mhd.gem_challenge_rmhd.skin_depth_norm = 1.0

sorder

str examples.dg.mhd.gem_challenge_rmhd.sorder = 'secondOrder'

spatial_solver

warpy examples.dg.mhd.gem_challenge_rmhd.spatial_solver

Initial value:

=  warpy.spatial_solvers.dg(name="dg",
                                          spatial_order=sorder,
                                          penalty_beta=penalty_beta,
                                          penalty_eta=penalty_eta,
                                          applications=apps,
                                          on_subdomains=['all'])

T

list examples.dg.mhd.gem_challenge_rmhd.T = [Ti, Te]

Te

float examples.dg.mhd.gem_challenge_rmhd.Te = B0**2/2.0/n0/(Ti_Te + 1.0)

theta

int examples.dg.mhd.gem_challenge_rmhd.theta = 1

Ti

float examples.dg.mhd.gem_challenge_rmhd.Ti = Ti_Te * Te

ti

warpy examples.dg.mhd.gem_challenge_rmhd.ti

Initial value:

=  warpy.host_actions.erk(name="rk",
                            scheme=torder,
                            spatial_solvers=[spatial_solver],
                            variable_adjusters=variable_adjusters)

Ti_Te

float examples.dg.mhd.gem_challenge_rmhd.Ti_Te = 1.0/theta

torder

str examples.dg.mhd.gem_challenge_rmhd.torder = 'RK4'

True

examples.dg.mhd.gem_challenge_rmhd.True

useLimiter

bool examples.dg.mhd.gem_challenge_rmhd.useLimiter = False

variable_adjusters

list examples.dg.mhd.gem_challenge_rmhd.variable_adjusters = []

writer

warpy examples.dg.mhd.gem_challenge_rmhd.writer

Initial value:

=  warpy.host_actions.writer(name='writer',
                                   ReadVars=[fluid,
                                             fluid_gradient,
                                             fluid_gradient2])

z_i

int examples.dg.mhd.gem_challenge_rmhd.z_i = 1