Fluid Mechanics for MAP/Dimensional Analysis

>back to Chapters of Fluid Mechanics for MAP

>back to Chapters of Microfluid Mechanics

Motivation

File:Motivation1.png

Consider the steady flow around a submerged sphere. Say we are interested in the drag force on the sphere.

Template:Center top$F_D=f(U_{\mathrm{\infty }},D,\rho ,\mu )$Template:Center bottom

Let's say 10Bold text measurement points are enough for one curve, which shows the effect of one parameter while keeping the others constant. In order to monitor the effect of all variables, $10^4$ separate measurements are needed. When each measurement takes 1/2 hr, all test series need 2.5 years.

By using dimensional analysis, we can reduce $F_D$ to an equivalent form:

Template:Center top$\frac{F_D}{\rho U_{\mathrm{\infty }}^2{D}^2}=g\left(\frac{\rho U_{\mathrm{\infty }}D}{\mu }\right)\text{or commonly written as}\underset{Dragcoefficient}{\underbrace{C_D}}=\underset{Reynoldsnumber}{\underbrace{g\left(\frac{\rho U_{\mathrm{\infty }}D}{\mu }\right)}}$Template:Center bottom

Both sides are dimensionless. Function $g$ is different than $f$ but it contains all necessary information.

Full scale models might be impossible or prohibitively expensive. There might be numerous number of parameters to be tested.

Figure:Commercial Aircraft

Figure:Model flight in windtunnel

Under which conditions these tests are realistic? Are the both flows dynamically similar?

“Dimensional analysis is a method which reduces number and complexity of experimental variables which affect a given physical phenomenon, using a sort of compacting technique.” (White 1979)

Dimensional analysis (DA) helps us to formulate dimensionless forms of governing equations and simplify them by the determination insignificant terms:

Template:Center top$\rho [\frac{\partial U_j}{\partial t}+U_i\frac{\partial U_j}{\partial x_i}]=-\frac{\partial P}{\partial x_j}+\frac{\partial }{\partial x_i}[\mu (\frac{\partial U_j}{\partial x_i}+\frac{\partial U_i}{\partial x_j})-\frac{2}{3}{\delta }_{ij}\mu \frac{\partial U_k}{\partial x_k}]\rho g_j$Template:Center bottom

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-Eu\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-\frac{\partial P^{*}}{\partial x_j^{*}}$Template:Center bottom

This is the non-dimensional form of Euler equation.

Benefits of dimensional analysis (DA)

One obvious goal of our efforts should be to obtain the most information from the fewest experiments. Dimensional analysis (DA) is an important tool that often helps us to achieve this goal.

The benefits of this dimensionless representation are:

1. We do not have to vary each variable: It is enough that we vary only Re and measure CD. As a result, we have enormous time and money saving.
2. Dimensional analysis (DA) provides “scaling laws” which can convert data from a cheap small model to provide information about an expensive large prototype.
3. DA helps our thinking and planning of an experiment or a theory.
   • It suggests dimensionless ways of writing equations.
   • One can find variables which can be discarded.

As a result, we can avoid to waste money and time while we use experimental and computational resources.

Derivation of dimensionless numbers

File:Motivation1.png

Buckingham Π theorem provides a method for computing sets of dimensionless numbers from the given variables, even if the form of the equation is still unknown.

Consider the steady flow around a sphere, which shows no separation. We will find out the dimensionless numbers governing this flow:

Step 1: List all the parameters involved. Let $n$ be the number of parameters.

Template:Center top$F_D,U_{\mathrm{\infty }},D,\rho ,\mu \to n=5$Template:Center bottom

Step 2: Select a set of primary dimensions, e.g. $MLtTq$ (mass, length, time, temperature, elec. charge or $FLtTq$ (Force, length, ...).

Template:Center top$M,L\text{and}t\text{are selected to be the primary dimensions.}$Template:Center bottom

Step 3: List the dimensions of all parameters. Let $r$ be the number of primary dimensions.

Template:Center top$F,U_{\mathrm{\infty }},D,\rho ,\mu \to \frac{ML}{t^2},\frac{L}{t},L,\frac{M}{L^3},\frac{M}{Lt}\to r=3\text{primary dimensions}$Template:Center bottom

Step 4: Select $r$ number of parameters, including all the primary dimensions. No two parameters should have the same net dimension differing by only a single exponent: for example do not include both $L$ and $L^4$.

Template:Center top$U_{\mathrm{\infty }},D,\rho \to m=r=3\text{repeating parameters}$Template:Center bottom

Step 5: Set up dimensional equations combining the parameters selected in Step 4 with each of the other parameters in turn to form $n-m$ dimensionless groups.

$n-m=2$ dimensionless groups will result. The dimensional equations are:

Template:Center top${\mathrm{\Pi }}_1=U_{\mathrm{\infty }}^a{D}^b{\rho }^{c}F={\left(\frac{L}{t}\right)}^a{\left(L\right)}^b{\left(\frac{M}{L^3}\right)}^c\left(\frac{ML}{t^2}\right)=M^0{L}^0{t}^0$Template:Center bottom

Template:Center top${\mathrm{\Pi }}_2=U_{\mathrm{\infty }}^d{D}^e{\rho }^f\mu ={\left(\frac{L}{t}\right)}^d{\left(L\right)}^e{\left(\frac{M}{L^3}\right)}^f\left(\frac{M}{Lt}\right)=M^0{L}^0{t}^0$Template:Center bottom

solving above equations delivers $a=-2,b=-2,c=-1,d=-1,e=-1\text{and}f=-1$. Hence

Template:Center top${\mathrm{\Pi }}_1=\frac{F}{\rho U_{\mathrm{\infty }}^2{D}^2}\text{and}{\mathrm{\Pi }}_2=\frac{\mu }{U_{\mathrm{\infty }}D\rho }=\frac{1}{Re}$Template:Center bottom

The functional relationship ${\mathrm{\Pi }}_1=f({\mathrm{\Pi }}_2)$ should be determined experimentally and/or numerically.

Application examples of non-dimensional numbers

Flow around a cylinder and a sphere

File:Flow around.png

In a very large Re range the wake is unsteady, quasi time-periodic.

Important dimensionless numbers are:

$Re=\frac{\rho Ud}{\mu }=\frac{Ud}{\nu }$ where $\nu =\frac{\mu }{\rho }$ is the kinematik viscosity.

$St=\frac{fd}{U}=\frac{wd}{2\pi U}$ where f is the dominant frequency of the wake structures.

$C_D=\frac{F_D}{0.5\rho U^{2}A}$ where A is the projected area to the flow.

$\frac{L}{d}$ is also important for flow around a cylinder.

File:Flow around 2a.png

File:Flow around 2b.png

Losses in pipe flows

The pressure drop along the pipe can be written as: $\mathrm{\Delta }P=g(D,L,e,V,\rho ,\mu )$ where $e$ is the average height of surface roughness and V is the average velocity. Template:Center top$\zeta =\frac{\mathrm{\Delta }P}{0.5\rho V^2}=\mathrm{\Phi }(Re,\frac{L}{D},\frac{e}{D})\to \zeta =\frac{L}{D}f(Re,\frac{e}{D})$Template:Center bottom where $f$ is the friction factor. Moody diagram

Velocity profiles for laminar (upper) and turbulent (lower) states at the same mass flow rate

Lift and drag forces on airfoils

File:Airfoil c.svg

File:Stallpoint.svg

File:Airfoil a.png

File:Airfoil b.png

Template:Center top$C_L=\frac{F_L}{0.5\rho U^{2}A}\to C_L=f(\alpha ,Re)$Template:Center bottom

Template:Center top$C_D=\frac{F_D}{0.5\rho U^{2}A}\to C_D=f(\alpha ,Re)$Template:Center bottom

Non-dimensionalization of Governing Equations

An isothermal and newtonian fluid flow is governed by the following equations

Conservation of Mass (continuity equation):

Template:Center top$\frac{\partial \rho }{\partial t}+\frac{\partial \rho U_i}{\partial x_i}=0$Template:Center bottom

Conservation of Momentum

Template:Center top$\rho [\frac{\partial U_j}{\partial t}+U_i\frac{\partial U_j}{\partial x_i}]=-\frac{\partial P}{\partial x_j}+\frac{\partial }{\partial x_i}[\mu (\frac{\partial U_j}{\partial x_i}+\frac{\partial U_i}{\partial x_j})-\frac{2}{3}{\delta }_{ij}\mu \frac{\partial U_k}{\partial x_k}]\rho g_j$Template:Center bottom

One can normalize each variable with a characteristic parameter, which has the same dimension

Template:Center top$x_i=L_c{x}_i^{*},t=t_c{t}^{*},U_j=U_c{U}_j^{*},\rho ={\rho }_c{\rho }^{*},g_j=g_c{g}_j^{*},\mu ={\mu }_c{\mu }^{*},P=\mathrm{\Delta }{P}_c{P}^{*}$Template:Center bottom

Inserting these into the continuity equation

Template:Center top$\frac{\partial {\rho }_c{\rho }^{*}}{\partial t_c{t}^{*}}+\frac{\partial {\rho }_c{\rho }^{*}{U}_c{U}_i^{*}}{\partial L_c{x}_i^{*}}=0\to \frac{{\rho }_c}{t_c}\frac{\partial {\rho }^{*}}{\partial t^{*}}+\frac{{\rho }_c{U}_c}{L_c}\frac{\partial {\rho }^{*}{U}_i^{*}}{\partial x_i^{*}}=0$Template:Center bottom

Dividing by $\frac{{\rho }_c{U}_c}{L_c}$

Template:Center top$\underset{St}{\underbrace{\frac{L_c}{t_c{U}_c}}}\frac{\partial {\rho }^{*}}{\partial t^{*}}+\frac{\partial {\rho }^{*}{U}_i^{*}}{\partial x_i^{*}}=0$Template:Center bottom

where St: Strouhal Number

Template:Center top$St=\frac{L_c}{t_c{U}_c}=\frac{\left(\frac{U_c}{t_c}\right)}{\frac{U_c^2}{L_c}}=\frac{localacceleration}{spatialacceleration}$Template:Center bottom

Similarly, we can apply this technique to make the momentum equation non-dimensional:

Template:Center top${\rho }_c{\rho }^{*}[\frac{U_c}{t_c}\frac{\partial U_j^{*}}{\partial t^{*}}+\frac{U_c^2}{L_c}{U}_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-\frac{\mathrm{\Delta }{P}_c}{L_c}\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{{\mu }_c{U}_c}{L_c^2}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+{\rho }_c{g}_c{\rho }^{*}{g}_j^{*}$Template:Center bottom

Dividing by $\frac{{\rho }_c{U}_c^2}{L_c}$

Template:Center top${\rho }^{*}[\underset{St}{\underbrace{\frac{L_c}{t_c{U}_c}}}\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=\underset{Eu}{\underbrace{-\frac{\mathrm{\Delta }{P}_c}{{\rho }_c{U}_c^2}}}\frac{\partial P^{*}}{\partial x_j^{*}}+\underset{\frac{1}{Re}}{\underbrace{\frac{{\mu }_c}{{\rho }_c{U}_c{L}_c}}}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\underset{\frac{1}{Fr}}{\underbrace{\frac{g_c{L}_c}{U_c^2}}}{\rho }^{*}{g}_j^{*}$Template:Center bottom

Eu: Euler Number $=\frac{\mathrm{\Delta }{P}_c}{{\rho }_c{U}_c^2}=\frac{staticpressure}{dynamicpressure}$

Re: Reynolds Number $=\frac{{\rho }_c{U}_c{L}_c}{{\mu }_c}=\frac{\frac{{\rho }_c{U}_c^2}{L_c}}{\frac{{\mu }_c{U}_c}{L_c^2}}=\frac{inertialforces}{viscousforces}$

Fr: Froude Number $=\frac{U_c^2}{g_c{L}_c}=\frac{\frac{{\rho }_c{U}_c^2}{L_c}}{{\rho }_c{g}_c}=\frac{inertialforces}{gravitationalforces}$

The resulting dimensionless conservation of mass and of momentum equations are:

Template:Center top$St\frac{\partial {\rho }^{*}}{\partial t^{*}}+\frac{\partial {\rho }^{*}{U}_i^{*}}{\partial x_i^{*}}=0$Template:Center bottom

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-Eu\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

These equation indicates that solution of similar problems for two problems, say A and B, can only be same when the dimensionless numbers are equal:

Template:Center top$S{t}_A=S{t}_B,E{u}_A=E{u}_B,R{e}_A=R{e}_B,F{r}_A=F{r}_B$Template:Center bottom

Furthermore, we should warrant the similarity of the boundary conditions and this requires geometric similarity.

Selection of characteristic variables and Simplifications to governing equations

Consider the momentum transport equation:

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-Eu\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

For very low St, i.e. when tc < Uc/Lc, in other words when local acceleration much smaller than spatial acceleration, the flow can be accepted to be steady. Thus the first term at LHS can be neglected, i.e. the momentum equation reads

Template:Center top${\rho }^{*}{U}_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}=-Eu\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

When we select: $\mathrm{\Delta }{P}_c={\rho }_c{U}_c^2\to Eu=1$

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

In one phase flows, this kind of selection of the characteristic pressure is common. Thus, one can see that the pressure gradient can not be neglected.

For flows in which cavitation occurs, one can use the vapor pressure of the fluid in the characteristic pressure term:

$\mathrm{\Delta }{P}_c=P-P_v\to Eu=Ca=CavitationNumber=\frac{P-P_v}{{\rho }_c{U}_c^2}$

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-Ca\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]+\frac{1}{Fr}{\rho }^{*}{g}_j^{*}$Template:Center bottom

The Froude number is important when gravity starts to play an important role, i.e. when buoyancy effects occurs due to density differences and, especially for multi-phase flows. For example, we deal with a two-phase flow, for which two fluids can be accepted to be immiscible, like water and air in the spill way of a dam. In order to introduce the effect of density into the Fr number, the characteristic gravity in the Fr number can be selected to be:

Template:Center top$g_c=g\frac{{\rho }_2-{\rho }_1}{{\rho }_2}$ where ${\rho }_2>{\rho }_1\to Fr=\left(\frac{1}{g({\rho }_2-{\rho }_1)}\right)\left(\frac{{\rho }_2{U}_c^2}{L_c}\right)=\frac{inertialforces}{bouyancyforces}$Template:Center bottom

Thus, when the density difference becomes negligible, Fr becomes very large and the momentum equation reduces to:

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-Eu\frac{\partial P^{*}}{\partial x_j^{*}}+\frac{1}{Re}\frac{\partial }{\partial x_i^{*}}[{\mu }^{*}(\frac{\partial U_j^{*}}{\partial x_i^{*}}+\frac{\partial U_i^{*}}{\partial x_j^{*}})-\frac{2}{3}{\delta }_{ij}{\mu }^{*}\frac{\partial U_k^{*}}{\partial x_k^{*}}]$Template:Center bottom

The Froude number can be important also for one phase flows with large density gradients. In order to introduce this effect, the characteristic gravity in the Fr number can be written by using the density gradient:

Template:Center top$g_c=g{\left(\frac{dp}{dx}\right)}_c\frac{L_c}{{\rho }_c}\to Fr=\left(\frac{1}{g{\left(\frac{dp}{dx}\right)}_c{L}_c}\right)\left(\frac{{\rho }_c{U}_c^2}{L_c}\right)=\frac{initialforces}{bouyancyforces}$Template:Center bottom

This form of the Fr number can be used in the field of convective heat transfer at which the characteristic velocities are very low such that buoyancy dominates the flow motion. In the mixing of fluids having different densities, this form can also be one of the important non-dimensional numbers.

Another simplification can be made: when the Re number of the flow is very high, the viscous stresses can be eliminated in regions where the velocity gradients are low. For example, for the cruising aircrafts the Re number is in the order of millions. Therefore, it is common approach, to use the inviscid form of Navier Stokes equations for flow simulations around cruising aircrafts, except in the vicinity of the aircraft body. For such a flow, we can neglect the gravitational effects, since Fr becomes very large. Moreover, Eu can be accepted to be 1 as shown before. Thus, the governing momentum equation reduces to the non-dimensional form of Euler equations:

Template:Center top${\rho }^{*}[St\frac{\partial U_j^{*}}{\partial t^{*}}+U_i^{*}\frac{\partial U_j^{*}}{\partial x_i^{*}}]=-\frac{\partial P^{*}}{\partial x_j^{*}}$Template:Center bottom

Similarity considerations

Once the variables are selected and the dimensional analysis is performed, the experimenter seeks to achieve similarity between the model tested and the prototype to be designed.

Flow conditions for a model test are completely similar if all relevant dimensionless parameters $({\mathrm{\Pi }}_1,{\mathrm{\Pi }}_2,...,{\mathrm{\Pi }}_n)$have the same corresponding values for the model and the prototype (White):

Template:Center top${\mathrm{\Pi }}_{1model}={\mathrm{\Pi }}_{1prototype},{\mathrm{\Pi }}_{2model}={\mathrm{\Pi }}_{2prototype},...,{\mathrm{\Pi }}_{nmodel}={\mathrm{\Pi }}_{nprototype}$Template:Center bottom

But this is easier said than done.

Instead of complete similarity, the engineering literature speaks of particular types of similarity:

• Geometric similarity

• Kinematic similarity

• Dynamic similarity

• Thermal similarity

Geometric similarity

A model and a prototype are geometrically similar if and only if all body dimensions in all three coordinates have the same linear-scale ratio.

All angles are preserved in geometric similarity. All flow directions are preserved. The orientations of model and prototype with respect to the surroundings must be identical.

File:Aircraft scaledown.svg

1:2 down-scaled aircraft model

Kinematic similarity

The motions of two systems are kinematically similar if homologous particles lie at homologous points at homologous times.

In other words, kinematic similarity requires that the model and prototype have dependent length-scale and time-scale ratios. This requires, for example, similarity in Re or Mach number. The result is that the velocity-scale ratio will be the same for both:

One exception, in which the length and the time-scales becomes independent is the frictionless low-speed flow, without any free-surface.

File:Similarity considerations2 2.png

Frictionless low speed flow

$D_{model}=\alpha D_{prototype}$

and

$V_{model}=\beta V_{prototype}$

For a frictionless two-phase flow, the Froude numbers should be same, thus

From:

Template:Center top$F{r}_{model}=\frac{V_{model}^2}{g{L}_{model}}=\frac{V_{prototype}^2}{g{L}_{prototype}}=F{r}_{prototype}$Template:Center bottom

and,

Template:Center top$L_{model}=\alpha L_{prototype}\to V_{model}=\sqrt{\alpha }{V}_{prototype}$Template:Center bottom

we obtain:

Template:Center top$\frac{L_{model}}{L_{protorype}}=\alpha$Template:Center bottom

and,

Template:Center top$\frac{T_{model}}{T_{prototype}}=\frac{L_{model}/V_{model}}{L_{prototype}/V_{prototype}}=\sqrt{\alpha }$Template:Center bottom

File:Similarity consideration 03.svg

Similarity consideration for wave effect modelling

Dynamic similarity

Dynamic similarity exists when the model and the prototype have the same length scale ratio, time-scale ratio, and force-scale (or mass-scale) ratio.

Mathematically, Newton’s law for any fluid particle requires that the sum of the pressure force, gravity force, and friction force equal the acceleration term, or inertia force:

Template:Center top$F_{inertia}=F_{body}+F_{pressure}+F_{viscous}$Template:Center bottom

The dynamic-similarity ensure that each of these forces will be in the same ratio and have equivalent directions between model and prototype.

Again geometric similarity is a first requirement. Then dynamic similarity exists, simultaneous with kinematic similarity, if the model and prototype force and pressure coefficients are identical. This is ensured if:

1. For compressible flow, the model and prototype Reynolds number and Mach number and specific-heat ratio are correspondingly equal.
2. For incompressible flow
   • With no free surface: model and prototype Reynolds numbers are equal.
   • With a free surface: model and prototype Reynolds number, Froude number, and (if necessary) Weber number and cavitation number are correspondingly equal.

Incomplete similarity

It is not possible to reach complete similarity in the laboratory for the wavy water flow around the ship

Not always, it is easy to obtain complete dynamic similarity.

Consider the drag force on the surface of the ship shown. To have complete dynamic similarity between test model and prototype:

$F{r}_{model}=F{r}_{prototype},R{e}_{model}=R{e}_{prototype}$

For a 100 times downscaled test model, the Fr number similarity requires:

Template:Center top$L_{model}=\frac{1}{100}{L}_{prototype}\to \frac{100V_{model}^2}{g{L}_{prototype}}=\frac{V_{prototype}^2}{g{L}_{prototype}}\to V_{model}=10V_{prototype}$Template:Center bottom

Thus, the Re number similarity can only be obtained at the same time

Template:Center top$R{e}_{model}=\frac{V_{model}{L}_{model}}{{\upsilon }_{model}}=\frac{10V_{prototype}{L}_{prototype}}{{\upsilon }_{model}100}=\frac{U_{prototype}{L}_{prototype}}{{\upsilon }_{prototype}}=R{e}_{prototype}$Template:Center bottom

$\to {\upsilon }_{model}=0.1{\upsilon }_{prototype}$

Only liquid which has less kinematic viscosity than water is mercury, and even with that fluid Re similarity can not be established. Since Froude number becomes more important than Re, in the experiments scaling is done using only Fr number.

An example study: Inhaler development

File:Inhaler development.png

Flow Visualization in Inhalers: Inhaler in the Water Tunnel

File:Inhaler development2.png

Experimental setup

The flow in the Handihaler is unsteady. Thus, visualization is preferred in order to understand the flow character within the inhaler. The up-scaled models in water are used to slow down the speed of the flow so that we can record streaklines with a high-speed camera.

Flow in the Handihaler® and Respimat®

File:Inhaler development3.png

After preliminary visualization experiments, scenarios were planned. The flow region was divided into sub-regions. Dye injection holes were built in these sub-regions. Around 100 visualizations were made for original and modified geometries of Handihaler and Respimat.

The mechanism of particle discharge from the capsule is discovered.

File:Inhaler development4.png

With the help of the prototypes we could easily study the effect of different flow and geometrical parameters on the unsteady flow within the Handihaler

The visualization is used to validate numerical simulations.

Experiments versus simulations

Complete design of inhalers can be made by the utilization of non-dimensional numbers.

File:Exp vs sim.png

Correlation between experiments and simulations

Conclusions

• Dimensional analysis and similitude are two important tools in experimental research.

• While analyzing a phenomenon one can use the dimensionless numbers to judge the importance of the governing effects.

• Similarity considerations should be carefully performed to obtain useful results, especially when similarity is not complete.