Section MultiScale RealUnits

UNITCODE RealUnits
DIMENSION 1

###### constants ######
LET r0 = 1.23157195294235
LET rho = 39.95/r0
PRINT rho
LET ELSIZE = 8*r0
LET NBATOMS = 5000
LET L = NBATOMS*r0
LET rcut = 2.2*r0
LET TSIZE = 100*r0
LET iTemp = 5
LET fTemp = 20
LET damp = 1000

###### geometries ######
GEOMETRY feGeom   BALL CENTER 0 RADII 0 L/2
GEOMETRY left	  CUBE BBOX -L/2 -L/2+TSIZE
GEOMETRY right	  CUBE BBOX L/2-TSIZE L/2

###### model declaration ######
ELAST_CODE AKANTU fe

###### filters ######
FILTER leftGeom GEOM INPUT fe GEOMETRY left
FILTER rightGeom GEOM INPUT fe GEOMETRY right

###### stimulations ######
STIMULATION init_temp TEMPERATURE INPUT fe TEMP iTemp ONESHOT 0
STIMULATION init_thermostat TEMPERATURE INPUT leftGeom FREQ 1 TEMP fTemp START 1001

###### computes ######
COMPUTE temp-left  EXTRACT INPUT leftGeom  FIELD temperature
COMPUTE temp-right EXTRACT INPUT rightGeom FIELD temperature
COMPUTE average-temp-left  REDUCE INPUT temp-left  OPERATOR AVERAGE
COMPUTE average-temp-right REDUCE INPUT temp-right OPERATOR AVERAGE


###### dumpers ######
DUMPER para PARAVIEW INPUT fe TEMPERATURE DISP VEL FORCE FREQ 1000
DUMPER temperature TEXT INPUT average-temp-left ADD_COMPUTE average-temp-right PRECISION 5 FREQ 1000

endSection


################################################################
#FE section
################################################################
Section AKANTU:fe RealUnits
DOMAIN_GEOMETRY feGeom
TIMESTEP 2
ELEM_SIZE ELSIZE
MASS 39.95
R0 r0
RCUT rcut
EPSILON 1.6567944e-21/6.94769e-21
SIGMA 1.1
HEAT_TRANSFER_MATERIAL_FILENAME heat_properties.dat
HEAT_TRANSFER
HEAT_BOUNDARY left
endSection