<ph f="cmex"> </ph><ph f="cmbx">On the interaction between quasilinear elastodynamics and the Navier-Stokes equations</ph>
E-mail address : coutand@math.ucdavis.edu E-mail address : shkoller@math.ucdavis.edu Daniel Coutand and Steve Shkoller Department of Mathematics University of California at Davis Davis, CA 95616
-
Abstract.
The interaction between a viscous fluid and an elastic solid is modeled by a system of parabolic and hyperbolic equations, coupled to one another along the moving material interface through the continuity of the velocity and traction vectors. We prove the existence and uniqueness (locally in time) of strong solutions in Sobolev spaces for quasilinear elastodynamics coupled to the incompressible Navier-Stokes equations. Unlike our approach in [5] for the case of linear elastodynamics, we cannot employ a fixed-point argument on the nonlinear system itself, and are instead forced to regularize it by a particular parabolic artificial viscosity term. We proceed to show that with this specific regularization, we obtain a time interval of existence which is independent of the artificial viscosity; together with a priori estimates, we identify the global solution (in both phases), as well as the interface motion, as a weak limit in srong norms of our sequence of regularized problems.
1 Introduction
We establish the existence and uniqueness in Sobolev spaces of strong solutions to the unsteady fluid-structure interaction problem consisting of a nonlinear large-displacement elastic solid coupled to a viscous incompressible Newtonian fluid. The fluid motion is governed by the incompressible Navier-Stokes equations, while the solid, which can be either compressible or incompressible, is modeled by the celebrated St. Venant-Kirchhoff constitutive law (although our method can be applied to more general quasilinear hyperelastic models).
The first fluid-solid interaction problems solved were for the case of a rigid body inside of a viscous flow in a bounded domain (see [7] , [12] , [4] , [13] ), and the case of a rigid body inside of a viscous flow in an infinite domain ([22] , [20] , [15] ). Later, the elastic body was modeled with the restriction of either a finite number of modes ([8] ) or with a hyperviscous type law for the solid ([2] , [10] ), essentially by the same type of Eulerian global variational methods developed in [7] . For the steady-state problem, which is elliptic in both phases, [11] solved the case of solid modeled as a St. Venant-Kirchhoff material. In [18] , an Eulerian approach was used for the case in which the solid is a visco-hyperelastic material, which is a regularization of a hyperbolic model of solid deformation.
With the exception of our recent well-posedness result for the case of a linear elastic solid in [5] , there are no known existence results for fluid-structure interaction when the solid is modeled by a standard second-order hyperbolic equation. This may be attributed to the difficulties associated with coupling a parabolic PDE for the fluid with a hyperbolic PDE for the solid through continuity of the velocity and traction vectors across the moving material interface. As we explained in [5] , an iteration scheme between fluid and solid phases fails to converge due to a regularity loss induced by the hyperbolic phase (this divergent behavior has been computationally noted as well in [14] ), and so we developed a method comprised of the following new ideas: first, a functional framework which scales in a hyperbolic fashion for both the fluid and solid phases. This scaling leads to additional compatibility conditions in the fluid phase (when compared to the use of the classical parabolic framework), and is absolutely crucial for obtaining consistent energy estimates.
Second, we developed a regularity theory founded upon central trace estimates for the velocity vector restricted to the interface, rather than traditional interior regularity arguments which do not work for our problem. Third, we were forced to bypass the use of the frozen (or constant) coefficient basic linear problem, which requires estimates on one more time derivative of the pressure function than the initial data allows, and created a new method wherein the solution was found as a limit of a sequence of penalized problems set in the Lagrangian framework.
The penalization scheme approximates the divergence-free constraint, whereas the Lagrangian framework alleviates the difficulties associated with the lack of a priori estimates in the solid phase for the frozen coefficient problem; this method indeed differs significantly from the classical methods used in fluid-fluid interface problems (see for instance [21] , [1] ).
The fundamental difficulty in extending our result to the case of nonlinear elasticity is the absence of any method of analysis for quasilinear elastodynamics which is compatible with the general scheme of [5] , involving a global Lagrangian variational formulation and the use of difference quotients to track the regularity of interface data. We remind the reader that unlike the analysis of elastostatic motion, direct inverse function theorem arguments cannot be applied directly to the case of quasilinear elastodynamics due to the fact that the perturbation term arising from the nonlinear operator is not an element of the appropriate function space for optimal regularity. Alternatively, one might attempt a fixed-point approach, wherein a portion of the nonlinear elasticity operator is viewed as a forcing function coming from a given velocity
, and then try to solve a linear problem for an unknown
. The difficulty in this approach stems from the fact that one has to find exact time derivatives of elastic energies for the forcing term associated to the elasticity operator, which is complicated by the inner-product of a term involving
and a term involving
. This difficulty is overcome in [6] , by a clever and essential use of the Dirichlet boundary condition in order to reformulate the problem in a non-standard way. As it turns out, the variety of known methods that have been used in the well-studied area of quasilinear elasticity, such as those in [6] and [17] for the Dirichlet boundary condition, or in [19] and [9] for the Neumann boundary condition, require a priori knowledge of the boundary data regularity, and are hence instrinsically incompatible with fluid-structure interaction analysis (in fact, the methods devised for Dirichlet conditions do not work for Neumann conditions and vice versa). Indeed, of these various methodologies, only [6] and [17] use a variational approach, the others employing either semi-group techniques as in the early work of [16] in the full space, or the technical paradifferential calculus in [19] for the two dimensional Neumann case.
In this paper, we develop a new method for quasilinear elastodynamics, variational in nature, that is compatible with our method in [5] ; we proceed in two steps.
First, we add a specific artificial viscosity to the solid phase which regularizes the system, thus converting our hyperbolic PDE into a parabolic one, transforming the fluid-structure interaction into a fluid-fluid interface-type problem for which existence and uniqueness of solutions is already known on a time interval that a priori shrinks to zero as the artificial viscosity
tends to zero. Second, and this is where the primary difficulty rests, we prove that our specific choice of parabolic smoothing renders the time interval (on which a unique solution exists) independent of
; furthermore, our a priori estimates allow us to construct a solution by weak convergence in strong norms. We note that the use of higher-order operators in the artificial viscosity term, while providing the necessary a priori control of the regularity of the moving interface, would not yield
-independent estimates which are essential here. Also, as our parabolic reguralization method is not specialized to any particular boundary condition, it thus provides a unified approach to the classical problem of quasilinear elastodynamics when the solid is not coupled to a fluid.
We now proceed to the formulation of our problem. Let
denote an open, bounded, connected and smooth domain with smooth boundary
which represents the fluid container in which both the solid and fluid move. Let
denote the closure of an open and bounded subset representing the solid body at each instant of time
with
denoting the fluid domain at each
. Note that in our analysis
is not necessarily connected, which allows us to handle the case of several elastic bodies moving in the fluid.
Remark 1.
If a function
is defined on all of
, we will denote
and
. This allows us to indicate from which phase the traces on
of various discontinuous terms arise, and also to specify functions that are associated with the fluid and solid phases.
For each
, we wish to find the location of these domains inside
, the divergence-free velocity field
of the fluid, the fluid pressure function
on
, the fluid Lagrangian volume-preserving configuration
, and the elastic Lagrangian configuration field
such that
, where
and
solves the Navier-Stokes equations in
:
| |
with
|
(1.1)
|
and
solves the elasticity equations on
with
and where the equations are coupled together by the continuity of the normal component of stress along the material interface
expressed in the Lagrangian representation on
as
and the continuity of particle displacement fields along
together with the initial conditions
, and
and the Dirichlet (no-slip) condition on the boundary
of the container
, where
is the kinematic viscosity of the fluid,
and
denote the Lamé constants of the elastic material,
is the outward unit normal to
and
is twice the rate of deformation tensor of
, given in coordinates by
. All Latin indices run through
, the Einstein summation convention is employed, and indices after commas denote partial derivatives.
We now briefly outline the proof. As the solid and fluid phases are naturally expressed in the Lagrangian and Eulerian framework, respectively, we begin by transforming the fluid phase into Lagrangian coordinates, leading us to the system of equations ( 3.2 ) and, as in [5] , we work in an hyperbolic framework in order to accomodate the dual nature of the problem (parabolic in the fluid and hyperbolic in the solid).
In order to solve ( 3.2 ), in Section 7 , we first add a particular form of artificial viscosity to the quasilinear hyperbolic equation in the solid, transforming the hyperbolic phase into a parabolic one; specifically, we add the term
, where
denotes the linearized (about the identity) elasticity operator and
is the material velocity.
We hence obtain an interface problem that is parabolic in nature in both phases, and can be thought as a fluid-fluid parabolic interface problem for which well-posedness is classical, save for the fact that both phases are required to scale in an hyperbolic fashion. The time interval of existence
for this parabolic system a priori shrinks to zero as
.
In Section 8 , we establish
-independent estimates on the solutions
to the regularized parabolic problem on the time interval
by identifying exact time derivatives of elastic energies, and establish regularity of the interface. A direct fixed-point approach for ( 3.2 ) does not appear to yield these exact time derivatives for the elastic energy, whereas the the regularized problem ( 7.1 ) does indeed lead to them. An essential key for obtaining estimates independent of
inside the solid is Lemma 6.1 . Whereas the trace estimates could be carried with other choices of artificial viscosity, we absolutely need the special choice made in our analysis in order to recover the regularity inside the solid independently of
. In particular, a different choice of a regularizing operator either of the same order such as
or of higher order such as
or
would not provide
-independent estimates.
In Section 9 , we then explain how our estimates allow the construction of solutions
on a time interval independent of
, still with energy estimates independent of
. The existence of a solution to ( 3.2 ) then follows by weak convergence as
.
Uniqueness is established in Section 11 in the same functional framework used for existence.
As our method seemingly requires more regularity on the initial data in the solid than it should, due to the artificial viscosity in the compatibility conditions, we explain in Section 12 how this extra regularity can be removed, thus leading to the result with optimal regularity.
Section 13 is dedicated to the case where the incompressibility constraint is added to the solid. The additional difficulty with respect to the compressible case comes from the fact that we control the velocity uniformly in
in function spaces which possess less regularity than in the fluid, whereas the pressure is controlled uniformly in the same regularity spaces in both phases. Also, we cannot use Lemma 6.1 in the most optimal form for the regularity of the pressure in the solid phase.
2 Notational simplification
Although a fluid with a Neumann (free-slip) boundary condition indeed obeys the constitutive law ( 1.1 ), we will replace for notational convenience ( 1.1 ) with
|
(2.1)
|
this amounts to replacing the energy
by
, which is not a problem mathematically due to the well-known Korn inequality. Henceforth, we shall take ( 2.1 ) as the fluid constitutive law.
3 Lagrangian formulation of the problem
In regards to the forcing functions, we shall use the convention of denoting both the fluid forcing
and the solid forcing
by the same letter
. Since
has to be defined in
(because of the composition with
), and
must be defined in
, we will assume that the forcing
is defined over the entire domain
.
Let
|
(3.1)
|
where
denotes the matrix of partial derivatives of
.
Clearly, the matrix
depends on
and we shall sometimes use the notation
to denote the formula ( 3.1 ).
Let
denote the Lagrangian or material velocity field,
is the Lagrangian pressure function (in the fluid), and
is the fluid forcing function in the material frame. Then, as long as no collisions occur between the solids (if there are initially more than one) or between a solid and
, the problem can be reformulated as
|
(3.2a)
|
|
(3.2b)
|
|
(3.2c)
|
|
(3.2d)
|
|
(3.2e)
|
|
(3.2f )
|
|
(3.2g)
|
|
(3.2h)
|
where
denotes the outward-pointing unit normal to
(pointing into the solid phase), and
Throughout the paper, all Greek indices run through
and all Latin indices run through
. Note that the continuity of the velocity along the interface is satisfied in the sense of traces on
by condition ( 3.2f ), whereas the continuity of the normal stress along the interface is represented by ( 3.2e ).
4 Notation and conventions
We begin by specifying our notation for certain vector and matrix operations.
-
We write the Euclidean inner-product between two vectors
and
as
, so that
.
-
The transpose of a matrix
will be denoted by
, i.e.,
.
-
We write the product of a matrix
and a vector
as
, i.e,
.
-
The product of two matrices
and
will be denoted by
, i.e.,
.
For
and
, we set
| |
with
defined with
replacing
.
In order to specify the initial data for the weak formulation, we introduce the space
which is endowed with the
scalar product.
The space of velocities,
, where the solution to ( 3.2 ) lives, is defined as the following separable Hilbert space:
|
(4.1)
|
endowed with its natural Hilbert norm
| |
| |
We also need the space
| |
endowed with its natural Hilbert norm
| |
We shall also need
-in-time control of certain norms of the velocity, which necessitates the use of the following closed subspace of
:
| |
| |
endowed with the following norm
| |
Finally, we will also make use of the space
| |
| |
endowed with its natural norm
| |
Remark 2.
Note that our functional framework does not make use of the third time derivative of the pressure
, even though we do use the third time derivative of velocity
; this functional framework is necessitated by the fact that the Dirichlet boundary condition together with the limited regularity of
does not allow us to obtain
with the appropriate regularity. Note also that we have added the
-in-time control of
in the definition of
mostly for a more convenient way to prove our theorems, rather than absolute necessity.
Throughout the paper, we shall use
to denote a generic constant, which may possibly depend on the coefficients
,
,
, or on the initial geometry given by
and
(such as a Sobolev constant or an elliptic constant). For the sake of notational convenience, we will also write
for
.
5 A first theorem
We now state our first theorem. We impose greater regularity requirements on the initial data than is optimal so as to avoid technical difficulties associated with a particular type of initial data regularization that would otherwise be necessitated.
We consider the case of optimal regularity on the initial data in Theorem 12.1 .
Theorem 5.1.
Let
be a bounded domain of class
, and let
be an open set (with a finite number
of connected components) of class
such that
and such that the distance between two distinct connected components of
(if there are multiple solid components) is greater than zero. Let us denote
. Let
,
,
be given. Let
|
(5.1a)
|
|
(5.1b)
|
Assume that the initial data satisfies
as well as the compatibility conditions
|
(5.2a)
|
| |
|
(5.2b)
|
| |
|
(5.2c)
|
| |
|
(5.2d)
|
where the time derivatives appearing in these equations and in the following ones are computed from any
satisfying
,
(
), and from any
satisfying
(
), the quantities
and
being defined as follows. First,
is defined by
|
(5.3a)
|
|
(5.3b)
|
|
(5.3c)
|
and
by
|
(5.4a)
|
|
(5.4b)
|
Note that
since
and
on
,
on
. We also have
defined by
| |
|
(5.5a)
|
| |
|
(5.5b)
|
|
(5.5c)
|
and
by
|
(5.6a)
|
|
(5.6b)
|
Finally,
is defined by
| |
|
(5.7a)
|
| |
|
(5.7b)
|
| |
|
(5.7c)
|
Then there exists
depending on
,
, and
, such that there exists a unique solution
of the problem ( 3.2 ). Furthermore,
.
Remark 3.
The remarks appearing in [
5]
at the end of Section 5 concerning the compatibility conditions and forcing functions for the linear elasticity case still hold in this setting with the necessary adjustments. In particular, we do not need the forcing functions to have the same regularity in both phases.
6 Preliminary result
In the remainder of the paper, we set
In our limit process as the artificial viscosity tends to zero, we will make use in a crucial way of the basic following result:
Lemma 6.1.
Let
and
be such that
and
|
(6.1)
|
Then, independently of
,
-
Proof.
Since
, let
be such that
If
, then the statement of the Lemma is satisfied. Now, let us assume that
. Let
be arbitrary. From ( 6.1 ), we infer that
| |
| |
From the definition of
we then infer that for any
,
which after division by
gives at the limit
:
which concludes the proof of the Lemma. □
Remark 4.
It should be clear that Lemma 6.1 applies to a more general class of linear operators than
.
7 The smoothed problem and its basic linear problem
As we described in the introduction, we cannot find an appropriate linear problem whose fixed-point provides a solution to ( 3.2 ). We are thus lead to introduce introduce the following (parabolic) regularization of ( 3.2 ), with the artificial viscosity
:
|
(7.1a)
|
|
(7.1b)
|
|
(7.1c)
|
|
(7.1d)
|
|
(7.1e)
|
|
(7.1f )
|
where
|
(7.2a)
|
|
(7.2b)
|
and
|
(7.3a)
|
|
(7.3b)
|
Solutions to ( 3.2 ) will be obtained as the limit as
of solutions to ( 7.1 ).
Suppose that
is given. Let
and let
be the quantity associated with
through ( 3.1 ).
We are concerned with the following time-dependent linear problem, whose fixed-point
provides a solution to ( 7.1 ):
|
(7.4a)
|
|
(7.4b)
|
|
(7.4c)
|
| |
|
(7.4d)
|
|
(7.4e)
|
|
(7.4f )
|
Remark 5.
The two forcing functions ( 7.3a ) are introduced for compatibility conditions at time
, allowing the solution of ( 7.4 ) to satisfy
and even to satisfy the same initial conditions as solutions to ( 3.2 ) would.
In the following, for the sake of notational convenience, we will denote by
a generic smooth function depending only on
(with the convention that
), by
a generic smooth function depending only on
and by
a generic smooth function depending only on
. Then, let
be defined by
|
(7.5a)
|
|
(7.5b)
|
where the time derivatives are computed with any
,
,
,
.
Let us now define
|
(7.6a)
|
|
(7.6b)
|
|
(7.6c)
|
By proceding as in [5] , we can establish the existence of a fixed-point to the system ( 7.1 ). This follows the lines of [5] by first approximating by a penalty scheme the divergence-free constraint in the fluid in our Lagrangian setting, and by performing a regularity analysis of the solution of ( 7.4 ) allowing the use of the Tychonoff fixed-point theorem. Given the estimates obtained in [5] , no new difficulty arise, since the parabolic artificial viscosity in the solid controls the forcing coming from the quasilinear part on a short time which is a priori shrinking to zero, and for this reason the proof is omitted here.
This leads us to the following
Lemma 7.1.
There exists
depending a priori on
and on a given expression of the type
, so that there exists a unique solution
of the regularized problem ( 7.1 ). Moreover,
.
In the next section we will study the limit of these solutions of the smoothed problems as
(this being problematic since the solutions to these regularized problems are a priori defined on a time interval shrinking to zero as
).
Moreover, the following variational equations (for
are satisfied for any test function
:
| |
| |
| |
|
(7.7)
|
together with the initial conditions
,
,
and
,
,
. Moreover for the third time differentiated problem in time, we also have that a.e. in
,
| |
| |
| |
| |
| |
| |
|
(7.8)
|
where, recall that
does not depend on the artificial viscosity
. The following result will be fundamental to our proof that the time interval of existence of solutions to ( 7.1 ) is in fact
-independent.
Lemma 7.2.
The mapping
is continuous on
.
-
Proof.
The continuity with respect to
of the terms of the type
is obvious, and since
(due to our artificial viscosity), so is the continuity of
. The only terms that remain are
and
.
In order to treat them, we will invoke the fact that due to our artificial viscosity in the solid, we in fact have
, which provides
. For the second time derivative of the pressure, we notice that from the variational form, true almost everywhere on
for any
,
| |
| |
| |
and the Lagrange multiplier Lemma 13 of [5] associated to the continuity results previously established, we have the continuity of
on
.
We now explain briefly why such a control on the fourth time derivative of
holds, and is possible only with the addition of the artificial viscosity in the solid. In particular, this norm cannot be controlled as
, which is not crucial for our purposes in any case. In order to understand the idea, we return to the level of the setting of the fixed-point argument, where we assume that
in an appropriate convex set of
is given, and search for a solution
of ( 7.4 ) by a Galerkin approximation on a penalized problem (for the pressure), in a way similar to [5] . The penalization parameter
is given, and we denote
, where
is solution of the Galerkin approximation at rank
, and where
is computed from
associated to the given
. Our interest will be with the first problem that appears in our methodology in [5] ; namely, the highest order time-differentiated problem is multiplied by
(which is permitted since it belongs to the appropriate finite dimensional space), and then integrate from
to
. We obtain
| |
| |
| |
| |
| |
leading us for a time small enough depending on the artificial viscosity
(but not on
and
) to an inequality of the type,
| |
| |
where
depends a priori on
. By proceding in a way inspired by our methodology in Section 9 of [5] , we can then prove that we have control, independently of
, on the first three norms. Taking the limit first as
and then as
, indeed provides us with
as announced. □
We note that this latter regularity property in the solid is only possible with the artificial viscosity
.
8 Estimate for the solutions of ( 7.4 ) independently of
In this section, we will denote
and denote the corresponding quantities
by
. In what follows,
is a given positive number to be made precise later when it will be chosen to be sufficiently small.
8.1 Energy estimate for
independently of
.
We are now going to use the regularity result
in the energy inequality ( 7.8 ) (which was established independently of the artificial viscosity), this time by interpolating and using the energy properties of the nonlinear elasticity operator, in order to get an estimate independent of the artificial viscosity.
Step 1. Let
. An integration by parts in time shows that
| |
| |
and thus with the properties of the Bochner integral in
,
| |
we deduce
| |
|
(8.1)
|
Step 2. Let
. Similarly,
| |
| |
By the same type of argument as for the previous step, we then get
| |
| |
| |
| |
| |
| |
and thus,
|
(8.2)
|
Step 3. Let
. By an integration by parts in time,
| |
| |
Similarly as before, we get
| |
| |
| |
| |
| |
| |
| |
and therefore
|
(8.3)
|
Step 4. Let
. By symmetry of
, we notice that
| |
| |
and thus,
| |
| |
| |
| |
| |
|
(8.4)
|
Now, Let
. By expanding the integrand with respect to the time derivatives and using the relation in
:
and estimates similar as in the previous steps, we find that
|
(8.5)
|
Step 5. By using ( 8.4 ) and ( 8.5 ) we find that
|
(8.6)
|
Step 6. By proceding in a way similar to [5] , except that we replace the constants
appearing there by appropriate powers of
, we find that the integrals set in the fluid domain are bounded by
| |
Step 7. Thus, from ( 7.8 ), and Steps 1-6, we then get on
:
| |
|
(8.7)
|
Step 8. The estimate of
in
, independently of
, will require some adjustments with respect to the methodology of [5] . To this end, we notice that we can apply a Lagrange multiplier Lemma similar to Lemma 13 of [5] , but corresponding to the case
, to the variational form true on
: for all
,
| |
| |
| |
| |
which provides for any
:
| |
| |
| |
By using ( 8.7 ) for the first four terms of the right-hand side of this inequality and remembering that the
norm of
is part of the norm
for the next two terms of this inequality, we get
| |
|
(8.8)
|
8.2 Estimate on
and
.
From the previous estimates, and the arguments that we will see hereafter for the case of
, we have
| |
| |
|
(8.9)
|
Similarly, we infer from ( 8.9 ) that
| |
| |
|
(8.10)
|
8.3 Estimate on
We will denote
,
, and
. We denote by
an
diffeomorphism from
into a neighborhood
of a point
such that
,
,
, with
. We consider a cut-off function
compactly supported in
, and equal to
in
.
With the use of test functions
(which is in
) in ( 7.7 ) for
, and by denoting
,
,
, we get after integrating by parts appropriately and letting
,
| |
| |
| |
| |
| |
| |
|
(8.11)
|
where
,
,
.
Remark 6.
Note that this limit process as
for the nonlinear elastic energy is possible because
(
) due to our artificial viscosity in the solid. Whereas we could also use difference quotients, it appears that the product rules are less cumbersome with the use of horizontal derivatives instead, which is permitted since we already know at this stage the regularity of
and
. Also, the limits on the right-hand side of ( 8.11 ) do not present any difficulty, given the regularity of the forcing functions and three integrations by parts with respect to horizontal variables.
Remark 7.
Since
is compactly supported in
, the integrals set on
,
,
do not depend on the extension that we chose for
,
or
, and simply represent a more convenient way to write these integrals.
Step 1. Let
. By using the
regularity of the coefficients
,
| |
| |
| |
| |
|
(8.12)
|
Step 2. Let
denoting the set of permutations of
,
| |
| |
| |
| |
| |
| |
| |
| |
From the regularity of
and the
regularity of
, we then infer
| |
| |
| |
with
|
(8.13)
|
By integrating by parts in time, we deduce
| |
| |
| |
| |
| |
| |
which implies in turn
| |
|
(8.14)
|
With
(
) denoting the canonical vectors of
, let
where
and
are arbitrarily fixed in
. We then have
with
| |
| |
We first notice that
| |
| |
Next, by writing
and
respectively in
and
, we obtain
| |
| |
| |
Next, we see that
| |
| |
and by the same type of arguments as for
,
| |
implying
|
(8.15)
|
Now, from the definition of a solution to the smoothed problem ( 7.1 ),
| |
| |
which implies with ( 8.15 )
| |
| |
Since this inequality also holds for any
, Lemma 6.1 provides
| |
| |
which with the estimate on
from the previous subsection leads to
| |
|
(8.16)
|
Step 3. From the estimates on
, and similar estimates that we could get in the fluid as in [5] , but this time by replacing
by appropriate powers of
, we then deduce that for all
,
| |
| |
| |
By the trace theorem, we then get
| |
| |
where
. By a finite covering argument, we then get
| |
|
(8.17)
|
From the estimate ( 8.10 ) on
and the trace estimate ( 8.17 ), we infer in a way similar to [5] by elliptic regularity arguments that
| |
| |
|
(8.18)
|
Similarly, from ( 8.16 ), and the trace estimate ( 8.17 ), elliptic regularity provides
| |
|
(8.19)
|
9 Time of existence independent of
From ( 8.7 ), ( 8.8 ), ( 8.9 ), ( 8.10 ), ( 8.18 ) and ( 8.19 ), we then have for any
,
| |
| |
The subscripts
in
,
,
mean that we no longer consider generic constants from now on.
Now, let
be such that
. For
small enough, for any
we have
|
(9.1)
|
where
. For conciseness, we will denote
and
.
Now for
fixed, let
, so that
Now let
, which does not depend on
, and let
.
From now on, we assume that
. We then have
which implies that
has three real roots
,
,
, with
. From the product
and
, we infer that
. From ( 9.1 ) and the continuity of
(established in Lemma 7.2 ) we then infer since
that we have
|
(9.2)
|
This implies that
,
,
,
,
,
,
,
, with a bound that depends only on the right-hand side of ( 9.2 ). The compatibility conditions for the smoothed problem ( 7.1 ) at
are also satisfied by definition of a solution, which means that we do not have any new term of the type of
,
or
associated to
to add to the already existing forcing terms coming from
.
We can thus build a solution of the smoothed problem ( 7.1 ) defined on
,
depending solely on the right-hand side of ( 9.2 ), that we will still denote
.
It is then readily seen that
and is a solution of the approximated problem ( 7.1 ) on
. If
, we have our solution defined on the
independent time interval
, with the
independent estimate ( 9.2 ). Otherwise, if
, we can also assume that
, which implies, in the same fashion as we got ( 9.1 ),
|
(9.3)
|
This implies in turn that
,
,
are in the same spaces as their respective counterparts at time
, with the same bound as well, since we could from ( 9.3 ) repeat the same argument leading to ( 9.2 ), this time on
. The compatibility conditions at
being also automatically satisfied, we can thus build a solution of the approximated problem ( 7.1 ) defined on
, the time of existence being the same as starting from
from the similarity of the bound that we obtain on
and their respective counterparts at time
. We will still denote this solution
. It is then readily seen that
and is a solution of the approximated problem on
. We then have in the same fashion as we got ( 9.1 ),
By induction, we then see that we get a solution
defined on
,satisfying the estimate
|
(9.4)
|
establishing the independence of the time of existence respectively to
, since
does not depend on
. In the following we will note
.
10 Existence for ( 3.2 )
-
Proof.
We can here choose to take
, and let
. By the bound ( 9.4 ) independent of
on
, we then have the existence of a weakly convergent subsequence of
in the reflexive Hilbert space
, to a limit that we call
, which also belongs to
and satisfies the estimate
The usual compactness theorems ensure at this stage that
is a solution of ( 3.2 ) on
. The smoothness of our solution ensures that the solids do not collide with each other (if there is more than one) or the boundary (for an eventually smaller time), which establishes the existence part of Theorem. 5.1 □
11 Uniqueness for ( 3.2 )
-
Proof.
Since we cannot use a contractive mapping scheme for our problem, we have to establish uniqueness separately. Let then
denote another solution of ( 3.2 ) in
. Then, taking
as a test function in the variational formulation of the difference between the systems ( 3.2 ) associated to each solution provides for
:
| |
| |
|
(11.1)
|
For the viscous term in the fluid, we write
which with the
control of
and
provides us with an estimate of the type (where
denotes once again a generic constant)
| |
|
(11.2)
|
Concerning the forcing term in the fluid, we first notice that if we still denote
as
,
leading us to
| |
| |
with
. We have
. Moreover
. We then have by invariance by homotopy of the Brouwer degree (for the parameter
)
which together with the regularity of
establishes that
and that
for almost all
. Thus,
which with the
control of
and
provides
Consequently,
| |
implying
|
(11.3)
|
Concerning the elastic term,
| |
with
| |
| |
| |
| |
where we have used the
control of
and
for the inequality.
Next, for the same reasons,
| |
| |
| |
| |
| |
| |
We then write for the second term on the right-hand side of the last equality
, to get by Korn's inequality
| |
Similarly,
| |
| |
Thus,
| |
|
(11.4)
|
Concerning the pressure term, with
and the
control of the pressure, we get
| |
|
(11.5)
|
In order to get the estimate of
in
, we have to introduce the time differentiated problem. By taking
in the variational formulation associated to the difference between the time differentiated systems, we obtain
| |
| |
|
(11.6)
|
For the fluid viscous term, we easily find with the
control of the first time derivative of the velocity that
|
(11.7)
|
Concerning the forcing term in the fluid, since
(with a similar formula for
), we then deduce in a way similar to the steps leading to ( 11.3 ) that
| |
|
(11.8)
|
For the elastic term, we can also essentially reproduce the arguments leading to ( 11.4 ), leading us to
| |
|
(11.9)
|
The pressure term will require more care since we want to avoid the introduction of
that the most direct method would lead to. To do so, we notice that
| |
with
| |
| |
| |
For
, we have in a way similar to ( 11.5 ),
| |
| |
For
, the
control of
provides us with
For
we have:
| |
| |
where we have used the relations
in
for the first integral.
By integrating by parts in time,
| |
| |
With the
control of
we have
| |
| |
| |
| |
The remaining terms are more delicate. We first have
| |
| |
|
(11.10)
|
The apparent problem here is that
is estimated in
in terms of
in
. Now, a bound of this quantity in
will require a bound of
in
. In order to get such an estimate, we will bound
in
by lower order terms in
. To do so, let us first estimate the trace of
on
by using the test function
in the difference between the variational problems satisfied by
and
. By proceeding as in Section 10 , we would then get an estimate of the type, where
is given:
| |
| |
| |
which by patching all the charts defining
leads to an estimate of
in
yielding by elliptic regularity:
| |
| |
| |
Thus, with a choice of
small enough, we have for
small enough and the use of Gronwall's inequality,
| |
| |
By using this estimate in ( 11.10 ), we then get for a time small enough
| |
| |
By putting together the estimates on
,
and
, we have
| |
|
(11.11)
|
Now, by considering the difference between the two variational forms satisfied repectively by
and
, and writing the difference between the pressure terms as
| |
the Lagrange multiplier Lemma 13 of [5] provides for all
:
| |
|
(11.12)
|
By putting together the estimates ( 11.1 )-( 11.12 ), we then obtain for
small enough an inequality of the type:
| |
which shows that
on
. Let
If
, we can repeat the same procedure with
replacing
, which would lead to uniqueness for
as well. Thus, we have
, which concludes the proof of the theorem. □
12 Optimal regularity on the initial data
We first remind some extensions and regularization results on domains:
Lemma 12.1.
Let
be a domain of class
. Then, there exists a linear and continuous operator
from
into
(for each
) such that
in
. Also, if the
norms of a family of domains stay bounded, the norms of the corresponding linear operators also stay bounded.
Lemma 12.2.
Since
is of class
, let
(
) be a collection of charts defining a neighborhood of its boundary. We note
Then, there exists a sequence of domains (
) of class
, so that
, and which are defined with a collection of charts
(
) so that
as
. We then denote the complementary of
in
by
and
. We also assume
large enough so that the different connected components of
(if there is more than one solid) do not intersect each other or the boundary of
. We denote
.
We now state the optimal regularity assumptions needed in our analysis, and explain the adjustements required to the previous proofs.
Theorem 12.1.
With the same assumptions as in Theorem 5.1 , except for the following concerning the regularity of the initial data:
|
(12.1a)
|
|
(12.1b)
|
the conclusion of Theorem 5.1 holds.
Remark 8.
We have chosen here to take different forcings for the fluid, that we still note
with the same assumptions as in Theorem 5.1 , and the solid, in order to stress out that the higher order regularity required indeed comes from the hyperbolic scaling of the Navier-Stokes equations. The somewhat not so natural condition
is made in order to get
associated to the condition
.
-
Proof.
The idea is to first regularize the domains, initial data and modify the forcings in an appropriate way, and then pass to the limit.
Given
with
, we define as usual
.
We first notice that
,
,
and
still have the same regularity in
as in Theorem 5.1 . We first define in
,
and
,
,
(which is permitted since
. We next define
in
,
|
(12.2a)
|
|
(12.2b)
|
|
(12.2c)
|
|
(12.2d)
|
where
denotes the unit normal exterior to
, and finally
by
where the time derivatives on the right-hand side are computed with the usual rules from
,
(
,
.
We next define in the solid
by
|
(12.3a)
|
|
(12.3b)
|
|
(12.3c)
|
where the right-hand sides of the previous boundary conditions come from the fluid regularization previously carried. Note also that
|
(12.4)
|
(with an estimate that may blow up as
) since
| |
| |
We can then define
in
by
| |
| |
| |
| |
with the same conventions as for the previous system for the time derivatives evaluated from
and
evaluated from
.
We then define in
,
| |
| |
| |
where the time derivatives on the right-hand side are evaluated with
,
. We also define the regularized forcing in the solid
| |
We then have
,
,
in
and
in
,
, with
|
(12.5a)
|
| |
|
(12.5b)
|
|
(12.5c)
|
|
(12.5d)
|
where
is a given polynomial expression of
and
. We briefly explain how those constants appear. For instance, for the first estimate of ( 12.5c ), we have by elliptic regularity on ( 12.3 ) that
is bounded by a sum of terms, one of which being
,
being a polynomial which does not depend on
. Next, still by elliptic regularity on ( 12.2 ), we have that
is bounded by a sum of terms such as
.
This particular term, by the properties of the convolution, is in turn bounded by
. This shows that a term of the type
appears in the sum of all terms bounding
. Since the other terms in the sum can be dealt with similarly, this explains our estimate ( 12.5c ).
For the pressures, we have
| |
|
(12.6)
|
Since the initial data
and forcings
,
,
are smooth enough to ensure the regularity properties ( 12.5 ), we then deduce that we have similarly as for theorem ( 5.1 ) the existence of a solution
of a system similar to ( 7.7 ) with
,
,
,
being replaced by their counterparts with an exponent
, and
,
,
being replaced by
,
,
(with the choice
) given by
| |
| |
| |
| |
| |
| |
where the time derivatives are computed with a velocity satisfying
,
and a pressure such that
,
. Note that by construction, the solutions
to these problems in
satisfy
,
,
,
. Next, we proceed to energy estimates similar to Section 8 . The bounds obtained are similar, except that this time the terms associated to
,
and
tend to zero as
. This is clear from the convergence results ( 12.5 ), ( 12.6 ) for the integral terms associated to the fluid.
The terms associated to the solid asymptotically tend to zero by properties of the convolution. For instance, with the notations of Section 8 , for
, we get after change of variables, an integration by parts in time, three integrations by parts in space:
| |
| |
and thus with our estimate ( 12.5c ), we have
| |
where
denotes the same type of space as
with
and
being replaced by their counterparts with an exponent
. This type of estimates thus shows that this term does not change the energy inequalities in Section 8 . We can thus reproduce the arguments of Section 9 , establishing that
can be defined over a time
independently of
, and that its norm in
depends solely on
and thus, thanks to the estimates ( 12.5a ), ( 12.5b ), ( 12.5d ), solely on
. We can then consider the sequence
which is bounded in a space similar as
, but defined on
and extract (with respect to
) a weakly convergent sequence in a space modified from
by replacing the condition
by
. By the classical compactness results, we next see that the weak limit
and is a solution of ( 3.2 ) with
as forcing and
. This solution is also unique in
. □
13 The case of incompressible elasticity
In this Section, we explain how to treat the supplementary difficulties appearing when the incompressibility constraint is added in the solid. This leads to the same system as ( 3.2 ), with the addition of the condition
a.e. in
and the addition of
on the left-hand side of ( 3.2d ) and the addition of
(the trace of
being from the solid phase in this new term) on the left-hand side of ( 3.2e ). We now state our result and explain how to overcome the additional difficulties related to this constraint.
We first update our functional frameworks. While
and
do not change,
and
become respectively
| |
| |
| |
| |
Remark 9.
Whereas the pressure in the solid satisfies
, it appears that the limit pressures
are controlled uniformly in the norm of
and seemingly not in these norms. Note also that whereas the velocity field is smoother in the fluid phase for the solution of our next theorem, the pressure field is actually smoother in the solid phase. Whereas our artificial viscosity smoothes the velocity field in the solid, it also interestingly makes the pressure in the solid for the regularized system less smooth than the one associated to the solution of the constrained problem, which is source of difficulties that we shall describe later.
We now state our result:
Theorem 13.1.
With the same regularity assumptions as in Theorem 12.1 and assuming that the compatibility conditions associated to our new system at
hold (for the sake of conciseness we do not state them here), the conclusion of Theorem 5.1 holds for the case where the incompressibility constraint is added to the solid part.
Furthermore,
.
-
Proof.
The extra regularity (with respect to the norm of
) on the pressure in the solid simply comes from the equation
which once the regularity for the solution
is known provides immediately the result. We now explain how to obtain a solution in
.
The beginning of the proof follows the same lines as for the compressible elasticity case. We first assume that the initial data satisfies the regularity assumptions of Theorem 5.1 , and define the same smoothed problem as ( 7.1 ) with the corresponding updates for the incompressibility constraint. We then define the same fixed point linear problem as ( 7.4 ) where the condition
in
is added (the
being computed from the given
) and add
on the left-hand side of ( 7.4c ) and
(the traces being taken from
) on the left-hand side of ( 7.4d ).
We then proceed as in [5] to construct a solution to this system by a penalty method (the penalty term being this time defined over
) and get the same type of regularity result. This provides us with a solution
, that we also denote by
, for the incompressible version of ( 7.1 ) on a time
shrinking to zero.
As for the compressible case,
is in
, and since our smoothed problem has a parabolic artificial viscosity, we also have for the velocity in the solid the regularity
(
) (with estimates that blow up as
). Thus,
with
endowed with the norm
We next proceed as in Section 8 to get energy estimates, that will be carried this time for the
dependent norm of
, independently of
on
, and for such a purpose it is important to keep the
factor in the definition of the norm. We could extend the sum to
, whereas it is not necessary.
As before, the first set of estimates has to be carried on the highest order time derivative. Our energy inequality ( 7.8 ) has the same form, except that the integrals over
where
appears has to be taken this time on
. The part over
is estimated as before. We now explain how to deal with the integrals set on
for the pressure, which indeed needs some justifications given that the velocity in the solid is not controlled uniformly in
in a space as smooth as the velocity in the fluid, while the pressure is controlled in the same type of spaces in both phases.
-
Estimates on
Here
denotes any time in
. The most difficult integrals set in
and associated to the incompressibility constraint in the solid are
and
, the others being either less difficult or similar to estimate.
Step 1. For
, if we denote
, we have
| |
| |
| |
|
(13.1)
|
|
(13.2)
|
where we have used the continuity of
in the sense of traces along
to bound the
norm of
by means of the
norm. Note that we have also used the fact that the
norm of
is in the definition of the norm of
. In order to get an estimate on this norm, we would proceed in a way similar as to get ( 8.8 ) in Section 9 .
Step 2. Concerning
, we have by integrating by parts in space
| |
since our artificial viscosity provides the regularity
and
(with estimates that may blow up as
).
The difficulty here comes from the second integral. Whereas as for
we can estimate the trace of
on
from the fluid, we have to take the norm of
in
, which is problematic given that the norm
contains only its
norm. In order to circumvent this difficulty, we notice that the same formula holds if we replace
by
(the extension to
of the velocity in the fluid). Since
on
, we have
on
, which implies:
| |
| |
and thus,
| |
| |
| |
| |
| |
| |
| |
The most difficult integral set at time
on
and containing
is
for which we apparently just have an estimate of the type
(without any small parameter in front). We now explain how to treat this difficulty.
Step 3. We first notice that
| |
If we could say that
is
controlled, the
control of
would give us a suitable bound for
. Whereas we have seen in the statement of our theorem that
for the limit solution is indeed in
, we cannot seemingly get such a bound on the approximate pressures
. In order to get around this, we introduce similarly as in the previous step the extension to the solid domain of the velocity in the fluid. Since a similar integration by parts formula holds when we replace
by
, we deduce
|
(13.3)
|
The easier term to estimate is
, for which we have for an arbitrary
:
| |
| |
| |
| |
| |
| |
|
(13.4)
|
For the first integral, the nonlinear elastodynamics equation in
provides
| |
leading us for
to (since
in virtue of
),
| |
|
(13.5)
|
The integrals on the first line of this equality do not give any trouble and can be estimated in the same fashion. For instance, we have for
| |
| |
|
(13.6)
|
Now, the difficult term to handle is
. We first write the divergence form
, and integrate by parts:
| |
| |
| |
leading us to
|
(13.7)
|
and thus by putting together ( 13.5 ), ( 13.6 ) and ( 13.7 ),
| |
|
(13.8)
|
|
(13.9)
|
Now, the apparent problem comes from the term
on
that should be taken in
, which is troublesome since the norm in
appropriate for our limit process only contains its
norm. In order to circumvent this, we notice that we also have, since
is at least as smooth as
in
,
| |
| |
|
(13.10)
|
leading us, since
on
, to
| |
| |
|
(13.11)
|
Thus, by using ( 13.3 ), ( 13.4 ) and ( 13.11 ), we have
| |
|
(13.12)
|
Thus, we finally arrive to estimates analogous to ( 8.7 ) and ( 8.8 ), with the right-hand side being of the same type as in ( 13.12 ).
-
Estimate on
and
.
With the same arguments as in the next subsection, we have for
:
| |
| |
| |
|
(13.13)
|
We now explain on the case of the highest space derivative how to obtain elliptic estimates independent of
, since the addition of the pressure term does not allow to use Lemma 6.1 directly in the present case.
-
Estimate on
in
.
-
Regularity of the trace of
. First, by proceeding as in Section 8 , and as for the case of the highest order time derivative, we get an estimate for the trace similar to ( 8.17 ), with a majorant of the same type as in ( 13.12 ). We explain hereafter how to handle the estimates related to the pressure in the solid in order to get this trace estimate since difficulties different than in the higher order time derivative problem appear in the higher order space derivative problem.
Step 1. Let
.
Then,
with
| |
| |
| |
For
, we first notice that for
, if
is the sign of the permutation between
and
if
are distinct, and is set to zero otherwise,
| |
| |
| |
| |
| |
| |
| |
where we have use
on the second equality and
on the third one. Thus,
| |
| |
| |
| |
| |
showing
| |
In order to estimate the remaining term, we notice by integrating by parts twice for
that
| |
| |
| |
Since
, we then infer
| |
| |
Now, if we note
, we also have for
a similar formula. Since
on
, we then have
| |
| |
leading us to
| |
| |
| |
| |
Step 2. We see by integrating by parts with respect to the direction
that we have
| |
Step 3. Next,
, where
| |
| |
We first have
| |
For
the divergence condition
on
implies
| |
which in turn provides,
| |
which concludes the estimates on the pressure terms in the solid, justifying why we obtain a trace estimate similar as ( 8.17 ) with a majorant of the type of the right-hand side of ( 13.12 ). Now, we turn our attention to the recovery of the regularity in the solid, which will need some justifications since we cannot directly apply Lemma 6.1 .
-
Regularity in the incompressible solid First, with the introduction of
and of
, the solution in
of
we have for the nonlinear elastodynamics
i.e.,
|
(13.14)
|
We now apply Lemma 6.1 to this equation, leading us to
| |
| |
and, with
, to
| |
|
(13.15)
|
We then want to use elliptic regularity on the system:
|
(13.16a)
|
|
(13.16b)
|
|
(13.16c)
|
where the trace on
is estimated as we explained in the previous subsection.
Now, for the divergence condition in
, we notice that:
| |
| |
For the apparently problematic first term on the right-hand side, we first notice that
| |
| |
which with the condition
, provides
and thus
| |
We then deduce that
|
(13.17)
|
Now, with ( 13.15 ) and ( 13.17 ), elliptic regularity on ( 13.16 ) provides for
a bound of the same type as the right-hand side of ( 13.12 ), with however the norms in
being replaced by the norms in
, due to the term
appearing on the right-hand side of ( 13.15 ), that we bound by
| |
| |
We now turn our attention to the pressure, that we just need to control in
. In order to do so, we notice from ( 13.14 ) that we have for
:
| |
which with the previous estimate on
and
shows that we have a bound on
of the same type as the right-hand side of ( 13.12 ), but where the norms in
are replaced by norms in
due to the estimate in
of
. Now, elliptic regularity on the system:
| |
| |
| |
provides after integrating in time an estimate for
with a bound similar as in ( 13.12 ), still with the norms in
being replaced by norms in
.
Thus, we obtain for
the same type of estimate as well. Given our estimate on
, this also implies the same type of majoration for
.
Thus, we are lead to
| |
| |
which leads as in Section 9 to the introduction of a polynomial, this time of degree 4, which does not bring any substantial change with respect to Section 9 .
Note that the addition of
does not create any difficulty since a small
is chosen at the same stage as
, and the conclusion is similar as in Section 9 from the continuity of
on
which is established in the same way as the continuity of
. We then infer that there is a time of existence of
for our smoothed problems, with a bound on
and thus on
independent of
. Existence follows then by weak convergence in
and uniqueness can be established similarly as for the compressible case in Section 11 .
□
Acknowledgments
DC was partially supported by the National Science Foundation under grant NSF ITR-0313370. SS was partially supported by the National Science Foundation under grants DMS-0105004 and NSF ITR-0313370. References
-
J.T. Beale, Large-time regularity of viscous surface waves, Arch. Rational Mech. Anal., 84 (1983/84), 307–352.
-
A. Chambolle, B. Desjardins, M.J. Esteban, C. Grandmont, Existence of weak solutions for an unsteady fluid-plate interaction problem, Preprint.
-
P.G. Ciarlet, Mathematical Elasticity, Volume I: Three-Dimensional Elasticity, North-Holland, Amsterdam, 1986.
-
C. Conca, J. San Martin, M. Tucsnak, Existence of solutions for the equations modelling the motion of a rigid body in a viscous fluid, Comm. Partial Differential Equations, 25 (2000), 1019-1042.
-
D. Coutand, S. Shkoller, On the motion of an elastic solid inside of an incompressible viscous fluid, to appear in Arch. Rational Mech. Anal., http://arxiv.org/abs/math.AP/0407365.
-
C.M. Dafermos, W.J. Hrusa, Energy methods for quasilinear hyperbolic initial boundary value problems. Applications to elastodynamics, Arch. Rational Mech. Anal., 87 (1985), 267-292.
-
B. Desjardins, M.J. Esteban, Existence of weak solutions for the motion of rigid bodies in a viscous fluid, Arch. Rational Mech. Anal., 146 (1999), 59–71.
-
B. Desjardins, M.J. Esteban, C. Grandmont, P. Le Tallec, Weak solutions for a fluid-structure interaction problem, Rev. Mat. Complut., 14 (2001), 523–538.
-
D.G. Ebin, S.R. Simanca, Deformations of incompressible bodies with free boundaries, Arch. Rational Mech. Anal., 103 (1992), 61-97.
-
F. Flori, P. Orenga, Fluid-structure interaction: analysis of a 3-D compressible model, Ann. Inst. H. Poincaré Anal. Non Linéaire, 17 (2000), 753-777.
-
C. Grandmont, Existence et unicite de solutions d'un probleme de couplage fluide-structure bidimensionnel stationnaire, C. R. Acad. Sci. Paris Sér. I., 326 (1998), 651–656.
-
C. Grandmont, Y. Maday, Existence for unsteady fluid-structure interaction problem, Math. Model. Numer. Anal., 34 (2000), 609–636.
-
M.D. Gunzburger, H-C. Lee, G.A. Seregin, Global existence of weak solutions for viscous incompressible flows around a moving rigid body in three dimensions,, J. Math. Fluid Mech., 2 (2000), 219–266.
-
M. Heil, Stokes flow in an elastic tube – a large-displacement fluid-structure interaction problem, Int. J. Numer. Meth. Fluids, 28 (1998), 243-265.
-
M. Hillairet, D. Serre, Chute stationnaire d'un solide dans un fluide visqueux incompressible le long d'un plan incline, Ann. Inst. H. Poincaré Anal. Non Linéaire, 20 (2003), 779-803.
-
T.J.R. Hughes, T. Kato, J.E. Marsden, Well-posed quasi-linear hyperbolic systems with applications to non-linear elastodynamics and general relativity, Arch. Rational Mech. Anal., 63 (1977), 273-294.
-
W.J. Hrusa, M. Renardy, An existence theorem for the Dirichlet problem in the elastodynamics of incompressible materials, Arch. Rational Mech. Anal., 102 (1988), 95-117.
-
C. Liu, N.J. Walkington, An Eulerian description of fluids containing visco-elastic particles, Arch. Rational Mech. Anal., 159 (2001), 229-252.
-
M. Sablé-Tougeron, Existence pour un probleme de l'elastodynamique Neumann non lineaire en dimension 2, Arch. Rational Mech. Anal., 101 (1988), 261-292.
-
D. Serre, Chute libre d'un solide dans un fluide visqueux incompressible: Existence, Japan J. Appl. Math., 4 (1987), 33–73.
-
V.A. Solonnikov, Solvability of the problem of evolution of a viscous incompressible fluid bounded by a free surface on a finite time interval, St. Petersburg Math. J., 3 (1992), 189–220.
-
H.F. Weinberger, Variational properties of steady fall in Stokes flow, J. Fluid Mech., 52 1972, 321–344.
E-mail address : coutand@math.ucdavis.edu E-mail address : shkoller@math.ucdavis.edu