Section MultiScale RealUnits

UNITCODE RealUnits
DIMENSION 2

###### constants ######
LET sqrt3 = 1.73205080756887729353
LET sqrt2 = 1.4142135623730950488
LET r0 = 4.0320*sqrt2/2
LET r1 = r0*sqrt3

LET L0 = r0*200
LET L1 = r1*40

LET TSIZE = 5*r0
LET iTemp = 5
LET fTemp = 20
LET tol = r0*1e-3

###### geometries ######
GEOMETRY feGeom	  CUBE BBOX -L0/2-tol L0/2+tol -L1/2-tol L1/2+tol
GEOMETRY left	  CUBE BBOX -L0/2-tol -L0/2+TSIZE/2 -L1/2-tol L1/2+tol
GEOMETRY right	  CUBE BBOX L0/2-TSIZE/2 L0/2+tol -L1/2-tol L1/2+tol

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

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

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

###### dumpers ######
DUMPER para PARAVIEW INPUT fe TEMPERATURE TEMPERATURE_VARIATION HEAT_RATE DISP VEL FORCE FREQ 100
DUMPER left PARAVIEW INPUT leftGeom TEMPERATURE TEMPERATURE_VARIATION HEAT_RATE DISP VEL FORCE ONESHOT 0
DUMPER right PARAVIEW INPUT rightGeom TEMPERATURE TEMPERATURE_VARIATION HEAT_RATE DISP VEL FORCE ONESHOT 0

COMPUTE ethermalFE ETHERMAL INPUT fFeGeom
DUMPER energy-mesh TEXT INPUT ethermalFE FREQ 100 PREFIX ./

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
DUMPER temperature TEXT INPUT average-temp-left ADD_COMPUTE average-temp-right PRECISION 5 FREQ 100

endSection

################################################################
#FE section
################################################################
Section AKANTU:fe RealUnits
DOMAIN_GEOMETRY feGeom
TIMESTEP 1.0
MESH_FILENAME mesh.msh
MATERIAL_FILENAME material.dat
PBC 0 1
#SCALE L0/40 L1/6
#SHIFT -L0/2 -L1/2

HEAT_TRANSFER
HEAT_TRANSFER_MATERIAL_FILENAME heat_properties.dat
HEAT_BOUNDARY left
endSection