MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_hb_function.fpp.f90 File Reference

Go to the source code of this file.

Modules

module  m_hb_function
 Herschel-Bulkley non-Newtonian viscosity: formula, shear rate, and mixture inverse-Re.

Functions/Subroutines

real(wp) function, public m_hb_function::f_compute_hb_viscosity (tau0, k_val, nn_val, mu_min_val, mu_max_val, shear_rate, hb_m_val)
 Papanastasiou-regularized Herschel-Bulkley viscosity.
real(wp) function, public m_hb_function::f_compute_shear_rate_from_components (d_xx, d_yy, d_zz, d_xy, d_xz, d_yz)
 Shear rate gamma_dot = sqrt(2 D_ij D_ij). Absent dims pass 0.
subroutine, public m_hb_function::s_compute_mixture_inv_re (alpha, shear_rate, res, re_out)
 Mixture inverse Reynolds (= 1/mu_mix) per direction (1=shear, 2=bulk) at one state. Newtonian fluids accumulate the legacy per-fluid alpha/Res term. For non-Newtonian fluids only the shear direction (i==1) uses the Herschel-Bulkley viscosity at the given shear rate; the bulk direction (i==2) contributes zero, since non-Newtonian bulk viscosity is not supported.