examples.dg.5-moment.shock_1d_sod Namespace Reference

Functions
def	run_name ()

Variables
str	torder = 'SSPRK3'
int	sord = 3
int	integral_calculation_stage = 1
	else :
str	sorder = ''
float	cfl = 0.8 / (2*sord - 1)
float	cfl_diffusion = 1.0
str	num_flux = "Roe"
bool	use_entropy_fix = True
int	dt = 1e-4
int	tend = .2
int	nout = 100
str	dt_controller_type = 'fixed'
int	num_elem = 100
str	gradient_method = "ldg"
None	penalty_beta = None
float	penalty_eta = 1.0
float	gamma = 1.4
float	mass = 1.0
float	charge = 1.0
int	nu_p_tau = 1.
int	jump_fac = 1
int	rho_L = 1
int	P_L = 1
int	v_L = 0
float	rho_R = 0.125 / jump_fac
float	P_R = 0.1 / jump_fac
int	v_R = 0
bool	use_pos_limiter = False
bool	use_moe_ross_limiter = False
int	moe_ross_alpha = 100
bool	use_art_diss = False
float	coeff_art_dens = 20.0
float	coeff_art_mom = 40.0
float	coeff_art_heat = 40.0
int	art_diss_diff_cap = 1e-4
int	dx_cut = 1e-3
bool	use_v_anom = False
bool	use_consv_floor_default = True
bool	use_floors_app = False
None	rho_min = None
None	P_min = None
warpy	mesh
list	fluid_comps = ['rho', 'px', 'py', 'pz', 'e']
warpy	fluid = warpy.variable(name='fluid', components=fluid_comps, basis_array_set=sorder)
warpy	fluid_analytic
bool	use_gradients = False
warpy	fluid_gradient
warpy	fluid_gradient2
warpy	ic_fluid
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
list	bc_apps = []
list	grad_bc_apps = []
list	analytic_solution_apps = []
list	integrator_apps = []
warpy	spatial_solver
warpy	time_integrator
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
warpy	writer
warpy	integrator_host_action
warpy	diag_writer
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Function Documentation

run_name()

def examples.dg.5-moment.shock_1d_sod.run_name

(

)

Variable Documentation

analytic_solution_apps

list examples.dg.5-moment.shock_1d_sod.analytic_solution_apps = []

apps

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

art_diss_diff_cap

int examples.dg.5-moment.shock_1d_sod.art_diss_diff_cap = 1e-4

bc_apps

list examples.dg.5-moment.shock_1d_sod.bc_apps = []

cfl

float examples.dg.5-moment.shock_1d_sod.cfl = 0.8 / (2*sord - 1)

cfl_diffusion

float examples.dg.5-moment.shock_1d_sod.cfl_diffusion = 1.0

charge

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

coeff_art_dens

float examples.dg.5-moment.shock_1d_sod.coeff_art_dens = 20.0

coeff_art_heat

float examples.dg.5-moment.shock_1d_sod.coeff_art_heat = 40.0

coeff_art_mom

float examples.dg.5-moment.shock_1d_sod.coeff_art_mom = 40.0

detect_nonscalar

examples.dg.5-moment.shock_1d_sod.detect_nonscalar

diag_writer

warpy examples.dg.5-moment.shock_1d_sod.diag_writer

Initial value:

=  warpy.host_actions.diagnostics_writer(
    name='diagnostics',
    write_interval=1e-3,
    diagnostics={integrator_host_action: ['integral_analytic']},
    output_file='diagnostics.csv')

dt

int examples.dg.5-moment.shock_1d_sod.dt = 1e-4

dt_controller

warpy examples.dg.5-moment.shock_1d_sod.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

dt_controller_type

str examples.dg.5-moment.shock_1d_sod.dt_controller_type = 'fixed'

dx_cut

int examples.dg.5-moment.shock_1d_sod.dx_cut = 1e-3

else

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

floors_app

None examples.dg.5-moment.shock_1d_sod.floors_app

Initial value:

=  warpy.apps.functions.five_moment.set_density_and_pressure_floors_single_species_5moment(
        name="set_density_and_pressure_floors_app",
        mass_density_min=rho_min,
        pressure_min=P_min,
        gamma=gamma,
        fluid=fluid,
        )

fluid

warpy examples.dg.5-moment.shock_1d_sod.fluid = warpy.variable(name='fluid', components=fluid_comps, basis_array_set=sorder)

fluid_analytic

warpy examples.dg.5-moment.shock_1d_sod.fluid_analytic

Initial value:

=  warpy.variable(
    name='fluid_analytic', components=fluid_comps, basis_array_set=sorder
    )

fluid_comps

list examples.dg.5-moment.shock_1d_sod.fluid_comps = ['rho', 'px', 'py', 'pz', 'e']

fluid_gradient

warpy examples.dg.5-moment.shock_1d_sod.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
        )

fluid_gradient2

warpy examples.dg.5-moment.shock_1d_sod.fluid_gradient2

Initial value:

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

gamma

float examples.dg.5-moment.shock_1d_sod.gamma = 1.4

gen_xdmf

examples.dg.5-moment.shock_1d_sod.gen_xdmf

grad_bc_apps

list examples.dg.5-moment.shock_1d_sod.grad_bc_apps = []

gradient_method

str examples.dg.5-moment.shock_1d_sod.gradient_method = "ldg"

ic_fluid

warpy examples.dg.5-moment.shock_1d_sod.ic_fluid

Initial value:

=  warpy.apps.functions.five_moment.euler1d_arbitrary_shocktube_5moment(
    name='ic_shock',
    gamma=gamma,
    fluid=fluid,
    density_left=rho_L,
    velocity_left=v_L,
    pressure_left=P_L,
    density_right=rho_R,
    velocity_right=v_R,
    pressure_right=P_R,
    xc=0.0,
    transition_width=0.0
    )

ic_ha

warpy examples.dg.5-moment.shock_1d_sod.ic_ha = warpy.host_actions.va_runner(name='ic_ha', variable_adjusters=ic_va)

ic_va

list examples.dg.5-moment.shock_1d_sod.ic_va = []

integral_calculation_stage

int examples.dg.5-moment.shock_1d_sod.integral_calculation_stage = 1

integrator_apps

list examples.dg.5-moment.shock_1d_sod.integrator_apps = []

integrator_host_action

warpy examples.dg.5-moment.shock_1d_sod.integrator_host_action

Initial value:

=  warpy.host_actions.integrator(
    name='mass_integrator',
    time_integrator=time_integrator,
    patch_process_integrators=[
        warpy.integrators.dg_integrator(
            name='integrator',
            spatial_order=sorder,
            applications=integrator_apps)
    ])

jump_fac

int examples.dg.5-moment.shock_1d_sod.jump_fac = 1

limiter_apps

list examples.dg.5-moment.shock_1d_sod.limiter_apps = []

mass

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

mesh

warpy examples.dg.5-moment.shock_1d_sod.mesh

Initial value:

=  warpy.mesh.block(
    Bounds=[-.5, .5],
    NumCells=[num_elem],
    NodeSets=['Left', 'Right'],
    NumLayers=1,
    basis_array_set=sorder
    )

moe_ross_alpha

int examples.dg.5-moment.shock_1d_sod.moe_ross_alpha = 100

nout

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

nu_p_tau

int examples.dg.5-moment.shock_1d_sod.nu_p_tau = 1.

num_elem

int examples.dg.5-moment.shock_1d_sod.num_elem = 100

num_flux

str examples.dg.5-moment.shock_1d_sod.num_flux = "Roe"

P_L

int examples.dg.5-moment.shock_1d_sod.P_L = 1

P_min

min examples.dg.5-moment.shock_1d_sod.P_min = None

P_R

float examples.dg.5-moment.shock_1d_sod.P_R = 0.1 / jump_fac

penalty_beta

None examples.dg.5-moment.shock_1d_sod.penalty_beta = None

penalty_eta

float examples.dg.5-moment.shock_1d_sod.penalty_eta = 1.0

rho_L

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

rho_min

min examples.dg.5-moment.shock_1d_sod.rho_min = None

rho_R

float examples.dg.5-moment.shock_1d_sod.rho_R = 0.125 / jump_fac

sim

warpy examples.dg.5-moment.shock_1d_sod.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'
    )

sord

int examples.dg.5-moment.shock_1d_sod.sord = 3

sorder

str examples.dg.5-moment.shock_1d_sod.sorder = ''

spatial_solver

warpy examples.dg.5-moment.shock_1d_sod.spatial_solver

Initial value:

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

tend

int examples.dg.5-moment.shock_1d_sod.tend = .2

time_integrator

warpy examples.dg.5-moment.shock_1d_sod.time_integrator

Initial value:

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

torder

str examples.dg.5-moment.shock_1d_sod.torder = 'SSPRK3'

True

examples.dg.5-moment.shock_1d_sod.True

use_art_diss

bool examples.dg.5-moment.shock_1d_sod.use_art_diss = False

use_consv_floor_default

bool examples.dg.5-moment.shock_1d_sod.use_consv_floor_default = True

use_entropy_fix

bool examples.dg.5-moment.shock_1d_sod.use_entropy_fix = True

use_floors_app

bool examples.dg.5-moment.shock_1d_sod.use_floors_app = False

use_gradients

bool examples.dg.5-moment.shock_1d_sod.use_gradients = False

use_moe_ross_limiter

bool examples.dg.5-moment.shock_1d_sod.use_moe_ross_limiter = False

use_pos_limiter

bool examples.dg.5-moment.shock_1d_sod.use_pos_limiter = False

use_v_anom

bool examples.dg.5-moment.shock_1d_sod.use_v_anom = False

v_L

int examples.dg.5-moment.shock_1d_sod.v_L = 0

v_R

int examples.dg.5-moment.shock_1d_sod.v_R = 0

variable_adjusters

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

writer

warpy examples.dg.5-moment.shock_1d_sod.writer

Initial value:

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