examples.dg.5-moment.tf_shock Namespace Reference

Variables
str	torder = 'TVDRK2'
str	sorder = 'secondOrder'
int	gamma = 5/3
float	mass = 1.0
int	Zi = 1
int	Ze = -1
int	skin_depth_norm = 1e-2
int	omega_p_norm = 1e3
int	rho_L = 1
int	P_L = 1/gamma
int	rho_R = 4
int	P_R = 4/gamma
int	Ai = 1
int	Ae = 1/1836
warpy	mesh
warpy	ions
warpy	electrons
warpy	field
warpy	ions_ic
warpy	elec_ic
list	applicators = []
list	apps = []
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=[ions, electrons, field])
warpy	spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps)
list	variable_adjusters = []
warpy	time_integrator = warpy.host_actions.erk(name="ti", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)
int	dt = 1e-6
int	nout = 100
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

Ae

int examples.dg.5-moment.tf_shock.Ae = 1/1836

Ai

int examples.dg.5-moment.tf_shock.Ai = 1

applicators

list examples.dg.5-moment.tf_shock.applicators = []

apps

list examples.dg.5-moment.tf_shock.apps = []

detect_nonscalar

examples.dg.5-moment.tf_shock.detect_nonscalar

dt

int examples.dg.5-moment.tf_shock.dt = 1e-6

dt_controller

warpy examples.dg.5-moment.tf_shock.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

elec_ic

warpy examples.dg.5-moment.tf_shock.elec_ic

Initial value:

=  warpy.functions.shock_tube_5moment(name='elec_shock_tube',
                                             density_left=Ae*rho_L,
                                             density_right=Ae*rho_R,
                                             pressure_left=P_L*.5,
                                             pressure_right=P_R*.5,
                                             gamma=gamma,
                                             center=[5,0,0])

electrons

warpy examples.dg.5-moment.tf_shock.electrons

Initial value:

=  warpy.variable(name='electrons',
                           components=['rho', 'px', 'py', 'pz', 'e'],
                           basis_array_set=sorder)

field

warpy examples.dg.5-moment.tf_shock.field

Initial value:

=  warpy.variable(name='field',
                       components=['Ex','Ey','Ez', 'Bx','By','Bz'],
                       basis_array_set=sorder)

gamma

int examples.dg.5-moment.tf_shock.gamma = 5/3

gen_xdmf

examples.dg.5-moment.tf_shock.gen_xdmf

ions

warpy examples.dg.5-moment.tf_shock.ions

Initial value:

=  warpy.variable(name='ion',
                      components=['rho', 'px', 'py', 'pz', 'e'],
                      basis_array_set=sorder)

ions_ic

warpy examples.dg.5-moment.tf_shock.ions_ic

Initial value:

=  warpy.functions.shock_tube_5moment(name='ion_shock_tube',
                                             density_left=Ai*rho_L,
                                             density_right=Ai*rho_R,
                                             pressure_left=P_L*.5,
                                             pressure_right=P_R*.5,
                                             gamma=gamma,
                                             center=[5,0,0])

mass

float examples.dg.5-moment.tf_shock.mass = 1.0

mesh

warpy examples.dg.5-moment.tf_shock.mesh

Initial value:

=  warpy.mesh.block(Bounds=[0, 10],
                        NumCells=[128],
                        NodeSets=['Left', 'Right'],
                        NumLayers=1,
                        basis_array_set=sorder)

nout

int examples.dg.5-moment.tf_shock.nout = 100

omega_p_norm

int examples.dg.5-moment.tf_shock.omega_p_norm = 1e3

P_L

int examples.dg.5-moment.tf_shock.P_L = 1/gamma

P_R

int examples.dg.5-moment.tf_shock.P_R = 4/gamma

rho_L

int examples.dg.5-moment.tf_shock.rho_L = 1

rho_R

int examples.dg.5-moment.tf_shock.rho_R = 4

sim

warpy examples.dg.5-moment.tf_shock.sim

Initial value:

=  warpy.dg_sim(name='tf_shock',
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[time_integrator],
                   writers=[writer],
                   time=[0,1],
                   dt_controller = dt_controller,
                   flexible_writeout=True,
                   write_steps=nout,
                   verbosity='info')

skin_depth_norm

int examples.dg.5-moment.tf_shock.skin_depth_norm = 1e-2

sorder

str examples.dg.5-moment.tf_shock.sorder = 'secondOrder'

spatial_solver

warpy examples.dg.5-moment.tf_shock.spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps)

time_integrator

warpy examples.dg.5-moment.tf_shock.time_integrator = warpy.host_actions.erk(name="ti", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)

torder

str examples.dg.5-moment.tf_shock.torder = 'TVDRK2'

True

examples.dg.5-moment.tf_shock.True

variable_adjusters

list examples.dg.5-moment.tf_shock.variable_adjusters = []

writer

warpy examples.dg.5-moment.tf_shock.writer = warpy.host_actions.writer(name='writer', ReadVars=[ions, electrons, field])

Ze

int examples.dg.5-moment.tf_shock.Ze = -1

Zi

int examples.dg.5-moment.tf_shock.Zi = 1