slider

<< Click to Display Table of Contents >>

Navigation:  Sample Problems > Applications > Heatflow >

slider

Previous pageReturn to chapter overviewNext page

{  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