Viscous Shear Behavior
The resistance to flow of a viscous fluid is described by the following relationship between deviatoric stress and strain rate
where is the deviatoric stress, is the deviatoric part of the strain rate, is the viscosity, and is the engineering shear strain rate.
Newtonian fluids are characterized by a viscosity that only depends on temperature, . In the more general case of non-Newtonian fluids the viscosity is a function of the temperature and shear strain rate:
where is the equivalent shear strain rate. In terms of the equivalent shear stress, , we have:
Non-Newtonian fluids can be classified as shear-thinning (or pseudoplastic), when the apparent viscosity decreases with increasing shear rate, and shear-thickening (or dilatant), when the viscosity increases with strain rate.
In addition to the Newtonian viscous fluid model, Abaqus/Explicit supports several models of nonlinear viscosity to describe non-Newtonian fluids: power law, Carreau-Yasuda, Cross, Herschel-Bulkley, Powell-Eyring, and Ellis-Meter. Other functional forms of the viscosity can also be specified in tabular format. In addition, in Abaqus/Explicit user subroutine VUVISCOSITY can be used.
In addition to the Newtonian viscous fluid model, Abaqus/Standard supports several models of nonlinear viscosity to describe non-Newtonian fluids for use with fluid pipe elements and fluid pipe connector elements: power law, Herschel-Bulkley, and Bingham plastic.
Newtonian
The Newtonian model is useful to model viscous laminar flow governed by the Navier-Poisson law of a Newtonian fluid, . Newtonian fluids are characterized by a viscosity that depends only on temperature, . You need to specify the viscosity as a tabular function of temperature when you define the Newtonian viscous deviatoric behavior.
In Abaqus/Standard this model is used to define the fluid behavior in fluid pipe and fluid pipe connector elements.
Power Law
The power law model is commonly used to describe the viscosity of non-Newtonian fluids. The viscosity is expressed as
where is the flow consistency index and is the flow behavior index. When , the fluid is shear-thinning (or pseudoplastic): the apparent viscosity decreases with increasing shear rate. When , the fluid is shear-thickening (or dilatant); and when , the fluid is Newtonian. Optionally, you can place a lower limit, , and/or an upper limit, , on the value of the viscosity computed from the power law.
In Abaqus/Standard this model is used to define the fluid behavior in fluid pipe and fluid pipe connector elements.
Carreau-Yasuda
The Carreau-Yasuda model describes the shear thinning behavior of polymers and is available only in Abaqus/Explicit. This model often provides a better fit than the power law model for both high and low shear strain rates. The viscosity is expressed as
where is the low shear rate Newtonian viscosity, is the infinite shear viscosity (at high shear strain rates), is the natural time constant of the fluid ( is the critical shear rate at which the fluid changes from Newtonian to power law behavior), and represents the flow behavior index in the power law regime. The coefficient is a material parameter. The original Carreau model is recovered when =2.
Cross
The Cross model is commonly used when it is necessary to describe the low-shear-rate behavior of the viscosity and is available only in Abaqus/Explicit. The viscosity is expressed as
where is the Newtonian viscosity, is the infinite shear viscosity (usually assumed to be zero for the Cross model), is the natural time constant of the fluid ( is the critical shear rate at which the fluid changes from Newtonian to power-law behavior), and is the flow behavior index in the power law regime.
Herschel-Bulkley
The Herschel-Bulkley model can be used to describe the behavior of viscoplastic fluids that exhibit a yield response. The viscosity is expressed as
Here is the “yield” stress and is a penalty viscosity to model the “rigid-like” behavior in the very low strain rate regime ( ), when the stress is below the yield stress, . With increasing strain rates, the viscosity transitions into a power law model once the yield threshold is reached, . The parameters and are the flow consistency and the flow behavior indexes in the power law regime, respectively.
In Abaqus/Standard this model is used to define the fluid behavior in fluid pipe and fluid pipe connector elements.
Bingham Plastic
The Bingham plastic model is a special case of the Herschel-Bulkley model. Bingham plastics correspond to the case where the flow behavior index in the power law regime, , is equal to one.
This model is available only in Abaqus/Standard for fluid pipe and fluid pipe connector elements.
Powell-Eyring
The Powell-Eyring model, which is derived from the theory of rate processes, is relevant primarily to molecular fluids but can be used in some cases to describe the viscous behavior of polymer solutions and viscoelastic suspensions over a wide range of shear rates. This model is available only in Abaqus/Explicit. The viscosity is expressed as
where is the Newtonian viscosity, is the infinite shear viscosity, and represents a characteristic time of the measured system.
Ellis-Meter
The Ellis-Meter model, available only in Abaqus/Explicit, expresses the viscosity in terms of the effective shear stress, , as:
where is the effective shear stress at which the viscosity is 50% between the Newtonian limit, , and the infinite shear viscosity, , and represents the flow index in the power law regime.
Tabular
In Abaqus/Explicit the viscosity can be specified directly as a tabular function of shear strain rate and temperature.
User-Defined
In Abaqus/Explicit you can specify a user-defined viscosity in user subroutine VUVISCOSITY (see VUVISCOSITY).