0% found this document useful (0 votes)
17 views2 pages

Power Law and Bingham Fluid Flow Analysis

This tutorial focuses on solving fluid dynamics problems involving power law and Bingham fluids over an inclined plane and in an annulus between concentric cylinders. It includes deriving velocity profiles, volumetric flow rates, and relationships involving yield stress, inclination angles, and film thickness. The tutorial emphasizes the use of Cauchy equations and the assumption of negligible fluid inertia in the analysis.

Uploaded by

parthrajput.iit
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
17 views2 pages

Power Law and Bingham Fluid Flow Analysis

This tutorial focuses on solving fluid dynamics problems involving power law and Bingham fluids over an inclined plane and in an annulus between concentric cylinders. It includes deriving velocity profiles, volumetric flow rates, and relationships involving yield stress, inclination angles, and film thickness. The tutorial emphasizes the use of Cauchy equations and the assumption of negligible fluid inertia in the analysis.

Uploaded by

parthrajput.iit
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Tutorial 8 - CLL110B

1. In tutorial 7a we solved the problem of a steady falling film of a Newtonian viscous


incompressible fluid over an inclined plane using the equations of motions. We now
consider the same problem for a power law fluid. In this case, however, the starting
point is the Cauchy equations:
𝜕𝒗
𝜌( + 𝒗 ⋅ 𝛁𝒗) = −𝛁 ⋅ 𝛕 + 𝜌𝒈
𝜕𝑡
You may use the below schematic as a reference. Consider a coordinate system
such that the positive 𝑥 -axis is along the downward incline, 𝑦-axis is orthogonal to it
with the incline plane at 𝑦 = 𝛿 and the free surface of the liquid at 𝑦 = 0. Solve for
only the 𝑥 momentum balance and obtain the fluid velocity profile when fluid inertia is
negligible. Given the 2D nature of the problem, you may forgo any dependence on 𝑧-
axis. A power law fluid has the following relationship between the stress and the
𝑑𝑣 𝑑𝑣 𝑛−1
velocity gradient: 𝜏 = −𝜇(𝛾̇ ) 𝑑 𝑦𝑥 , where 𝜇(𝛾̇ ) = 𝑚 | 𝑑 𝑦𝑥 | . Obtain the
expression for the volumetric flow rate of the fluid.

2. Solve the above question for a Bhingham fluid, which obeys the constitutive
equation:
𝑑 𝑣𝑥 𝑑𝑣
𝜏 = −𝜇(𝛾̇ ) , where 𝜇(𝛾̇ ) = ∞ (𝑜𝑟 𝑑 𝑦𝑥 = 0) ; 𝜏 ≤ 𝜏0 𝑎𝑛𝑑 𝜇(𝛾̇ ) = 𝜇0 +
𝑑𝑦
𝜏0
𝑑𝑣𝑥 ; 𝜏 ≥ 𝜏0 .
|𝑑 𝑦|

Here, 𝜏0 is the yield stress specified for a Bingham fluid.


𝑑 𝑣𝑥
a) For the choice of the coordinate system, is
𝑑𝑦
> 0 or < 0?
b) Obtain the velocity profile when the fluid exhibits yielding character at a
2𝛿
height 𝐻 = above the inclined plane.
3
c) How is the yield stress 𝜏0 related to inclination angle 𝛼 and film thickness
𝛿?
d) What is the maximum velocity of the fluid, and where?
e) Plot the velocity profile.
N.B. Use the same coordinate system as used in the above problem.

3. In the flow through annulus between two vertically oriented concentric cylinders,
where the outer cylinder is rotating with an angular velocity Ω and inner cylinder
is stationary, assume that the intervening fluid is a power-law fluid. Obtain the
velocity profile, assuming that fluid inertia is negligible – that is the left-hand side
of the above written Cauchy equations is 0. Assume that 𝑣𝑟 = 𝑣𝑧 = 0, pressure
𝑃 ≠ 𝑃(𝜃), and 𝑣𝜃 = 𝑣𝜃 (𝑟). You’ll first need to write down the relevant stress
equation from appendix B5. In what limit do you get the expression for the
Newtonian fluid – show it explicitly?

Common questions

Powered by AI

The velocity profile for a power law fluid in an annulus between two vertically oriented concentric cylinders where the outer cylinder rotates is derived using the stress balance equation from the modified Cauchy equation, where fluid inertia is negligible \(( \alpha = 0)\). Assuming \( v_r = v_z = 0 \) and \( v_\theta = v_\theta(r) \), and that the pressure \( P \) is not a function of \( \theta \), the stress relationship \( \tau = -\mu(\dot{\gamma}) \frac{dv_\theta}{dr} \) is utilized. The integration yields \( v_\theta(r) = C_1 \times r^{(n+1)/n} + C_2 \), with constants determined by boundary conditions. The Newtonian fluid limit is obtained as \( n \to 1 \), resulting in a linear velocity profile \( v_\theta(r) = C_3 \times \ln(r) + C_4 \), reflecting the classical solution.

The maximum velocity for a Bingham fluid flowing down an inclined plane occurs at the free surface of the fluid layer, where shearing is at a minimum. The velocity profile derived from integrating the Bingham fluid constitutive equations illustrates that as height \( y \) increases towards the free surface (\( y = 0 \)), shear effects reduce, achieving maximum velocity. Calculating this involves taking the derived velocity equation \( v_x(y) = \frac{\tau_0}{\mu_0} (y - H) + \frac{\rho g \sin(\alpha)}{\mu_0} (\frac{y^2}{2}) \) and evaluating it at \( y = 0 \), considering the yield stress \( \tau_0 \) and other parameters \( \mu_0, \alpha, H \).

The yield stress \( \tau_0 \) for a Bingham fluid flowing down an inclined plane is related to the inclination angle \( \alpha \) and the film thickness \( \delta \) by considering the balance between gravitational forces and the yield condition. The critical yield stress is \( \tau_0 = \rho g \delta \sin(\alpha) \). This expression implies that the yield stress is directly proportional to the fluid density \( \rho \), gravitational acceleration \( g \), film thickness \( \delta \), and the sine of the inclination angle \( \alpha \), reflecting the dependence on both gravitational forces and geometric configuration.

The velocity profile for a power law fluid flowing down an inclined plane, neglecting fluid inertia, is derived using the x-momentum balance from the Cauchy equation: \( \rho( \partial v / \partial t + v \cdot \nabla v) = - \nabla \cdot \tau + \rho g \). Assuming a steady state and negligible inertia, this simplifies to \( - \nabla \cdot \tau + \rho g = 0 \). For a power law fluid, the stress \( \tau \) is related to the velocity gradient by \( \tau = -\mu( \dot{\gamma} ) \frac{dv_x}{dy} \), where \( \mu(\dot{\gamma}) = m| \frac{dv_x}{dy}|^{n-1} \). Integrating the resulting equation and applying boundary conditions at the free surface and film thickness, the velocity profile is obtained as \( v_x(y) = A \cdot (y/\delta)^n \), where \( A \) is a constant derived from the balance equations.

To calculate the velocity profile of a power law fluid in the annular flow between rotating cylinders, assumptions include negligible fluid inertia (\( \alpha = 0 \)), axisymmetric flow conditions (\( v_r = 0, v_z = 0 \)), and constant pressure (\( P \neq P(\theta) \)). These assumptions simplify the Cauchy equations, enabling the focus on shear effects in \( v_\theta = v_\theta(r) \) caused by cylinder rotation. Justification stems from the dominance of shearing over inertial forces in narrow gaps and slow rotational speeds, where viscosity overcomes inertial effects. These approximations serve analytical tractability while preserving core flow characteristics, giving useful insights for practical applications like polymer flow in similar geometries.

The choice of the coordinate system in fluid dynamics affects both the mathematical expression and the physical interpretation of velocity gradients. In the case of Bingham fluids flowing down an inclined plane, setting the positive x-axis along the incline ensures that the velocity gradient \( \frac{dv_x}{dy} \) aligns with the direction of gravity-driven flow. Conventionally, this results in negative gradients as the velocity decreases orthogonally to the incline plane from air-fluid interface downwards. If \( \frac{dv_x}{dy} > 0 \), it indicates a non-physical situation within the given setup, disrupting the correct mathematical application and interpretation of stress and shear within the governing equations.

The velocity profile for a Bingham fluid flowing down an inclined plane can be graphically represented by plotting \( v_x(y) \) against \( y \), illustrating the piecewise nature of Bingham fluid flow. Below the yield height \( H \), the fluid behaves as a solid \( v_x(y)=0 \), while above, it flows with the velocity increasing towards the free surface. Such a plot highlights regions of different flow behavior and clarifies the marked transition across the yield point, substantiating analytical predictions. This graphical depiction aids in visualizing the impact of viscosity, yield stress and inclination on flow, revealing critical parameters that control flow regimes.

The volumetric flow rate \( Q \) for a power law fluid flowing over an inclined plane is derived by integrating the velocity profile along the depth of the fluid. Given the velocity profile \( v_x(y) = A (y/\delta)^n \), the volumetric flow rate is \( Q = \int_0^\delta v_x(y) \ dy \). Performing the integration results in \( Q = A \delta \frac{\delta^{n}}{n+1} \). Here, \( A \) is a constant derived from the momentum balance equations, and \( n\) is the flow behavior index of the power law fluid.

Bingham fluids differ from power law fluids in that they exhibit a yield stress \( \tau_0 \) below which the material behaves as a rigid body. When flow occurs, the constitutive equation for Bingham fluids is \( \tau = \mu( \dot{\gamma} ) \frac{dv_x}{dy} \), with \( \mu(\dot{\gamma}) = \mu_0 + \tau_0/| \frac{dv_x}{dy}| \) for \( \tau \geq \tau_0 \). For a Bingham fluid flowing down an incline, the velocity profile considering yield stress is obtained by integrating the stress equation. For \( \tau > \tau_0 \), the fluid flows, and the velocity profile is \( v_x(y) = \frac{\tau_0}{\mu_0} (y - H) + \frac{\rho g \sin(\alpha)}{\mu_0} (\frac{y^2}{2}) \), where \( H = \frac{2\delta}{3} \) is the height above which the fluid yields.

The stress equation for a power-law fluid in concentric annular flows states \( \tau = -\mu(\dot{\gamma}) \frac{dv_\theta}{dr} \) with \( \mu(\dot{\gamma}) = m| \frac{dv_\theta}{dr}|^{n-1} \). By integrating this relation, the velocity profile is \( v_\theta(r) = C_1 \cdot r^{1+(n-1)/n} + C_2 \). The Newtonian fluid limit is achieved at \( n = 1 \), simplifying \( \mu(\dot{\gamma}) \) to the constant dynamic viscosity \( \eta \), and the velocity profile becomes linear, \( v_\theta(r) = C_3 \cdot \ln(r) + C_4 \). The explicit limit is verified by substituting \( n = 1 \) into the constitutive and velocity equations, verifying that they reduce to classical Newtonian linear profiles familiar in laminar flow between rotating surfaces.

You might also like