WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE

WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper focuses on the two-dimensional Benjamin-Bona-Mahony and Benjamin-Bona-Mahony-Burgers equations with a general flux function. The aim is at the global (in time) well-posedness of the initial-and boundary-value problem for these equations defined in the upper half-plane. Under suitable growth conditions on the flux function, we are able to establish the global well-posedness in a Sobolev class. When the initial- and boundary-data become more regular, the corresponding solutions are shown to be classical. In addition, the continuous dependence on the data is also obtained.. Introduction This paper is concerned with the two-dimensional (2D) Benjamin-Bona-Mahony- Burgers equation of the form u t + div (φ(u)) = ν u + ν 2 u t in Ω (, T ), (.a) u = g on Ω {t = }, (.b) u = h on Ω (, T ), (.c) where Ω R 2 denotes a smooth domain, u = u(x, t) is a scalar function, φ(u) is a vector-valued flux function, and ν and ν 2 > are real parameters. (.a) is a natural generalization of the D Benjamin-Bona-Mahony equation u t + u x + uu x u xxt =, (.2) which governs the unidirectional propagation of D long waves with small amplitudes. Therefore (.a) with ν = is sometimes called the generalized Benjamin-Bona- Mahony (GBBM) equation while (.a) with ν > the Benjamin-Bona-Mahony- Burgers (GBBM-B) equation. (.a) is also a regularization of the scalar conservation law u t + div(φ(u)) =. (.3) 2 Mathematics Subject Classification. 35B45, 35B65, 76B5, 76S5. Key words and phrases. Generalized Benjamin-Bona-Mahony equation, Benjamin-Bona-Mahony- Burgers equation, Buckley-Leverett equation, initial-boundary value problem, global well-posedness.

2 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN In addition, (.a) has also been derived to model the two-phase fluid flow in a porous medium, as in the oil recovery. In fact, (.a) is a special case of the well-known Buckley-Leverett equation u t + div(φ(u)) = div{h(u) (J(u) τu t )}, (.4) where u denotes the saturation of water, the functions φ, H and J are related to the capillary pressure and the permeability of water and oil [3]. (.a) follows from (.4) by linearizing the static capillary pressure J(u) and H around a constant state. Attention here will be focused on the case when Ω =, the upper half-plane. The aim is at the global well-posedness of (.) with inhomogeneous boundary data, namely h. One motivation behind this study is to rigorously validate the laboratory experiments involving water waves generated by a wavemaker mounted at the end of a water channel. We are able to prove the global existence and uniqueness of the mild and classical solutions to (.). In addition, a continuous dependence result is also obtained. Our main theorems can be stated as follows. Theorem. (Existence and uniqueness). Let ν, ν 2 >, and Ω =. Let T >. Suppose that (g, h) H 2 () C ([, T ]; H 2 (R)), and the flux φ C 2 (R, R 2 ) satisfies the conditions φ() = φ () = and φ L (R) C. (.5) Then (.) admits a unique mild solution u C([, T ]; H 2 ()). If we further assume that (g, h) C 2,α loc (R2 +) C ([, T ]; C 2,α loc (R)) for some < α <, then the mild solution is in fact a classical solution. Theorem.2 (Continuous dependence on data). Let ν, ν 2 >, and Ω =. Suppose that φ C 2 (R, R 2 ) satisfies the condition (.5). Then the mild solution obtained in Theorem. depends continuously on the initial datum g and boundary datum h. If we further assume that φ C 3 (R, R 2 ) satisfies φ L (R) C, (.6) then the same result also holds for the classical solution. It is worth remarking that Theorems. and.2 hold with either ν = or with ν > and do not rely on the regularization due to the Burgers dissipation. We briefly review related well-posedness results and then explain the main difficulties in proving Theorems. and.2. There is a very large literature on the global well-posedness for the D BBM equation on the whole line (see, e.g. [3, 4]). Extensive results

2D GBBM EQUATION 3 have also been obtained on the global well-posedness for the initial- and boundaryvalue problem of the BBM equation posed on the half-line (see, e.g. [2, 5, 6, 7, 8, 9,,, 2, 7, 9]). In particular, in the well-known articles [2, 3], the existence of classical solutions and their continuous dependence on the specified data were investigated. While current analytic results for the multi-dimensional BBM or BBM- Burgers equations only dealt with the existence of mild solutions on either the whole space or bounded domains with homogeneous boundary data (see, e.g. [, 5, 8]). The results presented here allow inhomogeneous boundary data, which correspond to the setup of a wavemaker mounted at the end of a channel in laboratory experiments. We emphasize that, our methods are also suitable for the corresponding initial value problem. Therefore, Theorems. and.2 are the complete extensions of the results in [2, 3] to the multi-dimensional case. The main difficulty in proving Theorem. is from two sources. First, the Green function for operator I in 2D is much more singular than the D case; and second, the inhomogeneous boundary data prevents us from obtaining a time-independent H upper bound, which very much simplifies the process of global-in-time estimates. To overcome the difficulties, we introduce a new function that assumes the homogeneous boundary data and rewrite the equation in an integral form through the Green function of the elliptic operator. In addition, we use the bootstrapping technique to obtain the classical solution of (.) instead of looking for the solution in classical spaces directly. The rest of this paper is divided into six sections. The first five sections deal with the case when ν = while the last section explains why the results for ν = can be extended to the case when ν >. Section 2 introduces a new function that assumes homogenous boundary data and converts (.) into an integral formulation in terms of this new function. Section 3 presents preliminary regularity estimates for the operator (I ) in the Sobolev space H 2 and in Hölder spaces. Section 4 proves that (.) has a unique local (in time) classical solution. We make use of the integral representation (2.5). Section 5 establishes the global existence and uniqueness of the local solution obtained in Section 4 by showing global bounds for the solution in H 2 and in Hölder spaces. Section 6 contains the continuous dependence results. The continuous dependence of the solution on the initial data and the boundary data is proven in two functional settings and the proof is lengthy. As aforementioned, Section 7 is devoted to the case when ν >.

4 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN 2. An alternative formulation In this section we set ν =. The case when ν > is handled in Section 7. This section provides an integral formulation of (.). In order to apply the standard elliptic theory in the functional framework of Sobolev spaces, we shall rewrite equation (.) with homogeneous boundary data. This is achieved by setting v(x, t) = u(x, t) h(x, t)e x 2, which satisfies (I )v t + div ( φ(v + he x 2 ) ) = h x x te x 2 in (, T ), (2.a) v = g on {t = }, (2.b) v = on (, T ). (2.c) where g(x) = g(x) h(x, )e x 2. (2.2) Denoting (I ) f as the unique solution to the elliptic equation (I )u = f in Ω, (2.3a) u = on Ω, (2.3b) we can formally write the solution v of (2.) via the integral representation v(x, t) = g(x) + (I ) ( {h x x (x, t) h x x (x, )}e x ) 2 t For short, we rewrite (2.4) as the form (I ) { div ( φ(v + he x 2 ) )} dτ. (2.4) v = Av = g + Bh + Cv, (2.5) where, for x and t, Bh(x, t) := (I ) ( {h x x (x, t) h x x (x, )}e x 2 ), Cv(x, t) := t (I ) { div ( φ(v + he x 2 ) )} dτ. 3. Preliminary results This section specifies the functional spaces and provides two preliminary estimates on the solutions to the elliptic equation (2.3). In the rest of this paper, we write C k ([, T ]; H l ()) for the space { u : [, T ] H l (R 2 k +) lim k u t t t k (t) k u } t k (t ) = t H l ( ) [, T ] j=

equipped with norm u C k t Hx l := max 2D GBBM EQUATION 5 t [,T ] j= k k u t k (t) H. l ( ) The spaces with the particular indices k =, and l =, 2 will be frequently used. For the simplicity of notation, when k =, we omit the super-index, that is, C([, T ]; H l ()) C ([, T ]; H l ()) and CtH x l Ct. We will also need Hl x the space { } C([, T ]; L p ()) u : [, T ] L p () lim u(t) u(t ) t t L p ( ) = equipped with norm u CtL p x = max t [,T ] u(t) L p ( ). Similar notation is used to define the space of the boundary data which is only defined on the real line R. We introduce { } C([, T ]; H 2 (R)) h : [, T ] H 2 (R) lim h(t) h(t ) H 2 (R) = t t equipped with norm h CtH 2 x = max t [,T ] h(t) H 2 (R). To study the classical solutions, we let C k,α (Ω) denote the space of k-times classically differentiable functions whose k-th derivatives are Hölder continuous with exponent α. The norm on C k,α (Ω) is given by u C k,α (Ω) = k sup j= x Ω where D j u denotes the j-th classical derivative of u. D j D k u(x) D k u(y) u(x) + sup x,y Ω x y α, To deal with the integral representation (2.4), we need some crucial estimates on the operator (I ). In particular, the bounds in the following propositions will be employed in the subsequent sections. Proposition 3.. Assume f L 2 (Ω). Then the Dirichlet problem (2.3) admits a unique solution u H 2 (Ω). Furthermore, u H 2 (Ω) C f L 2 (Ω). If f is instead in a Hölder space, then we have the following Hölder s estimates for the solution of (2.3).

6 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN Proposition 3.2. Assume that Ω R 2 is a smooth domain. Assume that f is in C,α loc (Ω) L2 (Ω) for some < α <. Then the solution u of the Dirichlet problem (2.3) lies in C 2,α loc (Ω) H2 (Ω). Furthermore, for any compact subsets Ω and Ω 2 of Ω with Ω 2 Ω, u C 2,α (Ω 2 ) C( f C,α (Ω ) + f L 2 (Ω )), where C > depends only on the distance between Ω 2 and Ω. Proof. By Proposition 3., we have u H 2 (Ω) and u H 2 (Ω) C f L 2 (Ω). Sobolev embedding theorem says that u C,α loc (Ω) and u C,α (Ω ) C u H 2 (Ω ) for any compact subset Ω. Thus, we get that u = u f C,α loc (Ω). It follows from Lemma 4.2 and Theorem 4.6 in [4] that u C 2,α loc (Ω) and, for any compact subsets Ω and Ω 2 of Ω with Ω 2 Ω, we have u C 2,α (Ω 2 ) C( f C,α (Ω ) + u C,α (Ω )) C( f C,α (Ω ) + u H 2 (Ω )) C( f C,α (Ω ) + f L 2 (Ω )), where C > depends only on the distance between Ω 2 and Ω. 4. Local-in-time existence This section proves that (.) has a unique local (in time) classical solution. We make use of the integral representation (2.5). Due to the difficulty of applying the contraction mapping principle in the setting of Hölder spaces, the proof is divided into two steps. The first step applies the contraction mapping principle to (2.5) in the setting of Sobolev spaces to obtain a unique local solution. The second step obtains the desired regularity of the local solution through a bootstrapping procedure. Lemma 4.. Let (g, h) H 2 () C([, T ]; H 2 (R)), and φ satisfy the condition (.5). Then there is S with < S T, depending only on g and h, such that (2.5) has a unique solution v C([, S]; H 2 ()). Proof. This local existence and uniqueness result is proven through the contraction mapping principle. More precisely, we show that A defined in (2.5) is a contraction map from B(, R) C([, S]; H 2 ()) to itself, where B(, R) denotes the closed ball centered at with radius R in C([, S]; H 2 ()). S and R will be specified later in the proof. It follows from (.5), (2.5), Proposition 3. and the mean value theorem

2D GBBM EQUATION 7 that Av CtH 2 x g H 2 + C h C th 2 x + CS div(φ(v + he x 2 )) CtL 2 x g H 2 + C h CtH 2 x + CS φ (v + he x 2 ) CtL x (v + he x 2 ) CtL 2 x (4.) and C + C 2 S( + R)R Av Aw 2 C thx 2 CS2 div(φ(v + he x 2 )) div(φ(w + he x 2 )) 2 C tl 2 x 2 ( CS 2 φ j(v + he x 2 )(v w) xj 2 j= + ( φ j(v + he x 2 ) φ (w + he x 2 ) ) (w + he x 2 ) xj 2 ) dx (4.2) CS 2 { φ j(v + he x 2 ) 2 C tl x v w 2 C th x } + φ (v + he x 2 ) 2 C tl (w + he x 2 ) 2 x C v tl 4 w 2 x C tl 4 x C 2 S 2 ( + R) 2 v w 2 C, thx 2 where v lies between the line segment joining v and w, C is a constant depending on g H 2 and h CtH x 2, and C 2 is a constant depending on h CtH x 2. Note that (4.2) implies A is a continuous map of C([, S]; H 2 ()) to itself. According to (4.), A maps B(, R) onto itself if R C + C 2 S( + R)R. (4.3) Hence, by (4.2), A is a contraction mapping of this ball if C 2 S( + R) <. These conditions will be met if we take R = 2C and find a positive value S > small enough such that C 2 S( + 2C ) 2. (4.4) Now, let v (x, t) = g(x) + Bh(x, t) and v n (x, t) = Av n (x, t) = v (x, t) + Cv n (x, t) for n. The contraction mapping principle gives that the sequence v n (x, t) converges in C([, S]; H 2 ()) to the unique solution v of (2.5) in the ball v CtH x 2 R. If the initial data g and the boundary data h are also Hölder, then the corresponding solution can also be shown to be Hölder. This is achieved through the Sobolev embeddings and a bootstrapping procedure.

8 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN Lemma 4.2. Assume g C 2,α loc (R2 +) H 2 () and h C ([, T ]; C 2,α loc (R) H2 (R)) for some < α <. Let φ satisfy the condition (.5). Then any solution v C([, T ]; H 2 ()) of (2.5) actually belongs to C([, T ]; H 2 ()) C ([, T ]; C 2,α loc (R2 +)). Proof. Since v C([, T ]; H 2 ()) and h C ([, T ]; H 2 (R)), we have Proposition 3. gives div(φ(v + he x 2 )) C([, T ]; H ()). (I ) {div(φ(v + he x 2 ))} C([, T ]; H 3 ()) C([, T ]; C,α ()). (4.5) On the other hand, h C ([, T ]; C 2,α loc (R)) implies Proposition 3.2 then yields In view of (4.5) and (4.6), we obtain h x x te x 2 C([, T ]; C,α loc (R2 +)). (I ) {h x x te x 2 } C([, T ]; C 2,α loc (R2 +)). (4.6) v t = (I ) {h x x te x 2 div(φ(v + he x 2 ))} C([, T ]; C,α loc (R2 +)), which implies that v(x, t) = g(x) + t v τ (x, τ)dτ = g(x) h(x, )e x 2 + t Using (4.7) and Proposition 3.2, we have v τ (x, τ)dτ C ([, T ]; C,α loc (R2 +)). (4.7) v t = (I ) {h x x te x 2 div(φ(v + he x 2 ))} C([, T ]; C 2,α loc (R2 +)), and hence v C ([, T ]; C 2,α loc (R2 +)). 5. Global-in-time existence This section shows that the local (in time) solution obtained in the previous section can be extended into a global one. This is achieved by establishing a global bound for v(t) H 2 under the condition that the flux φ obeys suitable growth condition. We start with a global H -bound.

2D GBBM EQUATION 9 Lemma 5.. Suppose (g, h) H 2 () C ([, S]; H 2 (R)), and φ satisfy the condition (.5). Then the solution v of (2.5) obtained in Lemma 4. satisfies the estimates and v H v 2 H g 2 H + CS h 2 C t H2 x ( + h C t H 2 x ) + C( + h C t H 2 x ) t v 2 H ds (5.) [ g 2 H + CS h 2 C t H2 x ( + h C t H 2 x )] /2e CS(+ h C t H 2 x ), (5.2) where C > is a constant depending only on φ. Proof. Multiplying (2.a) by v and integrating over, we get d (v 2 + v 2 ) dx = {h x x 2 dt te x 2 div(φ(v + he x 2 ))}v dx = h x x te x 2 v dx div(vφ(v + he x 2 )) dx + φ(v + he x 2 ) v dx. Since v = on, div(vφ(v + he x 2 )) dx = Now let Φ C (R; R 2 ) satisfy Φ = φ and Φ() =. Then φ(v + he x 2 ) v dx = = = φ(v + he x 2 ) (v + he x 2 ) dx div(φ(v + he x 2 )) dx (5.3) vφ(v + he x 2 ) n ds =. (5.4) φ(v + he x 2 ) (he x 2 ) dx φ(v + he x 2 ) (he x 2 ) dx Φ(h) n ds φ(v + he x 2 ) (he x 2 ) dx. By the Sobolev embedding H 2 (R) L q (R) for all q 2, h C ([, S]; H 2 (R)) C ([, S]; L q (R)). (5.5) Employing the mean value theorem together with (.5) and the properties of Φ, we obtain Hence Φ(h) n ds = Φ(h) Φ() C h 3. Φ(h) Φ() ds C h 3 C tl 3 x C h 3 C t H2 x. (5.6)

Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN Applying the mean value theorem and (.5) again, we have φ(v + he x 2 ) C v + he x 2 2. As a consequence, φ(v + he x 2 ) (he x 2 ) dx v + he x 2 2 L h 4 H C( v 2 L 4 + h 2 L 4 x ) h H x (5.7) C( v 2 H + h 2 C t H2 x ) h C t H 2 x. From (5.3)-(5.7), we can conclude that d dt v 2 H 2 h x x te x 2 v dx + C( v 2 H + h 2 C ) h t H2 x C t Hx 2 that is, 2 h C t Hx 2 v L 2 + C( v 2 H + h 2 C ) h t H2 x C t Hx 2 C( v 2 H + h 2 C )( + h t H2 x C t Hx 2 ), v 2 H g 2 H + CS h 2 C ( + h t H2 x C t Hx 2 ) + C( + h C t Hx 2 ) where C depends only on φ. Gronwall s inequality gives t v 2 H ds, v H { g 2 H + CS h 2 C t H2 x ( + h C t H 2 x )} /2 e CS(+ h C t H2 x ). Now we derive the H 2 -estimates based on the H -estimates we just obtained. Lemma 5.2. Suppose (g, h) H 2 () h C ([, S]; H 2 (R)), and φ satisfy the condition (.5). estimate Then the solution v of (2.5) obtained in Lemma 4. satisfies the { v H 2 C( + S) /2 exp CS( + S) /2 e CS}, where C > is a constant depending only on g, h and φ. Proof. Multiplying (2.a) by u and then integrating on, we have d ( v 2 + v 2 )dx 2 dt = h x x te x 2 v dx div(φ(v + he x 2 )) v dx. By the mean value theorem and (.5), div(φ(v + he x 2 )) φ (v + he x 2 ) (v + he x 2 ) C v + he x 2 (v + he x 2 ). (5.8) (5.9)

Thus, Hölder s inequality gives d ( v 2 + v 2 )dx 2 h x x dt t L 2 x v L 2 2D GBBM EQUATION + C v + he x 2 L 4 (v + he x 2 ) L 4 v L 2. By the Sobolev embedding H () L 4 () and Young s inequality, d dt ( v 2 L 2 + v 2 L 2 ) C h 4 C t H2 x + C( + v CtH x + h C t H2 x ) v 2 H 2, where C > depends only on φ; that is, v 2 L 2 + v 2 L 2 C( g 2 H + S h 4 2 C ) + C( + v t H2 x CtH + h x C ) t H2 x Combining (5.) and (5.), we obtain v 2 H 2 C{ g 2 H 2 + S( + h 4 C t H2 x )} + C( + v CtH x + h Ct ) H2 x Applying (5.2) to (5.), we have t t v 2 H 2 ds. t v 2 H C( + S) + C( + S) /2 e CS v 2 2 H ds, 2 v 2 H 2 ds. where C depends only on g, h, and φ. Therefore, by Gronwall s inequality, v H 2 C( + S) /2 exp {CS( + S) /2 e CS} which concludes the proof of the lemma. (5.) (5.) 6. Continuous dependence of the solution on data This section is devoted to proving Theorem.2. That is, we establish the desired continuous dependence. For the sake of clarity, we will divide the rest of this section into two subsections. The first subsection proves the continuous dependence in the regularity setting of H 2 while the second subsection focuses on the continuous dependence in the intersection space of H 2 and a Hölder class. The precise statements can be found in the lemmas below. 6.. Continuous dependence in H 2. Let L m denote the mapping that takes the data g and h to the corresponding solutions of (.). By Theorem. we have L m : X m = H 2 () C ([, T ]; H 2 (R)) C([, T ]; H 2 ()). Since H 2 () and C ([, T ]; H 2 (R)) are Banach spaces, the space X m equipped with the usual product topology is also a Banach space.

2 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN Lemma 6.. Suppose that φ C 2 (R, R 2 ) satisfies the condition (.5). Then L m is continuous. Proof. Let (g i, h i ) X m and u i = L m (g i, h i ) be the mild solution of (.) corresponding to the initial data g i and the boundary data h i, i =, 2. Set v i = u i h i e x 2, i =, 2. Then v i satisfies the following initial-boundary value problem: (v i ) t (v i ) t (h i ) x x te x 2 + div (φ(v i + h i e x 2 )) =, v i (x, ) = g i (x) h i (x, )e x 2 := g i (x), ( v i (x, ), t ) =. Define w = v v 2. Then w satisfies: { w t w t h x x te x 2 + div (φ(v + h e x 2 )) div (φ(v 2 + h 2 e x 2 )) =, w(x, ) = g(x), w ( (x, ), t ) =, (6.) where g = g g 2 and h = h h 2. In addition, we derive that w satisfies the following integral equation: w(x, t) = g(x) + (I ) ( {h x x (x, t) h x x (x, )}e x 2 ) t (I ) { div(φ(v + h e x 2 )) } dτ. (6.2) Given ε >. Suppose that the distance between (g, h ) and (g 2, h 2 ) in X m is small enough such that (a) g H 2 ε, (b) h C t H 2 x ε. Taking H 2 norm on both sides of (6.2) and using Proposition 3., we derive Since w H 2 g H 2 + h C t H 2 x + C t div(φ(v + h e x 2 )) L 2 dτ. div(φ(v + h e x 2 )) = φ (v + h e x 2 ) (v + h e x 2 ) φ (v 2 + h 2 e x 2 ) (v 2 + h 2 e x 2 ) (6.3) = {φ (v + h e x 2 ) φ (v 2 + h 2 e x 2 )} (v + h e x 2 ) (6.4) + φ (v 2 + h 2 e x 2 ) (w + he x 2 ), the mean value theorem and condition (.5) yield div(φ(v + h e x 2 )) C{ w + he x 2 (v + h e x 2 ) + v 2 + h 2 e x 2 (w + he x 2 ) }. (6.5)

Applying (6.5) and Hölder s inequality to (6.3), we obtain w H 2 g H 2 + h C t H 2 x + C + C t 2D GBBM EQUATION 3 t w + he x 2 L 4 (v + h e x 2 ) L 4dτ v 2 + h 2 e x 2 L 4 (w + he x 2 ) L 4dτ. By Sobolev s inequality, w H 2 g H 2 + h C t Hx 2 + C t g H 2 + C( + T ) h C t H 2 x + C w + he x 2 H 2dτ t w H 2dτ, (6.6) where C depends only on v, v 2, h, h 2, and φ. Then Gronwall s inequality gives w CtH x 2 { g H 2 + h Ct }e CT. (6.7) H2 x Note that w = v v 2 = u u 2 (h h 2 )e x 2 = L m (g, h ) L m (g 2, h 2 ) he x 2. Therefore, by (6.7), L m (g, h ) L m (g 2, h 2 ) CtH 2 x w C th 2 x + h C th 2 x { g H 2 + h C t H 2 x }e CT e CT ε. 6.2. Continuous dependence in the intersection of H 2 and a Hölder space. This subsection proves the continuous dependence in the setting of the intersection of H 2 and a Höler space. First, we introduce the metrics on the spaces C 2,α loc (Ω) and C ([, T ]; C 2,α loc (Ω)), where Ω can be R2 + or R. Let {Ω i } i= be an increasing sequence of compact subsets of Ω satisfy (i) Ω i Ω i+ for all i N, (ii) Ω i = Ω. i= For a function f C ([, T ]; C 2,α loc (Ω)), we define ρ i (f) = f C ([,T ];C 2,α (Ω i )) for i N. Then {ρ i } forms a family of seminorms on C 2,α loc (Ω). For f, f 2 C ([, T ]; C 2,α loc (Ω)), we define d(f, f 2 ) = 2 i ρ i (f f 2 ) + ρ i (f f 2 ). i= Then d is a metric on C ([, T ]; C 2,α loc (Ω)). It is clear that f k f with respect to d if and only if ρ i (f k f) for all i. The topology induced by the metric d is

4 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN independent of the choice of the sequence {Ω i } i=. A metric on the space C2,α loc (Ω) can be defined similarly if we replace the seminorm above by for f C 2,α loc (Ω) and i N. ρ i (f) = f C 2,α (Ω i ) Now we fix a sequence {Ω i } i= of compact convex sets in R2 + such that the conditions (i) and (ii) hold. Let I i denote the projection of Ω i to x -axis. Then {I i } i= forms a sequence of compact sets in R satisfying (i) and (ii). As stated above, the sequences {Ω i } i= and {I i} i= induce metrics d on C 2,α loc (R2 +), d 2 on C ([, T ]; C 2,α loc (R)), and d 3 on C ([, T ]; C 2,α loc (R2 +)) respectively. In Theorem., we get that for a given pair (g, h) X c := [H 2 () C 2,α loc (R2 +)] C ([, T ]; H 2 (R) C 2,α loc (R)) of initial and a boundary data, then (.) admits a unique classical solution u Y := C([, T ]; H 2 ()) C ([, T ]; C 2,α loc (R2 +)). If we let L c denote the mapping that takes the pair (g, h) into the corresponding classical solution u, then L c : X c Y. (6.8) We remark that if (M, d M ) and (N, d N ) are two metric spaces, then the product space M N is a metric space with metric d M N (x, y) = d M (x, y) + d N (x, y) for x, y M N, and their intersection M N is a metric space with metric d M N (x, y) = d M (x, y) + d N (x, y) for x, y M N. We can immediately conclude that L c is a continuous mapping in (6.8) if we prove that i L c and j L c are both continuous where i and j are the natural inclusions of Y into C([, T ]; H 2 ()) and C ([, T ]; C 2,α loc (R2 +)) respectively. Lemma 6.2. Suppose that φ C 3 (R, R 2 ) satisfies the conditions (.5)-(.6). Then L c is continuous. Proof. By the discussions before the statement of this lemma, it suffices to prove L c : X c C([, T ]; H 2 ()) and L c : X c C ([, T ]; C 2,α loc (R2 +)) (6.9) are both continuous. Comparing the metrics of the spaces X m and X c, we can easily get the continuity of the mapping L c : X c C([, T ]; H 2 ()) from Lemma 6.. In this proof, we focus on showing that the second mapping of (6.9) is sequentially continuous.

2D GBBM EQUATION 5 Let (g i, h i ) X c and u i = L c (g i, h i ) be the classical solution of (.) corresponding to the initial data g i and the boundary data h i, i =, 2. Let w, g, and h be defined as in the proof of Lemma 6.. Then w satisfies the initial-boundary value problem (6.), the integral equation (6.2), and the estimate (6.7). Let Ω, Ω, and Ω be any given three compact convex sets in with Ω Ω Ω and let I and I be the projections of Ω and Ω to the x -axis respectively. First, we take C,α (Ω ) norm on both sides of (6.2) and use Sobolev s inequality to obtain w C,α (Ω ) g C,α (Ω ) + (I ) ( {h x x (x, t) h x x (x, )}e x ) 2 C 2,α (Ω ) t + (I ) { div(φ(v + h e x 2 )) } C,α ( ) dτ g C,α (Ω ) + (I ) ( {h x x (x, t) h x x (x, )}e x 2 ) C 2,α (Ω ) (6.) + t (I ) { div(φ(v + h e x 2 )) } H 3 ( ) dτ. Lemma (3.2) implies (I ) ( {h x x (x, t) h x x (x, )}e x 2 ) C 2,α (Ω ) { C {h x x (x, t) h x x (x, )}e x 2 C,α (Ω) + {h x x (x, t) h x x (x, )}e x } 2 L 2 (Ω) (6.) C ( h C 2,α (I) + h H 2 (I)), where C depends only on the distance between Ω and Ω. Employing Proposition 3. and applying mean value theorem together with (.5) to (6.4), we derive that (I ) { div(φ(v + h e x 2 )) } H 3 ( ) C div(φ(v + h e x 2 )) H ( { ) C w + he x 2 W,4 (R 2 ) (v + h + e x 2 ) L 4 ( ) + w + he x 2 L ( ) (v + h e x 2 ) H ( ) + v 2 + h 2 e x 2 W,4 ( ) (w + he x 2 ) L 4 ( ) + v 2 + h 2 e x 2 L ( ) (w + he x 2 ) H ( ) }. Thus, Sobolev s inequality gives (I ) { div(φ(v + h e x 2 )) } H 3 ( ( ) ) C w H 2 ( ) + h H 2 (R), (6.2)

6 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN where C depends only on v, v 2, h, h 2, and φ. Combining the estimates (6.)-(6.2), we get w C,α (Ω ) g C,α (Ω ) + C ( ) h C ([,T ];C 2,α (I)) + h C([,T ];H 2 (I)) ( ) + CT w C([,T ];H 2 ( )) + h C([,T ];H 2 (R)) g C,α (Ω ) + C h C ([,T ];C 2,α (I)) + C( + T ) h C([,T ];H 2 (R)) + CT w C([,T ];H 2 ( )). Next, by taking C 2,α (Ω ) norm on both sides of (6.2), we have (6.3) w C 2,α (Ω ) g C 2,α (Ω ) + (I ) ( {h x x (x, t) h x x (x, )}e x 2 ) C 2,α (Ω ) (6.4) + t (I ) { div(φ(v + h e x 2 )) } C 2,α (Ω ) dτ. We use Proposition 3.2 to deduce that (I ) { div(φ(v + h e x 2 )) } C 2,α (Ω ) ( C div(φ(v + h e x 2 )) C,α (Ω + div(φ(v + h e x 2 )) ) L. 2 (Ω ) (6.5) For the estimate of the last term in the right hand side of (6.5), we use the proof of Lemma 6. to obtain div(φ(v + h e x 2 )) L 2 (Ω ) C w + he x 2 H 2 ( ). (6.6) In view of (6.4) and the convexity of Ω, div(φ(v + h e x 2 )) C,α (Ω ) φ (v + h e x 2 ) φ (v 2 + h 2 e x 2 ) (v C + h (Ω ) e x 2 ) C,α (Ω ) + φ (v 2 + h 2 e x 2 ) C (w + he x 2 ) (Ω ) C,α (Ω ). Since φ C 3 (R, R 2 ) satisfies the conditions (.5)-(.6), φ (v 2 + h 2 e x 2 ) C (Ω ) C v 2 + h 2 e x 2 C (Ω ) C and φ (v + h e x 2 ) φ (v 2 + h 2 e x 2 ) C (Ω ) φ (v + h e x 2 ) φ (v 2 + h 2 e x 2 ) (v C(Ω + h ) e x 2 ) C(Ω ) + φ (v 2 + h 2 e x 2 ) C(Ω (w + he x 2 ) ) C(Ω ) C w + he x 2 C,α (Ω ),

2D GBBM EQUATION 7 where C depends only on v, v 2, h, h 2, and φ. Thus, we have div(φ(v + h e x 2 )) C,α (Ω ) C w + he x 2 C,α (Ω ) C g C,α (Ω ) + C h C ([,T ];C 2,α (I)) (6.7) + C( + T ) h C([,T ];H 2 (R)) + CT w C([,T ];H 2 ( )), where we used (6.3) in the last inequality. The estimates (6.), (6.4)-(6.7) yield w C ([,T ];C 2,α (Ω )) C( + T ) g C 2,α (Ω ) + C( + T ) h C ([,T ];C 2,α (I)) It follows from (6.7) that + C( + T ) 2 h C([,T ];H 2 (R)) + CT ( + T ) w C([,T ];H 2 ( )). w C ([,T ];C 2,α (Ω )) e CT { g C 2,α (Ω ) + g H 2 ( ) + h C ([,T ];C 2,α (I)) + h C([,T ];H 2 (R)) }. (6.8) Finally, let {(g k, h k )} k= be a sequence of X c that converges to (g, h ) in X c. Suppose that u k = L c (g k, h k ), k N {} be the corresponding classical solutions of (.) with respect to the initial data g k and the boundary data h k. Set, for k N {}, { v k = u k h k e x 2, g k = g k h k e x 2. We define, for k N, w k = v k v, g k, = g k g, h k, = h k h. Since {(g k, h k )} k= converges to (g, h ) in X c, we have (a) d ( g k, g ), (b) g k g H 2 ( ), (c) d 2 (h k, h ), (d) h k h C([,T ];H 2 (R)) ; or equivalently, (a ) g k, C 2,α (Ω i ) for all i N, (b ) g k, H 2 ( ), (c ) h k, C ([,T ];C 2,α (I i )) for all i N, (d ) h k, C([,T ];H 2 (R)). For any fixed i N, applying the estimate (6.8), we get { w k C ([,T ];C 2,α (Ω i )) e CT g k, C 2,α (Ω i+ ) + g k, H 2 ( ) } + h k, C ([,T ];C 2,α (I i+ )) + h k, C([,T ];H 2 (R)),

8 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN and hence L c (g k, h k ) L c (g, h ) C ([,T ];C 2,α (Ω i )) which implies that w k C ([,T ];C 2,α (Ω i )) + h k, C ([,T ];C 2,α (Ω i )) { e CT g k, C 2,α (Ω i+ ) + g k, H 2 ( ) } + h k, C ([,T ];C 2,α (I i+ )) + h k, C([,T ];H 2 (R)), d 3 (L c (g k, h k ), L c (g, h )) which is equivalent to that the mapping L c : X c C ([, T ]; C 2,α loc (R2 +)) is sequentially continuous. 7. Results for the GBBM-Burgers equation The purpose of this section is to generalize the above results to the 2D GBBM- Burgers equation, that is, equation (.) with ν =. The results established in previous sections also hold for the GBBM-Burgers equation. The proofs of Theorems. and.2 for the case when ν = are essentially the same as those for the case when ν =. In fact, as in the case when ν =, we rewrite equation (.) as (I )v t + div ( φ(v + he x 2 ) ) = h x x te x 2 in (, T ), (7.a) v = g on {t = }, (7.b) v = on (, T ), (7.c) where g is again given by (2.2). (7.a) can be modified as v + v t = (I ) {v + h x x te x 2 div(φ(v + he x 2 ))} (7.2) which suggest that t v(x, t) = e t g(x) + e (t s) (I ) {v + h x x te x 2 div(φ(v + he x 2 ))}ds = e t g(x) + (I ) [ ] h x x (x, t)e x 2 e t h x x (x, )e x 2 (7.3) + t e (t s) (I ) {v h x x e x 2 div(φ(v + he x 2 ))}ds. Based on the fact that e (t s) is bounded by for s [, t], exactly the same procedure of proving the existence of a unique solution (using the contraction mapping principle) for the case ν = can be applied to yield the results corresponding to Theorems. and.2.

2D GBBM EQUATION 9 Acknowledgments. AC was supported by the National Science Council (NSC) of Taiwan under grant -25-M-8-9-MY3, MH was supported by NSC under grant -25-M-8-5, YL was supported by NSC under grant 2-25-M-8-, JW was supported by NSF under grant DMS 2953, and JY was supported by NSC under grant -25-M-26-2. References [] J. Avrin and J. A. Goldstein, Global existence for the Benjamin-Bona-Mahony equation in arbitrary dimensions, Nonlinear Anal. 9 (985), pp. 86-865. [2] J. L. Bona and P. J. Bryant, A mathematical model for long waves generated by wavemakers in non-linear dispersive systems, Proc. Camb. Phil. Soc. 73 (973), pp. 39-45. [3] T. B. Benjamin, J. L. Bona, and J. J. Mahony, Model equations for long waves in nonlinear dispersive system, Phil. Trans. R. Soc. London., Ser A 272 (972), pp. 47-78. [4] J. L. Bona and H. Chen, Well-posedness for regularized nonlinear dispersive wave equations, Discrete Contin. Dyn. Syst. 23 (29), pp. 253275. [5] J. L. Bona, H. Chen, S. M. Sun, and B.-Y. Zhang, Approximating initial-value problems with two-point boundary-value problems: BBM-equation, Bull. Iranian Math. Soc. 36 (2), pp. -25. [6] J. L. Bona, H. Chen, S. M. Sun, and B.-Y. Zhang, Comparison of quarter-plane and two-point boundary value problem: the BBM-equation, Discrete Contin. Dyn. Syst. 3 (25), pp. 92-94. [7] J. L. Bona and V. A. Dougalis, An initial- and boundary-value problem for a model equation for propagation of long waves, J. Math. Anal. Appl. 75 (98), pp. 53-522. [8] J. L. Bona and L. Luo, More results on the decay of solutions to nonlinear, dispersive wave equations, Discrete Contin. Dyn. Syst. (995), pp. 5-93. [9] J. L. Bona and L. Luo, Initial-boundary value problems for model equations for the propagation of long waves. Evolution equations (Baton Rouge, LA, 992), 63, 65V94, Lecture Notes in Pure and Appl. Math., 68, Dekker, New York, 995. [] J. L. Bona, W. G. Pritchard, and L. R. Scott, An evaluation of a model equation for water waves, Phil. Trans. R. Soc. London., Ser A 32 (98), PP. 457-5. [] J. L. Bona and N. Tzvetkov, Sharp well-posedness results for the BBM equation, Discrete Contin. Dyn. Syst. 23 (29), pp. 24-252. [2] J. L. Bona and J. Wu, Temporal growth and eventual periodicity for dispersive wave equations in a quarter plane, Discrete Contin. Dyn. Syst. 23 (29), pp. 4-68. [3] S. E. Buckley and M. C. Leverett, Mechanism of fluid displacements in sands, Trans. AIME 46 (942), pp. 7-6. [4] D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer- Verlag Berlin Heidelberg New York, 2. [5] J. A. Goldstein and B. J. Wichnoski, On the Benjamin-Bona-Mahony equation in higher dimensions, Nonlinear Anal. 4 (98), pp. 665-675. [6] J.-M. Hong, J. Wu, and J.-M. Yuan, The generalized Buckley-Leverett and the regularized Buckley-Leverett equations, J. Math. Phys. 53 (22), 537, 24 pp. [7] J.-M. Hong, J. Wu, and J.-M. Yuan, A new solution representation for the BBM equation in a quarter plane and the eventual periodicity, Nonlinearity 22 (29), pp. 927-944. [8] G. Karch, Asymptotic behaviour of solutions to some pseudoparabolic equations, Math. Methods Appl. Sci. 2 (997), pp. 27-289. [9] Y. Martel, F. Merle, and T. Mizumachi, Description of the inelastic collision of two solitary waves for the BBM equation, Arch. Rational Mech. Anal. 96 (2), pp. 57-574. Ying-Chieh Lin, C. H. Arthur Cheng, and John M. Hong Department of Mathematics National Central University

2 Y.-C. LIN, C. H. A. CHENG, J. M. HONG, J. WU, AND J.-M. YUAN Chung-Li, Taiwan 32, ROC E-mail address: linyj@math.ncu.edu.tw; cchsiao@math.ncu.edu.tw; jhong@math.ncu.edu.tw Jiahong Wu Department of Mathematics Oklahoma State University 4 Mathematical Sciences Sillwater, OK 7478, USA E-mail address: jiahong@math.okstate.edu Juan-Ming Yuan Department of Financial and Computational Mathematics Providence University Taichung, Taiwan 433, ROC E-mail address: jmyuan@pu.edu.tw