Skip to content

Commit fdf5c21

Browse files
committed
fix: restore original PHYSICS_DOCS entry for non-Newtonian viscosity
1 parent 1b1f6b3 commit fdf5c21

1 file changed

Lines changed: 13 additions & 4 deletions

File tree

toolchain/mfc/case_validator.py

Lines changed: 13 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -147,10 +147,19 @@
147147
"explanation": "Reynolds numbers must be positive. Not supported with model_eqns = 1.",
148148
},
149149
"check_non_newtonian": {
150-
"title": "Non-Newtonian Viscosity",
151-
"category": "Numerical Schemes",
152-
"math": r"\mu = \frac{\tau_0}{\dot\gamma}(1 - e^{-m\dot\gamma}) + K\dot\gamma^{n-1}",
153-
"explanation": ("Herschel-Bulkley non-Newtonian viscosity model. Requires viscous=T, K > 0, nn > 0, tau0 >= 0, mu_min >= 0, mu_max > mu_min, hb_m > 0."),
150+
"title": "Non-Newtonian Viscosity (Herschel-Bulkley)",
151+
"category": "Physics Models",
152+
"math": r"\mu(\dot\gamma) = \frac{\tau_0}{\dot\gamma}(1 - e^{-m\dot\gamma}) + K\dot\gamma^{n-1}",
153+
"explanation": (
154+
"Herschel-Bulkley model with Papanastasiou regularization for "
155+
"non-Newtonian viscosity. The yield stress tau0 introduces a "
156+
"threshold below which the fluid resists deformation; the "
157+
"Papanastasiou parameter m smooths the transition (higher m = "
158+
"sharper yield). K is the consistency index, nn the flow behavior "
159+
"index (nn < 1: shear-thinning, nn = 1: Bingham, nn > 1: "
160+
"shear-thickening). Viscosity is clamped to [mu_min, mu_max] for "
161+
"numerical stability. Requires viscous = T."
162+
),
154163
},
155164
# Feature Compatibility
156165
"check_mhd": {

0 commit comments

Comments
 (0)