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

###### geometries ######
GEOMETRY feGeom     BALL CENTER 0 RADII 0 L/2
GEOMETRY feBoundary BALL CENTER 0 RADII L/2-10*r0 L/2
###### model declaration ######
ELAST_CODE AKANTU fe

###### computes ######
COMPUTE positions0 EXTRACT INPUT fe FIELD position0
COMPUTE Tgradient PYTHON INPUT positions0 FILENAME wl

###### stimulations ######
STIMULATION Tgradient FIELD INPUT fe FIELD temperature COMPUTE Tgradient

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

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 feBoundary
endSection