Loading web-font TeX/Math/Italic
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ secondInvariant_2D()

static double pism::secondInvariant_2D ( const Vector2d U_x,
const Vector2d U_y 
)
inlinestatic

This uses the definition of squared second invariant from Hutter and several others, namely the output is D^2 = \frac 1 2 D_{ij} D_{ij} where incompressibility is used to compute D_{zz}.

This is the approximation of the full second invariant corresponding to the shallow shelf approximation. In particular, we assume that u and v are depth-independent ( u_z = v_z = 0 ) and neglect horizontal derivatives of the vertical velocity ( w_x = w_y = 0 ).

Definition at line 45 of file FlowLaw.hh.

References pism::Vector2d::u, and pism::Vector2d::v.

Referenced by pism::inverse::IP_SSAHardavForwardProblem::apply_jacobian_design(), pism::inverse::IP_SSAHardavForwardProblem::apply_jacobian_design_transpose(), pism::stressbalance::compute_2D_stresses(), pism::stressbalance::SSAFDBase::compute_nuH_cfbc(), pism::stressbalance::SSAFDBase::compute_nuH_everywhere(), and pism::stressbalance::SSAFEM::PointwiseNuHAndBeta().