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

◆ D2()

static double pism::stressbalance::D2 ( double  u_x,
double  u_y,
double  u_z,
double  v_x,
double  v_y,
double  v_z 
)
inlinestatic

This function computes \(D^2\) defined by

\[ 2D^2 = D_{ij} D_{ij}\]

or

\[ D^2 = \frac{1}{2}\,\left(\frac{1}{2}\,(v_{z})^2 + (v_{y} + u_{x})^2 + (v_{y})^2 + \frac{1}{2}\,(v_{x} + u_{y})^2 + \frac{1}{2}\,(u_{z})^2 + (u_{x})^2\right) \]

(note the use of the summation convention). Here \(D_{ij}\) is the strain rate tensor. See StressBalance::compute_volumetric_strain_heating() for details.

Parameters
u_x,u_y,u_zpartial derivatives of \(u\), the x-component of the ice velocity
v_x,v_y,v_zpartial derivatives of \(v\), the y-component of the ice velocity
Returns
\(D^2\), where \(D\) is defined above.

Definition at line 440 of file StressBalance.cc.

Referenced by pism::Time::calendar_year_start(), pism::stressbalance::StressBalance::compute_volumetric_strain_heating(), and pism::Time::year_fraction().