HeatConduction

This is the test for the Integrator::HeatConduction integrator. All BCs are Dirichlet. Low faces have constant nonzero temperature, high faces have constant zero temperature. Domain is initialized with zero temperature. Sample result below in 3D:

2d-serial

	Two-dimensional
	Serial
	Validated using check script
	./bin/alamo-2d-g++ tests/HeatConduction/input

2d-serial-coverage

	Two-dimensional
	Serial
	Not validated
	./bin/alamo-2d-g++ tests/HeatConduction/input stop_time="0.0002"

3d-parallel

	Three-dimensional
	Parallel (4 procs)
	Validated using check script
	7.66s (beaker) 21.17s (waldorf)
	mpiexec -np 4 ./bin/alamo-3d-g++ tests/HeatConduction/input

3d-serial-coverage

	Three-dimensional
	Serial
	Not validated
	./bin/alamo-3d-g++ tests/HeatConduction/input stop_time="0.0002"

Input file (../../tests/HeatConduction/input)

#@
#@ [2d-serial]
#@ dim=2
#@ 
#@ [2d-serial-coverage]
#@ dim=2
#@ coverage=true
#@ args=stop_time=0.0002
#@ check=false
#@
#@ [3d-parallel]
#@ dim=3
#@ nprocs=4
#@ benchmark-beaker=7.66
#@ benchmark-waldorf=21.17
#@
#@ [3d-serial-coverage]
#@ dim=3
#@ coverage=true
#@ args=stop_time=0.0002
#@ check=false

alamo.program = heat

plot_file     = tests/HeatConduction/output

# Simulation length
timestep = 0.0001
stop_time = 0.05

# AMR parameters
amr.plot_int = 10
amr.max_level = 3
amr.n_cell = 8 8 8
amr.blocking_factor = 1
amr.regrid_int = 10
amr.grid_eff = 1.0

# Specify geometry and unrefined mesh
geometry.prob_lo = 0 0 0
geometry.prob_hi = 1 1 1
geometry.is_periodic= 0 0 0

# Criterion for mesh refinement
heat.alpha = 1.0
heat.refinement_threshold = 0.01

# Specify initial conditions
ic.type = sphere
ic.sphere.center = 0.5 0.5 0.5
ic.sphere.radius = 0.25
ic.sphere.inside = 0.0



# Boundary conditions
bc.temp.constant.type.xhi = dirichlet
bc.temp.constant.type.xlo = dirichlet
bc.temp.constant.type.yhi = dirichlet
bc.temp.constant.type.ylo = dirichlet
bc.temp.constant.type.zhi = dirichlet
bc.temp.constant.type.zlo = dirichlet
bc.temp.constant.val.xhi  = 0.0
bc.temp.constant.val.xlo  = 1.0
bc.temp.constant.val.yhi  = 0.0
bc.temp.constant.val.ylo  = 1.0
bc.temp.constant.val.zhi  = 0.0
bc.temp.constant.val.zlo  = 1.0