Non-Newtonian Flow
For the Newtonian fluid, the viscous stress tensor is proportional to the strain rate tensor, and the constitutive relation is written as:
using the compressible and incompressible formulations. Here denotes the strain-rate tensor defined by:
Its magnitude, the shear rate, is:
where the contraction operator “:” is defined by
For non-Newtonian fluids, viscous stress versus shear rate curve is nonlinear or does not pass through the origin. The Laminar Flow interfaces provides five predefined inelastic non-Newtonian constitutive models: Power law, Carreau, Bingham-Papanastasiou, Herschel-Bukley-Papanastasiou, and Casson-Papanastasiou.
Power Law
The Power law model is an example of a generalized Newtonian model. It prescribes
(4-18)
where m and n are scalars that can be set to arbitrary values. For n > 1, the power law describes a shear thickening (dilatant) fluid. For n < 1, it describes a shear thinning (pseudoplastic) fluid. A value of n equal to one gives the expression for a Newtonian fluid.
Equation 4-18 predicts an infinite viscosity at zero shear rate for n < 1. This is however never the case physically. Instead, most fluids have a constant viscosity for shear rates smaller than 102 s1 (Ref. 21). Since infinite viscosity also makes models using Equation 4-18 difficult to solve, COMSOL Multiphysics implements the Power law model as
(4-19)
where is a lower limit for the evaluation of the shear rate magnitude. The default value for is 102 s1, but can be given an arbitrary value or expression using the corresponding text field.
Carreau Model
The Carreau model defines the viscosity in terms of the following four-parameter expression
(4-20)
where λ is a parameter with the unit of time, μ0 is the zero shear rate viscosity, μinf is the infinite shear-rate viscosity, and n is a dimensionless parameter. This expression is able to describe the viscosity for most stationary polymer flows.
BINGHAM-PAPANASTASIOU
Viscoplastic fluid behavior is characterized by existence of the yield stress τy — a limit which must be exceeded before significant deformation can occur. To model the stress-deformation behavior of viscoplastic materials, different constitutive equations have been propose. The Bingham plastic model is written as
where μp is the plastic viscosity.
To allow computation in both yielded and unyielded region, the Papanastasiou continuous regularization for the viscosity function is used:
(4-21)
where mp is a scale that controls the exponential growth of stress.
HERSCHEL-BULKLEY-PAPANASTASIOU
The Herschel-Bulkley -Papanastasiou model combines the effects of the Power law and Bingham behavior:
(4-22)
where m and n are scalars that can be set to arbitrary values. For a value of n equal to one, the Bingham-Papanastasiou model is recovered. This equation is valid in both yielded and unyielded regions. The exponent mp controls the smoothness of the viscosity function.
CASSON-PAPANASTASIOU
The Casson-Papanastasiou model combines Casson eqaution with Papanastasiou regularization:
(4-23)
The equation is valid in both yielded and unyielded regions.