landau_damping 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
str	pert_strength = "weak"
float	perturbation = 0.01
	else :
float	thermal_velocity = 1.0
float	density = 1.0
list	v_space_num_elements = [40]
list	min_v = [-5.0]
list	dv = [0.25]
int	num_physical_elements = 40
int	L = 2 * np.pi
float	k = 0.5
warpy	mesh
	Generate a mesh -> vlasov.
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-5
	Simulation (dg)
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")
f	sim_name = f'landau_damping_{num_physical_elements}x{v_space_num_elements[0]}_{pert_strength}'
warpy	config_args = warpy.load_configs()
warpy	case_loc = config_args['test_rundir'](f'kinetics/landau_damping/{now}', '')
list	max_vx = min_v[0] + dv[0] * v_space_num_elements[0]
f	sim_param
warpy	sim
	gen_xdmf
	Run Simulation ... 1 now = datetime.datetime.now()
	detect_nonscalar

Variable Documentation

applications_moments

list landau_damping.applications_moments = []

apps_kinetic

list landau_damping.apps_kinetic

Initial value:

=  [warpy.apps.vlasovmaxwell.vlasovmaxwell(name = 'vlasov',
                                 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 landau_damping.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 landau_damping.case_loc = config_args['test_rundir'](f'kinetics/landau_damping/{now}', '')

cfl

float landau_damping.cfl = 0.235

charge

int landau_damping.charge = 1

charge_to_mass

int landau_damping.charge_to_mass = charge/mass

config_args

warpy landau_damping.config_args = warpy.load_configs()

density

float landau_damping.density = 1.0

detect_nonscalar

landau_damping.detect_nonscalar

distribution

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

Initial value:

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

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

distribution_ic_va

warpy landau_damping.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 landau_damping.dt = 1e-5

Simulation (dg)

dt_controller

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

dv

list landau_damping.dv = [0.25]

else

landau_damping.else :

field

warpy landau_damping.field

Initial value:

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

fluid

warpy landau_damping.fluid

Initial value:

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

gen_xdmf

landau_damping.gen_xdmf

Run Simulation ... 1 now = datetime.datetime.now()

ic_ha

warpy landau_damping.ic_ha

Initial value:

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

k

float landau_damping.k = 0.5

L

int landau_damping.L = 2 * np.pi

mass

int landau_damping.mass = 1

max_vx

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

maxwell_ic

warpy landau_damping.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 landau_damping.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 landau_damping.mesh

Initial value:

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

Generate a mesh -> vlasov.

min_v

list landau_damping.min_v = [-5.0]

now

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

num_physical_elements

int landau_damping.num_physical_elements = 40

omega_p_norm

float landau_damping.omega_p_norm = 1.0

periodic_boundary_conditions

bool landau_damping.periodic_boundary_conditions = True

pert_strength

str landau_damping.pert_strength = "weak"

perturbation

float landau_damping.perturbation = 0.01

sim

warpy landau_damping.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],
                   # time = [0,dt*write_steps],                   
                   dt_controller = dt_controller,
                   flexible_writeout = False,
                   write_steps = write_steps,
                   verbosity = 'info',
                   # num_partitions=4
)

sim_name

f landau_damping.sim_name = f'landau_damping_{num_physical_elements}x{v_space_num_elements[0]}_{pert_strength}'

sim_param

landau_damping.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 landau_damping.skin_depth_norm = 1.0

sorder

str landau_damping.sorder = 'thirdOrder'

spatial_solver_kinetic

warpy landau_damping.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 landau_damping.spatial_solver_physical

Initial value:

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

t_final

int landau_damping.t_final = 80

temporal_solver

warpy landau_damping.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 landau_damping.thermal_velocity = 1.0

torder

str landau_damping.torder = 'RK4'

v_space_num_elements

list landau_damping.v_space_num_elements = [40]

variable_adjusters

list landau_damping.variable_adjusters = []

Variable adjusters.

write_steps

int landau_damping.write_steps = 1000

writer

warpy landau_damping.writer

Initial value:

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

Variable writer.