ElasticSoftVoid

2d-serial

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

2d-parallel

	Two-dimensional
	Parallel (2 procs)
	Validated using check script
	mpiexec -np 2 ./bin/mechanics-2d-g++ tests/ElasticSoftVoid/input

Input file (../../tests/ElasticSoftVoid/input)

#@ [2d-serial]
#@ exe = mechanics
#@ dim = 2
#@ nprocs = 1
#@ check = true
#@ timeout = 120
#@
#@ [2d-parallel]
#@ exe = mechanics
#@ dim = 2
#@ nprocs = 2
#@ check = true
#@ timeout = 120
#@

alamo.program = mechanics
alamo.program.mechanics.model = finite.neohookean

type = static
timestep = 0.1
stop_time = 0.1

amr.plot_int = 1
amr.max_level = 1
amr.n_cell = 32 32 8
amr.blocking_factor = 4
amr.max_grid_size = 16
amr.regrid_int = -1
amr.grid_eff = 1.0
amr.node.all = 1
amr.cell.all = 1

explicitmesh.on = 1
explicitmesh.lo1 = 8 8 0
explicitmesh.hi1 = 23 23 0

geometry.prob_lo = 0 0 0
geometry.prob_hi = 1 1 1
geometry.is_periodic = 0 0 0

ic.type = ellipse
ic.ellipse.a = 0.22 0.22 0.22
ic.ellipse.x0 = 0.5 0.5 0.5
ic.ellipse.eps = 0.04

nmodels = 2
model1.mu = 1
model1.kappa = 1
model2.mu = 26000
model2.kappa = 70000

rhs.type = expression
rhs.expression.region0 = "-(x-0.5)/0.04*exp(-((sqrt((x-0.5)^2+(y-0.5)^2)-0.22)/0.04)^2)"
rhs.expression.region1 = "-(y-0.5)/0.04*exp(-((sqrt((x-0.5)^2+(y-0.5)^2)-0.22)/0.04)^2)"

solver.conservative_face_flux = 1
solver.line_search = 1
solver.resync_coeffs = 1
solver.nriters = 3
solver.nrtolerance = 0
solver.max_iter = 1500
solver.pre_smooth = 2
solver.post_smooth = 2
solver.bottom_solver = bicgstab
solver.verbose = 3

max_coarsening_level = 2
zero_out_displacement = 1
tol_rel = 2.0e-5
tol_abs = 1.0e-8
ref_threshold = 100

bc.type = constant
bc.constant.type.xlo = disp disp
bc.constant.type.xhi = disp trac
bc.constant.type.ylo = trac disp
bc.constant.type.yhi = disp disp
bc.constant.type.xloylo = disp disp
bc.constant.type.xloyhi = disp disp
bc.constant.type.xhiylo = disp disp
bc.constant.type.xhiyhi = disp disp