shear_estimate_2d Function

private pure function shear_estimate_2d(nvrstk, mean, deviation, grad) result(maxshear)

Estimate the magnitude of shear terms.

We use an estimate for mu(dv/dx) from the deviations and derivative estimates of the conservative variables: dm/dx = dv/dx * rho + drho/dx * v dv/dx = (dm/dx - drho/dx * v)/rho dv/dx < (max(dm/dx) - max(drho/dx) * max(v))/min(rho) dv/dx < (max(dm/dx) * min(rho) - max(drho/dx) * max(m) ) / min(rho)*2

Arguments

Type	Intent	Optional	Attributes		Name
class(atl_navierStokes_type),	intent(in)			::	nvrstk	Description of the equation
real(kind=rk),	intent(in)			::	mean(:)	The mean of each state variable.
real(kind=rk),	intent(in)			::	deviation(:)	Estimation of maximal deviation of each state.
real(kind=rk),	intent(in)			::	grad(:)	Estimation of maximal gradient of each state.

Return Value real(kind=rk)

Resulting estimate for maximal shear.