Example : Heatflow Script

{
Script:     SLIDER.PDE
This problem represents a cross section of a wood-frame sliding window.
(submitted by Elizabeth Finleyson, Lawrence Berkeley Lab.)
}

TITLE  "NFRC Wood Slider"

VARIABLES

Temp

DEFINITIONS

K = 0.97 {Thermal Conductivity}
B1 = 1.34 {Film coefficients interior wood}
B2 = 1.41 { Interior glass}
B3 = 5.11 { Exterior glass}
Tin = 70.0 {Ambient Temperature Inside}
Tout= 0.0 { Outside}

EQUATIONS

dx(K*dx(Temp)) + dy(K*dy(Temp)) = 0

BOUNDARIES

REGION 1 {Defines the maximum extent of the system (wood) }

start(6.813,1.813)
natural(Temp) = B1*(Tin - Temp)
line to (6.813,3.001) to (6.344,3.001) to (6.344,3.323)
line to (6.183,3.323) to (6.183,4.885) to (5.988,4.885)
line to (5.988,5.104) to (5.678,5.104)

natural(Temp) = B2*(Tin - Temp)
line to (5.678,7.604)

natural(Temp) = 0.0
line to (5.153,7.604)

natural(Temp) = B3*(Tout- Temp)
line to (5.153,5.104) to (5.012,5.104) to (5.012,4.889)
line to (4.871,4.889) to (4.871,3.323) to (4.248,3.323)
line to (4.248,2.845) to (3.233,2.845) to (3.233,3.323) to (2.906,3.323)
line to (2.906,3.001) to (2.250,3.001) to (2.250,2.501) to (1.156,2.501)

natural(Temp) = 0.0
line to (1.156,1.813) to close

REGION 2 K = 1.18 {Rigid PVC}

start(6.516,2.800)
line to (6.516,2.845) to (6.344,2.845) to (6.344,3.323)
line to (5.737,3.323) to (5.737,3.278) to (6.017,3.278)
line to (6.017,2.845) to (5.002,2.845) to (5.002,3.278)
line to (5.317,3.278) to (5.317,3.323)
line to (4.248,3.323) to (4.248,2.845) to (3.233,2.845)
line to (3.233,3.323) to (2.906,3.323) to (2.906,2.845)
line to (2.547,2.845) to (2.547,2.800) to close

REGION 3 K = 0.59 {air cavity overlays}

start(4.293,2.845)
line to (4.957,2.845) to (4.957,3.278) to (4.293,3.278) to close

REGION 4 k = 0.31

start(2.951,2.800)
line to (3.188,2.800) to (3.188,3.278) to (2.951,3.278) to close

REGION 5 k = 0.51

start(2.547,2.501)
line to (3.188,2.501) to (3.188,2.800) to (2.547,2.800) to close

REGION 6 k = 0.81

start(5.002,2.845)
line to (6.017,2.845) to (6.017,3.278) to (5.002,3.278) to close

REGION 7 k = 0.39

start(5.317,3.278)
line to (5.737,3.278) to (5.737,3.551) to (5.317,3.551) to close

REGION 8 k = 0.31

start(6.062,2.800)
line to (6.299,2.800) to (6.299,3.278) to (6.062,3.278) to close

REGION 9 k = 0.41

start(6.062,2.501)
line to (6.516,2.501) to (6.516,2.800) to (6.062,2.800) to close

REGION 10 k = 2.5 {Silicon sealant}

start(5.133,4.573)
line to (5.153,4.573) to (5.153,5.104) to (5.133,5.104) to close

REGION 11 k = 2.5

start(5.678,4.573)
line to (5.698,4.573) to (5.698,5.104) to (5.678,5.104) to close

REGION 12 k = 6.93 {glass layers}

start(5.153,4.573)
line to (5.678,4.573) to (5.678,7.604) to (5.153,7.604) to close

REGION 13 k = 2.5 {eurythane spacer seal}

start(5.278,4.573)
line to (5.553,4.573) to (5.553,4.771) to (5.278,4.771) to close

REGION 14 k = 18.44 {spacer}

start(5.278,4.771)
line to (5.553,4.771) to (5.553,5.012) to (5.278,5.012) to close

REGION 15 k = 0.32 {gas gap}

start(5.278,5.012)
line to (5.553,5.012) to (5.553,7.604) to (5.278,7.604) to close

REGION 16 k = 0.21 {frame fill}

start(3.188,2.501)
line to (6.062,2.501) to (6.062,2.800) to (3.188,2.800) to close

REGION 17 k = 0.28 {spacer air gap}

start(5.133,4.479)
line to (5.698,4.479) to (5.698,4.573) to (5.133,4.573) to close

MONITORS

contour(Temp)

PLOTS

contour(Temp)
contour(Temp) zoom(4.6,4.2,1.8,1.8)
elevation(Temp) from (5.416,1.813) to (5.416,7.604)
vector((K*(-dx(Temp))),(K*(-dy(Temp)))) as "HEAT FLUX"

END