examples.dg.diffusion.diffusion2d Namespace Reference

Variables
str	torder = 'SSPRK3'
int	sorder = 3
float	cfl_factor = 0.8
float	cfl = cfl_factor/(2*sorder - 1)
float	diffusivity = 1.0
list	center_bump = [0, 0, 0]
pi	width_bump = pi
int	min_bump = 0
int	max_bump = 1
float	t_tot = 1.0
int	nframe = 10
int	dt = 1e-4
str	dt_calc = "stability"
list	sorders = ['firstOrder','secondOrder','thirdOrder','fourthOrder']
list	sorder_str = sorders[sorder - 1]
pi	xmax = pi
pi	xmin = -pi
pi	ymax = pi
pi	ymin = -pi
warpy	mesh
str	gradient_method = "ldg"
None	penalty_beta = None
float	penalty_eta = 1.0
warpy	fluid = warpy.variable(name='fluid',components=['rho'],basis_array_set=sorder_str)
warpy	fluid_gradient = warpy.variable(name='fluid_gradient',components=['rho_x','rho_y','rho_z'],basis_array_set=sorder_str,explodable=False)
warpy	fluid_gradient2 = warpy.variable(name='fluid_gradient2',components=['rho_x','rho_y','rho_z'],basis_array_set=sorder_str,explodable=False)
warpy	ic_fluid
warpy	ic_va
warpy	ic_ha
list	apps
list	bc_apps = []
list	grad_bc_apps = []
list	variable_adjusters = []
warpy	writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid, fluid_gradient, fluid_gradient2])
warpy	dg
warpy	temporal_solver
warpy	dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)
	else :
warpy	sim
	local_path
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

apps

list examples.dg.diffusion.diffusion2d.apps

Initial value:

=  [warpy.apps.diffusion(name='diffusion',
                             variable=fluid,
                             gradient_method=gradient_method,
                             gradient_variable=fluid_gradient,
                             gradient_variable_interior_penalty_internal_approximation=fluid_gradient2,
                             diffusivity = diffusivity)]

bc_apps

list examples.dg.diffusion.diffusion2d.bc_apps = []

center_bump

list examples.dg.diffusion.diffusion2d.center_bump = [0, 0, 0]

cfl

float examples.dg.diffusion.diffusion2d.cfl = cfl_factor/(2*sorder - 1)

cfl_factor

float examples.dg.diffusion.diffusion2d.cfl_factor = 0.8

detect_nonscalar

examples.dg.diffusion.diffusion2d.detect_nonscalar

dg

warpy examples.dg.diffusion.diffusion2d.dg

Initial value:

=  warpy.spatial_solvers.dg(
    name='dg',
    spatial_order=sorder_str,
    applications=apps,
    penalty_beta=penalty_beta,
    penalty_eta=penalty_eta,
    cfl=cfl)

diffusivity

float examples.dg.diffusion.diffusion2d.diffusivity = 1.0

dt

int examples.dg.diffusion.diffusion2d.dt = 1e-4

dt_calc

str examples.dg.diffusion.diffusion2d.dt_calc = "stability"

dt_controller

warpy examples.dg.diffusion.diffusion2d.dt_controller = warpy.dt_calc.stability_dt(init_dt=dt)

else

examples.dg.diffusion.diffusion2d.else :

fluid

warpy examples.dg.diffusion.diffusion2d.fluid = warpy.variable(name='fluid',components=['rho'],basis_array_set=sorder_str)

fluid_gradient

warpy examples.dg.diffusion.diffusion2d.fluid_gradient = warpy.variable(name='fluid_gradient',components=['rho_x','rho_y','rho_z'],basis_array_set=sorder_str,explodable=False)

fluid_gradient2

warpy examples.dg.diffusion.diffusion2d.fluid_gradient2 = warpy.variable(name='fluid_gradient2',components=['rho_x','rho_y','rho_z'],basis_array_set=sorder_str,explodable=False)

gen_xdmf

examples.dg.diffusion.diffusion2d.gen_xdmf

grad_bc_apps

list examples.dg.diffusion.diffusion2d.grad_bc_apps = []

gradient_method

str examples.dg.diffusion.diffusion2d.gradient_method = "ldg"

ic_fluid

warpy examples.dg.diffusion.diffusion2d.ic_fluid

Initial value:

=  warpy.apps.functions.arbitrary.bump(
    name = 'ic_fluid',
    center = center_bump,
    width = width_bump,
    minimum = min_bump,
    maximum = max_bump,
    variable = fluid,
    variable_components = None,
)

ic_ha

warpy examples.dg.diffusion.diffusion2d.ic_ha

Initial value:

=   warpy.host_actions.va_runner(name = 'ic_ha',
                                      variable_adjusters = [ic_va])

ic_va

warpy examples.dg.diffusion.diffusion2d.ic_va

Initial value:

=  warpy.variable_adjusters.function_evaluation.function_evaluator(
    name = 'ic_va',
    priority = 0,
    spatial_order = sorder_str,
    on_subdomains = ['all'],
    applications = [ic_fluid])

local_path

examples.dg.diffusion.diffusion2d.local_path

max_bump

int examples.dg.diffusion.diffusion2d.max_bump = 1

mesh

warpy examples.dg.diffusion.diffusion2d.mesh

Initial value:

=  warpy.mesh.block(Bounds=[xmin, xmax, ymin, ymax],
                        NumCells=[20,20],
                        # NumCells=[4,4],                        
                        NodeSets=['Left', 'Right','Bottom','Top'],
                        NumLayers=1,
                        PeriodicBoundaries=None,
                        basis_array_set=sorder_str,block_type='symmetric')

min_bump

int examples.dg.diffusion.diffusion2d.min_bump = 0

nframe

int examples.dg.diffusion.diffusion2d.nframe = 10

penalty_beta

float examples.dg.diffusion.diffusion2d.penalty_beta = None

penalty_eta

float examples.dg.diffusion.diffusion2d.penalty_eta = 1.0

sim

warpy examples.dg.diffusion.diffusion2d.sim

Initial value:

=  warpy.dg_sim(name='diffusion2d_example',
                   meshes=[mesh],
                   initial_conditions=[ic_ha],
                   temporal_solvers=[temporal_solver],
                   writers=[writer],
                   time=[0,t_tot],
                   dt_controller=dt_controller,
                   flexible_writeout=False,
                   write_steps=nframe,
                   verbosity='info',
                   enable_fpe=True,
                   num_partitions=1)

sorder

int examples.dg.diffusion.diffusion2d.sorder = 3

sorder_str

list examples.dg.diffusion.diffusion2d.sorder_str = sorders[sorder - 1]

sorders

list examples.dg.diffusion.diffusion2d.sorders = ['firstOrder','secondOrder','thirdOrder','fourthOrder']

t_tot

float examples.dg.diffusion.diffusion2d.t_tot = 1.0

temporal_solver

warpy examples.dg.diffusion.diffusion2d.temporal_solver

Initial value:

=  warpy.host_actions.erk(
    name='rk',
    scheme=torder,
    spatial_solvers=[dg],
    variable_adjusters=variable_adjusters)

torder

str examples.dg.diffusion.diffusion2d.torder = 'SSPRK3'

True

examples.dg.diffusion.diffusion2d.True

variable_adjusters

list examples.dg.diffusion.diffusion2d.variable_adjusters = []

width_bump

pi examples.dg.diffusion.diffusion2d.width_bump = pi

writer

warpy examples.dg.diffusion.diffusion2d.writer = warpy.host_actions.writer(name='writer',ReadVars=[fluid, fluid_gradient, fluid_gradient2])

xmax

pi examples.dg.diffusion.diffusion2d.xmax = pi

xmin

pi examples.dg.diffusion.diffusion2d.xmin = -pi

ymax

pi examples.dg.diffusion.diffusion2d.ymax = pi

ymin

pi examples.dg.diffusion.diffusion2d.ymin = -pi