Slab1D

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