examples.dg.mhd.dai_woodward_shocktube Namespace Reference

Variables
str	torder = 'RK4'
str	sorder = 'secondOrder'
float	gamma = 5.0/3.0
bool	useLimiter = True
bool	always_limit = True
str	numerical_flux_type = "Roe"
float	cfl = 0.8
warpy	mesh = warpy.mesh.block(Bounds=[-.5, .5],NumCells=[512],NodeSets=['leftWall', 'rightWall'],NumLayers=1,basis_array_set=sorder)
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	daiwoodward_ic = warpy.functions.daiwoodward_mhd(name='daiwoodward_ic',gamma=gamma,transition_width=0.0)
list	ic_applicators = [warpy.applicator(spatial_order=sorder,fun=daiwoodward_ic,var=fluid,components=None,spatial_scheme='Nodal')]
list	apps = []
list	bc_apps = []
list	bc_apps_copyout = []
list	bc_apps_dirichlet = []
list	rho = [1.08,1.0]
list	pressure = [0.95,1.0]
list	vx = [1.2,0.0]
list	vy = [0.01,0.0]
list	vz = [0.5,0.0]
list	Bx = [4.0/math.sqrt(4.0*math.pi),4.0/math.sqrt(4.0*math.pi)]
list	By = [3.6/math.sqrt(4.0*math.pi),4.0/math.sqrt(4.0*math.pi)]
list	Bz = [2.0/math.sqrt(4.0*math.pi),2.0/math.sqrt(4.0*math.pi)]
list	px = [0,0]
list	py = [0,0]
list	pz = [0,0]
list	e = [0,0]
warpy	writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid,fluid_gradient])
warpy	spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps, on_subdomains=['all'],cfl=cfl)
list	variable_adjusters = []
warpy	ti = warpy.host_actions.erk(name="rk", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)
warpy	dt_controller = warpy.dt_calc.stability_dt(init_dt=1e-5)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

always_limit

bool examples.dg.mhd.dai_woodward_shocktube.always_limit = True

apps

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

bc_apps

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

bc_apps_copyout

list examples.dg.mhd.dai_woodward_shocktube.bc_apps_copyout = []

bc_apps_dirichlet

list examples.dg.mhd.dai_woodward_shocktube.bc_apps_dirichlet = []

Bx

list examples.dg.mhd.dai_woodward_shocktube.Bx = [4.0/math.sqrt(4.0*math.pi),4.0/math.sqrt(4.0*math.pi)]

By

list examples.dg.mhd.dai_woodward_shocktube.By = [3.6/math.sqrt(4.0*math.pi),4.0/math.sqrt(4.0*math.pi)]

Bz

list examples.dg.mhd.dai_woodward_shocktube.Bz = [2.0/math.sqrt(4.0*math.pi),2.0/math.sqrt(4.0*math.pi)]

cfl

float examples.dg.mhd.dai_woodward_shocktube.cfl = 0.8

component_names

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

daiwoodward_ic

warpy examples.dg.mhd.dai_woodward_shocktube.daiwoodward_ic = warpy.functions.daiwoodward_mhd(name='daiwoodward_ic',gamma=gamma,transition_width=0.0)

detect_nonscalar

examples.dg.mhd.dai_woodward_shocktube.detect_nonscalar

dt_controller

warpy examples.dg.mhd.dai_woodward_shocktube.dt_controller = warpy.dt_calc.stability_dt(init_dt=1e-5)

e

list examples.dg.mhd.dai_woodward_shocktube.e = [0,0]

fluid

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

fluid_gradient

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

gamma

float examples.dg.mhd.dai_woodward_shocktube.gamma = 5.0/3.0

gen_xdmf

examples.dg.mhd.dai_woodward_shocktube.gen_xdmf

gradient_component_names

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

ic_applicators

list examples.dg.mhd.dai_woodward_shocktube.ic_applicators = [warpy.applicator(spatial_order=sorder,fun=daiwoodward_ic,var=fluid,components=None,spatial_scheme='Nodal')]

mesh

warpy examples.dg.mhd.dai_woodward_shocktube.mesh = warpy.mesh.block(Bounds=[-.5, .5],NumCells=[512],NodeSets=['leftWall', 'rightWall'],NumLayers=1,basis_array_set=sorder)

numerical_flux_type

str examples.dg.mhd.dai_woodward_shocktube.numerical_flux_type = "Roe"

pressure

list examples.dg.mhd.dai_woodward_shocktube.pressure = [0.95,1.0]

px

list examples.dg.mhd.dai_woodward_shocktube.px = [0,0]

py

list examples.dg.mhd.dai_woodward_shocktube.py = [0,0]

pz

list examples.dg.mhd.dai_woodward_shocktube.pz = [0,0]

rho

list examples.dg.mhd.dai_woodward_shocktube.rho = [1.08,1.0]

sim

warpy examples.dg.mhd.dai_woodward_shocktube.sim

Initial value:

=  warpy.dg_sim(name='daiwoodward_july_14_2nd_order_testing_mixed_glm',
                   meshes=[mesh],
                   initial_conditions=ic_applicators,
                   temporal_solvers=[ti],
                   writers=[writer],
                   time=[0.,0.2],
                   # time=[0.,0.002],                   
                   # time=[0.,1.e-5*1000],
                   dt_controller=dt_controller,
                   flexible_writeout=False,
                   write_steps=100,
                   verbosity='info')

sorder

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

spatial_solver

warpy examples.dg.mhd.dai_woodward_shocktube.spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps, on_subdomains=['all'],cfl=cfl)

ti

warpy examples.dg.mhd.dai_woodward_shocktube.ti = warpy.host_actions.erk(name="rk", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)

torder

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

True

examples.dg.mhd.dai_woodward_shocktube.True

useLimiter

bool examples.dg.mhd.dai_woodward_shocktube.useLimiter = True

variable_adjusters

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

vx

list examples.dg.mhd.dai_woodward_shocktube.vx = [1.2,0.0]

vy

list examples.dg.mhd.dai_woodward_shocktube.vy = [0.01,0.0]

vz

list examples.dg.mhd.dai_woodward_shocktube.vz = [0.5,0.0]

writer

warpy examples.dg.mhd.dai_woodward_shocktube.writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid,fluid_gradient])