arxiv: v2 [math.na] 11 Dec 2016

Transcription

1 Noname manuscript No. will be inserted by te editor Sallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical metods Hendrik Ranoca arxiv: v [mat.na] Dec 06 Received: date / Accepted: date Abstract For te first time, a general two-parameter family of entropy conservative numerical fluxes for te sallow water equations is developed and investigated. Tese are adapted to a varying bottom topograpy in a well-balanced way, i.e. preserving te lake-at-rest steady state. Furtermore, tese fluxes are used to create entropy stable and well-balanced split-form semidiscretisations based on general summation-by-parts SBP operators, including Gauß nodes. Moreover, positivity preservation is ensured using te framework of Zang and Su Maximum-principle-satisfying and positivity-preserving igorder scemes for conservation laws: survey and recent developments, 0. In: Proceedings of te Royal Society of London A: Matematical, Pysical and Engineering Sciences, Te Royal Society, vol 67, pp Terefore, te new two-parameter family of entropy conservative fluxes is enanced by dissipation operators and investigated wit respect to positivity preservation. Additionally, some known entropy stable and positive numerical fluxes are compared. Furtermore, finite volume subcells adapted to nodal SBP bases wit diagonal mass matrix are used. Finally, numerical tests of te proposed scemes are performed and some conclusions are presented. Keywords Skew-symmetric sallow water equations Summation-by-parts Split-form Entropy stability Well-balancedness Positivity preservation Matematics Subject Classification 000 MSC 65M70 MSC 65M60 65M06 65M Hendrik Ranoca TU Braunscweig, Pockelsstraße, 06 Braunscweig, Germany Tel.: , Fax: ranoca@tu-bs.de

2 Hendrik Ranoca Introduction For te first time, a two-parameter family of entropy conservative and wellbalanced numerical fluxes for te sallow water equations is investigated, resulting in genuinely ig-order semidiscretisations tat are bot entropy stable and positivity preserving. Tese semidiscretisations of te sallow water equations in one space dimension are based on summation-by-parts SBP operators, see inter alia te review articles by Svärd and Nordström 0; Fernández et al 0 and references cited terein. Tis setting of SBP operators originates in te finite difference FD setting, but can also be used in polynomial metods as nodal discontinuous Galerkin DG Gassner 0 or flux reconstruction / correction procedure via reconstruction Ranoca et al 06. Entropy stability as long been known as a desirable stability property for conservation laws. Here, te semidiscrete setting of Tadmor 97, 00 will be used. Oter desirable stability properties for te sallow water equations are te preservation of non-negativity of te water eigt and te correct andling of steady states, especially te lake-at-rest initial condition, resulting in wellbalanced metods. References for numerical metods for te sallow water equations can be found in te review article of Xing and Su 0 and references cited terein. Tis article extends te entropy conservative split form of Gassner et al 06b; Wintermeyer et al 06 to a new two-parameter family of wellbalanced and entropy conservative splittings. Moreover, SBP bases not including boundary nodes are considered and corresponding semidiscretisations are developed in one spatial dimension. To te autor s knowledge, tis general two-parameter family of numerical fluxes as not appeared in te literature before. Furtermore, positivity preservation in te framework of Zang and Su 0 is considered. Terefore, te new fluxes are investigated regarding positivity and finite volume FV subcells are introduced. Additionally, some known entropy stable and positivity preserving numerical fluxes are compared. Te extension to general SBP bases can be seen as some kind of positive negative result : It is possible to get te desired properties of te scemes using Gauß nodes, including a iger accuracy compared to Lobatto nodes, but tis comes at te cost of complicated correction terms. Terefore, altoug a similar extension to two-dimensional unstructured and curvilinear grids can be conjectured to exist, it is expected to be even more complex and tus not suited for ig-performance production codes. However, te combination of entropy stability and positivity preservation for te sallow water equations using ig-order SBP metods as not been considered before. Te results can be expected to be extendable to igperformance codes using Lobatto nodes and te flux-differencing form of Fiser and Carpenter 0 on two-dimensional unstructured and curvilinear grids if tensor product bases on quadrilaterals are used.

3 Entropy stable, well-balanced, positivity preserving metods for sallow water At first, some analytical properties of te sallow water equations are reviewed in section and te existing split form of Gassner et al 06b; Wintermeyer et al 06 is described in section. Afterwards, a new two-parameter family of entropy conservative numerical fluxes for te sallow water equations wit constant bottom topograpy is developed in section and extended to a varying bottom in section 5. Te corresponding semidiscretisation using general SBP bases is designed in section 6. Te positivity preserving framework of Zang and Su 0 is introduced to tis setting and numerical fluxes based on te entropy conserving scemes are investigated wit respect to positivity preservation in section 7. Additionally, some known fluxes are presented. An extension of te idea to use FV subcells to te setting of nodal SBP bases wit diagonal mass matrix is proposed in section and numerical experiments are presented in section 9. Finally, te results are summed up in section 0 and some conclusions and directions of furter researc are presented. Review of some properties of te sallow water equations Te sallow water equations in one space dimension are v 0 t + v x v + =, }{{} g g x b }{{}}{{} =u =fu =su,x were is te water eigt, v its speed, v te discarge, b describes te bottom topograpy, and g is te gravitational constant. In te following, some well known results tat will be used in te remainder of tis work are presented. As described inter alia by Boucut 00, Section., Dafermos 00, Section., Fjordolm et al 0; Wintermeyer et al 06, te entropy / total energy U = v + g + gb = u + gu + gu b, is strictly convex for positive water eigts > 0. Tus, in tis case, te associated entropy variables w = U g + b u = v v and conserved variables u can be used intercangeably. Wit te entropy flux u F = v + g v + gbv, smoot solutions satisfy t U + x F = 0, and te entropy inequality t U + x F 0 will be used as an additional admissibility criterion for weak solutions. If te bottom topograpy b 0 is constant, te entropy variables are g w = v. v

4 Hendrik Ranoca In tis case, te flux expressed in terms of te entropy variables b 0 is f uw w + w g w = w + w g w + g w + w = w w + w g g w + w w + 5, w and te flux potential b 0 is given by ψ = g w w + g w w + g w5 = g v, 6 fulfilling ψ w = f uw. Finally, te entropy Jacobian v w u = u w g g =, u w = g + u u u u v g + v g u u u g + = v v v 7 can be expressed by using a scaling of te eigenvectors in te form proposed by Bart 999, Teorem as w u = RR T, R = g v g v +, g were te columns of R are eigenvectors of te flux Jacobian f u. Tis scaling as also been used inter alia by Fjordolm et al 0, Section.. 5, Review of an existing split-form SBP metod In order to fix some notation, present te general setting and motivate te extensions developed in tis work, some existing results will be reviewed. A general SBP SAT semidiscretisation is obtained by a partition of te domain into disjoint elements. On eac element, te solution is represented in some basis, mostly nodal bases. Tese cells are mapped to a standard element for te following computations. Tere, te symmetric and positive definite mass matrix M induces a scalar product, approximating te L scalar product. Te derivative is represented by te matrix D. Interpolation to te two point boundary of te cell interval is performed via te restriction operator R and evaluation of te values at te rigt boundary minus values at te left boundary is conducted by te boundary matrix B = diag,. Togeter, tese operators fulfil te summation-by-parts SBP property M D + D T M = R T B R, 9 mimicking integration by parts on a discrete level Ω u xv + Ω xuv = uv Ω. Here, te notation of Ranoca et al 06, 05 as been used. Ten, similar to strong form discontinuous Galerkin metods, te semidiscretisation

5 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 can be written as te sum of volume terms, surface terms, and numerical fluxes at te boundaries. Gassner et al 06b proposed as semidiscretisation of te sallow water equations wit continuous bottom topograpy b in te setting of a discontinuous Galerkin spectral element metod DGSEM using Lobatto-Legendre nodes in eac element, tat can be generalised to diagonal norm SBP operators wit nodal bases including boundary nodes. Wintermeyer et al 06 extended tis setting to two space dimensions, curvilinear grids and discontinuous bottom topograpies. In one space dimension on a linear grid, tis semidiscretisation can be written as t = D v M R T B f num R v, t v = D v + v D v + v D v g D g D b M R T B f num R v gr v were e k is te k-t unit vector and for cell i + g { } [b] e 0 + g { } + [b] + e, 0 { } = i,0 + i,p, { } + = i,p + i+,0, [b] = b i,0 b i,p, [b] + = b i+,0 b i,p. Here, i,0, i,p are te values of at te first and last node 0, p in cell i, respectively. Using f num = { } {v }, f num v = { } {v } + g { }, as numerical surface flux, were {a } = a +a+, te resulting sceme. conserves te mass in general and te discarge for a constant bottom topograpy,. conserves te total energy wic is used as entropy,. andles te lake-at-rest stationary state correctly, i.e. it is conservative, stable and well-balanced, as proved by Wintermeyer et al 06, Teorem. Te split form discretisation as been recast into te flux differencing framework of Fiser et al 0; Fiser and Carpenter 0 using te translations provided by Gassner et al 06a, Lemma. Te resulting volume fluxes are f vol = {v }, f vol v = {v } {v } + g { } g { }, i.e. te volume terms in 0 can be rewritten using te differentiation matrix D as p vol ik, were fik = f vol u i, u k. k=0 D ikf vol

6 6 Hendrik Ranoca Te split form in 0 corresponds to te entropy conservative fluxes f vol. f num corresponds to a splitting, too, but excanging f num for f vol does not yield a well-balanced metod. Similarly, excanging f vol for f num as surface flux does not work properly, since te numerical flux and source discretisation ave to be coupled properly in order to result in a well-balanced discretisation, as described in section 5. Te remaining parts of tis paper are dedicated to te investigation of te following questions. Are tere oter entropy conservative fluxes tan f vol, f num and corresponding split forms? Section. Are tere oter discretisations of g x b tat can be used to get a wellbalanced sceme, respecting te lake-at-rest stationary state? Section 5. Can te split forms be used for a nodal SBP metod witout boundary nodes, e.g. for Gauß nodes? Section 6. Are tere entropy conservative / stable and positivity preserving numerical fluxes tat can be used to apply te bound preserving framework of Zang and Su 0? If so, is te resulting metod still entropy stable and well-balanced? Section 7 Entropy conservative fluxes and split forms for vanising bottom topograpy b 0 In tis section, several numerical fluxes and associated split forms of te sallow water equations wit constant bottom topograpy b 0 will be considered. Te numerical flux f num u, u + as to be consistent, i.e. f num u, u = fu. In order to be entropy conservative, te condition [w] f num u, u + = [ψ] of Tadmor 97, 00 as to be fulfilled in a semidiscrete setting. Similarly, if [w] f num u, u + [ψ], te numerical flux is entropy stable, since it contains more dissipation tan an entropy conservative flux. Here, [a] = a + a.. A two-parameter family of entropy conservative numerical fluxes In order to investigate te entropy conservation condition, te jump of te entropy variables w and te flux potential ψ ave to be expressed in a common set of variables. Here, te entropy variables w will be used. Tus, te jump of te flux potential ψ = g w w + g w w + g w5 6 can be written as [ψ] = g [w w ] + g [w w ] + g [w5 ]. 5

7 Entropy stable, well-balanced, positivity preserving metods for sallow water 7 Using a discrete analogue of te product rule [ab] = a + b + a b = a + + a b + b + a + a b + + b = {a }[b] + [a] {b }, 6 te first jump term [w w ] can be written in two different ways. Weigting bot variants wit weigts a R first line and a second line, respectively, results in [ww ] a = {w }[w ] + {w }[w ] = {w }[w ] + {w } {w }[w ] 7 a = {w w }[w ] + {w }[w w ] = {w w }[w ] + {w } [w ] + {w } {w }[w ] = + a {w } {w } + a {w w } [w ] + a {w } + a {w } [w ], for a R. Similarly, te second jump term [w w ] can be expressed in four different ways as [w w ] a = {w }[w ] + {w }[w ] = {w }[w ] + {w } {w }[w ] + {w } {w } [w ] a = {w }[w w ] + {w w }[w ] = {w } {w }[w ] + {w } {w }[w ] + {w w } {w }[w ] a = {w }[w w ] + {w w }[w ] = {w } {w }[w ] + {w } {w } [w ] + {w w }[w ] = {w }[w w ] + {w w }[w ] = {w } [w w ] + {w w } {w }[w ] + {w w }[w ] = {w } [w ] + {w } {w } [w ] + {w w } {w }[w ] + {w w }[w ] = a {w } + a + a {w } {w } + a a a {w } [w ] + a + a {w } {w } + + a a + a {w } {w } + a + a a {w w } {w } + a a {w w } [w ], were a, a, a R. However, expanding te terms in te last expression results in a {w } + a + a {w } {w } a + a + a {w } 9 = a + a + a w+ w+w w + w + w,

8 Hendrik Ranoca and a + a {w } {w } + a a + a {w } {w } 0 a +a + a {w w } {w } a + a {w w } = a + a + a w + w+ + w + w + w w+ w w. Tus, te expression depends only on a + a + a and can be simplified by setting a = a = 0 to [w w ] = a {w } + a {w } [w ] + a {w } {w } + + a {w } {w } + a {w w } {w } + a {w w } [w ] Finally, te last jump term can be expressed as [w 5 ] a 5 = {w }[w ] + {w }[w ] = {w }[w ] + {w } {w } [w ] a = 6 {w }[w ] + {w }[w ] = {w }[w ] + {w } [w ] + {w } {w }[w ] = {w }[w ] + {w } [w ] + {w } {w }[w ] + {w } {w }[w ] = {w }[w ] + {w }[w ] = = {w } {w }[w ] + {w } {w }[w ] + {w } [w ] a 5 + a 6 {w } + + a 5 a 6 {w } {w } + a 6 {w } + a 5 a 6 {w } {w } + a 5 a 6 {w } [w ], were a 5, a 6 R. However, tese parameters a 5, a 6 are also redundant, since te last expression can be simplified as [w 5 ] = w+ + w+w + w+w + w + w + w [w ] = {w } + {w } {w } [w ]. Inserting tese forms in te condition [w] f num u, u + = [ψ] for an entropy conservative flux, f a,a w, w + = + a {w } {w } + a {w w } + a g g g {w } + a g {w },

9 Entropy stable, well-balanced, positivity preserving metods for sallow water 9 f a,a v w, w + = a g {w } + a {w } + a g g {w } {w } + + a g + a g {w } {w } + a g {w w } {w } {w w } + g {w } + g {w } {w }, turns out to be a two-parameter family of entropy conservative numerical T fluxes expressed in terms of te entropy variables w = g v, v for vanising bottom topograpy b 0. Expanding te terms in entropy and primitive variables, respectively, results after direct but tedious calculations in f a,a = a w + w + + w w + + a g g + + a 6g f a,a v = + a g and + 7 a 6g w+ + w + a 6g w + + w + a w + w g w + w + + w w w + w + w w + 5 w+w + w + w, + + a 6g + w +w + w + w w + w g g + w+ + w g +w + w+w + w + w + w, w + w + w w+ + v + + v + + a + v + v + 6 v+ v+v v + v + v, f a,a = a + a + a 6g f a,a v = + a g a g + a + a a + a + v + v+ + +v + v 6 + a + a v+ v g +v + v. Tis proves + v + v + v+ + v Lemma Te two-parameter family of numerical fluxes f a,a, expressed also as 5 and 6, is a family of consistent and entropy conservative numerical fluxes for te sallow water equations wit vanising bottom topograpy b 0. To te autor s knowledge, tis family of entropy conservative fluxes as not been considered en bloc in te literature before. Instead, only te fluxes corresponding to tree distinct coices of te parameters a, a ave been proposed, as described in te Remarks and in te following section.

10 0 Hendrik Ranoca. Relations wit known fluxes and metods Te crucial ingredient to obtain te entropy conservative fluxes in Lemma as been te expression of bot te entropy variables w and te flux potential ψ as polynomials in te same variables. Terefore, it may be conjectured, tat if suc a condition is complied wit, tere are entropy conservative fluxes expressed in terms of tese variables, corresponding to a split form as described inter alia by Gassner et al 06a and in section 6. Remark Te general entropy conservative flux of Tadmor 97, Equation.6a, obtained by integration in pase space, can be recovered by te family of entropy conservative fluxes. Indeed, for f uw as in 5, f u sw + sw + ds 7 0 = w + w + + w +w + w w + + w w g w + + w +w + w +w + w, and 0 f v u sw + sw + ds = w+ + w + w + 9w + w+ + 6w + w + w + w + w g + w + w w+ + 6w w + w + 9w w + w+ + w+w + w+w + w + w + w. Comparing tis wit te numerical fluxes f a,a 5, it can be seen tat Tadmor s flux as above is recovered by setting a = a =. Remark Using entropy variables, tere does not seem to be a special coice of te parameters a, a for te numerical flux 5. However, te expression in primitive variables, v reveals tat te coice a = a is special, since iger order terms in te velocity tat are consistent wit zero are removed. In tis case, te one-parameter family of entropy conservative fluxes can be written as f a = a {v } + + a { } {v }, f a v = a {v } {v } + + a or { } {v } + g a { } + a g { }, 9 f a = a + v + + v + + a + v + v +, 0

11 Entropy stable, well-balanced, positivity preserving metods for sallow water f a v = + a g + + a a + v + v+ g + + a 6 + +v + v + v + v. + v+ + v Remark Te numerical surface and volume fluxes f vol, f num of Gassner et al 06b and Wintermeyer et al 06 are members of tis family wit parameters a = and a =, respectively. Terefore, tis oneparameter family of entropy conservative fluxes 9 can also be seen as linear combinations of te two fluxes and wit coefficients summing up to one. Remark 5 Te derivations of section. can also be conducted using primitive variables, v instead of te entropy variables w. In tis case, only te oneparameter family of entropy conservative numerical fluxes 9, 0 can be obtained. Remark 6 As proved by Fiser and Carpenter 0, Teorem., a twopoint entropy conservative flux as f a,a can be used to construct a igorder spatial discretisation for diagonal-norm SBP operators wit nodal basis including boundary nodes. Gassner et al 06a, Lemma provided some examples for analogous split forms and numerical fluxes given by simple products of mean values. Analogously, using a diagonal-norm SBP derivative operator D i.e. an SBP derivative operator D, were te corresponding norm / mass matrix M is diagonal wit nodal basis including boundary nodes, te split form corresponding to te flux f a,a expressed via primitive variables 6 is for te component given by = = = [ VOL a,a p k=0 ]i = p k=0 D i,k f a,a u i, u k a D i,k i v i + k v k + + a i v k + k v i p D i,k k=0 [ + a + a 6g vi vi v k v i vk + vk a k v k + + a i v k + k v i + a + a g v k v i v k v i v k a D v + + a D v + v D + a + a ] D v v D v v D v. g i

12 Hendrik Ranoca Here, some summands ave been dropped, because te derivative is exact for constants, i.e. D = 0, resulting in p k=0 D i,k = 0. Te first two terms form a consistent discretisation of x v = x v+ x v for smoot solutions. Te tird term is consistently zero, since x v = x v v + v x v by te product rule for smoot solutions. Similarly, te v component can be computed via = = = [ VOL a,a v p D i,k k=0 p D i,k k=0 [ + a ]i = p k=0 + a g D i,k f a,a v u i, u k i + k a + a iv i v k + a + kv i v k + a g i k i v i + k v k + a + a g + a + a 6 v i v i v k + v k g k + a g i k a + a 5 k vk + a + a i vk + k vi + iv i v k + a + a 6g vk vi vk gd + a g D a + a 5 D v + a + a D v + v D + ] D v v D v i, + a + a 6g + kv i v k v D v + v D v i vk + k vi were again p k=0 D i,k = 0 as been used. Te first two terms form a consistent discretisation of g x, te tree following terms are consistent wit x v and te last two terms are consistently zero. Te entropy conservation follows from te general result of Fiser and Carpenter 0, Teorem. and will also be investigated in more detail for Gauß nodes and oter SBP bases in section 6. Remark 7 It is also possible to start wit a general ansatz of te split form and use conditions for consistency, conservation, and entropy stability similar to section 6, in order to determine te coefficients. Tis yields te same twoparameter family of fluxes and corresponding split forms, but is muc more tedious.

13 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 Adding well-balanced source discretisations In tis section, te discretisation of te source term g x b in te sallow water equations will be investigated. It sould be consistent, stable, and well-balanced, if combined wit te remaining semidiscretisation derived in te previous section. 5. Connections between finite volume and SBP SAT scemes A general semidiscretisation of a conservation law t u + x fu = 0 wit a polynomial SBP metod using te notation of Ranoca et al 06 can be written as t u = VOL + SURF M R T B f num, were VOL contains te volume terms consistent wit x fu, possibly using some split form, f num is te numerical surface flux, and SURF contains additional surface terms, consistent wit te difference of te flux values fu at te boundaries, i.e. M R T B R f in te simplest case, but additional terms may also appear, especially if a nodal basis witout boundary nodes is used, see also section 6. If te polynomial degree p is set to zero, te volume terms vanis, since te derivative is exact for constants, i.e. D = 0. Additionally, since te extra surface terms SURF are a consistent evaluation of te difference of boundary values, tey vanis, too, because tis difference is zero for constants. Terefore, tis metod reduces to a simple finite volume metod. If te cell i is of size x and te flux f num between te cells i and k is denoted as fi,k num, te FV metod can be written as t u i = fi,i+ num fi,i num, x and is determined solely by te numerical flux f num used at te boundaries. On te oter and, using te teory of Fiser and Carpenter 0, a finite volume metod wit entropy conservative flux f num can be used to construct te volume terms VOL, if a nodal SBP basis including boundary nodes is used. In tis case, since te evaluation at te boundary is exact and commutes wit nonlinear operations, te surface terms are simply M R T B R f. Tis strong correspondence between SBP scemes and FV metods will be used in te following sections to extend results from one area to te oter and vice versa. 5. Extended numerical fluxes and entropy conservation In order to incorporate source terms in a finite volume metod, te additional contributions can be incorporated into te numerical fluxes, resulting in extended numerical fluxes as described inter alia in te monograp by Boucut

14 Hendrik Ranoca 00, Section and references cited terein. Here, te extended flux for te discarge v will ave te form f num,ext i,k = f num i,k + S i,k, 5 were fi,k num is a usual entropy conservative and symmetric numerical flux of te problem witout source terms and S i,k describes te source terms and is in general not symmetric. Remark Setting te polynomial degree p to zero, te surface flux fv num and te surface source terms of te SBP SAT semidiscretisation of Wintermeyer et al 06 can be written as te FV metod of Fjordolm et al 0, wic can be written using te extended numerical flux f num v i,k + g { } i,k [b] i,k, f num v i,k = { } i,k {v } i,k + g { } i,k, 6 were {a } i,k = ai+a k, [a] i,k = a k a i. Rewriting te FV evolution equation by adding f i f i = 0 motivated by te form of SBP SAT metods and using extended numerical fluxes yields t u i = f num,ext i,i+ f i f num,ext i,i f i. 7 x Terefore, te rate of cange of te entropy U can be calculated as w i f num,ext i,i+ f i w i f num,ext i,i f i t U i =w i t u i = x [ = w i x ] f num,ext i,i+ f i + F i [ ] w i f num,ext i,i f i + F i. Tus, adding te contributions of te rigt and side of cell i and te left and side of cell i + yields after multiplication wit x [ w i+ f num,ext i+,i f i+ =w i+ w i f num i,i+ + w i+ S i+,i w i S i,i+, [ ] w i f num,ext i,i+ f i + F i [w i f i + F i ] }{{} =ψ i+ ] + F i+ [w i+ f i+ + F i+ ] }{{} =ψ i 9 were te extended flux f num,ext 5 as been inserted and te symmetry of te numerical flux f num as been used. Assume now tat te numerical flux fi,i+ num is cosen as an entropy conservative one, fulfilling [w] i,i+ fi,i+ num = [ψ] i,i+, for vanising bottom T topograpy b 0. Here, te entropy variables are w = g v, v, since T b 0. In te general case, te entropy variables are w = g + b v, v,

15 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 resulting in [w] i,i+ f num i,i+ = [ψ] i,i+ + gf num i,i+ [b] i,i+. Tus, te contribution of one boundary to te rate of cange of te entropy 9 is Tis proves gf num i,i+ [b] i,i+ + w i+ S i+,i w i S i,i+! = 0. 0 Lemma 9 If te source discretisation S i,k in te extended numerical flux 5 is cosen suc tat te expression 0 is zero for an entropy conservative numerical flux f num for te sallow water equations wit vanising bottom topograpy b 0, ten te resulting sceme is entropy conservative for general bottom topograpy. To te autor s knowledge, tis general result as not been formulated before. Instead, special source discretisations ave been cosen and similar calculations adapted to te specific discretisation ave been performed. Remark 0 Te source discretisation of Fjordolm et al 0, Lemma. and Wintermeyer et al 06 results in te extended numerical flux 6 wit source terms S i,k = g { } i,k [b] i,k. Tus, inserting tis and teir numerical flux f num i,k = { } i,k {v } i,k into 0 results in g { } i,i+ {v } i,i+ [b] i,i+ + v i+ g { } i,i+ [b] i+,i v i g { } i+,i [b] i,i+ =g { } i,i+ {v } i,i+ [b] i,i+ g { } i,i+ v i + v i+ fulfilling te condition of Lemma 9. [b] i,i+ = 0, Remark Te discretisation of te volume terms for te discarge of Gassner et al 06b and Wintermeyer et al 06 can be formulated using te flux difference form and an extended numerical flux f vol v i,k + g i [b] i,k, f vol v i,k = {v } i,k {v } i,k + g { } i,k g { } i,k. Indeed, since D = 0, te discretisation of g x + g x b is given by =g =g p k=0 p k=0 D i,k g { } i,k g { } i,k + g i [b] i,k D i,k i + k p D i,k i k + i b k b i = g i k=0 i + k + ib k b i p D i,k k + b k g i b i k=0 } {{ } =[D +b] i p k=0 D i,k }{{} =[D ] i=0.

16 6 Hendrik Ranoca In te same way te entropy conservative numerical fluxes and can be combined to get te one-parameter family of entropy conservative fluxes f a 9 for te sallow water equations wit vanising bottom topograpy b 0, te extended numerical fluxes 6 and can be combined to get entropy conservative extended numerical fluxes for te sallow water equations wit general bottom topograpy. For te two-parameter family, te source discretisation S i,k in te extended numerical flux 5 as to be adapted to te additional terms wit a in order to fulfil te condition 0 of Lemma 9. Since te two-parameter flux 6 for contains an additional term a+a 6g v + v+v v + v + v compared to te one-parameter flux 0, te new source term S i,k can be written as te sum of te source term g a i + +a k b k b i for te one-parameter flux 0 and an additional source term S i,k, obeying gf num i,k [b] i,k + w k S k,i w i S i,k = a + a vi vi v k v i vk + vk! b k b i + v k Sk,i v i Si,k = 0. 6 Tis can be rewritten using vi v i v k v i vk + v k = v i + v k v i v k. Tus, coosing S i,k = a+a 6 v k v i b k b i results in te desired equality a + a vi vi v k v i vk + vk b k b i + v k Sk,i v i Si,k 5 6 = a + a v i + v k v i v k b k b i + v k v i v k b i b k 6 v i v k v i b k b i =0. 5. Well-balancedness preserving te lake-at-rest steady state Te extended numerical flux 6 of Fjordolm et al 0 and Wintermeyer et al 06 preserves te lake-at-rest steady state + b const, v = 0 in a finite volume metod. Similarly, te extended numerical flux is well-balanced, since for +b const and v = 0 g f vol v i,,i+ + g i [b] i,i+ g f vol v i,,i + g i [b] i,i = { } i,i+ { } i,i+ + i [b] i,i+ { } i,i + { } i,i i [b] i,i = i i+ + ib i+ i i ib i = 0. 6 Comparing te two-parameter family of numerical fluxes f a,a 6 wit te one-parameter family f a 0, te additional parameter a contributes

17 Entropy stable, well-balanced, positivity preserving metods for sallow water 7 only to terms containing te velocity v. Tus, it is irrelevant for well-balancing, since tese terms vanis for te lake-at-rest initial condition. Togeter wit te results of te previous subsection, tis proves Lemma Te two-parameter family f a,a i,k = a + a + a 6g f a,a v i,k = + a g i + k i v i + k v k + + a i v k + k v i vi vi v k v i vk + vk, + a g i k a + a a + a i vk + k vi + iv i v k + kv i v k 6 + a + a g + g a v i v i v k + v k i + + a k b k b i, i vi + k vk + a + a v k v i b k b i 6 7 is a family of entropy conservative and well-balanced extended numerical fluxes including contributions of te source term g x b for te sallow water equations wit general bottom topograpy b. To te autor s knowledge, tis general family of entropy conservative and well-balanced extended numerical fluxes as only been considered for te special coices a = a, a {, } before. Remark Again, a one-parameter family is given by te special coice a = a as f a i,k = a i v i + k v k + + a i v k + k v i, v i,k = + a g i + k + a g i k + a i vi + k vk a i vk + k vi + iv i v k + kv i v k 6 f a + g a i + + a k b k b i. Remark Te volume terms corresponding to te two-parameter family of extended fluxes 7 are VOL a,a VOL a,a v = a D v + + a + a + a g = + a gd + a D v + v D D v v D v v D v, g D a + a 5 D v 9

18 Hendrik Ranoca were = p k=0 + a + a D v + v D + v D v + v D v + a + a D v v D v + a g D b 6g + + a g D b b D + a + a D bv b D v v D bv + v D b + bv D v D i,k g a i + + a k b k b i 50 p D i,k k=0 a + a + a v k v i b k b i 6 g i b k + + a g k b k b i vkb k v i v k b k + vi b k vkb i + v i b i v k + a + a [ a = g D b + + a + a + a g D b b D D bv v D bv + v D b b D v + bv D v as been used. Te corresponding surface terms using nodal bases including boundary nodes are simply given by SURF = M R T B R v, SURF v = M R T B R v + gr. 5 Te numerical fluxes used for te volume discretisation and as surface flux may be combined arbitrarily, as done by Gassner et al 06a; Wintermeyer et al 06, were tey ave used f, as volume flux and f, as surface flux. Tus, a general semidiscretisation is of te form ] i, wit a, a, b, b R. t = VOL a,a t v = VOL a,a v + SURF M R T B f b,b, + SURF v M R T B f b,b, v 5 Remark 5 Te volume terms 9 vanis for te lake-at-rest initial condition + b const, v = 0. Since v = 0, tis is immediately clear for VOL a,a. For te discarge term VOL a,a v, using + b = + bi, + b const = 0, + a D + a D + a D b + + a D b b D 5

19 Entropy stable, well-balanced, positivity preserving metods for sallow water 9 = + a = + a D + b + b = a D + b = 0. + a D + a D + a D + a D b + a b D D b + a b D 6 Extension to general SBP bases In tis section, an extension of te previous result to a nodal DG metod using Gauß nodes instead of Lobatto nodes or more general SBP bases will be investigated. Altoug te volume terms 9 ave been derived in section 5 wit te assumption of a diagonal-norm SBP basis including boundary nodes, tey can be easily transferred to te setting of a general SBP basis. If te multiplication operators are self-adjoint wit respect to te scalar product induced by M, e.g. for a nodal basis wit diagonal mass matrix, ten te same volume terms 9 can be used. Oterwise, some multiplication operators a ave to be replaced by teir M -adjoints a = M a T M, as proposed by Ranoca et al 05. Tis results in te volume terms VOL a,a = a D v + + a + a + a g VOL a,a D v + v D D v v D v v D v v = + a gd + a g D a + a 5 D v + a + a D v + v D + v D v + v D v + a + a D v v D v + a g D b 6g + + a g D b b D + a + a, D bv b D v v D bv + v D b + bv D v 5 However, te surface terms 5 also ave to be adapted to a general basis. Often, te split form of te volume terms is described as some correction for te product rule tat does not old discretely. However, as described by Ranoca et al 06, it is te multiplication tat is not correct on a discrete level, resulting in an invalid product rule. Moreover, if no boundary nodes are included in te basis, tis inexactness also as to be compensated in te surface terms. Tus, in te same spirit as te split form of te volume terms can be seen as corrections to inexact multiplication, some kind of correction as to be used for te interpolation to te boundaries..

20 0 Hendrik Ranoca Investigating conservation across elements, te time derivatives of te conserved variables 5 are multiplied wit T M, corresponding to integration over an element. Tis yields for te volume terms 5 T M VOL a,a T M D v + + a T M D v + v D 55 = a SBP + a + a T M g = a + a + a g T M D v + + a D v v D v v D v T M D v + v T M D T M D v v T M D v v T M D v a = T R T B R v + + a T R T B R v + a + a T R T B R v v T R T B R v. g Here, =, te SBP property M D = R T B R D T M, and D = 0 ave been used. If multiplication and restriction to te boundary commute, tese volume terms are simply T R T B R v and yield te desired integral form. Similarly, T M VOL a,a v 56 g T M D + a g T M D a + a 5 T M D v = + a + a + a + a + a T M 6g + + a g T M + a + a T M T M D v + v D + T M v D v + v D v D v v D v + a g T M D b D b b D D bv b D v v D bv + v D b + bv D v = + a g T M D + a g T M D a + a 5 T M D v + a + a T M D v + v T M D + v T M D v + v T M D v + a + a T M D v v T M D v + a g T M D b 6g + + a g T M D b b T M D + a + a T M D bv b T M D v a + a v T M D bv v T M D b bv T M D v

21 Entropy stable, well-balanced, positivity preserving metods for sallow water SBP + a = g T R T B R + a g T R T B R T R T B R v + a + a a + a 5 + vt R T B R v + a + a { + g T M D b + + a g + a + a 6g a + a T R T B R v [ T R T B R v v T R T B R v ] } { } v T M D bv v T M D b bv T M D v [ ] T R T B R b b T R T B R [ T R T B R bv b T R T B R v ]. Here, te terms in squared brackets [ ] vanis if restriction to te boundary and multiplication commute, i.e. for a basis using Lobatto nodes. However, for oter bases using e.g. Gauß nodes, tese contributions are not zero in general. Te first term in curly brackets { } is a consistent discretisation of te source term g x b. Te second term in curly brackets { } vanises, if te product rule is valid, e.g. for constant bottom topograpy b. However, for general bottom topograpy, it is not of te desired form for te source influence g x b and it migt be better to set it to zero by coosing a = a, i.e. only te one-parameter family instead of te two-parameter family. Tese surface terms obtained in 55 and 56 ave to be balanced by te surface terms of te SBP SAT semidiscretisation 5 in order to get te desired result T M t = T R T B f b,b, T M t v = T R T B f b,b v + consistent contribution of g x b, 57 leading to a conservative sceme. Turning to stability, te approximation of t U = U u t u = w t u, w = influenced by te volume terms is given by g + b v, v T, 5 w T M VOL a,a + w T M VOL a,a v 59 = a g + b T M D v + + a g + b T M D v + v D + a + a + b T M D v v D v v D v a v T M D v + a v T M D v + v D a + a v T M D v v D v v D v 6g

22 Hendrik Ranoca + + a + a + a + a + a v T M 6g + + a gv T M + a + a v T M = a gv T M D + a gv T M D a + a 5 v T M D v v T M D v + v D + vt M v D v + v D v D v v D v + a gv T M D b D b b D D bv b D v v D bv + v D b + bv D v T + b M D v g + b T M D v + + a g + + a g v + bv T M D + a + a a + a a a + a 6g + a + a 6g a + a 5 v + bv T M D v a + a v T M D v + a + b T M D v v + bv T M D v v T M D v + a v T M D v T M D v + a + a v T M D v 6g v T M D v + + a gv T M D + a gv T M D v T M D v + a + a v T M D v + a + a v T M D + vt M D v + vt M D v + a + a 6g + a + a a + a + a + a v T M D v a + a v T M D v g gv T M D b + + a gv T M D b gbv T M D + a + a v T M D bv bv T M D v a + a v T M D bv v T M D b + a + a bv T M D v = a g T M D v + v T M D + a g b T M D v + v T M D b + + a g T M D v + v T M D + + a g b T M D v + v T M D b + 0 bv T M D + a + a T M D v + v T M D + a + v T M D v + v T M D v

23 Entropy stable, well-balanced, positivity preserving metods for sallow water SBP + 5 a a v T M D v + v T M D v a + a v T M D v + v T M D v 6g + a + a 6g + a + a a + a + a + a a = + + a + a + v T M D v + v T M D v bv T M D v + v T M D bv bv T M D v + v T M D bv b T M D v + v T M D b g T R T B R v + a a + a 6g + a + a gb T R T B R v + + a gv T R T B R gv T R T B R b + a + a T R T B R v v T R T B R v + 5 a a v T R T B R v + a + a b T R T B R v. v T R T B R v + a + a v T R T B R v 6g v T R T B R bv a + a bv T R T B R v If multiplication and restriction to te boundary commute, tese terms simplify to T R T B R F, were F = g v + gbv + v is te entropy flux. Tese surface contributions resulting from te volume terms ave to be balanced by te surface terms w T M SURF + w T M SURF v, in order to get an estimate of te form [w] f num [ψ] for te entropy cange influenced by one boundary node if te bottom topograpy is continuous across elements. Tat is, te simple interpolations M R T B R v, M R T B R v, M R T B R, 60 in te surface terms 5 for te metod including boundary nodes ave to be adapted. Te following combination of surface term structures proposed by Ranoca et al 06; Ortleb 06 will be investigated SURF a,a 6 =b M R T B R v + b M R T B R R v + b M R T B R v + b v M R T B R + c g M R T B R v + c g M R T B R vr v + c g v M R T B R v + c g v M R T B R v,

24 Hendrik Ranoca SURF a,a v =d M R T B R v + d M R T B R R v + d M R T B R vr v + d M R T B R R v + d 5 v M R T B R v + d 6 v M R T B R v + d 7 v M R T B R + d M R T B R v + e gm R T B R + e gm R T B R + e g M R T B R + k g M R T B R v + k g M R T B R vr v + k g M R T B R v + k g M R T B R v R v + k 5 g M R T B R v + k 6 g v M R T B R v + k 7 g v M R T B R vr v + k g v M R T B R v + k 9 g v M R T B R v + k 0 g v M R T B R v + k g v M R T B R v + l M R T B R bv + l M R T B R br v + l M R T B R bvr v + l M R T B R br v + l 5 b M R T B R v + l 6 v M R T B R b + l 7 bv M R T B R v + l v M R T B R bv + l 9 b M R T B R v + l 0 v M R T B R br v + m gm R T B R b + m gm R T B R br + m g M R T B R b + m gb M R T B R v M R T B f a,a + M R T B f a,a R v, were b i, c i, d i, e i, k i, l i, m i R are free parameters tat ave to be determined. Considering conservation for, te relevant conditions are obtained by multiplying te surface terms wit T M. T M SURF a,a 6 =b T R T B R v + b T R T B R R v + b T R T B R v + b v T R T B R + c g T R T B R v + c g T R T B R vr v + c g vt R T B R v + c g vt R T B R v =b T R T B R v + b + b + b T R T B R R v + c g T R T B R v + c + c + c T R T B R vr v.

25 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 Here, some manipulations as T R T B R v = R v R L v L = T R T B R R v proposed by Ranoca et al 06 ave been used. Tus, comparison wit 55 yields te conditions b = a c = a + a Similarly, for v,, b + b + b = + a, 6, c + c + c = a + a. T M SURF a,a v 6 =d T R T B R v + d T R T B R R v + d T R T B R vr v + d T R T B R R v + d 5 v T R T B R v + d 6 v T R T B R v + d 7 v T R T B R + d T R T B R v + e g T R T B R + e g T R T B R + e g T R T B R + k g T R T B R v + k g T R T B R vr v + k g T R T B R v + k g T R T B R v R v + k 5 g T R T B R v + k 6 g vt R T B R v + k 7 g vt R T B R vr v + k g vt R T B R v + k 9 g vt R T B R v + k 0 g vt M R T B R v + k g vt M R T B R v + l T R T B R bv + l T R T B R br v + l T R T B R bvr v + l T R T B R br v + l 5 b T R T B R v + l 6 v T R T B R b + l 7 bv T R T B R v + l v T R T B R bv + l 9 b T R T B R v + l 0 v T R T B R br v + m g T R T B R b + m g T R T B R br + m g T R T B R b + m gb T R T B R. Analogously, comparing tis wit 56 results in te conditions d = 5 a a, d + d 7 + d = a + a, 65 d + d 5 + d 6 =, d =0, e = + a, e + e = a, k = a + a, k + k 6 + k =0, 6 k + k 9 = a + a, k + k 7 + k 0 =0, 6 k 5 + k =0, l = a + a,

26 6 Hendrik Ranoca l + l 5 + l 6 = a + a, l + l 7 + l =0, l + l 9 + l 0 =0, m = + a, m + m + m = + a. Considering stability, te surface terms 6 yield w T M SURFa,a + w T M SURFa,a v 66 =b g T R T B R v + b g T R T B R R v + b g T R T B R v + b gv T R T B R + c T R T B R v + c T R T B R vr v + c v T R T B R v + c v T R T B R v + b gb T R T B R v + b gb T R T B R R v + b gb T R T B R v + b gbv T R T B R + c b T R T B R v + c b T R T B R vr v + c bv T R T B R v + c bv T R T B R v b v T R T B R v b v T R T B R R v b v T R T B R v b v T R T B R c g vt R T B R v c g vt R T B R vr v c g vt R T B R v c g vt R T B R v + d v T R T B R v + d v T R T B R R v + d v T R T B R vr v + d v T R T B R R v + d 5 v T R T B R v + d 6 v T R T B R v + d 7 v T R T B R + d v T R T B R v + e gv T R T B R + e gv T R T B R + e gv T R T B R + k g vt R T B R v + k g vt R T B R vr v + k g vt R T B R v + k g vt R T B R v R v + k 5 g vt R T B R v + k 6 g vt R T B R v + k 7 g vt R T B R vr v + k g vt R T B R v + k 9 g vt R T B R v + k 0 g vt M R T B R v + k g vt M R T B R v + l v T R T B R bv + l v T R T B R br v + l v T R T B R bvr v + l v T R T B R br v + l 5 bv T R T B R v + l 6 v T R T B R b + l 7 bv T R T B R v + l v T R T B R bv + l 9 bv T R T B R v + l 0 v T R T B R br v + m gv T R T B R b + m gv T R T B R br + m gv T R T B R b

27 Entropy stable, well-balanced, positivity preserving metods for sallow water 7 + m gbv T R T B R vt R T B f a,a + vt R T B f a,a R v = b + b + e g T R T B R v + b g T R T B R R v + b + e g T R T B R v + c b + d 7 T R T B R v + c b + d T R T B R vr v + c b + d 5 v T R T B R v + c b + d + d 6 v T R T B R v + b + m gb T R T B R v + b + m gb T R T B R R v + b + m gb T R T B R v + b + m gbv T R T B R + c + l 6 b T R T B R v + c + l + l 0 b T R T B R vr v + c + l 5 + l bv T R T B R v + c + l + l 7 bv T R T B R v + c c + k 6 + k 9 g vt R T B R v + c + k + k 7 g vt R T B R vr v + c + k + k g vt R T B R v + d v T R T B R vr v + d v T R T B R R v + d v T R T B R v + e gv T R T B R + k + k 0 g vt R T B R vr v + k + k g vt R T B R v R v + k 5 g vt R T B R v + l + l 9 v T R T B R bvr v + l v T R T B R br v vt R T B f a,a + vt R T B f a,a R v. Comparing tis wit 59 yields te conditions b + b + e = a, b + e =, 67 b + e = + a, c b + d 7 = a + a, c b + d =0, c b + d 5 + d = a +, c b + d + d 6 = 5 a a, b + m = a, b + m =0, b + m = a +,

28 Hendrik Ranoca c c b + m =0, c + l 6 = a + a, c + l + l 0 =0, c + l 5 + l = a + a, c + l + l 7 = a + a, l =0, + k 6 + k 9 = a + a, 6 c + k + k 7 =0, c + k + k = a + a, 6 d =0, d =0, k + k 0 =0, k + k =0, k 5 =0 l + l 9 =0, l =0. Solving te linear system given by 6, 65, and 67 wit SymPy SymPy Development Team 06 results in te free parameters m, k 9, k 0, k, l 0 for any given parameters a, a : b = a +, b = a + m +, 6 b =0, b = m, c = a + a, c = a a k 0 k 9 +, c =k 0 k + k 9, c =k, d = a a + 5, d = a + a + k 0 + k 9 + m, d =0, d =0, d 5 =k +, d 6 = k, d 7 = m, d = k 0 k 9, e = a +, e = a m +, e =m, k = a 6 + a 6, k = k 0, k = a 6 a 6 k 9 +, k =0, k 5 =0, k 6 =k 0 k, k 7 = k 0, k =0, l = a + a, l =l 0, l = a + a + k 0 + k 9 l 0, l =0, l 5 =l 0 a + a k 0 k 9,

29 Entropy stable, well-balanced, positivity preserving metods for sallow water 9 l 6 =0, l 7 = k, l =k l 0, l 9 = l 0, m = a +, m = a m, m =0. Tis proves te following Teorem 6 For a, a, α, α R, using a general SBP operator, te semidiscretisation t = VOL a,a t v = VOL a,a v + SURF a,a + SURF a,a v M R T B f α,α, M R T B f α,α, v 69 wit volume terms 5 and surface terms 6, were te parameters are cosen according to 6 wit free parameters m, k 9, k 0, k, l 0 R,. conserves te total mass. Additionally, it conserves te total momentum v, if te bottom topograpy is constant. Oterwise, te rate of cange is consistent wit te source term g x b.. conserves te total entropy / energy U.. andles te lake-at-rest condition correctly. Tat is, tis semidiscretisation is conservative across elements, stable entropy conservative, and well-balanced. Remark 7 For an SBP operator using a nodal basis wit diagonal mass matrix M, te volume terms can be equivalently expressed in te flux difference form corresponding to te extended numerical fluxes 7 in Lemma. Remark For te special case a = a, i.e. te one-parameter family, te volume terms 5 can be considerably simplified. Accordingly, te ansatz 6 can be simplified by setting te coefficients c i, k i, and l i to zero. In tis case, only m remains as free parameter. In tis case, te surface terms 6 become

30 0 Hendrik Ranoca SURF a, a = a M R T B R v + a + + m M R T B R R v m v M R T B R, SURF a, a v = a M R T B R v + a + + m M R T B R R v + v M R T B R v m v M R T B R + a + gm R T B R + a m gm R T B R + m g M R T B R + a + gm R T B R b a + + m gm R T B R br + m gb M R T B R v M R T B f a,a + M R T B f a,a R v. 70 Here, te coice m = 0 results in te fewest terms. Additionally, coosing a = cancels some oter terms and results in te skew-symmetric form of Gassner et al 06b. Remark 9 In tis one-dimensional setting using a linear grid, te surface correction terms are more complicated tan te corresponding terms for Burgers equation Ranoca et al 06. Tus, altoug an extension to Gauß nodes on two-dimensional curvilinear grids may seem possible, it would result in a vast amount of surface and volume terms. As described in an earlier version by Wintermeyer et al 05, te use of te flux differencing form is clearly superior for te implementation. Altoug te usage of Gauß nodes may yield improved accuracy, te additional correction terms render te metod inferior compared to te usage of Lobatto nodes and te flux differencing form. Remark 0 Te two-parameter family of fluxes 5 as been derived in entropy variables w and translated to primitive variables, v 6 in te following calculations. Hence, it may seem natural to consider a splitting similar to 5 and 6 but expressed using entropy instead of primitive variables. Te volume terms can be obtained similarly to 5, but te surface terms are more delicate. A general ansatz similar to 6 witout te correction v M R T B f a,a + M R T B f a,a R v as not been successful, i.e. te resulting linear system was not solvable. However, it migt be possible to get te desired result using anoter ansatz for te surface terms. 7 Numerical surface fluxes and positivity preservation In tis section, several numerical surface fluxes will be presented, wic can be used to get entropy stable, positivity preserving, and well-balanced scemes.

31 Entropy stable, well-balanced, positivity preserving metods for sallow water Terefore, te framework of positivity preservation by Zang and Su 0 will be presented and adapted to te setting of a nodal SBP metod. 7. Positivity preservation Te argumentation of Zang and Su 0 can be summarised as. Since te metod sould be conservative, ensure a non-negative mean value of in eac cell.. If te eigt becomes negative somewere but te mean value is nonnegative, use a suitable limiter to enforce te non-negativity as needed. In a pure finite volume framework using an explicit Euler step, in cell i + i = i t f num u i, u i+ f num u i, u i. 7 x If tere are numerical fluxes suc tat under a suitable CFL condition t c x te non-negativity can be guaranteed, tis can be extended to te DG setting as proved by Zang and Su 0 and references cited by tem as follows. Coose q + Lobatto-Legendre nodes x k in te element i. Since Lobatto quadrature wit weigt ω k at x k is exact for polynomials of degree q, ensure q p, i.e. q p. If te finite volume metod 7 using te numerical flux f num is positivity preserving under te CFL condition t c x, te DG metod is positivity preserving under te scaled CFL condition t ω q c x, if te values of i,k at te quadrature nodes x 0,..., x are nonnegative. Ten, after te Euler step, te mean value in eac cell i is ensured to be non-negative. Applying te simple linear scaling limiter i, min i 0, i := θ i + θ i = i + θ i i, θ = i, min i min i < 0, i 7 of Liu and Oser 996, Section. ensures min i 0. Here, min i can be computed as te minimum of te polynomial i of degree p in te complete cell, or as te minimum of i at te nodes x 0,..., x q used to guarantee te nonnegativity of te cell mean. Applying tis limiter in eac Euler sub-step, tis procedure can be extended to SSP metods consisting of convex combinations of Euler steps. In a nodal collocation framework suc as nodal DG or FD, it would be more natural to enforce te non-negativity of te water eigt at te collocation nodes ξ 0,..., ξ p in te standard element. Tus, if Lobatto nodes are used, te framework of Zang and Su 0 described above extends directly, were te non-negativity of te mean value + i can be guaranteed under te possibly worse CFL condition t ω p c x, were p q, wen te limiter

32 Hendrik Ranoca 7 is applied wit min i = min { i ξ 0,..., i ξ p }. Alternatively, to get a better CFL condition t ω q c x, te water eigt can be interpolated to te oter nodes x 0,..., x q and te limiter 7 can be applied wit min i = min { i x 0,..., i x q, i ξ 0,..., i ξ p }. Similarly, if Gauß nodes are considered, only one possibility is obvious: Interpolate to suitable Lobatto nodes x 0,..., x q and use te limiter 7, again wit te coice of te minimum as minimum over solution nodes and interpolation points min i = min { i x 0,..., i x q, i ξ 0,..., i ξ p }. If tis limiter sould be coupled wit an entropy conservative / stable metod, a natural question is weter te limiter 7 is entropy stable. Tis is indeed true, since for a convex entropy U U i =θu i + θu i U convex θu i + θu i = θu i + θu i Jensen θu i + θu i = U i, 7 were te monotonicity of te mean value, te convexity of U, and Jensen s inequality ave been used. Tis proves Lemma Te scaling limiter 7 is entropy stable for θ [0, ]. 7. Entropy conservative fluxes for constant bottom topograpy b Here, te one-parameter entropy conservative numerical flux f a = a {v }+ +a { } {v } 9 will be considered. Inserting tis in a FV evolution equation 7 and dropping te bars, + i = i t f a x u i, u i+ f a u i, u i 7 = i t a i v i + i+ v i+ + + a i + i+ v i + v i+ x = i [ t + i t x a i v i + i v i + a i + i ] [ + a t a v i+ v i i+ x x [ a v i + + a ] v i. v i + v i v i+ + + a ] v i Considering non-negative water eigt, if i = 0, i, i+ > 0, v i = 0, a v i < 0, a v i+ > 0, te eigt + i becomes negative for t > 0. If only positive water eigt sould be considered, using te same conditions as before but i > 0, te new water eigt can be guaranteed to be positive, but only under a CFL condition depending on i wit allowed t 0 as i 0.

33 Entropy stable, well-balanced, positivity preserving metods for sallow water Lemma Te one-parameter family of entropy conservative fluxes 9 for te sallow water equations wit constant bottom topograpy b is not positivity preserving under a CFL condition not blowing up as 0. If te full two-parameter family 6 of entropy conservative fluxes is considered, terms of te form vi± v i± v i v i± vi + v i ave to be added. Since tese terms may be of arbitrary size and sign and do not contain any multiple of te water eigt, tey can render te water eigt negative. Tus, tis family is not positivity preserving, too. 7. Adding dissipation for constant bottom topograpy b Classically, adding dissipation to some kind of central flux f cent would result in f cent P [u], were P is positive semi-definite. However, tis migt not result in an entropy stable sceme. Terefore, te amount of dissipation is better cosen as proportional to te jump of entropy variables [w] instead of [u], as used inter alia by Bart 999; Fjordolm et al 0; Gassner et al 06b. In te simplest case, P = λ I wit λ 0. Ten λ[u] λ u w [w]. Since w u = u w, w u = U u, and U is convex, tis dissipation matrix is positive semi-definite. Tus, te resulting numerical flux f num = f a λ wu[w] 75 is entropy stable, were f a is te one-parameter entropy conservative flux 9 and w u is a suitable positive semi-definite approximation of te entropy Jacobian w u, e.g. te Jacobian evaluated at some mean value u. For te sallow water equations cf. 7, g + u w = v v v g g v, w u = v g v g Tus, adding dissipation in te form 75, te following additional contribution as to be added to te rigt and side of 7 for te entropy variables T w = g v, v, i.e. if te bottom topograpy b is continuous across cell boundaries t [ ] λ x i+ wu i+ w i+ w i λ i wu i w i w i 77 = t λ i+ i+ v i+ x g i + v i g + v i+ g v i+ v i + t λ i i + v i x g + i v i g v i g v i v i. Te additional positive values of i± are crucial to obtain te positivity preserving property under a suitable CFL condition, if λ is large enoug. Te

34 Hendrik Ranoca negative values of i are weigted by t and can tus be bounded by te positive terms in 7 if t is small enoug. Te remaining terms containing only te velocity but not te eigt of te water may be problematic. However, if te Jacobian u w is evaluated at te aritmetic mean value, v i+ = {v } i,i+ = v i+ + v i v i+ v i+ v i = v i+ v i, 7 and similarly for v i. Tus, tese additional terms vanis and cf. equation 7 [ + i = i t + a x + i+ t x + i t x [ λi+ [ λi v i+ v i λ i+ + λ i a + a ] v i+ + a v i ] v i + + a v i. ] t x 79 Hence, for λ i± max { v i, v i± }, te coefficients of i± are non-negative, and under te CFL condition t x + a v i+ + v i + λ i+ + λ i, 0 te new eigt + i is non-negative. Note tat tis CFL condition does not blow up as i 0. Wit tis estimate, te coice a = seems to be optimal in order to get te least restrictive CFL condition. Again, considering te full two-parameter family 6 instead of te one-parameter family 0 results in additional terms of order v not containing any contribution of te water eigt. Tus, tese terms are of arbitrary size and sign and destroy te positivity preservation as in Lemma Te one-parameter local Lax-Friedrics type flux 75 is entropy stable, if λ u w is positive semi-definite. Additionally, it is positivity preserving under te CFL condition 0, if u w = u w { }, {v }, λ max { v, v + }, and te bottom topograpy b is continuous across cell boundaries. { In te implementation, λ = max v + g, v + + } g + is cosen, as in te classical local Lax-Friedrics flux. However, if te bottom topograpy b is discontinuous across cell boundaries, additional terms ave to be considered, since te entropy variables are T w = g + b v, v. Tese terms are t x λ i+ b i+ b i + t x λ i b i + b i.

35 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 Adding tese terms to te rigt and side of equation 79, te CFL condition 0 gets lost. If all eigts are positive, it is possible to guarantee + i 0, but te corresponding CFL condition blows up as 0, since b i± b i may be of arbitrary size and as to be balanced by positive contributions of te terms. Wit te coice of u w as in Lemma, te dissipation term becomes λ uw {u } [w] {v } g g g [ + b] [v ] [v] = λ {v } g { } + = λ {v }[ + b] = λ [ + b] [v] + {v }[b] {v } g [ + b] {v } g [v] + {v } g {v }, g [v] + { }[v] + [v] {v } g [v] were te product rule 6 as been used. Tus, if b is continuous across cell boundaries, [b] = 0 and te dissipation term is simply te classical local Lax-Friedrics dissipation term λ [u]. Remark Te same numerical flux can also be obtained as a Rusanov type flux. Coosing a dissipation approximately as f u [u] u f w u[w] and using te scaling of Bart 999, Teorem for te eigenvectors, resulting in u f = R Λ R and w u = RR T, were Λ is te diagonal matrix of eigenvalues of f u, yields f u [u] R Λ R T [w]. Tus, te Rusanov coice of dissipation wit Λ = λ I, were λ > 0 is te largest eigenvalue, yields exactly te same Lax Friedrics type dissipative flux 75. Remark 5 A Roe type dissipation operator can be constructed by coosing Λ = diag λ i, were λi are te eigenvalues of f u. However, some attempts to prove te preservation of positivity ave not been successful. 7. Classical numerical fluxes for constant bottom topograpy b Not only te local Lax Friedrics type numerical flux 75 in te previous section is positivity preserving and entropy stable, but also its classical variant f num = {f } λ [u]. Tis can be establised using general results of Frid 00, 00; Boucut 00, but also by direct calculation. Using te FV update procedure 7, te water eigt after one time step using te local Lax Friedrics flux becomes + i = i t f num u i, u i+ f num u i, u i 5 x

36 6 Hendrik Ranoca = i t i+ v i+ + i v i λ i+ x i+ i = i [ t x iv i + i v i λ i + λ ] i+ + λ i i i λ i+ + v i+ λ i i+ + i v i. Tus, if λ i± > v x i± and t < λ i +λ, te new water eigt + i+ i is nonnegative, if te previous water eigts i, i, i+ are non-negative. Te entropy stability in te semidiscrete setting wit vanising bottom topograpy b can be establised by [w] f num [ψ] 6 = g [] [v ] {v } λ [] + [v] {v } + g { } λ [v] g [ v] =g {v }[] {v } {v }[v] λ g [] + λ {v }[][v] + {v }[v] + g { }[v] λ { }[v] λ {v }[][v] g [ v] = λ g [] λ { }[v] + g {v }[] {v } {v }[v] + {v }[v] g { } {v }[], were te product rule 6 as been used. Direct calculation yields [w] f num [ψ] 7 = g + λ v + v }{{} + λ v + v + v }{{}}{{}}{{} 0 }{{} λ v v + v }{{}}{{}}{{} , if +, 0 and λ v +, v. Tis proves Lemma 6 Te local Lax-Friedrics flux used in a simple FV metod for te sallow water equations wit constant bottom topograpy b is positivity preserving, if te water eigt is non-negative and te CFL condition t λ i + λ i+ x 0

37 Entropy stable, well-balanced, positivity preserving metods for sallow water 7 is fulfilled. Additionally, it is entropy stable in te semidiscrete sense, if te water eigts are non-negative +, 0 and λ v +, v. { In te implementation, λ = max v + g, v + + } g + is cosen. Remark 7 Using a Godunov sceme wit exact solution of te Riemann problem results in an entropy stable and positivity preserving sceme, since te exact solution as tese properties. However, some nonlinear root finding algoritm as to be applied to compute te exact solution of te Riemann problem. Terefore, tis will not be considered ere in more detail. Te Riemann problem for te sallow water equations wit vanising bottom topograpy is described in detail inter alia by Holden and Risebro 00, Capter 5. Remark Harten et al 9 proposed an approximate Riemann solver using only one intermediate state, known as HLL Riemann solver, using estimates of te slowest and fastest wave speeds s, s +. If tese estimates are lower and upper bounds of te wave speeds in te solution of te Riemann problem, tis flux is entropy stable, as remarked by Harten et al 9. Additionally, by te rigt coice of wave speed estimates as in te famous HLLE version for Euler equations proposed by Einfeldt 9, te numerical flux is positivity preserving, similar to te results for gas dynamics establised by Einfeldt et al 99, since te intermediate state in te approximating solution of te Riemann problem satisfies tis condition. 7.5 Suliciu relaxation solver for constant bottom topograpy b Te Suliciu relaxation solver described by Boucut 00, Section. for te sallow water equations as been implemented wit some tecnical modifications to allow vanising water eigt. Tis numerical flux is entropy stable and positivity preserving under a CFL condition not blowing up as Kinetic solver for constant bottom topograpy b Using a kinetic approac, Pertame and Simeoni 00 proposed a finite volume metod for te sallow water equations wit general bottom topograpy b tat as te tree desired properties, i.e. tat is entropy stable, positivity preserving, and well-balanced under a suitable CFL condition not blowing up as 0. However, te corresponding numerical flux as to be evaluated by some quadrature if te bottom topograpy varies. Tis will not be used ere. However, in te case of a constant topograpy, all integrals can be evaluated analytically.

38 Hendrik Ranoca 7.7 Hydrostatic reconstruction approac for general bottom topograpy b Te ydrostatic reconstruction as been introduced by Audusse et al 00 as a means to extend a numerical flux for te sallow water equations wit constant bottom topograpy b to varying b, preserving good properties of te numerical flux, notably entropy stability and positivity preservation. In addition, te resulting metod is well-balanced for te lake-at-rest initial condition. In te context of extended numerical fluxes incorporating te source term, te flux using ydrostatic reconstruction can be described as follows: Compute at first te limited values i = max { 0, i + b i max {b i, b k } }, k = max { 0, k + b k max {b i, b k } } and use te extended numerical fluxes wit new arguments f num i, v i, k, v k, f num v i, v i, k, v k + g i i, 9 for te water eigt and te discarge v to compute te rate of cange in cell i influenced by cell k. Tis results in a consistent and well-balanced numerical flux tat is positivity preserving and entropy stable, if te given fluxes f num, fv num ave tese properties for te sallow water equations wit constant bottom b. However, tis ydrostatic reconstruction as some disadvantages for some combinations of bottom slope, mes size, and water eigt, as described by Delestre et al 0, at least if used for a first order FV sceme. Remark 9 By an approac based on relaxation, Berton and Calons 06 constructed an approximate Riemann solver tat as te tree desired properties, i.e. tat is entropy stable, positivity preserving, and well-balanced under a suitable CFL condition. However, te existence of some parameters and a suitable time step are based on asymptotic arguments and can terefore not be implemented directly. However, te sallow water equations are derived based on te assumption of low variations in te bottom topograpy. Hence, a discretisation of it tat is continuous across elements seems to be natural. Finite volume subcells Altoug te analysis of te previous sections suggests tat te semidiscretisation of Teorem 6 wit appropriate positivity preserving and entropy stable fluxes of section 7 and te positivity preserving limiter of Zang and Su 0, described also in section 7., is stable, tere are problems at wetdry fronts in te practical implementation. Tese problems can be andled by some appropriate limiting strategy, e.g. TVB limiters used by Xing et al 00 or te slope limiter used by Duran and Marce 0. However, since te ig order of te approximation is lost in tese cases, te approac of finite volume subcells used similarly by Meister and Ortleb 06 in te context of te sallow water equations will be pursued. Applications of finite volume

39 Entropy stable, well-balanced, positivity preserving metods for sallow water 9 subcells ave also been proposed inter alia by Huerta et al 0; Dumbser et al 0; Sonntag and Munz 0. Additionally, given te interpretation of SBP metods wit diagonal norm as subcell flux differencing metods by Fiser et al 0; Fiser and Carpenter 0, te usage of FV subcells seems to be quite natural. In order to use finite volume subcells to compute te time derivative, te general procedure can be described as follows:. Decide, weter te ig-order discretisation or FV subcells of first order sould be used.. Project te polynomial of degree p onto a piecewise constant solution.. Compute te classical FV time derivative. As a detector to use FV subcells, te water eigt in te element or adjacent elements will be used, as described in section 9. Te projection in step is done for a diagonal-norm nodal SBP basis simply by taking subcells of lengt M i,i wit corresponding value u i. Tis is not an exact projection for te polynomial u in general, but is very simple and fits to te subcell flux differencing framework of Fiser et al 0; Fiser and Carpenter 0. It as also been used by Sonntag and Munz 0 in te context of te Euler equations. Tus, for an SBP SAT semidiscretisation t u = VOL + SURF M R T B f num, 90 te surface terms SURF are set to zero and te numerical flux f num is computed using te outer values u 0, u p instead of a iger order interpolation R u for nodal bases including boundary points, tis makes no difference. Te volume terms VOL are computed via FV subcells as [VOL] 0 = f 0, num, M 0,0 [VOL] p = f num p,p M p,p, 9 [VOL] i = f num i,i+ f num i,i M i,i for i {,..., p }. 9 Terefore, te numerical flux terms M R T B f num and te volume terms VOL form togeter a finite volume discretisation of te subcells. 9 Numerical tests In tis section, some numerical tests of te proposed scemes will be performed. For te time integration, te tird-order, tree stage SSP metod SSPRK, given by Gottlieb and Su 99 will be used, see also Gottlieb et al 0.

40 0 Hendrik Ranoca 9. Well-balancedness and entropy conservation In tis test case, te well-balancedness and entropy conservation properties of te two-parameter family of fluxes 6 and corresponding volume 5 and surface terms 6, were te parameters are cosen according to 6 and te free parameters m, k 9, k 0, k, l 0 ave been set to zero, are investigated. Te domain [, ] is equipped wit periodic boundary conditions and te solution is evolved in te time interval [0, ] using 000 steps of te SSPRK, metod. Te bottom topograpy and initial condition are given by bx = sin πx, 0x = bx, v 0 x = 0 9 for te lake-at-rest test case. Using N = 0 p+ = 5 elements wit polynomials of degree p = 7, represented using Gauß nodes, te simulations ave been { } performed for a, a + k 0, 0 k 60 wit gravitational constant g =. Te results, visualised in Figure sow{ te excellent well-balancedness of te metods. Te maximum errors max 0, } v v0 computed at te nodes are of order 0 using Float6 i.e. double precision in Julia 0..7 Bezanson et al 0. Coosing te initial condition 0 x =, v 0 x = 0, 9 +b is no longer constant, but te solution remains smoot at first. Computing again until t = wit te same parameters as before, te loss of entropy is visualised in Figure. Using 000 steps of SSPRK,, te relative entropy dissipation U U0 / U0 is of order 0 0, wit variations of order 0 for different parameters a, a. Tis loss of entropy is caused by te time integrator, as can be seen by refining te time step. Using 000 steps, te relative dissipation is order order 0 wit variations of order 0. Te smoot solutions for a =, a = a = are plotted in Figure. Te influence of te numerical surface flux is visualised in Figure. Tere, te number of degrees of freedom N p + as been kept constant, wile te polynomial degree p varies between 0 first order FV sceme and 5. As can be seen tere, te entropy conservative flux f, is indeed entropy conservative, wile te entropy stable fluxes are a bit dissipative. Te dissipation increases from te Suliciu flux Boucut 00, Section. over te kinetic flux Pertame and Simeoni 00 to te local Lax-Friedrics flux. All tree fluxes ave been implemented using te ydrostatic reconstruction of Audusse et al 00 described in section 7.7. In te finite volume setting p = 0, te dissipation is of order 0 and decreases wit increasing polynomial degree p. For p, te curves for te dissipative fluxes become visually indistinguisable and for p 5 tey coincide wit te entropy conservative flux f, for tis smoot solution.

41 Entropy stable, well-balanced, positivity preserving metods for sallow water a = a e.0.7 a 0 0 a { Fig. : Maximum norm error max 0, } v v 0 of solutions computed using te entropy conservative fluxes wit varying parameters a, a for te lake-at-rest initial condition 9. a 0 a = a 0 a e.55e a 0 a = a 0 a e 0..e a 000 steps. b 000 steps. Fig. : Relative entropy dissipation U U0 / U0 of solutions computed using te entropy conservative fluxes wit varying parameters a, a for te initial condition 9.

42 Hendrik Ranoca +b b +b b x a Water eigt. v v 0 v x b Discarge v. Fig. : Solutions computed using te entropy conservative fluxes wit parameters a =, a = a = for te initial condition Lake at rest wit emerged bump Here, te lake-at-rest initial condition of SWASHES Delestre et al 0, Section.. bx = { x 0, if < x <, 0, else, 0 x = max { 0., bx } bx, v 0 x = 0, 95 will be used in te domain [0, 5] wit periodic boundary conditions for simulations in te time interval [0, ] wit gravitational constant g = 9.. If no FV subcells are used, te result strongly depends on te resolution of te sore and needs in general some additional dissipation to be stable near te wet-dry front. However, activating FV subcells if te water eigt at some node in te element is smaller tan 0 5, te simulation is stable. Tese results are sown in Figure 5 for N = 0 elements of polynomials of degree p = 5 and te local Lax-Friedrics flux wit ydrostatic reconstruction as numerical flux. { Te maximum error norm max 0, } v v0 computed at te nodes is of order of magnitude 0 6 for varying parameters { } a, a + k 0, 0 k 60 used for te volume terms 5. Again, Gauß nodes and corresponding surface terms ave been used, were te additional free parameters ave been set to zero. Additionally, te water eigt for te coice a =, a = a = is visualised tere.

43 Entropy stable, well-balanced, positivity preserving metods for sallow water Ut U0 Ut U0 Ut U0 / U0 / U0 / U f, Suliciu kinetic LLF t a p = 0, N = 0 p+ = e f, Suliciu kinetic LLF t c p =, N = 0 p+ = e f, Suliciu kinetic LLF t e p =, N = 0 p+ = 0. Ut U0 Ut U0 Ut U0 / U0 / U0 / U e f, Suliciu kinetic LLF t b p =, N = 0 p+ = 60. f, Suliciu kinetic LLF t d p =, N = 0 p+ = e 0 f, Suliciu kinetic LLF t f p = 5, N = 0 p+ = 5. Fig. : Relative entropy dissipation U U0 / U0 of solutions computed using different surface fluxes for te initial condition 9 wit varying degree p and number of elements N = 0 p+.

44 Hendrik Ranoca a = a e b +b b a 0. +b a x a Maximum error. b Water eigt for a =, a =. { Fig. 5: Maximum norm error max 0, } v v0 of solutions computed using te entropy conservative fluxes wit varying parameters a, a for te lake-at-rest wit emerged bump initial condition 95 and solution at t = for a =, a =. 9. Moving water equilibrium wit varying bottom b Here, a moving water equilibrium of te sallow water equations wit gravitational constant g = 9. given by v m = const, is considered. Te bottom topograpy is v + g + b E = const 96 bx = { cos 0πx + +, if 0. < x < 0., 0, else, 97 and te initial condition is computed by solving te second equation of 96 for, inserting v = v m =, E = mg/ + g Audusse et al 05. = m. Two initial conditions m =, E = 5 and = are considered, similar to { Computing te maximum error max 0, } v v0 at te nodes yields identical results for bot initial conditions wit polynomial } degrees p {0,..., 9}, and parameters a, a { + k0, 0 k 60 for te volume terms, wile te local Lax-Friedrics flux as been used as numerical flux. Te domain is divided into N = 0 elements using Gauß nodes. Tese results wit m =, E = 5 are sown in Figure 6. As can be seen tere, te coice a = a is optimal for tis problem, wile te coice of a does not seem to be critical. Tis can be explained by te additional term in v for a a in te numerical flux 6 and te corresponding volume terms 5. Te results for m =, E = mg/ are visually indistinguisable.

45 Entropy stable, well-balanced, positivity preserving metods for sallow water 5 a 0 a = a 0 a min a,a u u p a Maximum error for p = 9. b Minimal errors for varying p. { Fig. 6: Maximum norm error max 0, } v v0 of solutions computed using te entropy conservative fluxes wit varying parameters a, a for moving water equilibrium 96 wit m =, E = 5 for polynomials of degree p = 9 and minimal values of te maximum errors over a, a. Additionally, te minimal values of te maximum error over te parameters a, a are plotted in Figure 6 for m =, E = 5. Te usual superior properties of odd polynomial degrees p as well as exponential convergence can be seen tere. 9. Dam break Here, te dam break problem wit dry domain and an analytical solution described by Delestre et al 0, Section.. will be considered. Te initial condition { 0.005, if x < 5, 0 x = v 0 x = 0, bx = 0, 9 0, else, is evolved in te domain [0, 0] until t = 6, and te gravitational constant is again g = 9.. Te results of a simulation using N = 00 elements wit polynomials of degree p = and te local Lax-Friedrics numerical flux are plotted in Figure 7. Here, FV subcells are used in a cell if te water eigt in te cell itself or adjacent cells is less ten 0 6, and te parameters are cosen as a =, a = a =. Motivated by te result of section 9., only te parameter a as been varied for tis problem, wile te parameter } a is fixed at a = a. te results for a { + k0, 0 k 60 are sown in Figures and 9. Tere, te L errors ave been computed exactly for te polynomials using Gauß nodes and te errors are computed at te same nodes. In tese experiments, Gauß nodes yield a lower error in te solutions, bot in L and and tis error is nearly independent of te parameter a