Regularization Model#

Regularization Model = <MODEL>

Description / Usage#

This card is used to specify the model for the regularization used in the generalized Newtonian viscosity.

Currently used for HERSCHEL_BULKLEY models

Choices for <MODEL> in HERSCHEL_BULKLEY fluids:

EPSILON

Epsilon regularization applied to both the power law term and the yield term (default for HERSCHEL_BULKLEY)

\[\mu = \mu_0 (\dot{\gamma} + \epsilon)^{n-1} + \frac{\tau_y}{(\dot{\gamma} + \epsilon)}\]
PAPANASTASIOU

Papanastasiou regularization applied to the yield parameter (default for HERSCHEL_BULKLEY_PAPANASTASIOU)

\[\mu = \mu_0 \dot{\gamma}^{n-1} + (1-exp(-f \dot{\gamma})) \frac{\tau_y}{\dot{\gamma}}\]
EPSILON_YIELD_ONLY

Epsilon regularization applied only to the yield term

\[\mu = \mu_0 (\dot{\gamma})^{n-1} + \frac{\tau_y}{(\dot{\gamma} + \epsilon)}\]
PAPANASTASIOU_EPSILON

Papanastasiou regularization applied to the yield term and an epsilon on the power law term

\[\mu = \mu_0 (\dot{\gamma}+ \epsilon)^{n-1} + (1-exp(-f \dot{\gamma})) \frac{\tau_y}{\dot{\gamma}}\]

Examples#

The following is a sample card that sets Papanastasiou regularization to a HERSCHEL_BULKLEY fluid with an epsilon on the power law term

Liquid Constitutive Equation = HERSCHEL_BULKLEY
...
Regularization Model         = PAPANASTASIOU_EPSILON
Yield Exponent               = CONSTANT 200
Epsilon Regularization       = CONSTANT 1e-12

Without the epsilon term

Liquid Constitutive Equation = HERSCHEL_BULKLEY
...
Regularization Model         = PAPANASTASIOU
Yield Exponent               = CONSTANT 200

Similarly applying only to the Yield term for epsilon regularization

Liquid Constitutive Equation = HERSCHEL_BULKLEY
...
Regularization Model         = EPSILON_YIELD_ONLY
Epsilon Regularization       = CONSTANT 1e-5

Or to both terms

Liquid Constitutive Equation = HERSCHEL_BULKLEY
...
# optional this is the default
Regularization Model         = EPSILON
Epsilon Regularization       = CONSTANT 1e-5

Technical Discussion#

See Description/Usage for Liquid Constitutive Equation