Power Law and Bingham Fluid Flow Analysis
Power Law and Bingham Fluid Flow Analysis
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.