examples.dg.open_bc.lobc_maxwell_2D 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_2D_3subdomain_i2t1e4_n140x10"
list	interior_subdomains = ['interior']
list	tran_subdomains = ['transition']
list	exterior_subdomains = ['exterior']
list	interior_to_transition_boundary = ['leftWallTran']
list	open_boundary = ['rightWallTran']
list	exterior_boundary = ['rightWall']
list	physical_tran_boundaries = ['lowerWallTran','upperWallTran']
list	physical_ext_boundaries = ['lowerWallExt','upperWallExt','rightWall']
list	transition_center = [-1.,0.,0.]
int	transition_width = 1.
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,
bool	cond_wall = True
None	periodic_boundaries = None
	else :
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
list	apps_bc_physics = []
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

app_flux_q

warpy examples.dg.open_bc.lobc_maxwell_2D.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_2D.app_flux_w

app_flux_wlacuna

examples.dg.open_bc.lobc_maxwell_2D.app_flux_wlacuna

apps_bc_physics

list examples.dg.open_bc.lobc_maxwell_2D.apps_bc_physics = []

cond_wall

bool examples.dg.open_bc.lobc_maxwell_2D.cond_wall = True

detect_nonscalar

examples.dg.open_bc.lobc_maxwell_2D.detect_nonscalar

dt

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

dt_controller

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

dt_final

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

else

examples.dg.open_bc.lobc_maxwell_2D.else :

exterior_boundary

list examples.dg.open_bc.lobc_maxwell_2D.exterior_boundary = ['rightWall']

exterior_subdomains

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

fluid_components

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

gen_xdmf

examples.dg.open_bc.lobc_maxwell_2D.gen_xdmf

host_action_lobc_mu

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

host_actions_lobc_pre_ti

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

ic_applicators

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

ic_applicators_lobc

warpy examples.dg.open_bc.lobc_maxwell_2D.ic_applicators_lobc

Initial value:

=  q_lobc_right.vars_setup(
    direction = [1, 0, 0],
    center = transition_center,
    transitionwidth = transition_width
)

interior_bc_apps

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

interior_subdomains

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

interior_to_transition_boundary

list examples.dg.open_bc.lobc_maxwell_2D.interior_to_transition_boundary = ['leftWallTran']

lacuna_time

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

mesh

warpy examples.dg.open_bc.lobc_maxwell_2D.mesh

Initial value:

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

mesh_file

str examples.dg.open_bc.lobc_maxwell_2D.mesh_file = "meshes/mesh_2D_3subdomain_i2t1e4_n140x10"

null_function

warpy examples.dg.open_bc.lobc_maxwell_2D.null_function

Initial value:

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

omega_p_norm

examples.dg.open_bc.lobc_maxwell_2D.omega_p_norm = 1.0

open_boundary

list examples.dg.open_bc.lobc_maxwell_2D.open_boundary = ['rightWallTran']

periodic_boundaries

list examples.dg.open_bc.lobc_maxwell_2D.periodic_boundaries = None

physical_ext_boundaries

list examples.dg.open_bc.lobc_maxwell_2D.physical_ext_boundaries = ['lowerWallExt','upperWallExt','rightWall']

physical_tran_boundaries

list examples.dg.open_bc.lobc_maxwell_2D.physical_tran_boundaries = ['lowerWallTran','upperWallTran']

q

warpy examples.dg.open_bc.lobc_maxwell_2D.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_2D.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 = physical_tran_boundaries# + physical_ext_boundaries
)

---

Lacuna-based open boundary condition (LOBC) setup

reintegration_time

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

sim

warpy examples.dg.open_bc.lobc_maxwell_2D.sim

Initial value:

=  warpy.dg_sim(
    name                     = 'lobc_maxwell_2D',
    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_2D.skin_depth_norm = 1.0

sorder

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

spatial_solvers

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

temporal_solver

warpy examples.dg.open_bc.lobc_maxwell_2D.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_2D.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,
    apps_bc_physics = apps_bc_physics if cond_wall else None
)

torder

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

tran_subdomains

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

transition_center

list examples.dg.open_bc.lobc_maxwell_2D.transition_center = [-1.,0.,0.]

transition_width

int examples.dg.open_bc.lobc_maxwell_2D.transition_width = 1.

True

examples.dg.open_bc.lobc_maxwell_2D.True

use_limiters

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

variable_adjusters

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

variable_name

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

write_steps

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

writer

warpy examples.dg.open_bc.lobc_maxwell_2D.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]
)