TITLE 'Heat flow through an Insulating layer in 1D'
COORDINATES
Cartesian1
VARIABLES
Phi { the temperature }
DEFINITIONS
K = 1 { default conductivity }
R = 0.5 { insulator thickness }
EQUATIONS
Div(-k*grad(phi)) = 0
BOUNDARIES
REGION 1 { the total domain }
START(-1) POINT VALUE(Phi)=0
LINE TO (1) POINT VALUE(Phi)=1
{ note: no 'close'! }
REGION 2 'blob' { the embedded layer }
k = 0.001
START (-R) LINE TO (R)
PLOTS
ELEVATION(Phi) FROM (-1) to (1)
END