two_stream_instability Namespace Reference

Variables
str	torder = 'RK4'
str	sorder = 'thirdOrder'
bool	periodic_boundary_conditions = True
float	skin_depth_norm = 1.0
float	omega_p_norm = 1.0
int	charge = 1
int	mass = 1
int	charge_to_mass = charge/mass
float	cfl = 0.235
float	perturbation = 0.01
float	thermal_velocity = 0.1
float	density = 1.0
np	velocity_separation = np.pi/2.0
list	v_space_num_elements = [160]
list	min_v = [-10.0]
list	dv = [0.125]
int	num_physical_elements = 80
int	L = 2 * np.pi
float	k = 0.5
warpy	mesh
warpy	distribution
warpy	field
warpy	fluid
warpy	distribution_ic
	Initial conditions-------------------------------— Step 1: define the function as an application:
warpy	maxwell_ic
warpy	distribution_ic_va
warpy	maxwell_ic_va
warpy	ic_ha
list	apps_kinetic
	Initial conditions-------------------------------—.
list	apps_physical
list	applications_moments = []
list	variable_adjusters = []
	Variable adjusters.
warpy	writer
	Variable writer.
warpy	spatial_solver_kinetic
	Spatial solver.
warpy	spatial_solver_physical
warpy	temporal_solver
	Temporal solver.
int	dt = 1e-1
	Simulation (dg) dt = 1e-5.
int	t_final = 80
warpy	dt_controller = warpy.dt_calc.stability_dt(dt)
int	write_steps = 1000
datetime	now = datetime.datetime.now().strftime("%Y_%m_%d")
str	sim_name = '80x160_3rdorder_rk4_t80_1000writes_vth_0p1_sep_pi_over_2'
warpy	config_args = warpy.load_configs()
warpy	case_loc = config_args['test_rundir'](f'kinetics/two_stream_instability/{now}', '')
list	max_vx = min_v[0] + dv[0] * v_space_num_elements[0]
f	sim_param
warpy	sim
	gen_xdmf
	Run Simulation ... 1.
	detect_nonscalar

Variable Documentation

applications_moments

list two_stream_instability.applications_moments = []

apps_kinetic

list two_stream_instability.apps_kinetic

Initial value:

=  [warpy.apps.vlasovmaxwell.vlasovmaxwell(name = 'advection_squares',
                                 pdf_variable = distribution,
                                 field_variable = field,
                                 skin_depth_norm = skin_depth_norm,
                                 charge_to_mass = charge_to_mass)]

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

Applications (consists of equation sets and boundary conditions)

apps_physical

list two_stream_instability.apps_physical

Initial value:

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

case_loc

warpy two_stream_instability.case_loc = config_args['test_rundir'](f'kinetics/two_stream_instability/{now}', '')

cfl

float two_stream_instability.cfl = 0.235

charge

int two_stream_instability.charge = 1

charge_to_mass

int two_stream_instability.charge_to_mass = charge/mass

config_args

warpy two_stream_instability.config_args = warpy.load_configs()

density

float two_stream_instability.density = 1.0

detect_nonscalar

two_stream_instability.detect_nonscalar

distribution

warpy two_stream_instability.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 two_stream_instability.distribution_ic

Initial value:

=  warpy.apps.functions.kinetics.vlasov_1d1v_two_stream(name = 'distribution_ic',
                                                                       distribution=distribution,
                                                         thermal_velocity = thermal_velocity,
                                                         density = density,
                                                         perturbation = perturbation,
                                                         wave_number = k,
                                                         velocity_separation = velocity_separation
)

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

distribution_ic_va

warpy two_stream_instability.distribution_ic_va

Initial value:

=  warpy.variable_adjusters.function_evaluation.phase_space_function_evaluator(name='distribution_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])

dt

int two_stream_instability.dt = 1e-1

Simulation (dg) dt = 1e-5.

dt_controller

warpy two_stream_instability.dt_controller = warpy.dt_calc.stability_dt(dt)

dv

list two_stream_instability.dv = [0.125]

field

warpy two_stream_instability.field

Initial value:

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

fluid

warpy two_stream_instability.fluid

Initial value:

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

gen_xdmf

two_stream_instability.gen_xdmf

Run Simulation ... 1.

ic_ha

warpy two_stream_instability.ic_ha

Initial value:

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

k

float two_stream_instability.k = 0.5

L

int two_stream_instability.L = 2 * np.pi

mass

int two_stream_instability.mass = 1

max_vx

list two_stream_instability.max_vx = min_v[0] + dv[0] * v_space_num_elements[0]

maxwell_ic

warpy two_stream_instability.maxwell_ic

Initial value:

=  warpy.apps.functions.maxwell.field_1d1v_landau_damping('maxwell_ic',
                                                                    fields=field,
                                                     omega_p_tau=omega_p_norm,
                                                     skin_depth_norm=skin_depth_norm,
                                                     density=density,
                                                     charge=charge,
                                                     perturbation = perturbation,
                                                     wave_number = k
                                                     )

maxwell_ic_va

warpy two_stream_instability.maxwell_ic_va

Initial value:

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

mesh

warpy two_stream_instability.mesh

Initial value:

=  warpy.mesh.block(Bounds = [0, 2*np.pi/k],
                        NumCells = [num_physical_elements],
                        NodeSets = ['Left', 'Right'],
                        NumLayers = 1,
                        PeriodicBoundaries = ['Left','Right'] if periodic_boundary_conditions else None,
                        basis_array_set = sorder)

min_v

list two_stream_instability.min_v = [-10.0]

now

datetime two_stream_instability.now = datetime.datetime.now().strftime("%Y_%m_%d")

num_physical_elements

int two_stream_instability.num_physical_elements = 80

omega_p_norm

float two_stream_instability.omega_p_norm = 1.0

periodic_boundary_conditions

bool two_stream_instability.periodic_boundary_conditions = True

perturbation

float two_stream_instability.perturbation = 0.01

sim

warpy two_stream_instability.sim

Initial value:

=  warpy.dg_sim(name = sim_name,
                   meshes = [mesh],
                   initial_conditions = [ic_ha],
                   temporal_solvers = [temporal_solver],
                   writers = [writer],
                   time = [0,t_final],
                   dt_controller = dt_controller,
                   flexible_writeout = False,
                   write_steps = write_steps,
                   verbosity = 'info',
                   # num_partitions=4
)

sim_name

str two_stream_instability.sim_name = '80x160_3rdorder_rk4_t80_1000writes_vth_0p1_sep_pi_over_2'

sim_param

two_stream_instability.sim_param

Initial value:

=  \
f'''case_loc={case_loc}
specific_case={sim_name}
start_frame=0
end_frame={write_steps}
domain=domain
variable={distribution.name()}
v_space_num_vx={v_space_num_elements[0]}
basis_name={sorder}
min_vx={min_v[0]}
max_vx={max_vx}
'''

skin_depth_norm

float two_stream_instability.skin_depth_norm = 1.0

sorder

str two_stream_instability.sorder = 'thirdOrder'

spatial_solver_kinetic

warpy two_stream_instability.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.

spatial_solver_physical

warpy two_stream_instability.spatial_solver_physical

Initial value:

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

t_final

int two_stream_instability.t_final = 80

temporal_solver

warpy two_stream_instability.temporal_solver

Initial value:

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

Temporal solver.

thermal_velocity

float two_stream_instability.thermal_velocity = 0.1

torder

str two_stream_instability.torder = 'RK4'

v_space_num_elements

list two_stream_instability.v_space_num_elements = [160]

variable_adjusters

list two_stream_instability.variable_adjusters = []

Variable adjusters.

velocity_separation

np two_stream_instability.velocity_separation = np.pi/2.0

write_steps

int two_stream_instability.write_steps = 1000

writer

warpy two_stream_instability.writer

Initial value:

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

Variable writer.