arxiv: v1 [math.na] 14 Jan 2019
|
|
- Ambrose Carroll
- 5 years ago
- Views:
Transcription
1 arxiv:9.423v [math.na] 4 Jan 29 Galerkin finite element methods for the Shallow Water equations over variable bottom G. Kounadis,2 and V.A. Dougalis,2 Department of Mathematics, National and Kapodistrian University of Athens, 5784 Zografou, Greece 2 Institude of Appled and Computational Mathematics, FORTH, 73 Heraklion, Greece Abstract We consider the one-dimensional shallow water equations (SW) in a finite channel with variable bottom topography. We pose several initialboundary-value problems for the SW system, including problems with transparent (characteristic) boundary conditions in the supercritical and the subcritical case. We discretize these problems in the spatial variable by standard Galerkin-finite element methods and prove L 2 -error estimates for the resulting semidiscrete approximations. We couple the schemes with the 4 th order-accurate, explicit, classical Runge-Kutta time stepping procedure and use the resulting fully discrete methods in numerical experiments of shallow water wave propagation over variable bottom topographies with several kinds of boundary conditions. We discuss issues related to the attainment of a steady state of the simulated flows, including the good balance of the schemes. Keywords: Shallow water equations, Standard Galerkin finite element method, error estimates, characteristic boundary conditions, variable bottom topography 2 MSC: 65M6, 65M2 Introduction In this paper we will consider standard Galerkin finite element approximations to the one-dimensional system of shallow water equations over a variable bottom that we write following [], as η t + (ηu) x + (βu) x =, u t + η x + uu x =. (SW) The system (SW) approximates the two-dimensional Euler equations of water wave theory and models two-way propagation of long waves of finite amplitude on the surface of an ideal fluid in a channel with a variable bottom. The variables in (SW) are nondimensional and unscaled; x R and t are proportional
2 to position along the channel and time, respectively. With the depth variable z taken to be positive upwards, the function η = η(x, t) is proportional to the elevation of the free surface from a level of rest corresponding to z = and u = u(x, t) is proportional to the horizontal velocity of the fluid at the free surface. The bottom of the channel is defined by the function z = β(x); it will be assumed that β(x) >, x R, and that the water depth η(x, t)+β(x) is positive for all x, t. It should be noted that there are several equivalent formulations of the system represented by (SW), some of which will be considered in section 3 of the paper. It is well known that given smooth initial conditions η(x, ) = η (x), u(x, ) = u (x), x R, and smooth bottom topography, the Cauchy problem for (SW) has smooth solutions, in general only locally in t. In this paper we will be concerned with numerical approximations of (SW) and suppose that its solution is sufficiently smooth so that the error estimates of section 2 hold. We will specifically consider three initial-boundary-value problems (ibvp s) for (SW), posed on the spatial interval [, ]: A simple ibvp with vanishing fluid velocity at the endpoints and two ibvp s with transparent (characteristic) boundary conditions, in the supercitical and subcritical flow cases, respectively. For these types of ibvp s there exists a well-posedness theory locally in t, cf. e.g. [2], [3], [4]. For the formulation and numerical solution of ibvp s with transparent boundary conditions see also [5], [6]. In section 2 we will specify in detail these ibvp s and summarize their well-posedness theory. The literature on the numerical solution of the shallow water equations is vast. We will just mention that in recent years there has been considerable interest in solving them numerically by Discontinuous Galerkin finite element methods and refer the reader to [7] and the recent surveys [8], [9], for an overview of issues related to the implementation of such methods in the presence of discontinuities and also in two space dimensions. In section 2 of the paper we consider the ibvp s previously mentioned, discretize them in space by the standard Galerkin finite element method, and prove L 2 -error estimates for the semidiscrete approximations assuming smooth solutions of the equations and extending results of [], [], to the variable bottom case. In section 3 we discretize the semidiscrete problem in the temporal variable using he classical fourth-order accurate, four-stage explicit Runge-Kutta method. The resulting fully discrete scheme is stable under a Courant number stability condition and its convergence has been analyzed for (SW) in the case of a horizontal bottom in [2]. We use this scheme in a series of numerical experiments simulating shallow water wave propagation over variable bottom topography and in the presence of absorbing (characteristic) boundary conditions up to the attainment of steady-state solutions. We also discuss issues of good balance, cf. [3], [7], of the standard Galerkin method applied to the shallow water equations written in balance-law form. In the sequel we denote, for integer m, by H m = H m (, ) the usual L 2 -based real Sobolev spaces of order m, and by m their norm. The space H = H (, ) will consist of the H functions that vanish at x =,. The inner product and norm on L 2 = L 2 (, ) will be denoted by (, ),, respectively, 2
3 while C m will be the m times continuously differentiable functions on [, ] The norms of L and of the L -based Sobolev space W, on (, ) will be denoted by,,, respectively. P r will be the space of polynomials of degree at most r. 2 Initial-boundary-value problems and error estimates In this section we will specify the initial-boundary-value problems (ibvp s) for the shallow water equations to be analyzed numerically, their Galerkin-finite element space discretizations and the properties of the attendant finite element spaces. We will then prove L 2 -error estimates for these discretizations assuming that the data and the solutions of the ibvp s are smooth enough for the purposes of the error estimation. 2. Semidiscretization of a simple ibvp with vanishing fluid velocity at the endpoints We consider first a simple ibvp for (SW) posed in the finite channel [, ]. let T > be given. We seek η = η(x, t), u = u(x, t), for x, t T, satisfying η t + (ηu) x + (βu) x =, u t + η x + uu x =, η(x, ) = η (x), u(x, ) = u (x), x u(, t) = u(, t) =, t T. x, t T, (2.) In [2] Petcu and Temam, using an equivalent form of (2.), established the existence-uniqueness of solutions (η, u) of (2.) in H 2 H 2 H for some T = T ( η 2, u 2 ) under the hypotheses that η H 2, and, say, β H 2, such that η (x) + β(x) >, x [, ], and u H 2 H. Moreover, it holds that η(x, t)+β(x) > for (x, t) [, ] [, T ], i.e. the water depth is always positive. (This property will be assumed in all the error estimates to follow in addition to the sufficient smoothness of η and u.) In order to solve (2.) numerically let = x < x 2 <... < x N+ = be a quasiuniform partition of [, ] with h := max i (x i+ x i ), and for integers k, r such that r 2, k r 2, consider the finite element spaces S h = {ϕ C k : ϕ [xj,x j+] P r, j N} and S h, = {ϕ S h : ϕ() = ϕ() = }. It is well known that given w H r, there exists χ S h such that w χ + h w χ Ch r w (r), (2.2a) and, in addition, if r 3, such that w χ 2 Ch r 2 w (r), (2.2b) 3
4 where C is a constant independent of h and w; a similar property holds in S h, provided w H r H. It follows from (2.2a), cf. [4], that if P is the L 2 -projection operator onto S h, then P w C w, w H, (2.3a) P w C w, w C, (2.3b) P w w L Ch r w (r), w C r, (2.3c) and that the analogous properties also hold for P, the L 2 -projection operator onto S h,. In addition, as a consequence of the quasiuniformity of the mesh, the inverse properties χ Ch χ, χ j, Ch (j+/2) χ, j =,, (2.4) hold for χ S h or χ S h,. The standard Galerkin semidiscretization of (2.) is defined as follows: Seek η h : [, T ] S h, u h : [, T ] S h,, such that for t [, T ] (η ht, ϕ) + ((η h u h ) x, ϕ) + ((βu h ) x, ϕ) =, ϕ S h, (u ht, χ) + (η hx, χ) + (u h u hx, χ) =, χ S h,, (2.5) with initial conditions η h () = P η, u h () = P u. (2.6) We will prove below that the semidiscrete approximations (η h, u h ) satisfy an L 2 -error bound of O(h r ). Is is well known that this order of accuracy cannot be improved in the case of the standard Galerkin finite element method for firstorder hyperbolic problems in the presence of general nonuniform meshes, [5], []; for uniform meshes better results are possible, cf. [] and the numerical experiments of section 3. Proposition 2.. Let (η, u) be the solution of (2.), assumed to be sufficiently smooth and satisfying β + η > in [, ] [, T ], where β C, β >. Let r 3 and h be sufficiently small. Then, the semidiscrete ivp (2.5) (2.6) has a unique solution (η h, u h ) for t [, T ], such that max ( η η h + u u h ) Ch r, (2.7) t T where, here and in the sequel, C will denote a generic constant independent of h. Proof. As the proof is similar to that of Proposition 2.2 in [], which is valid in the case of horizontal bottom (β(x) = ), we will only indicate the steps where the two proofs differ. We let ρ := η P η, θ := P η η h, σ := u P u, ξ := P u u h. While the solution exists we have (θ t, φ) + (β(ξ x + σ x ), φ) + (β x (ξ + σ), φ) + ((ηu) x (η h u h ) x, φ) =, φ S h, (2.8) (ξ t, χ) + (θ x + ρ x, χ) + (uu x u h u hx, χ) =, χ S h,. (2.9) 4
5 Taking φ = θ in (2.8) and integrating by parts we have d 2 dt θ 2 + ([(β + η)ξ] x, θ) = (βσ x, θ) (β x σ, θ) ((ησ) x, θ) ((uρ) x, θ) ((uθ) x, θ) + ((ρσ) x, θ) + ((θσ) x, θ) + ((ρξ) x, θ) + ((θξ) x, θ). (2.) In view of (2.6), we conclude by continuity that there exists a maximal temporal instance t h > such that (η h, u h ) exist and ξ x for t t h. Suppose that t h < T. Using the approximation and inverse properties of S h and S h,, we may then estimate the various terms in the r.h.s. of (2.) for t [, t h ] in a similar way as in [], since β C, and conclude that for t [, t h ] d 2 dt θ 2 (γ, θ x ) C(h r θ + θ 2 + ξ 2 ), (2.) where we have put γ := (β + η)ξ. We turn now to (2.9) in which we take χ = P γ = P [(β + η)ξ]. t t h it follows that (ξ t, γ) + (θ x, P γ) = (ρ x, P γ) ((uξ) x, P γ) ((uσ) x, P γ) Arguing now as in [], since β C, noting that For + ((σξ) x, P γ) + (σσ x, P γ) + (ξξ x, P γ). (2.2) ((uξ) x, P γ) = ((uξ) x, P γ γ) + (u x (β + η), ξ 2 ) 2 ([β + η)u x], ξ 2 ), and using a well-known superapproximation property of S h, to estimate the term P γ γ: we get P γ γ = P [(β + η)ξ] (β + η)ξ Ch ξ, ((uξ) x, P γ) Ch ξ ξ + C ξ 2 C ξ 2. With similar estimates as in [], using the hypothesis that ξ x for t t h, we conclude from this inequality and (2.2) that for t t h (ξ t, (β + η)ξ) + (θ x, P γ) C(h r ξ + ξ 2 ). (2.3) Adding now (2.2) and (2.3) we obtain d 2 dt θ 2 + (ξ t, (β + η)ξ) + (θ x, P γ γ) C[h r ( θ + ξ ) + θ 2 + ξ 2 ]. But, since β = β(x), we have (ξ t, (β + η)ξ) = 2 dt ((β + η)ξ, ξ) 2 (η tξ, ξ). Therefore, for t t h d 2 dt [ θ 2 + ((β + η)ξ, ξ)] C[h r ( θ + ξ ) + θ 2 + ξ 2 ], for a constant C independent of h and t h. Since β + η >, the norm ((β + η), ) /2 is equivalent to that of L 2 uniformly for t [, T ]. Hence, Gronwall s inequality and (2.6) yield for a constant C = C(T ) θ + ξ Ch r for t t h. (2.4) d 5
6 We conclude from (2.4), using inverse properties, that ξ x Ch r 5/2 for t t h, and, since r 3, if h is taken sufficiently small, we see that t h is not maximal. Hence we may take t h = T and (2.7) follows from (2.4). The hypothesis that r 3 seems to be technical, as numerical experiments indicate that (2.7) apparently holds for r = 2 as well, cf. []. 2.2 Semidiscretization of an ibvp with absorbing (characteristic) boundary conditions in the supercritical case We consider now the shallow water equations with variable bottom with transparent (characteristic) boundary conditions. First we examine the supercritical case: For (x, t) [, ] [, T ] we seek η = η(x, t) and u = u(x, t) satisfying the ibvp η t + (βu) x + (ηu) x =, u t + η x + uu x =, x, t T, (2.5) η(x, ) = η (x), u(x, ) = u (x), x, η(, t) = η, u(, t) = u, t T, where β C, η, u are given functions on [, ] and η, u constants such that β(x) + η >, u >, u > β(x) + η, x [, ]. The ibvp (2.5) was studied by Huag et al., [3], in the more general case of the presence of a lateral component of the horizontal velocity depending on x only (nonzero Coriolis parameter). In the simpler case of (2.5), we assume that (η, u ) is a suitable constant solution of (2.5) and that η (x), u (x) are sufficiently smooth initial conditions close to (η, u ) and satisfying appropriate compatibility relations at x =. Then, as is proved in [3], given positive constants c, α, ζ, and ζ, there exists a T > and a sufficiently smooth solution (η, u) of (2.5) satisfying for (x, t) [, ] [, T ] the strong supercriticality properties u 2 (β + η) c 2, u α, ζ (β + η) ζ. (2.6a) (2.6b) (2.6c) For the purposes of the error estimation to follow we will assume in addition that the solution of (2.5) satisfies a strengthened supercriticality condition of the following form: There exist positive constants a, and b, such that for (x, t) [, ] [, T ] β + η b, u 2a, β + η (u a)(u 2a 3 ). (2.7a) (2.7b) (2.7c) 6
7 Obviously (2.7a) and (2.7b) imply that u β + η. It is not hard to see that (2.7c) follows from (2.6a) (2.6c) if e.g. α is taken sufficiently small and c sufficiently large. We also remark here that in the error estimates to follow (2.7c) will be needed only at x = for t [, T ]. We will approximate the solution of (2.5) in a slightly transformed form. We let η = η η, ũ = u u and rewrite (2.5) as an ibvp for η and ũ with homogeneous boundary conditions. Dropping the tildes we obtain the system η t + u η x + (β + η )u x + (ηu) x + (u + u )β x =, u t + η x + u u x + uu x =, η(x, ) = η (x) η, u(x, ) = u (x) u, x, η(, t) =, u(, t) =, t T. In terms of the new variables (2.7a) (2.7c) become x, t T, (2.8) β + η + η b, u + u 2a, β + η + η (u + u a)(u + u 2a 3 ). (2.9a) (2.9b) (2.9c) In the rest of this subsection, for integer k, let C k = {v C k [, ] : v() = }, and H k+ = {v H k+ (, ) : v() = }. Using the hypotheses of section 2. on the finite element space discretization we define Sh = {φ C r 2 : φ P [xj,x j+] r, j N} and P the L 2 projection operator onto Sh. Note that (2.2) (2.4) also hold on Sh mutatis mutandis. The standard Galerkin semidiscretization of (2.8) is defined as follows: We seek η h, u h, : [, T ] Sh such that for t T (η ht, φ) + (u η hx, φ) + ((β + η )u hx, φ) + ((η h u h ) x, φ) + ((u h + u )β x, φ) =, φ Sh, (2.2) with (u ht, φ) + (η hx, φ) + (u u hx, φ) + (u h u hx, φ) =, φ Sh, (2.2) η h () = P (η ( ) η ), u h () = P (u ( ) u ). (2.22) The boundary conditions implied by the choice of Sh are no longer exactly transparent, but they are highly absorbing as will be seen in the numerical experiments of Section 3. Proposition 2.2. Let (η, u) be the solution of (2.8), and assume that the hypotheses (2.9a) (2.9c) hold, that r 3, and h is sufficiently small. Then the semidiscrete ivp (2.2) (2.22) has a unique solution (η h, u h ) for t T satisfying max ( η(t) η h(t) + u(t) u h (t) ) Ch r. (2.23) t T 7
8 Proof. Let ρ = η P η, θ = P η η h, σ = u P u, ξ = P u u h. After choosing a basis for Sh, it is straightforward to see that the semidiscrete problem represents an ivp for an ode system which has a unique solution locally in time. While this solution exists, it follows from (2.2) (2.22) and the pde s in (2.8), that (θ t, φ) + (u (ρ x + θ x ), φ) + ((β + η )(σ x + ξ x ), φ) + ((ηu η h u h ) x, φ)+ ((σ + ξ)β x, φ) =, φ Sh, (ξ t, φ) + (ρ x + θ x, φ) + (u (σ x + ξ x ), φ) + (uu x u h u hx, φ) =, φ Sh Proceeding as in the proof of Proposition 2. of [], which is valid for a horizontal bottom, we obtain from the above in the case of variable bottom that (θ t, φ) + (u θ x, φ) + (γ x, φ) + ((uθ) x, φ) ((θξ) x, φ) = (R, φ), φ Sh, (2.24) (ξ t, φ) + (θ x, φ) + (u ξ x, φ) + ((uξ) x, φ) (ξξ x, φ) = (R 2, φ), φ Sh, (2.25) where γ = (β + η + η)ξ and R = u ρ x + (β + η )σ x + σβ x + (ησ) x + (uρ) x (ρσ) x (ρξ) x (θσ) x, (2.26) R 2 = ρ x + u σ x + (uσ) x (σξ) x σσ x. (2.27) Putting φ = θ in (2.24), using integration by parts, and suppressing the dependence on t we have d 2 dt θ 2 (γ, θ x ) + 2 (u + u())θ 2 () + (β() + η + η())ξ()θ() 2 ξ()θ2 () = 2 (u xθ, θ) + 2 (ξ xθ, θ) (R, θ) (2.28) Take now φ = P γ = P [(β + η + η)ξ] in (2.25) and get (ξ t, γ) + (θ x, γ) + (u ξ x, γ) + ((uξ) x, γ) (ξξ x, γ) = (R 3, P γ γ) (R 2, P γ), (2.29) where R 3 = θ x + u ξ x + (uξ) x ξξ x. (2.3) Integration by parts in various terms in (2.29) gives (ξ t, γ)+(θ x, γ)+ 2 (u +u())(β()+η +η())ξ 2 () 3 (β()+η +η())ξ 3 () where = (R 4, ξ) (R 3, P γ γ) (R 2, P γ), (2.3) R 4 = 2 u (β x + η x )ξ 2 u x(β + η + η)ξ + 2 u(β x + η x )ξ 3 (β x + η x )ξ 2. (2.32) 8
9 Adding now (2.28) and (2.3) we obtain d 2 dt where [ θ 2 + ((β + η + η)ξ, ξ) ] + ω = 2 (η tξ, ξ) 2 (u xθ, θ) + 2 (ξ xθ, θ) (R, θ) + (R 4, ξ) (R 3, P γ γ) (R 2, P γ), (2.33) ω = 2 (u + u())θ 2 () + 2 (u + u())(β() + η + η())ξ 2 () + (β() + η + η())ξ()θ() 2 ξ()θ2 () 3 (β() + η + η())ξ 3 (). (2.34) In view of (2.22), by continuity we conclude that there exists a maximal temporal instance t h > such that (η h, u h ) exist and ξ x a for t t h. Suppose that t h < T. Then, since ξ ξ x, it follows from (2.34) that for t [, t h ] ω 2 (u + u() a)θ 2 () + 2 (β() + η + η()) ( ) u + u() 2a 3 ξ 2 () ( ) ( ) µ λ θ() + (β() + η + η())ξ()θ() = 2 (θ(), ξ())t, (2.35) λ λν ξ() where µ = u + u() a, λ = β() + η + η(), ν = u + u() 2a 3. The hypotheses (2.9a) (2.9b) give that < µ < ν, λ >. It is easy to see then that the matrix in (2.35) will be positive semidefinite precisely when (2.9c) holds. Hence, (2.35) implies that ω. We now estimate the various terms in the right-hand side of (2.33) for t t h. As in the proof of Proposition 2. of [] adapted in the case of a variable β(x) C and using an appropriate variable-β superapproximation property to estimate P γ γ. We finally obtain from (2.33) and the fact that ω, that for t t h it holds that [ θ 2 + ((β + η + η)ξ, ξ) ] Ch r ( θ + ξ ) + C( θ 2 + ξ 2 ), d dt where C is a constant independent of h and t h. By (2.9a) the norm ((β + η + η), ) /2 is equivalent to that of L 2 uniformly for t [, T ]. Hence, Gronwall s inequality and the fact that θ() = ξ() = yield for a constant C = C(T ) θ + ξ Ch r for t t h. (2.36) We conclude from the inverse properties that ξ x Ch r 5/2 for t t h, and, since r 3, if h is taken sufficiently small, t h is not maximal. Hence we may take t h = T and (2.23) follows from (2.36). 2.3 Semidiscretization in the case of absorbing (characteristic) boundary conditions in the subcritical case We finally consider the shallow water equations with variable bottom in the presence of transparent (characteristic) boundary conditions in the subcritical case. In this case, instead of the variable η, we will use the total height of the 9
10 water, H = β + η. For (x, t) [, ] [, T ] we seek H = H(x, t) and u = u(x, t) satisfying the ibvp H t + (Hu) x =, u t + H x + uu x = β x, x, t T, (2.37) H(x, ) = H (x), u(x, ) = u (x), x, u(, t) + 2 H(, t) = u + 2 H, t T, u(, t) 2 H(, t) = u 2 H, t T, where H, u are given functions on [, ] and H, u constants such that H > and u 2 < H. Implicit in the formulation of the boundary conditions in (2.37) is that outside the spatial domain [, ] u and H are equal to constants u, H, respectively. The ibvp (2.37) in a slightly different but equivalent form was studied by Petcu and Temam, [4], under the hypotheses that for some constant c > it holds that u 2 H c 2 and that the initial conditions H (x) and u (x) are sufficiently smooth and satisfy the condition (u (x)) 2 H (x) c 2 and suitable compatibility relations at x = and x =. Under these assumptions one may infer from the theory of [4] that there exists a T > such that a sufficiently smooth solution (H, u) of (2.37) exists for (x, t) [, ] [, T ] with the properties that H is positive and the strong supercriticality condition u 2 H c 2, (2.38) holds for (x, t) [, ] [, T ]. Here we will assume that the solution satisfies a stronger subcriticality solution; specifically that for some constant c > it holds that u + H c, u H c, (2.39a) and for (x, t) [, ] [, T ] that u + H c, u H c. (2.39b) In this section we will approximate the solution of (2.37) after transforming the system in diagonal form. We write the system of pde s in (2.37) as ( ) ( ) ( ) Ht Hx + A = (2.4) u t u x β x ( ) u H where A =. The matrix A has eigenvalues λ u = u+ H, λ 2 = u H, (note that (2.39b) implies that λ c and λ 2 c in [, ] [, T ]), with associated eigenvectors X = ( H, ) T, X2 = ( H, ) T. If S is the matrix with columns X, X 2 it follows from (2.4) that ( ) ( ) ( ) ( ) S Ht λ + S Hx = S. (2.4) u t λ 2 u x β x
11 ( ) If we try to define now functions v, w on [, ] [, T ] by the equations S Ht = ( ) ( ) ( ) u t vt, S Hx vx =, we see that these equations are consistent and their w t u x w x solutions are given by v = 2 u + H + c v, w = 2 u H + c w, for arbitrary constants c v, c w. Choosing the constants c v, c w so that v(, t) =, w(, t) =, and using the boundary conditions in (2.37) we get v = 2[ u u + 2 ( )] [ H δ, w = 2 u u 2 ( )] H δ (2.42) where δ = H. The original variables H, u are given in terms of v and w by the formulas H = ( 2 (v w) + δ ) 2, u = v + w + u (2.43) Since λ = u + H = u + δ + 3v + w, λ 2 = u H = u δ + v + 3w 2 2 we see that the ibvp (2.37) becomes ( ) ( ) ( ) ( ) vt u + δ + + 3v+w 2 vx w t u δ + v+3w = 2 w 2 β x, x x, t T. v(x, ) = v (x), w(x, ) = w (x), x, v(, t) =, w(, t) =, t T, (2.44) (2.45) where v (x) = 2 [u (x) u + 2( H (x) δ )], w (x) = 2 [u (x) u 2( H (x) δ )]. Under our hypotheses (2.45) has a unique solution in [, ] [, T ] which will be assumed to be smooth enough for the purposes of the error estimation that follows. Given a quasiuniform partition of [, ] as in section 2., in addition to the spaces defined there, let for integer k C k = {f C k [, ] : f() = }, H k+ = {f H k+ (, ), f() = }, and, for integer r 2, S h = {φ C r 2 : φ [xj,x j+] P r, j N}. Note that the analogs of the approximation and inverse properties (2.2), (2.4) hold for Sh as well, and that the estimates in (2.3) are also valid for the L 2 projection P onto Sh, mutatis mutandis. The (standard) Galerkin semidiscretization of (2.45) is then defined as follows: Seek v h : [, T ] Sh, w h : [, T ] Sh, such that for t [, T ] (v ht, φ) + ((u + δ )v hx, φ) (v hv hx, φ) + 2 (w hv hx, φ) = 2 (β x, φ), φ Sh, (2.46) (w ht, χ) + ((u δ )w hx, χ) (w hw hx, χ) + 2 (v hw hx, χ) = 2 (β x, χ), χ Sh, (2.47) with v h () = P (v ), w h () = P (w ). (2.48)
12 The boundary conditions induced by the finite element spaces and the discrete variational formulation (2.46) (2.48) are no longer exactly transparent; they are highly absorbent nevertheless as will be checked in numerical experiments in Section 3. The main result of this section is Proposition 2.3. Let (v, w) be the solution of (2.45) and assume that the hypotheses (2.39a) (2.39b) hold, that r 3, and that h is sufficiently small. Then the semidiscrete ivp (2.46) (2.48) has a unique solution (v h, w h ) for t T that satisfies max ( v v h + w w h ) Ch r. (2.49) t T If (H, u) is the solution of (2.39) and we define then H h = [ 2 (v h w h ) + δ ] 2, u h = v h + w h + u, (2.5) max ( H H h + u u h ) Ch r. t T Proof. With the notation that we have introduced incorporating the variable β(x), it can be seen that the proof is entirely analogous to that of Proposition 3. of [] mutatis mutandis and will consequently be omitted. (Note that the source terms involving β x in the right-hand sides of (2.45), (2.46), (2.47) will cancel in the variational error equations.) 3 Numerical experiments In this section we present results of numerical experiments that we performed solving numerically the shallow water equations using standard Galerkin finite element space discretizations like the ones analyzed in the previous section. The semidiscrete schemes were discretized in the temporal variable by the classical, explicit, 4-stage, 4 th -order Runge-Kutta scheme (RK4), unless otherwise indicated. The resulting fully discrete scheme is stable and fourth-order accurate in time provided a Courant-number stability condition of the form k h α is imposed; here k denotes the (uniform) time step. In the case of a horizontal bottom the convergence of this scheme for the ibvp (2.) was analyzed in [2] and used in numerical experiments for the absorbing b.c. ibvp s (2.5) and (2.37) in []. In section 3. below we use this fully discrete scheme to study computationally various issues related to the discretization of the ibvp s with absorbing (characteristic) b.c. s considered in sections 2.2 and 2.3. In section 3.2 we write the shallow water equations in the form of a balance law and study various issues of the numerical solution of this model with Galerkin-finite element methods, including questions of good balance of the schemes. Since the numerical method 2
13 simulates only smooth solutions, initial conditions and bottom topographies were taken to be of small amplitude to ensure that no discontinuities developed within the time frame of the experiments. 3. Absorbing (characteristic) boundary conditions In the numerical experiments of this section we use the Galerkin finite element method with continuous, piecewise linear functions for the space discretization of the numerical solution of the ibvp s with absorbing (characteristic) boundary conditions considered in sections 2.2 and 2.3. The theoretical error estimates in Propositions 2.2 and 2.3 require at least piecewise quadratic elements, i.e. r 3, and predict L 2 -error bounds of O(h r ) for quasiuniform meshes. The results of numerical experiments shown in the sequel suggest that the method works with piecewise linear functions (i.e. r = 2) as well, and in this case the L 2 errors for a uniform mesh are of O(h 2 ). In the supercritical case, in order to find the numerical convergence rates of the scheme (2.2) (2.2) we consider an ibvp with η =, u = 3 and a bottom function and exact solution given for x (, ) by β(x) =.4 exp( (x.5) 2 ), η(x, t) = x exp( xt) + η, u(x, t) = ( x cos(πx)) exp(2t) + u. (3.) (The initial conditions and an appropriatee right-hand side were computed from these formulas.) The problem was solved with a uniform mesh with h = /N and k = h/. The L 2 errors and rates of convergence at T = are shown in Table. η u N L 2 error rate L 2 error rate 4.322e e e e e e e e e e Table : L 2 errors and rates of convergence at T =, r = 2, supercritical case, (3.), h = /N, k/h = / In the case of a subcritical flow we consider an ibvp with η =, u =, 3
14 and bottom function and exact solution given for x (, ) by b(x) =.4 exp( (x.5) 2 ), η(x, t) = (x + ) exp( xt), u(x, t) = (2x + cos(πx) ) exp(t) + xa(t) + ( x)b(t), where A(t) = 2 + η(, t) + u 2 + η, B(t) = 2 + η(, t) + u η. (3.2) (The initial conditions and an appropriate right-hand side were computed by these formulas). The problem was solved by the scheme (2.46) (2.48), (2.5), with h = /N and k = h/. The L 2 errors and rates of convergence for the variables η and u at T = are shown in Table 2. η u N L 2 error rate L 2 error rate e e e e e e e e e e Table 2: L 2 errors and rates of convergence at T =, r = 2, subcritical case, (3.2), h = /N, k/h = / It is clear that Tables and 2 suggest that the L 2 convergence rates are optimal in the case of piecewise linear elements on a uniform mesh. In order to check further the accuracy of the numerical schemes we consider in the supercritical case a problem with a variable bottom having a single hump, and constant initial conditions on (, ) given by β(x) =.4 exp( (x.5) 2 ), η (x) = η =, u (x) = u = 3, (3.3) that we integrate numerically using h = /4, k = h/3. In Figure we show some profiles of the temporal evolution of the numerical solution up to t =.5. The data given by (3.3) and the boundary conditions generate a wave moving to the right and sensing the effect of the variable bottom which is centered at x =.5. There are no spurious oscillations reflected from the boundary x = as the wave exits. By t =.5 the solution has attained a steady state shown in (d). The steady state of such flows is straightforward to determine analytically. Its profile η = η(x), u = u(x) satisfies the equations ((β + η)u) x =, ( η + 2 u2) x =, (3.4) 4
15 (a) η for t = u for t = (b) η for t = u for t = (c) η for t = u for t = (d) η for t = u for t =.5 Figure : Evolution with data (3.3), supercritical case, r = 2, h = /4, k = h/3 5
16 from which using the boundary conditions at x =, we see that u is given in terms of η by u = u (η + β()), (3.5a) η + β where η is the physically acceptable solution of the cubic equation (η + β) 2 ( η η 2 u2 ) + 2 u2 (η + β()) 2 =. (3.5b) (For the analysis of the solutions of the steady-state problem, cf. [6]). We checked the ability of the code to preserve steady-state solutions by taking the profile computed analytically from (3.5) for this problem as initial condition and integrating up to t =.6. The difference between the final profile and the L 2 projection of the analytical initial condition was of O( 9 ) in L 2 for both components when h = /4, k = h/. In Figure 2 we show instances of the temporal evolution up to the attainment of steady state (in (2d)) of the supercritical flow generated with h = /4, k = h/3, by η =, u = 3 and bottom topography and initial conditions given on [, ] by β(x) =.4 exp( (x.75) 2 ), η (x) =.5 exp( 4(x.25) 2 ) + η, u (x) =. exp( 4(x.25) 2 ) + u. (3.6) The variable initial profile gives rise to a wavetrain that moves to the right, interacts with the bottom and exits without spurious oscillations leaving behind the steady state that depends only on η, u and β. We now present some analogous results in the subcritical case. We used the fully discrete scheme with spatial discretization given by (2.46) (2.48), (2.5); the variables depicted in the figures are the approximations of η and u. The spatial discretization was effected on [, ] with piecewise linear functions on a uniform mesh with h = /2; the time-stepping procedure was RK4 as usual with k = h/. In the first example we took η =, u = and β(x) =.4 exp( (x.5) 2 ), η (x) = η, u (x) = u. (3.7) The ensuing evolution of the solution is shown in Figure 3. The generated wave interacts with the bottom and forms pulses that exit without artificial oscillations at both ends of the boundary; the steady-state solution may be found analytically as before. When used as initial condition, its L 2 projection differed from the numerical solution at t = 2 by an L 2 -error of O( 8 ) for this example. An example of subcritical flow with variable initial conditions is shown in Figure 4, where we took η = u =, and β(x) =.4 exp( (x.75) 2 ), η (x) =.5 exp( 4(x.5) 2 ) + η, (3.8) u (x) =. exp( 4(x.5) 2 ) + u 6
17 (a) η for t = u for t = (b) η for t = u for t = (c) η for t = u for t = (d) η for t = u for t =.6 Figure 2: Evolution with data (3.6), supercritical case, r = 2, h = /4, k = h/3 7
18 (a) η for t = u for t = (b) η for t = u for t = (c) η for t = u for t = (d) η for t = u for t = 2 Figure 3: Evolution with data (3.7), subcritical case, r = 2, h = /2, k = h/ 8
19 (a) η for t = u for t = (b) η for t = u for t = (c) η for t = u for t = (d) η for t = 3 (steady state) u for t = 3 (steady state) Figure 4: Evolution with data (3.8), subcritical case, r = 2, h = /2, k = h/ 9
20 and integrated with h = /2, k = h/. A two-way wavetrain emerges and attains steady-state by t = 3. We also tested the code in a few examples of the shallow water equations with absorbing (characteristic) boundary conditions, written in dimensional form, i.e. as η t + ((β + η)u) x =, x L, t T, (3.9) u t + gη x + uu x =, with initial conditions η(x, ) = η (x), u(x, ) = u (x), x L, and analogous characteristic boundary conditions in the super- and subcritical cases. (The Riemann invariants are now u ± g(β + η), g is the acceleration of gravity taken as 9.82 m/s 2, and the bottom is at z = β(x). If the bottom is horizontal it is located at z = h ; in the general case h will be a typical depth.) As an example of supercritical flow we considered a numerical experiment similar to the one described in Section 8.2 of [5]. Let β be the trapezoidal profile given by δ (x L2 ) cκ κ/2 + cκ, if cκ x L/2 κ/2, δ, if κ/2 x L/2 κ/2, β(x) = δ (x L2 ) (3.) cκ κ/2 cκ, if κ/2 x L/2 cκ,, otherwise, where L = 6 m, δ = 5 m, k = L/. The bottom was located at z = β(x), where β(x) = h β(x), h = m, and the problem (3.9) was solved with characteristic boundary conditions and initial conditions η (x) = η = and u (x) = u, where the constant u was varied in order to give flows with different Froude numbers F r = u / gh. We solved (3.9) (3.) numerically with piecewise linear elements and RK4 on a uniform mesh with h = m, k = s. Some profiles of the steady state of the free surface η and the associated bottom function β(x) for various Froude numbers and values of the parameter c are shown in Figure 5. As expected the eventual maximum value of η decreases as F r increases; the results are consistent with those of [5]. In an example of a dimensional subcritical flow we modified the profile given in 5. of [5] in order to avoid discontinuity formation. Thus, the initial η-profile was rounded and its amplitude decreased. Let β be defined by δ β(x) = 2 + δ [ ] π(x L/2) 2 cos, if κ x L 2 < κ, (3.), otherwise, where L = 6 m, δ = 5 m, k = L/. The bottom was taken at z = β(x), where β(x) = h β(x), h = 4 m, and the problem (3.9) was solved with[ characteristic boundary conditions with η = u = and η (x) =.2 ε h exp 5 8 (x 3L/2)/) 2], x L, where ε =.2. 2
21 t = 8sec, c = 3, Fr = 2.2, max = m 6 t = 36sec, c = 2, Fr = 3.3, max = 568.5m m m 5 6 t = 36sec, c = 2, Fr = 4.4, max = m 6 t = 2sec, c =.5, Fr = 4.54, max = m m m 5 Figure 5: Supercritical flows over a trapezoidal bottom, (3.9) (3.), r = 2, h = m, k = s. (Upper figures: steady-state η(x); lower figures: β(x).) The evolution of the η-profiles is shown in Figure 6 up to the attainment of the steady state η = u =. The results resemble qualitatively those of [5]. 3.2 Shallow water equations in balance-law form In this section we consider the numerical solution by the standard Galerkin finite element method of the shallow water equations written in balance-law form (i.e. in conservation-law form with a source term), as d t + (du) x =, (du) t + ( du d2) x = β (x)d, (3.2) where d = η+β is the water depth assumed as always to be positive; the variables in (3.2) are nondimensional. Is is straightforward to see that the system (3.2) is equivalent to (SW) since d. In the sequel we will consider the periodic initial-value problem for (3.2) on the spatial interval [, ] and assume that it has sufficiently smooth solutions for t [, T ], provided that β is smooth and - periodic. We will discretize the problem in space on a uniform or quasiuniform mesh {x i } in [, ] and seek approximations d h, u h of d, u, respectively, in the finite element space S h,p = {φ Cp k : φ [xj,x j+] P r, all j}, where as usual r, k are integers such that r 2, k r 2, and Cp k are the k times continuously differentiable, periodic functions on [, ]. The semidiscrete 2
22 2 t = seconds 2 t = 26 seconds t = 47 seconds 2 t = 85 seconds t = 2 seconds 2 t = 278 seconds t = 4 seconds Figure 6: Subcritical flow over a hump, (3.9), (3.), r = 2, h = m, k = s, η profiles 22
23 approximations satisfy (d ht, φ) + ((d h u h, φ) =, ((d h u h ) t, φ) + ((d h u 2 h + 2 d2 h) x, φ) = (β d h, φ), φ S h,p, t T, (3.3) d h () = P d, u h () = P u, (3.4) where d, u are the initial conditions of d and u and P is now the L 2 projection operator onto S h,p. (The second equation in (3.3) is advanced in time for the variable v h = d h u h and u h is recovered as v h /d h.) In the case of a uniform mesh it is expected that the L 2 errors of the semidiscrete solution will be of O(h r ) while, for a quasiuniform mesh, of O(h r ), cf. []. We verified these rates of accuracy in numerical experiments using C linear, C 2 cubic and C 4 quintic splines (i.e. spaces S h,p with r = 2, 4, and 6, respectively) on uniform and nonuniform spatial meshes, coupled with explicit Runge Kutta schemes of third, fourth, and sixth order of accuracy, respectively. The fully discrete methods were stable under Courant number restrictions. We note that in order to preserve the optimal order of accuracy, say in the case of a uniform mesh, one has to compute the integrals that occur in the finite element equations using, on each subinterval [x i, x i+ ], an s-point Gauss quadrature rule with s r. For example, in the case of a cubic spline spatial discretization, a 3-point Gauss rule is sufficient. It is interesting to examine whether the method (3.3) (3.4) preserves the still water solution η =, u =, e.g. of the periodic i.v.p. for the shallow water equations in the form (3.2). Discretizations that approximate accurately this solution are called well balanced, cf. e.g. [3], and [7] and its references. (It is easy to check that the standard Galerkin semidiscretization of e.g. the periodic ivp for (SW), i.e. for the shallow water equations in their nonconservative form, is trivially well-balanced, since it satisfies η h (x, t) = α, α constant, u h (x, t) = for all t and x [, ], provided η h (x, ) = α, u h (x, ) =. So, our attention is turned to the periodic ivp for (3.2) and its standard Galerkin semidiscretization (3.3) (3.4).) For this purpose, since d = η+β, assume that (suppressing the x-dependence in the variables), d h () = P β, u h () = in (3.4), and ask whether there exist time-independent solutions of (3.3) (3.4) that approximate well the steady state solution d = β, u = of the continuous problem. Taking u h = in (3.3) we see that a steady-state solution d h must satisfy (d h d hx, φ) = (d h β, φ), for all φ in S h,p, from which it is evident that the source term β should be replaced by some approximations β h S h,p thereof. Moreover for the equation (d h d h,x, φ) = (d h β h, φ) to hold for φ S h,p, (this will imply that d h = β h, i.e. good balance), it is necessary that the integrals on each subinterval [x i, x i+ ] that contribute to these L 2 inner products should be evaluated exactly. Since both integrands are polynomials of degree at most 3r 4 on each subinterval, if an s-point Gauss quadrature rule is used (recall that such a rule is exact for polynomials of degree at most 2s ), then it should hold that s 3 2 (r ). For example, in the case of cubic splines (r = 4), a 5-point Gauss rule must be used. Therefore, although a 3-point Gauss is enough to preserve the optimal- 23
24 order O(h 4 ) L 2 -error estimate, good balance of the solution with cubic splines requires that a 5-point Gauss rule be used. This is confirmed by the results of the following experiment. We solve the periodic ivp for (3.2) on [, ] by (3.3) (3.4) using cubic splines for the spatial discretization on a uniform mesh and taking β(x) =.3 exp( (x.5) 2 ), h =.2, k =., u h () =, d h () = P β. Table 3 shows the error d h () d h () (where d h () = d h T = ) in the L 2 and L norms when the analytical formula of β or β h = P β is taken in the source term, and a 3- or a 5-point Gauss rule is used. It is evident that when b h = P β and a 5-point Gauss quadrature rule is used, the scheme is well d h () β in source term s (-point Gauss rule) d h () d h () d h () d h () P β analytical formula 3.89e e-4 P β β h = P β 3.224e e-6 P β β h = P β e-5.24e-4 Table 3: Treatment of source terms and effect of quadrature in (3.3) (3.4), cubic splines and RK4, h =.2, k =., T = balanced to roundoff and there is no influence of the time-stepping error. It should be noted that similar results were found when d h () and β h were taken as the cubic spline interpolant of β at the nodes, and when piecewise smooth bottom profiles, e.g. like a parabolic perturbation of β = supported in the interval of [, ], were considered. Acknowledgement This work was partially supported by the project Innovative Actions in Environmental Research and Development (PErAn) (MIS 52358), implemented under the Action for the strategic development of the Research and Technological sector funded by the Operational Program Competitiveness, and Innovation (NSRF 24-22) and cofinanced by Greece and the EU (European Regional Development Fund). References [] D. H. Peregrine, Equations for water waves and the approximation behind them, in: R. Meyer (Ed.), Waves on Beaches and Resulting Sediment Transport, Academic Press, New York, 972, pp [2] M. Petcu, R. Temam, The one dimensional Shallow Water equations with Dirichlet boundary conditions on the velocity, Discrete Contin. Dyn. Syst. Ser. S 4 () (2)
25 [3] A. Huang, M. Petcu, R. Temam, The one-dimensional supercritical shallowwater equations with topography, Annals of the University of Bucharest (Mathematical Series) 2 (LX) (2) [4] M. Petcu, R. Temam, The one-dimensional shallow water equations with transparent boundary conditions, Mathematical Methods in the Applied Sciences 36 (5) (23) [5] M.-C. Shiue, J. Laminie, R. Temam, J. Tribbia, Boundary value problems for the shallow water equations with topography, J. Geophys. Res. 6 (C25) (2) 22. doi:.29/2jc635. [6] J. Nycander, A. M. Hogg, L. M. Frankcombe, Open boundary conditions for nonlinear channel flow, Ocean Modelling 24 (3-4) (28) 8 2. [7] Y. Xing, X. Zhang, C.-W. Shu, Positivity-preserving high order wellbalanced discontinuous Galerkin methods for the shallow water equations, Advances in Water Resources 33 (2) (2) [8] J. Qiu, Q. Zhang, Stability, error estimate and limiters of discontinuous Galerkin methods, in: R. Abgrall, C.-W. Shu (Eds.), Handbook of Numerical Analysis, Vol. 7, Elsevier, 26, pp [9] Y. Xing, Numerical Methods for the Nonlinear Shallow Water Equations, in: R. Abgrall, C.-W. Shu (Eds.), Handbook of Numerical Analysis, Vol. 8, Elsevier, 27, pp [] D. Antonopoulos, V. Dougalis, Error estimates for the standard Galerkinfinite element method for the shallow water equations, Mathematics of Computation 85 (299) (26) [] D. Antonopoulos, V. Dougalis, Galerkin-finite element methods for the shallow water equations with characteristic boundary conditions, IMA Journal of Numerical Analysis 37 () (27) [2] D. Antonopoulos, V. Dougalis, G. Kounadis, On the standard Galerkin method with explicit RK4 time stepping for the Shallow Water equations, arxiv preprint arxiv:8.8. [3] A. Bermudez, M. E. Vázquez, Upwind methods for hyperbolic conservation laws with source terms, Computers & Fluids 23 (8) (994) [4] J. Douglas, T. Dupont, L. Wahlbin, Optimal L error estimates for Galerkin approximations to solutions of two-point boundary value problems, Mathematics of Computation 29 (3) (975) [5] T. Dupont, Galerkin methods for first order hyperbolics: an example, SIAM Journal on Numerical Analysis (5) (973) [6] D. D. Houghton, A. Kasahara, Nonlinear shallow fluid flow over an isolated ridge, Communications on Pure and Applied Mathematics 2 () (968)
arxiv: v2 [math.na] 28 Feb 2016
NOTES ON GALERKIN-FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS WITH CHARACTERISTIC BOUNDARY CONDITIONS. D.C. ANTONOPOULOS AND V.A. DOUGALIS arxiv:57.829v2 [math.na] 28 Feb 26 Abstract. We consider
More informationarxiv: v2 [math.na] 25 Nov 2014
NOTES ON ERROR ESTIMATES FOR THE STANDARD GALERKIN-FINITE ELEMENT METHOD FOR THE SHALLOW WATER EQUATIONS. D.C. ANTONOPOULOS AND V.A. DOUGALIS arxiv:1403.5699v [mat.na] 5 Nov 014 Abstract. We consider a
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
More informationA modified Galerkin/finite element method for the numerical solution of the Serre-Green-Naghdi system
A modified Galerkin/finite element method for the numerical solution of the Serre-Green-Naghdi system Dimitrios Mitsotakis, Costas Synolakis, Mark Mcguinness To cite this version: Dimitrios Mitsotakis,
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationWell-balanced central finite volume methods for the Ripa system
Well-balanced central finite volume methods for the Ripa system R. Touma a C. Klingenberg b a Lebanese American University, Beirut, Lebanon b Würzburg University, Würzburg, Germany This paper is dedicated
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationCOMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS
COMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS, BUYANG LI, AND CHRISTIAN LUBICH Abstract. We analyze fully implicit and linearly
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationSpatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations
Spatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations N. Chalmers and L. Krivodonova March 5, 014 Abstract We apply the discontinuous Galerkin finite element method
More informationA POSTERIORI ERROR ESTIMATES FOR THE BDF2 METHOD FOR PARABOLIC EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR THE BDF METHOD FOR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND PANAGIOTIS CHATZIPANTELIDIS Abstract. We derive optimal order, residual-based a posteriori error estimates for
More informationNumerical methods for weak solution of wave equation with van der Pol type boundary conditions
Numerical methods for weak solution of wave equation with van der Pol type boundary conditions Jun Liu a, Yu Huang b, Haiwei Sun c, Mingqing Xiao a, a Department of Mathematics, Southern Illinois University,
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationNumerical resolution of discontinuous Galerkin methods for time dependent. wave equations 1. Abstract
Numerical resolution of discontinuous Galerkin methods for time dependent wave equations Xinghui Zhong 2 and Chi-Wang Shu Abstract The discontinuous Galerkin DG method is known to provide good wave resolution
More informationParallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain
Parallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain Keying Ma, 1 Tongjun Sun, 1 Danping Yang 1 School of Mathematics, Shandong University, Jinan 50100, People s Republic
More informationarxiv: v1 [math.na] 3 Jan 2017
ERROR ESTIMATES FOR GALERKIN APPROXIMATIONS OF THE SERRE EQUATIONS DIMITRIOS ANTONOPOULOS, VASSILIOS DOUGALIS, AND DIMITRIOS MITSOTAKIS arxiv:171.668v1 [math.na] 3 Jan 217 Abstract. We consider the Serre
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationStability and Shoaling in the Serre Equations. John D. Carter. March 23, Joint work with Rodrigo Cienfuegos.
March 23, 2009 Joint work with Rodrigo Cienfuegos. Outline The Serre equations I. Derivation II. Properties III. Solutions IV. Solution stability V. Wave shoaling Derivation of the Serre Equations Derivation
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationOn finite time BV blow-up for the p-system
On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationA well-balanced scheme for the shallow-water equations with topography or Manning friction.
A well-balanced scheme for the shallow-water equations with topography or Manning friction. Victor Michel-Dansac a,, Christophe Berthon a, Stéphane Clain b, Françoise Foucher a,c a Laboratoire de Mathématiques
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationPeriodic Solutions of the Serre Equations. John D. Carter. October 24, Joint work with Rodrigo Cienfuegos.
October 24, 2009 Joint work with Rodrigo Cienfuegos. Outline I. Physical system and governing equations II. The Serre equations A. Derivation B. Justification C. Properties D. Solutions E. Stability Physical
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationFinite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water
Journal of Scientific Computing ( 006) DOI: 0.007/s095-005-9060-x Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water Alina Chertock, Alexander Kurganov, and Guergana Petrova
More informationDRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS
DRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS Joseph W. Jerome Department of Mathematics Northwestern University Evanston, IL 60208 Abstract
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationA re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws
A re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c,3 a Department of Applied Mathematics and National
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
More informationNotes on RKDG methods for shallow-water equations in canal networks
Notes on RKDG methods for shallow-water equations in canal networks Maya Briani, Benedetto Piccoli, Jing-Mei Qiu December 9, 015 Abstract PDE models for network flows are used in a number of different
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More information