Variables
str	torder = 'RK4'

str	sorder = 'secondOrder'

bool	useLimiter = True

int	limiter_alpha = 5

float	cfl = 0.8

str	imhd_numerical_flux_type = "HLLD"

float	penalty_beta = -0.0

float	penalty_eta = 0.0

str	gradient_method = "ldg"

float	perturbation = 0.1

float	floor_density = 1.0E-2

float	floor_pressure = 1.0E-2

float	skin_depth_norm = 1.0

None	nu_p_norm = None

int	gamma = 5./3.

int	z_i = 1

int	A_i = 1

int	theta = 1

int	di_L = 1.

int	Lx = 8math.pidi_L

int	Ly = 4math.pidi_L

float	sheet_width = 0.5

float	n0 = 1.0

float	n_inf = 0.2

float	B0 = 1.0

float	Ti_Te = 1.0/theta

float	Te = B0**2/2.0/n0/(Ti_Te + 1.0)

float	Ti = Ti_Te * Te

list	T = [Ti, Te]

warpy	mesh

list	component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']

warpy	fluid

list	gradient_component_names = [prefix+"_"+suffix for prefix in component_names for suffix in ['x','y','z']]

warpy	fluid_gradient

warpy	fluid_gradient2

warpy	gem_ic

list	applicators

list	apps = []

list	dirichlet_apps = []

list	neumann_apps = []

list	limiter_apps = []

warpy	writer

warpy	spatial_solver

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

	local_path

	gen_xdmf

	True

	detect_nonscalar

Variable Documentation

◆ A_i

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

◆ applicators

list examples.dg.mhd.gem_challenge_hall2.applicators

Initial value:

=  [warpy.applicator(spatial_order=sorder,
                                fun=gem_ic,
                                var=fluid,
                                spatial_scheme='Nodal',
                                components=component_names[0:8])]

◆ apps

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

◆ B0

float examples.dg.mhd.gem_challenge_hall2.B0 = 1.0

◆ cfl

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

◆ component_names

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

◆ detect_nonscalar

examples.dg.mhd.gem_challenge_hall2.detect_nonscalar

◆ di_L

int examples.dg.mhd.gem_challenge_hall2.di_L = 1.

◆ dirichlet_apps

list examples.dg.mhd.gem_challenge_hall2.dirichlet_apps = []

◆ dt_controller

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

◆ floor_density

float examples.dg.mhd.gem_challenge_hall2.floor_density = 1.0E-2

◆ floor_pressure

float examples.dg.mhd.gem_challenge_hall2.floor_pressure = 1.0E-2

◆ fluid

warpy examples.dg.mhd.gem_challenge_hall2.fluid

Initial value:

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

◆ fluid_gradient

warpy examples.dg.mhd.gem_challenge_hall2.fluid_gradient

Initial value:

=  warpy.variable(name='fluid_gradient',
                                components=gradient_component_names,
                                basis_array_set=sorder)

◆ fluid_gradient2

warpy examples.dg.mhd.gem_challenge_hall2.fluid_gradient2

Initial value:

=  warpy.variable(name='fluid_gradient2',
                                 components=gradient_component_names,
                                 basis_array_set=sorder,
                                 explodable=False)

◆ gamma

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

◆ gem_ic

warpy examples.dg.mhd.gem_challenge_hall2.gem_ic

Initial value:

=  warpy.functions.gem_hallmhd(name='gem_ic',
                                     gamma=gamma,
                                     mass=A_i,
                                     temperature=Ti+Te,
                                     sheet_width=sheet_width,
                                     number_density=n0,
                                     number_density_infinity=n_inf,
                                     B0=B0,
                                     Lx=Lx,
                                     Ly=Ly,
                                     perturbation=perturbation)

◆ gen_xdmf

examples.dg.mhd.gem_challenge_hall2.gen_xdmf

◆ gradient_component_names

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

◆ gradient_method

str examples.dg.mhd.gem_challenge_hall2.gradient_method = "ldg"

◆ imhd_numerical_flux_type

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

◆ limiter_alpha

int examples.dg.mhd.gem_challenge_hall2.limiter_alpha = 5

◆ limiter_apps

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

◆ local_path

examples.dg.mhd.gem_challenge_hall2.local_path

◆ Lx

int examples.dg.mhd.gem_challenge_hall2.Lx = 8*math.pi*di_L

◆ Ly

int examples.dg.mhd.gem_challenge_hall2.Ly = 4*math.pi*di_L

◆ mesh

warpy examples.dg.mhd.gem_challenge_hall2.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-Lx/2.,Lx/2., -Ly/2.,Ly/2.],
                        # NumCells=[512,256],                                                
                        # NumCells=[256,128],                        
                        # NumCells=[128,64],
                        # NumCells=[64,32],                        
                        NumCells=[32,16],
                        # NumCells=[16,8],
                        # NumCells=[8,4],                                                
                        NodeSets=['Left','Right', 'Lower','Upper'],
                        NumLayers=1,
                        basis_array_set=sorder,
                        PeriodicBoundaries=['Left','Right'],
                        block_type="symmetric")

◆ n0

float examples.dg.mhd.gem_challenge_hall2.n0 = 1.0

◆ n_inf

float examples.dg.mhd.gem_challenge_hall2.n_inf = 0.2

◆ neumann_apps

list examples.dg.mhd.gem_challenge_hall2.neumann_apps = []

◆ nu_p_norm

None examples.dg.mhd.gem_challenge_hall2.nu_p_norm = None

◆ penalty_beta

float examples.dg.mhd.gem_challenge_hall2.penalty_beta = -0.0

◆ penalty_eta

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

◆ perturbation

float examples.dg.mhd.gem_challenge_hall2.perturbation = 0.1

◆ sheet_width

float examples.dg.mhd.gem_challenge_hall2.sheet_width = 0.5

◆ sim

warpy examples.dg.mhd.gem_challenge_hall2.sim

Initial value:

=  warpy.dg_sim(name='test_no_glm',   
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[ti],
                   writers=[writer],
                   time=[0,4],
                   # time=[0,40.0],
                   dt_controller=dt_controller,
                   write_steps=100,
                   verbosity='debug'
                   )

◆ skin_depth_norm

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

◆ sorder

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

◆ spatial_solver

warpy examples.dg.mhd.gem_challenge_hall2.spatial_solver

Initial value:

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

◆ T

list examples.dg.mhd.gem_challenge_hall2.T = [Ti, Te]

◆ Te

float examples.dg.mhd.gem_challenge_hall2.Te = B0**2/2.0/n0/(Ti_Te + 1.0)

◆ theta

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

◆ Ti

float examples.dg.mhd.gem_challenge_hall2.Ti = Ti_Te * Te

◆ ti

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

◆ Ti_Te

float examples.dg.mhd.gem_challenge_hall2.Ti_Te = 1.0/theta

◆ torder

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

◆ True

examples.dg.mhd.gem_challenge_hall2.True

◆ useLimiter

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

◆ variable_adjusters

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

◆ writer

warpy examples.dg.mhd.gem_challenge_hall2.writer

Initial value:

=  warpy.host_actions.writer(name='writer',
                                   ReadVars=[fluid,
                                             fluid_gradient])

◆ z_i

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

Variables

Variable Documentation

◆ A_i

◆ applicators

◆ apps

◆ B0

◆ cfl

◆ component_names

◆ detect_nonscalar

◆ di_L

◆ dirichlet_apps

◆ dt_controller

◆ floor_density

◆ floor_pressure

◆ fluid

◆ fluid_gradient

◆ fluid_gradient2

◆ gamma

◆ gem_ic

◆ gen_xdmf

◆ gradient_component_names

◆ gradient_method

◆ imhd_numerical_flux_type

◆ limiter_alpha

◆ limiter_apps

◆ local_path

◆ Lx

◆ Ly

◆ mesh

◆ n0

◆ n_inf

◆ neumann_apps

◆ nu_p_norm

◆ penalty_beta

◆ penalty_eta

◆ perturbation

◆ sheet_width

◆ sim

◆ skin_depth_norm

◆ sorder

◆ spatial_solver

◆ T

◆ Te

◆ theta

◆ Ti

◆ ti

◆ Ti_Te

◆ torder

◆ True

◆ useLimiter

◆ variable_adjusters

◆ writer

◆ z_i