Heat Equation Excel Solver

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 59

Created by Ron Hugo, 26 SEP 2015

Updated 22 OCT 2018 to correct for an error in radiation / convection boundary condition.

INTERIOR NODE

INTERIOR NODE with GENERATION

Generation Rate [W/m^3] q_dot =


Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =

CONVECTIVE BOUNDARY 1

Conv. Heat Transfer Coefficient [W/m^2 K] h_1 =


Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Biot Number Bi = h_1 Dx / k #DIV/0!


Freestream Temperature (C) T_conv =
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right convective surface, or inside left surface

T_conv, h_1
k, q_dot

Left convective surface, or inside right surface

T_conv, h_1

k, q_dot
Top convective surface, or inside bottom surface

T_conv, h_1

k, q_dot

Bottom convective surface, or inside top surface

k, q_dot

T_conv, h_1

Top right convective corner

T_conv, h_1

T_conv, h_1
k, q_dot
Top left convective corner

T_conv, h_1

k, q_dot
T_conv, h_1
Bottom right convective corner
T_conv, h_1
k, q_dot

T_conv, h_1

Bottom left convective corner


T_conv, h_1
k, q_dot
k, q_dot

T_conv, h_1

Interior top left convective corner

k, q_dot

T_conv, h_1

Interior top right convective corner

k, q_dot

T_conv, h_1

Interior bottom left convective corner

T_conv, h_1

k, q_dot

Interior bottom right convective corner

T_conv, h_1

k, q_dot

CONVECTIVE BOUNDARY 2

Conv. Heat Transfer Coefficient [W/m^2 K] h_2 =


Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Biot Number Bi = h_2 Dx / k #DIV/0!


Freestream Temperature (C) T_conv =
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right convective surface, or inside left surface

T_conv, h_2
k, q_dot

Left convective surface, or inside right surface

T_conv, h_2

k, q_dot

Top convective surface, or inside bottom surface

T_conv, h_2

k, q_dot

Bottom convective surface, or inside top surface

k, q_dot

T_conv, h_2

Top right convective corner

T_conv, h_2
T_conv, h_2

T_conv, h_2
k, q_dot
Top left convective corner

T_conv, h_2

k, q_dot
T_conv, h_2
Bottom right convective corner
T_conv, h_2
k, q_dot

T_conv, h_2

Bottom left convective corner


T_conv, h_2
k, q_dot

T_conv, h_2

Interior top left convective corner

k, q_dot

T_conv, h_2

Interior top right convective corner

k, q_dot
T_conv, h_2

Interior bottom left convective corner

T_conv, h T_conv, h_2

k, q_dot

Interior bottom right convective corner

T_conv, h_2

k, q_dot

INSULATED BOUNDARY

Thermal Conductivity [W/m K] k=


Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right insulated surface, or inside left surface

k, q_dot

Left insulated surface, or inside right surface

k, q_dot
Top insulated surface, or inside bottom surface

k, q_dot
Bottom insulated surface, or inside top surface

k, q_dot

Top right insulated corner

k, q_dot
Top left insulated corner

k, q_dot
Bottom right insulated corner

k, q_dot

Bottom left insulated corner


k, q_dot

Interior top left insulated corner

k, q_dot

Interior top right insulated corner

k, q_dot

Interior bottom left insulated corner

k, q_dot

Interior bottom right insulated corner

k, q_dot

CONSTANT HEAT FLUX WITH CONVECTIVE BOUNDARY


Conv. Heat Transfer Coefficient [W/m^2 K] h_3 =
Surface Heat Flux [W/m^2] q'' =
Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Heat Flux Premultiplier Dx / (2 k) = #DIV/0!


Biot Number Bi = h_3 Dx / k #DIV/0!
Freestream Temperature (C) T_conv =
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right outer surface, or inside left surface

k, q_dot
T_conv, h_3
q''

Left outer surface, or inside right surface

k, q_dot
T_conv, h_3
q''

Top outer surface, or inside bottom surface

T_conv, h_3

q''

k, q_dot

Bottom outer surface, or inside top surface

k, q_dot
T_conv, h

q''
T_conv, h_3

Top right corner

T_conv, h_3

q''

T_conv, h_3

k, q_dot

Top left corner

T_conv, h_3

q''

T_conv, h_3
k, q_dot

Bottom right corner

k, q_dot

T_conv, h_3

q''

T_conv, h_3

Bottom left corner

k, q_dot
T_conv, h_3

q''
q''

T_conv, h_3

Interior top left corner

k, q_dot

q''
T_conv, h_3

Interior top right corner

k, q_dot

q''
T_conv, h_3

Interior bottom left corner

q''
T_conv, h_3

k, q_dot

Interior bottom right corner

q''

T_conv, h_3
k, q_dot

CONSTANT HEAT FLUX BOUNDARY

Surface Heat Flux [W/m^2] q'' =


Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =
Heat Flux Premultiplier Dx / (2 k) = #DIV/0!
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right outer surface, or inside left surface

k, q_dot

q''

Left outer surface, or inside right surface

k, q_dot

q''

Top outer surface, or inside bottom surface

q''

k, q_dot

Bottom outer surface, or inside top surface

k, q_dot

q''
q''

Top right corner

q''

k, q_dot

Top left corner

q''

k, q_dot

Bottom right corner

k, q_dot

q''

Bottom left corner

k, q_dot

q''

Interior top left corner

k, q_dot

q''
Interior top right corner

k, q_dot

q''

Interior bottom left corner

q''

k, q_dot

Interior bottom right corner

q''

k, q_dot

COMPOSITE SOLIDS 1

Internal Thermal Conductivity [W/m K] k_int =


External Thermal Conductivity [W/m K] k_ext =
Internal Generation Rate [W/m^3] q_dot_int =
External Generation Rate [W/m^3] q_dot_ext =
Mesh Size [m] Dx = Dy =
Conv. Heat Transfer Coefficient [W/m^2 K] h=
Freestream Temperature (C) T_conv =

k_ext / (2*(k_int + k_ext)) = #DIV/0!


k_int / (2*(k_int + k_ext)) = #DIV/0!
Dx*Dx / (4*(k_int + k_ext)) = #DIV/0!

k_ext / (k_int + 3 k_ext) = #DIV/0!


k_int / (2*(k_int + 3 k_ext)) = #DIV/0!
Dx*Dx / (4*(k_int + 3 k_ext)) = #DIV/0!

k_int/(2*(k_int + k_ext + h Dx)) = #DIV/0!


k_ext/(2*(k_int + k_ext + h Dx)) = #DIV/0!
h Dx T_conv / (k_int + k_ext + h Dx) = #DIV/0!
Dx*Dx / (4*(k_int + k_ext + h Dx)) = #DIV/0!

Top composite surface

k_ext,
q_dot_ext

k_int,
q_dot_int

Right composite surface


k_ext,
k_int, q_dot_ext
q_dot_int

Left composite surface


k_ext,
q_dot_ext k_int,
q_dot_int

Bottom composite surface

k_int,
q_dot_int

k_ext,
q_dot_ext

Top right composite corner


k_int k_ext
q_dot_int q_dot_ext

Top left composite corner

k_ext
q_dot_ext

k_int
q_dot_int

Bottom left composite corner

k_ext
q_dot_ext

k_int
q_dot_int

Bottom right composite corner


k_int
q_dot_int

k_ext
q_dot_ext

Top comp. conv. surface, internal material to the right

T_conv, h

k_ext k_int
q_dot_ext q_dot_int

Top comp. conv. surface, internal material to the left

T_conv, h

k_int k_ext
q_dot_int q_dot_ext
Left comp. conv. surface, internal material at bottom
k_ext
q_dot_ext

T_conv, h k_int
q_dot_int
Left comp. conv. surface, internal material at top

k_int
q_dot_int

T_conv, h
k_ext
q_dot_ext

Right comp. conv. surface, internal material at bottom


k_ext
q_dot_ext

k_int
q_dot_int T_conv, h

Right comp. conv. surface, internal material at top

k_int
q_dot_int

T_conv, h
k_ext
q_dot_ext

Bottom comp. conv. surface, internal material to the left

k_int k_ext
q_dot_int q_dot_ext
T_conv, h

Bottom comp. conv. surface, internal material to the right

k_ext k_int
q_dot_ext q_dot_int

T_conv, h

COMPOSITE SOLIDS 2

Internal Thermal Conductivity [W/m K] k_int =


External Thermal Conductivity [W/m K] k_ext =
Internal Generation Rate [W/m^3] q_dot_int =
External Generation Rate [W/m^3] q_dot_ext =
Mesh Size [m] Dx = Dy =
Conv. Heat Transfer Coefficient [W/m^2 K] h=
Freestream Temperature (C) T_conv =

k_ext / (2*(k_int + k_ext)) = #DIV/0!


k_int / (2*(k_int + k_ext)) = #DIV/0!
Dx*Dx / (4*(k_int + k_ext)) = #DIV/0!

k_ext / (k_int + 3 k_ext) = #DIV/0!


k_int / (2*(k_int + 3 k_ext)) = #DIV/0!
Dx*Dx / (4*(k_int + 3 k_ext)) = #DIV/0!

k_int/(2*(k_int + k_ext + h Dx)) = #DIV/0!


k_ext/(2*(k_int + k_ext + h Dx)) = #DIV/0!
h Dx T_conv / (k_int + k_ext + h Dx) = #DIV/0!
Dx*Dx / (4*(k_int + k_ext + h Dx)) = #DIV/0!

Top composite surface

k_ext,
q_dot_ext

k_int,
q_dot_int
Right composite surface
k_ext,
k_int, q_dot_ext
q_dot_int

Left composite surface


k_ext,
q_dot_ext k_int,
q_dot_int

Bottom composite surface

k_int,
q_dot_int

k_ext,
q_dot_ext

Top right composite corner

k_int k_ext
q_dot_int q_dot_ext
Top left composite corner

k_ext
q_dot_ext

k_int
q_dot_int

Bottom left composite corner


k_ext
q_dot_ext

k_int
q_dot_int

Bottom right composite corner


k_int
q_dot_int

k_ext
q_dot_ext

Top comp. conv. surface, internal material to the right

T_conv, h

k_ext k_int
q_dot_ext q_dot_int
Top comp. conv. surface, internal material to the left

T_conv, h

k_int k_ext
q_dot_int q_dot_ext

Left comp. conv. surface, internal material at bottom


k_ext
q_dot_ext

T_conv, h k_int
q_dot_int
Left comp. conv. surface, internal material at top

k_int
q_dot_int
k_int
q_dot_int

T_conv, h
k_ext
q_dot_ext

Right comp. conv. surface, internal material at bottom


k_ext
q_dot_ext

k_int
q_dot_int T_conv, h

Right comp. conv. surface, internal material at top

k_int
q_dot_int

T_conv, h
k_ext
q_dot_ext

Bottom comp. conv. surface, internal material to the left

k_int k_ext
q_dot_int q_dot_ext

T_conv, h

Bottom comp. conv. surface, internal material to the right

k_ext k_int
q_dot_ext q_dot_int

T_conv, h
RADIATION AND CONVECTION BOUNDARY 1

Surface Emissivity e=
Surroundings Temperature [K] T_surr =
Stefan-Boltzmann Constant [W/m^2 K^4] s = 5.669E-08
Conv. Heat Transfer Coefficient [W/m^2 K] h=
Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Biot Number Bi = h Dx / k #DIV/0!


Freestream Temperature (C) T_conv =
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!
2k + hDx 0

Right convective surface, or inside left surface

T_surr
k, q_rad
e,
q_dot

T_conv, h
Left convective surface, or inside right surface
T_surr

q_rad

k,
e,
q_dot
T_conv, h
Top convective surface, or inside bottom surface
T_surr
q_rad
T_conv, h

k, e, q_dot

Bottom convective surface, or inside top surface

k, e, q_dot
k, e, q_dot

T_conv, h
q_rad
T_surr

Top right convective corner


T_surr
T_conv, h q_rad

k, e, q_dot T_conv, h

Top left convective corner


T_surr
q_rad T_conv, h

k, e, q_dot
T_conv, h
Bottom right convective corner
T_conv, h
k, e, q_dot

q_rad
T_conv, h
T_surr
Bottom left convective corner
T_conv, h
k, e, q_dot

q_rad
T_conv, h
T_surr

Interior top left convective corner

k, e,
q_dot
q_rad

T_conv, h T_surr
Interior top right convective corner

k, e,
q_dot

q_rad

T_surr T_conv, h
Interior bottom left convective corner
T_surr
q_rad

T_conv, h
k, e,
q_dot

Interior bottom right convective corner


T_surr
q_rad k, e,
q_dot

T_conv, h

RADIATION AND CONVECTION BOUNDARY 2

Surface Emissivity e=
Surroundings Temperature [K] T_surr =
Stefan-Boltzmann Constant [W/m^2 K^4] s = 5.669E-08
Conv. Heat Transfer Coefficient [W/m^2 K] h=
Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Biot Number Bi = h Dx / k #DIV/0!


Freestream Temperature (C) T_conv =
Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!
2k + hDx 0

Right convective surface, or inside left surface

T_surr
k, q_rad
e,
q_dot

T_conv, h
Left convective surface, or inside right surface
T_surr

q_rad

k,
e,
q_dot
T_conv, h
Top convective surface, or inside bottom surface
T_surr
q_rad
T_conv, h

k, e, q_dot

Bottom convective surface, or inside top surface

k, e, q_dot

T_conv, h
q_rad
T_surr

Top right convective corner


T_surr
T_conv, h q_rad

k, e, q_dot T_conv, h
k, e, q_dot

Top left convective corner


T_surr
q_rad T_conv, h

k, e, q_dot
T_conv, h
Bottom right convective corner
T_conv, h
k, e, q_dot

q_rad
T_conv, h
T_surr
Bottom left convective corner
T_conv, h
k, e, q_dot

q_rad
T_conv, h
T_surr

Interior top left convective corner

k, e,
q_dot

q_rad

T_conv, h T_surr
Interior top right convective corner

k, e,
q_dot

q_rad

T_surr T_conv, h
Interior bottom left convective corner
T_surr
q_rad
q_rad

T_conv, h
k, e,
q_dot

Interior bottom right convective corner


T_surr
q_rad k, e,
q_dot

T_conv, h

RADIATION BOUNDARY

Surface Emissivity e=
Surroundings Temperature [K] T_surr =
Stefan-Boltzmann Constant [W/m^2 K^4] s = 5.669E-08
Thermal Conductivity [W/m K] k=
Mesh Size [m] Dx = Dy =
Generation Rate [W/m^3] q_dot =

Freestream Temperature (C) T_conv =


Generation Term Premultiplier Dx^2 / (2 k) = #DIV/0!

Right surface, or inside left surface

k,
e, q_rad T_surr
q_dot

Left surface, or inside right surface

q_rad
T_surr k,
e,
q_dot
k,
e,
q_dot

Top convective surface, or inside bottom surface


T_surr

q_rad

k, e, q_dot

Bottom convective surface, or inside top surface

k, e, q_dot

q_rad

T_surr

Top right convective corner

T_surr

q_rad

k, e, q_dot

Top left convective corner

T_surr q_rad

k, e, q_dot

Bottom right convective corner

k, e, q_dot

q_rad

T_surr
Bottom left convective corner
k, e, q_dot

q_rad

T_surr

Interior top left convective corner

k, e,
q_dot

q_rad

T_surr
Interior top right convective corner

k, e,
q_dot

q_rad

T_surr
Interior bottom left convective corner
T_surr
q_rad

k, e,
q_dot

Interior bottom right convective corner


T_surr
q_rad k, e,
q_dot
ection boundary condition.

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
0

0
0

Err:522

Err:522

Err:522

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!
0

Err:522
Err:522

Err:522

Err:522

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
#DIV/0!

#DIV/0!

#DIV/0!

#DIV/0!
25 25 25 25 25 25
25 25 25 25 25 25
25 25 25 25 25 25
25 25 25 25 25 25
25 25 25 25 25 25

You might also like