PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ 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::SSAFD::compute_nuH_staggered(), pism::stressbalance::SSAFD::compute_nuH_staggered_cfbc(), and pism::stressbalance::SSAFEM::PointwiseNuHAndBeta().