🔏
brainless
  • What is this?
  • Responsibility
  • Changelog
  • meta
    • Sharing
      • Inspirations
      • Workflows
      • Social Media
    • Geography
      • Life
      • Death
        • Family Death
    • Research
      • Project Index
      • 3D Printing
      • Photogrammetry
      • Drone Building
    • External Websites
    • [unreleased]
      • [Template]
      • [TEMP] CL V0.0.5 or whatever
      • Skincare
      • Travel
      • Working and Staying Busy
      • Stride
      • Funeral Playlist
      • Notes and Ideas
      • Boredom
      • Four Noble Truths of "Thermo"
      • Respect
      • Work
  • STE[A]M
    • [guide]
    • Science
      • Materials Modeling
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
          • index
          • Carbon Nanotubes
    • Technology
      • Computer Science
        • Commands
      • Photogrammetry
      • Quantum Computing
      • Computers
      • Programs
        • Matlab and Octave
        • Audacity
        • Google Chrome
          • Websites
            • Google Suite Sites
            • Github
              • Version Control
            • Product Hunt
            • Twitter
            • Youtube
              • Channels
            • Vimeo
          • Extensions
            • Dark Reader
            • Vimium
        • [miscellaneous]
          • Octave
          • PureRef
          • git
          • gnu stow
          • mermaid.js
        • Excel
        • Blender
        • LaTeX
        • Sublime: Text Editor
        • Spotify
        • VLC Media Player
      • Android and iOS
      • Operating Systems
        • macOS
          • mackup
        • Unix
          • folder structure
        • Windows
          • App Installation
          • Meshroom
          • Corsair Utility Engine
      • 3D Printing
    • Engineering
      • Accreditation
        • Fundamentals of Engineering
        • Professional Engineering
      • Continuum Mechanics
        • Fluid Mechanics
          • Incompressible Flow
            • corona final
            • indexhw4
            • indexhw3
            • index
            • hw2
            • hw1
          • Syllabus Description
          • Lecture Slides
          • Student Notes
            • Dynamic or Kinematic Viscosity
          • Assignments
            • Homeworks
              • Homework 1
              • Homework 2
              • Homework 3
              • Homework 4
              • Homework 5
            • Vortex Project
        • Solid Mechanics
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
        • Incompresible Flow
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
      • Experimental Mechanics
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
      • Finite Element Methods
        • Intro to Finite Elements
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
        • Fundamentals of FEM
          • Syllabus Description
            • index
          • Lecture Slides
          • Student Notes
          • Assignments
            • Project
              • index
              • Untitled
            • Homework 1
            • Homework 4
            • index
      • Heat Transfer
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
          • homework
            • hw10
            • q9
            • q8
            • q7
            • hw7
            • hw6
            • q5
            • q3
            • 1 ec
          • Discussions
            • d11
            • d10
            • d9
            • d8
            • d6
            • d4
            • d3
          • Project Notes
      • Machine Dynamics
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
    • Art
      • Color Theory
      • Origami
        • FolderMath
          • Surveying Origami Math
          • Represent a Folded Object
          • Creating a Crease Pattern
          • Making the Folds
          • Simulating Folding Origami
          • List of Resources
            • Codes
            • Papers, Programs, and Inspirations
    • Mathematics
      • Complex Numbers
        • What is i^i?
      • Analytic Hierarchy Process
      • Probability
      • Conway's Game of Life
      • Metallic Numbers
      • Cauchy's Formula for Repeated Integration
      • Wavelet Transform
      • Laplace Tidal Equation
      • Alternating Summation of Ones
      • Constants
      • Bad Maths
      • Calculus
        • Syllabus Description
        • Miscellaneous
  • Thoughts
    • Marksmanship
      • Archery
    • Schooling
    • ...and Ideas?
      • Perceived Time and Learning
      • Content Comprehension
    • Comics and Games
      • Rubik's Cube
      • Dungeons and Dragons
      • Beyond-All-Reason
      • Sekiro: Shadows Die Twice
      • Super Smash Bros
        • Project M
        • Project +
      • League of Legends
      • Satisfactory
    • Literature and Art
      • Books
      • Reading is Hard
      • Various Words and Phrases
      • Poems
      • Interviews
      • Quotes
        • Phrases
      • Jokes
      • ASCII Art
    • Shows and Films
      • Cowboy Bebop
      • My Hero Academia
      • Sword of the Stranger
    • Working and Life Balance
  • Projects
Powered by GitBook
On this page
  1. STE[A]M
  2. Engineering
  3. Heat Transfer
  4. Assignments
  5. homework

hw10

PrevioushomeworkNextq9

Last updated 4 years ago

Chris Nkinthorn $\texttt{20190805}$

Prompt : (100 Points) Consider a 2-dimensional (2D) $[x, y]$, steady, incompressible, buoyant flow :

​ a. Write down the governing Navier-Stokes equations

​ b. Express the boundary layer equations in terms of Prandtl $\text{Pr}$, Reynolds $\text{Re}$, Grashof $\text{Gr}$, and Eckert $\text{Ec}$ numbers

​ c. Develop a corresponding set of boundary layer equations [momentum, energy], where the momentum $[\delta]$ and thermal thicknesses $\delta_T$ are not the same.

​ d. If Ec is small and if the convection terms are assumed of order the conduction terms, infer a relation for T as a function of other non-dimensional parameters.

Solution:

Part a: Write down the governing Navier-Stokes (N-S) equations, using the given assumptions for 2D $[x, y]$ flow: steady or constant with respect to time and incompressible or constant with respect to density. We include the buoyant term in the N-S equation as:

ρ[V⃗∙∇]V⃗=−∇p+μ∇2V⃗+ρgβΔT⏞Bouyant.\rho[\vec{V} \bullet \nabla] \vec{V} = -\nabla p+\mu \nabla^{2} \vec{V}+ \overbrace{\rho g \beta \Delta T}^\texttt{Bouyant}.ρ[V∙∇]V=−∇p+μ∇2V+ρgβΔT​Bouyant​.

Momentum conservation is dependent on three factors: the pressure, viscosity, and the aforementioned buoyancy. Buoyancy is a function of the gravitational acceleration, the local pressure, coupled with the Laplacian of the temperature distribution, per Fourier’s law.

Additionally to derive the boundary layer equation, the energy conservation expression is required. By the inclusion of specific heat capacity, the conservation for energy expression is seen to be a factor the material by the thermal conductivity, $k$, and viscous heat dissipation, $\Phi$, as:

ρCp[V⃗∙∇]T=k∇2T+Φ.% viscous dissipation fuction \rho C_{p}[\vec{V} \bullet \nabla] T=k \nabla^{2} T+\Phi.ρCp​[V∙∇]T=k∇2T+Φ.

Part b: To express the boundary layer (BL) equations for momentum and energy in terms of Prandtl $\text{Pr}$, Reynolds $\text{Re}$, Grashof $\text{Gr}$, and Eckert $\text{Ec}$ numbers, we examine each BL individually. First, the momentum BL which depends on $\text{Gr}$ is premultiplied by $\frac{U_\infty^2}{L}$:

Similarly, expressing the thermal energy boundary layer involves the same process of premultiplication on each side of the equation.

Here, the coupling to the momentum boundary layer comes from the heat capacity of the fluid, the temperature and velocity distributions, as well as the Prandtl number. The Eckert number is found by definition as the ratio of advective transfer and the total possible heat dissipation.

Part c: To develop a corresponding set of boundary layer equations [momentum, energy], where the momentum $[\delta]$ and thermal thicknesses $[\delta_T]$ are not the same, order of magnitude analysis is employed. Let variation in the $x$ direction correspond to the free stream velocity $u$ normalized to $1$ and $\delta$ correspond to the high of the boundary layer on length scale much smaller than the free stream velocity.

The reciprocal second positional derivative in the normal direction contributes much more than the component in the flow direction. This means that the expression simplifies:

Similarly, analysis of the variation in the flow component normal to the boundary surface yields:

This leads to the conclusion that the pressure variation in the normal direction is very much negligible. At least in so far in comparison to the variation along the streamwise direction. Having investigated the vector components of momentum, additional order-of-magnitude analysis of the scalar energy expression yields:

Again, the relative contribution of second derivatives is found and the streamwise direction neglected, so the simplification produces:

Part d: Additionally, if $\text{Ec}$ is small and the convection terms are assumed to be on the order of the conduction terms, to infer a relation for $\delta_{T}$ as a function of other non-dimensional parameters, we reexamine the relationship previously established, after order of magnitude simplification:

so that application of ideal gas law in the far stream (noted with the subscript $\infty$) allows for correlation between the specific heat capacities for constant volume $c_v$ and constant pressure $c_p$. Algebraically, we can correlate with respect the ratio of these heat capacities, $\gamma$ .

Additionally, find the definition of the material parameter of the speed of sound of the fluid as:

These, in conjunction with the order of magnitude analysis on the Eckert number which is the ratio of convective heat transfer to the possible diffusive transport. The term $\Delta T$, represents local deviation from the reference temperature, $\Delta T_{ref}$.

Applying this to the form given in viscous boundary layer expression $\delta$ to the thermal BL $\delta_T$, produces:

ρ[V⃗f∙∇]V⃗=−∇p+μ∇2V⃗+ρgβΔTρU∞2L[V⃗∙∇]V⃗=−ρU∞2L∇p+μU∞L2∇2V⃗+ρgβΔT[ΔT]refρ[V⃗∙∇]V⃗=−∇p+1Re⁡∇2V⃗+GrRe⁡2\rho[\vec{V}f \bullet \nabla] \vec{V}=-\nabla p+\mu \nabla^{2} \vec{V}+\rho g \beta \Delta T\\ \rho \frac{U_{\infty}^{2}}{L}[\vec{V} \bullet \nabla] \vec{V}=-\rho \frac{U_{\infty}^{2}}{L} \nabla p+\mu \frac{U_{\infty}}{L^{2}} \nabla^{2} \vec{V}+\rho g \beta \Delta T[\Delta T]_{r e f}\\ \boxed{ \rho[\vec{V} \bullet \nabla] \vec{V}=-\nabla p+\frac{1}{\operatorname{Re}} \nabla^{2} \vec{V}+\frac{G r}{\operatorname{Re}^{2}} }ρ[Vf∙∇]V=−∇p+μ∇2V+ρgβΔTρLU∞2​​[V∙∇]V=−ρLU∞2​​∇p+μL2U∞​​∇2V+ρgβΔT[ΔT]ref​ρ[V∙∇]V=−∇p+Re1​∇2V+Re2Gr​​
ρCp[V⃗∙∇]T=k∇2T+ΦρCpU∞ΔTL[V⃗∙∇]T=kΔTL2∇2T+U∞2L2ΦρCp[V⃗∙∇]T=kCpμ∙μρU∞L∇2T+[μρU∞L]U∞2CpΔT[Φμ]ρCp′[V⃗∙∇]T=1Pr⁡Re∇2T+EcRe⁡Φ\rho C_{p}[\vec{V} \bullet \nabla] T=k \nabla^{2} T+\Phi\\ \rho C_{p} \frac{U_{\infty} \Delta T}{L}[\vec{V} \bullet \nabla] T=k \frac{\Delta T}{L^{2}} \nabla^{2} T+\frac{U_{\infty}^{2}}{L^{2}} \Phi\\ \rho \mathcal{C}_{p}[\vec{V} \bullet \nabla] T=\frac{k}{C_{p} \mu} \bullet \frac{\mu}{\rho U_{\infty} L} \nabla^{2} T+\left[\frac{\mu}{\rho U_{\infty} L}\right] \frac{U_{\infty}^{2}}{C_{p} \Delta T}\left[\frac{\Phi}{\mu}\right]\\ \boxed{ \rho C_{p}^{\prime}[\vec{V} \bullet \nabla] T=\frac{1}{\operatorname{Pr} \mathrm{Re}} \nabla^{2} T+\frac{E c}{\operatorname{Re}} \Phi }ρCp​[V∙∇]T=k∇2T+ΦρCp​LU∞​ΔT​[V∙∇]T=kL2ΔT​∇2T+L2U∞2​​ΦρCp​[V∙∇]T=Cp​μk​∙ρU∞​Lμ​∇2T+[ρU∞​Lμ​]Cp​ΔTU∞2​​[μΦ​]ρCp′​[V∙∇]T=PrRe1​∇2T+ReEc​Φ​
[uux+vuy]=−pxρ+1Re⁡[uxx+uyy]+GrRe⁡21∙11+δ∙1δ=1+δ[11+1δ2]⏞simplify+δ2Gr\left[u u_{x}+v u_{y}\right]=-\frac{p_{x}}{\rho}+\frac{1}{\operatorname{Re}}\left[u_{x x}+u_{y y}\right]+\frac{G r}{\operatorname{Re}^{2}}\\ 1 \bullet \frac{1}{1}+\delta \bullet \frac{1}{\delta}=1+ \delta \overbrace{ \left[\frac{1}{1}+\frac{1}{\delta^{2}}\right] }^\text{simplify} +\delta^{2} G r[uux​+vuy​]=−ρpx​​+Re1​[uxx​+uyy​]+Re2Gr​1∙11​+δ∙δ1​=1+δ[11​+δ21​]​simplify​+δ2Gr
[uux+vuy]=ueuex+1Re⁡[uyy]+GrRe⁡2\boxed{\left[u u_{x}+v u_{y}\right]=u_{e} u_{e x}+\frac{1}{\operatorname{Re}}\left[u_{y y}\right]+\frac{G r}{\operatorname{Re}^{2}}}[uux​+vuy​]=ue​uex​+Re1​[uyy​]+Re2Gr​​
[uvx+vvy]=−pyρ+v[vxx+vyy]1∙δ1+δ∙δδ⏞LHS=1δ+δ2[δ1+δδ2]⏞⋘LHSpy≅0\left[u v_{x}+v v_{y}\right] = -\frac{p_{y}}{\rho}+v\left[v_{x x}+v_{y y}\right]\\ \overbrace{ 1 \bullet \frac{\delta}{1}+\delta \bullet \frac{\delta}{\delta} }^\text{LHS} = \frac{1}{\delta}+ \overbrace{ \delta^{2}\left[\frac{\delta}{1}+\frac{\delta}{\delta^{2}}\right] }^{\lll\text{LHS}} \\ \boxed{p_{y} \cong 0}[uvx​+vvy​]=−ρpy​​+v[vxx​+vyy​]1∙1δ​+δ∙δδ​​LHS​=δ1​+δ2[1δ​+δ2δ​]​⋘LHS​py​≅0​
ρCp[uTx+vTy]=Ec[upx+vpy]+1Pr⁡Re⁡[Txx+Tyy]+EcRe⁡2Φ1∙11+δ∙1δT=δ2Pr⁡[11+1δT2]⏞simplify\rho \mathcal{C}_{p}\left[u T_{x}+v T_{y}\right]=E c\left[u p_{x}+v p_{y}\right]+\frac{1}{\operatorname{Pr} \operatorname{Re}}\left[T_{x x}+T_{y y}\right]+\frac{E c}{\operatorname{Re}^{2}} \Phi\\ 1 \bullet \frac{1}{1}+\delta \bullet \frac{1}{\delta_{T}}=\frac{\delta^{2}}{\operatorname{Pr}} \overbrace{ \left[\frac{1}{1}+\frac{1}{\delta_{T}^{2}}\right] }^\text{simplify}ρCp​[uTx​+vTy​]=Ec[upx​+vpy​]+PrRe1​[Txx​+Tyy​]+Re2Ec​Φ1∙11​+δ∙δT​1​=Prδ2​[11​+δT2​1​]​simplify​
ρCp[uTx+vTy]=Ec[upx+vpy]+1Pr⁡Re⁡[Tyy]+EcRe⁡2[uy]2.\boxed{ \rho C_{p}\left[u T_{x}+v T_{y}\right]=E c\left[u p_{x}+v p_{y}\right]+\frac{1}{\operatorname{Pr} \operatorname{Re}}\left[T_{y y}\right]+\frac{E c}{\operatorname{Re}^{2}}\left[u_{y}\right]^{2}.}ρCp​[uTx​+vTy​]=Ec[upx​+vpy​]+PrRe1​[Tyy​]+Re2Ec​[uy​]2.​
ρCp[uTx+vTy]=Ec[upx+vpy]+1Pr⁡Re⁡[Tyy]+EcRe⁡2[uy]2\rho \mathcal{C}_{p}\left[u T_{x}+v T_{y}\right]=E c\left[u p_{x}+v p_{y}\right]+\frac{1}{\operatorname{Pr} \operatorname{Re}}\left[T_{y y}\right]+\frac{E c}{\operatorname{Re}^{2}}\left[u_{y}\right]^{2}ρCp​[uTx​+vTy​]=Ec[upx​+vpy​]+PrRe1​[Tyy​]+Re2Ec​[uy​]2
p∞ρ∞=RT∞=(cp−cv)T∞=γ−1γcpT∞\frac{p_{\infty}}{\rho_{\infty}}=R T_{\infty}=\left(c_{p}-c_{v}\right) T_{\infty}=\frac{\gamma-1}{\gamma} c_{p} T_{\infty}ρ∞​p∞​​=RT∞​=(cp​−cv​)T∞​=γγ−1​cp​T∞​
c∞=γp∞ρ∞=(γ−1)cpT∞c_{\infty}=\frac{\gamma p_{\infty}}{\rho_{\infty}}=(\gamma-1) c_{p} T_{\infty}c∞​=ρ∞​γp∞​​=(γ−1)cp​T∞​
Ec=U∞2cp(ΔTref)=U∞2cpT∞T∞(ΔTref)=(γ−1)U∞2c∞2T∞(ΔTrσf)=(γ−1)M2T∞(ΔTref)∼M2T∞(ΔTref)∼0.01T∞(ΔTref)\begin{align} E c &=\frac{U_{\infty}^{2}}{c_{p}\left(\Delta T_{r e f}\right)}=\frac{U_{\infty}^{2}}{c_{p} T_{\infty}} \frac{T_{\infty}}{\left(\Delta T_{r e f}\right)}\\ & =(\gamma-1) \frac{U_{\infty}^{2}}{c_{\infty}^{2}} \frac{T_{\infty}}{\left(\Delta T_{r \sigma f}\right)}=(\gamma-1) M^{2} \frac{T_{\infty}}{\left(\Delta T_{r e f}\right)}\\ &\sim M^{2} \frac{T_{\infty}}{\left(\Delta T_{r e f}\right)} \sim 0.01 \frac{T_{\infty}}{\left(\Delta T_{r e f}\right)} \end{align}Ec​=cp​(ΔTref​)U∞2​​=cp​T∞​U∞2​​(ΔTref​)T∞​​=(γ−1)c∞2​U∞2​​(ΔTrσf​)T∞​​=(γ−1)M2(ΔTref​)T∞​​∼M2(ΔTref​)T∞​​∼0.01(ΔTref​)T∞​​​​
1∙11+δ∙1δT=δ2Pr⁡[11+1δT2]→δT≅1Pr⁡1/2Re1/21 \bullet \frac{1}{1}+\delta \bullet \frac{1}{\delta_{T}}=\frac{\delta^{2}}{\operatorname{Pr}}\left[\frac{1}{1}+\frac{1}{\delta_{T}^{2}}\right]\rightarrow \delta_{T} \cong \boxed{\frac{1}{\operatorname{Pr}^{1 / 2} \mathrm{Re}^{1 / 2}}} \\1∙11​+δ∙δT​1​=Prδ2​[11​+δT2​1​]→δT​≅Pr1/2Re1/21​​