we are concerned with underlying connections between fluids , elasticity , isometric embedding of riemannian manifolds , and the existence of wrinkled solutions of the underlying nonlinear partial differential equations .
one of the main purposes of this paper is to develop such connections for the case of two spatial dimensions and examine whether the continuum mechanical equations can be mapped into a corresponding geometric framework and the inherent direct application of the theory of isometric embeddings and the gauss - codazzi equations .
another motivation for such study is to explore the possibility whether the geometric theory can serve an avenue for addressing the admissibility criteria for nonlinear conservation laws in continuum mechanics .    in recent years ,
a theory of _ wild _ solutions to both the incompressible and compressible euler equations in two and higher spatial dimensions has been developed by de lellis , szkelyhidi , and others in @xcite and the references cited therein .
the approach is based on the analogy with the highly irregular ( static ) solutions of the isometric embedding of a two - dimensional riemannian manifold into three - dimensional euclidean space given by the nash - kuiper theorem @xcite .
specifically , the analogy arises because of the applicability of gromov s h - principle and convex integration to both problems ( _ cf_. @xcite ) .
this suggests that perhaps the initial value problem in fluid dynamics and the embedding problem in differential geometry would be more than just two analogous issues , and could in fact be mapped one to the other . a first thought on
this issue would suggest that the question is not even meaningful : the fluid problem is dynamic and the embedding problem is static .
thus , if the linkage is to make any sense at all , we must think of the embedding problem as dynamic and derive the equations of a time evolving two - dimensional riemannian manifold .
of course , conceptually it is easy to visualize this time evolving two - dimensional riemannian manifold as a two - dimensional surface moving in three - dimensional space that can be seen in such an everyday phenomenon as the vibration of the surface of a drum .
thus , to make the link , we must derive the equations of a new type of geometric flow and then interpret the consequences of this fluid - geometric duality .
in fact , once we have mapped a solution of the euler equations onto the evolving surface , it is rather easy to see that the dynamic metric @xmath0 is a short metric in the sense of the nash - kuiper theorem @xcite with respect to a metric associated with a developable surface . hence , an immediate consequence of our theory is that the evolving manifold can be approximated by wrinkled manifolds ( _ i.e. _ , manifolds with discontinuous second derivatives ) and thus gives some indication of the existence of _ wild _ solutions of the dual fluid problem
. a simple mental picture of the geometric image of the fluid problem would be the motion of @xmath1 wrinkles on a piece of paper .
moreover , the _ wild _ solutions of the geometric problem are completely time - reversible and are in reality a sequence of the nash - kuiper solutions of the embedding problem with the metric given by @xmath2 , where @xmath3 corresponds to a developable surface that is time - independent and has non - vanishing mean curvature .
the implication of this fact is immediate : the geometric nash - kuiper _ wild _ solutions are time - reversible .
thus , the question raised in @xcite as to which is the correct admissibility criterion to choose the relevant solution from the infinite number of non - unique solutions to the euler equations , has at least a partial answer .
namely , no dynamic admissibility criterion such as the energy inequality , entropy inequality , entropy rate can serve the purpose , since all the inequalities become the identities for such solutions .
the only possible useful criteria must be meaningful for time reversible fluid flow , such as energy minimization , artificial viscosity @xcite , or viscosity - capillarity @xcite .
the equations for geometric flow are abundant : the einstein equations of general relativity and the ricci flow equations are two of the better known examples . in both of these problems ,
the metric for a riemannian manifold becomes the dynamic unknown .
in the theory developed in this paper , the same situation arises : a dynamic metric @xmath0 is our unknown along with the second fundamental form for the evolving manifold .
furthermore , just as in the case of the einstein equations ( _ cf .
_ @xcite ) , the initial data must be consistent with the problem , that is , initially an embedded manifold does exist and the map is also consistent with the divergence free condition on the velocity for the incompressible fluid case . for the einstein equations , this consistency of the initial data yields the einstein constraint equations , while in our case a system of constraint equations is also required . in this paper
we show that our constraint equations have local analytic solutions ; moreover , there is a velocity @xmath4 that defines both a solution of the incompressible euler equations and an evolving two - dimensional riemannian manifold isometrically immersed in @xmath5 .
a similar result is given for the compressible case , as well as for neo - hookean elasticity .
the most important issue is the physical meaning of the results mentioned above . in short ,
we emphasize the comments which have appeared in szkelyhidi @xcite , bardos - titi - wiedermann @xcite , bardos - titi @xcite , and bardos - lopes filho - niu - nussenzveig lopes - titi @xcite .
the appearance of the _ wild _ solutions is due essentially to the application of the euler equations on the initial data that form a discontinuous shear , _
i.e. _ , a vortex sheet .
hence , the euler equations , which have no small scales built into their theory as opposed to the compressible or incompressible navier - stokes equations , have been used in a case where they should not be expected to apply . to this
, we must add the proviso that as again mentioned in de lellis - szkelyhidi @xcite ; also see elling @xcite : these _ wild _ solutions could be a demonstration of fluid turbulence .
the results given in this paper show that , if this is the case , then this _ wild _ fluid turbulent behavior is mirrored by the solutions produced by the non - smooth moving wrinkled surfaces produced by the nash - kuiper theorem @xcite .
this paper consists of nine sections after this brief introduction . in  2 , we recall the euler equations for an inviscid incompressible fluid and the gauss - codazzi equations for the isometric embedding of a two - dimensional riemannian manifold into three - dimensional euclidean space . in
3 , as discussed above , we derive the equations of the geometric flow , as well as the constraint conditions on the initial data . in  4 , we prove the solvability of the constraint equations for the initial data .
we emphasize that the constraint conditions require the initial fluid velocity _ not _ to be a shear flow .
this assumption pairs nicely with a theorem in de lellis - szeklyhidi @xcite that non - smooth shear flow initial data for both the compressible and incompressible euler equations are _ wild _ data and yield non - unique solutions to the cauchy problem for both the compressible and incompressible euler equations . in fact
our computation suggests that this is the only _ wild _ initial data .
on the other hand , we note that , for the degenerate case , when the fluid motion is a shear flow , we still have a metric that provides the desired map , namely , metric @xmath3 . in  5 , we state and
prove our main result : evolving from the initial data , there exists a metric @xmath0 of geometric flow equations which yields a solution of both the euler equations and the equations describing the evolving isometrically immersed riemannian manifold @xmath6 . in  6
, we continue the discussion of the initial data issue and show that , for the shear flow initial data in the hypotheses of lemma [ l41 ] , the symbol of the underlying system of second order partial differential equations vanishes so that the equations for @xmath0 are degenerate ; however , as we just commented , metric @xmath3 suffices in this singular case . in  7
, we give our principal result : the evolving manifold arising from the euler equations can be approximated by wrinkled @xmath7 manifolds @xmath8 , @xmath9 .
in fluid variables , these wrinkled manifolds correspond to discontinuous shear flows , which thus show conceptually that the discontinuous shear flows are dense in the set of solutions of the incompressible euler equations .
furthermore , these wrinkled solutions can be arranged as a time - sequence of solutions which render the initial value problem for the evolving manifold to have an infinite number of constant energy solutions . in  8 ,
a short discussion is provided to show that many of our earlier results for the incompressible euler equations carry over to the compressible case .
based on the knowledge that has been obtained from the fluid equations , we show in  9 how the case of general continuum mechanics can be placed in our mechanics - geometry framework .
as an illustrative example , we demonstrate results for elastodynamic motion of a neo - hookean solid . finally ,
our last section ,  10 , provides a discussion of admissibility criteria for the ( incompressible and compressible ) euler equations which asserts that , for the multidimensional euler equations , no dynamic admissibility condition would eliminate _ wild _ solutions , and hence the only meaningful one must be the one which eliminates or at least reduces the number of wrinkles of our dual geometric problem . furthermore , these results suggest that more refined continuum mechanical theories relying on micro - structure could play a key role in choosing admissible solutions .
we start with some basic geometric equations and notations for subsequent developments . for more details , see han - hong @xcite and the references cited therein .
let @xmath6 be a two - dimensional riemannian manifold with @xmath10 denoting a point on the manifold , @xmath11 .
the unit normal vector @xmath12 to the manifold is given by @xmath13 and the second fundamental form is @xmath14 where @xmath15 , @xmath16 , and @xmath17 with @xmath18 for @xmath19 .
we will use the alternative version of the second fundamental form @xmath20 and recall that @xmath21 where @xmath22 is the gauss curvature , and @xmath23 is the riemann curvature tensor .
furthermore , for notational simplicity , we henceforth drop the  ' " superscript in the alternative version of the second fundamental form .
then the codazzi equations are @xmath24,\\[1 mm ] \gamma_{12}^{1 } = \frac1{2\det { { \bf g}}}\left[g_{22}(\partial_{2}g_{11})-g_{12}(\partial_{1}g_{22})\right],\\[1 mm ] \gamma_{22}^{1 } = \frac1{2\det { { \bf g}}}\left[g_{22}(2\partial_{2}g_{21}-\partial_{1}g_{22})-g_{12}(\partial_{2}g_{22})\right],\\[1 mm ] \gamma_{11}^{2 } = \frac1{2\det { { \bf g}}}\left[-g_{12}(\partial_{1}g_{11})+g_{11}(2\partial_{1}g_{12}-\partial_{2}g_{11})\right],\\[1 mm ] \gamma_{12}^{2 } = \frac1{2\det { { \bf g}}}\left[-g_{12}(\partial_{2}g_{11})+g_{11}(\partial_{1}g_{22})\right],\\[1 mm ] \gamma_{22}^{2 } = \frac1{2\det { { \bf g}}}\left[-g_{12}(2\partial_{2}g_{21}-\partial_{1}g_{22})+g_{11}(\partial_{2}g_{22})\right ] , \end{cases } \ee and the riemann curvature tensor is given by \begin{equation}\label{2.4-a } r_{iljk}=g_{lm}\big(\partial_{k}\gamma_{ij}^{m } -\partial_{j}\gamma_{ik}^{m}+\gamma_{ij}^{n}\gamma_{nk}^{m } -\gamma_{ik}^{n}\gamma_{nj}^{m}\big),\ ] ] where we have used the einstein summation convention that the repeated indices are implicitly summed over in the terms , which will also be used from now on .
in particular , we have @xmath25    with , we have gauss s theorem egregium for the gauss curvature . a convenient form is given by brioschi s formula : [ e24 ]    = & 12(*g*)^2    - _
22g_11 + 2_12g_12- _ 11g_22 & _ 1g_11 & 2_1g_12- _ 2g_11 + 2_2g_12-_1g_22 & 2g_11 & 2g_12 + _ 2g_22 & 2g_12 & 2g_22     + & -12(*g*)^2    0 & _ 2g_11 & _ 1g_22 + _ 2g_11 & 2g_11 & 2g_12 + _ 1g_22 & 2g_12 & 2g_22    .
we recall the fundamental theorem of surface theory states that the solvability of the gauss - codazzi equations is a necessary and sufficient condition for the existence of an isometric embedding , _
i.e. _ , a simply connected surface @xmath26 which satisfies @xmath11
. a convenient reference for the smooth version of the fundamental theorem is do carmo @xcite , while a non - smooth version can be found in mardare @xcite .
the equations for the balance of linear momentum are [ e25 ]    _ 1(u^2+p)+_2(uv)=-_tu , + _ 1(uv)+_2(v^2+p)=-_tv ,    where the constant density @xmath27 is taken .
the condition of incompressibility for constant density is then given by the equation : [ e26 ] _ 1u+_2v=0 .
in addition , taking the divergence of equations and using the incompressibility condition , we have [ e27 ] _ 11(u^2 ) + 2_12(uv)+_22(v^2 ) = -p .      just as in chen - slemrod - wang @xcite , it is convenient to write the geometric equations in fluid variables
. set [ e28 ] l = v^2+p , m =- uv , n = u^2+p .
then the gauss equation becomes @xmath28 that is , @xmath29 where @xmath30 .
this quadratic equation then gives [ e29 ] p=-12q^212 .
this means that @xmath31 must be required .
as it will be seen in the analysis below , this condition is always satisfied .
we have just shown that @xmath32 can be written in fluid variables .
now we can write  the fluid variables in terms of the geometric variables @xmath32 . to do this , simply substitute formula into to find @xmath33 write @xmath34 to see that @xmath35 , which yields that @xmath36
. then we see @xmath37 and , using @xmath38 , @xmath39 this shows that the `` + '' sign must be chosen in the above formulas so that [ e210 ]    u^2=12(-(l - n)+ ) , + v^2=12((l - n)+ ) .
note that the condition : @xmath40 , with @xmath41 , is equivalent to @xmath42 that is , @xmath43 which is always satisfied .
thus , we have shown that @xmath44 are determined by @xmath32 , since @xmath45
we now construct a _ dual _ solution , which simultaneously satisfies the incompressible euler equations and the gauss - codazzi equations of isometric embeddings .    as before ,
the gauss - codazzi equations are [ e31 ]    _ 1n-_2 m = -_22 ^ 1l+2_12 ^ 1m-_11 ^ 1n , + _
1m-_2l=_22 ^ 2l-2_12 ^ 2m+_11 ^ 2n ,    and [ e32 ] ln - m^2=. hence , for any @xmath4 to be a solution of the euler equations  , we must have [ e33 ]    _ tu=_22 ^ 1l-2_12 ^ 1m+_11 ^ 1n , + _
^ 2l-2_12 ^ 2m+_11 ^ 2n ,    and [ e34 ] _ 1u+_2v=0 . taking the divergence of and using , we have [ e35 ] _ 1(_22 ^ 1l-2_12 ^ 1m+_11 ^ 1n ) + _ 2(_22 ^ 2l-2_12 ^ 2m+_11 ^ 2n)=0 . for convenience , define [ e36 ]    & u(l , m , n):=(12)^12 , + & v(l , m , n):=(12)^12 ,    so that [ e37 ] ( u , v)=(u(l , m , n),v(l , m , n ) ) . the other choice of @xmath46 can be handled similarly .    in summary
, we have the two evolution equations and four closure relations  and .
since @xmath47 and @xmath22 are the functions of @xmath48 through  , we obtain the six equations  and for the six unknowns @xmath48 .
in this section , we exposit the constraint equations on the initial data and the consequences of their solvability .
our first result is    [ theorem41 ] assume that the initial data @xmath48 at @xmath49 satisfy the five constraint equations : @xmath50 which mean that the initial data are consistent with both the incompressibility of the fluid and that the evolving manifold is initially indeed a riemannian manifold .
then , if the system of six equations  and in the six unknowns @xmath48 is satisfied , it produces simultaneously a solution of both the gauss - codazzi equations and the incompressible euler equations .
since @xmath51 , we see from the first two evolution equations and that @xmath52 since the initial data satisfy the constraint equations , we conclude that , for @xmath53 , @xmath54    next , since the gauss - codazzi equations are satisfied for all @xmath53 , the fluid representation for @xmath55 allows to write the codazzi equations as @xmath56 and the gauss equation as @xmath57 by , we see that the balance of linear momentum equations is also satisfied .
the proof is complete .
we next examine the solvability of the initial data system .
note that  imply [ e46 ] _
1(_1n-_2 m ) + _ 2(-_1m+_2l)=0 , and the initial data system can be written as @xmath58 we can reverse the above computation . if  are satisfied , we take the divergence of the left - hand sides of  and employ to yield .    from , we find that and  becomes an undetermined system of five equations in the six unknowns @xmath48 .
our existence result for the initial data satisfying the constraint equations reads as follows :    [ l41 ] let an analytic divergence free velocity @xmath59 be prescribed in a neighborhood of a point @xmath60 such that @xmath61 at this point .
set @xmath62 .
on @xmath63 ( respectively @xmath64 , prescribe the analytic initial data : @xmath65 , and @xmath66 satisfying the ordinary differential equation : @xmath67 with initial condition @xmath68 at @xmath69 so that @xmath70 at @xmath71 ( respectively , @xmath72 and a similar ordinary differential equation and initial data ) .
then the initial data system  has a local analytic solution @xmath73 .
we divide the proof into eight steps .    in fluid variables ,
three of our equations are @xmath74 next , prescribe the velocity to make this a determined system . if a divergence free velocity @xmath59 is prescribed , then is immediately solvable for @xmath75 under the standard regularity assumptions on @xmath59 , and hence defines the gauss curvature @xmath22 .
thus , @xmath55 are known with no dependence on metric @xmath0 .
now , determine @xmath0 by solving the three equations : @xmath76    we use  to solve for @xmath77 and @xmath78 .
simply use to write these two equations  as [ e421 ]    & 1*g * + & = _
1n-_2 m + 12*g*l + & -m + 12*g * n ,    [ e422 ]    & 1*g * + & = -_1m+_2l + 12*g*l + & - m + 12*g*n .    in the matrix form , we have    1*g *    g_12n & -g_22l + -g_11n & g_12l    _ 1g_12 + _ 2g_21    =    _
1n-_2 m + -_1m+_2l    + 1*g *    g_1 + g_2    ,    where    g_1=&l -m + & + n , + g_2=&l -m + & + n .    the inverse of the coefficient matrix is @xmath79 \frac{g_{11}}{l } & \frac{g_{12}}{l } \end{bmatrix},\ ] ] which gives [ e423 ]    _ 1g_12 + _ 2g_21    =    & + &    _ 1n-_2m+1*g*g_1 + -_1m+_2l+1*g*g_2    .    from , the equality of cross partials gives us the additional consistency equation : [ e424 ]    & _ 2((_1n-_2m+1*g*g_1 ) + ( -_1m+_2l+1*g*g_2 ) ) + & = _
1((_1n-_2m+1*g*g_1 ) + ( -_1m+_2l+1*g*g_2 ) ) .
the gauss curvature equation has the form : [ e425 ] = 1*g*+l.o.t .
substitution of into yields the two second - order equations for @xmath80 .
now compute the highest order terms of the differential operators associated with ( [ e423 ] ) at the origin with @xmath70 : @xmath81 \partial_{2}g_{21}\end{bmatrix}&= \begin{bmatrix } 0 & \frac{1}{n } \\[1.5 mm ] \frac{1}{l } & 0\end{bmatrix}\begin{bmatrix } \mathbf{\ } \partial _ { 1}n-\partial _ { 2}m-{\frac12}\partial_{1}g_{22}l-\partial_{2}g_{11}m+{\frac12}\partial _ { 1}g_{11}n \\[1.5 mm ] \mathbf{-}\partial _ { 1}m+\partial _ { 2}l+{\frac12}\partial_{2}g_{22}l-\partial _ { 1}g_{22}m+-{\frac12}\partial _ { 2}g_{11}n\end{bmatrix}\\[2 mm ] & = \begin{bmatrix } \frac{1}{n}(\mathbf{-}\partial _ { 1}m+\partial _ { 2}l+{\frac12}\partial_{2}g_{22}l-\partial_{1}g_{22}m-{\frac12}\partial_{2}g_{11}n ) \\[1.5 mm ] \frac{1}{l}\mathbf{\ ( } \partial_{1}n-\partial_{2}m-{\frac12}\partial _ { 1}g_{22}l-\partial _ { 2}g_{11}m+{\frac12}\partial _ { 1}g_{11}n)\end{bmatrix}. \end{split}\ ] ] then becomes [ e426 ]    = & + + & + + l.o.t .
and equation becomes @xmath82    in the matrix form ,  are written as [ e427 ]    14 & -14 + -12 & 12
_ 11    g_11 + g_22    +    - & - + & -    _ 12    g_11 + g_22    +    -14 & 14 + -12 & 12    _ 22    g_11 + g_22    = l.o.t .
make the change of independent variables : @xmath83 so that , for any @xmath84 , @xmath85    then we have [ e427-a ] a    g_11 + g_22    + a    g_11 + g_22    + ( + 1 )    & - + -1 & 1    g_11 + g_22    = l.o.t . where the coefficient matrix @xmath86 is @xmath87 -{\frac12}\frac{n}{l } & { \frac12}\end{bmatrix}+\begin{bmatrix } -\frac{m}{2l } & -\frac{m}{2n } \\[1.5 mm ]
\frac{m}{l } & -\frac{m}{n}\end{bmatrix}+\begin{bmatrix } -{\frac14 } & { \frac14}\frac{l}{n } \\[1.5 mm ] -{\frac12 } & { \frac12}\frac{l}{n}\end{bmatrix}= \begin{bmatrix } { \frac14}\frac{n}{l}-\frac{m}{2l}-{\frac14 } & -{\frac14}-\frac{m}{2n}\text { } + { \frac14}\frac{l}{n } \\[1.5 mm ] -{\frac12}\frac{n}{l}+\frac{m}{l}-{\frac12 } & { \frac12}-\frac{m}{n}+{\frac12}\frac{l}{n}\end{bmatrix}\ ] ] with determinant equal to @xmath88    if the term : @xmath89 , then @xmath90 .
thus , the conditions that @xmath91 and @xmath92 make this coefficient matrix non - singular . since @xmath93 in fluid variables ,
we need @xmath61 and curvature @xmath22 determined by the poisson equation with @xmath94 .    of course ,
when @xmath95 , @xmath96 , and @xmath97 , then @xmath98 .
clearly , if pressure @xmath75 determined by the poisson equation is positive , then the condition : @xmath91 , is automatically satisfied .
then line @xmath99 is non - characteristic for the equation in .
since @xmath100 , gives a partial differential equation for @xmath66 in the direction normal to the initial data line @xmath101 : @xmath102    we now determine the initial data for system and on line @xmath63 that is non - characteristic for the system .
since @xmath103 the initial conditions : @xmath104 also give us @xmath105 this yields @xmath106 then , from the formula : @xmath107 , , and @xmath108 \partial_{2}g_{21}\end{bmatrix}= \begin{bmatrix } \frac{g_{12}}{n } & \frac{1}{n } \\[1.5 mm ] \frac{1}{l } & \frac{g_{12}}{l } \end{bmatrix}\begin{bmatrix } \mathbf{\ } \partial_{1}n-\partial_{2}m \\[1.5 mm ] \mathbf{-}\partial_{1}m+\partial_{2}l \end{bmatrix},\ ] ] we obtain the ordinary differential equation on line @xmath63 . this differential equation and the initial condition : @xmath68 at @xmath109
determine the initial data for @xmath66 on line @xmath63 .
therefore , we have specified analytic initial data on the non - characteristic line @xmath63 for system and .
finally , we see that the cauchy - kowalewski theorem delivers a local analytic solution @xmath73 of system and with the analytic data on the non - characteristic line @xmath63 , where we have used the analyticity of @xmath75 as the solution to our poisson equation and the analyticity of @xmath22 from the formula : @xmath110 .
the positivity condition on @xmath75 is irrelevant and can always be satisfied by simply adding to any @xmath75 solving the poisson equation a sufficiently large positive constant .
this completes the proof .
notice that the shear flow data @xmath111 are divergence free , no matter how smooth or in fact how irregular @xmath112 is . in this case ,
@xmath113 , so that lemma [ l41 ] does not apply to yield a metric @xmath0 . however , in this case , we have the following trivial result .    [ l42 ] if @xmath113 in any open set of @xmath114 , then there exists a developable surface @xmath115 with corresponding metric @xmath3 given by @xmath116 with @xmath117 similarly , the corresponding developable surface with corresponding metric @xmath3 can be obtained for @xmath118 as another shear flow .
in particular , @xmath119 and @xmath120 for @xmath121 .
this metric @xmath3 yields a solution of our initial data system .
furthermore , this solution must be a shear flow @xmath122 , or @xmath123 , owing to the incompressibility .
consider the surface : @xmath124 then @xmath125 so that @xmath126 and @xmath127    furthermore , we have @xmath128 so that @xmath129    at @xmath130 , @xmath131 set @xmath132 so that @xmath133 .    then , at @xmath134 , @xmath135 which implies that @xmath136    with @xmath137 prescribed by the solution of @xmath138 we have constructed the developable surface , while the solution of is determined by @xmath139    similarly , we can obtain the corresponding results for @xmath118 as another shear flow .    [ r41 ] notice the very important fact that the shear flow initial data given by lemma [ l42 ] require only the square integrability of the shear functions @xmath140 and @xmath141
.    [ r42 ] the assumption of analyticity is made only for convenience .
in fact , if one follows the presentation of deturck - yang @xcite for a similar problem where the theory of systems of real principal type is exploited , it is quite evident that the analyticity could be replaced by @xmath142 .
now we prove our main existence theorem .    [ t51 ] for the initial data constructed in lemma [ l41 ] ,
there exists a solution of system  and , locally in space - time , which satisfies both the incompressible euler equations  and a propagating two - dimensional riemannian manifold isometrically immersed in @xmath5 .    for the given fluid data , by the cauchy - kowalewski theorem
, there exists a local space - time solution @xmath143 of the incompressible euler equations  with @xmath144 , and .
define the gauss curvature via the relation : @xmath145 so that @xmath146 .
thus system  and is satisfied if the following system can be solved :    _ tu=_22 ^ 1l-2_12 ^ 1m+_11 ^ 1n , + _
tv=_22 ^ 2l-2_12 ^ 2m+_11 ^ 2n , + * g*=g_2m(_2_11^m-_1_12^m + _ 11^n_n2^m-_12^n_n1^m )    for the time evolving metric @xmath0 , where    u^2=12 , + v^2=12 , + q^2=.    this is the same system for @xmath0 that has been solved in lemma [ l41 ] .
this completes the proof .
[ c51 ] let @xmath4 be a solution , local analytic in space - time , of the incompressible euler equations
, satisfying the initial condition @xmath61 at a point @xmath147 . then there is an evolving metric @xmath0 so that solution @xmath4 of the euler equations also defines a time evolving two - dimensional riemannian manifold @xmath6 immersed in @xmath5 .
[ c52 ] conversely , let @xmath0 be an isometrically embeddable manifold in @xmath5 with second fundamental form @xmath148,m[{{\bf g}}],n[{{\bf g}}])$ ] that is a solution of the gauss - codazzi equations for an as yet undetermined metric @xmath0 ) and @xmath149 .
identify the fluid variables @xmath44 by our relations : [ 5.1-a ]    u^2= , + v^2= , + p^2+pq^2= , q^2=u^2+v^2 ,    so that @xmath150 , v[{{\bf g } } ] , p[{{\bf g}}])$ ] are functionals of the unknown metric field @xmath0 . then the equations : [ 5.2-a ]    _ tu = _ 22 ^ 1l-2_12 ^ 1m+_11 ^ 1n , + _ tv = _ 22 ^ 2l-2_12 ^ 2 m + _ 11 ^ 2n , + _
11(u^2)+2_12(uv)+_22(v^2 ) = -p    define a system of three equations in the three unknowns @xmath73 .
if this implicit system for @xmath0 has a solution , then the three functions @xmath44 satisfy the incompressible euler equations .
substitute the first two equations in into the codazzi equations .
we see that the balance of linear momentum is satisfied , and the third equation in implies , via the balance of linear momentum , the incompressibility condition : @xmath151    the difficulty in employing corollary [ c52 ] would be that the equations for @xmath0 are only known implicitly .
notice that equations
satisfied by @xmath0 can be written as [ e51 ]    _ t u=_22 ^ 1(v^2+p)+2_12 ^ 1(uv)+_11 ^ 1(u^2+p ) , + _
t v=_22 ^ 2(v^2+p)+2_12 ^ 2(uv)+_11 ^ 2(u^2+p ) , + _
1(_22 ^ 1(v^2+p)+2_12 ^ 1(uv)+_11 ^ 1(u^2+p ) ) + + _
2(_22 ^ 2(v^2+p)+2_12 ^ 2(uv)+_11 ^ 2(u^2+p))=0 .    here
@xmath4 are defined by in terms of @xmath55 which in turn depend on @xmath0 , and @xmath75 is given by @xmath152 from .
thus , this is a system of three equations for the three components of @xmath0 .
the difficulty here is that the dependence of @xmath59 on @xmath0 is implicit
. it would be interesting to identify some suitable ways to check their solvability .
as we noted in remark [ r41 ] , the shear flow data do not satisfy the hypothesis : @xmath61 in lemma [ l41 ] .
hence , this is a particularly interesting case to study .
namely , if the fluid data of @xmath4 are @xmath153 , the question is whether a riemannian manifold @xmath6 locally immersed in @xmath5 can be determined . for this data
, we find immediately that @xmath154 so that @xmath75 is harmonic .
furthermore , @xmath155 , and our constraint equations  reduce to the equations : [ e61 ]    & _
1p=-_22 ^ 1p-_11 ^ 1(u^2(x_2)+p ) , + & _
2p=-_22 ^ 2p-_11 ^ 2(u^2(x_2)+p ) ,    and [ e62 ]    pu^2(x_2)+p^2=&12(*g*)^2    -_22g_11 + 2_12g_12-_11g_22 & _ 1g_11 & 2_1g_12-_2g_11 + 2_2g_12-_1g_22 & 2g_11 & 2g_12 + _ 2g_22 & 2g_12 & 2g_22     + & -12(*g*)^2    0 & _ 2g_11 & _ 1g_22 + _ 2g_11 & 2g_11 & 2g_12 + _ 1g_22 & 2g_12 & 2g_22    .    as in
 [ s4 ] ,  give equation that is now [ e63 ]    & -14 + - & 12    _ 11    g_11 + g_22    +    -14 & + -12 &    _ 22    g_11 + g_22    = l.o.t . the symbol of the above differential operator is given by    & (    & -14 + - & 12    _
1 ^ 2 +    -14 & + -12 &    _ 2 ^ 2 ) + & =    _ 1 ^ 2 -14_2 ^ 2 & -14_1 ^ 2+_2 ^ 2 + -_1 ^ 2 - 12_2 ^ 2 & 12_1 ^ 2 + _ 2 ^ 2     + & = 0 .    thus , for the shear flow
initial data , our system of partial differential equations for @xmath0 is degenerate . on the other hand , we know from lemma [ l42 ] that , if a shear flow is given in an open neighborhood of space - time , then metric @xmath3 provides the desired map .
now that the existence of a locally analytic riemannian manifold with metric @xmath0 has been established in corollary [ c51 ] , we wish to know the local form of the metric . to do this
, we adapt a new local system of normal coordinates @xmath156 , @xmath157 ( _ cf .
_ @xcite and @xcite ) so that the origin of the old coordinate system remains the origin of our new coordinate system and @xmath158 .
we can easily compute @xmath0 to the leading order in the normal coordinates by the formula : @xmath159 that is ,    & g_11(x_1^,x_2^,t ) = 1 + 13r_1221x_2^x_2^+h.o.t .
= 1 - 13 ( x_2^)^2+h.o.t . , + & g_22(x_1^,x_2^,t ) = 1 + 13r_2112x_1^x_1^+h.o.t .
= 1 - 13 ( x_1^)^2+h.o.t . , + & g_12(x_1^,x_2^,t ) = 13r_1212x_1^x_2^+h.o.t .
= 13 x_1^x_2^+h.o.t ..    hence , as a quadratic form , @xmath160    however , the matrix @xmath161 -{\frac12}x_{1}^{\prime \prime } x_{2}^{\prime \prime } & ( x_{1}^{\prime \prime})^2\end{bmatrix}\ ] ] has determinant @xmath162 and hence is positive definite away from @xmath163    furthermore , on line @xmath164 , the quadratic form is @xmath165 thus , we have proved that map @xmath166 locally induced by metric @xmath0 is shorter than that induced by @xmath3 or , in language of the nash - kuiper theorem @xcite , embedding @xmath167 induced by @xmath0 is a _
short embedding_. we state this in the following lemma :    [ l71 ] the metric , @xmath168 , given in theorem [ t51 ] and corollary [ c51 ] is locally a short embedding that satisfies @xmath169 in the sense of quadratic forms .
we then recall the nash - kuiper theorem as given in @xcite .
[ t71 ] let @xmath170 be a smooth , compact riemannian manifold of dimension @xmath171 .
assume @xmath172 is in @xmath142 .
then    if @xmath173 , any short embedding can be approximated by isometric embeddings of class @xmath142 ( cf .
nash @xcite and gromov @xcite ) ;    if @xmath174 then any short embedding can be approximated in @xmath175 by isometric embeddings of class @xmath176 ( cf .
nash @xcite , kuiper @xcite ) .
we note that part ( ii ) of theorem [ t71 ] has been extended for an analytic metric @xmath172 by borisov @xcite to obtain the greater regularity @xmath177 , and the condition of analyticity has been weakened by conti - delellis - szelelyhidi @xcite .
we are now interested in case ( ii ) with @xmath178 , and @xmath179 so that the analyticity of the metric is immediate and @xmath9 .
thus , an immediate consequence of the results by nash - kuiper @xcite , borisov @xcite , and conti - delellis - szelelyhidi @xcite is the following theorem .
[ t72 ] the short embedding @xmath180 established by theorem [ t51 ] and corollary [ c51 ] can be approximated in @xmath175 , by wrinkled embeddings @xmath181 in @xmath7 , @xmath9 , locally in space - time . in elementary terms , the geometric image of our locally analytic solution to the incompressible euler equations , i.e. , our surface propagating in time can be approximated in @xmath1 , @xmath9 , by @xmath7 wrinkled manifolds .
[ t74 ] there are an infinite number of ( non @xmath182 ) evolving manifolds for metric @xmath179 , all with the same initial wrinkled data at @xmath49 .
furthermore , for all these solutions , the energy remains constant : @xmath183 where the manifold is defined by the @xmath1solution of @xmath184 which trivially satisfies @xmath185 so that @xmath186 is in @xmath187 in time .    on any sequence of time interval , just switch back and forth from any of the infinite choices of wrinkled solutions .
specifically , let @xmath188 $ ] be a time interval in which there exists a unique analytic solution to the cauchy initial value problem for the euler equations . by theorem [ t72
] , this solution is the pre - image of a short embedding @xmath180 and can be approximated by wrinkled embedding @xmath181 .
divide the interval @xmath188 $ ] into sub - intervals @xmath189 , [ t_1 , t_2 ] , \dots , [ t_{n-1 } , t_n]$ ] with @xmath190 and @xmath191 .
take any sequence @xmath192 .
then , by theorem [ t72 ] , we have a sequence of wrinkled solution @xmath193 such that @xmath194 now define the wrinkled solution @xmath195 on @xmath188 $ ] by @xmath196    fix @xmath197 as the cauchy data , but allow @xmath198 and @xmath199 to vary for @xmath200 . in this way
, we have produced an infinite number of wrinkled solutions satisfying the same cauchy data .
then the solutions are in @xmath1 in space on each sub - interval @xmath188 $ ] for time , so that they are in @xmath201 ; c(\omega))$ ] , but not necessarily in @xmath202 ; c(\omega))$ ] .    to compute the energy , use @xmath203    in the results given in theorems [ t71][t74 ] , by construction , embedding @xmath204 is the geometric image of a smooth solution of the euler equations . on the other hand , the nash - kuiper @xmath1 embeddings have not been shown to be the image of solutions in any sense of the euler equations .
we now give a definition of weak solutions to the euler equations
, which provides this link .
let @xmath181 be a nash - kuiper solution , @xmath205 , @xmath206 , of the isometric embedding equations : @xmath207 in an open domain .
since @xmath208 , it is in @xmath209 . by equality of cross partials
, we have @xmath210 in the sense of distributions , where the second partials of @xmath181 lie in @xmath211 .
the motivation behind this definition is as follows : for nash - kuiper solutions , we know that the normal to the second fundamental form is not well defined so that the gauss - codazzi equations that are the geometric image of the euler equations for smooth flows are also not well defined
. therefore , we retreat this one step back and use the equality of cross partials which is the starting point for the derivation of the gauss - codazzi equations .
[ weak_shear_embedding ] we define an embedding @xmath212 , @xmath213 , as a _ weak euler shear solution inducing embedding _ if ( [ cross_partials ] ) is satisfied in the sense of distributions .    certainly , the nash - kuiper embeddings satisfy this definition .
the arguments we have used for the incompressible case also carry over to the compressible euler equations .
then the equations for the balance of linear momentum are [ e81 ]    _ 1(u^2+p)+_2(uv)=-_t(u ) , + _
1(uv)+_2(v^2+p)=-_t(v ) ,    the equation for the conservation of mass is given by the equation : [ e82 ] _ 1(u)+_2(v)=-_t .
set [ e83 ] l = v^2+p , m =- uv , n = u^2+p , and the gauss equation becomes @xmath214 _
i.e. _ , [ e84 ] p^2+pq^2= , where @xmath30 . for the compressible case
, we take [ e85 ] p= , 1 . substitute into to obtain [ e86 ] ( ) ^2 + q^2=.    for simplicity , let us first take the isothermal case @xmath215 so that [ e87 ] ^2(1 + q^2)= , that is , [ e88 ] = . from , @xmath216 as the density
is determined explicitly as a function of @xmath22 and @xmath217 . on the other hand
, @xmath216 must satisfy the mass balance equation so that [ e89 ] _
1(u)+_2(v ) = -_t ( ) .
we see that , for given @xmath59 , is a scalar conservation law for the gauss curvature @xmath22 .
this leads us to the following elementary lemma :    [ l81 ] for a given smooth solution @xmath218 of the compressible euler equations  , the initial value problem for the gauss curvature @xmath22 given by with initial data @xmath219 at @xmath49 has a global smooth solution @xmath22 in space - time , which satisfies .
notice that equation is a linear transport equation for @xmath220 of conservation form .
then the result holds simply by using the standard local existence - uniqueness theorem for the conservation law and the direct computation @xmath221 this equation has the unique solution @xmath222    the general @xmath223 case follows from noting that the left - hand side of is a monotone increasing function of @xmath216 so that , for any @xmath224 , there exists a function @xmath225 that solves .
the regularity of this solution is seen by differentiation of with respect to @xmath226 .
we can thus state the general version of lemma [ l81 ] as    [ l82 ] for a given smooth solution @xmath218 of the compressible euler equations  , the initial value problem for the gauss curvature for with initial data at @xmath49 has a global solution in space - time satisfying .
again , we can also write the fluid variables in terms of the geometric variables @xmath55 . to do this , simply substitute with @xmath216 given by the solution to into @xmath227 , and @xmath228 to find @xmath229 write @xmath230 to see @xmath231 , which yields @xmath232 then @xmath233,\ ] ] and , with @xmath38 , @xmath234.\ ] ] this tells us to choose sign @xmath235 in the above formulas so that [ e810 ]    u^2= , + v^2= ,    which are the desired relations
.    now the compressible analogs of theorem [ t51 ] and corollaries [ c51][c52 ] follow by the same arguments we have employed in  [ s5 ] .
[ t81 ] for locally analytic initial data @xmath218 satisfying and @xmath236 , there exists a local solution in space - time of system  which satisfies both the compressible euler equations  and a propagating two - dimensional riemannian manifold isometrically immersed in @xmath5 .
[ c81 ] let @xmath218 be a local analytic solution in space - time of the compressible euler equations  satisfying the initial condition : @xmath236 , at a point @xmath237 .
then there is an evolving metric @xmath0 so that solution @xmath218 of the compressible euler equations  also defines a time evolving two - dimensional riemannian manifold @xmath238 immersed in @xmath5 with the gauss curvature @xmath22 determined .
theorems [ t72][t74 ] remain unchanged , except for the fact that metric @xmath3 for the wrinkled manifold now corresponds to the special incompressible solution of the compressible euler equations  with @xmath239 .
in particular , theorem [ t72 ] now gives the smooth solutions of the geometric image of the compressible euler equations
being approximated by the wrinkled solutions that correspond weak shear solutions of the _ incompressible _ euler equations .
the construction of the wrinkled solutions in this case is the same as done in  [ s7 ] and produces the solutions of the incompressible euler equations . a wrinkled solution of the compressible euler equations  is impossible since it would correspond to a vacuum via ( [ e84 ] ) .
we note that the results given in  7 and this section on the existence of _ wild _ weak shear solutions have been given in terms of a local system of normal coordinates , rather than the original cartesian coordinates @xmath240 .
this means for example that the _ wild _ weak shear solutions could in fact correspond to different local coordinates , _
e.g. _ , polar coordinates .
the choice of local cartesian coordinates used earlier is only a convenience and the euler equations written in polar coordinates would suffice .
motivated by our results for fluid dynamics , we now consider solutions in two - dimensional general continuum mechanics .
denote by @xmath241 the ( symmetric ) cauchy stress tensor , and assume that fields @xmath4 , @xmath242 , and @xmath216 are consistent with some specific constitutive equation for a body and satisfy the balances of mass and the linear momentum ( the satisfaction of the balance of angular momentum is automatic ) .
the equations for the balance of linear momentum in the spatial representation are [ e101 ]    _ 1(u^2-t_11)+_2(uv - t_12 ) = -_t(u ) , + _
1(uv - t_12)+_2(v^2-t_22 ) = -_t(v ) ,    and the balance of mass is @xmath243 or @xmath244 where @xmath245 is the density of the body in the reference configuration , and @xmath246 is the deformation gradient of the current configuration with respect to this reference .      denote the geometric dependent variables as [ lmn ] l = v^2-t_22,n = u^2- t_11,m= -uv + t_12 , so that the gauss equation becomes @xmath247    under the assumption that @xmath248 at least on some initial time interval , such a solution of continuum mechanics corresponds to a positive gauss curvature for the corresponding isometric embedding problem . using @xmath55 as time - dependent data in ( [ e31])([e32 ] ) and expressing these three equations in terms of the metric components @xmath73 and their derivatives by using ( [ e23])([e24 ] )
, we have a system of three partial differential equations for the components of the metric . using the required analog of corollary [ c81 ]
, then we have    let @xmath249 be a local analytic solution in space - time of system  at a point satisfying the condition : @xmath250 then there is an evolving metric @xmath0 so that the continuum mechanical solution defines a time evolving two - dimensional riemannian manifold @xmath6 immersed in @xmath251 .
we note the ability to solve for the evolving metric @xmath0 in the fluid case relied on the initial data for the off - diagonal term @xmath252 in the second fundamental form being non - vanishing . in the cases of the incompressible and compressible fluids ,
this term was simply @xmath93 .
however , for general continuum mechanics , @xmath253 that is , the expression for @xmath252 has an additional contribution .
theorems [ t72][t74 ] remain unchanged , except for the fact that metric @xmath3 for the wrinkled manifold now corresponds to special _ steady _ smooth solutions of the equations of two - dimensional general continuum mechanics .
we now identify these solutions .
we want to define solutions to the mechanical equations  in continuum mechanics from a smooth , degenerate isometric embedding problem for which @xmath254 , are satisfied
. then , making the association analogous to that used in lemma 4.2 , we have [ 10.8-a ]    uv - t_12 = 0 , + u^2 - t_11 = z(x_2,t ) , + v^2 - t_22 = 0 ,    where function @xmath255 is defined from the embedding , and    _ 1(uv - t_12 ) + _ 2 ( v^2 - t_22 ) = 0 , + _ 1 ( u^2 - t_11 ) + _ 2 ( uv - t_12 ) = 0    are satisfied by the mechanical fields that are being defined .
these would be a consistent set of the fields satisfying the balances of linear momentum and mass if the following constraints :    _ t ( u ) = _ t ( v ) = 0 , + _ 1 ( u ) + _ 2 ( v ) = -_t    are satisfied , _ i.e. _ , these geometric solutions define continuum mechanical solutions with steady momenta .
this is easily done by noting that the conditions imply @xmath256 with the solution @xmath257 where @xmath258 and @xmath259 are arbitrary time - independent functions of the spatial variables , which are required to be so chosen that @xmath260 .
we then define @xmath59 by integrating the pointwise ordinary differential equations :    u_t = - _ 1 u , + v_t = - _ 1 v.    with @xmath261 in hand , we define the stress components from to obtain a class of mechanical solutions in general continuum mechanics from the smooth degenerate isometric embedding problem .
it is to be noted that , for such continuum mechanical solutions to be realizable for a specific material , the mechanical fields as defined have to be shown to be consistent with a constitutive equation for the stress for that material .
alternatively , by seeking solutions to the system [ steady ]    _ 1(uv - t_12 ) + _ 2 ( v^2 - t_22 ) = 0 , + _ 1 ( u^2 - t_11 ) + _ 2 ( uv - t_12 ) = 0 , + _ 1 ( u ) + _ 2(v ) = -_t ,    we can define solutions to the balance of mass and the _ steady equations of balance of linear momentum_.    the solutions to system do not necessarily constitute exact , steady solutions of the balance of linear momentum , i.e. , @xmath262 and @xmath263 may not evaluate to @xmath264 from such motions , much in the spirit of quasi - static evolutions in solid mechanics and stokes flow in fluid mechanics .
as is understood , such solutions are typically interpreted in an asymptotic sense when the velocities are assumed to `` equilibrate '' on a much faster time - scale than the evolution of driving boundary conditions or forcing , i.e. , @xmath265 is assumed to be the `` slow '' time scale and the right - hand sides of the first two equations actually have the terms : @xmath266 and @xmath267 , with @xmath268 , so that we deal with the singular perturbation with a first - order approximation .
the goal in this section is to display the image of the isometric embedding problem in the smooth degenerate case in nonlinear elastodynamics for a neo - hookean material .
we choose the constitutive equation for the cauchy stress of the compressible material as @xmath269 where @xmath246 is the deformation gradient from a fixed reference configuration of the body .
the generic point on the fixed reference is denoted by @xmath270 , the motion as @xmath271 , and @xmath272 . for the sake of simplicity
, we ignore a multiplicative scalar function of the invariants of @xmath273 that would ensure that the stress at the reference configuration would vanish , even in this simplest frame - indifferent elastic constitutive assumption .
_ we now seek special solutions to the steady equations . after obtaining any such solution consistent with the posed steady problem
, we will further check which of them also constitute ( steady ) solutions to the equations of balance of linear momentum_.    thus , it needs to be demonstrated that the conditions : [ nh_gov_eq ]    uv - b_12 = 0 , + v^2 - b_22 = 0 , + u^2 - b_11 = z(x_2,t )    are satisfied , and the balance of mass is satisfied by _ a - posteriori _ setting @xmath274 with @xmath275 calculated from the obtained solution to ( [ nh_gov_eq ] ) , where @xmath245 is the mass density distribution in the reference configuration .
we always require that @xmath276 .
the components of @xmath273 are given as @xmath277 thus , the equations required to be satisfied by a motion consistent with these constraints would be [ elastic_syst ]    = + , + ( ) ^2 = ( ) ^2 + ( ) ^2 , + ( ) ^2 = ( ) ^2 + ( ) ^2 + z ( x_2(x , t),t ) .
we define @xmath278 , and @xmath279 for subsequent use .
if solutions exist to the above system , then they also satisfy the following conditions : the satisfaction of the first equation of ( [ elastic_syst ] ) by a solution satisfying the second and third is equivalent to @xmath280 which implies that @xmath281 , or else @xmath282 which is not acceptable by hypothesis .
thus , the steady , neo - hookean image of the degenerate isometric embedding problem may be stated
as the system for functions @xmath283 : [ es2 ]    ( ) ^2 = |x_2|^2 , + ( ) ^2 = |x_1|^2 - ( ) ^2 , + z(x_2 , t ) |x_2|^2 = - ( j(x_1 , x_2))^2 ,    with the caveat that solutions to satisfy @xmath284 in order to generate steady solutions consistent with , we need to choose only the @xmath285 sign .    when restricting attention to the class of solutions to the whole system that are scale - invariant in @xmath283 , _
i.e. _ , those solutions that remain solutions if the dependent and independent variables are scaled by the same constant , we note that the first two equations in ( [ es2 ] ) remain invariant under rescaling the independent and dependent variables by @xmath286 .
similarly , @xmath287 is independent of @xmath286 .
we now assume that the given function @xmath288_.    [ elast_example ] consider a shearing motion of the form : [ ex_ansatz ]    & x_1 = x_1 + f(x_2,t ) , + & x_2 = x_2 t.    then @xmath289 , @xmath290 , @xmath291 , and @xmath292 .
thus , @xmath293 , @xmath294 , and @xmath295 .
then the first equation of ( [ es2 ] ) is identically satisfied , the second requires @xmath296 and the third is identically satisfied by our choice of @xmath297 .
thus , the traveling waves of the form : @xmath298 define solutions to the steady , neo - hookean , degenerate isometric embedding problem , for each sign of @xmath299 in the second equation of ( [ ex_ansatz ] ) .    in summary
, we have shown that the elasticity analog of theorem [ t72 ] for fluids holds , with the identification that @xmath166 is one of the possible smooth degenerate ( _ i.e. , _
@xmath3 ) isometric embedding images of the neo - hookean motion and the corresponding _
@xmath181 are the weak euler shear solution inducing embeddings in the sense of definition [ weak_shear_embedding ] _ , where the shear would be with respect to the reference @xmath300configuration .    for a spatially uniform density field on the reference configuration ,
it is easy to check that the equations of balance of linear momentum for the neo - hookean constitutive assumption we are considering ( the first piola - kirchhoff stress is given by @xmath301 ) reduces to the linear , second order system : @xmath302 thus , for the assumed ansatz ( [ ex_ansatz ] ) , the image of the degenerate isometric embedding problem has produced special _
exact _ non - steady , solutions to the balance of linear momentum .
clearly , these are not steady solutions on the reference configuration in general .
thus , it is an interesting question to check whether any of these are _ exact _ steady solutions on the current configuration .    1 .
the solutions generated in example [ elast_example ] correspond to @xmath303 .
let us consider each in turn . to check steadiness , we consider an arbitrarily fixed point @xmath304 in space .
we have to show that the velocities of material points that occupy it at different times is the same , since the density is constant , @xmath305 , everywhere for this example .
let the image in the reference configuration of point @xmath306 at time @xmath265 be @xmath307 for @xmath308 , @xmath309 which implies @xmath310 + now @xmath311 then @xmath312 thus , we indeed have a steady spatial velocity and momenta field .
the way to physically understand it is the following : let @xmath313 be specified .
the points on line @xmath314 have material points sitting on them at time @xmath265 with referential @xmath315coordinate given by @xmath316 .
however , for such referential points , the horizontal velocity @xmath317 corresponds to @xmath318 hence , this is a very curious situation where the picture is completely steady on the current configuration ; however , there is unsteady wave propagation with piola - kirchhoff shear stress waves propagating on the reference configuration .
now we consider the other solution @xmath319 .
@xmath320 so that @xmath321 then we have unsteadiness in the current configuration as well as the reference . note
that this case corresponds to choosing the wrong sign in solutions to to correspond to solutions of .
however , if we now choose the ansatz for the motion to be @xmath322 then it can be checked that again @xmath323 are both solutions to system , but now only solution @xmath324 works ( for obvious reasons by following the previous argument ) : @xmath325 then again @xmath326 thus , we have a steady field on the current configuration , while , on the reference , stress / velocity / position waves are moving from the bottom to the top .    while simple , it is important to appreciate that the generated exact solution in this extremely simple example can produce smooth analogs ( with continuous deformation gradient ) of the traveling wave profile shown in fig .
[ f1 ]        with arbitrarily small positive values of @xmath327 and @xmath328 . when @xmath327 is comparable to @xmath328 , this is a sequence of phase boundaries separating domains ; when @xmath329 , then the regions spanned within width @xmath328 are slip zones or shear bands separating undeformed blocks .
the greater freedom that elastodynamics provides over elastostatics in producing microstructure can now be appreciated . while such deformations in elastostatics can only be produced with a multi - well energy function @xcite
, such microstructures ( for @xmath330 ) in elastodynamics occur in this case of the simplified neo - hookean material ( with the energy density @xmath331 up to a multiplying constant ) , and persist for all times . indeed , multiple low - energy states corresponding to large total strains and static microstructures . ] are observed facts , and their representation consists of an important physical ingredient of solid mechanics .
our intent here is to simply point out the greater freedom of displaying microstructure without any length scale in elastodynamics . by considering the case for ever smaller @xmath332 , it is also now perhaps easy to intuitively see why the wrinkled deformations of the isometric embedding problem may have a limiting status in elastodynamics , as the parameters @xmath327 and/or @xmath328 tend to zero .
furthermore , it is an interesting question whether our elastodynamic model ( [ es2 ] ) corresponding to the degenerate isometric embedding also displays solutions that represent static and dynamic terminating lines of phase boundaries ( generalized disclinations ) and shear bands ( dislocations ) , ingredients of nature that may help to further understand physically the nature of the wrinkled embeddings from geometry . finally , it is important to recognize and accept that microstructure in solids is not about infinite refinement and , as is well - known , such scale - free deformations , no matter how exciting from the analytical point of view , are notoriously difficult to deal with practically in modeling endeavors ( following nature s guide , in some sense ) .
thus , our considerations seem to point to the need for models that can represent both static , quasi - static ( _ i.e. _ , evolving microstructure in the absence of material inertia ) , as well as fully dynamic microstructures with in - built length scales and accounting for the dissipation due to microstructure evolution .
suppose that we have a smooth , time - dependent two - dimensional manifold immersed in @xmath251
. then @xmath73 and @xmath32 can be defined from the manifold as in  [ s2 ] .
consider the initial data to be available on @xmath218 denoted by @xmath333 pointwise . then define the following quantities @xmath334 through the ordinary differential equations : [ momenta ]    _ t m_u = _ 22 ^ 1l-2_12 ^ 1m+_11 ^ 1n , + _
t m_v = _ 22 ^ 2l-2_12 ^ 2m+_11 ^ 2n ,    with initial condition on @xmath335 specified as @xmath336 and on @xmath337 as @xmath338 . with fields @xmath334 available ,
solve @xmath339 with initial condition on @xmath216 to be @xmath245 .
then , using this @xmath216 field along with @xmath340 and @xmath55 as time - dependent data , solve the pointwise ordinary differential equations for fields @xmath59 given by ( [ momenta ] ) with @xmath334 replaced by @xmath341 , with the initial condition on @xmath59 being @xmath342 , respectively . clearly
, then we have @xmath343 however , then ( [ bal_mass ] ) implies that the constructed fields @xmath344 satisfy the balance of mass .
finally , defining the stress components from ( [ lmn ] ) by using @xmath55 and the constructed fields @xmath218 as the data , and noting ( [ momenta ] ) and ( [ e31 ] ) , we find that every smooth time evolving two - dimensional manifold immersed in @xmath251 defines a solution of the balance laws of two - dimensional general continuum mechanics .
it requires the imposition of further constraints to obtain the motions within this class that are consistent with the constitutive response of any specific material .    for a given constitutive relation ,
the cauchy stress is determined from the velocity field @xmath4 and density @xmath216 .
the case for ( incompressible and compressible ) inviscid fluids has been exposited , and the case of elasticity requires the deformation gradients to be calculated from the deformation field and to be inferred from velocity @xmath4 via time integration .
hence , since @xmath345 have been determined by our argument above , the ability to satisfy the constitutive relations requires that metric @xmath0 be consistent with these extra constraints .
let @xmath246 be the deformation gradient field of the mechanical body from a fixed stress - free elastic reference configuration .
let @xmath346 denote the velocity gradient field .
for simplicity , assume a local constitutive equation for the stress of the form : @xmath347 we also use notation @xmath348 and @xmath349 for the known functions for the christoffel symbols and the gauss curvature .
a time - dependent set of mechanical fields @xmath350 and a time - dependent metric field @xmath0 related by the constitutive assumption @xmath351 are consistent in the sense that the balance of linear momentum and balance of mass are satisfied , and the metric is isometrically embedded in @xmath251 at each instant of time , provided that the following system of constrained partial differential equations of evolution are satisfied : @xmath352 \partial_k f_{ij } v_k - ( { { \bf l}}_v)_{ik}f_{kj } = -\partial_t f_{ij},\\[1 mm ] \partial_{1}\big(\rho u^{2}-\hat t_{11}({{\bf l}}_v,\rho , { { \bf f } } , { { \bf g}})\big ) + \partial_{2}\big(\rho uv- \hat t_{12}({{\bf l}}_v,\rho , { { \bf f } } , { { \bf g}})\big )   =   -\partial_{t}(\rho u),\\[1 mm ] \partial_{1}\big(\rho uv - \hat t_{12}({{\bf l}}_v,\rho , { { \bf f } } , { { \bf g}})\big )   + \partial_{2}\big(\rho v^{2 } - \hat t_{22}({{\bf l}}_v,\rho , { { \bf f } } , { { \bf g}})\big)=-\partial _ { t}(\rho v ) , \end{cases}\ ] ] and @xmath353 \partial_{1 } m -\partial_{2}l    = \gamma_{22}^{2}({{\bf g}},\nabla { { \bf g } } ) l -2\gamma_{12}^{2}({{\bf g}},\nabla { { \bf g } } ) m + \gamma_{11}^{2}({{\bf g}},\nabla { { \bf g } } ) n,\\[1 mm ]   ln - m^2 = \kappa({{\bf g}},\nabla { { \bf g } } , \nabla^2 { { \bf g}}),\\[1 mm ] l \mathbf{e}^1 \otimes \mathbf{e}^1 + m \left(\mathbf{e}^1 \otimes \mathbf{e}^2 + \mathbf{e}^2 \otimes \mathbf{e}^1 \right ) + n \mathbf{e}^2 \otimes \mathbf{e}^2 = \mathbf{f } \left ( \rho \mathbf{u } \otimes \mathbf{u } - \mathbf{t}\left(\mathbf{l}_v , \rho , \mathbf{f } , \mathbf{g } \right ) \right ) , \end{cases}\ ] ] where @xmath354 is a tensor - valued function of its tensorial argument , and @xmath355 is the dual basis corresponding to the natural basis on the surface given by @xmath356 .
the last tensorial equation consists of three independent equations .
the equations in are to be considered as the constraints that determine the family of pairs of first and second fundamental forms of embedded manifolds consistent with a given mechanical state at any give time .
conversely _ , the evolution of the mechanical fields following the first four equations must be constrained to the _ manifold _ defined by in the state - space of spatial @xmath357 fields .
one may consider eliminating variables @xmath55 from the equations in to obtain three constraint equations for the three components of the metric field .
abstractly , one may think of eliminating all of the equations in and replacing @xmath0 in the mechanical set of the first four equations as a spatially non - local term in the mechanical fields representing a solution of .
notice that this remark applies as well to the steady problem of continuum mechanics ( where the right - hand sides of the third and fourth equations of are assumed to be @xmath264 ) and the time - dependence of boundary conditions ( or body forces that have been assumed to vanish here for simplicity ) drive the evolution of the mechanical problem .
we note that the earlier sections of this paper dealing with the equations of incompressible and compressible fluid dynamics , and neo - hookean elastodynamics are special cases of the above system where the constitutive equation for the cauchy stress are independent of metric @xmath0 .
we close by discussing some broad implications and possible extensions of the presented work .
since the dynamics of the wrinkled solutions shadowing the incompressible or compressible euler equations are completely reversible nash - kuiper solutions corresponding to the metric for the developable surface , the usual irreversible entropy admissibility criteria become useless as selection criteria .
this tells us that the only hope of selecting a unique physically meaningful solution would be chosen without recourse to time evolution , _
e.g. _ , artificial viscosity or energy minimization .
for example , an energy minimization which penalizes second derivatives of @xmath166 would prefer the affine initial data over data with folds , so that the _ wild _ initial data would be ruled out .
similarly , szkelyhidi @xcite , bardos - titi - wiedermann @xcite , bardos - titi @xcite , and bardos - lopes filho - niu - nussenzveig lopes - titi @xcite have shown that the viscosity criterion also eliminates the _ wild _ solutions for the euler equations .
the condition of incompressibility and the poisson equation for the pressure remain as @xmath359 & \partial_{11}^{2}(u^{2})+2\partial_{12}(uv)+\partial_{22}(v^{2 } ) = -\delta p. \label{e93}\end{aligned}\ ] ] if we review all the previous arguments made for the inviscid euler equations leading up to and including
[ s8 ] , we see all the conclusions we have made regarding the euler equations hold true for the navier - stokes equations modulo one crucial point . for the euler equations , metric @xmath3 provides the map from the steady shear , @xmath360 , a solution of the euler equations to the gauss - codazzi equations .
however , for a shear solution of the navier - stokes equations , the right - hand sides of must vanish .
therefore , instead of the fluid pre - image being a steady shear , the fluid pre - image must satisfy the diffusion equation : @xmath361 or @xmath362 thus , the pre - image is smooth so that , by , we have a smooth second fundamental form and a smooth embedding @xmath166 .
however , this contradicts the non-@xmath182 property of our nash - kuiper wrinkled solution .
thus , the only possibility is that , for the navier - stokes equations , there is no fluid pre - image of the nash - kuiper wrinkled solutions .
one may be tempted to discount the physical relevance of the nash - kuiper theorem outlined here and in the recent work in @xcite . on the other hand , as noted in chen - glimm @xcite and others ( _ cf_. @xcite ) , viscous fluid turbulence arises with the imposition of an external force , which overcomes the viscous dissipation .
hence , it may be that the addition of an external force would indeed bring us back to the nash - kuiper - gromov turbulence scenario .
furthermore , the dynamics of defect microstructure like dislocations , phase and grain boundaries , triple junctions , and point defects in crystalline solids furnish a compelling physical argument for accepting / developing physically rigorous and practically computable models that account for the representation of microstructure , necessarily then not of infinite refinement and with a modicum of uniqueness in the predicted evolution of their fields . in analogy with a crumpled piece of paper that does not produce infinitely fine terminated folds and neither unfolds itself back to its original flat state ,
perhaps the correct notion of admissibility is to move to augmented physical models of continuum defect dynamics involving extra kinematics representing the microscopic , and hence smoothed , dynamics of discontinuity surfaces , their terminating lines , and point singularities of the fields of the original macroscopic model ( like nonlinear elasticity and navier - stokes ) , while accounting for the energetics of these defects and the dissipation produced owing to their motion .
such partial differential equation - based augmentations of nonlinear elasticity theory have begun to emerge , _
e.g. _ , acharya - fressengeas @xcite , along with their interesting predictions of soliton - like dynamical behavior of nonsingular defects ( zhang et al .
@xcite ) and their collective behavior .
[ [ a - plausible - minimal - model - for - internally - stressed - elastic - materials ] ] a plausible minimal model for internally stressed elastic materials ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    the constrained evolution system ( [ constrained_system ] ) , written on the deforming configuration of an elastic body being tracked in lagrangian fashion , appears to pose an interesting model for internally stressed elastic bodies whose constitutive response in terms of the deformation gradient and a metric representing a stress - free state is known : @xmath363 \partial_{1}\hat t_{11}({{\bf f } } , { { \bf g}})+\partial_{2}\hat t_{12}({{\bf f } } , { { \bf g}})=\rho\ , d_t u,\\[1 mm ] \partial_{1}\hat t_{12}({{\bf f } } , { { \bf g}})+\partial_{2}\hat t_{22}({{\bf f } } , { { \bf g}})= \rho \,d_t v , \end{cases}\ ] ] and @xmath364 \partial_{1}m -\partial_{2}l = \gamma _ { 22}^{2}({{\bf g}},\nabla { { \bf g } } ) l -2\gamma_{12}^{2}({{\bf g}},\nabla { { \bf g } } ) m + \gamma_{11}^{2}({{\bf g}},\nabla { { \bf g } } ) n,\\[1 mm ]   ln - m^2 = \kappa({{\bf g}},\nabla { { \bf g } } , \nabla^2 { { \bf g}}),\\[1 mm ]   l \mathbf{e}^1 \otimes \mathbf{e}^1   + m \left(\mathbf{e}^1 \otimes \mathbf{e}^2 + \mathbf{e}^2 \otimes \mathbf{e}^1 \right )   +
n \mathbf{e}^2 \otimes \mathbf{e}^2   = \mathbf{\hat{f } } \left (   \mathbf{t}\left(\mathbf{f } , \mathbf{g } \right ) \right ) , \end{cases}\ ] ] where @xmath365 , represent the spatial derivatives on the current configuration , @xmath366 represents the material time derivative operator , @xmath367 is a tensor - valued function of its tensorial argument , and @xmath355 is the dual basis corresponding to the natural basis on the surface given by @xmath356 .
the last tensorial equation consists of three independent equations .
first of all , we note that , on dimensional grounds , the identification of the second fundamental form with mechanical objects with physical dimensions of stress implies from the codazzi equations that the metric is physically dimensionless and the christoffel symbols have dimensions of reciprocal lengths .
thus , the metric may be interpreted as describing strain .
of course , this identification also implies that a material parameter with physical units of ( @xmath368 is required to make the gauss curvature equation dimensionally consistent , while introducing a length - scale into the traditional elastic problem of internal stress .    next , we consider the equations in ( [ elast_constrained_system - a ] ) for given @xmath246 .
considering , for the moment , the situation when the constitutive equation is independent of @xmath0 , this becomes a question of determining the metric , given the second fundamental form , such that an embedding exists in @xmath251 , which is the opposite of the isometric embedding problem that may be interpreted as the question of determining the second fundamental form , given a metric @xmath0 . at any rate , as some of our results show , this problem has a solution in many instances .
it is also perhaps reasonable to expect that the situation does not change drastically even when the constitutive equation depends on @xmath0 , and is not degenerate in the sense that , for almost all @xmath246 , there are many solutions for @xmath0 ( this would necessitate the evolution equations for @xmath0 ) .
based on this premise , the requirement of an embedding of @xmath0 in @xmath251 assumes a physical status replacing a separate constitutive equation for the evolution of @xmath0 that would be required for the mechanical problem otherwise .
moreover , much like in systems displaying relaxation oscillations , the geometry and stability of states on the constraint manifold can , on occasion , lead to interesting dynamical behavior of the mechanical fields .
the above model appears to have the possibility of being considered as a minimal model for the statics and dynamics of soft and biological materials ( _ e.g. _ , efrati et al @xcite , jin et al @xcite , and ambrosi et al .
@xcite ) depending only on the knowledge of the elastic response of the material and the material constant needed to define the gauss curvature equation .    *
acknowledgments*. a. acharya acknowledges the support of the rosi and max varon visiting professorship at the weizmann institute of science , rehovot , israel , and was also supported in part by grants nsf - cmmi-1435624 , nsf - dms-1434734 , and aro w911nf-15 - 1 - 0239 .
chen s research was supported in part by the uk epsrc science and innovation award to the oxford centre for nonlinear pde ( ep / e035027/1 ) , the uk epsrc award to the epsrc centre for doctoral training in pdes ( ep / l015811/1 ) , and the royal society  wolfson research merit award ( uk ) .
m. slemrod was supported in part by simons collaborative research grant 232531 .
d. wang was supported in part by nsf grant dms-1312800 .
a. acharya and c. fressengeas , continuum mechanics of the interaction of phase boundaries and dislocations in solids . in : _
differential geometry and continuum mechanics _
, pp . 123165 , eds .
chen et al .
( eds . ) , springer proceedings in mathematics and statistics , * 137 * , 2015 .
d. ambrosi , g.  a. ateshian , e.  m. arruda , s.  c. cowin , j. dumais , a. goriely , g.  a. holzapfel , j.  d. humphrey , r. kemkemer , e. kuhl , j.  e. olberding , l.  a. taber , and k. garikipati , perspectives on biological growth and remodeling .
_ j. mech .
* 59 * ( 2011 ) , 863883 .          c.  w. bardos , m.  c. lopes filho , d.  niu , h.  j. nussenzveig lopes , and e.  s. titi , stability of two - dimensional viscous incompressible flows under three - dimensional perturbations and inviscid symmetry breaking .
_ siam j. math .
_ * 45 * ( 2013 ) , 18711885 .
j.  f. borisov , on the question of parallel displacement on a smooth surface and the connection of space forms of smooth surfaces with their intrinsic geometries .
_ vestnik leningrad .
_ * 15 * , 19 ( 1960 ) , 127129 .
s. conti , c. de lellis , and l. szkelyhidi , the @xmath371-principle and rigidity for @xmath7 isometric embeddings . in : _
nonlinear partial differential equations _ , pp
. 83116 , abel symp .
* 7 * , springer : heidelberg , 2012 .
d. deturck and d. yang , local existence of smooth metrics with prescribed curvature . in : _
nonlinear problems in geometry ( mobile , ala . , 1985 )
_ , pp . 3743 , _ contemp .
_ * 51 * , ams : providence , ri , 1986 .
x. zhang , a. acharya , n. j. walkington , and j. bielak , a single theory for some quasi - static , supersonic , atomic , and tectonic scale applications of dislocations .
* 84 * ( 2015 ) , 145195 .