examples.dg.5-moment.shock_tube Namespace Reference

Variables
str	torder = 'RK4'
str	sorder = 'secondOrder'
int	gamma = 5/3
float	mass = 1.0
float	charge = 1.0
None	nu_p_tau = None
int	rho_L = 1
int	P_L = 1/gamma
int	rho_R = 4
int	P_R = 4/gamma
float	penalty_beta = 0.0
float	penalty_eta = 0.0
warpy	mesh = warpy.mesh.block(Bounds=[0, 10],NumCells=[16],NodeSets=['Left', 'Right'],NumLayers=1,basis_array_set=sorder)
warpy	fluid
warpy	fluid_gradient
warpy	shock_tube_function
list	applicators = []
list	apps = []
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid, fluid_gradient])
	else :
warpy	spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps,penalty_beta=penalty_beta,penalty_eta=penalty_eta)
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.stability_dt(init_dt=dt)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

applicators

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

apps

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

charge

float examples.dg.5-moment.shock_tube.charge = 1.0

detect_nonscalar

examples.dg.5-moment.shock_tube.detect_nonscalar

dt

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

dt_controller

warpy examples.dg.5-moment.shock_tube.dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)

else

examples.dg.5-moment.shock_tube.else :

fluid

warpy examples.dg.5-moment.shock_tube.fluid

Initial value:

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

fluid_gradient

warpy examples.dg.5-moment.shock_tube.fluid_gradient

Initial value:

=  warpy.variable(name='fluid_gradient',
                                    components=[prefix+"_"+suffix for prefix in fluid.components() for suffix in ['x','y','z']],
                                    basis_array_set=sorder)

gamma

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

gen_xdmf

examples.dg.5-moment.shock_tube.gen_xdmf

mass

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

mesh

warpy examples.dg.5-moment.shock_tube.mesh = warpy.mesh.block(Bounds=[0, 10],NumCells=[16],NodeSets=['Left', 'Right'],NumLayers=1,basis_array_set=sorder)

nout

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

nu_p_tau

None examples.dg.5-moment.shock_tube.nu_p_tau = None

P_L

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

P_R

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

penalty_beta

float examples.dg.5-moment.shock_tube.penalty_beta = 0.0

penalty_eta

float examples.dg.5-moment.shock_tube.penalty_eta = 0.0

rho_L

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

rho_R

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

shock_tube_function

warpy examples.dg.5-moment.shock_tube.shock_tube_function

Initial value:

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

sim

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

Initial value:

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

sorder

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

spatial_solver

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

time_integrator

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

torder

str examples.dg.5-moment.shock_tube.torder = 'RK4'

True

examples.dg.5-moment.shock_tube.True

variable_adjusters

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

writer

warpy examples.dg.5-moment.shock_tube.writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid, fluid_gradient])