examples.dg.advection.advection_test Namespace Reference

Variables
str	torder = 'RK1'
str	sorder = 'thirdOrder'
list	velocity1 = [1.0,0.0,0.0]
list	velocity2 = [-1.0,0.0,0.0]
bool	periodic_boundary_conditions = True
bool	use_limiters = True
warpy	mesh = warpy.mesh.block(Bounds=[-0.5, 0.5],NumCells=[128],NodeSets=['Left', 'Right'],NumLayers=2,PeriodicBoundaries=['Left','Right'] if periodic_boundary_conditions else None,basis_array_set=sorder)
warpy	fluid1 = warpy.variable(name='fluid1',components=['rho'],basis_array_set=sorder)
warpy	fluid2 = warpy.variable(name='fluid2',components=['rho'],basis_array_set=sorder)
warpy	fluid_ic1 = warpy.functions.heaviside(name='fluid_ic1',neg=0.5,pos=0,direction=[1,0,0],center=[0,0,0])
warpy	fluid_ic2 = warpy.functions.heaviside(name='fluid_ic2',neg=0,pos=0.5,direction=[1,0,0],center=[0,0,0])
list	applicators = []
list	apps = []
list	bc_apps = []
warpy	writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid1,fluid2])
warpy	spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps,on_subdomains=['all'])
list	variable_adjusters = []
warpy	temporal_solver = warpy.host_actions.erk(name='rk', scheme=torder,spatial_solvers=[spatial_solver],variable_adjusters=variable_adjusters)
int	dt = 1e-2
warpy	dt_controller = warpy.dt_calc.stability_dt(dt)
int	write_steps = 250
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

applicators

list examples.dg.advection.advection_test.applicators = []

apps

list examples.dg.advection.advection_test.apps = []

bc_apps

list examples.dg.advection.advection_test.bc_apps = []

detect_nonscalar

examples.dg.advection.advection_test.detect_nonscalar

dt

int examples.dg.advection.advection_test.dt = 1e-2

dt_controller

warpy examples.dg.advection.advection_test.dt_controller = warpy.dt_calc.stability_dt(dt)

fluid1

warpy examples.dg.advection.advection_test.fluid1 = warpy.variable(name='fluid1',components=['rho'],basis_array_set=sorder)

fluid2

warpy examples.dg.advection.advection_test.fluid2 = warpy.variable(name='fluid2',components=['rho'],basis_array_set=sorder)

fluid_ic1

warpy examples.dg.advection.advection_test.fluid_ic1 = warpy.functions.heaviside(name='fluid_ic1',neg=0.5,pos=0,direction=[1,0,0],center=[0,0,0])

fluid_ic2

warpy examples.dg.advection.advection_test.fluid_ic2 = warpy.functions.heaviside(name='fluid_ic2',neg=0,pos=0.5,direction=[1,0,0],center=[0,0,0])

gen_xdmf

examples.dg.advection.advection_test.gen_xdmf

mesh

warpy examples.dg.advection.advection_test.mesh = warpy.mesh.block(Bounds=[-0.5, 0.5],NumCells=[128],NodeSets=['Left', 'Right'],NumLayers=2,PeriodicBoundaries=['Left','Right'] if periodic_boundary_conditions else None,basis_array_set=sorder)

periodic_boundary_conditions

bool examples.dg.advection.advection_test.periodic_boundary_conditions = True

sim

warpy examples.dg.advection.advection_test.sim

Initial value:

=  warpy.dg_sim(name='advection_heaviside_test',
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[temporal_solver],
                   writers=[writer],
                   time=[0,dt*write_steps],
                   dt_controller=dt_controller,
                   flexible_writeout=False,
                   write_steps=write_steps,
                   verbosity='info')

sorder

str examples.dg.advection.advection_test.sorder = 'thirdOrder'

spatial_solver

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

temporal_solver

warpy examples.dg.advection.advection_test.temporal_solver = warpy.host_actions.erk(name='rk', scheme=torder,spatial_solvers=[spatial_solver],variable_adjusters=variable_adjusters)

torder

str examples.dg.advection.advection_test.torder = 'RK1'

True

examples.dg.advection.advection_test.True

use_limiters

bool examples.dg.advection.advection_test.use_limiters = True

variable_adjusters

list examples.dg.advection.advection_test.variable_adjusters = []

velocity1

list examples.dg.advection.advection_test.velocity1 = [1.0,0.0,0.0]

velocity2

list examples.dg.advection.advection_test.velocity2 = [-1.0,0.0,0.0]

write_steps

int examples.dg.advection.advection_test.write_steps = 250

writer

warpy examples.dg.advection.advection_test.writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid1,fluid2])