examples.dg.open_bc.lobc_maxwell_1D Namespace Reference

Variables
str	torder = 'RK2'
str	sorder = 'thirdOrder'
int	dt = 1e-3
int	dt_final = 8
warpy	dt_controller = warpy.dt_calc.fixed_dt(dt)
int	write_steps = 100
str	mesh_file = "meshes/mesh_1D_3subdomain-n501"
list	interior_subdomains = ['interior']
list	tran_subdomains = ['transition']
list	exterior_subdomains = ['exterior']
list	interior_to_transition_boundary = ['tLeft']
list	open_boundary = ['tRight']
list	exterior_boundary = ['Right']
list	all_tran_boundaries = interior_to_transition_boundary + open_boundary
float	skin_depth_norm = 1.0
float	omega_p_norm = 1.0
bool	use_limiters = False
str	variable_name = 'field'
list	fluid_components = ['Ex','Ey','Ez','Bx','By','Bz']
float	lacuna_time = 1.0,
float	reintegration_time = 3.0,
warpy	mesh
warpy	q
warpy	null_function
list	ic_applicators = []
warpy	app_flux_q
list	interior_bc_apps = []
list	spatial_solvers = []
list	variable_adjusters = []
warpy	temporal_solver
warpy	q_lobc_right
warpy	ic_applicators_lobc
warpy	host_action_lobc_mu = q_lobc_right.calc_dmu()
	app_flux_w
	app_flux_wlacuna
warpy	temporal_solver_lobc
warpy	host_actions_lobc_pre_ti = q_lobc_right.create_pre_ti_hostactions()
warpy	writer
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

all_tran_boundaries

list examples.dg.open_bc.lobc_maxwell_1D.all_tran_boundaries = interior_to_transition_boundary + open_boundary

app_flux_q

warpy examples.dg.open_bc.lobc_maxwell_1D.app_flux_q

Initial value:

=  warpy.apps.maxwell.flux(
    name = 'app_flux_q',
    field = q,
    skin_depth_norm = skin_depth_norm,
    omega_p_norm = omega_p_norm
)

app_flux_w

examples.dg.open_bc.lobc_maxwell_1D.app_flux_w

app_flux_wlacuna

examples.dg.open_bc.lobc_maxwell_1D.app_flux_wlacuna

detect_nonscalar

examples.dg.open_bc.lobc_maxwell_1D.detect_nonscalar

dt

int examples.dg.open_bc.lobc_maxwell_1D.dt = 1e-3

dt_controller

warpy examples.dg.open_bc.lobc_maxwell_1D.dt_controller = warpy.dt_calc.fixed_dt(dt)

dt_final

int examples.dg.open_bc.lobc_maxwell_1D.dt_final = 8

exterior_boundary

list examples.dg.open_bc.lobc_maxwell_1D.exterior_boundary = ['Right']

exterior_subdomains

list examples.dg.open_bc.lobc_maxwell_1D.exterior_subdomains = ['exterior']

fluid_components

list examples.dg.open_bc.lobc_maxwell_1D.fluid_components = ['Ex','Ey','Ez','Bx','By','Bz']

gen_xdmf

examples.dg.open_bc.lobc_maxwell_1D.gen_xdmf

host_action_lobc_mu

warpy examples.dg.open_bc.lobc_maxwell_1D.host_action_lobc_mu = q_lobc_right.calc_dmu()

host_actions_lobc_pre_ti

warpy examples.dg.open_bc.lobc_maxwell_1D.host_actions_lobc_pre_ti = q_lobc_right.create_pre_ti_hostactions()

ic_applicators

list examples.dg.open_bc.lobc_maxwell_1D.ic_applicators = []

ic_applicators_lobc

warpy examples.dg.open_bc.lobc_maxwell_1D.ic_applicators_lobc

Initial value:

=  q_lobc_right.vars_setup(
    direction = [1, 0, 0],
    center = [2.5,0.,0.],
    transitionwidth = .25
)

interior_bc_apps

list examples.dg.open_bc.lobc_maxwell_1D.interior_bc_apps = []

interior_subdomains

list examples.dg.open_bc.lobc_maxwell_1D.interior_subdomains = ['interior']

interior_to_transition_boundary

list examples.dg.open_bc.lobc_maxwell_1D.interior_to_transition_boundary = ['tLeft']

lacuna_time

float examples.dg.open_bc.lobc_maxwell_1D.lacuna_time = 1.0,

mesh

warpy examples.dg.open_bc.lobc_maxwell_1D.mesh

Initial value:

=  warpy.mesh.arbitrary(
    MeshFile        = str(mesh_file + '.inp'), 
    NumLayers       = 1, 
    PeriodicBoundaries = None,
    basis_array_set =sorder
    )

mesh_file

str examples.dg.open_bc.lobc_maxwell_1D.mesh_file = "meshes/mesh_1D_3subdomain-n501"

null_function

warpy examples.dg.open_bc.lobc_maxwell_1D.null_function

Initial value:

=  warpy.functions.set_to(
    name            = 'null_func',
    value           = 0.
)

omega_p_norm

examples.dg.open_bc.lobc_maxwell_1D.omega_p_norm = 1.0

open_boundary

list examples.dg.open_bc.lobc_maxwell_1D.open_boundary = ['tRight']

q

warpy examples.dg.open_bc.lobc_maxwell_1D.q

Initial value:

=  warpy.variable(
    name            = variable_name,
    components      = fluid_components,
    subdomains      = interior_subdomains + tran_subdomains,
    virtual_subdomains = exterior_subdomains,
    basis_array_set = sorder
)

q_lobc_right

warpy examples.dg.open_bc.lobc_maxwell_1D.q_lobc_right

Initial value:

=  warpy.helpers.lobc(
    interior_variable = q,
    sorder = sorder,
    lacuna_time = lacuna_time,
    reintegration_time = reintegration_time,
    interior_subdomains = interior_subdomains,
    tran_subdomains = tran_subdomains,
    exterior_subdomains = exterior_subdomains,
    interior_to_transition_boundary = interior_to_transition_boundary,
    open_boundary = open_boundary,
    exterior_boundary = exterior_boundary,
    other_tran_subdomain_boundaries = None    
)

---

Lacuna-based open boundary condition (LOBC) setup

reintegration_time

float examples.dg.open_bc.lobc_maxwell_1D.reintegration_time = 3.0,

sim

warpy examples.dg.open_bc.lobc_maxwell_1D.sim

Initial value:

=  warpy.dg_sim(
    name                     = 'lobc_maxwell_1D',
    meshes                   = [mesh],
    initial_conditions       = ic_applicators + ic_applicators_lobc,
    extra_start_host_actions = [host_action_lobc_mu],
    temporal_solvers         = [temporal_solver, temporal_solver_lobc],
    pre_ti_host_actions      = host_actions_lobc_pre_ti, 
    writers                  = [writer], 
    time                     = [0, dt_final],
    dt_controller            = dt_controller,
    flexible_writeout        = False,
    write_steps              = write_steps,
    verbosity                = 'info'
)

skin_depth_norm

examples.dg.open_bc.lobc_maxwell_1D.skin_depth_norm = 1.0

sorder

str examples.dg.open_bc.lobc_maxwell_1D.sorder = 'thirdOrder'

spatial_solvers

list examples.dg.open_bc.lobc_maxwell_1D.spatial_solvers = []

temporal_solver

warpy examples.dg.open_bc.lobc_maxwell_1D.temporal_solver

Initial value:

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

temporal_solver_lobc

warpy examples.dg.open_bc.lobc_maxwell_1D.temporal_solver_lobc

Initial value:

=  q_lobc_right.create_temporal_solver(
    torder = torder,
    apps_flux_interior_var = [app_flux_q],
    apps_source_interior_var = None,
    apps_flux_aux_vars = [app_flux_w, app_flux_wlacuna],
    apps_source_aux_vars = None,
)

torder

str examples.dg.open_bc.lobc_maxwell_1D.torder = 'RK2'

tran_subdomains

list examples.dg.open_bc.lobc_maxwell_1D.tran_subdomains = ['transition']

True

examples.dg.open_bc.lobc_maxwell_1D.True

use_limiters

bool examples.dg.open_bc.lobc_maxwell_1D.use_limiters = False

variable_adjusters

list examples.dg.open_bc.lobc_maxwell_1D.variable_adjusters = []

variable_name

str examples.dg.open_bc.lobc_maxwell_1D.variable_name = 'field'

write_steps

int examples.dg.open_bc.lobc_maxwell_1D.write_steps = 100

writer

warpy examples.dg.open_bc.lobc_maxwell_1D.writer

Initial value:

=  warpy.host_actions.writer(
    name = 'writer',
    ReadVars = [q, q_lobc_right.w, q_lobc_right.wlacuna, q_lobc_right.mu, q_lobc_right.dmu]
)