18.11.2014 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!