Variables
str	torder = 'RK4'

str	sorder = 'thirdOrder'

float	skin_depth_norm = 1.0

int	gamma = 5./3.

int	z_i = 1

int	A_i = 1

int	theta = 1

bool	useLimiter = True

float	limiter_alpha_coefficient = 500.0

bool	fixed_dt = False

str	verbosity = 'debug'

str	imhd_numerical_flux_type = "HLLD"

float	penalty_beta = -0.5

float	penalty_eta = 0.0

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	briowu_ic = warpy.functions.briowu_mhd(name='briowu_ic',gamma=gamma)

list	ic_applicators = [warpy.applicator(spatial_order=sorder,fun=briowu_ic,var=fluid,components=None,spatial_scheme='Nodal')]

list	apps = []

list	bc_apps_dirichlet = []

list	rho = [1.0,0.125]

list	px = [0.0,0.0]

list	py = [0.0,0.0]

list	pz = [0.0,0.0]

list	Bx = [0.75, 0.75]

list	By = [1.0, -1.0]

list	Bz = [0.0, 0.0]

list	pressure = [1.0,0.1]

list	e = [0,0]

list	bc_apps_neumann = []

list	limiter_apps = []

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, penalty_beta=penalty_beta, penalty_eta=penalty_eta, on_subdomains=['all'])

list	variable_adjusters = []

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

float	dt = 1.0e-9

warpy	dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)

warpy	sim

	gen_xdmf

	True

	detect_nonscalar

Variable Documentation

◆ A_i

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

◆ apps

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

◆ bc_apps_dirichlet

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

◆ bc_apps_neumann

list examples.dg.mhd.brio_wu_shocktube_hall.bc_apps_neumann = []

◆ briowu_ic

warpy examples.dg.mhd.brio_wu_shocktube_hall.briowu_ic = warpy.functions.briowu_mhd(name='briowu_ic',gamma=gamma)

◆ Bx

list examples.dg.mhd.brio_wu_shocktube_hall.Bx = [0.75, 0.75]

◆ By

list examples.dg.mhd.brio_wu_shocktube_hall.By = [1.0, -1.0]

◆ Bz

list examples.dg.mhd.brio_wu_shocktube_hall.Bz = [0.0, 0.0]

◆ component_names

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

◆ detect_nonscalar

examples.dg.mhd.brio_wu_shocktube_hall.detect_nonscalar

◆ dt

float examples.dg.mhd.brio_wu_shocktube_hall.dt = 1.0e-9

◆ dt_controller

warpy examples.dg.mhd.brio_wu_shocktube_hall.dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)

◆ e

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

◆ fixed_dt

bool examples.dg.mhd.brio_wu_shocktube_hall.fixed_dt = False

◆ fluid

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

◆ fluid_gradient

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

◆ gamma

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

◆ gen_xdmf

examples.dg.mhd.brio_wu_shocktube_hall.gen_xdmf

◆ gradient_component_names

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

◆ ic_applicators

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

◆ imhd_numerical_flux_type

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

◆ limiter_alpha_coefficient

float examples.dg.mhd.brio_wu_shocktube_hall.limiter_alpha_coefficient = 500.0

◆ limiter_apps

list examples.dg.mhd.brio_wu_shocktube_hall.limiter_apps = []

◆ mesh

warpy examples.dg.mhd.brio_wu_shocktube_hall.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-.5, .5],
                        # NumCells=[512],
                        # NumCells=[256],                        
                        NumCells=[128],                        
                        # NumCells=[16],                        
                        NodeSets=['leftWall', 'rightWall'],
                        # PeriodicBoundaries=['leftWall','rightWall'],
                        NumLayers=1,
                        basis_array_set=sorder)

◆ penalty_beta

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

◆ penalty_eta

float examples.dg.mhd.brio_wu_shocktube_hall.penalty_eta = 0.0

◆ pressure

list examples.dg.mhd.brio_wu_shocktube_hall.pressure = [1.0,0.1]

◆ px

list examples.dg.mhd.brio_wu_shocktube_hall.px = [0.0,0.0]

◆ py

list examples.dg.mhd.brio_wu_shocktube_hall.py = [0.0,0.0]

◆ pz

list examples.dg.mhd.brio_wu_shocktube_hall.pz = [0.0,0.0]

◆ rho

list examples.dg.mhd.brio_wu_shocktube_hall.rho = [1.0,0.125]

◆ sim

warpy examples.dg.mhd.brio_wu_shocktube_hall.sim

Initial value:

=  warpy.dg_sim(name='primitive_alpha500',                   
                   meshes=[mesh],
                   initial_conditions=ic_applicators,
                   temporal_solvers=[ti],
                   writers=[writer],
                   # time=[0.,0.2],
                   time=[0.,3.0e-4],
                   dt_controller=dt_controller,
                   write_steps=100,                   
                   # write_steps=20000,
                   verbosity=verbosity,
                   num_partitions=1)

◆ skin_depth_norm

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

◆ sorder

str examples.dg.mhd.brio_wu_shocktube_hall.sorder = 'thirdOrder'

◆ spatial_solver

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

◆ theta

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

◆ ti

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

◆ torder

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

◆ True

examples.dg.mhd.brio_wu_shocktube_hall.True

◆ useLimiter

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

◆ variable_adjusters

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

◆ verbosity

str examples.dg.mhd.brio_wu_shocktube_hall.verbosity = 'debug'

◆ writer

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

◆ z_i

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

Variables