THE PRIMITIVE EQUATIONS ON THE LARGE SCALE OCEAN UNDER THE SMALL DEPTH HYPOTHESIS

Transcription

1 DISCREE AND CONINUOUS Website: DYNAMICAL SYSEMS Volume 9, Number, January 00 pp. 97 HE PRIMIIVE EQUAIONS ON HE LARGE SCALE OCEAN UNDER HE SMALL DEPH HYPOHESIS Changbing Hu, Roger emam,, Mohammed Ziane he Institute for Scientific Computing & Applied Mathematics, Indiana University, Rawles Hall, Bloomington, IN Laboratoire d Analyse Numérique, Université Paris-Sud, Bâtiment 45, 9405 Orsay,France Department of Mathematics University of Southern California, Los Angeles, CA Abstract. In this article we study the global existence of strong solutions of the Primitive Equations (PEs) for the large scale ocean under the small depth hypothesis. he small depth hypothesis implies that the domain occupied by the ocean is a thin domain, its thickness parameter is the aspect ratio between its vertical and horizontal scales. Using and generalizing the methods developed in [],[4], we establish the global existence of strong solutions for initial data and volume and boundary forces, which belong to large sets in their respective phase spaces, provided is sufficiently small. Our proof of the existence results for the PEs is based on precise estimates of the dependence of a number of classical constants on the thickness of the domain. he extension of the results to the atmosphere or the coupled ocean and atmosphere or to other relevant boundary conditions will appear elsewhere.. Introduction. In this article we are concerned with the Primitive Equations (PEs for brevity) governing the motion and state of the ocean in a three dimensional thin domain (thin in a sense to be made precise later on). he primary purpose of this article is to establish the existence and uniqueness of global strong solutions of the PEs, i.e. for arbitrary time, for a broad class of data which includes the most physically relevant ones. he PEs of the ocean are derived from the Navier Stokes equations, with Coriolis force, coupled with the thermodynamic equation and the diffusion equation for the salinity by taking into account both the Boussinesq and hydrostatic approximations; see e.g. [7], [9], []. In a series of articles by Lions, emam and Wang [,,, 4, 5], the PEs of the ocean, atmosphere and their interaction, are studied from the mathematical viewpoint (existence of weak solutions, long time behavior, etc.), as well as from the numerical viewpoint. In [] emam and Ewald derived the maximum principles of the PEs of atmosphere by changing the temperature variable to a potential temperature variable. However, despite these advances in the study of the PEs, the mathematical theory for the PEs is far from being complete. For example, the uniqueness of weak solutions and the global existence of strong solutions for the PEs remain unsolved. his article presents a step in this direction 99 Mathematics Subject Classification. 5Q5, 5Q0, 76D0. Key words and phrases. Nonlinear partial differential equations, Primitive equations, Oceanography, Atmosphere sciences. Current address of CH: Department of Mathematics, exas A& M University, College Station, X

2 98 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE taking advantage of the small shape ratio of the ocean, the vertical dimension being much smaller than the horizontal ones. As it is well-known (see e.g. [7], [9], [,,, 4, 5]), there are essentially two characteristics used in simplifying the general hydrodynamics and thermodynamics equations of both the atmosphere and ocean. he first one is that, for large scale geophysical flows, the ratio between the vertical and horizontal scales is very small, and thus the vertical momentum equation of the atmosphere or the ocean can be well approximated by the hydrostatic equation; this hypothesis leads to the PEs; one of our objectives in this article is to further explore the mathematical consequences of this hypothesis. Another important characteristics is the rotating effect, i. e., the Coriolis force, see e.g. [], [6], [7] and the references therein. In this article we will treat the case where the Rossby number is of order, O(), concentrating on the small depth characteristic. In the mathematical literature, there is a large amount of work devoted to the study of partial differential equations in thin domains, see e.g. [4, 5], [8], [9], [8], [0, ], and [, 4]. In this article, we will consider the PEs in thin domains following the methods developed by two of the authors in [, 4]. We choose to work with the PEs in the actual domains. Compared to the work in [, 4], the PEs possess some specific difficulties to circumvent, namely: (i) the physical domain has corners and a non-flat bottom topography, (ii) the non-local constraint (incompressibility condition) and the integral expression of the vertical velocity in PEs, which leads to a strong non-linear term, i.e. ( 0 ) divv dz. (.) x (iii) non-homogeneous boundary conditions and the mixed boundary condition on the surface of the ocean which are not considered in []. he principal linear part of the PEs is a (stationary elliptic) Stokes type problem which we call the GFD-Stokes problem. In [6], [7], one of the authors, established the H regularity of the solutions of the GFD-Stokes problem under various boundary conditions. In the companion article [6], we extend this result to the more realistic boundary conditions considered in this article. In order to overcome the difficulties (i)-(iii), we also generalize the anisotropic Sobolev-type inequalities from [] to fit in our problem, (see [6]), and devote great efforts to the study of the strong nonlinear term, see (.). In this article, we obtain the global existence of the strong solutions for the PEs with a large class of initial data and boundary conditions (acting as external forcing in the PEs), provided is small enough; see the main theorem in Section. Let us mention a very recent work on PEs in [], in which the local existence of the strong solutions of the PEs is proved for domains of the class studied in [6]. Our results in this paper also imply the local existence of strong solutions in the case is of order (see Remark 4.). We present an outline of this article: in Section we briefly recall the mathematical setting of the PEs, the corresponding boundary conditions and their nondimensional form, their weak formulation and the main result of this article. In Section we recall some inequalities on v,, S from [6]. In Section 4 we derive a priori energy estimates for the PEs which will guarantee the large time existence of solutions. In Section 5 we prove the existence of the strong solutions of the PEs

3 PRMIIVE EQUAIONS IN HIN DOMAINS 99 with homogenized boundary conditions. In Section 6 we shift the solution obtained in Section 5 to its original counterpart, and prove the main theorem... he Brief primitive reminder equations of the PEs. of thewe ocean first recall and their full mathematical set of equations setting. which describe the motion and state of the ocean, obtained by taking into account both the Boussinesq and hydrostatic approximations, and which are called the Primitive Equations of the large scale ocean, PEs for brevity, (see e.g. [9], [7], [],[]): + (v )v + w + Ω cos θk v + gradp = µ v + ν v t ρ 0 x, (.) p = ρg, (.) divv + w = 0, (.) + (v ) + w = µ t x + ν x + Q, (.4) S S S + (v )S + w = µ t x S S + ν S x, (.5) ρ = ρ 0 ( β ( 0 ) + β S (S S 0 )), (.6) Here g is the gravity acceleration, µ, ν are the effective molecular dissipations in the horizontal and vertical directions, µ, ν, µ S, ν S reflect similarly the heat and salinity diffusions; ρ 0, 0 and S 0 are reference values for the density, the temperature and the salinity; Ω is the angular velocity of the earth, θ (0 θ π) is the colatitude of the earth. Q is the heating source mainly through radiation (sun and clouds). he unknown functions are the horizontal velocity v and vertical velocity w, the temperature, the salinity S, the density function of the ocean ρ and the pressure p. hroughout this article, we use,, div to denote the two dimensional gradient, Laplacian and divergence operators on the horizontal plane, and use, and div for the corresponding D differential operators. he domain occupied by the ocean is M h R : M h = {(x, x ) Γ i : h(x, x ) < x < 0}, where Γ i is the interface between the ocean and the atmosphere, and h : Γ i R is the depth function, which is positive and sufficiently smooth. hroughout this article we assume that Γ i R is convex and smooth. he boundary of M h consists of M h = Γ i Γ l Γ b, where Γ b is the bottom of the ocean, and Γ l is the lateral boundary of the ocean. he boundary conditions for the PEs read (see e.g. []): At the interface of the ocean and atmosphere, Γ i : + α v v = g v, w = 0, (.7) S + α x = g, = 0. At the bottom of the ocean and at the lateral boundary, Γ b Γ l : v = 0, w = 0, n = S = 0. (.8) n

4 00 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE Remark.. With the same technique as those developed in this article, we could treat the PEs in thin domains with other boundary conditions. As indicated in Section, one advantage of the above-mentioned boundary conditions is that a Poincaré-type inequality holds for v,, S (with small constants for v and ). Other boundary conditions for which the Poincaré inequality does not hold necessitate different methods, and will be considered elsewhere. Now we briefly recall the non dimensional form of the PEs, following []. However, and departing here from [], we want to introduce and emphasize in this work the role of the small depth assumption for the domain occupied by the ocean. o this end, let U, 0, ρ 0 be reference values respectively of the horizontal velocity, temperature and density. We also consider L to be the reference value of the horizontal length and H to be the reference value of the vertical length. Let = H/L be the small parameter which will be used throughout this article. hen we set { v = v U w = w U = 0 S = S 0 S ρ = ρ 0 ρ p = ρ 0 U p, (.9) t = Lt /U, h = h/l = h/h, (x, x, x ) = L(x, x, x ), (.0) Q = (U /L)Q = µ = ν Re LU Re LU = µ = ν Rt LU Rt LU = µ S = ν S Rs LU Rs LU b = gl U β = β 0 β S = β S S 0 Ro = U LΩ, (.) α v = α v/l α = α /L g = Ug /Lg v = (U/L)g v g = ( 0 /L)g. (.) Here Re i (i=,) are Reynolds type numbers, and Rt i and Rs i (i=,) are the nondimensional eddy diffusion coefficients; Ro is the Rossby number, which compares the angular velocity of the earth to the dynamical behavior of the ocean. In the present work, we will keep the Rossby number constant. For the Coriolis force we could as well consider the β approximation, see e.g. [9]: f = f 0 + β 0 x where f 0 = cos θ 0, β 0 = a sin θ 0, a is the radius of the earth, θ 0 is the mean colatitude. hen, by easy computation, we obtain the following non-dimensional

5 PRMIIVE EQUAIONS IN HIN DOMAINS 0 form of the PEs, in which all the primes have been dropped: + (v )v + w + f t Ro k v + gradp = v + v Re Re x, (.) p = bρ, (.4) div v + w = 0, (.5) + (v ) + w = + t Rt Rt x + Q, (.6) S S + (v )S + w = S + S t Rs Rs x, (.7) ρ = β ( ) + β S (S ). (.8) Introducing and then dropping the primes, we denote the non-dimensional domain occupied by the ocean, by = {(x, x ) ( h(x, x ), 0)}, (x,x ) Γ i where = Γ i Γ b Γ l is now the boundary of the domain in the non-dimensional space. hroughout this article, we will assume that h is independent of, where is the small parameter representing the thickness of the domain. In the final result we will assume that h is a positive constant but in several parts of the article, we will make the following assumptions concerning h: there exist positive constants h, h, h, h such that h C (Γ i ), 0 < h h h, h x i h, h (.9) L x i x j h. L i= i,j= hen we can write the boundary conditions of the non-dimensional problem as follows: S + α v v = g v, w = 0, + α x = g, n = 0 on Γ i, (v, w) = 0, n = S (.0) n = 0 on Γ b Γ l, where α v, g v, α, g are given functions, α v 0, α 0. Remark.. In [], by taking a different scale in the vertical direction, the authors obtained the non-dimensional PEs on a domain independent of the small parameter. However the parameter is hidden in the definition of some physical constants such as Re, Rt, and Rs, see [] for the details. he influence of this small parameter was not taken into account in [], and, on the contrary, we want here to emphasize the role of. he vertical velocity w and the pressure p are two diagnostic variables that we determine by integration of (.4) and (.5), taking into account the boundary conditions. Hence, for w, we obtain w(t; x, x, x ) = w(v)(t; x, x, x ) = div 0 x v(t; x, x, z) dz, (.)

6 0 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE and 0 div v(t; x, x, z) dz = 0. (.) h For the pressure, we introduce a new unknown function p s, the pressure of the sea water on the surface of the ocean hen p s : Γ i R, p s (t; x, x ) = p(t; x, x, 0). 0 p(t; x, x, x ) = p s (t; x, x ) + bρ(t; x, x, z) dz. x (.) With these notations the primitive equations consist of (.)-(.8), in which p in (.) is given by (.), and w in (.), (.6) and (.7) is given by (.). he boundary conditions are (.0). he main unknown functions is now the prognostic function u = (v,, S) and p s. Remark.. his formulation of the PEs was first introduced in [], []. In this new formulation, p s plays the role of the Lagrange multiplier corresponding to the non-local constraint (.). In [5], the author worked directly on the original PEs with anisotropic Sobolev spaces for the vertical velocity, and proved in particular the existence of strong solutions for the PEs for small time under the constant surface pressure assumption, p s = P. Remark.4. In the following argument, we show that the average of the salinity S over is zero if the initial salinity has zero average. i. e., S dx = 0. (.4) Indeed integrating (.7), we obtain [ d S dx + (v )S + w(v) S ] dt [ dx + S + ] S Rs Rs x dx = 0. Using the interface and boundary conditions (.0), we can deduce easily that the second and third integrals vanish and hence d S dx = 0, dt herefore the average value of S is constant and is denoted by S: (volume of ) S = S dx = S t=0 dx. his simply expresses the conservation of the total amount of salt in the sea water. In (.7) we actually replace S by S S. We can easily see that this translation does not change the equation and the boundary conditions for S. After this translation, we have the natural constraint (.4) for S. Relation (.4) makes the Poincaré inequality available for S, see (.7).

7 PRMIIVE EQUAIONS IN HIN DOMAINS 0.. Mathematical setting of the PEs. In this section we first define the function spaces suitable for the mathematical setting of the problem, see [] for the details. We denote by H s ( ), s R, the usual Sobolev spaces constructed on L ( ), and L ( ) = (L ( )), H s ( ) = (H s ( )). We also denote by H0(M s ) the closure in the space H s ( ) of Cc ( ), the space of infinitely differentiable functions with compact support in. Motivated by the boundary condition for the velocity field in the PEs, we define C l+b,0( ) = {v C ( ) : v is zero near Γ l Γ b }. We now define the function spaces for the unknown function u = (v,, S). Let hen we define and V = {v C l+b,0( ) : div 0 h v dz = 0}. H = the closure of V in L ( ), H = L ( ), H = L ( )/R, H = H H H, V = the closure of V for the H norm, V = H ( ), V = H ( )/R, V = V V V. he norms and inner products for the spaces H and H i, (i =, ) are the L ones, denoted by (, ) and. he inner products and norms for the spaces V are given by ((v, v )) = v v dx, v = v dx, v, v V. For V and V we take the usual inner products and norms, i. e., ((, )) = ( + ) dx, = ( + ) dx., V. and ((S, S )) = S S dx S = S dx, S, S V. For simplicity, we use also (, ) and to denote the norm and inner product in H, ((, )) and in V, i. e. for the later, { ((u, u )) = ((v, v )) + ((, )) + ((S, S )) u = ((u, u)), for u = (v,, S), u = (v,, S ). Remark.5. Here we propose different definitions for the inner products and norms for V, V and V based on the observation that the functions v and S in V and V respectively satisfy a Poincaré inequality; this will be made clear in Section. For in V, the classical Poincaré inequality does not hold.

8 04 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE Now by the Riesz representation theorem, we can identify the dual space H of H (H i, i =,, ) with H (respectively H i with H i, i =,, ), i. e., H = H. hen we have V H = H V where the two inclusions are compact and injective. he following lemma characterizes the relationship between H and its orthogonal complement. Lemma.. Let H H = {v L ( ) : div be the orthogonal complement of H in L ( ). hen 0 h v(x, x, z) dz = 0, n h 0 h v dz = 0 on Γ i }, (.5) H = {v L ( ) : v = grad φ, φ H (Γ i )}, (.6) V = {v H ( ) : div 0 h v dz = 0, v = 0 on Γ l Γ b }. (.7) where n h is the horizontal unit normal (in the horizontal plane) on the boundary of Γ i. Proof. Notice first that, in (.5), we give a slight correction to [] on the space H, so we prove this lemma in detail. We first prove (.6); the space in the right hand side of (.6) is clearly orthogonal to V (in L ( )) and is thus included in H. o prove the opposite inclusion, we first notice that if v is orthogonal to V in L ( ) (v H ), then v is independent of z, / z = 0. Indeed for any Ψ in Cc ( ), φ = Ψ/ z is in V ; hence v, Ψ/ z = / z, Ψ = 0, i. e. / z = 0. Now any distribution on which is independent of z defines a distribution on Γ i : for Ψ Cc (Γ i ), we define v, Ψ Γi as v, Ψθ 0 M, where θ 0 is a real valued C function with compact support in ( h, 0) and such that + θ 0(z)dz =. For this definition we observe that v, Ψθ 0 M is independent of the choice of θ 0 ; indeed if θ 0, θ 0 are two such functions, then θ = θ 0 θ 0 is such that + θ(z)dz = 0, so that θ = σ/ z with σ C with compact support in ( h, 0) and thus Ψθ is in V. Finally if Ψ is in Cc (Γ i ), and divψ = 0, then Ψθ 0 is in V, for θ 0 as above, since div 0 h Ψθ 0 dz = divψ = 0. hus v, Ψ Γi = v, Ψθ 0 Γi = 0 for all such Ψ s and by the classical De Rham heorem (see e.g. []), v is the gradient of a distribution φ on Γ i. Finally, by Deny-Lions theorem (see the same references), the distribution φ has its derivatives in L (Γ i ) (v = gradφ), and thus φ H (Γ i ). his completes the proof of (.6). hen we prove (.5). First notice, from the classical theory of Navier Stokes equations, that for v L ( ), with 0 h v dz L (Γ i ) and div 0 h v dz L (Γ i ), then n h 0 h v dz is well defined and belongs to H ( Γ i ); so the constraint appearing in the right hand side of (.5) makes sense. Denote by Ḣ the right hand side of (.5); the inclusion H Ḣ is trivial, since V Ḣ, and Ḣ is a closed subspace of L ( ). Denote by Ḧ the orthogonal complement of H in Ḣ, and assume that v Ḧ. hen v H Ḣ, and by (.6), there exists a

9 PRMIIVE EQUAIONS IN HIN DOMAINS 05 φ H (Γ i ) such that v = grad φ L (Γ i ); it should be noticed that v is a D vector depending on (x, x, x ), and furthermore div n h 0 h 0 h grad φ dz = div(h grad φ) = 0, v dz = n h (h grad φ) = h φ n h = 0. (.8) By the classical elliptic theory, (.8) implies that grad φ = 0, i. e., v = 0, which yields immediately Ḧ = {0}, and hence H = Ḣ; the proof is complete. For (.7) we then follow the same methods as in [], using (.5) and (.6). hroughout this article, we will denote by, the semi-norm in H ( ) or H ( ), for any u H ( ) or H ( ), u, = u x i x j i,j= We now define the bilinear forms a : V V R, a i : V i V i R, i =,,, and the corresponding linear operators A : V V, A i : V i V i, i =,,, by setting: and a (v, ṽ) = A v, ṽ = [ Re v ṽ + Re ] ṽ a (, ) = A, [ = Rt + ] Rt dx + α v v ṽdx dx, Re Γ i dx + α Rt a (S, S) = A S, S [ = S Rs S + S S ] dx. Rs a(u, ũ) = Au, ũ = A v, ṽ + A, + A S, S, for u = (v,, S), and ũ = (ṽ,, S) in V. Γ i dx dx, One of the main objectives of [6] is to establish, under suitable hypotheses, the following characterization of the operators A i and their domains, D(A ) = {v H ( ) V : + α v v = 0 on Γ i }, D(A ) = { H ( ) V : + α x = 0 on Γ i }, (.9) D(A ) = {S H S ( ) V : n = 0 on }.

10 06 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE and [ A v = P v + Re ] v Re x, v D(A ), [ A = + ] Rt Rt x, D(A ), [ A S = S + ] S Rs Rs x, S D(A ). (.0) Here P is the orthogonal projection from L ( ) onto H ; the regularity of the operator A has been studied in [6] in a different context; one of the main issues in [6] is to derive the similar results for our current boundary conditions and the dependence on of the norms of the operators A i. he linear operators A i, i =,,, which are isomorphisms from V i onto V i, can be seen as unbounded, selfadjoint linear operators on H i ; they are positive operators, and they admit compact inverses, so that the fractional powers of A i can be defined; A i, i =,,, especially will be frequently used, and A i, i =,, provides an equivalent norm on V i, i =,, respectively. Hereafter we shall alternatively use the following notations, justified by a simple integration by parts, we have A v = a (v, v) = v + Re Re + α v v L Re (Γ i) v V, A = a (, ) = + Rt Rt A S = a (S, S) = S + S Rs Rs + α Rt L (Γ i) V, S V. We postpone the discussion on the coercivity and continuity of the bilinear forms a i, i =,, and a until some inequalities on v,, S are established, see Proposition 4. in Section 4. In relation with the nonlinear terms appearing in the equations, we define the following trilinear functionals and associated operators: and finally b (v, ṽ, u ) = B (v, ṽ), v = b (v,, ) = B (v, ), = b (v, S, S ) = B (v, S), S = b(u, ũ, u ) = B(u, ũ), u ) [ (v )ṽ + w(v) ṽ z [ (v ) + w(v) [ (v )S + w(v) S = b (v, ṽ, v ) + b (v,, ) + b (v, S, S ), ] v dx, ] dx, ] S dx, for u = (v,, S), ũ = (ṽ,, S) and u = (v,, S ) in V.

11 PRMIIVE EQUAIONS IN HIN DOMAINS 07 hese trilinear terms possess the following anti-symmetry properties: b(u, u, u ) = b (v, v, v ) = b (,, ) = b (S, S, S ) = 0 (.) which play a crucial role in the proof of the global existence of weak solutions, see [], []. It should be noted that the PEs possess some highly nonlinear terms such as w(v) which need particular attention. Concerning the Coriolis term in the equations, we define a bilinear functional e : H H R and the associated linear operator E : H H by setting e(v, v ) = (E(v), v ) = f (k v) v dx, v, v H. Ro It is easy to see that e(v, v) = 0, v H. We are now in a position to recall the weak formulation of the PEs and to state the main results in this article... Weak formulation of the PEs and the main results. With the above notations, the weak formulation of the PEs is given by: Find u L (0, t ; H) L (0, t ; V ), t > 0, such that: d dt (u, u ) + a(u, u ) + b(u, u, u ) + e(v, v ) = g v, v Γi + g, Γi + (Q, ), u = (v,, S ) V, u(0, x) = u 0 (x) H. (.) We recall the existence of weak solutions for the PEs which is due to Lions, emam and Wang in []: heorem.. For any t positive, there exists at least one solution for problem (.), u = (v,, S) (not necessarily unique) defined on (0, t ), and such that u C([0, t ); H w ), u t L (0, t ; (V H ( )) ), where H w is the space H equipped with the weak topology and (V H ( )) is the dual space of V H ( ). In this article we are concerned with the global strong solutions (.), namely: { Find u L (0, t ; H ( )) L (0, t ; H ( )), t > 0, (.) such that (.) holds, where u 0 is now given in V. he terms in the right-hand side of (.) come by integration by parts from the boundary conditions on v and. he weak formulation will be sufficient to consider the weak solutions of the PEs. Since we are interested in the strong solutions of the PEs, it is necessary to homogenize the boundary terms, then the integral terms over the boundary will be replaced by some integration over the whole domain ; see (4.)-(4.6). It should be noted that the uniqueness of the weak solution is an open problem. o state the main theorem, let us introduce an arbitrary monotone positive function R 0 () from (0, ) into R + satisfying lim 0 R 0() = 0. (.4)

12 08 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE We are now in position to state the main theorem. heorem.. (Main result) he assumptions are those above; in particular is convex and the depth h is constant, and let g v H 0(Γ i ), g H (Γ i ), Q L ( ). Finally assume that the initial data and boundary conditions satisfy A v 0 + A 0 + A S 0 + Q + [ g v 4 L (Γ i) + g 4 L (Γ i) + g 4 L (Γ i) ] R 0(). (.5) hen there exists an 0 depending on Re i, Rs i, Rt i, α v, α, Ro, i =, and on the function R 0, and such that for 0 < 0, there exists a unique global strong solution which satisfy (.), i. e., for any t > 0, u L (0, ; H ( )) L (0, t ; H ( )), Moreover we have the following estimate on the solution u: A u(t) σr 0(), t 0, where σ > (independent of ) will be determined below (see (5.)). Remark.6 he convexity of and h constant are used in [6] to show (.9)- (.0) including the dependence on of the norms of the operators A i. However heorem. is valid for any domain and any set of boundary conditions for which the analogue of (.9)-(.0) has been proven, e.g. the Dirichlet boundary conditions as in [6]. his remark extends also to the result in Remark 4. (local existence of strong solutions for = ). he extension of the results of [6] (and hence of the results in heorem. and Remark 4. in this article), to the case where is not convex and h is not constant, will appear elsewhere. Remark.7. he assumptions (.5) are physically relevant, they allow large initial data for which the vertical variation (for v 0, 0, S 0, Q ), is O(). Indeed such functions can be treated as independent of z so that e.g. A v 0 = O( ). Hence v 0 independent of z can have a norm as large as C R 0 (), and we observe that R 0 () = α, α > 0, is suitable for (.4). A similar remark holds for functions slowly varying in the z-direction and for other data. Remark.8. he hypothesis that Γ i is convex is not natural, and is due to Propositions. and. in [6]. Note that heorem. above and all what follows extend without any modification to all domains for which Propositions. and. of [6] is valid.. Functional inequalities. In this section we recall from [6] a number of inequalities on the functions v,, S, satisfying the boundary conditions proposed in the first section. Since we are interested in the nonlinear partial differential equations in thin domains, we need to know the exact dependence, with respect to the thickness of the domain, of the constants appearing the inequalities. As it is well known, in the classical theory of Sobolev spaces, the isotropic feature is emphasized. In our problem, one of the main features is the natural anisotropy of the domain between different space directions. In [], the authors derived anisotropic versions of the classical Sobolev-type inequalities, such as the Poincaré inequality, Agmon s inequality and Ladyzhenkaya s inequality, which have proved to be successful in the study of the Navier Stokes equations in thin domains. In the current work, we

13 PRMIIVE EQUAIONS IN HIN DOMAINS 09 need some inequalities with explicit dependence of the constants on the thickness of the domain. We present these inequalities without proof, their proof is given in a separate article [6]. Here and after, c 0 will stand for some numerical constant, and C some constant depending on some non-dimensional physical quantities appearing in the PEs, namely, α v, α, Re i, Rs i, Rt i, i =,, and Ro. he constant c 0 is independent of the small parameter and may vary at each occurrence; C is a constant which possibly depend on. First we describe some inequalities that will be used for a linearized version of the primitive equation, the goal of which is to homogenize the boundary conditions and the initial conditions for the velocity and the temperature. For v satisfying the following boundary conditions: v = 0 on Γ b Γ l, + α v v = g v on Γ i. (.) he Poincaré-type inequalities: v, v H ( ), (.) x i v v x, for i =,, v H ( ), (.) i + α v v dx dx v Γ i x + 4 g v, v H ( ). (.4) he Grisvard-Iooss inequality: i,j= ( v x i x j ( c 0 he Agmon inequality: ) dx + α v i= v dx dx Γ i ) v x dx + c 0 ( g v L (Γ + g i) v L (Γ ). i) j (.5) v L c 0 v + c 0 ( gv L (Γ i) + g v L (Γ i)). (.6) he Ladyzhenskaya inequality: v L 6 () c 0 v, v V, (.7) v L6 () c 0 A v, v H ( ). (.8) Now we recall some inequalities concerning a function which satisfies the following boundary conditions on : n = 0 on Γ l Γ b, + α = g on Γ i. (.9) he classical Poincaré inequality does not follow directly. he following inequality will play a role similar to Poincaré s inequality.

14 0 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE i= he Poincaré-type inequalities: c 0 L (Γ + c i {0}) 0, c 0 + c 0. x i i= + α he Grisvard-Iooss inequality: ( ) dx + α i,j= x i x j c 0 ( (.0) Γ i dx dx + α g L (Γ i). (.) i= We have some further inequalities on : and x i Γ i dx dx ) x dx + c 0 g L (Γ + c i) 0 g L (Γ. i) i (.) C g + c 0, c 0 + c 0 g L (Γ i) + c 0 g L (Γ i), (.) + α dx dx + α g L (Γ i) Γ i ) c 0 ( g L (Γ i) + g L (Γ i) + x + α g i i= c 0 + c 0 g L (Γ i) + C g L (Γ i). (.4) Finally we have the following version of the Ladyzhenskaya inequality for : L 6 c 0 6 L (Γ i {0}) + c 0. (.5) We now present some inequalities on a function S satisfying S dx = 0, together with the following boundary condition on : S n = 0 on. (.6) he Poincaré-type inequalities: S c 0 S, he Grisvard-Iooss inequality: i,j= ( S x i x j ) dx c 0 S c 0 ( i= i= S x i. (.7) ) S x dx, S H ( ). (.8) i

15 PRMIIVE EQUAIONS IN HIN DOMAINS We conclude this section with the Ladyzhenskaya inequality for S: S L 6 c 0 S, S L 6 c 0 S. (.9) 4. A priori estimates. Before we start to derive a priori estimates for the PEs, we need to homogenize the boundary conditions on v,. For that purpose, we introduce v,, that are solutions of the following problems: t ( v + v x ) + p = 0 in, 0 div v dz = 0 in Γ i, h (4.) v = 0 on Γ l Γ b, + α v v = g v on Γ i. x v (, 0) = v 0, and ( + t x ) = 0, in n = 0 on Γ l Γ b, + α x = g on Γ i, (, 0) = 0. (4.) he regularity of solutions of the stationary problems associated to (4.) and (4.) has been investigated in [6], namely, v L (0, t ; H ( )) C([0, t ]; H ( )), L (0, t ; H ( )) C([0, t ]; ), t. herefore, one has, thanks to (.)-(.4), t 0 v (t) c 0 v 0 + c 0 [ g v + g v ]. A v ds c 0 t [ g v + g v ] + v 0 + c 0 g v. We set v = v + v, = + and p = p + p, and after dropping the primes, the new system for v,, S, p reads + (v )v + w(v) + f t Ro k v + gradp s (4.) + (v )v + (v )v + w(v) + w(v ) 0 + b x grad( β + β S S) dz Re v Re v x = f v, + (v ) + w t Rt Rt x (4.4) + (v ) + (v ) + w(v ) + w(v) = f x + Q, S S + (v )S + w(v) + (v )S + w(v ) S (4.5) t S S = 0, Rs Rs div 0 h x v dz = 0. (4.6)

16 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE where f v = (v )v + w(v ) + f Ro k v (4.7) 0 + b grad( β ) dz, x f = (v ) + w(v ). (4.8) he homogenized system is supplemented with the following homogeneous boundary conditions on : and initial data + α v v = 0, v = 0, + α = 0, n = S n = 0 on Γ b Γ l, v(, 0) = 0, (, 0) = 0, S(, 0) = S 0. S n = 0 on Γ i, (4.9) In this section, f v and f are assumed to be in L ( ). In Section 6, we will give estimates on f v and f needed to prove the main theorem in this article. Remark 4. In [], the authors assumed that div g v = 0, and they solved an ordinary differential equation to obtain v (the notation is different in []). In this work, since v is the solution of the time-dependent GFD-Stokes problem (4.), it is not necessary to assume this divergence-free condition on g v. For the readers convenience, we list in the following proposition some inequalities for the new functions v, and S satisfying the homogeneous boundary conditions (4.9). Proposition 4.. For v D(A ): C v A v C v, (4.0) v C A v C A v, (4.) v L C A v, (4.) For D(A ): v L 6 c 0 v C A v, (4.) v x i x j c 0 v C A v. (4.4) i,j= C A C, (4.5) C A C A, (4.6) L 6 C 6 A, (4.7) x i x j C A. (4.8) i,j=

17 PRMIIVE EQUAIONS IN HIN DOMAINS For S D(A ): A S = Rs S (4.9) S C A S C A S (4.0) S L 6 C A S, (4.) S x i x j C A S. (4.) i,j= Proof. Most of the inequalities presented here follow from the previous section by making the boundary conditions homogeneous. he new one are inequalities (4.4), (4.8) and (4.); their proofs are given in [6]. he following inequality will also be frequently used ( Γ i ( 0 ) ) 4 q dz dx dx h ( ) 4 q 4 4 dx, q L 4 ( ). (4.) We start with the a priori estimates for v. We multiply equation (4.) by A v, then integrate in x over. Noticing that (gradp s, A v) = 0, we are led to d dt A v + A v + (v )v A v dx + w(v) A v dx + (v )v A v dx + (v )v A v dx M + w(v) M A v dx + w(v ) A v dx f 0 + Ro k va v dx + b A v grad( β + β S S) = (f v, A v). x (4.4) In what follows we will write estimates for each term appearing in (4.4). (v )v A v dx C v A v, (4.5) (v )v A v dx v L v L 6 A v c 0 v / 6 A v + c 0 A v v, v / A v A v (4.6)

18 4 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE (v )v A v dx v L v A v C [ / A v + g v + g v ] v A v By Poincaré s inequality (4.0), we have 6 A v + c 0 [ / A v + g v + g v ] v 6 A v + c 0 A v v + c 0 [ gv + g v ] v. (4.7) ( ) f Ro k v A v dx f Ro v A v f Ro v A v. (4.8) Now we handle the most problematic term mentioned in the Introduction: w(v) A v dx w(v) A v dx ( 0 ) (divv) dz h A v dx ( 0 ) ( (divv) 0 ( dz Γ i h h }{{} =I, L 4 ( 0 (A v) dz h ) } {{ } L I L 4 I L 4 A v C v A v. ) dz) }{{} =I, L 4 dx dx (by Lemma 4. and Lemma 4.4 below) (4.9) Remark 4.. By making full use of the specific structure of the integrand, we obtain the above estimate on the trilinear term which can be controlled by the principal linear operator A. A rough computation without attention to this fact leads to the following w(v) A v dx w(v) L 4 A v L 4 ( ) w(v) 4 w(v) A v (w v) C v A v 5, and with this estimate, the nonlinear term is too strong to be dominated by the linear operator A. We continue the estimate, apply the same technique as for

19 PRMIIVE EQUAIONS IN HIN DOMAINS 5 (4.9), and we find M w(v) A v dx I L 4 So that Γ i ( 0 h c 0 v v A v A v, ( ) ) dz dx dx 4 A v (4.0) w(v) A v dx 6 A v + c 0 v v A v (4.) In the same way as in (4.0), w(v ) A v dx c 0 v v A v [ A v + [ ] ] g v + gv 6 A v + c 0 v v [ A v + [ g v + g v ]]. (4.) By the Cauchy-Schwarz inequality, we have [ 0 ] ( 0 ) A v ( β + β S S)dz dx A v ( + S )dz dx h h A v ( ( 0 h A v ( + S ). Using again the Cauchy-Schwarz inequality, we obtain ( + S )dz) dx (4.) (f v, A v) f v A v ( f v + A v ). (4.4) We collect (4.4)-(4.4) and obtain the final inequality: d ] dt A v + [ c 0 v A v f v + c 0 ( + S ) + c 0 [ A v + g v + g v + f ] Ro + v A v v with ϕ (t) + ψ (t) v. ϕ (t) = f v + c 0 ( + S ), ψ (t) = c 0 [ A v + + g v + g v + f ] Ro + v A v (4.5) Remark 4..(Local existence of strong solutions)we note here that inequality (4.5) is enough to establish the global well-posedness for small data and the local well-posedness. Initially the velocity v is zero and therefore its H norm stays small for a small time say on (0, t ). let t > 0, so that c 0 / v (t) on (0, t ). hen we have d dt A v ϕ (t) + ψ (t) A v. )

20 6 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE herefore, since v = 0 at t = 0, the following inequality hold for t (0, t ) ( t ) t A v (t) exp ψ (s) ds ϕ (s) ds 0 he integrability of ϕ and ψ follows from the existence of weak solutions and the linearity of the equations satisfied by v. hus choosing t, so that the right hand side in the above inequality is less than c implies the local existence of strong 0, solutions. he following lemma specifies the estimates on I, I used in (4.9). Lemma 4.. For v D(A), ( 0 w(v) (x, x, x ) (div v) dz h and 0 ) = I (4.5) I L4 (Γ i) c 0 v A v. (.40) Proof. Inequality (4.5) is a direct consequence of the Cauchy-Schwarz inequality. o prove (.40), we first recall Ladyzhenskaya s inequality in -D [0]: I L 4 (Γ i) c 0 I L (Γ i) I L (Γ i). We need to bound I L (Γ i) in terms of A v. o this end, we consider I x i, i =,. I x i = ( 0 ) [ (div v) 0 dz h h (div v) (div v) dz + (div v) x Γb i ( h )]. x i o control the boundary term, we need Lemma 4. below. Note div v + w = 0 (continuity equation), w = 0 on Γ i Γ b, so by Lemma 4. below, we obtain (div v) 0 Γb div v (div v) dz. (4.6) h We infer from (4.6) by using Hölder s inequality ( ) 0 [ ] I 0 (div v) dz + x i x i h which implies (.40). he lemma is proved. h [ ] (div v) dz. Lemma 4.. Given a < b, given f H 0 (a, b) H (a, b), then f satisfies f L (a,b) f L (a,b) f L (a,b). (4.7) Proof. If f C 0(a, b) C (a, b), using the existence of ξ (a, b) such that f (ξ) = 0, we easily obtain (4.7) for such a function. he lemma then follows by density. Lemma 4.4. Denote I = ( 0 h ( ) ) dz for v D(A ), then I L4 (Γ i) c 0 v A v.

21 PRMIIVE EQUAIONS IN HIN DOMAINS 7 Proof. Applying the same technique as in Lemma 4., we are led to estimate I x i, i =,, I x i = ( 0 h ( ) dz) [ 0 h ( ) ( v dz + Γb h ) ]. x i x i he main task here is to bound the boundary term by suitable integration on v. Using the boundary condition for v, we obtain ( ) ( Γb = + α v v) Γb 0 ( ) ( ) = + α v v + α v v dz h 0 ( ) h 0 ( ) + α v v dz + α v v + α v v dz h ( ( L + α v v) + c 0 v L L ( h,0) x L ( h,0) + α v v) ( h,0) ( h,0) ( c 0 L + α v v). L ( h,0) ( h,0) he remaining part of the proof is easy and will be omitted. We now derive the a priori estimates for the temperature. Multiply the equation (4.4) with A, then integrate over, to obtain A + A + (v ) A dx + w(v) A dx + (v ) A dx + w(v) M A dx + (v ) A dx + w(v ) A dx = (f x + Q, A ). (4.8) As we have done for v, we will estimate each term in (4.8) separately. First by Agmon s inequality for v, we obtain (v ) A dx v L L 6 A (4.9) c 0 v A

22 6 A + c 0 v A v A (4.4) 8 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE w(v) A dx w(v) A dx ( 0 ) (divv) dz h A dx ( 0 ) ( (divv) 0 ( dz Γ i h h }{{} =I, L 4 ( 0 (A ) dz h ) } {{ } L ) dz) }{{} =I, L 4 dx dx I L 4 I L 4 A (by Lemma 4. and Lemma 4.5 below) c 0 v A v A 6 A + c 0 v A v. (4.40) In the same fashion, by Lemma 4., w(v) M A dx I L4 (Γ x i) Γ i ( 0 h c 0 v A v A ( ) ) dz dx dx 4 A hen by Agmon s inequality, (v ) A dx v L A c 0 A v A 6 A + c 0 A v. (4.4) Similarly, (v ) A dx v L A 6 A + c 0 [ A v + g v + g v ]. (4.4)

23 PRMIIVE EQUAIONS IN HIN DOMAINS 9 ) ( Let I = Γ (div v i ) dz, then by Hölder s inequality, Lemma 4. and Lemma 4.5 below, we have w(v ) A dx I L 4 I L 4 A Finally c 0 v A v A A, (by (4.8)) 6 A + c 0 v A v (4.44) (f + Q, A ) f + Q A f + Q + 4 A. (4.45) Putting (4.8)-(4.45) together, we find d ] dt A + [ c 0 v ] A f + Q + + c 0 [ v A v A + A v ] + [ + c 0 v A v + A v + v A v ] + c 0 [ g v + g v ] ϕ (t) + ψ (t) A, where ϕ (t) = f + Q + c 0 [ v A v A + A v ] ψ (t) = c 0 [ v A v + A v + v A v + g v + g v ] Now we prove the estimate on I = which we used before. ( 0 h Lemma 4.5. Let D(A ), and denote ( 0 ) ( I 4 = 0 dz, I 4 = h hen we have the following estimates: ( h (4.46) ) ) dz (appearing in (4.40)), ( ) + α x dz). I L4 (Γ i) c 0 A, (4.47) I 4 L 4 (Γ i) c 0 A, (4.48) I 4 L4 (Γ i) c 0 A. (4.49) Proof. We first prove (4.48). By utilization of the Ladyzhenskaya inequality in -D, we have I 4 L 4 (Γ i) c 0 I 4 / L (Γ I i) 4 / L (Γ. i) It suffices to prove I 4 L (Γ i) c 0 A,

24 0 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE for some constant c 0. We compute, for i =,, I 4 x i = ( 0 [ 0 h ( ) + α x dz) ( ) ( ) ( ) ] + α x + α x i x dz + + α x h. x i Γb h ( aking into account the boundary condition ( ) 0 ( ) + α x = + α h x x α Γb + α x x α From(4.50)(4.5), we have for i =,, ( I4 x i ) (x, x ) c 0 ( x i which implies (4.48). We now prove (4.49); by (4.) We deduce from (4.7) L ( h,0) ( ) + α x L ( h,0) ( 0 I 4 (x, x ) 4 4 dz h I 4 L4 () L6 () ) Γi + α = 0, we obtain ( ) + α x dz ( (4.50) ) L + α x. ( h,0) (4.5) + ) c 0 4 ( 6 A ) 4 ( 6 A ) 4 c 0 4 A c 0 A A,. ( ) ) + α x, L ( h,0) so (4.49) is proved. We obtain (4.47) by adding (4.48) and (4.49) together, and the lemma is proved. Now we turn to the a priori estimates for the salinity S. Multiply (4.5) by A S, then integrate both sides over, to obtain d dt A S + A S + (v ) A S dx + w(v) S A S dx + (v ) A S dx + w(v ) S A S dx = 0. (4.5)

25 PRMIIVE EQUAIONS IN HIN DOMAINS By the Agmon inequality, we have (v )S A S dx v L S A S c 0 A v S A S 8 A S + c 0 A S A v. w(v) S A S dx M w(v) S A S dx ( 0 ) (divv) dz S h A S dx ( 0 ) ( (divv) 0 ( S dz Γ i h h }{{} =I, L 4 ( 0 (A S) dz h ) } {{ } L ) dz) }{{} =I 5, L 4 dx dx I L 4 I 5 L 4 A S, (by Lemma 4.6 below) c 0 v A v S A S 8 A S + c 0 v A v A S. (4.5) (4.54) (v )S A S dx v L S L A S v L S A S 8 A S + c 0 S [ A v + g v + g v ]. (4.55) Applying the same technique as in (4.9), (4.40) and (4.54), we find w(v ) S A S dx 8 A S + c 0 v A v A S. (4.56) From (4.5)-(4.56), one is led to d dt A S + A S c 0 A S [ A v + A v + g v + g v ] + c 0 A S [ v A v + v A v (4.57) ] ψ (t) A S, with ϕ (t) = c 0 [ A v + A v + g v + g v + v A v + v A v ].

26 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE he following lemma yields the desired estimate for I 5 = ) ) S dz appearing in (4.54). ( 0 h ( Lemma 4.6. For S D(A ), S I 5 L 4 (Γ i) c 0 A S. (4.58) Proof. his lemma can be proved in the same fashion as Lemma 4.5 by combining the Ladyzhenskaya inequality in -D and the boundary condition for S, we will just give a sketch of the proof. For i =,, [ S ( ) S I 5 x i = ) ( S x i dz + ( S ) S By taking advantage of the boundary condition for S, Γi = 0: ( ) S 0 ( ) ( S ) S dz. Γb h x Γ b ( h ) ]. x α he remaining part can be completed by Hölder s inequality and Ladyzhenskaya s inequality. So (4.58) is justified. We summarize all estimates for v, S, into the following proposition. Proposition 4.7. Assume u = (v,, S) is the solution of (4.)-(4.5) with initial data u 0 = (0, 0, S 0 ) V. hen we have the energy inequalities: d ] dt A v + [ c 0 v A v f v + c 0 ( + S ) (4.59) + c 0 [ A v + ( g v + g v ) + f ] Ro + v A v v d ] dt A + [ c 0 v A f + Q + (4.60) + c 0 [ v A v A + A v ] + + c 0 [ v A v + A v + v A v ] + c 0 [ gv + g v ]. d dt A S + A S c 0 A S [ A v + A v ( g v + g v ) ] (4.6) + c 0 A S [ v A v + v A v ] We conclude this section with a remark. Remark 4.4 Assume c 0 v is small. he second terms in (4.59) and (4.60) will be positive, the right hand side of (4.59) can be controlled using the inequalities valid for weak solutions and if f v and g v are small, say of order, then, the right hand side will be of order and therefore v will stay of order and the estimates can be closed to obtain the global boundedness of v and the integrability (in time) of A v. hen, without any further assumptions on and S, their H norms will stay bounded for all time.

27 PRMIIVE EQUAIONS IN HIN DOMAINS 5. Existence of global in time strong solutions. In this section we will prove the existence of global in time strong solutions to the PEs under small data depending on the thickness. It is noteworthy that the local/global strong solution for the PEs is so far an open problem due to the lack of the estimates on the vertical component w. his work will give a partial answer to this question when the domain is thin. Recall that R 0 () is a monotone positive function satisfying o proceed, observe first lim 0 R 0() = 0. (5.) Also we assume A u0 + Q + f v + g v + g v + f R 0(). (5.) R () = c 0 v + c 0 + c 0 v 4 v 4, + c 0 4 v v, + c 0 4 4, (5.) 4. Notice that all terms in R () appear in the left-hand side of (4.59)-(4.6). Fixed σ > to be specified later on (σ > e 4 ), then take such that for 0 <, 0 σr 0() 4. (5.4) Set u = (v,, 0). hen the weak formulation for (4.)-(4.6) can be rewritten as follows: o find u L (0, t ; V ) L (0, t ; D(A)) such that d dt (u, u ) + a(u + u, u ) + b(u + u, u + u, u ) +e(v, v ) = (Q, ), u V, u(0, x) = u 0 (x) V. We are in a position to state the main theorem in this section. heorem 5.. Assume that (5.)- (5.4) hold. hen there exists a unique strong solution u to the PEs such that u L (0, ; V ) L ((0, t ; D(A)), u t L (0, t ; H), Furthermore the solution u satisfies the following estimate: A u(t) σr 0(), t 0. Proof. We will prove this theorem by using the Galerkin method. Since this method is standard, we will just present an outline of the proof. We divide the proof into 4 steps. Step. Galerkin approximate solution: Let e j = (e v j, e j, es j ), j =,,, be the eigenvectors of the principle operator A corresponding to the eigenvalues λ j, j =,,, such that Ae j = λ j e j, λ j.

28 4 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE hen the {e j } j= define a complete orthogonal basis in V. Furthermore, let P m the orthogonal projector on the first eigenvalues e,..., e m. Now we look for approximate solutions u m = m j= g jme j which satisfy, with u m = (v m, m, S m ): d dt (u m, e j ) + a(u m + u, e j ) + b(u m + u, u m + u, e j ) +e(v m, e v j ) = (Q, e j ), for j =,,, m, (u m (0), e j ) = (u 0, e j ) for j =,,, m. Solving for u m amounts to solving a system of differential equations for the g jm, j =,,, m. It is easy to see that the above a priori estimates apply as well to u m, namely {u m } is a sequence bounded in L (0, ; V ) L (0, t ; D(A)). We are then able to pass to the limit using the standard compactness methods using some additional estimate on (u m ) t. Step. Global solution for the Galerkin Approximation. By the general theory of ODE, the system of differential equations for g jm admits continuous solutions. For σ >, there exists t σ > 0 such that A um (t) σr 0(), t [0, t σ ]. (5.5) We claim that t σ =. Assume the contrary, i. e. t σ <, then let t = t σ be the first time such that A um (t σ ) = σr 0(). (5.6) For any t [0, t σ ], We rewrite (4.59)-(4.6) as follows d dt A v m + A v m f v + 4 σr 0() + ψ (t) A v m, (5.7) d dt A m + A m A v m + f + Q + ψ (t) A m, (5.8) d dt A S m + A S m A S m A v m. (5.9) we have, using the estimates on the linear equation satisfied by v and, for small enough: t 0 ψ (s) ds c 0 tr 0+R 0, t 0 ψ (s) ds c 0 tr 0+R 0, t 0 ψ (s) ds c 0 tr 0+R 0. First we add (5.7) and (5.8) and then use Poincaré s inequalities (4.6) to obtain d dt ( A v m + A m ) + c 0 ( A v m + A m ) f v + f + Q + 4 σr 0() + ψ A v m + ψ A m. aking into consideration the zero initial data for v and, the Gronwall inequality yields, for t [0, t σ ) and small so that R0 c0 A v m + A m t 0 t ( A v m + A m )dt 0 [ f v + f ] ds + Q ) + c 0 t 0 4 σr 0(), (5.0) ( f v + f ) ds + c 0 R 0t. (5.)

29 PRMIIVE EQUAIONS IN HIN DOMAINS 5 For the salinity S m, by using the Poincaré s inequality (4.0), we have from (5.9) d dt A S m ( 4 A v m ) A S m. he Gronwall inequality yields A S m A S 0 e t 0 ( 4 Avm ) dτ herefore from (5.0) and (5.), we obtain A S 0 e A v 0 e ( f v 4 σr 0 ())t. A v m + A m + A S m A v 0 + A 0 + A S 0 e A v 0 hen at t = t σ (5.6) gives σr 0() e 4 R 0 () + c 0 + c 0 ( f v + f + Q ) + c 0 (5.) 4 σr 0(). 4 σr 0(). (5.) We take σ > e 4, and a contradiction occurs, so t σ =. Step. Passage to limit and the existence of solution: o pass to the limit, we need to show that {(u m ) t } bounded in L (0, t ; H), t > 0, which can be proved by repeating some of the estimates in the previous section, we omit them. So far we infer that there exists u m L (0, ; V ) L (0, t ; D(A)) and subsequence of u m (still denoted by u m ), such that u m u weakly in L (0, t ; V ), (5.4) u m u weakly in L (0, t ; D(A)), (5.5) (u m ) t u t weakly in L (0, t ; H). (5.6) By Aubin compactness lemma, we have u m u strongly in L (0, t ; V ). Since the PEs bear much similarity with the Navier Stokes equations, to prove that u is a local strong solution of the PEs, only the vertical convection term needs more care, which is presented below.. Lemma 4.. Assume u m u in the sense of (5.4)-(5.6), for any vector function v with components in C((0, t ) ) and scalar function Q in C((0, t ) ), then we have t w(v m ) t m v dxdt w(v) v dxdt, 0 0 t w(v m ) t m Q dxdt w(v) Q dxdt, 0 0 t w(v m ) S t m Q dxdt w(v) S Q dxdt. 0 0 Step 4. Uniqueness of the strong solution: o finish the proof, it remains to confirm the uniqueness of the solution obtained. Recall that the uniqueness of strong solutions of the -D Navier Stokes equation is well known; the main difference of the PEs from D NS equations is the term w(v). In D NS equations, w H, however all we obtain for the PEs is w L. Let u = (v,, S ), u =

30 6 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE (v,, S ) be two solutions of (.)-(.8) satisfying the non-homogeneous boundary condition (4.9). Denote δu = (δv, δ, δs) = (v v,, S S ), δρ = ρ ρ, δp = p p. hus δu satisfies the following equations: (δv) t + (v )(δv) + w(v) (δv) + ((δv) )v (5.7) + w(δv) + f Ro k (δv) + grad(δp) s 0 + b grad(δρ) dz = (δv) + (δv) x Re Re x, (δ ) + (v )(δ ) + w(v) (δ ) + ((δv) ) (5.8) t + w(δv) = (δ ) + (δ ), Rt Rt (δs) t x + (v )(δs) + w(v) (δs) + ((δv) )S (5.9) + w(δv) S = (δs) + (δs) Rs Rs x, 0 div (δv) dz = 0, (5.0) h δρ = ρ 0 ( β (δ ) + β S (δs)). (5.) supplemented with the homogeneous boundary condition: (δv) + α v (δv) = 0, (δ ) (δv) = 0, n (δ ) + α (δ ) = 0, = (δs) n = 0 on Γ b Γ l. (δs) n = 0 on Γ i, (5.) ake the inner product in H on both side of (5.7) with δv, noticing the skewsymmetry identity for the trilinear term (.), we have d dt δv + δv + α δv dx dx + ((δv) )v (δv) Γ i {0} + w(δ) 0 (δv)dx + (δv) grad ( β x (δ ) + β S (δs)) dx = 0. x Now we make some brief estimates. ((δv) )v (δv) δv L 4 v 4 δv + C δv v 4.

31 PRMIIVE EQUAIONS IN HIN DOMAINS 7 and and w(δv) (δv)dx (δv) grad 0 Γ i ( 0 ) (divv) dz } h {{ } L ( 0 ) (δv) dz h }{{} L 4 ( 0 ( ) dz) h }{{} =I, L 4 dx dx δv A v δv δv 4 δv + C A v δv. x ( β (δ ) + β S (δs)) dx δv ( δ + δs ) 4 ( δv + δ + δs ). we arrive at d dt δv + δv + α v δv dx dx Γ i ) C ( v 4 + A v δv + 4 ( δ + δs ). Similar estimates can be carried over to and S, we have the following d dt δ + δ + α δ dx dx Γ ( i ) 4 v + C A + 4 δ. and d dt δs + δs (5.) (5.4) ( ) S 4 δv + C A S + S 4 δ. (5.5) Adding (5.)-(5.5), we have ( d dt u + u C v 4 + A v + A + 4 ) + S A S + S 4 u = G(t) u, where G(t) = v 4 + A v + A S A S + S 4.

32 8 CHANGBING HU, ROGER EMAM, MOHAMMED ZIANE By estimates (5.5), we have t 0 G(τ)dτ <. hen Gronwall s inequality yields u = 0, so the uniqueness of the solution (.)- (.8) is confirmed. Finally the theorem is proved. A direct consequence from the process of above proof is given below. Corollary 5.. Assume u = (v,, S) is the solution obtained in the previous theorem, then there exists a constant λ such that S(t) S 0 e λt, which implies S behaves trivially in the long time. 6. Proof of the main theorem. In this part we will prove the main theorem stated in Section. Basically the core part of the proof has been done in the previous section. Recall that we obtain the existence of the solution of the PEs under assumptions in terms of u,, in this part we will transform the results into the physical boundary condition, i. e., g v, g. oward this end, we first present a lemma which gives estimates on v, in term of g v, g. Lemma 6.. Let v, be solutions of (4.) (4.), assume g v H (Γ i ), g H (Γ i ), then (v, ) L (0, t ; H ( )) L (0.t ; H ( )). (6.) Moreover (v, ) satisfies the following estimates: v (t) + v L (Γ i)(t) c 0 g v L (Γ i), (6.) (t) + L (Γ i) (t) c 0 g L (Γ i), (6.) t 0 t 0 v, ds c 0 t gv L (Γ i) + c 0 t 4 gv L (Γ i), (6.4), ds c 0 t g L (Γ i) + c 0 t g L (Γ i). (6.5) Proof. Inequalities (6.) and (6.) can be obtained from the weak formulation of (4.) (4.) and Lax-Milgram theorem. Inequality (6.4) is the main result in heorem. in [6]. Inequality (6.5) is just the Grisvard-Iooss inequality for. Lemma 6.. f v and f defined in (4.7) and (4.8) satisfy f v c 0 4 v 5 4 v, ( v ) 4 x + v, (6.6) ( + c 0 4 v 4 v, ) + v 4, + c0 v + c 0, f c 0 v 4 4, + c 0 v v 4, 4,. (6.7)