The velocity and shear stress versus radial position are obtained for the laminar flow of a power-law fluid in a pipe. Pipe radius and applied pressure gradients can be set by the user. If you choose a power-law exponent, n, equal to 1 then a Newtonian fluid is recovered. Dilatant and pseudo-plastic fluids are obtained for n>1 and n<1, respectively. A non-Newtonian fluid has a viscosity that changes with the applied shear force.
For a Newtonian fluid (such as water), the viscosity is independent of how fast you are stirring it, but for a non-Newtonian fluid the viscosity is dependent. It gets easier or harder to stir faster for different types of non-Newtonian fluids. Different constitutive equations, giving rise to various models of non-Newtonian fluids, have been proposed in order to express the viscosity as a function of the strain rate. In power-law fluids, n is the power-law exponent and kappa is the power-law consistency index. Dilatant fluids correspond to the case where the exponent in the power-law constitutive equation is positive while pseudo-plastic fluids are obtained when n<1. We see that viscosity decreases with strain rate for n<1, which is the case for pseudo-plastic fluids, also called shear-thinning fluids. On the other hand, dilatant fluids are shear-thickening. If n=1, one recovers the Newtonian fluid behavior.
The user of the GUI application should be aware that good data for the power law exponent and consistency index must be used in order for the result to be correct.
For a treatment using Mathematica, please visit the following links:
Housam Binous (2022). GUI application for pipe flow of a power-law fluid (https://www.mathworks.com/matlabcentral/fileexchange/17910-gui-application-for-pipe-flow-of-a-power-law-fluid), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!