examples.dg.5-moment.forward_facing_step Namespace Reference

Variables
str	torder = 'TVDRK2'
str	sorder = 'secondOrder'
float	cfl = 0.9
float	gamma = 1.4
float	limiter_alpha_coefficient = 0.0
bool	use_limiter = True
float	density_min = 0.01
float	pressure_min = 0.01
warpy	mesh
warpy	fluid
warpy	shock_tube_function
warpy	shock_tube_function_va
list	applicators = []
list	solver_apps = []
list	boundary_condition_apps = []
float	rhoin = gamma
float	pxin = 3.0*rhoin
float	pyin = 0.0
float	pzin = 0.0
float	pressurein = 1.0
float	energyin = pressurein/(gamma-1) + 1.0/2.0*(pxin**2 + pyin**2 + pzin**2)/rhoin
list	limiter_apps = []
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid])
warpy	spatial_solver
list	variable_adjusters = []
warpy	time_integrator = warpy.host_actions.erk(name="ti", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)
int	dt = 1e-6
int	nout = 100
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=0.00134634557913)
warpy	sim
datetime	now = datetime.datetime.now()
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

applicators

list examples.dg.5-moment.forward_facing_step.applicators = []

boundary_condition_apps

list examples.dg.5-moment.forward_facing_step.boundary_condition_apps = []

cfl

float examples.dg.5-moment.forward_facing_step.cfl = 0.9

density_min

float examples.dg.5-moment.forward_facing_step.density_min = 0.01

detect_nonscalar

examples.dg.5-moment.forward_facing_step.detect_nonscalar

dt

int examples.dg.5-moment.forward_facing_step.dt = 1e-6

dt_controller

warpy examples.dg.5-moment.forward_facing_step.dt_controller = warpy.dt_calc.fixed_dt(init_dt=0.00134634557913)

energyin

float examples.dg.5-moment.forward_facing_step.energyin = pressurein/(gamma-1) + 1.0/2.0*(pxin**2 + pyin**2 + pzin**2)/rhoin

fluid

warpy examples.dg.5-moment.forward_facing_step.fluid

Initial value:

=  warpy.variable(name='fluid',
                       components=['rho', 'px', 'py', 'pz', 'e'],
                       basis_array_set=sorder)

gamma

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

gen_xdmf

examples.dg.5-moment.forward_facing_step.gen_xdmf

limiter_alpha_coefficient

float examples.dg.5-moment.forward_facing_step.limiter_alpha_coefficient = 0.0

limiter_apps

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

mesh

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

Initial value:

=  warpy.mesh.arbitrary(
    # MeshFile        = "fstepA001.neu",
    # MeshFile        = "FS_949.neu",
    # MeshFile        = "forward3A00004.neu",
    # MeshFile        = "iman_step.inp",
    MeshFile        = "iman_step_coarse.inp",        
    NumLayers       = 1,
    # PeriodicBoundaries=['Wall', 'Inflow', 'Outflow'],
    basis_array_set =sorder
    )

nout

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

now

datetime examples.dg.5-moment.forward_facing_step.now = datetime.datetime.now()

pressure_min

float examples.dg.5-moment.forward_facing_step.pressure_min = 0.01

pressurein

float examples.dg.5-moment.forward_facing_step.pressurein = 1.0

pxin

float examples.dg.5-moment.forward_facing_step.pxin = 3.0*rhoin

pyin

float examples.dg.5-moment.forward_facing_step.pyin = 0.0

pzin

float examples.dg.5-moment.forward_facing_step.pzin = 0.0

rhoin

float examples.dg.5-moment.forward_facing_step.rhoin = gamma

shock_tube_function

warpy examples.dg.5-moment.forward_facing_step.shock_tube_function

Initial value:

=  warpy.apps.functions.euler2d_forward_facing_step_5moment(name='rp1',
                                                                               gamma=gamma,
                                                                               fluid=fluid)

shock_tube_function_va

warpy examples.dg.5-moment.forward_facing_step.shock_tube_function_va

Initial value:

=  warpy.variable_adjusters.function_evaluation.function_evaluator(name="rp1_va",
                                                                                         priority=0,
                                                                                         spatial_order=sorder,
                                                                                         on_subdomains=['all'],
                                                                                         applications=[shock_tube_function])

sim

warpy examples.dg.5-moment.forward_facing_step.sim

Initial value:

=  warpy.dg_sim(name='test',                   
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[time_integrator],
                   writers=[writer],
                   # time=[0,dt*nout],                   
                   time=[0,4.0],
                   # time=[0,0.25],                   
                   dt_controller = dt_controller,
                   flexible_writeout=True,
                   write_steps=nout,
                   verbosity='debug',
                   enable_fpe=True,
                   num_partitions=1)

solver_apps

list examples.dg.5-moment.forward_facing_step.solver_apps = []

sorder

str examples.dg.5-moment.forward_facing_step.sorder = 'secondOrder'

spatial_solver

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

Initial value:

=  warpy.spatial_solvers.dg(name="dg",
                                          spatial_order=sorder,
                                          applications=solver_apps,
                                          cfl = cfl)

time_integrator

warpy examples.dg.5-moment.forward_facing_step.time_integrator = warpy.host_actions.erk(name="ti", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)

torder

str examples.dg.5-moment.forward_facing_step.torder = 'TVDRK2'

True

examples.dg.5-moment.forward_facing_step.True

use_limiter

bool examples.dg.5-moment.forward_facing_step.use_limiter = True

variable_adjusters

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

writer

warpy examples.dg.5-moment.forward_facing_step.writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid])