WARPXM v1.10.0
Loading...
Searching...
No Matches
examples.dg.mhd.gem_challenge_imhd2 Namespace Reference

Variables

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

Variable Documentation

◆ A_i

int examples.dg.mhd.gem_challenge_imhd2.A_i = 1

◆ always_limit

bool examples.dg.mhd.gem_challenge_imhd2.always_limit = True

◆ applicators

list examples.dg.mhd.gem_challenge_imhd2.applicators = [warpy.applicator(spatial_order=sorder,fun=gem_ic,var=fluid,spatial_scheme='Nodal')]

◆ apps

list examples.dg.mhd.gem_challenge_imhd2.apps = []

◆ B0

float examples.dg.mhd.gem_challenge_imhd2.B0 = 1.0

◆ component_names

list examples.dg.mhd.gem_challenge_imhd2.component_names = ['rho', 'px', 'py', 'pz', 'e', 'Bx', 'By', 'Bz']

◆ detect_nonscalar

examples.dg.mhd.gem_challenge_imhd2.detect_nonscalar

◆ di_L

int examples.dg.mhd.gem_challenge_imhd2.di_L = 1.

◆ dirichlet_apps

list examples.dg.mhd.gem_challenge_imhd2.dirichlet_apps = []

◆ dt_controller

warpy examples.dg.mhd.gem_challenge_imhd2.dt_controller = warpy.dt_calc.stability_dt(init_dt=1e-5)

◆ fluid

warpy examples.dg.mhd.gem_challenge_imhd2.fluid = warpy.variable(name='fluid',components=component_names,basis_array_set=sorder)

◆ fluid_gradient

warpy examples.dg.mhd.gem_challenge_imhd2.fluid_gradient = warpy.variable(name='fluid_gradient', components=gradient_component_names,basis_array_set=sorder)

◆ gamma

int examples.dg.mhd.gem_challenge_imhd2.gamma = 5./3.

◆ gem_ic

warpy examples.dg.mhd.gem_challenge_imhd2.gem_ic = warpy.functions.gem_hallmhd(name='gem_ic',gamma=gamma,mass=A_i,temperature=Ti+Te,sheet_width=sheet_width,number_density=n0,number_density_infinity=n_inf,Bx=B0,Lx=Lx,Ly=Ly)

◆ gen_xdmf

examples.dg.mhd.gem_challenge_imhd2.gen_xdmf

◆ gradient_component_names

list examples.dg.mhd.gem_challenge_imhd2.gradient_component_names = [prefix+"_"+suffix for prefix in component_names for suffix in ['x','y','z']]

◆ Lx

int examples.dg.mhd.gem_challenge_imhd2.Lx = 8*math.pi*di_L

◆ Ly

int examples.dg.mhd.gem_challenge_imhd2.Ly = 4*math.pi*di_L

◆ mesh

warpy examples.dg.mhd.gem_challenge_imhd2.mesh = warpy.mesh.block(Bounds=[-Lx/2.,Lx/2., -Ly/2.,Ly/2.],NumCells=[64,32],NodeSets=['Left','Right', 'Lower','Upper'],NumLayers=1,basis_array_set=sorder,PeriodicBoundaries=['Left','Right'],block_type="symmetric")

◆ n0

float examples.dg.mhd.gem_challenge_imhd2.n0 = 1.0

◆ n_inf

float examples.dg.mhd.gem_challenge_imhd2.n_inf = 0.2

◆ neumann_apps

list examples.dg.mhd.gem_challenge_imhd2.neumann_apps = []

◆ sheet_width

float examples.dg.mhd.gem_challenge_imhd2.sheet_width = 0.5

◆ sim

warpy examples.dg.mhd.gem_challenge_imhd2.sim = warpy.dg_sim(name='gem_imhd_may28_testing',meshes=[mesh],initial_conditions=applicators,temporal_solvers=[ti],writers=[writer],time=[0,4.0],dt_controller=dt_controller,write_steps=100,verbosity='debug')

◆ skin_depth_norm

float examples.dg.mhd.gem_challenge_imhd2.skin_depth_norm = 1.0

◆ sorder

str examples.dg.mhd.gem_challenge_imhd2.sorder = 'secondOrder'

◆ spatial_solver

warpy examples.dg.mhd.gem_challenge_imhd2.spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps, on_subdomains=['all'])

◆ T

list examples.dg.mhd.gem_challenge_imhd2.T = [Ti, Te]

◆ Te

float examples.dg.mhd.gem_challenge_imhd2.Te = B0**2/2.0/n0/(Ti_Te + 1.0)

◆ theta

int examples.dg.mhd.gem_challenge_imhd2.theta = 1

◆ Ti

float examples.dg.mhd.gem_challenge_imhd2.Ti = Ti_Te * Te

◆ ti

warpy examples.dg.mhd.gem_challenge_imhd2.ti = warpy.host_actions.erk(name="rk", scheme=torder, spatial_solvers=[spatial_solver], variable_adjusters=variable_adjusters)

◆ Ti_Te

float examples.dg.mhd.gem_challenge_imhd2.Ti_Te = 1.0/theta

◆ torder

str examples.dg.mhd.gem_challenge_imhd2.torder = 'RK4'

◆ True

examples.dg.mhd.gem_challenge_imhd2.True

◆ useLimiter

bool examples.dg.mhd.gem_challenge_imhd2.useLimiter = True

◆ v

int examples.dg.mhd.gem_challenge_imhd2.v = 2. * skin_depth_norm / (sheet_width * B0)

◆ V

list examples.dg.mhd.gem_challenge_imhd2.V = T[0]*v

◆ variable_adjusters

list examples.dg.mhd.gem_challenge_imhd2.variable_adjusters = []

◆ writer

warpy examples.dg.mhd.gem_challenge_imhd2.writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid])

◆ z_i

int examples.dg.mhd.gem_challenge_imhd2.z_i = 1