examples.dg.mhd.pn_test Namespace Reference

Variables
str	torder = 'RK1'
str	sorder = 'firstOrder'
float	cfl = 0.5
int	gamma = 5/3
str	numerical_flux_type = "Rusanov"
float	qe = 1.602177e-19
float	mp = 1.672622e-27
float	mu0 = 1.256637e-6
float	eps0 = 8.854187e-12
int	c_light = 1/(mu0*eps0)**.5
int	L = 1
float	ompt = 4.251e3
float	omct = 13.89
int	skin_depth_norm = 1/omct
float	n0 = mp*omct**2/(mu0*L**2*qe**2)
float	tau = ompt*L/(omct*c_light)
float	b0 = mp*c_light*omct**2/(qe*L*ompt)
float	e0 = b0**2/mu0
float	T0 = e0/(n0*qe)
float	nu_p_norm = 1.0
int	atom = 1
int	Zi = 1
int	Ze = -1
int	Zn = 0
int	Ai = 1
int	Ae = 1
int	An = 1 + Ae
float	phi_ion_eV = 13.6
float	phi_ion = phi_ion_eV/T0
warpy	mesh
list	component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']
warpy	plasma
warpy	ic_fun = warpy.functions.custom_mhd(name='ic_fun',skin_depth_norm=skin_depth_norm,gamma=gamma)
list	ic_applicators
list	apps = []
list	bc_apps = []
float	rhoin = 0.1
float	pxin = 0.0
float	pyin = 0.0
float	pzin = 0.0
float	pressurein = 0.1
float	energyin = pressurein/(gamma-1)
warpy	spatial_solver
list	var_adj = []
warpy	writer = warpy.host_actions.writer(name='writer',ReadVars=[plasma])
warpy	time_integrator
int	dt = 1e-4
int	nout = 30
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

Ae

int examples.dg.mhd.pn_test.Ae = 1

Ai

int examples.dg.mhd.pn_test.Ai = 1

An

int examples.dg.mhd.pn_test.An = 1 + Ae

apps

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

atom

int examples.dg.mhd.pn_test.atom = 1

b0

float examples.dg.mhd.pn_test.b0 = mp*c_light*omct**2/(qe*L*ompt)

bc_apps

list examples.dg.mhd.pn_test.bc_apps = []

c_light

int examples.dg.mhd.pn_test.c_light = 1/(mu0*eps0)**.5

cfl

float examples.dg.mhd.pn_test.cfl = 0.5

component_names

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

detect_nonscalar

examples.dg.mhd.pn_test.detect_nonscalar

dt

int examples.dg.mhd.pn_test.dt = 1e-4

dt_controller

warpy examples.dg.mhd.pn_test.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

e0

float examples.dg.mhd.pn_test.e0 = b0**2/mu0

energyin

float examples.dg.mhd.pn_test.energyin = pressurein/(gamma-1)

eps0

float examples.dg.mhd.pn_test.eps0 = 8.854187e-12

gamma

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

gen_xdmf

examples.dg.mhd.pn_test.gen_xdmf

ic_applicators

list examples.dg.mhd.pn_test.ic_applicators

Initial value:

=  [warpy.applicator(spatial_order=sorder,
                                   fun=ic_fun,
                                   var=plasma,components=None,
                                   spatial_scheme='Nodal')]

ic_fun

warpy examples.dg.mhd.pn_test.ic_fun = warpy.functions.custom_mhd(name='ic_fun',skin_depth_norm=skin_depth_norm,gamma=gamma)

L

int examples.dg.mhd.pn_test.L = 1

mesh

warpy examples.dg.mhd.pn_test.mesh

Initial value:

=  warpy.mesh.block(Bounds=[0.0, 1.0, -0.5, 0.5],
                        # NumCells=[300,60],
                        # NumCells=[150,30],
                        # NumCells=[30,6],
                        NumCells=[50,50],
                        # NumCells=[2000,400],
                        NodeSets=['leftWall','rightWall',
                                  'bottomWall','topWall'],
                        NumLayers=1,
                        PeriodicBoundaries = ['bottomWall', 'topWall'],
                        # PeriodicBoundaries = ['leftWall', 'rightWall'],
                        # PeriodicBoundaries = ['leftWall','rightWall',
                        #                       'bottomWall', 'topWall'],
                        basis_array_set=sorder,
                        block_type="symmetric")

mp

float examples.dg.mhd.pn_test.mp = 1.672622e-27

mu0

float examples.dg.mhd.pn_test.mu0 = 1.256637e-6

n0

float examples.dg.mhd.pn_test.n0 = mp*omct**2/(mu0*L**2*qe**2)

nout

int examples.dg.mhd.pn_test.nout = 30

nu_p_norm

float examples.dg.mhd.pn_test.nu_p_norm = 1.0

numerical_flux_type

str examples.dg.mhd.pn_test.numerical_flux_type = "Rusanov"

omct

float examples.dg.mhd.pn_test.omct = 13.89

ompt

float examples.dg.mhd.pn_test.ompt = 4.251e3

phi_ion

float examples.dg.mhd.pn_test.phi_ion = phi_ion_eV/T0

phi_ion_eV

float examples.dg.mhd.pn_test.phi_ion_eV = 13.6

plasma

warpy examples.dg.mhd.pn_test.plasma

Initial value:

=  warpy.variable(
    name='plasma',components=component_names,basis_array_set=sorder)

pressurein

float examples.dg.mhd.pn_test.pressurein = 0.1

pxin

float examples.dg.mhd.pn_test.pxin = 0.0

pyin

float examples.dg.mhd.pn_test.pyin = 0.0

pzin

float examples.dg.mhd.pn_test.pzin = 0.0

qe

float examples.dg.mhd.pn_test.qe = 1.602177e-19

rhoin

float examples.dg.mhd.pn_test.rhoin = 0.1

sim

warpy examples.dg.mhd.pn_test.sim

Initial value:

=  warpy.dg_sim(name='pnt00e',
                   meshes=[mesh],
                   initial_conditions=ic_applicators,
                   temporal_solvers=[time_integrator],
                   writers=[writer],
                   time=[0.,1.0],
                   dt_controller=dt_controller,
                   write_steps=nout,
                   verbosity='debug',
                   enable_fpe=True)

skin_depth_norm

int examples.dg.mhd.pn_test.skin_depth_norm = 1/omct

sorder

str examples.dg.mhd.pn_test.sorder = 'firstOrder'

spatial_solver

warpy examples.dg.mhd.pn_test.spatial_solver

Initial value:

=  warpy.spatial_solvers.dg(name="dg",
                                          spatial_order=sorder,
                                          applications=apps,
                                          cfl=cfl)

T0

float examples.dg.mhd.pn_test.T0 = e0/(n0*qe)

tau

float examples.dg.mhd.pn_test.tau = ompt*L/(omct*c_light)

time_integrator

warpy examples.dg.mhd.pn_test.time_integrator

Initial value:

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

torder

str examples.dg.mhd.pn_test.torder = 'RK1'

True

examples.dg.mhd.pn_test.True

var_adj

list examples.dg.mhd.pn_test.var_adj = []

writer

warpy examples.dg.mhd.pn_test.writer = warpy.host_actions.writer(name='writer',ReadVars=[plasma])

Ze

int examples.dg.mhd.pn_test.Ze = -1

Zi

int examples.dg.mhd.pn_test.Zi = 1

Zn

int examples.dg.mhd.pn_test.Zn = 0