Grashof number explained

In fluid mechanics (especially fluid thermodynamics), the Grashof number (after Franz Grashof) is a dimensionless number which approximates the ratio of the buoyancy to viscous forces acting on a fluid. It frequently arises in the study of situations involving natural convection and is analogous to the Reynolds number .^[1]

Definition

Heat transfer

Free convection is caused by a change in density of a fluid due to a temperature change or gradient. Usually the density decreases due to an increase in temperature and causes the fluid to rise. This motion is caused by the buoyancy force. The major force that resists the motion is the viscous force. The Grashof number is a way to quantify the opposing forces.^[2]

The Grashof number is:

Gr _{L} = (g \beta (T _{s} - T _{infty} ) L ^{3})/(\nu ^{2}) for vertical flat platesGr _{D} = (g \beta (T _{s} - T _{infty} ) D ^{3})/(\nu ^{2}) for pipes and bluff bodies

where:

• is gravitational acceleration due to Earth
• is the coefficient of volume expansion (equal to approximately for ideal gases)
• is the surface temperature
• is the bulk temperature
• is the vertical length
• is the diameter
• is the kinematic viscosity.

The and subscripts indicate the length scale basis for the Grashof number.

The transition to turbulent flow occurs in the range for natural convection from vertical flat plates. At higher Grashof numbers, the boundary layer is turbulent; at lower Grashof numbers, the boundary layer is laminar, that is, in the range .

Mass transfer

There is an analogous form of the Grashof number used in cases of natural convection mass transfer problems. In the case of mass transfer, natural convection is caused by concentration gradients rather than temperature gradients.

$$\backslash mathrm\_c\; =\; \backslash frac$$

where

$$\backslash beta^*\; =\; -\backslash frac\; \backslash left\; (\backslash frac\; \backslash right)\_$$

and:

• is gravitational acceleration due to Earth
• is the concentration of species at surface
• is the concentration of species in ambient medium
• is the characteristic length
• is the kinematic viscosity
• is the fluid density
• is the concentration of species
• is the temperature (constant)
• is the pressure (constant).

Relationship to other dimensionless numbers

The Rayleigh number, shown below, is a dimensionless number that characterizes convection problems in heat transfer. A critical value exists for the Rayleigh number, above which fluid motion occurs.

$$\backslash mathrm\_\; =\; \backslash mathrm\_\backslash mathrm$$

The ratio of the Grashof number to the square of the Reynolds number may be used to determine if forced or free convection may be neglected for a system, or if there's a combination of the two. This characteristic ratio is known as the Richardson number . If the ratio is much less than one, then free convection may be ignored. If the ratio is much greater than one, forced convection may be ignored. Otherwise, the regime is combined forced and free convection.

Ri = (Gr)/(Re ^{2}) \gg 1 \implies ignoreforcedconvection Ri = (Gr)/(Re ^{2}) ≈ 1 \implies combinedforcedandfreeconvection Ri = (Gr)/(Re ^{2}) \ll 1 \implies ignorefreeconvection

Derivation

The first step to deriving the Grashof number is manipulating the volume expansion coefficient,

\beta as follows.

$$\backslash beta\; =\; \backslash frac\backslash left(\backslash frac\backslash right)\_p\; =\backslash frac\backslash left(\backslash frac\backslash right)\_p$$

The

v in the equation above, which represents specific volume, is not the same as the v in the subsequent sections of this derivation, which will represent a velocity. This partial relation of the volume expansion coefficient, \beta , with respect to fluid density, \rho , given constant pressure, can be rewritten as

$$\backslash rho=\backslash rho\_0\; (1-\backslash beta\; \backslash Delta\; T)$$

where:

\rho _{0} is the bulk fluid density\rho is the boundary layer density\Delta T = (T - T _{0)} , the temperature difference between boundary layer and bulk fluid.

There are two different ways to find the Grashof number from this point. One involves the energy equation while the other incorporates the buoyant force due to the difference in density between the boundary layer and bulk fluid.

Energy equation

This discussion involving the energy equation is with respect to rotationally symmetric flow. This analysis will take into consideration the effect of gravitational acceleration on flow and heat transfer. The mathematical equations to follow apply both to rotational symmetric flow as well as two-dimensional planar flow.

$$\backslash frac(\backslash rho\; u\; r\_0^)+(\backslash rho\; v\; r\_0^)=0$$

where:

s is the rotational direction, i.e. direction parallel to the surfaceu is the tangential velocity, i.e. velocity parallel to the surfacey is the planar direction, i.e. direction normal to the surfacev is the normal velocity, i.e. velocity normal to the surfacer _{0} is the radius.

In this equation the superscript is to differentiate between rotationally symmetric flow from planar flow. The following characteristics of this equation hold true.

n = 1: rotationally symmetric flown = 0: planar, two-dimensional flowg is gravitational acceleration

This equation expands to the following with the addition of physical fluid properties:

$$\backslash rho\backslash left(u\; \backslash frac\; +\; v\; \backslash frac\backslash right)\; =\; \backslash frac\backslash left(\backslash mu\; \backslash frac\backslash right)\; -\; \backslash frac\; +\; \backslash rho\; g.$$

From here we can further simplify the momentum equation by setting the bulk fluid velocity to 0 (

u = 0 ).

$$\backslash frac=\backslash rho\_0\; g$$

This relation shows that the pressure gradient is simply a product of the bulk fluid density and the gravitational acceleration. The next step is to plug in the pressure gradient into the momentum equation.

$$\backslash left(u\; \backslash frac+v\; \backslash frac\backslash right)=\backslash nu\; \backslash left(\backslash frac\backslash right)+g\backslash frac=\backslash nu\; \backslash left(\backslash frac\backslash right)-\backslash frac\; g\; \backslash beta\; (T-T\_0)$$

where the volume expansion coefficient to density relationship

\rho - \rho _{0} = - \rho _{0} \beta (T - T _{0)} found above and the kinematic viscosity relationship \nu = (\mu)/(\rho) were substituted into the momentum equation.$$u\backslash left(\backslash frac\backslash right)+v\; \backslash left(\backslash frac\backslash right)=\backslash nu\; \backslash left(\backslash frac\backslash right)-\backslash fracg\; \backslash beta(T\; -\; T\_0)$$

To find the Grashof number from this point, the preceding equation must be non-dimensionalized. This means that every variable in the equation should have no dimension and should instead be a ratio characteristic to the geometry and setup of the problem. This is done by dividing each variable by corresponding constant quantities. Lengths are divided by a characteristic length,

L _{c} . Velocities are divided by appropriate reference velocities, V , which, considering the Reynolds number, gives V = (Re _{L} \nu)/(L _{c}) . Temperatures are divided by the appropriate temperature difference, (T _{s} - T _{0)} . These dimensionless parameters look like the following:s ^{* = (s)/(L _{c})} ,y ^{*} = (y)/(L _{c}) ,u ^{* = (u)/(V)} ,v ^{*} = (v)/(V) , andT ^{* = ((T - T _{0)})/((T _{s} - T _{0)})} .

The asterisks represent dimensionless parameter. Combining these dimensionless equations with the momentum equations gives the following simplified equation.

u ^{*} (\partial u ^{*})/(\partial s ^{*}) + v ^{*} (\partial u ^{*})/(\partial y ^{*}) =- \left[ (\rho _{0} g \beta(T _{s} - T _{0)L} ^{3} ^{c})/(\rho\nu ^{2} Re ^{2} _{L}) \right] T ^{* + (1)/(Re _{L})} (\partial ^{2} u ^{*})/(\partial {y ^{*}) ^{2}} =- \left( (\rho _{0})/(\right)\left[ \rho) (g \beta(T _{s} - T _{0)L} ^{3} ^{c})/(\nu ^{2}) \right] (T ^{*})/(+ Re ^{2} _{L}) (1)/(Re _{L}) (\partial ^{2} u ^{*})/(\partial {y ^{*}) ^{2}}

where:

T _{s} is the surface temperatureT _{0} is the bulk fluid temperatureL _{c} is the characteristic length.

The dimensionless parameter enclosed in the brackets in the preceding equation is known as the Grashof number:

Gr = (g \beta(T _{s-T} _{0)L} ^{3} ^{c})/(\nu ^{2}) .

Buckingham π theorem

Another form of dimensional analysis that will result in the Grashof number is known as the Buckingham π theorem. This method takes into account the buoyancy force per unit volume,

F _{b} due to the density difference in the boundary layer and the bulk fluid.

$$F\_b\; =\; (\backslash rho\; -\; \backslash rho\_0)\; g$$

This equation can be manipulated to give,

$$F\_b\; =\; -\backslash beta\; g\; \backslash rho\_0\; \backslash Delta\; T.$$

The list of variables that are used in the Buckingham π method is listed below, along with their symbols and dimensions.

Variable	Symbol	Dimensions
Significant length	L	L
Fluid viscosity	\mu	(M)/(Lt)
Fluid heat capacity	c _{p}	(Q)/(MT)
Fluid thermal conductivity	k	(Q)/(LtT)
Volume expansion coefficient	\beta	(1)/(T)
Gravitational acceleration	g	(L)/(t ^{2})
Temperature difference	\Delta T	T
Heat transfer coefficient	h	(Q)/(L ^{2} tT)

With reference to the Buckingham π theorem there are dimensionless groups. Choose,

\mu, , and \beta as the reference variables. Thus the \pi groups are as follows:\pi _{1} = L ^{a} \mu ^{b} k ^{c} \beta ^{d} g ^{e} c _{p} ,\pi _{2} = L ^{f} \mu ^{g} k ^{h} \beta ^{i} g ^{j} \rho ,\pi _{3} = L ^{k} \mu ^{l} k ^{m} \beta ^{n} g ^{o} \Delta T ,\pi _{4} = L ^{q} \mu ^{r} k ^{s} \beta ^{t} g ^{u} h .

Solving these

\pi groups gives:\pi _{1} = (\mu(c _{p)})/(k) = Pr ,\pi _{2} = (l ^{3} g \rho ^{2})/(\mu ^{2}) ,\pi _{3} = \beta \Delta T ,\pi _{4} = (h L)/(k) = Nu

From the two groups

\pi _{2} and \pi _{3,} the product forms the Grashof number:\pi _{2} \pi _{3= (\beta g \rho ^{2} \Delta T L ^{3})/(\mu ^{2})} = Gr .

Taking

\nu = (\mu)/(\rho) and \Delta T = (T _{s} - T _{0)} the preceding equation can be rendered as the same result from deriving the Grashof number from the energy equation.Gr = (\beta g \Delta T L ^{3})/(\nu ^{2})

In forced convection the Reynolds number governs the fluid flow. But, in natural convection the Grashof number is the dimensionless parameter that governs the fluid flow. Using the energy equation and the buoyant force combined with dimensional analysis provides two different ways to derive the Grashof number.

Physical Reasoning

It is also possible to derive the Grashof number by physical definition of the number as follows:

$$\backslash mathrm\; =\; \backslash frac\; =\backslash frac=\backslash frac=\backslash frac$$

However, above expression, especially the final part at the right hand side, is slightly different from Grashof number appearing in literature. Following dimensionally correct scale in terms of dynamic viscosity can be used to have the final form.

$$\backslash mathrm\; =\; \backslash rho\; V\; L$$Writing above scale in Gr gives;

$$\backslash mathrm\; =\; \backslash frac$$Physical reasoning is helpful to grasp the meaning of the number. On the other hand, following velocity definition can be used as a characteristic velocity value for making certain velocities nondimensional.

$$\backslash mathrm\; =\backslash frac$$

Effects of Grashof number on the flow of different fluids

In a recent research carried out on the effects of Grashof number on the flow of different fluids driven by convection over various surfaces.^[3] Using slope of the linear regression line through data points, it is concluded that increase in the value of Grashof number or any buoyancy related parameter implies an increase in the wall temperature and this makes the bond(s) between the fluid to become weaker, strength of the internal friction to decrease, the gravity to becomes stronger enough (i.e. makes the specific weight appreciably different between the immediate fluid layers adjacent to the wall). The effects of buoyancy parameter are highly significant in the laminar flow within the boundary layer formed on a vertically moving cylinder. This is only achievable when the prescribed surface temperature (PST) and prescribed wall heat flux (WHF) are considered. It can be concluded that buoyancy parameter has a negligible positive effect on the local Nusselt number. This is only true when the magnitude of Prandtl number is small or prescribed wall heat flux (WHF) is considered. Sherwood number, Bejan Number, Entropy generation, Stanton Number and pressure gradient are increasing properties of buoyancy related parameter while concentration profiles, frictional force, and motile microorganism are decreasing properties.

Further reading

• Book: Cengel, Yunus A.. Heat and Mass Transfer: A Practical Approach. 3rd . Boston. McGraw Hill. 2003.
• Book: Ernst R. G. Eckert. Eckert. Ernst R. G.. Drake. Robert M.. Analysis of Heat and Mass Transfer. New York. McGraw Hill. 1972.
• Book: Jaluria, Yogesh. Natural Convection Heat and Mass Transfer. New York: Pergamon Press. 1980.
• Book: Welty, James R.. Fundamentals of Momentum, Heat, and Mass Transfer. New York. John Wiley & Sons. 1976.

Notes and References

1. Book: Incropera, Frank. Fundamentals of Heat and Mass Transfer. limited. Wiley. 2007. Hoboken, NJ. 408, 599, 629. 6th. 288958608. 9780471457282.
2. Book: Bird. R. Byron. Warren E.. Stewart. Edwin N.. Lightfoot. Transport Phenomena. limited. New York. J. Wiley. 2002. 318, 359. 2nd. 471520548. 9780471410775.
3. Shah. Nehad Ali. Animasaun. I.L.. Ibraheem. R.O.. Babatunde. H.A.. Sandeep. N.. Pop. I.. Scrutinization of the effects of Grashof number on the flow of different fluids driven by convection over various surfaces. Journal of Molecular Liquids. 249. 2018. 980–990. 0167-7322. 10.1016/j.molliq.2017.11.042.