Difference between revisions of "Heat transfer"

From UMaine SECS Numerical Modeling Laboratory
Jump to: navigation, search
 
(25 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
''Created by Clara Deck, Mariama Dryak and He Feng on 2019-03-14. Because this page has some high-resolution animations, it is best viewed on a reliable internet connection.''
 +
 +
== '''Introduction''' ==
 +
 +
Welcome to our heat transfer module! Here we will explain the basic principles of heat transfer using some applied examples of heat transfer. During our discussion of heat transfer in a lava lamp, we introduce you to the different components of the heat transfer equation as applied to our lava lamp and then bring the different components of this equation all together at the end. This heat transfer equation and other critical mathematical terms like the Peclét and Rayleigh Numbers are then applied to examples of heat transfer in a room and ice stream start-up and shut-down.
  
 
== '''Heat Transfer in a Lava Lamp''' ==
 
== '''Heat Transfer in a Lava Lamp''' ==
Line 21: Line 26:
 
**velocity of the fluid (v)
 
**velocity of the fluid (v)
 
**thermal gradient (∇T)
 
**thermal gradient (∇T)
 +
  
  
Line 37: Line 43:
 
The Peclét Number is a dimensionless number which indicates whether conduction or advection dominates the system.
 
The Peclét Number is a dimensionless number which indicates whether conduction or advection dominates the system.
 
   Pe = uL/κ
 
   Pe = uL/κ
where u = flow velocity, L = characteristic length Height of lava lamp, κ = thermal diffusivity
+
where u = flow velocity, L = characteristic length (height of lava lamp), κ = thermal diffusivity
 
*Low peclet regime- dominated by heat transport by conduction
 
*Low peclet regime- dominated by heat transport by conduction
 
*High peclet regime- dominated by heat transport by advection
 
*High peclet regime- dominated by heat transport by advection
Line 43: Line 49:
 
With the movement of wax in the upward direction away from the heat source at the bottom of the container and into a cooler material the wax density increases again (more than liquid). This wax density increase causes the wax to sink back down to the bottom of the container towards the heat source. This process is called convection.
 
With the movement of wax in the upward direction away from the heat source at the bottom of the container and into a cooler material the wax density increases again (more than liquid). This wax density increase causes the wax to sink back down to the bottom of the container towards the heat source. This process is called convection.
  
[[File:lava5.png|230px|Image: 230 pixels]] [[File:isotherm3.png|180px|Image: 180 pixels]]
+
[[File:lava5.png|230px|Image: 230 pixels]] [[File:isotherm3.png|185px|Image: 185 pixels]]
  
  
Line 63: Line 69:
  
 
At some critical Rayleigh (Racr) number (dependent on the system), the convection regime shifts.
 
At some critical Rayleigh (Racr) number (dependent on the system), the convection regime shifts.
 +
 +
== '''Indoor Natural Convection in Different Rayleigh Numbers''' ==
 +
 +
'''Natural Convection'''
 +
 +
Natural convection is a condition in which there is '''no external driving force''' but the fluid is still in motion. In natural convection, the fluid surrounding the heat source '''receives heat''' and becomes '''less dense''' and '''rises''' by thermal expansion. Then the cooler fluid moves around to replace it. The cooler fluid is then heated, and the process continues, which creats convection currents.
 +
 +
Natural convection plays an important role in earth sciences, such as atmospheric sciences, oceans, and the earth's mantle. Indoor heating, which we often used, is to use this principle to form complex indoor heat transfer mechanism. In this module, I use a finite volume numerical method to analyze the fluid temperature field, flow field distribution characteristic under different Rayleigh numbers, in order to study the indoor natural convection and heat transfer process.
 +
 +
[[File:Indoorconvection.png|350px|Image: 350 pixels]]
 +
 +
'''Governing equations'''
 +
 +
Continuity equation (Law of conservation of mass in fluid mechanics):
 +
 +
  ∂ρ/∂t+∇∙(ρ'''v''')=0
 +
 +
Navier-Stokes equation (Newton's Second Law of Motion-Force and Accelerat in fluid mechanics):
 +
 +
  m'''a'''='''F'''
 +
 +
  Inertial forces = viscous forces +body forces - pressure gradient
 +
 +
 +
  ρ (d'''v''')/dt  = μ∇² '''v''' + ∆ρg - ∇P
 +
 +
Because temperature affects fluid density, buoyancy changes the term of force in the N-S equation.
 +
 +
  ρ (d'''v''')/dt  = μ∇² '''v''' + ∆ρgα(T-T<sub>0</sub>) - ∇P
 +
 +
'''v''' : vector of fluid velocity;
 +
 +
'''a''': acceleration;
 +
 +
m: mass;
 +
 +
'''F''': force exerted on a body;
 +
 +
P: pressure;
 +
 +
T: temperature;
 +
 +
T<sub>0</sub>: a reference temperature;
 +
 +
g: gravity acceleration;
 +
 +
ρ: reference density;
 +
 +
μ: dynamic viscosity;
 +
 +
α: coefficient of volumetric thermal expansion.
 +
 +
 +
'''Model setup'''
 +
 +
We defined a 1m*1m square cavity filled with the fluid we defined above (reference pressure p=0, reference temperature T<sub>0</sub>=0). The four sides are no-sliping Wall. The right side is constant temperature (T{<sub>0</sub>=0), and the lower right wall and bottom edge 0.1m*0.1m are heat sources (T=1). Pressure Point Constraint with the four corners set to p=0.We use this model to simulate the natural convection of indoor enclosed space caused by heating (modified from Buoyancy Flow in Free Fluids by COMSOL Multiphysics).
 +
 +
[[File:Model-INC.png|400px|Image: 400 pixels]] [[File:Mesh.png|400px|Image: 400 pixels]]
 +
 +
I examine a wide range of scenarios sequencing through different values of Ra using an Auxiliary sweep with continuation selected for the Ra parameter in the Stationary study step (Ra=[100 1e3 1e4 1e5 1e6]). With the increase of Ra, the importance of viscous force decreases.
 +
 +
For the convenience of numerical simulation, I defined a fluid that does not exist in reality, whose fluid parameters are shown in the following table (Pr: Prandtl number, the ratio of kinematic viscosity to thermal diffusivity.):
 +
 +
  Pr=(μC<sub>p</sub>)/k
 +
 +
μ: dynamic viscosity;
 +
 +
C<sub>p</sub>: specific heat capacity of the fluid;
 +
 +
k = thermal diffusivity
 +
 +
 +
{| class="wikitable"
 +
|-
 +
! Name
 +
! Value
 +
|-
 +
| ρ
 +
| 1
 +
|-
 +
| μ
 +
| 1
 +
|-
 +
| k
 +
| 1
 +
|-
 +
| Pr
 +
| 0.71
 +
|-
 +
| Ra
 +
| Changed as setting
 +
|}
 +
 +
 +
So '''numerically''' ('''ignoring units''') :
 +
 +
  Cp=Pr
 +
 +
body force in the y-direction for the momentum equation:
 +
 +
  F=(Ra/Pr)(T-T<sub>0</sub>)
 +
 +
'''Result and Discussion'''
 +
 +
We analyze our results from the aspects of velocity field and temperature field. Free convective flow is completely caused by temperature difference and the velocity field is generated by the temperature field and in turn affects the temperature field.
 +
 +
'''''Time to stable state'''''
 +
 +
We found that with the increase of Ra, the time of heat transfer to the relatively stable state within the space became shorter.
 +
 +
Ra=100 (in ~250ms)
 +
 +
[[File:Velocity100.gif|500px|Image: 500 pixels]] [[File:Temp100.gif|500px|Image: 500 pixels]]
 +
 +
Ra=1e4 (in ~130ms)
 +
 +
[[File:Velocity1e4.gif|500px|Image: 500 pixels]] [[File:Temp1e4.gif|500px|Image: 500 pixels]]
 +
 +
Ra=1e6 (in ~70ms)
 +
 +
[[File:Velocity1e6.gif|500px|Image: 500 pixels]] [[File:Temp1e6.gif|500px|Image: 500 pixels]]
 +
 +
'''''Velocity field'''''
 +
 +
With the increase of Ra, the flow velocity of the fluid in stable state also increases significantly. Convection is near the boundary layer near the hot (cold) wall. Due to the influence of the heat source area on the left, the fluid flows upward along the left wall surface to the top of the square cavity. Due to the limitation of the height of the thermal insulation wall surface, the flow line turns to the right. Due to the influence of the right cold wall, the fluid flows down the right cold wall to the bottom of the square cavity. Due to the limitation of the height of the adiabatic wall, the streamline turns left and flows to the left hot wall. The whole flow field flows in a clockwise direction. When Ra= 100 or 1e3, there is a vortex at the center of the flow field. With the increase of Rayleigh number, the vortices in the flow field are divided into two symmetrical vortices and move towards the hot (cold) wall continuously.
 +
 +
[[File:RaU.png|1100px|Image: 1100 pixels]]
 +
 +
'''''Temperature field'''''
 +
 +
Convective heat transfer is the result of the combined action of fluid flow and heat conduction between fluid and molecules. When Ra = 100 or 10e3 (relatively low), the isotherm is uniformly arched outward with the heat source as the center, and the temperature gradient near the heat source is large. With the increase of Ra, the heat transfer changes from heat conduction to heat convection, and the isotherm gradually bends and deforms, forming a thin boundary layer near the cold wall and the heat source.
 +
 +
[[File:RaT.png|1100px|Image: 1100 pixels]]
  
 
=='''Heat Transfer in Ice Stream Shut Down and Start (based on Joughin and Alley, 2011)'''==
 
=='''Heat Transfer in Ice Stream Shut Down and Start (based on Joughin and Alley, 2011)'''==
Line 74: Line 213:
 
Watch the video below to walk through the steps of ice stream shut-down and speed-up, and see how this relates to the heat transfer equation introduced above.  
 
Watch the video below to walk through the steps of ice stream shut-down and speed-up, and see how this relates to the heat transfer equation introduced above.  
  
{{#ev:youtube|"https://youtu.be/JZWqx4qiGJ0"}}
+
{{#ev:youtube|"https://youtu.be/Sae8Rh-oLX0"}}
  
 
'''Further reading:'''
 
'''Further reading:'''
Line 86: Line 225:
 
*Lipovsky, B., and E. Dunham (2016) Tremor during ice-stream stick slip, Cryosphere, 10(1), pp. 385–399. https://www.the-cryosphere.net/10/385/2016/tc-10-385-2016.html
 
*Lipovsky, B., and E. Dunham (2016) Tremor during ice-stream stick slip, Cryosphere, 10(1), pp. 385–399. https://www.the-cryosphere.net/10/385/2016/tc-10-385-2016.html
 
*Price, S.F., Bindschadler, R.A., Hulbe, C.L. and Joughin, I.R. (2001) Post-stagnation behavior in the upstream regions of Ice Stream C, West Antarctica. Journal of Glaciology, 47(157), pp. 283-294. https://doi.org/10.3189/172756501781832232
 
*Price, S.F., Bindschadler, R.A., Hulbe, C.L. and Joughin, I.R. (2001) Post-stagnation behavior in the upstream regions of Ice Stream C, West Antarctica. Journal of Glaciology, 47(157), pp. 283-294. https://doi.org/10.3189/172756501781832232
 +
 +
== '''Summary''' ==
 +
 +
This teaching module covered conduction, advection, and convection as heat transfer mechanisms. We can quantitatively characterize these processes in the heat transfer equation explained above. In the lava lamp, the transition from a conduction to advection dominated system occurs when temperature driven density reaches a critical low, allowing the ‘lava’ to rise through the surrounding fluid. The dimensionless Peclet Number can help us understand whether conduction or advection is driving most of the heat transfer. Though the same principles apply, behavior is a little different when considering the ice stream stop and start behavior. The regime shift between conduction and advection is greatly dependent on the phase change from ice to liquid water at the ice stream bed. The film of meltwater which allows advection to proceed is a product of conduction of heat from the bed. Convection, described above in the contexts of the lava lamp and heat circulation through an enclosed room, includes both conduction and advection as important mechanisms of heat transfer, but requires a spatial cycling of heat flow. The Rayleigh Number is a useful tool to evaluate the drivers and resistors of convection in a given system, and to determine whether convection will tend to grow or decay. Please check out the included links for more information or examples of heat transfer in the real world!

Latest revision as of 19:30, 19 March 2019

Created by Clara Deck, Mariama Dryak and He Feng on 2019-03-14. Because this page has some high-resolution animations, it is best viewed on a reliable internet connection.

Introduction

Welcome to our heat transfer module! Here we will explain the basic principles of heat transfer using some applied examples of heat transfer. During our discussion of heat transfer in a lava lamp, we introduce you to the different components of the heat transfer equation as applied to our lava lamp and then bring the different components of this equation all together at the end. This heat transfer equation and other critical mathematical terms like the Peclét and Rayleigh Numbers are then applied to examples of heat transfer in a room and ice stream start-up and shut-down.

Heat Transfer in a Lava Lamp

When you flip on the light in a lava lamp, heat from the bulb works slowly to warm up the wax (or ‘lava’). This kind of heat transfer is called conduction, or movement of heat from one place to another. Heat energy flows from high heat (light bulb) to low heat (wax), and the gradient is the driver of conductive heat transfer.

Image: 220 pixels Image: 220 pixels Image: 190 pixels Image: 170 pixels


Until things start moving, conduction is the dominant heat transfer process in our lava lamp system.

  • Conduction, defined as κ∇²T depends on:
    • thermal diffusivity of the material through which heat is moving (κ)
    • thermal gradient (∇²T)

The density of the wax is inversely related to temperature, meaning that at higher temperatures, the wax is less dense. When the density becomes lower than the density of the liquid in the lava lamp, the wax will tend to float upward.

Image: 215 pixels Image: 200 pixels Image: 180 pixels

Wax reaches critical density (less than liquid) and tends to float upward. That movement of wax, and thus heat, introduces velocity into the system and is a process that is referred to as advection.

  • Advection, defined as v∇T depends on:
    • velocity of the fluid (v)
    • thermal gradient (∇T)


Heat Transfer Equation:

Heat transfer within the lava lamp can be represented by an equation. The heat transfer equation brings these aforementioned processes together to describe how heat changes through time (∂T/∂t). It is defined as:

 ∂T/∂t = κ∇²T + v∇T + A
 ∂T/∂t = conduction + advection + production

The production term (A) represents internally generated heat energy and we do not consider it in this lava lamp scenario. It becomes important in systems such as glaciers, when internal strain heating is an important source of heat.

Peclét Number:

The Peclét Number is a dimensionless number which indicates whether conduction or advection dominates the system.

 Pe = uL/κ

where u = flow velocity, L = characteristic length (height of lava lamp), κ = thermal diffusivity

  • Low peclet regime- dominated by heat transport by conduction
  • High peclet regime- dominated by heat transport by advection

With the movement of wax in the upward direction away from the heat source at the bottom of the container and into a cooler material the wax density increases again (more than liquid). This wax density increase causes the wax to sink back down to the bottom of the container towards the heat source. This process is called convection.

Image: 230 pixels Image: 185 pixels


Rayleigh Number:

The Rayleigh Number is a dimensionless number used to describe whether convection grows or decays.

 Ra ≡ (gρΔTαd³)/μκ
  • Convection depends on driving forces and resisting forces.
    • Driving forces (when these terms are dominant convection grows):
      • Acceleration due to gravity (g)
      • Density (ρ)
      • The temperature difference between the bottom and top of the convection cell (ΔT)
      • The volume coefficient of thermal expansion (α)
      • The height of the convection cell (d)
    • Resisting forces (when these terms are dominant convection decays):
      • Viscosity (μ)
      • Thermal diffusivity (κ)

At some critical Rayleigh (Racr) number (dependent on the system), the convection regime shifts.

Indoor Natural Convection in Different Rayleigh Numbers

Natural Convection

Natural convection is a condition in which there is no external driving force but the fluid is still in motion. In natural convection, the fluid surrounding the heat source receives heat and becomes less dense and rises by thermal expansion. Then the cooler fluid moves around to replace it. The cooler fluid is then heated, and the process continues, which creats convection currents.

Natural convection plays an important role in earth sciences, such as atmospheric sciences, oceans, and the earth's mantle. Indoor heating, which we often used, is to use this principle to form complex indoor heat transfer mechanism. In this module, I use a finite volume numerical method to analyze the fluid temperature field, flow field distribution characteristic under different Rayleigh numbers, in order to study the indoor natural convection and heat transfer process.

Image: 350 pixels

Governing equations

Continuity equation (Law of conservation of mass in fluid mechanics):

 ∂ρ/∂t+∇∙(ρv)=0

Navier-Stokes equation (Newton's Second Law of Motion-Force and Accelerat in fluid mechanics):

 ma=F

 Inertial forces = viscous forces +body forces - pressure gradient


 ρ (dv)/dt  = μ∇² v + ∆ρg - ∇P

Because temperature affects fluid density, buoyancy changes the term of force in the N-S equation.

  ρ (dv)/dt  = μ∇² v + ∆ρgα(T-T0) - ∇P

v : vector of fluid velocity;

a: acceleration;

m: mass;

F: force exerted on a body;

P: pressure;

T: temperature;

T0: a reference temperature;

g: gravity acceleration;

ρ: reference density;

μ: dynamic viscosity;

α: coefficient of volumetric thermal expansion.


Model setup

We defined a 1m*1m square cavity filled with the fluid we defined above (reference pressure p=0, reference temperature T0=0). The four sides are no-sliping Wall. The right side is constant temperature (T{0=0), and the lower right wall and bottom edge 0.1m*0.1m are heat sources (T=1). Pressure Point Constraint with the four corners set to p=0.We use this model to simulate the natural convection of indoor enclosed space caused by heating (modified from Buoyancy Flow in Free Fluids by COMSOL Multiphysics).

Image: 400 pixels Image: 400 pixels

I examine a wide range of scenarios sequencing through different values of Ra using an Auxiliary sweep with continuation selected for the Ra parameter in the Stationary study step (Ra=[100 1e3 1e4 1e5 1e6]). With the increase of Ra, the importance of viscous force decreases.

For the convenience of numerical simulation, I defined a fluid that does not exist in reality, whose fluid parameters are shown in the following table (Pr: Prandtl number, the ratio of kinematic viscosity to thermal diffusivity.):

 Pr=(μCp)/k

μ: dynamic viscosity;

Cp: specific heat capacity of the fluid;

k = thermal diffusivity


Name Value
ρ 1
μ 1
k 1
Pr 0.71
Ra Changed as setting


So numerically (ignoring units) :

 Cp=Pr

body force in the y-direction for the momentum equation:

 F=(Ra/Pr)(T-T0)

Result and Discussion

We analyze our results from the aspects of velocity field and temperature field. Free convective flow is completely caused by temperature difference and the velocity field is generated by the temperature field and in turn affects the temperature field.

Time to stable state

We found that with the increase of Ra, the time of heat transfer to the relatively stable state within the space became shorter.

Ra=100 (in ~250ms)

Image: 500 pixels Image: 500 pixels

Ra=1e4 (in ~130ms)

Image: 500 pixels Image: 500 pixels

Ra=1e6 (in ~70ms)

Image: 500 pixels Image: 500 pixels

Velocity field

With the increase of Ra, the flow velocity of the fluid in stable state also increases significantly. Convection is near the boundary layer near the hot (cold) wall. Due to the influence of the heat source area on the left, the fluid flows upward along the left wall surface to the top of the square cavity. Due to the limitation of the height of the thermal insulation wall surface, the flow line turns to the right. Due to the influence of the right cold wall, the fluid flows down the right cold wall to the bottom of the square cavity. Due to the limitation of the height of the adiabatic wall, the streamline turns left and flows to the left hot wall. The whole flow field flows in a clockwise direction. When Ra= 100 or 1e3, there is a vortex at the center of the flow field. With the increase of Rayleigh number, the vortices in the flow field are divided into two symmetrical vortices and move towards the hot (cold) wall continuously.

Image: 1100 pixels

Temperature field

Convective heat transfer is the result of the combined action of fluid flow and heat conduction between fluid and molecules. When Ra = 100 or 10e3 (relatively low), the isotherm is uniformly arched outward with the heat source as the center, and the temperature gradient near the heat source is large. With the increase of Ra, the heat transfer changes from heat conduction to heat convection, and the isotherm gradually bends and deforms, forming a thin boundary layer near the cold wall and the heat source.

Image: 1100 pixels

Heat Transfer in Ice Stream Shut Down and Start (based on Joughin and Alley, 2011)

Stop and start of ice streams is dependent on:

  • Thinning (steepens basal temperature gradients)
  • Basal freezing (puts on the brakes for a system)
  • Thickening (build up traps geothermal heat)
  • Ice stream activation & the cycle is repeated

Watch the video below to walk through the steps of ice stream shut-down and speed-up, and see how this relates to the heat transfer equation introduced above.

Further reading:

Different theories exist for ice stream shut down, and here we present varying hypothesis. Please note that the behavior we model in our video is based upon the hypothesis of Joughin and Alley (2011).

Summary

This teaching module covered conduction, advection, and convection as heat transfer mechanisms. We can quantitatively characterize these processes in the heat transfer equation explained above. In the lava lamp, the transition from a conduction to advection dominated system occurs when temperature driven density reaches a critical low, allowing the ‘lava’ to rise through the surrounding fluid. The dimensionless Peclet Number can help us understand whether conduction or advection is driving most of the heat transfer. Though the same principles apply, behavior is a little different when considering the ice stream stop and start behavior. The regime shift between conduction and advection is greatly dependent on the phase change from ice to liquid water at the ice stream bed. The film of meltwater which allows advection to proceed is a product of conduction of heat from the bed. Convection, described above in the contexts of the lava lamp and heat circulation through an enclosed room, includes both conduction and advection as important mechanisms of heat transfer, but requires a spatial cycling of heat flow. The Rayleigh Number is a useful tool to evaluate the drivers and resistors of convection in a given system, and to determine whether convection will tend to grow or decay. Please check out the included links for more information or examples of heat transfer in the real world!