1d1v_double_rarefaction Namespace Reference

Variables
str	torder = 'TVDRK2'
str	sorder = 'secondOrder'
float	cfl = 0.1
float	gamma = 1.4
float	skin_depth_norm = 1.0
float	omega_p_norm = 1.0
float	charge = 0.0
int	mass = 1
float	charge_to_mass = charge/mass
float	thermal_velocity = 1.0
float	density = 1.0
list	v_space_num_elements = [40]
list	min_v = [-10]
list	dv = [0.5]
int	num_physical_elements = 128
bool	periodic_boundary_conditions = False
int	nu_p_tau = 1E-6
int	nu = 1
warpy	mesh
warpy	distribution
warpy	field
warpy	fluid
warpy	distribution_ic
	Initial conditions-------------------------------— Step 1: define the function as an application:
warpy	ic_va
warpy	ic_ha
list	apps_kinetic = []
	Initial conditions-------------------------------—.
list	apps_physical
list	applications_moments = []
list	applications_2nd_moments = []
list	bc_apps_dirichlet_field = []
list	Ex = [0.0,0.0]
list	Ey = [0.0,0.0]
list	Ez = [0.0,0.0]
list	Bx = [0.0,0.0]
list	By = [0.0,0.0]
list	Bz = [0.0,0.0]
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=[distribution, fluid, field])
warpy	spatial_solver_kinetic
warpy	spatial_solver_physical
list	variable_adjusters = []
warpy	time_integrator
	Variable adjusters variable_adjusters = [].
int	dt = 1e-4
int	nout = 100
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
warpy	sim
datetime	now = datetime.datetime.now()
	category
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

applications_2nd_moments

list 1d1v_double_rarefaction.applications_2nd_moments = []

applications_moments

list 1d1v_double_rarefaction.applications_moments = []

apps_kinetic

list 1d1v_double_rarefaction.apps_kinetic = []

Initial conditions-------------------------------—.

Applications (consists of equation sets and boundary conditions)

apps_physical

list 1d1v_double_rarefaction.apps_physical

Initial value:

=  [warpy.apps.maxwell.flux(name='maxwell',
                                         field=field,
                                         skin_depth_norm=skin_depth_norm,
                                         omega_p_norm=omega_p_norm)]

bc_apps_dirichlet_field

list 1d1v_double_rarefaction.bc_apps_dirichlet_field = []

Bx

list 1d1v_double_rarefaction.Bx = [0.0,0.0]

By

list 1d1v_double_rarefaction.By = [0.0,0.0]

Bz

list 1d1v_double_rarefaction.Bz = [0.0,0.0]

category

1d1v_double_rarefaction.category

cfl

float 1d1v_double_rarefaction.cfl = 0.1

charge

float 1d1v_double_rarefaction.charge = 0.0

charge_to_mass

float 1d1v_double_rarefaction.charge_to_mass = charge/mass

density

float 1d1v_double_rarefaction.density = 1.0

detect_nonscalar

1d1v_double_rarefaction.detect_nonscalar

distribution

warpy 1d1v_double_rarefaction.distribution

Initial value:

=  warpy.variable(name = 'distribution',
                              components = ['f'],
                              basis_array_set = sorder,
                              phase_space_basis_array_set = sorder,
                              v_space_num_elements = v_space_num_elements)

distribution_ic

warpy 1d1v_double_rarefaction.distribution_ic

Initial value:

=  warpy.apps.functions.kinetics.kinetics1d1v_double_rarefaction(name='double_rarefaction_problem_1d1v',
                                                                                distribution=distribution,
                                                                                mass=mass)

Initial conditions-------------------------------— Step 1: define the function as an application:

dt

int 1d1v_double_rarefaction.dt = 1e-4

dt_controller

warpy 1d1v_double_rarefaction.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

dv

list 1d1v_double_rarefaction.dv = [0.5]

Ex

list 1d1v_double_rarefaction.Ex = [0.0,0.0]

Ey

list 1d1v_double_rarefaction.Ey = [0.0,0.0]

Ez

list 1d1v_double_rarefaction.Ez = [0.0,0.0]

field

warpy 1d1v_double_rarefaction.field

Initial value:

=  warpy.variable(name = 'field',
                       components = ['Ex','Ey','Ez','Bx','By','Bz'],
                       basis_array_set = sorder)

fluid

warpy 1d1v_double_rarefaction.fluid

Initial value:

=  warpy.variable(name = 'fluid',
                       components = ['n','nvx','nvy','nvz','nT'],
                       basis_array_set = sorder,
                       explodable=False)

gamma

float 1d1v_double_rarefaction.gamma = 1.4

gen_xdmf

1d1v_double_rarefaction.gen_xdmf

ic_ha

warpy 1d1v_double_rarefaction.ic_ha

Initial value:

=  warpy.host_actions.va_runner(name='ic_ha',
                                     variable_adjusters=[ic_va])

ic_va

warpy 1d1v_double_rarefaction.ic_va

Initial value:

=  warpy.variable_adjusters.function_evaluation.phase_space_function_evaluator(name='ic_va',
                                                                                    priority=0,
                                                                                    spatial_order=sorder,
                                                                                    phase_space_spatial_order=sorder,
                                                                                    v_space_num_elements=v_space_num_elements,
                                                                                    min_v=min_v,
                                                                                    dv=dv,
                                                                                    on_subdomains=['all'],
                                                                                    applications=[distribution_ic])

mass

int 1d1v_double_rarefaction.mass = 1

mesh

warpy 1d1v_double_rarefaction.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-0.5, 0.5],
                        NumCells=[num_physical_elements],
                        NodeSets=['Left', 'Right'],
                        PeriodicBoundaries = ['Left','Right'] if periodic_boundary_conditions else None,                        
                        NumLayers=1,
                        basis_array_set=sorder)

min_v

list 1d1v_double_rarefaction.min_v = [-10]

nout

int 1d1v_double_rarefaction.nout = 100

now

datetime 1d1v_double_rarefaction.now = datetime.datetime.now()

nu

int 1d1v_double_rarefaction.nu = 1

nu_p_tau

int 1d1v_double_rarefaction.nu_p_tau = 1E-6

num_physical_elements

int 1d1v_double_rarefaction.num_physical_elements = 128

omega_p_norm

float 1d1v_double_rarefaction.omega_p_norm = 1.0

periodic_boundary_conditions

bool 1d1v_double_rarefaction.periodic_boundary_conditions = False

sim

warpy 1d1v_double_rarefaction.sim

Initial value:

=  warpy.dg_sim(name='test',                                                         
                   meshes=[mesh],
                   initial_conditions=[ic_ha],
                   temporal_solvers=[time_integrator],
                   writers=[writer],
                   time=[0,0.15],
                   # time=[0,dt*nout],                   
                   dt_controller = dt_controller,
                   flexible_writeout=True,
                   write_steps=nout,
                   verbosity='info')

skin_depth_norm

float 1d1v_double_rarefaction.skin_depth_norm = 1.0

sorder

str 1d1v_double_rarefaction.sorder = 'secondOrder'

spatial_solver_kinetic

warpy 1d1v_double_rarefaction.spatial_solver_kinetic

Initial value:

=  warpy.spatial_solvers.dg_kinetic(name="dg_kinetic",
                                                  spatial_order = sorder,
                                                  phase_space_spatial_order = sorder,
                                                  v_space_num_elements = v_space_num_elements,
                                                  min_v = min_v,
                                                  dv = dv,
                                                  applications = apps_kinetic,
                                                  on_subdomains = ['all'],
                                                  cfl = cfl)

spatial_solver_physical

warpy 1d1v_double_rarefaction.spatial_solver_physical

Initial value:

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

thermal_velocity

float 1d1v_double_rarefaction.thermal_velocity = 1.0

time_integrator

warpy 1d1v_double_rarefaction.time_integrator

Initial value:

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

Variable adjusters variable_adjusters = [].

torder

str 1d1v_double_rarefaction.torder = 'TVDRK2'

True

1d1v_double_rarefaction.True

v_space_num_elements

list 1d1v_double_rarefaction.v_space_num_elements = [40]

variable_adjusters

list 1d1v_double_rarefaction.variable_adjusters = []

writer

warpy 1d1v_double_rarefaction.writer = warpy.host_actions.writer(name='writer', ReadVars=[distribution, fluid, field])