WARPXM v1.10.0
Loading...
Searching...
No Matches
examples.dg.5-moment.ein_bm Namespace Reference

Variables

str torder = 'RK1'
 
str sorder = 'firstOrder'
 
float cfl = 0.5
 
float density_min = 0.01
 
float pressure_min = 0.01
 
float qe = 1.602177e-19
 
float mp = 1.672622e-27
 
float mu0 = 1.256637e-6
 
float eps0 = 8.854187e-8
 
int c_light = 1/(mu0*eps0)**.5
 
int L = 1
 
float ompt = 4.251e3
 
float omct = 13.89
 
float c_tilde = ompt/omct
 
int skin_depth_norm = 1/omct
 
float n0 = mp*omct**2/(mu0*L**2*qe**2)
 
float tau = ompt*L/(omct*c_light)
 
float b0 = mp*c_light*omct**2/(qe*L*ompt)
 
float e0 = b0**2/mu0
 
float T0 = e0/(n0*qe)
 
int atom = 1
 
int gamma = 5/3
 
int Zi = 1
 
int Ze = -1
 
int Zn = 0
 
int Ai = 2
 
int Ae = 1/10
 
int An = Ai
 
float phi_ion_eV = 13.6
 
float phi_ion = phi_ion_eV/T0
 
float phi_eff_eV = 33.0
 
float phi_eff = phi_eff_eV/T0
 
float lx = 1.0
 
int nx = 300
 
int ny = nx
 
float dx = lx/nx
 
warpy mesh
 
warpy ions
 
warpy electrons
 
warpy neutrals
 
warpy field
 
tuple ic_fun_i
 
tuple ic_fun_e
 
tuple ic_fun_n
 
tuple ic_fun_EM
 
list applicators = []
 
list solver_apps = []
 
list bc_apps = []
 
warpy spatial_solver
 
list var_adj = []
 
warpy writer
 
warpy time_integrator
 
float t_tot = 0.03
 
int dt = 1e-6
 
int nout = 50
 
warpy dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
 
warpy sim
 
 gen_xdmf
 
 True
 
 detect_nonscalar
 

Variable Documentation

◆ Ae

int examples.dg.5-moment.ein_bm.Ae = 1/10

◆ Ai

int examples.dg.5-moment.ein_bm.Ai = 2

◆ An

int examples.dg.5-moment.ein_bm.An = Ai

◆ applicators

list examples.dg.5-moment.ein_bm.applicators = []

◆ atom

int examples.dg.5-moment.ein_bm.atom = 1

◆ b0

float examples.dg.5-moment.ein_bm.b0 = mp*c_light*omct**2/(qe*L*ompt)

◆ bc_apps

list examples.dg.5-moment.ein_bm.bc_apps = []

◆ c_light

int examples.dg.5-moment.ein_bm.c_light = 1/(mu0*eps0)**.5

◆ c_tilde

float examples.dg.5-moment.ein_bm.c_tilde = ompt/omct

◆ cfl

float examples.dg.5-moment.ein_bm.cfl = 0.5

◆ density_min

float examples.dg.5-moment.ein_bm.density_min = 0.01

◆ detect_nonscalar

examples.dg.5-moment.ein_bm.detect_nonscalar

◆ dt

int examples.dg.5-moment.ein_bm.dt = 1e-6

◆ dt_controller

warpy examples.dg.5-moment.ein_bm.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

◆ dx

float examples.dg.5-moment.ein_bm.dx = lx/nx

◆ e0

float examples.dg.5-moment.ein_bm.e0 = b0**2/mu0

◆ electrons

warpy examples.dg.5-moment.ein_bm.electrons
Initial value:
1= warpy.variable(name='electrons',
2 components=['rho','px','py','pz','e'],
3 basis_array_set=sorder)
Definition: variable.py:4

◆ eps0

float examples.dg.5-moment.ein_bm.eps0 = 8.854187e-8

◆ field

warpy examples.dg.5-moment.ein_bm.field
Initial value:
1= warpy.variable(name='field',
2 components=['Ex','Ey','Ez','Bx','By','Bz'],
3 basis_array_set=sorder)

◆ gamma

int examples.dg.5-moment.ein_bm.gamma = 5/3

◆ gen_xdmf

examples.dg.5-moment.ein_bm.gen_xdmf

◆ ic_fun_e

tuple examples.dg.5-moment.ein_bm.ic_fun_e
Initial value:
1= (
3 gamma=gamma,
4 ic_type=1,
5 var_type=1,
6 mass=Ae))
Function which sets IC for 5-moment electron-ion-neutral problem.
Definition: ein_ic.py:4

◆ ic_fun_EM

tuple examples.dg.5-moment.ein_bm.ic_fun_EM
Initial value:
1= (
2 warpy.functions.ein_ic_EM_5moment(name='ic_EM',ic_type=1))
Function which sets IC for EM fields in 5-moment electron-ion-neutral problem.
Definition: ein_ic_EM.py:4

◆ ic_fun_i

tuple examples.dg.5-moment.ein_bm.ic_fun_i
Initial value:
1= (
3 gamma=gamma,
4 ic_type=1,
5 var_type=2,
6 mass=Ai))

◆ ic_fun_n

tuple examples.dg.5-moment.ein_bm.ic_fun_n
Initial value:
1= (
3 gamma=gamma,
4 ic_type=1,
5 var_type=3,
6 mass=An))

◆ ions

warpy examples.dg.5-moment.ein_bm.ions
Initial value:
1= warpy.variable(name='ions',
2 components=['rho','px','py','pz','e'],
3 basis_array_set=sorder)

◆ L

int examples.dg.5-moment.ein_bm.L = 1

◆ lx

float examples.dg.5-moment.ein_bm.lx = 1.0

◆ mesh

warpy examples.dg.5-moment.ein_bm.mesh
Initial value:
1= warpy.mesh.block(Bounds=[0.0, lx, -0.5, 0.5],
2 NumCells=[nx,ny],
3 NodeSets=['leftWall','rightWall',
4 'bottomWall','topWall'],
5 NumLayers=1,
6 PeriodicBoundaries = ['bottomWall', 'topWall'],
7 basis_array_set=sorder,
8 block_type="symmetric")
block mesh generator
Definition: mesh.py:87

◆ mp

float examples.dg.5-moment.ein_bm.mp = 1.672622e-27

◆ mu0

float examples.dg.5-moment.ein_bm.mu0 = 1.256637e-6

◆ n0

float examples.dg.5-moment.ein_bm.n0 = mp*omct**2/(mu0*L**2*qe**2)

◆ neutrals

warpy examples.dg.5-moment.ein_bm.neutrals
Initial value:
1= warpy.variable(name='neutrals',
2 components=['rho','px','py','pz','e'],
3 basis_array_set=sorder)

◆ nout

int examples.dg.5-moment.ein_bm.nout = 50

◆ nx

int examples.dg.5-moment.ein_bm.nx = 300

◆ ny

int examples.dg.5-moment.ein_bm.ny = nx

◆ omct

float examples.dg.5-moment.ein_bm.omct = 13.89

◆ ompt

float examples.dg.5-moment.ein_bm.ompt = 4.251e3

◆ phi_eff

float examples.dg.5-moment.ein_bm.phi_eff = phi_eff_eV/T0

◆ phi_eff_eV

float examples.dg.5-moment.ein_bm.phi_eff_eV = 33.0

◆ phi_ion

float examples.dg.5-moment.ein_bm.phi_ion = phi_ion_eV/T0

◆ phi_ion_eV

float examples.dg.5-moment.ein_bm.phi_ion_eV = 13.6

◆ pressure_min

float examples.dg.5-moment.ein_bm.pressure_min = 0.01

◆ qe

float examples.dg.5-moment.ein_bm.qe = 1.602177e-19

◆ sim

warpy examples.dg.5-moment.ein_bm.sim
Initial value:
1= warpy.dg_sim(name='bm02a_1proc',
2 meshes=[mesh],
3 initial_conditions=applicators,
4 temporal_solvers=[time_integrator],
5 writers=[writer],
6 time=[0,t_tot],
7 dt_controller = dt_controller,
8 flexible_writeout=True,
9 write_steps=nout,
10# verbosity='debug',
11 verbosity='info',
12 enable_fpe=True,
13 num_partitions=1)
Discontinuous finite element RK simulation.
Definition: dg_sim.py:11

◆ skin_depth_norm

int examples.dg.5-moment.ein_bm.skin_depth_norm = 1/omct

◆ solver_apps

list examples.dg.5-moment.ein_bm.solver_apps = []

◆ sorder

str examples.dg.5-moment.ein_bm.sorder = 'firstOrder'

◆ spatial_solver

warpy examples.dg.5-moment.ein_bm.spatial_solver
Initial value:
1= warpy.spatial_solvers.dg(name="dg",
2 spatial_order=sorder,
3 applications=solver_apps,
4 cfl = cfl)

◆ T0

float examples.dg.5-moment.ein_bm.T0 = e0/(n0*qe)

◆ t_tot

float examples.dg.5-moment.ein_bm.t_tot = 0.03

◆ tau

float examples.dg.5-moment.ein_bm.tau = ompt*L/(omct*c_light)

◆ time_integrator

warpy examples.dg.5-moment.ein_bm.time_integrator
Initial value:
1= warpy.host_actions.erk(name="ti", scheme=torder,
2 spatial_solvers=[spatial_solver],
3 variable_adjusters=var_adj)
Explicit Runge-Kutta temporal solver Note: Dormand45 currently will not work correctly with limiters ...
Definition: erk.py:5

◆ torder

str examples.dg.5-moment.ein_bm.torder = 'RK1'

◆ True

examples.dg.5-moment.ein_bm.True

◆ var_adj

list examples.dg.5-moment.ein_bm.var_adj = []

◆ writer

warpy examples.dg.5-moment.ein_bm.writer
Initial value:
1= warpy.host_actions.writer(name='writer',
2 ReadVars=[ions, electrons, neutrals, field])
Writes out a list of variables.
Definition: writer.py:4

◆ Ze

int examples.dg.5-moment.ein_bm.Ze = -1

◆ Zi

int examples.dg.5-moment.ein_bm.Zi = 1

◆ Zn

int examples.dg.5-moment.ein_bm.Zn = 0