EshelbyDynamics

This test demonstrates the use of the dynamics solver using the Mechanics integrator. It is similar to the Eshelby test problem: the inclusion is given an eigenstrain and compressed. This creates a stress wave at the interface into and out of the inclusion. Viscosity is used to damp out spurious oscillations, and so the solution converges fairly quickly to the static solution.

../_images/movie1.gif

2D-serial-4levels

Two-dimensional

Serial

Validated using check script

42s (beaker)

 

./bin/mechanics-2d-g++ tests/EshelbyDynamics/input  stop_time="0.5"
Input file (../../tests/EshelbyDynamics/input)
#@
#@  [2D-serial-4levels]
#@  exe    = mechanics
#@  dim    = 2
#@  nprocs = 1
#@  args   = stop_time=0.5
#@  check  = true
#@  benchmark-beaker = 42
#@

alamo.program = mechanics
alamo.program.mechanics.model = affine.isotropic

plot_file		    = tests/EshelbyDynamics/output

type=dynamic

# this is not a time integration, so do
# exactly one timestep and then quit
timestep		    = 0.0005
stop_time		    = 10.0

# amr parameters
amr.plot_dt		    = 0.01
amr.max_level		    = 4
amr.n_cell		    = 64 64 64
amr.blocking_factor	    = 8
amr.regrid_int		    = 100
amr.grid_eff		    = 1.0
amr.cell.all                = 1

# use an explicit mesh (i.e. no adaptive meshing)
explicitmesh.on               = 0
explicitmesh.lo1              = 16 16 16
explicitmesh.hi1              = 47 47 47
explicitmesh.lo2              = 48 48 48
explicitmesh.hi2              = 79 79 79
explicitmesh.lo3              = 112 112 112
explicitmesh.hi3              = 145 145 145

# geometry
geometry.prob_lo	    = -8 -8 -8 
geometry.prob_hi	    = 8 8 8
geometry.is_periodic	    = 0 0 0

# ellipse configuration
ic.type        = ellipse
ic.ellipse.a   = 1.0 0.75 0.5  # ellipse radii
ic.ellipse.x0  = 0 0 0 # location of ellipse center
ic.ellipse.eps = 0.1 # diffuse boundary

# elastic moduli
nmodels = 2
model1.E = 210 
model1.nu = 0.3
model1.F0  = 0.001 0 0 0 0.001 0 0 0 0.001 # eigenstrain
model2.E = 210 
model2.nu = 0.3
model2.F0  = 0 0 0 0 0 0 0 0 0 # eigenstrain

ref_threshold = 1000000000

viscous.mu_dashpot = 0.1
viscous.mu_newton  = 1.0