examples.dg.mhd.shock_2d_brio-wu Namespace Reference

Functions
def	run_name ()

Variables
str	torder = 'SSPRK3'
int	sord = 4
str	sorder_str = ''
int	cfl = .8 / (2*sord - 1)
float	cfl_diffusion = 1.0
str	num_flux = "Rusanov"
int	dt = 1e-4
int	tend = .1
int	nout = 100
int	num_elem = 16
bool	use_bc_flux_diagnostic = True
str	gradient_method = "ldg"
None	penalty_beta = None
float	penalty_eta = 1.0
int	gamma = 2
float	mass = 1.0
float	charge = 1.0
int	nu_p_tau = 1.
int	omega_p_tau = 1.
int	skin_depth_norm = 1
float	lnlam = 10.0
float	jump_fac = 0.125
int	rho_L = 1
int	P_L = 1
int	v_L = 1
int	By_L = 1
float	Bx = 0.75
float	rho_R = 0.125 / jump_fac
int	P_R = .1 / jump_fac
int	v_R = -1
int	By_R = -1
str	mhd_str = 'imhd'
bool	use_resistivity = 'rmhd' else False
str	resistivity_type = "spitzer-vacuum"
int	max_resistive_diffusivity = 1e-4
float	eta_max = 0.01
float	eta_min = eta_max / 100.0
int	vacuum_diffusivity = max_resistive_diffusivity
float	eta_vacuum = eta_max
int	n_vacuum = 0
bool	use_pos_limiter = True
bool	use_moe_ross_limiter = False
int	moe_ross_alpha = 100
bool	use_viscosity = False
int	visc = -5.e-5
int	kappa = visc
bool	use_art_diss = False
float	coeff_art_dens = 20.0
float	coeff_art_mom = 40.0
float	coeff_art_heat = 40.0
abs	art_diss_diff_cap = abs(visc)
None	dx_cut = None
bool	use_v_anom = False
None	rho_min = None
None	P_min = None
bool	use_imhd_floor = True
bool	use_floors_app = False
list	boundary_names = ['Left', 'Right', 'Bottom', 'Top']
warpy	mesh
warpy	plasma
bool	use_gradients = False
warpy	plasma_gradient
warpy	plasma_gradient2
dict	bc_flux_vars = {}
	name
	bc
	components
	explodable
warpy	ic_plasma
list	ic_va = []
warpy	ic_ha = warpy.host_actions.va_runner(name='ic_ha', variable_adjusters=ic_va)
list	variable_adjusters = []
list	apps = []
list	limiter_apps = []
warpy	floors_app
	else :
list	bc_apps = []
list	bc_apps_dirichlet = []
list	rho = [rho_L, rho_R]
list	px = [mass * v_L, mass * v_R]
list	py = [0.0, 0.0]
list	pz = [0.0, 0.0]
list	By = [By_L, By_R]
list	Bz = [0.0, 0.0]
list	pressure = [P_L, P_R]
list	e = [0, 0]
list	vals = [rho, px, py, pz, e, Bx, By, Bz]
list	grad_bc_apps = []
warpy	spatial_solver
warpy	time_integrator
warpy	dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)
dict	diagnostics_dict = {}
warpy	integ_rho
warpy	diag
int	diag_writer_interval = 1.e-3
warpy	diag_writer
list	write_vars = [plasma]
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=write_vars)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Function Documentation

run_name()

def examples.dg.mhd.shock_2d_brio-wu.run_name

(

)

Variable Documentation

apps

list examples.dg.mhd.shock_2d_brio-wu.apps = []

art_diss_diff_cap

abs examples.dg.mhd.shock_2d_brio-wu.art_diss_diff_cap = abs(visc)

bc

examples.dg.mhd.shock_2d_brio-wu.bc

bc_apps

list examples.dg.mhd.shock_2d_brio-wu.bc_apps = []

bc_apps_dirichlet

list examples.dg.mhd.shock_2d_brio-wu.bc_apps_dirichlet = []

bc_flux_vars

dict examples.dg.mhd.shock_2d_brio-wu.bc_flux_vars = {}

boundary_names

list examples.dg.mhd.shock_2d_brio-wu.boundary_names = ['Left', 'Right', 'Bottom', 'Top']

Bx

list examples.dg.mhd.shock_2d_brio-wu.Bx = 0.75

By

list examples.dg.mhd.shock_2d_brio-wu.By = [By_L, By_R]

By_L

int examples.dg.mhd.shock_2d_brio-wu.By_L = 1

By_R

int examples.dg.mhd.shock_2d_brio-wu.By_R = -1

Bz

list examples.dg.mhd.shock_2d_brio-wu.Bz = [0.0, 0.0]

cfl

int examples.dg.mhd.shock_2d_brio-wu.cfl = .8 / (2*sord - 1)

cfl_diffusion

float examples.dg.mhd.shock_2d_brio-wu.cfl_diffusion = 1.0

charge

float examples.dg.mhd.shock_2d_brio-wu.charge = 1.0

coeff_art_dens

float examples.dg.mhd.shock_2d_brio-wu.coeff_art_dens = 20.0

coeff_art_heat

float examples.dg.mhd.shock_2d_brio-wu.coeff_art_heat = 40.0

coeff_art_mom

float examples.dg.mhd.shock_2d_brio-wu.coeff_art_mom = 40.0

components

examples.dg.mhd.shock_2d_brio-wu.components

detect_nonscalar

examples.dg.mhd.shock_2d_brio-wu.detect_nonscalar

diag

warpy examples.dg.mhd.shock_2d_brio-wu.diag

Initial value:

=  warpy.host_actions.unique_variable_diagnostic(
            f'diagnostic_{bc}',
            ReadVars=[var],
            BasisArraySet=sorder_str,
            stage=0
            )

diag_writer

warpy examples.dg.mhd.shock_2d_brio-wu.diag_writer

Initial value:

=  warpy.host_actions.diagnostics_writer(
    name='diagnostics',
    write_interval=diag_writer_interval,
    diagnostics=diagnostics_dict,
    output_file='diagnostics.csv'
    )

diag_writer_interval

int examples.dg.mhd.shock_2d_brio-wu.diag_writer_interval = 1.e-3

diagnostics_dict

dict examples.dg.mhd.shock_2d_brio-wu.diagnostics_dict = {}

dt

int examples.dg.mhd.shock_2d_brio-wu.dt = 1e-4

dt_controller

warpy examples.dg.mhd.shock_2d_brio-wu.dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)

dx_cut

None examples.dg.mhd.shock_2d_brio-wu.dx_cut = None

e

list examples.dg.mhd.shock_2d_brio-wu.e = [0, 0]

else

examples.dg.mhd.shock_2d_brio-wu.else :

eta_max

float examples.dg.mhd.shock_2d_brio-wu.eta_max = 0.01

eta_min

float examples.dg.mhd.shock_2d_brio-wu.eta_min = eta_max / 100.0

eta_vacuum

float examples.dg.mhd.shock_2d_brio-wu.eta_vacuum = eta_max

explodable

examples.dg.mhd.shock_2d_brio-wu.explodable

floors_app

None examples.dg.mhd.shock_2d_brio-wu.floors_app

Initial value:

=  warpy.apps.functions.mhd.set_density_and_pressure_floors_mhd(
        name="set_density_and_pressure_floors_app",
        mass_density_min=rho_min,
        pressure_min=P_min,
        gamma=gamma,
        fluid=plasma,
        )

gamma

int examples.dg.mhd.shock_2d_brio-wu.gamma = 2

gen_xdmf

examples.dg.mhd.shock_2d_brio-wu.gen_xdmf

grad_bc_apps

list examples.dg.mhd.shock_2d_brio-wu.grad_bc_apps = []

gradient_method

str examples.dg.mhd.shock_2d_brio-wu.gradient_method = "ldg"

ic_ha

warpy examples.dg.mhd.shock_2d_brio-wu.ic_ha = warpy.host_actions.va_runner(name='ic_ha', variable_adjusters=ic_va)

ic_plasma

warpy examples.dg.mhd.shock_2d_brio-wu.ic_plasma

Initial value:

=  warpy.apps.functions.mhd.briowu_mhd(
    name='briowu_ic',
    fluid=plasma,
    gamma=gamma,
    density_left=rho_L,
    density_right=rho_R,
    vx_left=v_L,
    vx_right=v_R,
    by_left=By_L,
    by_right=By_R,
    pressure_left=P_L,
    pressure_right=P_R,
    center=[0, 0, 0],
    transition_width=0.0
    )

ic_va

list examples.dg.mhd.shock_2d_brio-wu.ic_va = []

integ_rho

warpy examples.dg.mhd.shock_2d_brio-wu.integ_rho

Initial value:

=  warpy.host_actions.integrator(
    name='rho_integrator',
    time_integrator=time_integrator,
    patch_process_integrators=[
        warpy.integrators.dg_integrator(
            name='integ_rho',
            spatial_order=sorder_str,
            applications=[
                warpy.apps.integrate.int_q_dv(
                    name='int_rho_dv',
                    variables=plasma,
                    stage=0,
                    # variables_components=['rho']
                    )
                ]
            )
        ]
    )

jump_fac

float examples.dg.mhd.shock_2d_brio-wu.jump_fac = 0.125

kappa

int examples.dg.mhd.shock_2d_brio-wu.kappa = visc

limiter_apps

list examples.dg.mhd.shock_2d_brio-wu.limiter_apps = []

lnlam

float examples.dg.mhd.shock_2d_brio-wu.lnlam = 10.0

mass

float examples.dg.mhd.shock_2d_brio-wu.mass = 1.0

max_resistive_diffusivity

int examples.dg.mhd.shock_2d_brio-wu.max_resistive_diffusivity = 1e-4

mesh

warpy examples.dg.mhd.shock_2d_brio-wu.mesh

Initial value:

=  warpy.mesh.block(
    Bounds=[-1.5, 1.5, -.5, .5],
    NumCells=[num_elem*3, num_elem],
    NodeSets=boundary_names,
    # PeriodicBoundaries=['Bottom', 'Top'],
    NumLayers=1,
    basis_array_set=sorder_str
    )

mhd_str

str examples.dg.mhd.shock_2d_brio-wu.mhd_str = 'imhd'

moe_ross_alpha

int examples.dg.mhd.shock_2d_brio-wu.moe_ross_alpha = 100

n_vacuum

int examples.dg.mhd.shock_2d_brio-wu.n_vacuum = 0

name

examples.dg.mhd.shock_2d_brio-wu.name

nout

int examples.dg.mhd.shock_2d_brio-wu.nout = 100

nu_p_tau

int examples.dg.mhd.shock_2d_brio-wu.nu_p_tau = 1.

num_elem

int examples.dg.mhd.shock_2d_brio-wu.num_elem = 16

num_flux

str examples.dg.mhd.shock_2d_brio-wu.num_flux = "Rusanov"

omega_p_tau

int examples.dg.mhd.shock_2d_brio-wu.omega_p_tau = 1.

P_L

int examples.dg.mhd.shock_2d_brio-wu.P_L = 1

P_min

min examples.dg.mhd.shock_2d_brio-wu.P_min = None

P_R

int examples.dg.mhd.shock_2d_brio-wu.P_R = .1 / jump_fac

penalty_beta

None examples.dg.mhd.shock_2d_brio-wu.penalty_beta = None

penalty_eta

float examples.dg.mhd.shock_2d_brio-wu.penalty_eta = 1.0

plasma

warpy examples.dg.mhd.shock_2d_brio-wu.plasma

Initial value:

=  warpy.variable(
    name='plasma',
    components=['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz'],
    basis_array_set=sorder_str
    )

plasma_gradient

warpy examples.dg.mhd.shock_2d_brio-wu.plasma_gradient

Initial value:

=  warpy.variable(
        name='plasma_gradient',
        components=[
            prefix + "_" + suffix
            for prefix in plasma.components()
            for suffix in ['x', 'y', 'z']
            ],
        basis_array_set=sorder_str
        )

plasma_gradient2

warpy examples.dg.mhd.shock_2d_brio-wu.plasma_gradient2

Initial value:

=  warpy.variable(
        name='plasma_gradient2',
        components=[
            prefix + "_" + suffix
            for prefix in plasma.components()
            for suffix in ['x', 'y', 'z']
            ],
        basis_array_set=sorder_str
        )

pressure

list examples.dg.mhd.shock_2d_brio-wu.pressure = [P_L, P_R]

px

list examples.dg.mhd.shock_2d_brio-wu.px = [mass * v_L, mass * v_R]

py

list examples.dg.mhd.shock_2d_brio-wu.py = [0.0, 0.0]

pz

list examples.dg.mhd.shock_2d_brio-wu.pz = [0.0, 0.0]

resistivity_type

str examples.dg.mhd.shock_2d_brio-wu.resistivity_type = "spitzer-vacuum"

rho

list examples.dg.mhd.shock_2d_brio-wu.rho = [rho_L, rho_R]

rho_L

int examples.dg.mhd.shock_2d_brio-wu.rho_L = 1

rho_min

min examples.dg.mhd.shock_2d_brio-wu.rho_min = None

rho_R

float examples.dg.mhd.shock_2d_brio-wu.rho_R = 0.125 / jump_fac

sim

warpy examples.dg.mhd.shock_2d_brio-wu.sim

Initial value:

=  warpy.dg_sim(
    name=run_name(),
    meshes=[mesh],
    initial_conditions=[ic_ha],
    temporal_solvers=[time_integrator],
    writers=[writer, diag_writer],
    time=[0, tend],
    dt_controller=dt_controller,
    flexible_writeout=True,
    write_steps=nout,
    verbosity='info'
    # verbosity='debug'
    )

skin_depth_norm

int examples.dg.mhd.shock_2d_brio-wu.skin_depth_norm = 1

sord

int examples.dg.mhd.shock_2d_brio-wu.sord = 4

sorder_str

str examples.dg.mhd.shock_2d_brio-wu.sorder_str = ''

spatial_solver

warpy examples.dg.mhd.shock_2d_brio-wu.spatial_solver

Initial value:

=  warpy.spatial_solvers.dg(
    name="dg",
    spatial_order=sorder_str,
    applications=apps,
    penalty_beta=penalty_beta,
    penalty_eta=penalty_eta
    )

tend

int examples.dg.mhd.shock_2d_brio-wu.tend = .1

time_integrator

warpy examples.dg.mhd.shock_2d_brio-wu.time_integrator

Initial value:

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

torder

str examples.dg.mhd.shock_2d_brio-wu.torder = 'SSPRK3'

True

examples.dg.mhd.shock_2d_brio-wu.True

use_art_diss

bool examples.dg.mhd.shock_2d_brio-wu.use_art_diss = False

use_bc_flux_diagnostic

bool examples.dg.mhd.shock_2d_brio-wu.use_bc_flux_diagnostic = True

use_floors_app

bool examples.dg.mhd.shock_2d_brio-wu.use_floors_app = False

use_gradients

bool examples.dg.mhd.shock_2d_brio-wu.use_gradients = False

use_imhd_floor

bool examples.dg.mhd.shock_2d_brio-wu.use_imhd_floor = True

use_moe_ross_limiter

bool examples.dg.mhd.shock_2d_brio-wu.use_moe_ross_limiter = False

use_pos_limiter

bool examples.dg.mhd.shock_2d_brio-wu.use_pos_limiter = True

use_resistivity

bool examples.dg.mhd.shock_2d_brio-wu.use_resistivity = 'rmhd' else False

use_v_anom

bool examples.dg.mhd.shock_2d_brio-wu.use_v_anom = False

use_viscosity

bool examples.dg.mhd.shock_2d_brio-wu.use_viscosity = False

v_L

int examples.dg.mhd.shock_2d_brio-wu.v_L = 1

v_R

int examples.dg.mhd.shock_2d_brio-wu.v_R = -1

vacuum_diffusivity

int examples.dg.mhd.shock_2d_brio-wu.vacuum_diffusivity = max_resistive_diffusivity

vals

list examples.dg.mhd.shock_2d_brio-wu.vals = [rho, px, py, pz, e, Bx, By, Bz]

variable_adjusters

list examples.dg.mhd.shock_2d_brio-wu.variable_adjusters = []

visc

int examples.dg.mhd.shock_2d_brio-wu.visc = -5.e-5

write_vars

list examples.dg.mhd.shock_2d_brio-wu.write_vars = [plasma]

writer

warpy examples.dg.mhd.shock_2d_brio-wu.writer = warpy.host_actions.writer(name='writer', ReadVars=write_vars)