AIDJEX Bulletin #40 - Polar Science Center - University of Washington
AIDJEX Bulletin #40 - Polar Science Center - University of Washington
AIDJEX Bulletin #40 - Polar Science Center - University of Washington
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
in which there is shearing but no extending since these lie at the origin<br />
<strong>of</strong> the horizontal axis that is a measure <strong>of</strong>,normal stretching.<br />
is shown as the angle from the x-axis to point c and appears as 2Cp in the<br />
Mohr's circle.<br />
axis and are at +(Cp-y).<br />
The angle Cp<br />
Points c and c' are symmetric about the principal stretching<br />
Geometrically, it is clear that when IDI I > DII<br />
the circle does not intersect the origin, and so there cannot be directions<br />
<strong>of</strong> no stretching. It is also clear that when lDIl = DII the two roots for<br />
4 coincide providing only one direction <strong>of</strong> no stretching. These special<br />
cases occur when the equations are elliptic and parabolic, respectively.<br />
The General Non-Normal Flow Rule Case<br />
When advection is included in the momentum equation, the eigenvectors<br />
are somewhat different from the simplest uncoupled system. Our attention<br />
returns to equations (59) and (60) to define the eigenvectors R a = 1, 2,<br />
-a<br />
We note first that the eigenvectors associated with the veZocity<br />
..., 4.<br />
characteristics (a = 3,4) are unchanged by including advection. These are<br />
defined by det C = 0 and for distinct roots (non-normal flow rule) require<br />
-a<br />
that det c 0. As before, we obtain R = 9. Then S satisfies (59) or<br />
-a -a -a<br />
(67) and is given by (78). The equation governing the solution along these<br />
characteristics remains the same (93) so that the velocity components still<br />
satisfy the simple relationships.<br />
as given in (94).<br />
The roots are also unchanged and remain<br />
A similar consideration <strong>of</strong> the equations governing solutions along<br />
stress characteristics shows changes. The stress characteristics (a = 1,2)<br />
are defined by det ?? = 0 and det # 0. Thus, from equation (60) we find<br />
-a<br />
the eigenvectors R to be unchanged from the uncoupled case and given by<br />
-a<br />
equations (73) and (74). But by (59) we find that S is no longer zero.<br />
-a<br />
Instead, we have<br />
since C is nonsingular. Substituting into equation (66) and evaluating<br />
--a<br />
all coefficients provide the equations governing solutions along these<br />
characteristics. Tt is seen that the forcing function and coefficient <strong>of</strong><br />
145