Heat Flow

{  SLIDER.PDE  

 This problem represents a cross section of a wood-frame sliding window.
     -- submitted by Elizabeth Finleyson, Lawrence Berkeley Labs

} 

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
 Temp : 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  

{Rigid PVC}
region 2      K = 1.18
  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  

{Air cavity overlays}
region 3     K = 0.59
  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  

{Silicon sealant}
region 10  k = 2.5
  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  

{Glass layers}
region 12  k = 6.93
  start(5.153,4.573)
  line to (5.678,4.573) to (5.678,7.604) to (5.153,7.604) to close  

{Eurythane spacer seal}
region 13  k = 2.5
  start(5.278,4.573)
  line to (5.553,4.573) to (5.553,4.771) to (5.278,4.771) to close  

{Spacer}
region 14  k = 18.44
  start(5.278,4.771)
  line to (5.553,4.771) to (5.553,5.012) to (5.278,5.012) to close  

{Gas gap}
region 15  k = 0.32
  start(5.278,5.012)
  line to (5.553,5.012) to (5.553,7.604) to (5.278,7.604) to close  

{Frame fill}
region 16  k = 0.21
  start(3.188,2.501)
  line to (6.062,2.501) to (6.062,2.800) to (3.188,2.800) to close  

{Spacer air gap}
region 17  k = 0.28
  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
grid(x,y)
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