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.
Called by
proc~~shear_estimate_2d~~CalledByGraph
proc~shear_estimate_2d
shear_estimate_2d
proc~atl_shearestimate_getpoint
atl_shearestimate_getPoint
proc~atl_shearestimate_getpoint->proc~shear_estimate_2d
proc~atl_shearestimate_getelement
atl_shearestimate_getElement
proc~atl_shearestimate_getelement->proc~shear_estimate_2d
Nodes of different colours represent the following:
Graph Key
Subroutine
Subroutine
Function
Function
Interface
Interface
Unknown Procedure Type
Unknown Procedure Type
Program
Program
This Page's Entity
This Page's Entity
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Variables
Type Visibility Attributes Name Initial
real(kind=rk),
private ::
rho_min
real(kind=rk),
private ::
m_max
real(kind=rk),
private ::
grad_mag