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*800
LET L1 = r1*10

LET TSIZE = 100*r0
LET iTemp = 5
LET fTemp = 20
LET tol = r0/2

###### 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 -L1/2-tol L1/2+tol
GEOMETRY right	  CUBE BBOX L0/2-TSIZE 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

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

###### dumpers ######
DUMPER para PARAVIEW INPUT fe TEMPERATURE TEMPERATURE_VARIATION HEAT_RATE DISP VEL FORCE FREQ 10

LET nSliceFEM = L0/TSIZE
FOR slice in range 0 ${nSliceFEM} 1
LET rxMin = -L0/2+TSIZE*slice
LET rxMax = -L0/2+TSIZE*(slice+1)
GEOMETRY feGeom${slice} CUBE BBOX rxMin-tol rxMax+tol -L1/2-tol L1/2+tol
FILTER geom${slice} GEOM INPUT fe GEOMETRY feGeom${slice}
COMPUTE temp${slice}  EXTRACT INPUT geom${slice}  FIELD temperature
COMPUTE average-temp${slice}  REDUCE INPUT temp${slice}  OPERATOR AVERAGE
DUMPER temperature${slice} TEXT INPUT average-temp${slice} PRECISION 5 FREQ 10
DUMPER zone${slice} PARAVIEW INPUT geom${slice} TEMPERATURE FREQ 10
DUMPER zone-text${slice} TEXT INPUT temp${slice} FREQ 10 PRECISION 5
endFOR

#DUMPER temperature TEXT INPUT average-temp7 ADD_COMPUTE average-temp7 PRECISION 5 FREQ 10


endSection

################################################################
#FE section
################################################################
Section AKANTU:fe RealUnits
DOMAIN_GEOMETRY feGeom
TIMESTEP 1
MESH_FILENAME mesh.msh
MATERIAL_FILENAME material.dat
PBC 0 0

HEAT_TRANSFER
HEAT_TRANSFER_MATERIAL_FILENAME heat_properties.dat
#HEAT_BOUNDARY left
endSection