8.6.1 Boundary Conditions Categories
- Page ID
- 749
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)The governing equations that were developed earlier requires some boundary conditions and initial conditions. These conditions described physical situations that are believed or should exist or approximated. These conditions can be categorized by the velocity, pressure, or in more general terms as the shear stress conditions (mostly at the interface). For this discussion, the shear tensor will be separated into two categories, pressure (at the interface direction) and shear stress (perpendicular to the area). A common velocity condition is that the liquid has the same value as the solid interface velocity. In the literature, this condition is referred as the "no slip'' condition. The solid surface is rough thus the liquid participles (or molecules) are slowed to be at the solid surface velocity. This boundary condition was experimentally observed under many conditions yet it is not universal true. The slip condition (as oppose to "no slip'' condition) exist in situations where the scale is very small and the velocity is relatively very small. The slip condition is dealing with a difference in the velocity between the solid (or other material) and the fluid media. The difference between the small scale and the large scale is that the slip can be neglected in the large scale while the slip cannot be neglected in the small scale. In another view, the difference in the velocities vanishes as the scale increases.
Fig. 8.13 Dimensional free surface describing \(\widehat{\pmb{n}}\) and \(\widehat{\pmb{t}}\).
Another condition which affects whether the slip condition exist is how rapidly of the velocity change. The slip condition cannot be ignored in some regions, when the flow is with a strong velocity fluctuations. Mathematically the "no slip'' condition is written as
\[ \label{dif:eq:noSlipCondV} \widehat{\mathbf{t}} \cdot \left( \pmb{U}_{fluid} - \pmb{U}_{boundary} \right) = 0 \]
where \(\widehat{\mathbf{n}}\) is referred to the area direction (perpendicular to the area see Figure ??. While this condition (61) is given in a vector form, it is more common to write this condition as a given velocity at a certain point such as
\[ \label{dif:eq:noSlipCond} U(\ell) = U_{\ell} \]
Note, the "no slip'' condition is applicable to the ideal fluid ("inviscid flows'') because this kind of flow normally deals with large scales. The 'slip' condition is written in similar fashion to equation (61) as
\[
\label{dif:eq:slipCondV}
\widehat{\mathbf{t}} \cdot \left( \pmb{U}_{fluid} - \pmb{U}_{boundary} \right) = f(Q,scale, etc)
\]
The condition (61) can be mathematically represented in another way for free surface
conditions. To make sure that all the material is accounted for in the control volume (does not cross the free surface), the relative perpendicular velocity at the interface must be zero. The location of the (free) moving boundary can be given as \(f(\widehat{\pmb{r}},t) =0\) as the equation which describes the bounding surface. The perpendicular relative velocity at the surface must be zero and therefore
\[ \label{dif:eq:perpendicularUSurface} \dfrac{Df}{Dt} = 0 \quad \mbox{ on the surface } f (\widehat{\pmb{r}},t) = 0 \]
This condition is called the kinematic boundary condition. For example, the free surface in the two dimensional case is represented as \(f(t,x,y)\). The condition becomes as
\[
\label{dif:eq:2perpendicularUSurface}
0 = \dfrac{\partial f}{\partial t} +
U_x\, \dfrac{\partial f}{\partial x} +
U_y\, \dfrac{\partial f}{\partial y}
\]
\[
\label{dif:eq:JumpStress}
\Delta \boldsymbol{\tau}^{(n)} = 0 =
\Delta {\boldsymbol{\tau}^{(n)}}_{\text{ upper surface}} -
\Delta {\boldsymbol{\tau}^{(n)}}_{\text{ lower surface}}
\]
where the index \((n)\) indicate that shear stress are normal (in the surface area).
If the surface is straight there is no jump in the shear stress. The condition with curved surface are out the scope of this book yet mathematically the condition is given as without explanation as
\[
\label{dif:eq:bcCurvedSurface}
\widehat{\pmb{n}}\cdot {\boldsymbol{\tau}^{(n)}} = \sigma \left( \dfrac{1\dfrac{}{}}{R_1} + \dfrac{1}{R_2}\right) \
\widehat{\pmb{t}}\cdot {\boldsymbol{\tau}^{(t)}} = - \widehat{\pmb{t}}\cdot \nabla\sigma
\]
where \(\widehat{\pmb{n}}\) is the unit normal and \(\widehat{\pmb{t}}\) is a unit tangent to the surface (notice that direction pointed out of the "center'' see Figure ??) and \(R_1\) and \(R_2\) are principal radii. One of results of the free surface condition (or in general, the moving surface condition) is that integration constant is unknown). In same instances, this constant is determined from the volume conservation. In index notation equation (67) is written
\[
\label{dif:eq:bcCurvedSurfaceIndex}
\tau_{ij}^{(1)}\,n_j + \sigma \, n_i \left( \dfrac{1}{R_1} + \dfrac{1}{R_2}\right) =
\tau_{ij}^{(2)}\,n_j
\]
where \(1\) is the upper surface and \(2\) is the lower surface. For example in one dimensional
\[
\label{dif:eq:nAndT}
\begin{array}{rl}
\widehat{\pmb{n}} = & \dfrac{\left(-f^\prime{}(x), 1 \right) }{\sqrt{1+ \left(f^\prime{}(x)\right)^2 }} \
\widehat{\pmb{t}} = & \dfrac{\left(1,f^\prime{}(x) \right) }{\sqrt{1+ \left(f^\prime{}(x)\right)^2 }}
\end{array}
\]
the unit vector is given as two vectors in \(x\) and \(y\) and the radius is given by equation (??). The equation is given by
\[
\label{dif:eq:1DfreeSurface}
\dfrac{\partial f}{\partial t} + U_x \dfrac{\partial f}{\partial x} = U_y
\]
The Pressure Condition
The second condition that commonality prescribed at the interface is the static pressure at a specific location. The static pressure is measured perpendicular to the flow direction. The last condition is similar to the pressure condition of prescribed shear stress or a relationship to it. In this category include the boundary conditions with issues of surface tension which were discussed earlier. It can be noticed that the boundary conditions that involve the surface tension are of the kind where the condition is given on boundary but no at a specific location.
Gravity as Driving Force
The body forces, in general and gravity in a particular, are the condition that given on the flow beside the velocity, shear stress (including the surface tension) and the pressure. The gravity is a common body force which is considered in many fluid mechanics problems. The gravity can be considered as a constant force in most cases (see for dimensional analysis for the reasons).
Shear Stress and Surface Tension as Driving Force
Fig. 8.14 Kerosene lamp.
If the fluid was solid material, pulling the side will pull all the material. In fluid (mostly liquid) shear stress pulling side (surface) will have limited effect and yet sometime is significant and more rarely dominate. Consider, for example, the case shown in Figure 8.14. The shear stress carry the material as if part of the material was a solid material. For example, in the kerosene lamp the burning occurs at the surface of the lamp top and the liquid is at the bottom. The liquid does not move up due the gravity (actually it is against the gravity) but because the surface tension.
Fig. 8.15 Flow in a kendle with a surfece tension gradient.
The physical conditions in Figure 8.14 are used to idealize the flow around an inner rode to understand how to apply the surface tension to the boundary conditions. The fluid surrounds the rode and flows upwards. In that case, the velocity at the surface of the inner rode is zero. The velocity at the outer surface is unknown. The boundary condition at outer surface given by a jump of the shear stress. The outer diameter is depends on the surface tension (the larger surface tension the smaller the liquid diameter). The surface tension is a function of the temperature therefore the gradient in surface tension is result of temperature gradient. In this book, this effect is not discussed. However, somewhere downstream the temperature gradient is insignificant. Even in that case, the surface tension gradient remains. It can be noticed that, under the assumption presented here, there are two principal radii of the flow. One radius toward the center of the rode while the other radius is infinite (approximatly). In that case, the contribution due to the curvature is zero in the direction of the flow (see Figure 8.15). The only (almost) propelling source of the flow is the surface gradient (\(\dfrac{\partial \sigma}{\partial n}\)).
Contributors and Attributions
Dr. Genick Bar-Meir. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or later or Potto license.