m_hb_function Module Reference

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

Functions/Subroutines
real(wp) function, public	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	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	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.

Detailed Description

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

Function/Subroutine Documentation

f_compute_hb_viscosity()

real(wp) function, public m_hb_function::f_compute_hb_viscosity	(	real(wp), intent(in)	tau0,
		real(wp), intent(in)	k_val,
		real(wp), intent(in)	nn_val,
		real(wp), intent(in)	mu_min_val,
		real(wp), intent(in)	mu_max_val,
		real(wp), intent(in)	shear_rate,
		real(wp), intent(in)	hb_m_val )

Papanastasiou-regularized Herschel-Bulkley viscosity.

Definition at line 330 of file m_hb_function.fpp.f90.

Here is the caller graph for this function:

f_compute_shear_rate_from_components()

real(wp) function, public m_hb_function::f_compute_shear_rate_from_components	(	real(wp), intent(in)	d_xx,
		real(wp), intent(in)	d_yy,
		real(wp), intent(in)	d_zz,
		real(wp), intent(in)	d_xy,
		real(wp), intent(in)	d_xz,
		real(wp), intent(in)	d_yz )

Shear rate gamma_dot = sqrt(2 D_ij D_ij). Absent dims pass 0.

Definition at line 369 of file m_hb_function.fpp.f90.

Here is the caller graph for this function:

s_compute_mixture_inv_re()

subroutine, public m_hb_function::s_compute_mixture_inv_re	(	real(wp), dimension(*), intent(in)	alpha,
		real(wp), intent(in)	shear_rate,
		real(wp), dimension(2,*), intent(in)	res,
		real(wp), dimension(2), intent(out)	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.

Definition at line 397 of file m_hb_function.fpp.f90.

Here is the call graph for this function:

Here is the caller graph for this function: