Add Herschel-Bulkley non-Newtonian viscosity model#1545

Merged

sbryngelson merged 36 commits into

MFlowCode:masterfrom

sbryngelson:feature/non-newtonian-viscosity

Jun 10, 2026

Commits on Jun 10, 2026

docs: add non-Newtonian (Herschel-Bulkley) viscosity design spec
sbryngelson
committed
docs: add non-Newtonian viscosity implementation plan
sbryngelson
committed
feat(params): register Herschel-Bulkley non-Newtonian fluid parameters
sbryngelson
committed
feat(params): default-initialize non-Newtonian fields in all targets
sbryngelson
committed
feat(params): MPI-broadcast non-Newtonian fields in all targets
sbryngelson
committed
feat(nn): add m_hb_function (HB formula, shear rate, mixture inv-Re)
sbryngelson
committed
feat(nn): add any_non_newtonian flag and GPU HB parameter arrays
sbryngelson
committed
feat(nn): shear-dependent viscosity in shared Cartesian viscous flux (HLL/HLLC)
sbryngelson
committed
test(nn): power-law Poiseuille example with analytic velocity-profile check
sbryngelson
committed
feat(nn): shear-dependent viscosity in cylindrical viscous flux and axis boundary
sbryngelson
committed
feat(nn): per-stencil-sample viscosity in IBM viscous stress tensor
sbryngelson
committed
feat(nn): bound viscous CFL by mu_max for non-Newtonian fluids
sbryngelson
committed
feat(nn): case-validator guards for non-Newtonian parameters
sbryngelson
committed
test(nn): regression cases for shear-thinning and shear-thickening fluids
sbryngelson
committed
test(nn): lid-driven cavity non-Newtonian example (Li et al. 2015)
sbryngelson
committed
docs(nn): document Herschel-Bulkley non-Newtonian viscosity
sbryngelson
committed
fix(nn): address final review — case-opt dimension guards, device-safe seq-routine dummies, doc accuracy
sbryngelson
committed
test(nn): Bingham Poiseuille example validating the yield-stress term (plug flow)
sbryngelson
committed
fix(nn): forbid unsupported mu_bulk (non-Newtonian bulk viscosity) in the case validator
sbryngelson
committed
Delete docs/superpowers/plans/2026-06-08-non-newtonian-viscosity.md
sbryngelson
committed
Discard changes to docs/superpowers/specs/2026-06-08-non-newtonian-viscosity-design.md
sbryngelson
committed
docs(nn): add READMEs to the validated non-Newtonian example cases
sbryngelson
committed
test(nn): shear-thickening and general Herschel-Bulkley Poiseuille analytic examples
sbryngelson
committed
fix(nn): validate HB parameter positivity and forbid bulk Re(2) for non-Newtonian fluids
sbryngelson
committed
refactor(nn): drop pure from m_hb_function device routines; clarify mixture docstring
sbryngelson
committed
docs(nn): fix cavity benchmark reference, HB citation, and bluntness wording
sbryngelson
committed
test(nn): golden files for non-Newtonian example regression tests
sbryngelson
committed
fix(validator): close non-Newtonian validation gaps (num_fluids default, stray HB params)
sbryngelson
committed
refactor(nn): remove dead hb_mu_bulk plumbing; clamp NN interface volume fractions
sbryngelson
committed
test(nn): fast regression cases for yield-stress branch and mixed NN/Newtonian fluids
sbryngelson
committed
docs(nn): fix cavity citation/contradiction, document validator constraints, drop dev-history notes
sbryngelson
committed
test(nn): IBM-walled power-law Poiseuille validating the IBM non-Newtonian path
sbryngelson
committed
test(nn): IBM non-Newtonian regression case and example golden
sbryngelson
committed
docs(nn): list non-Newtonian rheology in README features and equations overview
sbryngelson
committed
fix(ci): skip truncated ibm_poiseuille_nn example test (under-resolved IB slabs are platform-marginal)
sbryngelson
committed
fix(ci): actually register the ibm_poiseuille_nn example skip (missed in 419ade9ba)
sbryngelson
committed