<TeXmacs|1.0.6>

<style|article>

<\body>
  <doc-data|<doc-title|EN221 Summary>>

  <section|Tensor Stuff>

  <\itemize>
    <item><em|Divergence>:

    <\equation*>
      <tabular|<tformat|<table|<row|<cell|\<nabla\>\<cdot\>\<b-u\>=\<partial\><rsub|i>u<rsub|i>>|<cell|<big|int><rsub|R>\<nabla\>\<cdot\>\<b-u\>=<big|int><rsub|\<partial\>R>\<b-u\><rsup|T>\<b-n\>\<mathd\>a,>>|<row|<cell|\<nabla\>\<cdot\>T=\<partial\><rsub|i>T<rsub|i
      j>\<b-e\><rsub|j>>|<cell|<big|int><rsub|R>\<nabla\>\<cdot\>T=<big|int><rsub|\<partial\>R>T<rsup|T>\<b-n\>\<mathd\>a.>>|<row|<cell|\<nabla\>\<otimes\>\<b-u\>=<with|mode|text|Jacobian>>|<cell|<big|int><rsub|R>\<nabla\>\<otimes\>\<b-u\>=<big|int><rsub|\<partial\>R>\<b-u\>\<otimes\>\<b-n\>\<mathd\>a,>>>>>
    </equation*>

    (matrix divergence: columns stay separate)

    <item><em|Box product>: <with|mode|math|[\<b-a\>,\<b-b\>,\<b-c\>]=\<b-a\>\<cdot\>(\<b-b\>\<wedge\>\<b-c\>)>

    <item><em|Levi-Civita tensor>:

    <\equation*>
      \<varepsilon\><rsub|i j k>=det<matrix|<tformat|<table|<row|<cell|\<delta\><rsub|i,1>>|<cell|\<delta\><rsub|j,1>>|<cell|\<delta\><rsub|k,1>>>|<row|<cell|\<delta\><rsub|i,2>>|<cell|\<delta\><rsub|j,2>>|<cell|\<delta\><rsub|k,2>>>|<row|<cell|\<delta\><rsub|i,3>>|<cell|\<delta\><rsub|j,3>>|<cell|\<delta\><rsub|k,3>>>>>>=[\<b-e\><rsub|i>,\<b-e\><rsub|j>,\<b-e\><rsub|k>]=<choice|<tformat|<table|<row|<cell|1>|<cell|(i
      j k)<with|mode|text| an even permut. of
      >(123),>>|<row|<cell|-1>|<cell|(i j k)<with|mode|text| an odd permut.
      of >(123),>>|<row|<cell|0>|<cell|<with|mode|text|if not>.>>>>>
    </equation*>

    <with|mode|math|\<b-e\><rsub|j>\<wedge\>\<b-e\><rsub|k>=\<varepsilon\><rsub|i
    j k>\<b-e\><rsub|i>>.

    <with|mode|math|det(\<b-a\>\<b-b\>\<b-c\>)=\<varepsilon\><rsub|i j
    k>a<rsub|i>b<rsub|j>c<rsub|k>>

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<varepsilon\><rsub|i j
      k>\<varepsilon\><rsub|i l m>>|<cell|=>|<cell|\<delta\><rsub|j
      l>\<delta\><rsub|k m>-\<delta\><rsub|j m>\<delta\><rsub|k
      l>,>>|<row|<cell|\<varepsilon\><rsub|i j k>\<varepsilon\><rsub|i j
      l>>|<cell|=>|<cell|2\<delta\><rsub|k
      l>,>>|<row|<cell|\<varepsilon\><rsub|i j k>\<varepsilon\><rsub|i j
      k>>|<cell|=>|<cell|6>>>>
    </eqnarray*>

    <item><em|Principal Invariants>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<Iota\><rsub|A>>|<cell|=>|<cell|\<lambda\><rsub|1>+\<lambda\><rsub|2>+\<lambda\><rsub|3>=([A\<b-a\>,\<b-b\>,\<b-c\>]+[\<b-a\>,A\<b-b\>,\<b-c\>]+[\<b-a\>,\<b-b\>,A\<b-c\>])/[\<b-a\>,\<b-b\>,\<b-c\>]=tr
      A,>>|<row|<cell|II<rsub|A>>|<cell|=>|<cell|\<lambda\><rsub|1>\<lambda\><rsub|2>+\<lambda\><rsub|2>\<lambda\><rsub|3>+\<lambda\><rsub|1>\<lambda\><rsub|3>=([A\<b-a\>,A\<b-b\>,\<b-c\>]+[A\<b-a\>,\<b-b\>,A\<b-c\>]+[\<b-a\>,A\<b-b\>,A\<b-c\>])/[\<b-a\>,\<b-b\>,\<b-c\>]=<frac|1|2>[tr<rsup|2>A-tr
      A<rsup|2>],>>|<row|<cell|III<rsub|A>>|<cell|=>|<cell|\<lambda\><rsub|1>\<lambda\><rsub|2>\<lambda\><rsub|3>=[A\<b-a\>,A\<b-b\>,A\<b-c\>]/[\<b-a\>,\<b-b\>,\<b-c\>]=det
      A.>>>>
    </eqnarray*>

    <item><em|Adjugate/Cofactor of a Tensor>:
    <with|mode|math|A<rsup|\<ast\>>(\<b-a\>\<wedge\>\<b-b\>)=(A\<b-a\>)\<wedge\>(A\<b-b\>)>
    <with|mode|math|\<Rightarrow\>> <with|mode|math|A<rsup|\<ast\>>=det
    A(A<rsup|-T>)>.

    <with|mode|math|\<partial\><rsub|t> det A(t)=det
    A*tr((\<partial\><rsub|t>A)A<rsup|-1>)>

    <item><em|Tensor Product>: <with|mode|math|<with|mode|text|TO>\<otimes\><with|mode|text|FROM>>

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<b-e\><rsub|i>\<otimes\>\<b-e\><rsub|j>>|<cell|=>|<cell|\<b-e\><rsub|i>\<b-e\><rsub|j><rsup|T>>>|<row|<cell|(\<b-u\>\<otimes\>\<b-v\>)\<b-a\>>|<cell|=>|<cell|\<b-u\>(\<b-v\>\<cdot\>\<b-a\>)>>|<row|<cell|(\<b-u\>\<otimes\>\<b-v\>)(\<b-w\>\<otimes\>\<b-x\>)>|<cell|=>|<cell|\<b-v\>\<cdot\>\<b-w\>(\<b-u\>\<otimes\>\<b-x\>)>>|<row|<cell|(\<b-u\>\<otimes\>\<b-v\>)A>|<cell|=>|<cell|\<b-u\>\<otimes\>(A<rsup|T>\<b-v\>)>>>>
    </eqnarray*>

    <item><em|Skewsymmetric matrices>: Rotation around axis
    <with|mode|math|<with|math-font-series|bold|\<Omega\>>> given by
    orthogonal matrix <with|mode|math|Q(t)>.

    <with|mode|math|<wide|\<b-x\>|\<dot\>>=<wide|Q|\<dot\>>\<b-x\>><with|mode|math|\<Rightarrow\>\<partial\><rsub|t>(Q<rsup|T>Q)=0>,\ 

    <with|mode|math|W=<wide|Q|\<dot\>>Q<rsup|T>>,
    <with|mode|math|W=-W<rsup|T>>. <with|mode|math|W\<b-x\>=<with|math-font-series|bold|\<Omega\>>\<wedge\>\<b-x\>>.
  </itemize>

  <section|Kinematics>

  <subsection|Static>

  <\itemize>
    <item><em|Reference> and <em|deformed configurations>.

    <item><em|Deformation gradient>: assumed regular. <with|mode|math|J=det
    F\<neq\>0>.

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<b-x\>(\<b-X\>)>|<cell|=>|<cell|\<b-X\>+\<b-u\>(\<b-X\>),>>|<row|<cell|F>|<cell|=>|<cell|\<nabla\><rsub|\<b-X\>>\<otimes\>\<b-x\>(\<b-X\>),>>|<row|<cell|F<rsup|-1>>|<cell|=>|<cell|\<partial\><rsub|x<rsub|j>>X<rsub|\<beta\>>\<b-E\><rsub|\<beta\>>\<otimes\>\<b-e\><rsub|i>>>>>
    </eqnarray*>

    <item><em|Isochoric>: <with|mode|math|J=1>.

    <item><em|Polar decomposition>:

    <\itemize>
      <item><with|mode|math|F=R*U>,

      <with|mode|math|F<rsup|T>F=U<rsup|2>>, <with|mode|math|R=F*U<rsup|-1>>.

      <item><with|mode|math|F=V*R>.
    </itemize>

    Features:

    <\itemize>
      <item>Is unique.

      <item><with|mode|math|R> is rotation of principal axes.

      <item><with|mode|math|R> average of all rotations.

      <item>Principal axes of <with|mode|math|V> are
      <with|mode|math|R\<b-u\><rsub|i>>.\ 

      <item><with|mode|math|\<sigma\>(V)=\<sigma\>(U)>.

      <item><with|mode|math|R=\<b-v\><rsub|k>\<otimes\>\<b-u\><rsub|k>>.

      <item><with|mode|math|F=\<lambda\><rsub|k>\<b-v\><rsub|k>\<otimes\>\<b-u\><rsub|k>>.
    </itemize>

    <item><em|Left/Right Cauchy-Green Deformation Tensor>:
    <with|mode|math|F*F<rsup|T>/F<rsup|T>F> SPD.

    <item><em|Strain>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|E>|<cell|=>|<cell|<frac|1|2>(F<rsup|T>F-Id)<space|1em>(<with|mode|text|Lagrangean>:<space|1em>\|\<mathd\>\<b-x\>\|<rsup|2>-\|\<mathd\>\<b-X\>\|<rsup|2>=2\<mathd\>\<b-X\>\<cdot\>E\<mathd\>\<b-X\>),>>|<row|<cell|E<rprime|'>>|<cell|=>|<cell|<frac|1|2>(Id-F<rsup|-T>F<rsup|-1>)<space|1em>(Eulerian:<space|1em>\|\<mathd\>\<b-x\>\|<rsup|2>-\|\<mathd\>\<b-X\>\|<rsup|2>=2\<mathd\>\<b-x\>\<cdot\>E<rprime|'>\<mathd\>\<b-x\>).>>>>
    </eqnarray*>

    <item><em|Stretch>:

    <\equation*>
      \<lambda\>(\<b-M\>)=(\<b-M\>\<cdot\>F<rsup|T>F\<b-M\>)<rsup|1/2>=\|U\<b-M\>\|.
    </equation*>

    Has local maxima and minima when <with|mode|math|\<b-M\>> is an
    eigenvector of <with|mode|math|U>.

    <item><em|Transformation of area elements>:

    <\equation*>
      \<b-n\>\<mathd\>a=F<rsup|\<ast\>>\<b-N\>\<mathd\>A
    </equation*>

    <item><em|Deformation gradient in cylindrical coordinates>: Given

    <\equation*>
      <matrix|<tformat|<table|<row|<cell|r>>|<row|<cell|\<theta\>>>|<row|<cell|z>>>>>=f(R,\<Theta\>,Z),
    </equation*>

    we have

    <\equation*>
      F=\<partial\><rsub|R>\<b-x\>\<otimes\>\<b-E\><rsub|R>+<frac|1|R>\<partial\><rsub|\<Theta\>>\<b-x\>\<otimes\>\<b-E\><rsub|\<Theta\>>+\<partial\><rsub|Z>\<b-x\>\<otimes\>\<b-E\><rsub|Z>.
    </equation*>

    Also expressible as <em|mixed tensor> from
    <with|mode|math|\<b-E\><rsub|(R,\<Theta\>,Z)>> to
    <with|mode|math|\<b-E\><rsub|(r,\<theta\>,z)>>:

    <\equation*>
      F=<matrix|<tformat|<table|<row|<cell|<phantom|<frac|r|R>*><frac|\<partial\>r|\<partial\>R>>|<cell|<frac|1|<with|color|blue|R>>*<frac|\<partial\>r|\<partial\>\<Theta\>>>|<cell|<phantom|<frac|r|R>*><frac|\<partial\>r|\<partial\>z>>>|<row|<cell|<frac|<with|color|orange|r>|1><frac|\<partial\>\<theta\>|\<partial\>R>>|<cell|<frac|<with|color|orange|r>|<with|color|blue|R>>*<frac|\<partial\>\<theta\>|\<partial\>\<Theta\>>>|<cell|<frac|<with|color|orange|r>|1><frac|\<partial\>\<theta\>|\<partial\>z>>>|<row|<cell|<phantom|<frac|r|R>*><frac|\<partial\>z|\<partial\>R>>|<cell|<frac|1|<with|color|blue|R>><frac|\<partial\>z|\<partial\>\<Theta\>>>|<cell|<phantom|<frac|r|R>*><frac|\<partial\>z|\<partial\>Z>>>>>>.
    </equation*>

    Caveat for mixed tensors: <with|mode|math|tr(F)\<neq\>F<rsub|i i>>.
    However <with|mode|math|det>, <with|mode|math|V>, <with|mode|math|U> as
    usual. Also works for spherical basis, but more complicated.
  </itemize>

  <subsubsection|Static Examples>

  <\itemize>
    <item><em|Pure shear>: <with|mode|math|F=\<lambda\>\<b-e\><rsub|1>\<otimes\>\<b-e\><rsub|1>+\<lambda\><rsup|-1>\<b-e\><rsub|2>\<otimes\>\<b-e\><rsub|2>>.

    <item><em|Simple shear>: <with|mode|math|F=Id+\<lambda\>\<b-e\><rsub|1>\<otimes\>\<b-e\><rsub|2>>.

    <item><em|Pure bending>:

    <\equation*>
      <matrix|<tformat|<table|<row|<cell|x>>|<row|<cell|y>>|<row|<cell|z>>>>>=<matrix|<tformat|<table|<row|<cell|(R-Y)sin
      \<alpha\>(x)>>|<row|<cell|R-(R-Y)cos
      \<alpha\>(x)>>|<row|<cell|Z>>>>>,<space|1em>J=(R-Y)\<alpha\><rprime|'>.
    </equation*>

    <item><em|Tension and torsion>:

    <\equation*>
      <matrix|<tformat|<table|<row|<cell|x>>|<row|<cell|y>>|<row|<cell|z>>>>>=<matrix|<tformat|<table|<row|<cell|<frac|X|<sqrt|\<lambda\>>>cos<left|(><frac|\<alpha\>|l>\<lambda\>Z<right|)>-<frac|Y|<sqrt|\<lambda\>>>sin<left|(><frac|\<alpha\>|l>\<lambda\>Z<right|)>>>|<row|<cell|<frac|X|<sqrt|\<lambda\>>>sin<left|(><frac|\<alpha\>|l>\<lambda\>Z<right|)>+<frac|Y|<sqrt|\<lambda\>>>cos<left|(><frac|\<alpha\>|l>\<lambda\>Z<right|)>>>|<row|<cell|\<lambda\>Z>>>>>.
    </equation*>

    <item><em|Turning a cylinder inside out>.
  </itemize>

  <subsection|Dynamic>

  <\itemize>
    <item><em|Steady motion>: <with|mode|math|\<partial\>/\<partial\>t
    \<b-v\>(\<b-x\>,t)=0>.

    <item><em|Material/Lagrangean POV>: focus on particle, expressions in
    terms of <with|mode|math|\<b-X\>> and
    <with|mode|math|t><with|mode|math|\<rightarrow\>>Solids.

    <item><em|Spatial/Eulerian POV>: focus on point in space, expressions in
    terms of <with|mode|math|\<b-x\>> and
    <with|mode|math|t><with|mode|math|\<rightarrow\>>Fluids.

    <item>Lines:

    <\itemize>
      <item><em|Path line>: Curve traced by a fixed particle.

      <item><em|Streamlines>: Field lines of velocity in Eulerian POV.
    </itemize>

    Both coincide under steady motion.

    <item><em|Material derivative>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|<wide|\<varphi\>|\<dot\>>>|<cell|=>|<cell|<frac|\<partial\>\<varphi\>|\<partial\>t>+\<nabla\><rsub|\<b-x\>>\<varphi\>\<cdot\>\<b-v\>,>>|<row|<cell|<wide|\<b-w\>|\<dot\>>>|<cell|=>|<cell|<frac|\<partial\>\<b-w\>|\<partial\>t>+(\<nabla\><rsub|\<b-x\>>\<otimes\>\<b-w\>)\<b-v\>,>>|<row|<cell|<wide|T|\<dot\>>>|<cell|=>|<cell|<frac|\<partial\>T|\<partial\>T>+(\<nabla\><rsub|\<b-x\>>\<otimes\>T)\<b-v\>.>>>>
    </eqnarray*>

    <item><em|Acceleration>: <with|mode|math|\<b-a\>=<wide|\<b-v\>|\<dot\>>>.

    <item><em|Velocity gradient>: <with|mode|math|L=\<nabla\><rsub|\<b-x\>>\<otimes\>\<b-v\><space|1em>\<Rightarrow\><space|1em><wide|F|\<dot\>>=L*F>
    (chain rule).

    <with|mode|math|F> requires a ``reference state'', <with|mode|math|L>
    does not.

    <item> <with|mode|math|\<mathd\><wide|\<b-x\>|\<dot\>>=<wide|F|\<dot\>>\<mathd\>\<b-X\>=L*F\<mathd\>\<b-X\>=L\<mathd\>\<b-x\>>.
    Assume <with|mode|math|\<mathd\>\<b-x\>=\<b-m\>\|\<mathd\>\<b-x\>\|>.

    <\eqnarray*>
      <tformat|<table|<row|<cell|<with|mode|text|<em|Strain
      rate>:><space|1em><frac|\|\<mathd\>\<b-x\>\|<rsup|\<bullet\>>|\|\<mathd\>\<b-x\>\|>>|<cell|=>|<cell|\<b-m\>\<cdot\>L\<b-m\>=\<b-m\>\<cdot\>D\<b-m\>>>|<row|<cell|<wide|\<b-m\>|\<dot\>>>|<cell|=>|<cell|L\<b-m\>-\<b-m\>(\<b-m\>\<cdot\>L\<b-m\>)>>>>
    </eqnarray*>

    <item><em|Stretch and Spin>: <with|mode|math|L=D+W>,
    <with|mode|math|D=D<rsup|T>>, <with|mode|math|W=-W<rsup|T>>.

    <with|mode|math|D<rsub|11>>: stretching rate of a line element along the
    1-direction

    <with|mode|math|D<rsub|12>>: (roughly) change in angle between the 1- and
    2-direction.

    Principal axes <with|mode|math|\<b-p\><rsub|i>> of <with|mode|math|D> are
    rigidly rotating about

    <\equation*>
      <with|math-font-series|bold|\<omega\>>=<frac|1|2>curl \<b-v\>
    </equation*>

    with <with|mode|math|W\<b-p\><rsub|i>=<with|math-font-series|bold|\<omega\>>\<times\>\<b-p\><rsub|i>>.

    <item><em|Vorticity>: <with|mode|math|curl
    \<b-v\>=2\<cdot\><with|mode|text|angular velocity>>. (Letter here is also
    <with|mode|math|<with|math-font-series|bold|\<omega\>>>.)

    <item><with|mode|math|<wide|J|\<dot\>>=J*tr L=J div \<b-v\>>.

    <item><em|Integrals over moving contours>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|<big|oint><rsub|C<rsub|t>>\<b-v\>\<cdot\>\<mathd\>\<b-x\>>|<cell|=>|<cell|<big|oint><rsub|C<rsub|R>>\<b-v\>(\<b-x\>,t)\<cdot\>F\<mathd\>\<b-X\>>>|<row|<cell|<frac|\<mathd\>|\<mathd\>t><big|oint><rsub|C<rsub|t>>\<b-v\>\<cdot\>\<mathd\>\<b-x\>>|<cell|=>|<cell|<big|oint><rsub|C<rsub|R>><left|(><wide|\<b-v\>|\<dot\>>(\<b-x\>,t)\<cdot\>F+\<b-v\>(\<b-x\>,t)\<cdot\>L*F<right|)>\<mathd\>\<b-X\>>>|<row|<cell|>|<cell|=>|<cell|<big|oint><rsub|C<rsub|t>><wide|\<b-v\>|\<dot\>>(\<b-x\>,t)+L<rsup|T>\<b-v\>(\<b-x\>,t)\<mathd\>\<b-x\>>>>>
    </eqnarray*>

    <item><em|Integrals over moving surfaces>: Similar, taking into account
    that <with|mode|math|F<rsup|\<ast\>>=J*F<rsup|-T>>.

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int><rsub|S<rsub|t>>\<b-u\>\<cdot\>\<b-n\>\<mathd\>s=<big|int><rsub|S<rsub|t>>(<wide|\<b-u\>|\<dot\>>+\<b-u\>*tr(L)-L\<b-u\>)\<cdot\>\<b-n\>\<mathd\>s.
    </equation*>

    <item><em|Integrals over moving volumes/Reynolds' Transport Theorem>:

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int><rsub|R<rsub|t>>\<varphi\>(\<b-x\>)\<mathd\>v=<big|int><rsub|R<rsub|t>><wide|\<varphi\>|\<dot\>>+\<varphi\>*tr(L)\<mathd\>v.
    </equation*>

    Observe that <with|mode|math|tr(L)=div \<b-v\>>, which is zero in the
    incompressible case.

    <item><em|Circulation>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|<big|oint><rsub|C<rsub|t>>\<b-v\>\<cdot\>\<mathd\>\<b-x\>>|<cell|=>|<cell|<big|int><rsub|S<rsub|t>>curl
      \<b-v\>\<cdot\>\<mathd\>\<b-s\>=<big|int><rsub|S<rsub|t>><with|math-font-series|bold|\<omega\>>\<cdot\>\<mathd\>\<b-s\>>>|<row|<cell|L<rsup|T>\<b-v\>>|<cell|=>|<cell|<frac|1|2>\<nabla\>v<rsup|2>>>|<row|<cell|0<space|1em><above|=|<with|mode|text|if
      circulation-preserving>><space|1em><frac|\<mathd\>|\<mathd\>t><big|oint><rsub|C<rsub|t>>\<b-v\>\<cdot\>\<mathd\>\<b-x\>>|<cell|=>|<cell|<big|oint><rsub|C<rsub|t>><wide|\<b-v\>|\<dot\>>(\<b-x\>,t)+L<rsup|T>\<b-v\>(\<b-x\>,t)\<mathd\>\<b-x\>>>|<row|<cell|>|<cell|=>|<cell|<big|oint><rsub|C<rsub|t>><wide|\<b-v\>|\<dot\>>(\<b-x\>,t)\<mathd\>\<b-x\>+<big|oint><rsub|C<rsub|t>><frac|1|2>\<nabla\>v<rsup|2>\<mathd\>\<b-x\>>>|<row|<cell|>|<cell|=>|<cell|<big|oint><rsub|C<rsub|t>>\<b-a\>(\<b-x\>,t)\<mathd\>\<b-x\>>>|<row|<cell|>|<cell|=>|<cell|<big|int><rsub|S<rsub|t>>curl
      \<b-a\>\<cdot\>\<mathd\>\<b-s\>>>>>
    </eqnarray*>

    If <with|mode|math|\<b-a\>=\<nabla\>\<psi\>>, then the motion is
    <em|circulation-preserving>.

    If circulation-preserving, then

    <\equation*>
      curl \<b-a\>=<wide|<with|math-font-series|bold|\<omega\>>|\<dot\>>+<with|math-font-series|bold|\<omega\>>tr(L)-L<with|math-font-series|bold|\<omega\>>=0.
    </equation*>

    Then consider the product rule on

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t>(J*F<rsup|-1><with|math-font-series|bold|\<omega\>>)=\<ldots\>=0
    </equation*>

    to find <em|Cauchy's vorticity formula>:

    <\equation*>
      <with|math-font-series|bold|\<omega\>>=<frac|1|J>F<with|math-font-series|bold|\<omega\>><rsub|ref>.
    </equation*>

    Field lines of vorticity are <em|vortex lines>.

    If the motion is circulation-preserving, these are material curves.
  </itemize>

  <section|Balance Laws and Field Equations>

  <\itemize>
    <item><em|Conservation of Mass>: Assumption:

    <\eqnarray*>
      <tformat|<table|<row|<cell|J\<rho\>>|<cell|=>|<cell|\<rho\><rsub|ref><space|1em>(referential).>>>>
    </eqnarray*>

    Therefore,

    <\eqnarray*>
      <tformat|<table|<row|<cell|<wide|\<rho\>|\<dot\>>J+\<rho\><wide|J|\<dot\>>>|<cell|=>|<cell|0>>|<row|<cell|<wide|\<rho\>|\<dot\>>J+\<rho\>*J*div
      \<b-v\>>|<cell|=>|<cell|0>>|<row|<cell|<wide|\<rho\>|\<dot\>>+\<rho\>*div
      \<b-v\>>|<cell|=>|<cell|0>>|<row|<cell|\<partial\><rsub|t>\<rho\>+div(\<rho\>\<b-v\>)>|<cell|=>|<cell|0.>>>>
    </eqnarray*>

    <item><em|Transport Theorem with density>:

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int><rsub|R<rsub|t>>\<rho\>\<Phi\>\<mathd\>v=<big|int><rsub|R<rsub|t>>\<rho\><wide|\<Phi\>|\<dot\>>\<mathd\>v.
    </equation*>

    <item><em|Linear Momentum>: <with|mode|math|M=\<rho\>\<b-v\>>

    <\itemize>
      <item><em|Stress vector>: <with|mode|math|\<b-t\><rsub|(\<b-n\>)>> is
      force/unit area.

      <item><em|Balance law>:

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t><big|int><rsub|R<rsub|t>>\<rho\>\<b-v\>=<big|int><rsub|R<rsub|t>>\<rho\><wide|\<b-v\>|\<dot\>>=<big|int><rsub|R<rsub|t>>\<rho\>\<b-b\>\<mathd\>v+<big|int><rsub|\<partial\>R<rsub|t>>\<b-t\><rsub|(\<b-n\>)>\<mathd\>a.
      </equation*>

      <item><em|Stress tensor/Cauchy's Theorem>:
      <with|mode|math|\<sigma\><rsup|T>\<b-n\>=\<b-t\><rsub|(\<b-n\>)>>.
      Derivation:

      <\itemize>
        <item><with|mode|math|\<b-t\><rsub|(-\<b-n\>)>=-\<b-t\><rsub|(\<b-n\>)>>
        by pillbox and balance law.

        <item>Tetrahedron argument: <with|mode|math|\<b-n\>> the general
        normal of the coordinate-system-boxed tetrahedron. Then other faces
        <with|mode|math|a<rsub|i>=a\<b-n\><rsub|i>>, where <with|mode|math|a>
        is the area of the complicated face. Volume <with|mode|math|h*a/3>.
        Apply <with|mode|math|1/a\<cdot\>>balance law, let
        <with|mode|math|h\<rightarrow\>0>. Assume continuity, derive linear
        dependence by assuming values are locally constant.
      </itemize>

      <item>Updated balance law:

      <\equation*>
        <big|int><rsub|R<rsub|t>>\<rho\>\<b-a\>=<big|int><rsub|R<rsub|t>>\<rho\>\<b-b\>+\<nabla\>\<cdot\>\<sigma\>
      </equation*>

      <item><em|Field equations>:

      <\eqnarray*>
        <tformat|<table|<row|<cell|\<rho\><rsub|ref><wide|\<b-x\>|\<ddot\>>>|<cell|=>|<cell|\<nabla\><rsub|X>\<cdot\>s+\<rho\><rsub|ref>\<b-b\><space|1em>(<with|mode|text|referential>)>>|<row|<cell|\<rho\>\<b-a\>>|<cell|=>|<cell|\<nabla\><rsub|\<b-x\>>\<cdot\>\<sigma\>+\<rho\>\<b-b\><space|1em>(<with|mode|text|spatial>)>>>>
      </eqnarray*>
    </itemize>

    <item><em|Nominal Stress/Conjugate stress>:
    <with|mode|math|s=J*F<rsup|-1>\<sigma\>>.
    (<with|mode|math|\<sigma\><rsup|T>\<b-n\>\<mathd\>\<b-a\>=s<rsup|T>\<b-N\>\<mathd\>\<b-A\>>
    can be directly verified.)

    Also called <em|Piola-Kirchoff stress>. <with|mode|math|s<rsup|T>> is
    <em|2nd Piola-Kirchoff stress>.

    <item><em|Angular Momentum>: <with|mode|math|H=\<rho\>\<b-x\>\<wedge\>\<b-v\>>

    <\itemize>
      <item><em|Non-polar material>: no contact torques.

      <item><em|Balance law>:

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t><big|int><rsub|R<rsub|t>>\<rho\>\<b-x\>\<wedge\>\<b-v\><above|=|(\<ast\>)><big|int><rsub|R<rsub|t>>\<rho\>\<b-x\>\<wedge\><wide|\<b-v\>|\<dot\>>=<big|int><rsub|R<rsub|t>>\<rho\>(\<b-x\>\<wedge\>\<b-b\>+\<b-c\>)\<mathd\>v+<big|int><rsub|\<partial\>R<rsub|t>>\<b-x\>\<wedge\>\<b-t\><rsub|(\<b-n\>)>\<mathd\>a
      </equation*>

      <with|mode|math|\<b-c\>> is body torque. Equality
      <with|mode|math|(\<ast\>)> follows because
      <with|mode|math|\<b-x\>>-derivatives vanish once
      <with|mode|math|\<wedge\>\<b-v\>> is applied.

      <item>Subsituting Cauchy's Theorem into the balance law gives

      <\eqnarray*>
        <tformat|<table|<row|<cell|<big|int><rsub|R<rsub|t>>\<b-x\>\<wedge\><left|(>\<nabla\><rsub|\<b-x\>>\<cdot\>\<sigma\>+\<rho\>\<b-b\><right|)>>|<cell|=>|<cell|<big|int><rsub|R<rsub|t>>\<rho\>(\<b-x\>\<wedge\>\<b-b\>)\<mathd\>v+<big|int><rsub|\<partial\>R<rsub|t>>\<b-x\>\<wedge\>\<sigma\>\<b-n\>\<mathd\>a>>|<row|<cell|<big|int><rsub|R<rsub|t>>\<b-x\>\<wedge\><left|(>\<nabla\><rsub|\<b-x\>>\<cdot\>\<sigma\><right|)>>|<cell|>|<cell|<big|int><rsub|\<partial\>R<rsub|t>>\<b-x\>\<wedge\>\<sigma\>\<b-n\>\<mathd\>a>>>>
      </eqnarray*>

      View in component form, apply Gau˙, derive
      <with|mode|math|\<varepsilon\><rsub|i j k>\<sigma\><rsub|j
      i>=0><with|mode|math|\<Rightarrow\>><with|mode|math|\<sigma\>=\<sigma\><rsup|T>>.

      <item><em|Field equations>:

      <\eqnarray*>
        <tformat|<table|<row|<cell|s<rsup|T>F<rsup|T>>|<cell|=>|<cell|F*s<space|1em>(<with|mode|text|referential>)>>|<row|<cell|\<sigma\><rsup|T>>|<cell|=>|<cell|\<sigma\><space|1em>(<with|mode|text|spatial>)>>>>
      </eqnarray*>
    </itemize>

    <item><em|Vector identities>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|(\<b-A\>\<cdot\>\<nabla\>)\<b-A\>>|<cell|=>|<cell|<frac|1|2>\<nabla\>\|\<b-A\>\|<rsup|2>+(\<nabla\>\<times\>\<b-A\>)\<times\>\<b-A\>>>|<row|<cell|(\<b-A\>\<cdot\>\<nabla\>)\<b-A\>>|<cell|=>|<cell|(\<nabla\>\<otimes\>\<b-A\>)\<b-A\>.>>>>
    </eqnarray*>

    Use these identities to rewrite the <with|mode|math|<wide|\<b-v\>|\<dot\>>>
    as <with|mode|math|\<nabla\>(\<b-v\><rsup|2>)> for irrotational flow.

    <item><em|Types of fluid flow>:

    <\itemize>
      <item><em|Inviscid>: <with|mode|math|\<sigma\>=-p*Id><with|mode|math|\<Rightarrow\>><with|mode|math|div
      \<sigma\>=-\<nabla\>p>.

      <item><em|Incompressible>: <with|mode|math|<wide|\<rho\>|\<dot\>>=0> or
      <with|mode|math|div \<b-v\>=0>.

      <item><em|Steady>: <with|mode|math|\<partial\><rsub|t>\<b-v\>=\<b-0\>>.

      <with|mode|math|<wide|\<rho\>|\<dot\>>=\<b-v\>\<cdot\>\<nabla\>\<rho\>>.

      <item><em|Irrotational>: <with|mode|math|<with|math-font-series|bold|\<omega\>>=\<b-0\>>
      or <with|mode|math|\<b-v\>=\<nabla\>\<varphi\>>.

      <item><em|Elastic>: <with|mode|math|p(\<rho\>)>

      <item><em|Ideal=incompressible>: <with|mode|math|div \<b-v\>=0>,
      <with|mode|math|J=1>.
    </itemize>

    <item><em|Rayleigh-Plesset equation>: Begin with deformation of spherical
    shell (with extent!), assume <with|mode|math|J\<equiv\>1>. Derive ODE.

    <item><em|Conservative potentials>: <with|mode|math|\<b-b\>=-\<nabla\>\<beta\>>

    <\itemize>
      <item>Elastic or ideal flow here is circulation preserving, i.e.
      <with|mode|math|\<b-a\>=-\<nabla\>something>.

      <\itemize>
        <item>Have

        <\eqnarray*>
          <tformat|<table|<row|<cell|\<b-a\>>|<cell|=>|<cell|-<frac|1|\<rho\>>\<nabla\>p(\<rho\>)+\<b-b\>>>|<row|<cell|>|<cell|=>|<cell|-<frac|1|\<rho\>>p<rprime|'>(\<rho\>)\<nabla\>\<rho\>-\<nabla\>\<beta\>.>>>>
        </eqnarray*>

        <item>Define

        <\eqnarray*>
          <tformat|<table|<row|<cell|\<varepsilon\>(\<rho\>)>|<cell|=>|<cell|<big|int><rsub|0><rsup|\<rho\>><frac|1|\<rho\><rprime|'>>p<rprime|'>(\<rho\><rprime|'>)\<mathd\>\<rho\><rprime|'>>>|<row|<cell|\<Rightarrow\>\<nabla\>\<varepsilon\>(\<rho\>)>|<cell|=>|<cell|\<varepsilon\><rprime|'>(\<rho\>)\<nabla\>\<rho\>>>>>
        </eqnarray*>

        <item>Therefore

        <\equation*>
          \<b-a\>=-\<nabla\>(\<varepsilon\>(\<rho\>)+\<beta\>).
        </equation*>

        <item>For ideal fluid substitute <with|mode|math|p/\<rho\><rsub|0>>
        for <with|mode|math|\<varepsilon\>>.
      </itemize>

      <item><em|Bernoulli's Theorem>:\ 

      <\itemize>
        <item>Flow irrotational (i.e. <with|mode|math|\<b-v\>=-\<nabla\>\<varphi\>>):

        <\equation*>
          \<nabla\><left|(>\<partial\><rsub|t>\<varphi\>+<frac|v<rsup|2>|2>+\<varepsilon\>(\<rho\>)+\<beta\><right|)>=\<b-0\>.
        </equation*>

        Proof: Just rewrite, obtaining <with|mode|math|v<rsup|2>/2> from
        second term of material derivative.

        <item>Flow steady:\ 

        <\equation*>
          <left|(><frac|v<rsup|2>|2>+\<varepsilon\>(\<rho\>)+\<beta\><right|)><rsup|\<bullet\>>=0,
        </equation*>

        i.e. this quantity is constant along streamlines.

        Proof: Exploit <with|mode|math|\<b-v\>\<cdot\><wide|\<b-v\>|\<dot\>>=\<b-v\>\<cdot\>\<nabla\>(\<b-v\><rsup|2>)>

        <item>Flow both irrotational and steady:

        <\equation*>
          \<nabla\><left|(><frac|v<rsup|2>|2>+\<varepsilon\>(\<rho\>)+\<beta\><right|)>=\<b-0\>.
        </equation*>
      </itemize>
    </itemize>

    <item><em|Acoustic wave equation>:

    <\itemize>
      <item>Assume <with|mode|math|\<rho\>=\<rho\><rsub|0>+\<delta\>\<rho\>>,
      <with|mode|math|\|\<b-v\>\|\<ll\>1>,
      <with|mode|math|\|\<nabla\>\<b-v\>\|\<ll\>1>.

      <item>Start with <with|mode|math|\<partial\><rsub|t>*(\<nabla\>\<cdot\>\<b-v\>)>,
      use cons. of. momentum without second order term, cons. of mass as
      <with|mode|math|\<partial\><rsub|t>\<rho\>+\<rho\><rsub|0>div(\<b-v\>)=0.>

      <item><with|mode|math|\<delta\>\<rho\><rsub|t
      t>=c<rsup|2>\<nabla\><rsup|2>\<delta\>\<rho\>>, with
      <with|mode|math|c=<sqrt|\<partial\><rsub|\<rho\>>p>>.
    </itemize>

    <item><em|Mach number>: assume steadiness <with|mode|math|\<b-b\>=0>, use
    <with|mode|math|\<b-v\>\<cdot\><wide|\<b-v\>|\<dot\>>> in terms of
    <with|mode|math|c<rsup|2>>.

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<b-v\>\<cdot\>(\<rho\>\<b-v\>)<rsup|\<bullet\>>>|<cell|=>|<cell|\<b-v\>\<cdot\><wide|\<b-v\>|\<dot\>>\<rho\><left|(>1-<wide*|<frac|\<b-v\><rsup|2>|c<rsup|2>>|\<wide-underbrace\>><rsub|m\<assign\>><right|)>>>>>
    </eqnarray*>

    <\itemize>
      <item>Supersonic nozzle <with|mode|math|m\<less\>1>,
      <with|mode|math|m\<gtr\>1>.
    </itemize>

    <item><em|Conservation of Energy>:

    <\itemize>
      <item><em|Balance law>:

      <\eqnarray*>
        <tformat|<table|<row|<cell|<frac|\<mathd\>|\<mathd\>t>K(R<rsub|t>)>|<cell|=>|<cell|-S(R<rsub|t>)+P(R<rsub|t>)>>|<row|<cell|<frac|\<mathd\>|\<mathd\>t><wide*|<frac|1|2><big|int>\<rho\>\<b-v\>\<cdot\>\<b-v\>\<mathd\>v|\<wide-underbrace\>><rsub|<with|mode|text|Kinetic
        energy> K(t)>>|<cell|=>|<cell|-<wide*|<big|int><rsub|R<rsub|t>>tr(\<sigma\>D)\<mathd\>v|\<wide-underbrace\>><rsub|<with|mode|text|Stress
        power> S(R<rsub|t>)>+<wide*|<big|int><rsub|R<rsub|t>>\<rho\>\<b-b\>\<cdot\>\<b-v\>+<big|int><rsub|\<partial\>R<rsub|t>>\<sigma\>\<b-n\>\<cdot\>\<b-v\>\<mathd\>a|\<wide-underbrace\>><rsub|<with|mode|text|Power
        supplied> P(R<rsub|t>)>>>>>
      </eqnarray*>

      \;

      Proof: Multiply Equation of Motion by <with|mode|math|\<b-v\>>,
      integrate by parts in the <with|mode|math|\<sigma\>> term.

      <item><em|Field equation>:

      <\equation*>
        <wide*|\<rho\><left|(><frac|1|2>v\<cdot\>v<right|)><rsup|\<bullet\>>|\<wide-underbrace\>><rsub|<with|mode|text|Kinetic
        Energy>>+<wide*|tr(\<sigma\>D)|\<wide-underbrace\>><rsub|<with|mode|text|Stress
        Power>>=<wide*|\<nabla\><rsub|\<b-x\>>\<cdot\>(\<sigma\>\<b-v\>)+\<rho\>\<b-b\>\<cdot\>\<b-v\>|\<wide-underbrace\>><rsub|<with|mode|text|Rate-of-working>>.
      </equation*>

      <item>Key words for more global energy conservation: internal energy
      <with|mode|math|U(R<rsub|t>)>, heat supply per unit mass
      <with|mode|math|H(R<rsub|t>)>, heat flux through material surface.

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t>{K+U}=P(R<rsub|t>)+H(R<rsub|t>).
      </equation*>

      Now, because there is a stress power loss above, there needs to be a
      gain here:

      <\equation*>
        <frac|\<mathd\>|\<mathd\>t>U(R<rsub|t>)=S(R<rsub|t>)+H(R<rsub|t>).
      </equation*>
    </itemize>

    <item><em|Jump conditions>: For the balance law

    <\equation*>
      <frac|\<mathd\>|\<mathd\>t><big|int><rsub|R<rsub|t>>\<rho\>\<pi\>=<big|int><rsub|R<rsub|t>>\<rho\>s+<big|int><rsub|\<partial\>R<rsub|t>>f<rsub|(\<b-n\>)>,
    </equation*>

    we get

    <\equation*>
      [\<rho\>V\<pi\>+f<rsub|(\<b-n\>)>]=0.
    </equation*>

    <with|mode|math|V<rsub|n>> interface speed,
    <with|mode|math|V=V<rsub|n>-\<b-v\>\<cdot\>\<b-n\>>.

    <tabular|<tformat|<cwith|2|2|1|-1|cell-tborder|0.5pt>|<cwith|1|-1|3|3|cell-lborder|0.5pt>|<cwith|1|-1|2|2|cell-lborder|0.5pt>|<table|<row|<cell|>|<cell|>|<cell|Mass>|<cell|Mom.>|<cell|A.Mom.>|<cell|Energy>>|<row|<cell|<with|mode|math|\<pi\>>>|<cell|quantity
    per unit mass>|<cell|1>|<cell|<with|mode|math|\<b-v\>>>|<cell|<with|mode|math|\<b-x\>\<wedge\>\<b-v\>>>|<cell|<with|mode|math|\<varepsilon\>+<frac|1|2>\<b-v\>\<cdot\>\<b-v\>>>>|<row|<cell|s>|<cell|supply
    of <with|mode|math|\<pi\>> per unit mass
    >|<cell|0>|<cell|<with|mode|math|\<b-b\>>>|<cell|<with|mode|math|\<b-x\>\<wedge\>\<b-b\>>>|<cell|<with|mode|math|\<b-b\>\<cdot\>\<b-v\>+r>>>|<row|<cell|<with|mode|math|f<rsub|(\<b-n\>)>>>|<cell|influx
    of <with|mode|math|\<pi\>> per unit area>|<cell|0>|<cell|<with|mode|math|\<b-t\><rsub|(\<b-n\>)>>>|<cell|<with|mode|math|\<b-x\>\<wedge\>\<b-t\><rsub|(\<b-n\>)>>>|<cell|<with|mode|math|\<b-t\><rsub|(\<b-n\>)>\<cdot\>\<b-v\>+h<rsub|(\<b-n\>)>>>>>>>

    so that for example

    <\eqnarray*>
      <tformat|<table|<row|<cell|[\<rho\>V]>|<cell|=>|<cell|0,>>|<row|<cell|[\<rho\>V\<b-v\>+\<b-t\><rsub|(\<b-n\>)>]>|<cell|=>|<cell|\<b-0\>.>>>>
    </eqnarray*>

    Or for material jumps: <with|mode|math|[\<b-t\><rsub|(\<b-n\>)>]=\<b-0\>>.

    Derivation:

    <\itemize>
      <item>Modification for moving boundary is

      <\equation*>
        -<big|int><rsub|jump surface>[\<rho\>\<pi\>]V<rsub|n>.
      </equation*>

      <item>Then use pillbox that flattens around surface.
    </itemize>

    Examples:

    <\itemize>
      <item>Free boundary: pressure must be continuous, because otherwise
      there's a finite force on something massless.
    </itemize>

    <item><em|Stokes waves>:

    <\itemize>
      <item>Assume <with|mode|math|\<b-v\>=\<nabla\>\<varphi\>>.

      <item>Conservation of mass <with|mode|math|\<nabla\><rsup|2>\<varphi\>=0>.

      <item>Bernoulli's equation

      <\equation*>
        \<partial\><rsub|t>\<varphi\>+<frac|v<rsup|2>|2>+<frac|p|\<rho\><rsub|0>>+\<beta\>=const
      </equation*>

      <item>BCs: <with|mode|math|z> depthward, <with|mode|math|z=\<eta\>>
      free surface

      <\itemize>
        <item><with|mode|math|\<varphi\><rsub|z>=0> at bottom

        <item><with|mode|math|<frac|\<mathd\>|\<mathd\>t>(z-\<eta\>)=0> at
        <with|mode|math|z=\<eta\>><with|mode|math|\<rightarrow\>><with|mode|math|\<partial\><rsub|t>\<varphi\>=\<partial\><rsub|t>\<eta\>>
        at <with|mode|math|z=0> (!).

        <item>pressure continuous at interface. Use Bernoulli's equation to
        rewrite as condition

        <\equation*>
          \<partial\><rsub|t>\<varphi\>+g\<eta\>=0<space|1em><with|mode|text|at>
          z=0.
        </equation*>
      </itemize>

      <item><em|Surface tension>: <with|mode|math|p<rsub|1>-p<rsub|2>=-\<gamma\>\<cdot\><with|mode|text|curvature>>.

      <item><em|Rayleigh-Taylor instability>: Large density over small
      density.

      <item><em|Kelvin-Helmholtz instability>: Wave formation.
    </itemize>
  </itemize>

  <section|Constitutive Laws>

  <\itemize>
    <item><em|Observer>: A reference frame/coordinate system w.r.t. which
    vectors and tensors are seen.

    <\equation*>
      \<b-x\><rsup|\<ast\>>=\<b-c\>(t)+Q(t)\<b-x\>
    </equation*>

    so, for example, <with|mode|math|F<rsup|\<ast\>>=Q*F>,
    <with|mode|math|J<rsup|\<ast\>>=J>, <with|mode|math|U<rsup|\<ast\>>=U>,
    <with|mode|math|R<rsup|\<ast\>>=Q*R>.

    <item><em|Objective fields>:

    <\eqnarray*>
      <tformat|<table|<row|<cell|\<varphi\><rsup|\<ast\>>(\<b-x\><rsup|\<ast\>>)>|<cell|=>|<cell|\<varphi\>(\<b-x\>)>>|<row|<cell|\<b-u\><rsup|\<ast\>>(\<b-x\><rsup|\<ast\>>)>|<cell|=>|<cell|Q\<b-u\>(\<b-x\>)>>|<row|<cell|A(\<b-x\><rsup|\<ast\>>)>|<cell|=>|<cell|Q*A(\<b-x\>)Q<rsup|T>>>>>
    </eqnarray*>

    Examples: <with|mode|math|D>, regions, normals,
    <with|mode|math|\<sigma\>>

    Non-examples: <with|mode|math|\<b-v\>=<wide|\<b-c\>|\<dot\>>+Q*\<b-v\>>,
    <with|mode|math|L=Q*L*Q<rsup|T>+<wide|Q|\<dot\>>Q<rsup|T>>,
    <with|mode|math|W>.

    <item><em|Constraint stress>:

    <\itemize>
      <item>Must be workless, i.e. <with|mode|math|tr(N*D)=0>

      <item>Constraint given as <with|mode|math|\<lambda\>(C)=0\<rightarrow\><wide|\<lambda\>|\<dot\>>=tr(\<lambda\><rsub|C><wide|C|\<dot\>>)=0>,
      where <with|mode|math|C=F<rsup|T>F>.

      <item><with|mode|math|<wide|C|\<dot\>>=2F<rsup|T>D*F><with|mode|math|\<Rightarrow\>><with|mode|math|N=\<alpha\>F\<lambda\><rsub|C>F<rsup|T>>.
    </itemize>

    <item><em|Fluid>: <with|mode|math|\<sigma\>=g(L)>.

    Cannot support shear stress at equilibrium. If <em|ideal>, also cannot
    support shear stress when in motion.

    <\itemize>
      <item><em|Objectivity>: <with|mode|math|\<sigma\><rsup|\<ast\>>=g(L<rsup|\<ast\>>)>.

      <\itemize>
        <item>Choose <with|mode|math|Q=Id>,
        <with|mode|math|<wide|Q|\<dot\>>=-W> to obtain that
        <with|mode|math|g(L)=g(D)>.
      </itemize>

      <item>Most general such <with|mode|math|g>:

      <\equation*>
        \<sigma\>(D)=\<alpha\>I+\<beta\>D+\<gamma\>D<rsup|2>,
      </equation*>

      with <with|mode|math|\<alpha\>>, <with|mode|math|\<beta\>>,
      <with|mode|math|\<gamma\>> functions of the invariants of
      <with|mode|math|D>.

      Proof: Cayley-Hamilton.

      <item><em|Incompressible fluid>:

      <\equation*>
        \<sigma\>=-p*Id.
      </equation*>

      <item><em|Ideal fluid>:

      <\equation*>
        \<sigma\>=-p(\<rho\>)Id
      </equation*>

      <item><em|Newtonian fluid>:

      <\equation*>
        \<sigma\>=-p(\<rho\>)Id+2\<mu\>D
      </equation*>

      <item><em|Navier-Stokes equation>:

      <\equation*>
        \<rho\>\<b-a\>=-\<nabla\>p+\<mu\>\<Delta\>\<b-v\>+\<rho\>\<b-b\>
      </equation*>

      plus conservation of mass.

      <item><em|Rescaling> <with|mode|math|<wide|x|~>=x/l>,
      <with|mode|math|<wide|\<b-v\>|~>=\<b-v\>/v>,
      <with|mode|math|p=p/(\<rho\><rsub|0>v<rsup|2>)>,
      <with|mode|math|<wide|t|~>=t/l>.

      Then <em|kinematic viscosity> is <with|mode|math|\<nu\>=\<mu\>/p>.

      <item><em|Reynolds number>: <with|mode|math|Re=l*v/\<nu\>>.

      High: Dominated by inertial effects.

      Low: Dominated by viscous effects.

      <item>No-slip BCs apply only for viscous fluids.

      <item><em|Wiggling plate>: Watch for emergence of a boundary layer.
    </itemize>

    <item><em|Solid>: <with|mode|math|\<sigma\>=f(F)>

    <\itemize>
      <item><em|Material Symmetry>: <with|mode|math|P\<in\>\<cal-S\>>, where
      <with|mode|math|\<cal-S\>> is the symmetry group of the material.

      <\equation*>
        \<sigma\>=f(F)=f(F*P)
      </equation*>

      <em|Isotropic Material>: <with|mode|math|\<cal-S\>=SO(3)>. Then choose
      <with|mode|math|P=R<rsup|T>><with|mode|math|\<Rightarrow\>><with|mode|math|\<sigma\>=f(F)=f(V)>.

      <item><em|Objectivity>: <with|mode|math|\<sigma\><rsup|\<ast\>>=f(V<rsup|\<ast\>>)>.

      Most general expression to satisfy this:

      <\equation*>
        \<sigma\>(V)=\<alpha\>Id+\<beta\>V+\<gamma\>V<rsup|2>,
      </equation*>

      with <with|mode|math|\<alpha\>>, <with|mode|math|\<beta\>>,
      <with|mode|math|\<gamma\>> functions of the invariants.

      <item><em|Lamé constant/Young's Modulus>: Linearization!

      <\eqnarray*>
        <tformat|<table|<row|<cell|F>|<cell|=>|<cell|Id+\<nabla\>\<b-u\>>>|<row|<cell|E>|<cell|=>|<cell|<frac|1|2><left|[>F<rsup|T>F-Id<right|]>\<approx\><frac|1|2><left|[>\<nabla\>\<b-u\>-(\<nabla\>\<b-u\>)<rsup|T><right|]>>>|<row|<cell|V>|<cell|\<approx\>>|<cell|Id+E>>|<row|<cell|R>|<cell|\<approx\>>|<cell|Id+<frac|1|2><left|[>\<nabla\>\<b-u\>-(\<nabla\>\<b-u\>)<rsup|T><right|]>>>>>
      </eqnarray*>

      Use these in

      <\equation*>
        \<sigma\>=c<rsub|0>tr V*Id+c<rsub|1>V+c<rsub|3>V<rsup|2>\<approx\>\<lambda\>tr(E)Id+2\<mu\>E,
      </equation*>

      where <with|mode|math|\<lambda\>>, <with|mode|math|\<mu\>> are the
      <em|Lamé constants>.

      <item><em|Strain energy per volume>:
      <with|mode|math|W(F)><with|mode|math|\<leftarrow\>>the usual way to
      specify constitutive relations for solids

      Then

      <\equation*>
        \<sigma\>=<frac|1|J>\<cdot\><wide*|<frac|\<partial\>W|\<partial\>F>|\<wide-underbrace\>><rsub|<frac|\<partial\>W|\<partial\>V>R>F<rsup|T>=<frac|1|J>\<cdot\><frac|\<partial\>W|\<partial\>V>V
      </equation*>

      Invoke objectivity: <with|mode|math|W(F)=W(U)>

      Invoke isotropy: <with|mode|math|W(F)=W(V)>

      <with|mode|math|\<Rightarrow\>><with|mode|math|W> depends only on
      invariants of <with|mode|math|V>.

      <with|mode|math|\<Rightarrow\>> <with|mode|math|W>'s principal axes
      line up with those of <with|mode|math|V>, i.e. <em|principal stresses>
      <with|mode|math|\<\|\|\>> <em|principal stretches>:

      <\equation*>
        \<sigma\><rsub|\<alpha\>>=<frac|1|J>\<lambda\><rsub|\<alpha\>><frac|\<partial\>W(\<lambda\><rsub|1>,\<lambda\><rsub|2>,\<lambda\><rsub|3>)|\<partial\>\<lambda\><rsub|\<alpha\>>>.
      </equation*>

      Incompressible:

      <\equation*>
        \<sigma\><rsub|\<alpha\>>=<frac|1|J>\<lambda\><rsub|\<alpha\>><frac|\<partial\>W(\<lambda\><rsub|1>,\<lambda\><rsub|2>,\<lambda\><rsub|3>)|\<partial\>\<lambda\><rsub|\<alpha\>>>-p.
      </equation*>

      Specifying <with|mode|math|W> in terms of
      <with|mode|math|B=F*F<rsup|T>>:

      <\equation*>
        \<sigma\>=<frac|2|J><left|(>III<rsub|B>W<rsub|III<rsub|B>>*Id+(W*<rsub|\<Iota\><rsub|B>>+\<Iota\><rsub|B>W<rsub|II<rsub|B>>)B-W<rsub|II<rsub|B>>B<rsup|2><right|)>,
      </equation*>

      where subscripts by <with|mode|math|\<Iota\><rsub|B>,II<rsub|B>,III<rsub|B>>
      mean partial derivatives.

      <item><em|Neo-Hookean material>:

      <\eqnarray*>
        <tformat|<table|<row|<cell|W>|<cell|=>|<cell|<frac|1|2>\<mu\><left|[>\<lambda\><rsub|1><rsup|2>+\<lambda\><rsub|2><rsup|2>+\<lambda\><rsub|3><rsup|2>-3-2ln(J)<right|]>+<frac|1|2>\<mu\><rprime|'>(J-1)<rsup|2>,>>|<row|<cell|<with|mode|text|unconstrained:><space|1em>\<sigma\><rsub|i>>|<cell|=>|<cell|\<mu\>(\<lambda\><rsub|i><rsup|2>-1)+\<mu\><rprime|'>J*(J-1),>>|<row|<cell|<with|mode|text|incompressible:><space|1em>\<sigma\><rsub|i>>|<cell|=>|<cell|\<mu\>\<lambda\><rsub|i><rsup|2>-p.>>>>
      </eqnarray*>

      <item>Solving a solids problem:

      <\itemize>
        <item>Calculate <with|mode|math|F> (Kinematics)

        <item>Calculate <with|mode|math|B=F*F<rsup|T>>

        <item>Calculate <with|mode|math|\<sigma\>>

        <item>Apply conservation of momentum in deformed configuration. Solve
        for unknowns <with|mode|math|\<b-x\>>, <with|mode|math|p>, using BCs.
      </itemize>
    </itemize>
  </itemize>
</body>

<\initial>
  <\collection>
    <associate|page-type|letter>
    <associate|par-first|0>
  </collection>
</initial>

<\references>
  <\collection>
    <associate|auto-1|<tuple|1|1>>
    <associate|auto-2|<tuple|2|1>>
    <associate|auto-3|<tuple|2.1|2>>
    <associate|auto-4|<tuple|2.1.1|2>>
    <associate|auto-5|<tuple|2.2|3>>
    <associate|auto-6|<tuple|3|3>>
    <associate|auto-7|<tuple|4|5>>
    <associate|auto-8|<tuple|5|8>>
  </collection>
</references>

<\auxiliary>
  <\collection>
    <\associate|toc>
      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|1<space|2spc>Tensor
      Stuff> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-1><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|2<space|2spc>Kinematics>
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-2><vspace|0.5fn>

      <with|par-left|<quote|1.5fn>|2.1<space|2spc>Static
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-3>>

      <with|par-left|<quote|3fn>|2.1.1<space|2spc>Static Examples
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-4>>

      <with|par-left|<quote|1.5fn>|2.2<space|2spc>Dynamic
      <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-5>>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|3<space|2spc>Balance
      Laws and Field Equations> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-6><vspace|0.5fn>

      <vspace*|1fn><with|font-series|<quote|bold>|math-font-series|<quote|bold>|4<space|2spc>Constitutive
      Laws> <datoms|<macro|x|<repeat|<arg|x>|<with|font-series|medium|<with|font-size|1|<space|0.2fn>.<space|0.2fn>>>>>|<htab|5mm>>
      <no-break><pageref|auto-7><vspace|0.5fn>
    </associate>
  </collection>
</auxiliary>