Namespaces
namespace	sheath_flux_bc

Variables
str	sheath_flux_bc.torder = 'SSPRK3'

str	sheath_flux_bc.sorder = 'thirdOrder'

float	sheath_flux_bc.cfl = 0.209

str	sheath_flux_bc.mesh_res = '32x192x32'

f	sheath_flux_bc.meshfile = f'sheath_mesh_{mesh_res}.inp'

warpy	sheath_flux_bc.mesh

float	sheath_flux_bc.density_min = 0.01

float	sheath_flux_bc.pressure_min = 0.01

float	sheath_flux_bc.gamma = 3.0

bool	sheath_flux_bc.use_limiter_mr = False

int	sheath_flux_bc.limiter_alpha_coefficient = 2.

bool	sheath_flux_bc.use_phmaxwell = False

float	sheath_flux_bc.phchi = 1.5

float	sheath_flux_bc.phgamma = phchi

argparse	sheath_flux_bc.parser = argparse.ArgumentParser(description='1d Sheath problem')

	sheath_flux_bc.type

	sheath_flux_bc.str

	sheath_flux_bc.nargs

	sheath_flux_bc.help

	sheath_flux_bc.default

argparse	sheath_flux_bc.args = parser.parse_args()

argparse	sheath_flux_bc.model = args.model

list	sheath_flux_bc.tf_sds = ['middle']

list	sheath_flux_bc.tf_vds = ['left','right']

list	sheath_flux_bc.kinetic_sds = ['left','right']

list	sheath_flux_bc.kinetic_vds = ['middle']

	sheath_flux_bc.else :

float	sheath_flux_bc.qe = 1.602177e-19

float	sheath_flux_bc.mp = 1.672622e-27

float	sheath_flux_bc.mu0 = 1.256637e-6

float	sheath_flux_bc.eps0 = 8.854187e-12

int	sheath_flux_bc.c_light = 1/(mu0eps0)*.5

tuple	sheath_flux_bc.L_debye = (eps0qe1/(1e19qe2))*0.5

tuple	sheath_flux_bc.L = L_debye

float	sheath_flux_bc.ompt = 1.0

float	sheath_flux_bc.omct = 3.265e-5

float	sheath_flux_bc.c_tilde = ompt/omct

int	sheath_flux_bc.skin_depth_norm = 1/omct

float	sheath_flux_bc.omega_p_norm = ompt

float	sheath_flux_bc.n0 = mpomct2/(mu0L*2qe**2)

float	sheath_flux_bc.tau = omptL/(omctc_light)

float	sheath_flux_bc.b0 = mpc_lightomct*2/(qeL*ompt)

float	sheath_flux_bc.e0 = b0**2/mu0

float	sheath_flux_bc.T0 = e0/(n0*qe)

float	sheath_flux_bc.Z_i = +1.0

float	sheath_flux_bc.Z_e = -1.0

float	sheath_flux_bc.mass_ratio = 1836.0

float	sheath_flux_bc.A_i = 1.0

float	sheath_flux_bc.A_e = A_i / mass_ratio

float	sheath_flux_bc.n_i = 1.0

float	sheath_flux_bc.n_e = 1.0

float	sheath_flux_bc.rho_i = A_i * n_i

float	sheath_flux_bc.rho_e = A_e * n_e

float	sheath_flux_bc.charge_to_mass_ion = Z_i/A_i

float	sheath_flux_bc.charge_to_mass_electron = Z_e/A_e

int	sheath_flux_bc.T_i = 10

int	sheath_flux_bc.T_e = 10

float	sheath_flux_bc.P_i = n_i * T_i

float	sheath_flux_bc.P_e = n_e * T_e

math	sheath_flux_bc.v_th_i = math.sqrt(P_i/A_i/n_i)

math	sheath_flux_bc.v_th_e = math.sqrt(P_e/A_e/n_e)

float	sheath_flux_bc.nu_p_tau_bgk = 1.0

float	sheath_flux_bc.vth_factor = +6.0

list	sheath_flux_bc.v_space_num_elements = [48]

list	sheath_flux_bc.min_v_ion = [-vth_factor*v_th_i]

list	sheath_flux_bc.dv_ion = [0.25* v_th_i]

list	sheath_flux_bc.min_v_electron = [-vth_factor*v_th_e]

list	sheath_flux_bc.dv_electron = [0.25* v_th_e]

list	sheath_flux_bc.fluid_comps = ['rho', 'px', 'py', 'pz', 'e']

list	sheath_flux_bc.field_comps = ['Ex', 'Ey', 'Ez', 'Bx', 'By', 'Bz']

warpy	sheath_flux_bc.ions

warpy	sheath_flux_bc.electrons

warpy	sheath_flux_bc.field

warpy	sheath_flux_bc.ion_maxwellian_distribution

warpy	sheath_flux_bc.electron_maxwellian_distribution

warpy	sheath_flux_bc.ion_positive_maxwellian_flux_moments

warpy	sheath_flux_bc.electron_positive_maxwellian_flux_moments

warpy	sheath_flux_bc.ion_negative_maxwellian_flux_moments

warpy	sheath_flux_bc.electron_negative_maxwellian_flux_moments

warpy	sheath_flux_bc.electron_distribution

warpy	sheath_flux_bc.ion_distribution

warpy	sheath_flux_bc.kinetic_electron_fluid

warpy	sheath_flux_bc.kinetic_ion_fluid

warpy	sheath_flux_bc.ion_positive_kinetic_flux_moments

warpy	sheath_flux_bc.electron_positive_kinetic_flux_moments

warpy	sheath_flux_bc.ion_negative_kinetic_flux_moments

warpy	sheath_flux_bc.electron_negative_kinetic_flux_moments

warpy	sheath_flux_bc.chi_metric_i

warpy	sheath_flux_bc.chi_metric_e

tuple	sheath_flux_bc.ic_fun_i

tuple	sheath_flux_bc.ic_fun_e

tuple	sheath_flux_bc.ic_fun_EM

warpy	sheath_flux_bc.electron_ic

warpy	sheath_flux_bc.ion_ic

warpy	sheath_flux_bc.ic_fun_i_va

warpy	sheath_flux_bc.ic_fun_e_va

warpy	sheath_flux_bc.ic_fun_EM_va

warpy	sheath_flux_bc.electron_ic_va

warpy	sheath_flux_bc.ion_ic_va

warpy	sheath_flux_bc.tf_ic_ha

warpy	sheath_flux_bc.field_ic_ha

warpy	sheath_flux_bc.distribution_ic_ha

list	sheath_flux_bc.ic_ha = []

list	sheath_flux_bc.bc_apps = []

list	sheath_flux_bc.values_in = [A_e*1e-4,0,0,0,1e-4/(gamma-1)]

list	sheath_flux_bc.electron_distribution_wall_bcs = []

list	sheath_flux_bc.ion_distribution_wall_bcs = []

list	sheath_flux_bc.maxwell_bcs = []

list	sheath_flux_bc.virtual_bc_apps_kinetic_ions = []

list	sheath_flux_bc.virtual_bc_apps_kinetic_electrons = []

list	sheath_flux_bc.applications_moments_ions = []

list	sheath_flux_bc.applications_moments_electrons = []

list	sheath_flux_bc.applications_2nd_moments_ions = []

list	sheath_flux_bc.applications_2nd_moments_electrons = []

list	sheath_flux_bc.applications_positive_maxwellian_flux_moments_ions = []

list	sheath_flux_bc.applications_negative_maxwellian_flux_moments_ions = []

list	sheath_flux_bc.applications_positive_maxwellian_flux_moments_electrons = []

list	sheath_flux_bc.applications_negative_maxwellian_flux_moments_electrons = []

list	sheath_flux_bc.applications_positive_kinetic_distribution_flux_moments_ions = []

list	sheath_flux_bc.applications_negative_kinetic_distribution_flux_moments_ions = []

list	sheath_flux_bc.applications_positive_kinetic_distribution_flux_moments_electrons = []

list	sheath_flux_bc.applications_negative_kinetic_distribution_flux_moments_electrons = []

list	sheath_flux_bc.applications_maxwellian_metric_ions = []

list	sheath_flux_bc.applications_maxwellian_metric_electrons = []

list	sheath_flux_bc.limiter_apps_electron = []

list	sheath_flux_bc.limiter_apps_ion = []

list	sheath_flux_bc.variable_adjusters_tf = []

list	sheath_flux_bc.variable_adjusters_kinetic = []

list	sheath_flux_bc.kinetic_to_tf_va = []

list	sheath_flux_bc.tf_to_kinetic_va_ions = []

list	sheath_flux_bc.tf_to_kinetic_va_electrons = []

warpy	sheath_flux_bc.ion_timestep_start_moment_taker
	add these same moment takers that can be done at the beginning of each timestep before writeout

warpy	sheath_flux_bc.electron_timestep_start_moment_taker

warpy	sheath_flux_bc.ion_timestep_start_second_moment_taker

warpy	sheath_flux_bc.electron_timestep_start_second_moment_taker

warpy	sheath_flux_bc.ion_timestep_start_chi_metric_moment_taker

warpy	sheath_flux_bc.electron_timestep_start_chi_metric_moment_taker

list	sheath_flux_bc.extra_timestep_start_writer = []

warpy	sheath_flux_bc.moment_takers_timestep_start_ha

warpy	sheath_flux_bc.fluid_ions_to_maxwellian_function

warpy	sheath_flux_bc.fluid_electrons_to_maxwellian_function

list	sheath_flux_bc.maxwellian_evaluation_timestep_start_va = []
	add these maxwellian conversions that can be done at the beginning of each timestep before writeout

warpy	sheath_flux_bc.maxwellian_evaluation_timestep_start_ha

list	sheath_flux_bc.solver_apps_tf = []

list	sheath_flux_bc.solver_apps_maxwell_flux = []

list	sheath_flux_bc.solver_apps_maxwell_tf_fluid_source = []

list	sheath_flux_bc.solver_apps_maxwell_kinetic_fluid_source = []

list	sheath_flux_bc.apps_kinetic_ions = []

list	sheath_flux_bc.apps_kinetic_electrons = []

warpy	sheath_flux_bc.spatial_solver_tf
	Spatial solver.

warpy	sheath_flux_bc.spatial_solver_kinetic_ions

warpy	sheath_flux_bc.spatial_solver_kinetic_electrons

warpy	sheath_flux_bc.spatial_solver_kinetic_field

list	sheath_flux_bc.ss_s = [spatial_solver_tf]

list	sheath_flux_bc.va_s = variable_adjusters_tf

warpy	sheath_flux_bc.time_integrator = warpy.host_actions.erk(name="ti", scheme=torder, spatial_solvers=ss_s, variable_adjusters=va_s)

list	sheath_flux_bc.write_vars = [ions, electrons, field]

warpy	sheath_flux_bc.writer

int	sheath_flux_bc.dt = 1e-4
	Simulation (dg) dt = 1e-6.

tuple	sheath_flux_bc.ompe = (1e19qe2/(eps0A_emp))*0.5

float	sheath_flux_bc.t_final = 20.0/ompe/tau

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

int	sheath_flux_bc.write_steps = 100

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

f	sheath_flux_bc.sim_name = f'sheath_fm_{v_space_num_elements[0]}_{mesh_res}_{model}'

warpy	sheath_flux_bc.config_args = warpy.load_configs()

warpy	sheath_flux_bc.case_loc = config_args['test_rundir'](f'hybrid/sheath/{now}', '')

f	sheath_flux_bc.sim_param

warpy	sheath_flux_bc.sim

	sheath_flux_bc.gen_xdmf
	Run Simulation ...

	sheath_flux_bc.detect_nonscalar

Namespaces

Variables