arxiv: v1 [math.na] 17 Jan 2019

Transcription

1 arxiv: v1 [mat.na 17 Jan 019 Teoretical an numerical stuies for energy estimates of te sallow water equations wit a transmission bounary conition M. Masum Murse 1,, Kouta Futai 1, Masato Kimura 3 an Hirofumi Notsu 3,4 1. Division of Matematical an Pysical Sciences, Kanazawa University, Kanazawa , Japan. Department of Matematics, University of Rajsai, Rajsai 605, Banglaes 3. Faculty of Matematics an Pysics, Kanazawa University, Kanazawa , Japan 4. Japan Science an Tecnology Agency, PRESTO, Kawaguci , Japan. s: mmmurse8@gmail.com, futai.k.175@stu.kanazawa-u.ac.jp, mkimura@se.kanazawa-u.ac.jp an notsu@se.kanazawa-u.ac.jp Abstract Energy estimates of te sallow water equations (SWEs) wit a transmission bounary conition are stuie teoretically an numerically. In te teoretical part, using a suitable energy, we begin wit eriving an equality wic implies an energy estimate of te SWEs wit te Diriclet an te slip bounary conitions. For te SWEs wit a transmission bounary conition, an inequality for te energy estimate is prove uner some assumptions to be satisfie in practical computation. Hence, it is recognize tat te transmission bounary conition is reasonable in te sense tat te inequality ols true. In te numerical part, base on te teoretical results, te energy estimate of te SWEs wit a transmission bounary conition is confirme numerically by a finite ifference meto (FDM). Te coice of a positive constant c 0 use in te transmission bounary conition is investigate aitionally. Furtermore, we present numerical results by a Lagrange Galerkin sceme, wic are similar to tose by te FDM. From te numerical results, it is foun tat te transmission bounary conition works well numerically. 1 Introuction Te sallow water equations (SWEs) are often use for te simulation of tsunami/ storm surge in te bay. In suc simulation tere are some bounaries in te open sea, see Figure 1. In a real situation, if wave propagates towars suc 1

2 bounaries in te open sea, ten tere soul not ave any reflection on tese bounaries. Terefore, in te simulation a special type of bounary conition soul be impose on tese bounaries. In tis stuy, following [6, we employ a transmission bounary conition on te bounaries in te open sea wic is capable to remove tese kin of artificial reflection. Many stuies ave been conucte on te preiction of te surges ue to tropical storms for te Bay of Bengal region covering te coast of Banglaes an te east coast of Inia, see, e.g., [ 4,8 15. Almost all of te stuies mentione ere ave employe a raiation type bounary conition for te bounaries in te open sea, wic is very similar to te transmission bounary conition use in [6. It is to be note ere tat tese papers mainly present numerical results witout matematical iscussion on te stability of te moel wit tis kin of bounary conition. Figure 1: Te Bay of Bengal an te coastal region of Banglaes Te SWEs can be consiere as a couple system of a pure convection equation for te function φ of total wave eigt an a simplifie Navier Stokes equation for te velocity u = (u 1, u ) T by averaging function values in x 3 - irection. It is known tat a bounary ata for φ is necessary on te so-calle inflow bounary, were u n < 0 is satisfie for te outwar unit normal vector n. We can easily know weter te Diriclet ata for φ is require or not on te Diriclet an te slip bounaries for u, since te sign of u n is known a priori. On te transmission bounary Γ T, owever, te bounary conition for u an φ is mysterious an problematic from bot computational an matematical view points. Te transmission bounary conition of te form η(x, t) u(x, t) = c(x) n(x) (1) φ(x, t) is often use on Γ T, were c(x) is a given positive value function an η(x, t) = φ(x, t) ζ(x) is te elevation from te reference eigt for a given ept function ζ. Let φ k an uk be te approximations of φk := φ(, t k ) an u k := u(, t k ), respectively, were t k := k t (k Z) for a time increment t. In our computation we get φ k+1 by using u k, an ten uk+1 by using te conition (1) as

3 te Diriclet bounary conition for u. But at te same time we soul consier (1) as te bounary ata for φ k+1 if te position is on inflow bounary, i.e., u k+1 n < 0. In fact, if u k+1 n < 0, we nee to give te value of φ k+1 wic is unknown. We ave an experience of computation of te SWEs wit te transmission bounary conition by a Lagrange Galerkin (LG) meto, were te LG meto is base on te time iscretization of te material erivative, φ k+1 (x) φ k (x u k (x) t). t Te position x u k (x) t is te so-calle upwin point of x wit respect to u k. In te computation a nearest bounary value of φ k is use if te upwin point places outsie te omain, an te LG meto works witout bounary ata for φ k+1 even if u k+1 n < 0. In te LG meto te problem on te transmission bounary seems to be solve numerically, but it is still problematic matematically. In tis paper, in orer to unerstan te transmission bounary conition matematically, we stuy te stability of te SWEs in terms of a suitable energy, an confirm te stability numerically by a finite ifference meto (FDM). Since te LG meto solves te problem implicitly by using te upwin point an is not suitable to unerstan te problem wit te transmission bounary conition matematically, we employ an FDM to avoi te special numerical treatment in te LG meto. It is to be note ere tat we can sow a (successful) energy estimate of te SWEs, wen only te Diriclet an te slip bounary conitions are employe, cf. Corollary 3.3-(ii), were suc iscussions ave been one uner te perioic bounary conition, e.g., [1,7. As far as we know, owever, tere is no matematical results on te energy estimate of te SWEs wit te transmission bounary conition. Te stability is consiere teoretically wit respect to te energy as follows. Introucing a suitable energy, we begin wit eriving an equality tat te time-erivative of te energy consists of four terms, were tree terms are line integrals over te bounary an te oter term is an integral over te wole omain wic is always non-positive, cf. Teorem 3.1. Since te tree line integrals vanis over te Diriclet an te slip bounaries, as a result, we ave te tree line integrals over te transmission bounary an te integral over te omain, cf. Corollary 3.3-(i). An energy estimate is obviously obtaine if tere is no transmission bounary, cf. Corollary 3.3-(ii). In aition, we obtain tat a sum of two line integrals over te transmission bounary is non-positive uner some conitions to be satisfie in real computations, cf. Teorem 3.4. Altoug, at present, te matematical results o not erive te stability estimate of te SWEs wit te transmission bounary conition irectly, we ave goo information an can stuy te stability numerically by using te teoretical results. In te latter alf of te paper, Sections 4 an 5, te energy estimates erive in Section 3 are stuie numerically. Tis paper organize as follows. Te problem is state matematically in Section. Te energy estimate is teoretically stuie in Section 3. Numerical 3

4 results by an FDM are sown to see tat te transmission bounary conition works well an tat te solution is stable in terms of te energy in Section 4. Here, te effect of a positive constant to be use in te transmission bounary conition is investigate. In Section 5, numerical results by an LG meto are sown to support tose by te FDM. Finally, conclusions are given in Section 6. Statement of te problem In tis section, we state te matematical problem to be consiere in tis paper. Let R be a boune omain an T a positive constant. We consier te problem; fin (φ, u) : [0, T R R suc tat φ + (φu) = 0 in (0, T ), t [ u ρφ t + (u )u () µ (φd(u)) + ρgφ η = 0 in (0, T ), φ = η + ζ in (0, T ), wit bounary conitions an initial conitions u = 0 on Γ D (0, T ), (3) (D(u)n) n = 0, u n = 0 on Γ S (0, T ), (4) u = c η φ n on Γ T (0, T ), (5) u = u 0, η = η 0 in, at t = 0, (6) were φ is te total eigt of wave, u = (u 1, u ) T is te velocity, η : [0, T R is te water level from te reference eigt, ζ(x) > 0 (x ) is te ept of water from te reference eigt, see Figure, D(u) := ( u + ( u) T ) / is te strain-rate tensor, n is te unit outwar normal vector, Γ := is te bounary of, we assume tat Γ consists of non-overlappe tree parts, Γ D, Γ S an Γ T, i.e., Γ = Γ D Γ S Γ T, Γ D Γ S =, Γ S Γ T =, Γ T Γ D =, te subscripts D, S, an T mean Diriclet, slip, an transmission bounaries, respectively, ρ > 0 is te ensity of water, µ > 0 is te viscosity, g > 0 is te acceleration ue to gravity, an c(x) := c 0 gζ(x) wit a positive constant c0. In te rest of paper, we assume ζ C 1 (). It is important to note ere tat te equations in () are erive in [5 by consiering one-layer viscous SWEs. 3 Energy estimate 4

5 x # x $ x " η u " φ u # Figure : Moel omain In tis section, we efine te total energy an stuy te stability of solutions to te problem state in Section in terms of te energy. For a solution of () te total energy E(t) at time t [0, T is efine by E(t) := E 1 (t) + E (t), (7) were E 1 (t) an E (t) are te kinetic an te potential energies efine by ρ ρg η E 1 (t) := φ u x, E (t) := x. Let symbols I i (t; Γ ), i = 1,..., 3, an I 4 (t; ), t [0, T, be integrations efine by I 1 (t; Γ ) := ρ φ u u n s, I (t; Γ ) := ρg φηu n s, Γ Γ I 3 (t; Γ ) := µ φ [ D(u)n u s, I 4 (t; ) := µ φ D(u) x. Tese are use in te rest of tis paper. Γ Teorem 3.1. Suppose tat a pair of smoot functions (φ, u) : [0, T R R satisfies (). Ten, we ave t E(t) = 3 I i (t; Γ ) + I 4 (t; ). (8) i=1 We prove Teorem 3.1 after preparing a lemma. Lemma 3.. For smoot functions φ : [0, T R an u : [0, T R, we ave te following. (i) ( ) ( ) φ u t (φu) + [(φu) u = t + (φu) u + φ + (u )u, t (ii) [(φu) u) ux = ( 1 φ u u ns + 1 [ (φu) u x. Γ 5

6 Proof. We prove (i). From te ientity, we ave [(φu) u = [ (φu)u + φ[(u )u, (9) φ (φu) + [(φu) u = t t u + φ u + [ (φu)u + φ[(u )u ( t ) ( ) φ u = t + (φu) u + φ + (u )u, t wic completes te proof of (i). We prove (ii). Denoting te left an sie of (ii) by J an using te integration by parts formula, we ave J = ([(φu) un) us [(φu) u : u x Γ = φ u u ns φ(u )u u x, (10) Γ were A : B = i,j=1 A ijb ij. We, terefore, ave J = [ (φu) u + (φ(u )u) u x (from (9)) = [ (φu) u x + φ u u n s J (from (10)), wic implies te esire result of (ii). Proof of Teorem 3.1. Differentiating (7) wit respect to t, we get Γ t E(t) = t E 1(t) + t E (t). (11) We compute t E 1(t) an t E (t) separately. Firstly, t E 1(t) is compute as follows. From Lemma 3.-(i) an te first equation of (), we ave wic implies [ u φ t + (u )u = (φu) + [(φu) u, t [ ρ t (φu) + [(φu) u µ [φd(u) + ρgφ η = 0. (1) Multiplying (1) by u an integrating wit respect to x over, we get [ [ [ ρ t (φu) u x + ρ [(φu) u u x µ (φd(u)) u x 6

7 + ρg φ η u x = 0. (13) From te equation (13) above an te next two ientities: [ [ ρ t (φu) u x + ρ [(φu) u u x ( φ = ρ t u + φ u ) t u x + ρ φ u u n s + ρ = ρ = t 1 ( [ ρ φ t u + φu u t φ u x + ρ Γ [ (φu) u x (from Lemma 3.-(ii)) ) x + ρ φ u u n s (from te first eq. of ()) Γ Γ φ u u n s = t E 1(t) I 1 (t; Γ ), [ µ (φd(u)) u x = µ φ [ D(u)n u s + µ φ D(u) x Γ = I 3 (t; Γ ) I 4 (t; ), we obtain t E 1(t) = I 1 (t; Γ ) + I 3 (t; Γ ) + I 4 (t; ) ρg η (φu) x. (14) Seconly, t E (t) is compute as follows: t E (t) = [ ρg η x t = ρg η η t x = ρg η φ x (from te tir eq. of ()) t = ρg η [ (φu) x (from te first eq. of ()) = ρg (ηφu)x + ρg η (φu) x = I (t; Γ ) + ρg η (φu) x. (15) Te result (8) follows by aing (14) an (15) an recalling (11). Corollary 3.3. (i) Suppose tat a pair of smoot functions (φ, u) : [0, T R R satisfies () wit (3)-(5). Ten, we ave t E(t) = 3 I i (t; Γ T ) + I 4 (t; ). (16) i=1 7

8 (ii) Furtermore, if Γ = Γ D Γ S an φ(x, t) > 0 ((x, t) [0, T ), we ave t E(t) = I 4(t; ) 0. (17) Proof. On Γ S, from te first equation of (4), tere exists a scalar function w : [0, T R suc tat D(u)n = w(x, t)n, wic implies [ D(u)n u = (wn) u = w(u n) = 0. Hence, te result (16) is establise from Teorem 3.1 wit (3) an (4). Wen Γ = Γ D Γ S, i.e., Γ T =, te ientity (16) erives (17). Teorem 3.4. Suppose tat a pair of smoot functions (φ, u) : [0, T R R satisfies () wit (3)-(5), an an inequality an tat tere exists α (0, 1) suc tat Ten, we ave te following estimates: φ(x, t) > 0, (x, t) Γ T [0, T, (18) η(x, t) αζ(x), x Γ T, t [0, T, (19) 0 < c 0 (1 α). (0) α I 1 (t; Γ T ) + I (t; Γ T ) 0, (1) in particular, t E(t) I 3(t; Γ T ). () Proof. From (16), we prove (1) wic implies (), since I 4 (t; ) is always non-positive. We ave i=1 [ I i (t; Γ T ) = ρ φ(u n) Γ T = ρ φc η [gη + 1 Γ T φ c 0gζ η [ = ρg cη Γ T gη + 1 u s 1 + c 0 φ s ζη (ζ + η) s. Let f(r) := r/(1+r). From f (r) = (1 r)/(1+r) 3, it ols tat f(r 1 ) f(r ), 1 r 1 r 1. Since 1 α η/ζ 1, we obtain f( α) f(η/ζ), i.e., α (1 α) ηζ (ζ + η), 8

9 wic implies tat i=1 { I i (t; Γ T ) ρg cη c } 1 0α Γ T (1 α) s 0 from te conition on c 0, i.e., 0 < c 0 /α (1 α). Remark 3.5. From Teorem 3.4, we can say tat te transmission bounary conition (5) is reasonable in te sense of (1) uner te conitions (18)-(0) to be satisfie in practical computation. Altoug te sign of te(t) is as yet unknown ue to I 3 (t; Γ T ), Teorem 3.4 as meaningful for te energy estimate of te SWEs wit te transmission bounary conition. In fact, we observe numerically tat I (t; Γ ) is ominant an 3 i=1 I i(t; Γ ) is negative, wile I 1 (t; Γ ) an I 3 (t; Γ ) may be positive, cf. Subsection 4.. Remark 3.6. Te conition (0) is not strict in te practical computation, were α an c 0 are cosen typically as, e.g., α = 0.01 an c 0 = 0.9 [6. Tese satisfy (0), since /α (1 α) Numerical results by a finite ifference sceme In tis section, we present numerical results by a finite ifference sceme for problem () (6) wit = (0, L) for a positive constant L, T = 100, ζ = a > 0, µ = 10 3, g = , ρ = 10 1, η 0 = c 1 exp( 100 x p ) (c 1 > 0, p ). Tese values are in km (lengt), kg (mass) an s (time). We set Γ S = for simplicity. We consier five cases of Γ T : (i) Γ T =, (ii) Γ T = Γ top, (iii) Γ T = Γ top Γ rigt {(L, L)}, (iv) Γ T = Γ top Γ rigt Γ left {(L, L)} {(0, L)}, (v) Γ T = Γ, for Γ top := {(x 1, L); 0 < x 1 < L}, Γ rigt := {(L, x ); 0 < x < L}, Γ left := {(0, x ); 0 < x < L}, an set Γ D = Γ \ Γ T. For te above cases (ii)-(v), c 0 = 0.9 is taken following [ A finite ifference sceme Let N N an t > 0 be given, an let := L/N an N T := T/ t, x i,j := (i, j) T R (i, j Z), := {x i,j ; i, j Z}, := {x i,j ; i, j Z}, Γ D := {x i,j Γ D ; i, j Z}, Γ T := {x i,j Γ T ; i, j Z}. Let u 0 : R an φ 0 : R be given approximate functions of u 0 an φ 0, respectively. Our finite ifference sceme is to fin {(φ k, uk )(x i,j); x i,j, k = 1,..., N T } 9

10 suc tat, for k = 0,..., N T 1, φ k+1 φ k t ρφ k + ( u k ) φ k + u k up ) ( u k+1 u k + (u k )u k t φk = 0 on, µ (φ k D (u k ) ) + ρgφ k η k = 0 on, φ k = η k + ζ on, u k+1 = u k+1 D on Γ D, u k+1 = c φk+1 φ k+1 ζ n on Γ T, (3) were = ( 1, ) T an up = ( up 1, up )T represent te (stanar) central an upwin (wit respect to u k ) ifference operators, respectively, an D (v ) := [( v )+( v ) T / : R sym for v : R. If te require points for te operators an up are not in, one sie ifference is use. 4. Numerical results for five cases of bounary settings Numerical simulations are carrie out by sceme (3) for L = 10, a = 1, u 0 = 0, c 1 = 0.01, p = (5, 5) T, N = 1, 000 an t = 0.05 (N T =, 000). Figure 3 sows color contours of η k for k = 0, 500, 1, 000, 1, 500 an, 000, wic correspon to times t = 0, 5, 50, 75 an 100, respectively, were (i)-(v) represent simulate results for te cases (i)-(v) state at te beginning of tis section. It can be clearly foun tat te artificial reflection is almost remove on te transmission bounaries for te cases (ii)-(v). 4.3 Numerical stuy of energy estimate In tis subsection, we stuy te stability of solutions to te problem ()-(6) numerically by sceme (3) in terms of te energy E(t) efine in (7). Te values of E(t k ) an I i (t k ; Γ ), i = 1,, 3, I 4 (t k ; ) are approximately compute by using solution {(u k, φk )}N T k=1 wit {ηk }N T k=1 of sceme (3) as E(t k ) E k := ρ φ k (x i,j ) u k (x i,j ) + ρg η (x k i,j ), x i,j x i,j I 1 (t k ; Γ ) I1 k := ρ (Π φ k ) Π u k (Π u k ) n s, Γ I (t k ; Γ ) I k := ρg (Π φ k )(Π η)(π k u k ) n s, N I 3 (t k ; Γ ) I3 k := µ ( Γ m=1 l T m ) (Π φ k )(D(Π u k )n) (Π u k ) s, l B m l L m l R m 10

11 (i) t = 0 t = 5 t = 50 t = 75 t = 100 (ii) t = 0 t = 5 t = 50 t = 75 t = 100 (iii) t = 0 t = 5 t = 50 t = 75 t = 100 (iv) t = 0 t = 5 t = 50 t = 75 t = 100 (v) t = 0 t = 5 t = 50 t = 75 t = 100 Figure 3: Color contours of η k by finite ifference sceme (3) for te five cases (i)-(v) iscusse in Subsection

12 E k 4 i=1 Ik i t E(tk ) I k i, i = 1,..., 4 (i) (ii) (iii) (iv) (v) Figure 4: Graps of E k (left), 4 i=1 Ik i t E(t) (center) an Ik i, i = 1,..., 4, (rigt) versus t = t k ( 0, k Z) for te five cases (i)-(v). 1

13 N I 4 (t k ; ) I4 k := µ (Π φ k )(x i 1/,j 1/ ) D (u k )(x i 1/,j 1/ ), i,j=1 were Π f C(; R) is te bilinear interpolation of f : R for f = φ, η, Π u = (Π u 1, Π u ) T C(; R ), te bounary is represente as Γ = N m=1 (lt m l B m l L m l R m) for line segments l T m, l B m, l L m, l R m efine by l T m := x m 1,N x m,n, l B m := x m 1,0x m,0, l L m := x 0,m 1x 0,m, l R m := x N,m 1x N,m, an te omain is represente as = N i,j=1 ω i 1/,j 1/ for ω i 1/,j 1/ := ((i 1), i) ((j 1), j) wit te area. Numerical simulations for te problem ()-(6) wit L = 1, a = 0.1, u 0 = 0, c 1 = 0.001, p = (0.5, 0.5) T are carrie out by sceme (3) wit t = for = L/N = 1/N, N = 400, 500, 800 an 1, 000. Te results are presente in Figure 4, were (i)-(v) in te figure represent te cases (i)-(v) at te beginning of tis section. Te graps of E k an 4 i=1 Ik i versus t = tk (k N) are presente in te left an center figures, respectively. Tere are four lines in eac figure, but te lines are almost overlappe in te cases of (ii)-(v). In te case of (i) te graps are qualitatively similar. Te rigt figures sow te graps of Ii k, i = 1,..., 4 versus t = tk (k N). Te maximum an minimum values of Ii k, i = 1,..., 4, on te time interval (0, T ) are presente in Table 1. From te numerical results presente in Figure 4, it can be foun tat te total energy is mainly ecreasing wit respect to time. In te case of (i), i.e., Γ = Γ D, we can see tat at te early perio te graps are increasing in te left figure, wile te values are small. From te center figures it can be clearly seen tat te sum 4 i=1 Ik i corresponing to te erivative of te total energy is always non-positive, wic confirms te stability of solutions to te moel numerically. From rigt figures an Table 1, it can be observe tat te value of I is ominating negatively over I 1 an I 3 so tat te sum 4 i=1 I i becomes non-positive always. 4.4 Coice of c 0 In tis subsection, we stuy te value of c 0 by solving te problem wit L = 1, a = 0.1, p = (0.5, 0.5) T an c 0 = 0.1, 0.,..., 1. an 1.5 by sceme (3) wit N = 400, t = (N T = 0, 000), an by computing { S(c k 0 ) := η(x k i,j ) x i,j { N T S (c 0 ) := t k=0 } 1 S(t k ; c 0 ) := { S k (c 0 ) } 1 S(c0 ) := { T 0 } 1 η k (x) x = η k L (), } 1 S(t k ; c 0 ) t = η L (0,T ;L ()). We set six ifferent cases for te initial values of η 0 an u 0 as follows. Case I: η 0 = c 1 exp ( 100 x (0.5, 0.5) ), u 0 = 0, 13

14 Table 1: Maximum an minimum values of Ii k, i = 1,..., 4, wit respect to te number of transmission bounaries Γ T Γ D I 1 I I 3 I Max Min Max Min Max Min Max Min Max Min Case II: η 0 = c 1 exp ( 00 x (0.5, 0.5) ), u 0 = 0, Case III: η 0 = c 1 exp ( 100 x (0, 0.5) ), u 0 = 0, Case IV: η 0 = c 1 exp ( 100 x ), u 0 = 0, Case V: η 0 = c 1 exp ( 100 x (0.5, 0.5) ), u 0 = 10 4 (1, 1) T, Case VI: η 0 = c 1 exp ( 100 x (0.5, 0.5) ), u 0 = 10 4 (1, 1) T, were c 1 = Since te artificial reflection soul be remove after te time te wave touces te transmission bounary, we fin a value of c 0 wic provies te minimum of S (c 0 ). Te results are presente in Table, from were it can be conclue tat for te case of zero initial velocity te suitable value of c 0 lies in [0.7, 1.0 an for te case of nonzero initial velocity we cannot say anyting yet. 5 Numerical results by an LG sceme In tis section, we present an LG sceme for te same problem escribe in Subsection 4.. Let T = {K} be a triangulation of, an M te so-calle P1 (piecewise linear) finite element space. We set Ψ := M for te water level η, an V (ψ ) := v M; v (P ) = c(p ) ψ (P ) ζ(p ) ψ (P ) v (Q) = 0, n(p ), P : noe on Γ T, Q : noe on Γ D 14

15 Table : c 0 an S (c 0 ) c 0 S (c 0 ) Case I Case II Case III Case IV Case V Case VI for te velocity u. Te LG sceme is to fin {(φ k, uk )}N T k=1 Ψ V suc tat, for k = 1,..., N T, ρ φ k φ k φ k 1 X k 1 1 t u k ũk 1 t +ρg γk 1 X k 1 1 ψ x = 0, ψ Ψ, v x + µ φ k D(u k ) : D(v ) x φ k η k v x = 0, v V, φ k = η k + Π FEM ζ, were X1 k (x) := x uk (x) t, γk : R is efine by ( X γ(x) k k := et 1 (x) ), x (4) te symbol represents te composition of functions, i.e., [v X1 k (x) := v (X1 k (x)), ΠFEM : C() M is te Lagrange interpolation operator, an { ψ (x), x, ψ (x) = ψ (P x ), x R \, were P x Γ is a nearest noal point from x. In eac step, firstly, φ k Ψ is obtaine from te first equation of sceme (4). Seconly, u k V is obtaine by using φ k from te secon equation. In te first equation of (4), te iea of mass conservative Lagrange Galerkin sceme [16 is employe. 15

16 A numerical simulation is carrie out by LG sceme (4) wit te same setting escribe in Subsection 4. except t, were it is set as t = Te results are sown in Figure 5. Te figures are similar to te ones in Figure 3 obtaine by finite ifference sceme (3) an support te results in Figure 3. (i) t = 0 t = 5 t = 50 t = 75 t = 100 (ii) t = 0 t = 5 t = 50 t = 75 t = 100 (iii) t = 0 t = 5 t = 50 t = 75 t = 100 (iv) t = 0 t = 5 t = 50 t = 75 t = 100 (v) t = 0 t = 5 t = 50 t = 75 t = 100 Figure 5: Color contours of η k by Lagrange Galerkin sceme (4) for te five cases (i)-(v) iscusse in Section 5. 16

17 6 Conclusions Energy estimates of te SWEs wit a transmission bounary conition ave been stuie matematically an numerically. For a suitable energy, we ave obtaine an equality tat te time-erivative of te energy is equal to a sum of tree line integrals an a omain integral in Teorem 3.1. Te teorem implies a (successful) energy estimate of te SWEs wit te Diriclet an te slip bounary conitions, cf. Corollary 3.3-(ii). After tat, an inequality for te energy estimate of te SWEs wit te transmission bounary conition as been prove in Teorem 3.4. In te proof, it as been sown tat a sum of two line integrals over te transmission bounary is non-positive uner some conitions to be satisfie in practical computation. Base on te teoretical results, te energy estimate of SWEs wit te transmission bounary conition as been confirme numerically by an FDM. It is foun tat te transmission bounary conition works well numerically an tat te transmission bounary conition reuces te energy rastically via te term I k. Te coice of a positive constant c 0 use in te transmission bounary conition as been investigate aitionally, an it as been observe tat te suitable value lies in [0.7, 1.0 in te case of zero initial velocity. Furtermore, we ave presente numerical results by an LG sceme, wic are similar to tose by te FDM. Completeness of te (teoretical) energy estimate of te SWEs wit te transmission bounary conition is a future work. Acknowlegements M.M.M. is supporte by MEXT scolarsip. Tis work is partially supporte by JSPS KAKENHI Grant Numbers JP16H0155, JP17H0857 an JP18H01135, JSPS A3 Foresigt Program, an JST PRESTO Grant Number JPMJPR16EA. References [1 D. Bresc an B. Desjarins, On te construction of approximate solutions for te D viscous sallow water moel an for compressible Navier Stokes moels, Journal e Matématiques Pures et Appliquées, 86 (006), [ P. K. Das, Preiction moel for storm surges in te Bay of Bengal, Nature, 39 (197), [3 S. K. Debsarma, Simulations of storm surges in te Bay of Bengal, Marine Geoesy, 3 (009), [4 B. Jons an A. Ali, Te numerical moeling of storm surges in te Bay of Bengal, Quarterly Journal of te Royal Meteorological Society, 106 (1980),

18 [5 H. Kanayama an H. Dan, A finite element sceme for two-layer viscous sallow water equations, Japan Journal of Inustrial an Applie Matematics, 3, No. (006), [6 H. Kanayama an H. Dan, Tsunami Propagation from te Open Sea to te Coast, Tsunami, Capter 4, IntecOpen, 016, [7 C. Lucas, Cosine Effect on Sallow Water Equations an Matematical Properties, Quarterly of Applie Matematics, American Matematical Society, 67 (009), [8 G. C. Paul an A. M. I. Ismail, Tie surge interaction moel incluing air bubble effects for te coast of Banglaes, Journal of te Franklin Institute, 349 (01), [9 G. C. Paul an A. M. I. Ismail, Contribution of offsore islans in te preiction of water levels ue to tie-surge interaction for te coastal region of Banglaes, Natural Hazars, 65, No.1 (013), [10 G. C. Paul, A. M. I. Ismail an M. F. Karim, Implementation of meto of lines to preict water levels ue to a storm along te coastal region of Banglaes, Journal of Oceanograpy, 70 (014), [11 G. C. Paul, M. M. Murse, M. R. Haque, M. M. Raman an A. Hoque, Development of a cylinrical polar coorinates sallow water storm surge moel for te coast of Banglaes, Journal of Coastal Conservation, 1, No.6 (017), [1 G. C. Paul, S Sentilkumar an R. Pria, Storm surge simulation along te Megna estuarine area: an alternative approac, Acta Oceanologica Sinica, 37, No.1 (018), [13 G. C. Paul, S Sentilkumar an R. Pria, An efficient approac to forecast water levels owing to te interaction of tie an surge associate wit a storm along te coast of Banglaes, Ocean Engineering, 148 (018), [14 G. D. Roy, A. B. Humayun Kabir, M. M. Manal an M. Z. Haque, Polar coorinate sallow water storm surge moel for te coast of Banglaes, Dynamics of Atmosperes an Oceans, 9 (1999), [15 G. D. Roy an A. B. H. M. Kabir, Use of neste numerical sceme in a sallow water moel for te coast of Banglaes, BRAC University Journal, 1, No.1 (004), [16 H. Rui an M. Tabata. A mass-conservative caracteristic finite element sceme for convection-iffusion problems, Journal of Scientific Computing, 43 (010),