[Click here for a PDF of this post with nicer formatting and figures if the post had any (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]

Disclaimer.

This problem set is as yet ungraded.

Problem Q1.

Background.

In fluid convection problems one can make several choices for characteristic velocity scales. Some choices are given below for example:

$U_1 = g \alpha d^2 \nabla T/\nu$
$U_2 = \nu/d$
$U_3 = \sqrt{ g \alpha d \nabla T }$
$U_4 = \kappa/ d$

where is the acceleration due to gravity, $\alpha = ({\partial {V}}/{\partial {T}})/V$ is the coefficient of volume expansion, length scale associated with the problem, $\nabla T$ is the applied temperature difference, $\nu$ is the kinematic viscosity and $\kappa$ is the thermal diffusivity.

For :
Observing that

$\begin{aligned}\left[ \frac{\partial {\mathbf{u}}}{\partial {t}} \right] = [ \nu \boldsymbol{\nabla}^2 \mathbf{u} ]\end{aligned} \hspace{\stretch{1}}(2.1)$

we must have

$\begin{aligned}[\nu] = \frac{1}{{[t] [\boldsymbol{\nabla}^2]}} = \frac{1}{{T}} L^2\end{aligned} \hspace{\stretch{1}}(2.2)$

We also find

$\begin{aligned}[\alpha] = \frac{1}{{[V]}} \left[ \frac{\partial {V}}{\partial {T}} \right] = \frac{1}{{[K]}},\end{aligned} \hspace{\stretch{1}}(2.3)$

so that

$\begin{aligned}[U_1] = \frac{L}{T \not{{T}}} \not{{\frac{1}{{K}}}} L^2 \not{{K}} \frac{\not{{T}}}{L^2} = \frac{L}{T}\end{aligned} \hspace{\stretch{1}}(2.4)$
For :
$\begin{aligned}[U_2] = \frac{L^2}{T} \frac{1}{{L}} = \frac{L}{T}.\end{aligned} \hspace{\stretch{1}}(2.5)$
For
$\begin{aligned}[U_3] = \sqrt{ \frac{L}{T^2} \frac{1}{{K}} L K } = \frac{L}{T}\end{aligned} \hspace{\stretch{1}}(2.6)$
For
According to http://scienceworld.wolfram.com/physics/ThermalDiffusivity.html, the thermal diffusivity is defined by

$\begin{aligned}\frac{\partial {T}}{\partial {t}} = \kappa \boldsymbol{\nabla}^2 T\end{aligned} \hspace{\stretch{1}}(2.7)$

so that

$\begin{aligned}[\kappa] = \frac{1}{{[t][\boldsymbol{\nabla}^2]}} = \frac{L^2}{T}\end{aligned} \hspace{\stretch{1}}(2.8)$

That gives us

$\begin{aligned}[U_4] = \frac{L^2}{T} \frac{1}{{L}} = \frac{L}{T}.\end{aligned} \hspace{\stretch{1}}(2.9)$

We’ve verified that all of these have dimensions of velocity.

Part 1. Statement. Check whether the dimensions match in each case above.

Solution.

Part 2. Statement. Pure liquid.

For pure liquid, say pure water at room temperature, one has the following estimates in cgs units:

$\begin{aligned}\alpha &\sim 10^{-4} \\ \kappa &\sim 10^{-3} \\ \nu &\sim 10^{-2}\end{aligned}$

For a $d \sim 1 \text{cm}$ layer depth and a ten degree temperature drop convective velocities have been experimentally measured of about $10^{-2}$ .

With $g \sim 10^{-3}$ , calculate the values of U_1 , U_2 , U_3 , and U_4 . Which ones of the characteristic velocities (U_1 , U_2 , U_3, U_4) do you think are suitable for nondimensionalising the velocity in Navier-Stokes/Energy equation describing the water convection problem?

We have

$\begin{aligned}U_1 &\sim 10^{-3} 10^{-4} (1)^2 10^1 \frac{1}{{10^{-2}}} = 10^{-4} \\ U_2 &\sim 10^{-2}/1 = 10^{-2} \\ U_3 &\sim \sqrt{ 10^{-3} 10^{-4} (1) 10^1 } = 10^{-3} \\ U_4 &\sim 10^{-3}/1 = 10^{-3}\end{aligned} \hspace{\stretch{1}}(2.10)$

Use of $U_2 = \nu/d$ gives the closest match to the measured characteristic velocity of $10^{-2}$ .

Solution.

Part 3. Statement. Mantle convection.

For mantle convection, we have

$\begin{aligned}\alpha &\sim 10^{-5} \\ \nu &\sim 10^{21} \\ \kappa &\sim 10^{-2} \\ d &\sim 10^8 \\ \nabla T &\sim 10^3,\end{aligned}$

and the actual convective mantle velocity is $10^{-8}$ . Which of the characteristic velocities should we use to nondimensionalise Navier-Stokes/Energy equations describing mantle convection?

Solution

Let’s compute the characteristic velocities again with the mantle numbers

$\begin{aligned}U_1 &\sim \frac{10^{-3} 10^{-5} 10^{16} 10^3 }{10^{21}} = 10^{-10} \\ U_2 &\sim \frac{10^{21}}{10^8} = 10^{13} \\ U_3 &\sim \sqrt{ 10^{-3} 10^{-5} 10^8 10^3 } \sim 10^1 \\ U_4 &\sim \frac{10^{-2}}{10^8} = 10^{-10}\end{aligned} \hspace{\stretch{1}}(2.14)$

Both U_1 and U_4 come close to the actual convective mantle velocity of $10^{-8}$ . Use of U_1 to nondimensionalise is probably best, since it has more degrees of freedom, and includes the gravity term that is probably important for such large masses.

Problem Q2.

Statement

Nondimensionalise N-S equation

$\begin{aligned}\rho \frac{\partial {\mathbf{u}}}{\partial {t}} + \rho (\mathbf{u} \cdot \boldsymbol{\nabla}) \mathbf{u} = - \boldsymbol{\nabla} p + \mu \boldsymbol{\nabla}^2 \mathbf{u} + \rho g \hat{\mathbf{z}}\end{aligned} \hspace{\stretch{1}}(3.18)$

where $\hat{\mathbf{z}}$ is the unit vector in the direction. You may scale:

velocity with the characteristic velocity ,
time with , where is the characteristic length scale,
pressure with $\rho U^2$ ,

Reynolds number $\text{Re} = R U \rho/ \mu$ and Froude number $\text{Fr} = g R/U$ .

Solution

Let’s start by dividing by $g \rho$ , to make all terms (most obviously the $\hat{\mathbf{z}}$ term) dimensionless.

$\begin{aligned}\frac{1}{{g}} \frac{\partial {\mathbf{u}}}{\partial {t}} + \frac{1}{{g}} (\mathbf{u} \cdot \boldsymbol{\nabla}) \mathbf{u} = - \frac{1}{{g \rho}} \boldsymbol{\nabla} p + \frac{\mu}{g \rho} \boldsymbol{\nabla}^2 \mathbf{u} + \hat{\mathbf{z}}.\end{aligned} \hspace{\stretch{1}}(3.19)$

Our suggested replacements are

$\begin{aligned}\mathbf{u} &= U \mathbf{u}' \\ \frac{\partial {}}{\partial {t}} &= \frac{U}{R} \frac{\partial {}}{\partial {t'}} \\ p &= \rho U^2 p' \\ \boldsymbol{\nabla} &= \frac{1}{{R}} \boldsymbol{\nabla}'.\end{aligned} \hspace{\stretch{1}}(3.20)$

Plugging these in we have

$\begin{aligned}\frac{U^2}{g R} \frac{\partial {\mathbf{u}'}}{\partial {t'}} + \frac{U^2}{g R} (\mathbf{u}' \cdot \boldsymbol{\nabla}') \mathbf{u}' = - \frac{\not{{\rho}} U^2}{g \not{{\rho}} R} \boldsymbol{\nabla}' p' + \frac{\mu U}{g \rho R^2} {\boldsymbol{\nabla}'}^2 \mathbf{u}' + \hat{\mathbf{z}}.\end{aligned} \hspace{\stretch{1}}(3.24)$

Making a $\text{Fr} = gR/U$ replacement, using the Froude number, we have

$\begin{aligned}\frac{U}{\text{Fr}} \frac{\partial {\mathbf{u}'}}{\partial {t'}} + \frac{U}{\text{Fr}} (\mathbf{u}' \cdot \boldsymbol{\nabla}') \mathbf{u}' = - \frac{U}{\text{Fr}} \boldsymbol{\nabla}' p' + \frac{\mu }{\text{Fr} \rho R} {\boldsymbol{\nabla}'}^2 \mathbf{u}' + \hat{\mathbf{z}}.\end{aligned} \hspace{\stretch{1}}(3.25)$

Scaling by $\text{Fr}/U$ we tidy things up a bit, and also allow for insertion of the Reynold’s number

$\begin{aligned}\frac{\partial {\mathbf{u}'}}{\partial {t'}} + (\mathbf{u}' \cdot \boldsymbol{\nabla}') \mathbf{u}' = - \boldsymbol{\nabla}' p' + \frac{1}{\text{Re}} {\boldsymbol{\nabla}'}^2 \mathbf{u}' + \frac{\text{Fr}}{U}\hat{\mathbf{z}}.\end{aligned} \hspace{\stretch{1}}(3.26)$

Observe that the dimensions of Froude’s number is that of velocity

$\begin{aligned}[\text{Fr}] = [g] T = \frac{L}{T},\end{aligned} \hspace{\stretch{1}}(3.27)$

so that the end result is dimensionless as desired. We also see that Froude’s number, characterizes the significance of the body force for fluid flow at the characteristic velocity. This is consistent with [1] where it was stated that the Froude number is used to determine the resistance of a partially submerged object moving through water, and permits the comparison of objects of different sizes (complete with pictures of canoes of various sizes that Froude built for such study).

Problem Q3.

Statement

In case of Stokes’ boundary layer problem (see class note) calculate shear stress on the plate y = 0 . What is the phase difference between the velocity of the plate $U(t) = U_0 \cos\omega t$ and the shear stress on the plate?

Solution

We found in class that the velocity of the fluid was given by

$\begin{aligned}u(y, t) = U_0 e^{-\lambda y} \cos(\lambda y - \omega t)\end{aligned} \hspace{\stretch{1}}(4.28)$

where

$\begin{aligned}\lambda = \sqrt{\frac{\omega}{2 \nu}}\end{aligned} \hspace{\stretch{1}}(4.29)$

Calculating our shear stress we find

$\begin{aligned}\mu \frac{\partial {u}}{\partial {y}} &= U_0 \lambda \mu e^{-\lambda y}\left(-1- \sin(\lambda y - \omega t)\right)\end{aligned}$

and on the plate ( y = 0 ) this is just

$\begin{aligned}{\left.{{\mu \frac{\partial {u}}{\partial {y}}}}\right\vert}_{{y = 0}} = U_0 \lambda \mu (-1 + \sin(\omega t)).\end{aligned} \hspace{\stretch{1}}(4.30)$

We’ve got a constant term, plus one that is sinusoidal. Observing that

$\begin{aligned}\cos x &= \text{Real}( e^{ix} ) \\ \sin x &= \text{Real}( -i e^{ix} ) = \text{Real}( e^{i (x - \pi/2)} ),\end{aligned} \hspace{\stretch{1}}(4.31)$

The phase difference between the non-constant portion of the shear stress at the plate, and the plate velocity $U(t) = U_0 \cos\omega t$ is just $-\pi/2$ . The shear stress at the plate lags the driving velocity by 90 degrees.

References

[1] Wikipedia. Froude number — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 27-March-2012]. http://en.wikipedia.org/w/index.php?title=Froude_number&oldid=479498080.

PHY454H1S Continuum mechanics. Problem Set 3. Velocity scaling, non-dimensionalisation, boundary layers.

Disclaimer.

Problem Q1.

Background.

Part 1. Statement. Check whether the dimensions match in each case above.

Solution.

Part 2. Statement. Pure liquid.

Solution.

Part 3. Statement. Mantle convection.

Solution

Problem Q2.

Statement

Solution

Problem Q3.

Statement

Solution

References

Trending Articles

Practice Sheet of Right form of verbs for HSC Students

Download: FK ft Shenky – Nakuyewa ”Prod by: Shenky”

How to win at Markstrat (Markstrat Tips and Tricks) – Vodites

Ominde Commission Report and Recommendations – Ominde Report of 1964

Bureau of Internal Revenue: Regional Offices (Directory)

GO 53 on Enhancement of Ex-gratia upto 5 Lakhs Toddy Tappers in Telangana

Cakewalk CA-2A Leveling Amplifier v2.0.1.97 WiN, v2.0.1.96 OSX Incl Keygen

Mp3 Download: Mdu - Kunjenjenjena

How the kill the job , when DTP request running for long hours.

Microsoft Intune から展開しているアプリのアップデートについて

18-year-old girl was beaten for half an hour by two Northampton men in 'an...

Car crash in Dunton Bassett leaves driver in critical condition

Macky 2, Two Others In Road Accident

Application log 00000000000000089514: Could not convert queue DLVST90CLNT

Detroit mafia: D’Anna Brothers agree to plea deal

Delivery block field greyed out using VA02

Muloraki Au

【個人撮影】スマホのプライベート映像♪「中に出さないで///」カラオケ屋での生ハメ撮りが流出ｗ【リベンジポルノ】＠PornHub

BREAKING NEWS: Diamond Platnumz Is Reported Dead After Ghastly Car Accident

FIAT 500 B0111 B0112