examples.dg.advection.advection_squares Namespace Reference

Variables
str	torder = 'RK1'
str	sorder = 'firstOrder'
list	velocity = [1.0,0.0,0.0]
bool	periodic_boundary_conditions = True
bool	use_limiters = False
int	number_velocity_components = 10
int	skin_depth = 1
int	charge_to_mass = 1
str	numerical_flux_type = "upwind"
float	cfl = 0.8
str	perturbation = "sinusoidal"
warpy	mesh
	Generate a mesh.
np	component_names = np.zeros(number_velocity_components)
	Create a variable for the fluid and one for the array Make components array.
warpy	fluid
warpy	field
warpy	fluid_ic
	Initial conditions function TODO: Make an initial condition function that creates phase space distributions Probable parameters: Number of components, thermal velocity, density, type of perturbation fluid_ic = warpy.functions.fourier(name='fluid_ic', sine_amplitudes=[1],sine_frequencies=[2*np.pi])
list	applicators
	Initial condition applicator.
list	apps
	Applications (consists of equation sets and boundary conditions)
list	bc_apps = []
warpy	writer
	Variable writer.
warpy	spatial_solver
	Spatial solver.
list	variable_adjusters = []
	Variable adjusters.
warpy	temporal_solver
	Temporal solver.
int	dt = 1e-3
	Simulation (dg)
float	dt_final = 2.0
warpy	dt_controller = warpy.dt_calc.fixed_dt(dt)
int	write_steps = 100
warpy	sim
	gen_xdmf
	Run Simulation ...
	detect_nonscalar

Variable Documentation

applicators

list examples.dg.advection.advection_squares.applicators

Initial value:

=  [warpy.applicator(spatial_order = sorder,
                                fun = fluid_ic,
                                var = fluid,
                                spatial_scheme = 'Nodal')]

Initial condition applicator.

apps

list examples.dg.advection.advection_squares.apps

Initial value:

=  [warpy.apps.vlasovmaxwell(name = 'advection',
                                 pdf_variable = fluid,
                                 field_variable = field,
                                 skin_depth = skin_depth,
                                 charge_to_mass = charge_to_mass,
                                 number_velocity_components = number_velocity_components,
                                 numerical_flux_type = numerical_flux_type)]

Applications (consists of equation sets and boundary conditions)

bc_apps

list examples.dg.advection.advection_squares.bc_apps = []

cfl

float examples.dg.advection.advection_squares.cfl = 0.8

charge_to_mass

int examples.dg.advection.advection_squares.charge_to_mass = 1

component_names

np examples.dg.advection.advection_squares.component_names = np.zeros(number_velocity_components)

Create a variable for the fluid and one for the array Make components array.

detect_nonscalar

examples.dg.advection.advection_squares.detect_nonscalar

dt

int examples.dg.advection.advection_squares.dt = 1e-3

Simulation (dg)

dt_controller

warpy examples.dg.advection.advection_squares.dt_controller = warpy.dt_calc.fixed_dt(dt)

dt_final

float examples.dg.advection.advection_squares.dt_final = 2.0

field

warpy examples.dg.advection.advection_squares.field

Initial value:

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

fluid

warpy examples.dg.advection.advection_squares.fluid

Initial value:

=  warpy.variable(name = 'fluid',
                       components = component_names,
                       basis_array_set = sorder)

fluid_ic

warpy examples.dg.advection.advection_squares.fluid_ic

Initial value:

=  warpy.functions.distribution(name = 'fluid_ic',
                                        number_velocity_components = number_velocity_components,
                                        thermal_velocity = thermal_velocity,
                                        density = density,
                                        perturbation = perturbation)

Initial conditions function TODO: Make an initial condition function that creates phase space distributions Probable parameters: Number of components, thermal velocity, density, type of perturbation fluid_ic = warpy.functions.fourier(name='fluid_ic', sine_amplitudes=[1],sine_frequencies=[2*np.pi])

gen_xdmf

examples.dg.advection.advection_squares.gen_xdmf

Run Simulation ...

mesh

warpy examples.dg.advection.advection_squares.mesh

Initial value:

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

Generate a mesh.

number_velocity_components

int examples.dg.advection.advection_squares.number_velocity_components = 10

numerical_flux_type

str examples.dg.advection.advection_squares.numerical_flux_type = "upwind"

periodic_boundary_conditions

bool examples.dg.advection.advection_squares.periodic_boundary_conditions = True

perturbation

str examples.dg.advection.advection_squares.perturbation = "sinusoidal"

sim

warpy examples.dg.advection.advection_squares.sim

Initial value:

=  warpy.dg_sim(name = 'advection',
                   meshes = [mesh],
                   initial_conditions = applicators,
                   temporal_solvers = [temporal_solver],
                   writers = [writer],
                   time = [0,dt_final],
                   dt_controller = dt_controller,
                   flexible_writeout = False,
                   write_steps = write_steps,
                   verbosity = 'info')

skin_depth

int examples.dg.advection.advection_squares.skin_depth = 1

sorder

str examples.dg.advection.advection_squares.sorder = 'firstOrder'

spatial_solver

warpy examples.dg.advection.advection_squares.spatial_solver

Initial value:

=  warpy.spatial_solvers.dg(name="dg",
                                          spatial_order = sorder,
                                          applications = apps,
                                          on_subdomains = ['all'],
                                          phase_space_basis_primitive_id = "WMPRIMITIVE_KINETIC_SQUARE",
                                          cfl = cfl)

Spatial solver.

temporal_solver

warpy examples.dg.advection.advection_squares.temporal_solver

Initial value:

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

Temporal solver.

torder

str examples.dg.advection.advection_squares.torder = 'RK1'

use_limiters

bool examples.dg.advection.advection_squares.use_limiters = False

variable_adjusters

list examples.dg.advection.advection_squares.variable_adjusters = []

Variable adjusters.

velocity

list examples.dg.advection.advection_squares.velocity = [1.0,0.0,0.0]

write_steps

int examples.dg.advection.advection_squares.write_steps = 100

writer

warpy examples.dg.advection.advection_squares.writer

Initial value:

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

Variable writer.