Skip to content
Ahmet Çelik
← Back to Study Guides

Final Formula Sheet

MECH302

Scope: Chapters 1, 2, 3, 5, 6, 7, 8, 9, 11. Excluded: Pool boiling. Purpose: Dense exam-use formula sheet. No explanations or worked examples.


Chapter 1 — Introduction

Heat Rate / Heat Flux

q=q˙A,q˙=qAq''=\frac{\dot{q}}{A},\qquad \dot{q}=q''A

Conduction

qx=kdTdxq_x''=-k\frac{dT}{dx} q=kT1T2Lq''=k\frac{T_1-T_2}{L} q˙=kAT1T2L\dot{q}=kA\frac{T_1-T_2}{L}

Convection

q=h(TsT)q''=h(T_s-T_\infty) q˙=hA(TsT)\dot{q}=hA(T_s-T_\infty)

Radiation

Eb=σTs4E_b=\sigma T_s^4 E=εσTs4E=\varepsilon\sigma T_s^4 q˙rad=εσA(Ts4Tsur4)\dot{q}_{\text{rad}}=\varepsilon\sigma A(T_s^4-T_{\text{sur}}^4) q˙abs=αGA\dot{q}_{\text{abs}}=\alpha GA σ=5.67×108 W/(m2K4)\sigma=5.67\times10^{-8}\ \text{W/(m}^2\cdot\text{K}^4)

Energy Balance

E˙inE˙out+E˙g=E˙st\dot{E}_{\text{in}}-\dot{E}_{\text{out}}+\dot{E}_g=\dot{E}_{st}

Steady, no generation:

E˙in=E˙out\dot{E}_{\text{in}}=\dot{E}_{\text{out}}

Chapter 2 — Introduction to Conduction

Thermal Diffusivity

α=kρcp\alpha=\frac{k}{\rho c_p}

Cartesian Heat Equation

2Tx2+2Ty2+2Tz2+q˙k=1αTt\frac{\partial^2 T}{\partial x^2}+\frac{\partial^2 T}{\partial y^2}+\frac{\partial^2 T}{\partial z^2}+\frac{\dot{q}'''}{k}=\frac{1}{\alpha}\frac{\partial T}{\partial t}

Cylindrical Heat Equation

1rr(rTr)+1r22Tϕ2+2Tz2+q˙k=1αTt\frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial T}{\partial r}\right) +\frac{1}{r^2}\frac{\partial^2 T}{\partial \phi^2} +\frac{\partial^2 T}{\partial z^2} +\frac{\dot{q}'''}{k} =\frac{1}{\alpha}\frac{\partial T}{\partial t}

Spherical Radial Heat Equation

1r2ddr(r2dTdr)+q˙k=1αTt\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dT}{dr}\right)+\frac{\dot{q}'''}{k}=\frac{1}{\alpha}\frac{\partial T}{\partial t}

Boundary Conditions

Specified temperature:

Ts=T0T_s=T_0

Specified heat flux:

  • k\frac{\partial T}{\partial n}=q_s”

Adiabatic/symmetry:

Tn=0\frac{\partial T}{\partial n}=0

Convection:

  • k\frac{\partial T}{\partial n}=h(T_s-T_\infty)

Radiation:

  • k\frac{\partial T}{\partial n}=\varepsilon\sigma(T_s^4-T_{\text{sur}}^4)

Interface:

TA=TBT_A=T_B
  • k_A\frac{\partial T_A}{\partial n}=-k_B\frac{\partial T_B}{\partial n}

Chapter 3 — One-Dimensional Steady Conduction

Thermal Resistance

q˙=ΔTRtotal\dot{q}=\frac{\Delta T}{R_{\text{total}}}

Plane wall:

Rcond=LkAR_{\text{cond}}=\frac{L}{kA}

Cylinder:

Rcond=ln(r2/r1)2πkLR_{\text{cond}}=\frac{\ln(r_2/r_1)}{2\pi kL}

Sphere:

Rcond=14πk(1r11r2)R_{\text{cond}}=\frac{1}{4\pi k}\left(\frac{1}{r_1}-\frac{1}{r_2}\right)

Convection:

Rconv=1hAR_{\text{conv}}=\frac{1}{hA}

Contact:

Rt,c=Rt,cAR_{t,c}=\frac{R_{t,c}''}{A}

Composite Plane Wall

q˙=T,1T,21h1A+LikiA+1h2A\dot{q}=\frac{T_{\infty,1}-T_{\infty,2}}{\frac{1}{h_1A}+\sum\frac{L_i}{k_iA}+\frac{1}{h_2A}}

Tube Wall with Convection

q˙=T,iT,o1hiAi+ln(ro/ri)2πkL+1hoAo\dot{q}=\frac{T_{\infty,i}-T_{\infty,o}}{\frac{1}{h_iA_i}+\frac{\ln(r_o/r_i)}{2\pi kL}+\frac{1}{h_oA_o}} Ai=2πriL,Ao=2πroLA_i=2\pi r_iL,\qquad A_o=2\pi r_oL

Overall Heat Transfer Coefficient

q˙=UAΔT\dot{q}=UA\Delta T 1UA=Rtotal\frac{1}{UA}=R_{\text{total}} 1UiAi=1hiAi+ln(ro/ri)2πkL+1hoAo\frac{1}{U_iA_i}=\frac{1}{h_iA_i}+\frac{\ln(r_o/r_i)}{2\pi kL}+\frac{1}{h_oA_o}

Uniform Heat Generation

Plane wall, half-thickness LL:

T(x)=Ts+q˙2k(L2x2)T(x)=T_s+\frac{\dot{q}}{2k}(L^2-x^2) Tmax=Ts+q˙L22kT_{\text{max}}=T_s+\frac{\dot{q}L^2}{2k} qs=q˙Lq_s''=\dot{q}L

Solid cylinder:

T(r)=Ts+q˙4k(ro2r2)T(r)=T_s+\frac{\dot{q}}{4k}(r_o^2-r^2) Tmax=Ts+q˙ro24kT_{\text{max}}=T_s+\frac{\dot{q}r_o^2}{4k} q=q˙πro2q'=\dot{q}\pi r_o^2

Fins

θ=TT,θb=TbT\theta=T-T_\infty,\qquad \theta_b=T_b-T_\infty m=hPkAcm=\sqrt{\frac{hP}{kA_c}} d2θdx2m2θ=0\frac{d^2\theta}{dx^2}-m^2\theta=0

Adiabatic tip:

qf=hPkAcθbtanh(mL)q_f=\sqrt{hPkA_c}\,\theta_b\tanh(mL)

Convective tip:

qf=Msinh(mL)+hmkcosh(mL)cosh(mL)+hmksinh(mL)q_f=M\frac{\sinh(mL)+\frac{h}{mk}\cosh(mL)}{\cosh(mL)+\frac{h}{mk}\sinh(mL)} M=hPkAcθbM=\sqrt{hPkA_c}\,\theta_b

Infinite fin:

qf=hPkAcθbq_f=\sqrt{hPkA_c}\,\theta_b

Corrected length:

Lc=L+AcPL_c=L+\frac{A_c}{P}

Fin efficiency:

ηf=qfhAfθb\eta_f=\frac{q_f}{hA_f\theta_b}

Overall surface efficiency:

ηo=1AfAt(1ηf)\eta_o=1-\frac{A_f}{A_t}(1-\eta_f) qt=ηohAtθbq_t=\eta_o hA_t\theta_b

Chapter 5 — Transient Conduction

Biot Number

Bi=hLckBi=\frac{hL_c}{k}

Lumped criterion:

Bi<0.1Bi<0.1

Lumped length:

Lc=VAsL_c=\frac{V}{A_s}

Lumped Capacitance

T(t)TTiT=exp(hAsρVcpt)\frac{T(t)-T_\infty}{T_i-T_\infty}=\exp\left(-\frac{hA_s}{\rho Vc_p}t\right) θθi=et/τ\frac{\theta}{\theta_i}=e^{-t/\tau} τ=ρVcphAs\tau=\frac{\rho Vc_p}{hA_s} q˙(t)=hAs[T(t)T]\dot{q}(t)=hA_s[T(t)-T_\infty] Q=ρVcp[TiT(t)]Q=\rho Vc_p[T_i-T(t)] Q0=ρVcp(TiT)Q_0=\rho Vc_p(T_i-T_\infty) QQ0=1θθi\frac{Q}{Q_0}=1-\frac{\theta}{\theta_i}

Fourier Number

Fo=αtL2plane wallFo=\frac{\alpha t}{L^2}\quad\text{plane wall} Fo=αtro2cylinder/sphereFo=\frac{\alpha t}{r_o^2}\quad\text{cylinder/sphere}

One-Term Approximation

Plane wall:

θθi=C1eζ12Focos(ζ1xL)\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}\cos\left(\zeta_1\frac{x}{L}\right)

Cylinder:

θθi=C1eζ12FoJ0(ζ1rro)\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}J_0\left(\zeta_1\frac{r}{r_o}\right)

Sphere:

θθi=C1eζ12Fosin(ζ1r/ro)ζ1r/ro\frac{\theta}{\theta_i}=C_1e^{-\zeta_1^2Fo}\frac{\sin(\zeta_1r/r_o)}{\zeta_1r/r_o}

Center temperature, all three:

θ0θi=C1eζ12Fo\frac{\theta_0}{\theta_i}=C_1e^{-\zeta_1^2Fo}

Multidimensional product:

(θθi)multiD=(θθi)1D\left(\frac{\theta}{\theta_i}\right)_{multi-D}=\prod\left(\frac{\theta}{\theta_i}\right)_{1D}

Chapter 6 — Introduction to Convection

q=h(TsT)q''=h(T_s-T_\infty) q˙=hA(TsT)\dot{q}=hA(T_s-T_\infty) Nu=hLckfNu=\frac{hL_c}{k_f} h=NukfLch=\frac{Nu\,k_f}{L_c} Re=VLcν=ρVLcμRe=\frac{VL_c}{\nu}=\frac{\rho VL_c}{\mu} Pr=να=μcpkfPr=\frac{\nu}{\alpha}=\frac{\mu c_p}{k_f} Gr=gβΔTLc3ν2Gr=\frac{g\beta\Delta TL_c^3}{\nu^2} Ra=GrPrRa=GrPr Tf=Ts+T2T_f=\frac{T_s+T_\infty}{2} hˉ=1L0Lhxdx\bar h=\frac{1}{L}\int_0^Lh_x\,dx

Chapter 7 — External Flow

Flat Plate

Rex=Uxν,ReL=ULνRe_x=\frac{U_\infty x}{\nu},\qquad Re_L=\frac{U_\infty L}{\nu} Rex,c5×105Re_{x,c}\approx5\times10^5 xc=Rex,cνUx_c=\frac{Re_{x,c}\nu}{U_\infty}

Laminar local:

Nux=0.332Rex1/2Pr1/3Nu_x=0.332Re_x^{1/2}Pr^{1/3}

Laminar average:

NuL=0.664ReL1/2Pr1/3\overline{Nu}_L=0.664Re_L^{1/2}Pr^{1/3}

Turbulent local:

Nux=0.0296Rex4/5Pr1/3Nu_x=0.0296Re_x^{4/5}Pr^{1/3}

Turbulent average from leading edge:

NuL=0.037ReL4/5Pr1/3\overline{Nu}_L=0.037Re_L^{4/5}Pr^{1/3}

Mixed average:

NuL=(0.037ReL4/5871)Pr1/3\overline{Nu}_L=(0.037Re_L^{4/5}-871)Pr^{1/3} h=NukfLch=\frac{Nu\,k_f}{L_c} q˙=hˉA(TsT)\dot{q}=\bar hA(T_s-T_\infty)

Cylinder

ReD=VDνRe_D=\frac{VD}{\nu}

Hilpert:

NuD=CReDmPr1/3\overline{Nu}_D=CRe_D^mPr^{1/3} hˉ=NuDkfD\bar h=\frac{\overline{Nu}_Dk_f}{D} q˙=hˉπDL(TsT)\dot{q}=\bar h\pi DL(T_s-T_\infty) q=hˉπD(TsT)q'=\bar h\pi D(T_s-T_\infty)

Chapter 8 — Internal Flow

Ac=πD24,As=πDLA_c=\frac{\pi D^2}{4},\qquad A_s=\pi DL um=m˙ρAcu_m=\frac{\dot{m}}{\rho A_c} ReD=ρumDμ=umDν=4m˙πDμRe_D=\frac{\rho u_mD}{\mu}=\frac{u_mD}{\nu}=\frac{4\dot{m}}{\pi D\mu} Dh=4AcPD_h=\frac{4A_c}{P}

Laminar/turbulent cutoff:

ReD2300Re_D\approx2300

Entrance lengths:

xfd,hD0.05ReDlaminar\frac{x_{fd,h}}{D}\approx0.05Re_D\quad\text{laminar} xfd,tD0.05ReDPrlaminar\frac{x_{fd,t}}{D}\approx0.05Re_DPr\quad\text{laminar} xfdD10turbulent estimate\frac{x_{fd}}{D}\approx10\quad\text{turbulent estimate}

Energy balance:

q˙=m˙cp(Tm,oTm,i)\dot{q}=\dot{m} c_p(T_{m,o}-T_{m,i})

Constant heat flux:

q˙=qsπDL\dot{q}=q_s''\pi DL Tm(x)=Tm,i+qsπDxm˙cpT_m(x)=T_{m,i}+\frac{q_s''\pi Dx}{\dot{m} c_p} TsTm=qshT_s-T_m=\frac{q_s''}{h}

Constant wall temperature:

TsTm,oTsTm,i=exp(hˉAsm˙cp)\frac{T_s-T_{m,o}}{T_s-T_{m,i}}=\exp\left(-\frac{\bar hA_s}{\dot{m} c_p}\right) q˙=hˉAsΔTlm\dot{q}=\bar hA_s\Delta T_{\text{lm}} ΔTlm=(TsTm,i)(TsTm,o)ln[(TsTm,i)/(TsTm,o)]\Delta T_{\text{lm}}=\frac{(T_s-T_{m,i})-(T_s-T_{m,o})}{\ln[(T_s-T_{m,i})/(T_s-T_{m,o})]}

Laminar fully developed:

NuD=3.66constant TsNu_D=3.66\quad\text{constant }T_s NuD=4.36constant qsNu_D=4.36\quad\text{constant }q_s''

Laminar developing:

NuD=3.66+0.0668(D/L)ReDPr1+0.04[(D/L)ReDPr]2/3\overline{Nu}_D=3.66+\frac{0.0668(D/L)Re_DPr}{1+0.04[(D/L)Re_DPr]^{2/3}}

Turbulent Dittus–Boelter:

NuD=0.023ReD0.8PrnNu_D=0.023Re_D^{0.8}Pr^n n=0.4 heating fluid,n=0.3 cooling fluidn=0.4\ \text{heating fluid},\qquad n=0.3\ \text{cooling fluid}

Chapter 9 — Free Convection

Tf=Ts+T2T_f=\frac{T_s+T_\infty}{2}

Ideal gas:

β=1Tf\beta=\frac{1}{T_f} GrL=gβTsTL3ν2Gr_L=\frac{g\beta|T_s-T_\infty|L^3}{\nu^2} RaL=GrLPrRa_L=Gr_LPr Nu=hˉLckfNu=\frac{\bar hL_c}{k_f}

Vertical plate:

NuL=[0.825+0.387RaL1/6[1+(0.492/Pr)9/16]8/27]2\overline{Nu}_L=\left[0.825+\frac{0.387Ra_L^{1/6}}{[1+(0.492/Pr)^{9/16}]^{8/27}}\right]^2

Vertical cylinder as vertical plate if:

DL35GrL1/4\frac{D}{L}\gtrsim\frac{35}{Gr_L^{1/4}}

Horizontal cylinder:

NuD=[0.60+0.387RaD1/6[1+(0.559/Pr)9/16]8/27]2\overline{Nu}_D=\left[0.60+\frac{0.387Ra_D^{1/6}}{[1+(0.559/Pr)^{9/16}]^{8/27}}\right]^2

Horizontal plate characteristic length:

Lc=AsPL_c=\frac{A_s}{P}

Hot up / cold down:

NuL=0.54RaL1/4(104RaL107)\overline{Nu}_L=0.54Ra_L^{1/4}\quad(10^4\lesssim Ra_L\lesssim10^7) NuL=0.15RaL1/3(107RaL1011)\overline{Nu}_L=0.15Ra_L^{1/3}\quad(10^7\lesssim Ra_L\lesssim10^{11})

Hot down / cold up:

NuL=0.27RaL1/4\overline{Nu}_L=0.27Ra_L^{1/4}

Convection:

q˙conv=hˉAs(TsT)\dot{q}_{\text{conv}}=\bar hA_s(T_s-T_\infty)

Radiation:

q˙rad=εσAs(Ts4Tsur4)\dot{q}_{\text{rad}}=\varepsilon\sigma A_s(T_s^4-T_{\text{sur}}^4)

Combined:

q˙total=q˙conv+q˙rad\dot{q}_{\text{total}}=\dot{q}_{\text{conv}}+\dot{q}_{\text{rad}}

Chapter 11 — Heat Exchangers

Ch=m˙hcp,h,Cc=m˙ccp,cC_h=\dot{m}_h c_{p,h},\qquad C_c=\dot{m}_c c_{p,c} q˙=Ch(Th,iTh,o)=Cc(Tc,oTc,i)\dot{q}=C_h(T_{h,i}-T_{h,o})=C_c(T_{c,o}-T_{c,i}) 1UA=1hiAi+Rf,i+ln(ro/ri)2πkL+Rf,o+1hoAo\frac{1}{UA}=\frac{1}{h_iA_i}+R_{f,i}+\frac{\ln(r_o/r_i)}{2\pi kL}+R_{f,o}+\frac{1}{h_oA_o} ΔTlm=ΔT1ΔT2ln(ΔT1/ΔT2)\Delta T_{\text{lm}}=\frac{\Delta T_1-\Delta T_2}{\ln(\Delta T_1/\Delta T_2)} q˙=UAFΔTlm\dot{q}=UAF\Delta T_{\text{lm}}

Parallel:

ΔT1=Th,iTc,i\Delta T_1=T_{h,i}-T_{c,i} ΔT2=Th,oTc,o\Delta T_2=T_{h,o}-T_{c,o}

Counterflow:

ΔT1=Th,iTc,o\Delta T_1=T_{h,i}-T_{c,o} ΔT2=Th,oTc,i\Delta T_2=T_{h,o}-T_{c,i} Cmin=min(Ch,Cc),Cmax=max(Ch,Cc)C_{\text{min}}=\min(C_h,C_c),\qquad C_{\text{max}}=\max(C_h,C_c) Cr=CminCmaxC_r=\frac{C_{\text{min}}}{C_{\text{max}}} q˙max=Cmin(Th,iTc,i)\dot{q}_{\text{max}}=C_{\text{min}}(T_{h,i}-T_{c,i}) ε=q˙q˙max\varepsilon=\frac{\dot{q}}{\dot{q}_{\text{max}}} NTU=UACminNTU=\frac{UA}{C_{\text{min}}} q˙=εCmin(Th,iTc,i)\dot{q}=\varepsilon C_{\text{min}}(T_{h,i}-T_{c,i})

Parallel effectiveness:

ε=1eNTU(1+Cr)1+Cr\varepsilon=\frac{1-e^{-NTU(1+C_r)}}{1+C_r}

Counterflow, Cr<1C_r<1:

ε=1eNTU(1Cr)1CreNTU(1Cr)\varepsilon=\frac{1-e^{-NTU(1-C_r)}}{1-C_re^{-NTU(1-C_r)}}

Counterflow, Cr=1C_r=1:

ε=NTU1+NTU\varepsilon=\frac{NTU}{1+NTU}

Phase change, Cr=0C_r=0:

ε=1eNTU\varepsilon=1-e^{-NTU}