WARPXM v1.10.0
Loading...
Searching...
No Matches
gem_challenge_imhd.py File Reference

Namespaces

namespace  examples
 
namespace  examples.dg
 
namespace  examples.dg.mhd
 
namespace  examples.dg.mhd.gem_challenge_imhd
 

Variables

str examples.dg.mhd.gem_challenge_imhd.torder = 'RK4'
 
str examples.dg.mhd.gem_challenge_imhd.sorder = 'thirdOrder'
 
bool examples.dg.mhd.gem_challenge_imhd.useLimiter = False
 
bool examples.dg.mhd.gem_challenge_imhd.always_limit = True
 
float examples.dg.mhd.gem_challenge_imhd.skin_depth_norm = 1.0
 
int examples.dg.mhd.gem_challenge_imhd.gamma = 5./3.
 
int examples.dg.mhd.gem_challenge_imhd.z_i = 1
 
int examples.dg.mhd.gem_challenge_imhd.A_i = 1
 
int examples.dg.mhd.gem_challenge_imhd.theta = 1
 
int examples.dg.mhd.gem_challenge_imhd.di_L = 1.
 
int examples.dg.mhd.gem_challenge_imhd.Lx = 8*math.pi*di_L
 
int examples.dg.mhd.gem_challenge_imhd.Ly = 4*math.pi*di_L
 
float examples.dg.mhd.gem_challenge_imhd.sheet_width = 0.5
 
float examples.dg.mhd.gem_challenge_imhd.n0 = 1.0
 
float examples.dg.mhd.gem_challenge_imhd.n_inf = 0.2
 
float examples.dg.mhd.gem_challenge_imhd.B0 = 1.0
 
float examples.dg.mhd.gem_challenge_imhd.Ti_Te = 1.0/theta
 
tuple examples.dg.mhd.gem_challenge_imhd.Te = (di_L / skin_depth_norm * B0)**2 / (2 * (n0-n_inf) * (1 + Ti_Te))
 
list examples.dg.mhd.gem_challenge_imhd.T = [Ti_Te * Te, Te]
 
int examples.dg.mhd.gem_challenge_imhd.v = 2. * skin_depth_norm / (sheet_width * B0)
 
list examples.dg.mhd.gem_challenge_imhd.V = T[0]*v
 
warpy examples.dg.mhd.gem_challenge_imhd.mesh = warpy.mesh.block(Bounds=[-Lx/2.,Lx/2., -Ly/2.,Ly/2.],NumCells=[4,2],NodeSets=['Left','Right', 'Lower','Upper'],NumLayers=2,basis_array_set=sorder,PeriodicBoundaries=['Left','Right'])
 
list examples.dg.mhd.gem_challenge_imhd.component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']
 
warpy examples.dg.mhd.gem_challenge_imhd.fluid = warpy.variable(name='fluid',components=component_names,basis_array_set=sorder)
 
list examples.dg.mhd.gem_challenge_imhd.gradient_component_names = [prefix+"_"+suffix for prefix in component_names for suffix in ['x','y','z']]
 
warpy examples.dg.mhd.gem_challenge_imhd.fluid_gradient = warpy.variable(name='fluid_gradient', components=gradient_component_names,basis_array_set=sorder)
 
warpy examples.dg.mhd.gem_challenge_imhd.gem_ic = warpy.functions.gem_hallmhd(name='gem_ic',gamma=gamma,mass=A_i,temperature=T[0],velocity=V,sheet_width=sheet_width,number_density=n0,number_density_infinity=n_inf,Bx=B0,Lx=Lx,Ly=Ly)
 
list examples.dg.mhd.gem_challenge_imhd.applicators = [warpy.applicator(spatial_order=sorder,fun=gem_ic,var=fluid,spatial_scheme='Nodal')]
 
list examples.dg.mhd.gem_challenge_imhd.apps = []
 
list examples.dg.mhd.gem_challenge_imhd.dirichlet_apps = []
 
list examples.dg.mhd.gem_challenge_imhd.neumann_apps = []
 
warpy examples.dg.mhd.gem_challenge_imhd.writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid])
 
warpy examples.dg.mhd.gem_challenge_imhd.spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps, on_subdomains=['all'])
 
list examples.dg.mhd.gem_challenge_imhd.variable_adjusters = []
 
warpy examples.dg.mhd.gem_challenge_imhd.ti = warpy.host_actions.erk(name="rk", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)
 
warpy examples.dg.mhd.gem_challenge_imhd.dt_controller = warpy.dt_calc.stability_dt(init_dt=1e-5)
 
warpy examples.dg.mhd.gem_challenge_imhd.sim = warpy.dg_sim(name='gem_imhd',meshes=[mesh],initial_conditions=applicators,temporal_solvers=[ti],writers=[writer],time=[0,40.0],dt_controller=dt_controller,write_steps=1000,verbosity='info', profile='icpus')
 
 examples.dg.mhd.gem_challenge_imhd.gen_xdmf
 
 examples.dg.mhd.gem_challenge_imhd.True
 
 examples.dg.mhd.gem_challenge_imhd.detect_nonscalar