About the motion of nonhomogeneous ideal incompressible fluids

About the motion of nonhomogeneous ideal incompressible fluids

Nonhnear Ano/ysic. Theory. Printed m Great Britain. Methods & Applications, Vol. 12, No. 1. pp. 43-50, 1988. 01988 ABOUT THE MOTION OF NONHOMOGE...

529KB Sizes 32 Downloads 60 Views

Nonhnear Ano/ysic. Theory. Printed m Great Britain.

Methods & Applications,

Vol.

12, No.

1. pp. 43-50,

1988. 01988

ABOUT THE MOTION OF NONHOMOGENEOUS INCOMPRESSIBLE FLUIDS

0362-546x/88 53.00+ .oO Pergamon Journals Ltd.

IDEAL

ALBERTO VALLI* Dipartimentodi Matematica, Universita di Trento, 38050 Povo (Trento), Italy

and Institute of Fundamental

WOJCIECH M. ZAJACZKOWSKI* Technological Research, Polish Academy of Sciences, 00-049 Warsaw, Poland

(Received 2 August 1985; received for publication 27 October 1986) Key

words

and phrases:

Nonhomogeneous

incompressible

fluids, Euler

equations,

method

of

characteristics.

1. INTRODUCTION IN

THIS

paper the following

System

of

eqUatiOnS

+ (u - V)u] + Vp = div v = 0 p,+v*vp=o u],,o = 0 PI,=0 = PO(X) bnlan = 0 p[u,

pf

in IO,T[ x n = Qr,

in QT, in QT, inS2, in &, on]O, T[ x as2=xT

(1.1) (1.2) (1.3) (1.4) (1.5)

(1.6)

is considered in a bounded domain &2C UP’with boundary a!2 Here u, = u * n, and n = n(x) is the unit outward vector normal to the boundary. The equations (l.l)-( 1.3) describe the motion of an ideal incompressible nonhomogeneous fluid where u is the velocity, p the pressure, p the density and f the external force field per unit mass. This problem has been considered by many authors: see for instance [14] under some restrictions for f; [4-51 in Holder’s spaces; [6] in Sobolev’s spaces of Hilbert type, and in c”; and some additional references about this problem and others related to it can be found in [4,5]. However we think that the proof presented here is simpler than the others already known; we show the existence of a solution by means of a constructive procedure (by proving that a suitable sequence of successive approximations converges), and moreover one has no difficulties arising from the geometry of the domain B and the dimension m, while the approaches which use as unknown the vorticity o = curl u need some additional procedures if & is not simply connected and X2 is not connected, or if m > 3. The idea of the proof is taken from [17,18], in which some results for the Euler equation for homogeneous fluids were proved. By following this approach, one has to solve two equations for the velocity u and the density p which are of the same type and can be integrated by means of the characteristics (hence they are reduced essentially to ordinary differential * Work partially supported by G.N.A.F.A.

of C.N.R. (Italy). 43

A. VALLI and W. M. ZAJACZKOWSKI

44

equations), and then one needs to consider two Neumann problems for an elliptic equation of second order: the first one gives the projection over the subspace of divergence free vectors tangential to the boundary, and the second is related to the pressure p. We recall that the use of a Neumann problem for controlling the behaviour of the pressure was used also in Temam’s paper [15], whose approach however is different from ours. Finally, we remark that uniqueness is proved by Graffi [lo] (see also [3]). We want to finish this introduction by remarking that if p,,(x) = p* > 0, then problem (l.l)(1.6) reduces to the usual Euler equations for homogeneous incompressible fluids. Hence our approach gives also a method for showing the existence of a (local) solution for these equations. This method is different from the others already known, for instance: the vorticity method (see [ll, 12, 161 for m = 2; [2] for m = 3); the method of finding geodesics over an infinite dimensional Riemannian manifold of diffeomorphisms (see [l, 91); the method which transforms the problem into an ordinary differential equation on a closed set of a Banach space (see [8]); the Galerkin method (see [15]); the recent abstract method of Kato and Lai [13] (in which it is observed that Lai’s method based on successive approximations require “a considerable amount of preparation”. We think on the contrary that our result in this sense is quite simple, however compare our linearization with equation (E’) in [13] and the related linearization in Lai’s thesis). 2. NOTATIONS

In this paper the following notations are used

I141,,r~ lI4lL”(o.T;w!(n)) = lbllLr.-,T7

Ilullw!(n) =

for1 E

N, 1 c r < m. For noninteger

1 we set

IlullWkdQ)

=

Ib\ll,r,dR

.

We shall assume also that O
-i;fp,

Finally, the Einstein convention

S pO(x) S sup pa = M2 < $03. R about summation over repeated indices is assumed.

3. EXISTENCE

(2.1)

OF SOLUTIONS

To prove the existence of solutions of problem (l.l)-( 1.6) we shall use the following method of successive approximations: p;+l + “k. v k+l = 0 inQ,, (3.1) in 9, > Pk+Yr=o = PaLI where vk is treated as a given function, such that div vk = 0, vk . nldR = 0; div(l/#+‘Vpk+‘) I@+’

i!g

= divf-

Djv,“D,v”

in(t) x Sz,

tE [0, ZJ,l

on{f} x as2,

= f - n + vfv,kD3ni

I

(3.2)

where vk and pk+ ’ are treated as given functions, and n is a regular extension of n to 3; ,,,;+I

+ vk . v&+1

= _ -

P Wk+lj,+)

= a(x)

1 k+l VP k+l

+f

inQr, in 9,

(3.3)

Motion of nonhomogeneous

ideal incompressible

45

fluids

where vk, pk+l and pk+’ are treated as given functions. The function wk+l constructed by (3.3) is such that in general div wk’I # 0 and wkcl - n\an # 0. Therefore the method of successive approximations is possible if we introduce the projection n determined by nWk+l

3

Wk+l

-

Vlpk+‘,

(3.4)

where qk’l is the solution of the problem Ag7k+l zz div wk+l drp

in S2

k+l

an

an

= wktl

0naQ

dR

and then we set vk” E zwkil. We assume also that v0 = - a. For getting easily some estimates for l/pk+‘, we can consider also the following problem 17f+’

+ vk

.v

rl

k+l

(3.5)

Tlk+%O = l/PO where rfk+’ E l/pkf’. Assume now that aS2 E C4, vk E C’([O, T]; W:(Q)), I > m, p. E W;(Q) and satisfies (2.1), = 0 on a!& Then the existence of the solutions to ~~~~~~)withdivu=Oin~andu.~l~~ problems (3.1) and (3.3) is proved by means of the method of characteristics. Moreover, by applying D’ to (3.1), (1~1G 2), by multiplying by IDypk+11r-2D)‘pk”1 and by integrating in 5;2 (and repeating the same procedure also for (3.3), and (3~9,) one gets easily

(3 *6)

+ c(lJrlk+‘l12.,t~k+1J13.r + Ilfl12.~)Jlwk+‘Ili,‘~ s ~llvkll*.illWk+‘II~,~ (here and in the sequel each constant will depend at most on 52, m, r, M1 and M2). Moreover, the representation of p k+l by means of the characteristics gives that Ml s pk+l(f, X) c MZ,

M;’ 6 qktl(t, X) G Mi’.

Remark also that pk+ * and vk+’ belong to C”([O, T]; W:(Q))) , since the characteristic Vk(t, s, x) belong to c([O, T] X [0, T] X W:(Q)) ( see, for instance, [8]). The same will be true also for I@+’ since we will see that Vp k+l E L’(0, T, W:(Q)). From the projection (3.4) we obtain also at once that

ll~k+‘l12,rasclJwk+‘I/z,r.

(3.7)

A. VALLIand W. M. ZAJACZKOWSKI

46

Now we shall consider problem (3.2). We write it in the short form div(qVp) =: Fz + 5

&

divf - DiV,D,Ui

= G s f * n + v,v,D,ni

in Q (3.8)

on aB,

and one verifies easily that fQF = _fanG. Hence there exists a solution of problem (3.8) since the necessary compatibility condition is satisfied (see, for instance, [4]). Then we can restrict ourselves to get an a priori estimate. For having the uniqueness of the solution, we assume also that (3.9)

p = 0.

I .Q rp[tiplying (3.8)1 by p, integrating O,2one gets

over Q and using (3.9) for estimating /~]]o.z in terms of (3.10)

IIVP~IO,, =z cWl10.~ + IlGllLZ). Writing the problem (3.8) in the form hp = -pVr).

Vp + pF

in s2, (3.11)

ap

on asi,

G=PG we have, by classical results about the Neumann

problem

l/&r d c(lbF - ~Vrl . VPIIL,+ Il~Gllz-~/r.r.m1 s 4llPVrt - VP!1.r +-~lp~l~,~+ I~PG~~~,~).

(3.12)

The last two terms in the right-hand side of (3.12) are estimated by cl~~~~~.~(l~Fll~.~ + 11G]Iz.,) (recall that W;(B) is an algebra for ST> m). The first term in the right-hand side of (3.12) can be estimated by

s cl~~~i~.rlirll~~.r(~~V~ilo.~ + llD2~hr). ~llPlI~.rll~ll~.rllVPIl1.~ Using the interpolation inequalities of Gaghardo-Nirenberg have that this last expression is estimated by

(see for instance f7], p. 312) we

4lPlI%r+ 4IlPl1I.rlIT1ll2.r)llVPllO,Z~ where

(3.13) Therefore,

from the above considerations

LEMMA 3.1. L’(0,

T;

Let

W#L?)),

as2Ec3, G E L’(0,

we have W:(Q)), r > m. Then

pEL”(O,T; T;

W;(Q)),

qEL”(O,1”;

there

W;(Q)),

exists a unique

FE

solution

Motion of nonhomogeneous

p E L’(0, T; Wj(9))

ideal incompressible

fluids

47

of problem (3.8), (3.9) such that =S~~~Y(lI~II~.~llrlll~,~> + ll~l1dll%

II&,

+ llGhr),

(3.14)

where a( +) is described by (3.13). Hence using the form of F we have proved the following. LEMMA 3.2.

Let &2 E C4, pk+’ E C”([O, 7j; W:(Q)), $+‘E c”([O, T’J; W?(Q)), UkE P([O, T]; WZ(sz)), f E L’(O,T; w:(Q)), r > m. Then there exists a unique solution of the problem (3.2) such that pk+’ E L’(0, T, W:(Q)), jn pk+’ = 0 a.e. in [O, ?‘j and Ilpk+%J g

~~~tllP”+‘Il~,~ll~k+‘ll*,~~ + lI~k+‘ll~,~l~llfl~,~ + Ibkll:.r>-

(3.15)

From (3.6), and (3.6), we have

IIP”+‘M2,r d IJ~dl2,~exptc~lull~,~,~,~), t E IO97%

Ilsk+‘~ot12,‘ s ll~O”112,~ e~P@l141Z,r,m,lh cE lo, II.

(3.16)

Using (3.15) and (3.16) in (3X5), we get (3.17)

where I-I(.) is a nondecreasing function of its arguments. Using (3.7) (for k) and integrating (3.17) with respect to t we get IIw~+~IIz,~ 6 exp(c~~w~~~*,~,~,~ )[W .) i’ (f +

ilflld dt + tl~il~.,]

(3.18)

0

Let cr > 1 and let fl~‘$,,~,r s o$&; then assuming that 1 d fo, to such that the right-hand side of (3.18) is smaller than o$z~~~,,,we obtain from (3.7)

II~k+lll*,r=G41Wk+111Z,r s c4l&,rt for each natural k. To prove the convergence

f E P, hd

of these sequences we have to consider the following problems: Rf+’

+ ok. VRk+’ + Ok- Vpk = 0,

.

(3.20) Rk+llt=O = 0, where

Rk+’

= pk+’

(3.19)

- pk, @kE gk - Uk-1, k

&

1;

E,k+‘+vk.VEk++ek-v?jk=O, Ek+ll,_O = 0,

(3.21) I

whereEk+l==~k+l-~k,k~l; Wf”

+ (vk * V)Wk+’ + (ok. V)Wk = -,p+lVpk

- ,,k+lvpk+l,

(3.22) Wk+’ Itso = 0,

A.

48 where

wktl

s

,,,k+l

_

,+,k

9

VALLI and W. M. ZAJACZKOWSKI

k+l - pk, k 2 1. Moreover we have P + Ek”‘Vpk) = -D,S,“D,u; - L&U;-‘&?;,

pk+l=

div(rlk+‘VPk+’

aPk+l rfk+l- an + p+’

(3.23)

= (efuf + uye#Qzj.

g

At last from (3.4) we get ok+1

=

nwk+le

(3.24)

By proceeding in the same way used for getting (3.6), from problems (3.20), (3.21) and (3.22) we have

where E depends on the right-hand side of (3.19). At last from (3.23) and (3.24) we get Il@%,r s ~llWkllLr.

llP”“%,F g ~(ll~k+~vPkll~,~+ Jv%,A Knowing that t s to, the above inequalities imply

Ilwk+‘(t)\ll,r S I5 1 4 of”/lf(r)11z,, dr) i’ Ilw”(r)lli,~ dr. ( I 0 After standard considerations it follows that the sequence {uk, wk, pk, pk} converges to a solution of the problems (3.1)-(3.3) for the limit functions u, W, p and p, and u = xw. It remains to show that for the limit functions we have u = W, that is div w = 0, w *nldn = 0. Taking the divergence of (3.3) and using (3.2) and (3.4) one obtains (div w), I- u - V div w = -D,u,D,D,g, and taking the normal component (W * FZ), +

By using the properties

of (3.3), projecting U a V(W

*n) =

in Qr

(3.25)

on as2and using (3.2) one has

U~~~~iDj~

(3.26)

on Cr.

of the function sp defined by (3.4), we have easily

$IldivwIlo.2s cll~4l0,&412.2 c cII~ullo.mWiv wlh2+ Ilw44/~.2,dR).

(3.27)

To obtain an estimate for w . r&o we introduce the following curves on aS2 %(C

r, Y) = u(t, V(G t7 Y>)

V(t, T7Y) = Y

Therefore

on CT,

(3.28)

on ask

we can write (3.26) in the form

-$qV(t, V(t,

t, Y)) = U,(f, V(t, t, Y))D~ni(Vt~y t7Y))DiqCfy

= K(f, V(C t, Y>>,

‘ttt ‘T Y)) (3.29)

Motion of nonhomogeneous where

r+ = we n.

ideal incompressible

fluids

49

From (3.29), by setting $(t, y; t) = r,$t, V(t, z, y)) and by remarking that tl@r/~o

s c ex~(c~~ll~ll~.,a.~,)llrcll~.~.~~

we have $11~1~1/?,2,&2s clHI*.r.~J~ exP(cr*fl~ll2.r,~,t,>flvcoli1/2.2.dR s 4412.~.Dc.10 exp(cf,llu112.r,~:,i,)(ltdiv40.2 +

lb - 4l~m~2 1.

(3.30)

On the other hand one has also that III41l12,2,a~ s c exP(ctOllUtjZ.,.~.t~)~~~~~t~*.~.~~ , hence from (3.27) and (3.30), knowing that div ~],,a = diva = 0 in 52 and w . n],=, = a . n = 0 on dft, we get that div w = 0 in R, w. ~&o = 0 on aS2. Therefore we have proved the following. 3.1. Let a& E C*, a~ W;(Q), diva = 0, a +nlao = 0, p. E W?(a), infa po > 0,fE L’(0, T; W:(n)), r 2 m. Then there exist r. E IO, T], to sufficiently small (see 0 E c”([O, toI; Wf(Q 1) n Cl([O, Gil; ante)), (3.19)), P E co(K)> toI; (3.1% such that (p, u, p) is a solution of (1.1) W(W) fl w (0, to; WQ)),P E w)? to; ww (1.6) in Q,,. Moreover one has

THEOREM

0 < itf

p.

If fE CO([O, 7J; Wf(S2)), then (p, u,p) is a classical solution.

c p(t, p E

x) c

sip

po < +m

c”(]O, to]; W:(Q))

V(t, X) E &.

and u E C’([O, &,I; ~~(~)~,

hence

Remark 3.1. This method can be used also for obtaining more regular solutions, for instance

by choosing the space Wi(sL) instead of W:(G),

for I> 1 + y, q > 1 (and &2 E C’+2).

Acknowfedgemenr-The authors would like to thank Professor Hugo Beirlo da Veiga for very fruitful conversations about the subject of this paper. REFERENCES 1. ARNOLD V. I., Sur la geomitrie differentietle des groups de Lie de dimension infmie et ses applications a l’hydrodynamique des fluides parfaits, Annfs Inst. Fourier Univ. Grenoble 16, 319-361 (1966). 2. BARDOSC. & FRISCH U., Finite-time regularity for bounded and unbounded ideal incompressible fluids using Holder estimates, Lecture Nores in Mathematics, 565, l-13, Springer, New York (1975). . 3. BEIRAODA VEIGA H. & VALL~A., On the motion of a non-homogeneous ideal incomnressibfe fluid in an external 1 force field, Rc. Semin. mat. Univ. Padova 59, 117-145 (1978). 4. BEIRAO DA VEIGA H. & VALLI A., On the Euler equations for non-homogeneous fluids (I), Rc. Semin. mat. Univ. Padova 43, 151-158 (1980). 5. BEIRAODA VEIGA H. % VALLI A., On the Euier equations for non-homogeneous Appl. 73, 338-350 (1980).

fluids (II), J. math. Anatysk

6. BEIRAODA VEIGA H. & VALLI A., Existence of C solutions of the Euler equations for non-homogeneous fluids, Communs partial d$$ Eqns 5, 95-107 (1980). 7. BESOV0. V., IL’INV. P. & NIKOL’SKIIS. M., Integral Representations of Functions and Imbedding Theorems, 1, Winston & Sons, Washington (1978); translated from Russian. 8. BOURGUIGNON I. P. & BREZIS H., Remarks on the Euler equation, J. funct. Analysis 15, 341-363 (1974). 9. EBIN D. G. & MARSDENJ. E., Groups of diffeomorphisms and the motion of an incompressible fluid, Arm. Math. 92, 102-163 (1970).

50

A. VALLI and W. M. ZAJACZKOWSKI

10. GRAFFI D., 11 teorema di unicita per i fluidi incompressibili, perfetti, eterogenei, Reutu Un. map argenr., 17, 7% 77 (1955). 11. JUDOVIC V. I., A two-dimensional problem of unsteady flow of an ideal incompressible fluid across a given domain, Am. math. Sot. Transl. 57. 277-304 (1966); in Russian in Mar. Sb. 64, 562-588 (1964). 12. KATO T., On classical solutions of the two-dimensional non-stationary Euler equation. Archs rarion. Mech. Analysis 25, 188-200 (1967). 13. KATO T. & LAI C. Y., Nonlinear evolution equations and the Euler flow, J. funct. Analysis 56. 15-28 (1984). 14. MARSDENJ. E., Well-posedness of the equations of a non-homogeneous perfect fluid. Communs partial diff- Eqns 1, 215-230 (1976). 15. TEMAMR. On the Euler equations of incompressible perfect fluids, J. funcr. Analysis 20. 32-43 (1975). 16. WOLIBNERW., Un theorbme sur I’existence du mouvement plan d’un fluide parfait, homogene, incompressible, pendant un temps infiniment long, Math. 2. 37, 698-726 (1933). 17. ZAJACZKOWSKI W. M., Local solvability of nonstationary leakage problem for ideal incompressible fluid, 2, Pacif. J. Math. 113, 229-255 (1984). 18. ZAJACZKOWSKIW. M., Solvability of the leakage problem for the hydrodynamic Euler equations in Sobolev spaces, IFTR Reports 21/1983, Warsaw (1983).