dory_guest_harris_instability.py File Reference

Namespaces
namespace	dory_guest_harris_instability

Variables
str	dory_guest_harris_instability.torder = 'RK4'
str	dory_guest_harris_instability.sorder = 'thirdOrder'
bool	dory_guest_harris_instability.periodic_boundary_conditions = True
bool	dory_guest_harris_instability.use_phmaxwell = False
float	dory_guest_harris_instability.cfl = 0.235
float	dory_guest_harris_instability.qe = 1.602177e-19
float	dory_guest_harris_instability.me = 9.10938356E-31
float	dory_guest_harris_instability.mu0 = 1.256637e-6
float	dory_guest_harris_instability.eps0 = 8.854187e-12
int	dory_guest_harris_instability.c_light = 1/(mu0*eps0)**.5
float	dory_guest_harris_instability.Z_i = +1.0
float	dory_guest_harris_instability.Z_e = -1.0
float	dory_guest_harris_instability.mass_ratio = 1836.0
float	dory_guest_harris_instability.A_i = 1836.0
float	dory_guest_harris_instability.A_e = A_i / mass_ratio
float	dory_guest_harris_instability.n_i = 1.0
float	dory_guest_harris_instability.n_e = 1.0
float	dory_guest_harris_instability.rho_i = A_i * n_i
float	dory_guest_harris_instability.rho_e = A_e * n_e
float	dory_guest_harris_instability.charge_to_mass_ion = Z_i/A_i
float	dory_guest_harris_instability.charge_to_mass_electron = Z_e/A_e
float	dory_guest_harris_instability.v_th_i = 1.0
float	dory_guest_harris_instability.v_th_e = 1.0
float	dory_guest_harris_instability.P_i = n_i * A_i * v_th_i**2
float	dory_guest_harris_instability.P_e = n_e * A_e * v_th_e**2
float	dory_guest_harris_instability.T_i = P_i / n_i
float	dory_guest_harris_instability.T_e = P_e / n_e
int	dory_guest_harris_instability.j_value = 6
int	dory_guest_harris_instability.omega_p_over_omega_c = 20
float	dory_guest_harris_instability.k_tilde = 3.15
float	dory_guest_harris_instability.epsilon_pert = 1.0E-4
math	dory_guest_harris_instability.v_perp0 = math.sqrt(2)
int	dory_guest_harris_instability.Omega_c = 1 / omega_p_over_omega_c
float	dory_guest_harris_instability.Bz_ic = 1.0
float	dory_guest_harris_instability.L = 1.0
float	dory_guest_harris_instability.omega_p_tau = 1.0
int	dory_guest_harris_instability.skin_depth_norm = omega_p_over_omega_c
int	dory_guest_harris_instability.omega_c_tau = 1/skin_depth_norm
float	dory_guest_harris_instability.n0 = me*omega_c_tau**2/(mu0*L**2*qe**2)
float	dory_guest_harris_instability.c_tilde = omega_p_tau * skin_depth_norm
int	dory_guest_harris_instability.v0 = c_light / (skin_depth_norm * omega_p_tau)
float	dory_guest_harris_instability.tau = L/v0
float	dory_guest_harris_instability.omega_p = omega_p_tau / tau
int	dory_guest_harris_instability.omega_c = omega_c_tau / tau
int	dory_guest_harris_instability.b0 = v0 * math.sqrt(me*n0*mu0)
int	dory_guest_harris_instability.p0 = b0**2/mu0
int	dory_guest_harris_instability.T0 = p0/n0
math	dory_guest_harris_instability.v_th0 = math.sqrt(T0/me)
math	dory_guest_harris_instability.rl0 = v_th0/omega_p
int	dory_guest_harris_instability.dp0 = c_light / omega_p
tuple	dory_guest_harris_instability.L_debye0 = (eps0*T0/(n0*qe**2))**0.5
tuple	dory_guest_harris_instability.L_debye = (eps0*T_e*T0/(n0*n_e*qe**2))**0.5
int	dory_guest_harris_instability.normalized_domain_length = 2*math.pi/k_tilde * v_perp0 / Omega_c
int	dory_guest_harris_instability.v_bound = 4
list	dory_guest_harris_instability.v_space_num_elements = [24,24]
int	dory_guest_harris_instability.dv = 2*v_bound / v_space_num_elements[0]
list	dory_guest_harris_instability.min_v_ion = [-v_bound, -v_bound]
list	dory_guest_harris_instability.dv_ion = [dv, dv]
list	dory_guest_harris_instability.min_v_electron = [-v_bound, -v_bound]
list	dory_guest_harris_instability.dv_electron = [dv,dv]
int	dory_guest_harris_instability.num_physical_elements = 48
list	dory_guest_harris_instability.domain = ['domain']
warpy	dory_guest_harris_instability.mesh
warpy	dory_guest_harris_instability.electron_distribution
warpy	dory_guest_harris_instability.ion_distribution
warpy	dory_guest_harris_instability.kinetic_electron_fluid
warpy	dory_guest_harris_instability.kinetic_ion_fluid
list	dory_guest_harris_instability.field_components = ['Ex','Ey','Ez','Bx','By','Bz']
warpy	dory_guest_harris_instability.field
warpy	dory_guest_harris_instability.electron_ic
	Initial conditions-------------------------------— Step 1: define the function as an application:
warpy	dory_guest_harris_instability.ion_ic
warpy	dory_guest_harris_instability.field_ic
warpy	dory_guest_harris_instability.electron_ic_va
warpy	dory_guest_harris_instability.ion_ic_va
warpy	dory_guest_harris_instability.field_ic_va
warpy	dory_guest_harris_instability.ic_ha
list	dory_guest_harris_instability.apps_kinetic_electrons = []
	Initial conditions-------------------------------—.
list	dory_guest_harris_instability.apps_physical = []
	dory_guest_harris_instability.else :
list	dory_guest_harris_instability.applications_moments_electrons = []
list	dory_guest_harris_instability.applications_2nd_moments_electrons = []
list	dory_guest_harris_instability.variable_adjusters = []
	Variable adjusters.
warpy	dory_guest_harris_instability.writer
	Variable writer.
warpy	dory_guest_harris_instability.spatial_solver_kinetic_electrons
warpy	dory_guest_harris_instability.spatial_solver_physical
warpy	dory_guest_harris_instability.temporal_solver
	Temporal solver.
tuple	dory_guest_harris_instability.ompe = (( n_e*n0 * qe**2)/(eps0*A_e*me))**0.5
	Simulation (dg)
int	dory_guest_harris_instability.t_final = 500/ompe/tau
int	dory_guest_harris_instability.write_steps = 100
int	dory_guest_harris_instability.dt = t_final / write_steps
warpy	dory_guest_harris_instability.dt_controller = warpy.dt_calc.stability_dt(dt)