WARPXM v1.10.0
Loading...
Searching...
No Matches
Variables
examples.dg.5-moment.flow_ei 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 = 1
 
int Ae = 1/10
 
int An = 1 + Ae
 
float phi_ion_eV = 13.6
 
float phi_ion = phi_ion_eV/T0
 
float lx = 1.0
 
int nx = 100
 
int ny = 100
 
float dx = lx/nx
 
warpy mesh
 
warpy ions
 
warpy electrons
 
warpy field
 
tuple ic_fun_i
 
tuple ic_fun_e
 
tuple ic_fun_EM
 
list applicators = []
 
list solver_apps = []
 
list bc_apps = []
 
float rhoin = 0.1
 
float pxin = 0.0
 
float pyin = 0.0
 
float pzin = 0.0
 
float pressurein = 0.1
 
float energyin = pressurein/(gamma-1)
 
int inject_t = 0
 
float inject_rho = 1.0
 
float inject_momy = 1.0
 
warpy spatial_solver
 
list var_adj = []
 
warpy writer
 
warpy time_integrator
 
float t_tot = 0.2
 
int dt = 2e-7
 
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.flow_ei.Ae = 1/10

◆ Ai

int examples.dg.5-moment.flow_ei.Ai = 1

◆ An

int examples.dg.5-moment.flow_ei.An = 1 + Ae

◆ applicators

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

◆ atom

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

◆ b0

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

◆ bc_apps

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

◆ c_light

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

◆ c_tilde

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

◆ cfl

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

◆ density_min

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

◆ detect_nonscalar

examples.dg.5-moment.flow_ei.detect_nonscalar

◆ dt

int examples.dg.5-moment.flow_ei.dt = 2e-7

◆ dt_controller

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

◆ dx

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

◆ e0

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

◆ electrons

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

◆ energyin

float examples.dg.5-moment.flow_ei.energyin = pressurein/(gamma-1)

◆ eps0

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

◆ field

warpy examples.dg.5-moment.flow_ei.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.flow_ei.gamma = 5/3

◆ gen_xdmf

examples.dg.5-moment.flow_ei.gen_xdmf

◆ ic_fun_e

tuple examples.dg.5-moment.flow_ei.ic_fun_e
Initial value:
1= (
2 warpy.functions.flow_init2_5moment(name='ic_flow_e',gamma=gamma,mass=Ae))
warpy.functions.flow_init2.flow_init2_5moment
Function which sets IC for simple 5-moment fluid flow problem.
Definition: flow_init2.py:4

◆ ic_fun_EM

tuple examples.dg.5-moment.flow_ei.ic_fun_EM
Initial value:
1= (
2 warpy.functions.flow_init_EM_5moment(name='ic_flow_EM'))
warpy.functions.flow_init_EM.flow_init_EM_5moment
Sets IC for EM fields in a simple 5-moment fluid flow problem.
Definition: flow_init_EM.py:4

◆ ic_fun_i

tuple examples.dg.5-moment.flow_ei.ic_fun_i
Initial value:
1= (
2 warpy.functions.flow_init2_5moment(name='ic_flow_i',gamma=gamma,mass=Ai))

◆ inject_momy

float examples.dg.5-moment.flow_ei.inject_momy = 1.0

◆ inject_rho

float examples.dg.5-moment.flow_ei.inject_rho = 1.0

◆ inject_t

int examples.dg.5-moment.flow_ei.inject_t = 0

◆ ions

warpy examples.dg.5-moment.flow_ei.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.flow_ei.L = 1

◆ lx

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

◆ mesh

warpy examples.dg.5-moment.flow_ei.mesh
Initial value:
1= warpy.mesh.block(Bounds=[0.0, lx, -0.5, 0.5],
2 # NumCells=[300,60],
3 # NumCells=[150,30],
4 # NumCells=[30,6],
5 NumCells=[nx,ny],
6 # NumCells=[2000,400],
7 # NodeSets=['leftWall','rightWall'],
8 NodeSets=['leftWall','rightWall',
9 'bottomWall','topWall'],
10 NumLayers=1,
11 PeriodicBoundaries = ['bottomWall', 'topWall'],
12 # PeriodicBoundaries = ['leftWall', 'rightWall'],
13 # PeriodicBoundaries = ['leftWall','rightWall',
14 # 'bottomWall', 'topWall'],
15 basis_array_set=sorder,
16 block_type="symmetric")
warpy.mesh.block
block mesh generator
Definition: mesh.py:87

◆ mp

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

◆ mu0

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

◆ n0

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

◆ nout

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

◆ nx

int examples.dg.5-moment.flow_ei.nx = 100

◆ ny

int examples.dg.5-moment.flow_ei.ny = 100

◆ omct

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

◆ ompt

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

◆ phi_ion

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

◆ phi_ion_eV

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

◆ pressure_min

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

◆ pressurein

float examples.dg.5-moment.flow_ei.pressurein = 0.1

◆ pxin

float examples.dg.5-moment.flow_ei.pxin = 0.0

◆ pyin

float examples.dg.5-moment.flow_ei.pyin = 0.0

◆ pzin

float examples.dg.5-moment.flow_ei.pzin = 0.0

◆ qe

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

◆ rhoin

float examples.dg.5-moment.flow_ei.rhoin = 0.1

◆ sim

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

◆ skin_depth_norm

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

◆ solver_apps

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

◆ sorder

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

◆ spatial_solver

warpy examples.dg.5-moment.flow_ei.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.flow_ei.T0 = e0/(n0*qe)

◆ t_tot

float examples.dg.5-moment.flow_ei.t_tot = 0.2

◆ tau

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

◆ time_integrator

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

◆ torder

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

◆ True

examples.dg.5-moment.flow_ei.True

◆ var_adj

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

◆ writer

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

◆ Ze

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

◆ Zi

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

◆ Zn

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