Two-waves PVM-WAF method for non-conservative systems
|
|
- Hubert Watkins
- 5 years ago
- Views:
Transcription
1 Two-waves PVM-WAF method for non-conservative systems Manuel J. Castro Díaz 1, E.D Fernández Nieto, Gladys Narbona Reina and Marc de la Asunción 1 1 Departamento de Análisis Matemático University of Málaga (Spain) Departmento de Matemática Aplicada University of Sevilla (Spain) Hyp 01, Padova 5-9 June 01
2 Outline 1 Introduction Path-conservative Roe-based schemes PVM methods WAF schemes PVMU-WAF method 3 Numerical tests 4 Conclusions
3 Model problem Let us consider the system w t + F(w) x + B(w) w x = G(w)H x, (1) where w(x, t) takes values on an open convex set O R N, F is a regular function from O to R N, B is a regular matrix function from O to M N N(R), G is a function from O to R N, and H is a function from R to R. By adding to (1) the equation H t = 0, the system (1) can be rewritten under the form where W is the augmented vector and W = W t + A(W) W x = 0, () [ w H ] Ω = O R R N+1
4 Model problem Let us consider the system w t + F(w) x + B(w) w x = G(w)H x, (1) where w(x, t) takes values on an open convex set O R N, F is a regular function from O to R N, B is a regular matrix function from O to M N N(R), G is a function from O to R N, and H is a function from R to R. By adding to (1) the equation H t = 0, the system (1) can be rewritten under the form where A(W) is the matrix whose block structure is given by: [ A(w) G(w) A(W) = 0 0 W t + A(W) W x = 0, () ], where A(w) = J(w) + B(w), being J(w) = F w (w).
5 Difficulties Main difficulties Non conservative products A(W) W x. Solutions may develop discontinuities and the concept of weak solution in the sense of distributions cannot be used. The theory introduced by DLM 1995 is used here to define the weak solutions of the system. This theory allows one to give a sense to the non conservative terms of the system as Borel measures provided a prescribed family of paths in the space of states. Derivation of numerical schemes for non-conservative systems: path-conservative numerical schemes (Parés 006). The eigenstructure of systems like bilayer Shallow-Water system or two-phase flow model of Pitman Le are not explicitly known: PVM and/or WAF schemes.
6 PC-Roe-based schemes I Let us consider path-conservative numerical schemes that can be written as follows: where w n+1 i = w n i t ( D + i 1/ x + i+1/) D, (3) x and t are, for simplicity, assumed to be constant; w n i is the approximation provided by the numerical scheme of the cell average of the exact solution at the i-th cell, I i = [x i 1/, x i+1/ ] at the n-th time level t n = n t. H i is the cell average of the function H(x).
7 PC-Roe-based schemes I Let us consider path-conservative numerical schemes that can be written as follows: w n+1 i = w n i t ( D + i 1/ x + i+1/) D, (3) D ± i+1/ = 1 ( F(wi+1) F(w i) + B i+1/ (w i+1 w i) G i+1/ (H i+1 H i) ) ± Q i+1/ (w i+1 w i A 1 i+1/ G i+1/(h i+1 H i)), (4) where A i+1/ = J i+1/ + B i+1/. Here, J i+1/ is a Roe matrix of the Jacobian of the flux F in the usual sense: B i+1/ (w i+1 w i) = J i+1/ (w i+1 w i) = F(w i+1) F(w i); B(Φ w(s; W i, W i+1)) Φw (s; Wi, Wi+1) ds; s G i+1/ (H i+1 H i) = G(Φ w(s; W i, W i+1)) ΦH (s; Wi, Wi+1) ds; 0 s Q i+1/ is a numerical viscosity matrix.
8 PC-Roe-based schemes I Let us consider path-conservative numerical schemes that can be written as follows: w n+1 i = w n i t ( D + i 1/ x + i+1/) D, (3) D ± i+1/ = 1 ( F(wi+1) F(w i) + B i+1/ (w i+1 w i) G i+1/ (H i+1 H i) ) ± Q i+1/ (w i+1 w i A 1 i+1/ G i+1/(h i+1 H i)), (4) Conservative systems If the system is conservative and F i+1/ = F(wi) + F(wi+1) 1 Q i+1/(w i+1 w i) is a conservative flux, where Q i+1/ is defined in terms of J i+1/, then D i+1/ = F i+1/ F(w i) D + i+1/ = F(wi+1) F i+1/.
9 PC-Roe-based schemes I Let us consider path-conservative numerical schemes that can be written as follows: w n+1 i = w n i t ( D + i 1/ x + i+1/) D, (3) D ± i+1/ = 1 ( F(wi+1) F(w i) + B i+1/ (w i+1 w i) G i+1/ (H i+1 H i) ) ± Q i+1/ (w i+1 w i A 1 i+1/ G i+1/(h i+1 H i)), (4) Different numerical schemes can be obtained for different definitions of Q i+1/
10 PC-Roe-based schemes II Roe scheme corresponds to the choice Q i+1/ = A i+1/, Lax-Friedrichs scheme: being Id the identity matrix. Lax-Wendroff scheme: Q i+1/ = x t Id, Q i+1/ = t x A i+1/, FORCE and GFORCE schemes: Q i+1/ = (1 ω) x t Id + ω t x A i+1/, with ω = 0.5 and ω = 1, respectively, being α the CFL parameter. 1+α
11 PVM methods We propose a class of finite volume methods defined by being P l (x) a polinomial of degree l, Q i+1/ = P l (A i+1/ ), P l (x) = l j=0 α i+1/ j x j, and A i+1/ a Roe matrix. That is, Q i+1/ can be seen as a Polynomial Viscosity Matrix (PVM). See also: P. Degond, P.F. Peyrard, G. Russo, Ph. Villedieu. Polynomial upwind schemes for hyperbolic systems. C. R. Acad. Sci. Paris 1 38, , 1999.
12 PVM methods We propose a class of finite volume methods defined by being P l (x) a polinomial of degree l, Q i+1/ = P l (A i+1/ ), P l (x) = l j=0 α i+1/ j x j, and A i+1/ a Roe matrix. That is, Q i+1/ can be seen as a Polynomial Viscosity Matrix (PVM). Q i+1/ has the same eigenvectors than A i+1/ and if λ i+1/ is an eigenvalue of A i+1/, then P l (λ i+1/ ) is an eigenvalue of Q i+1/.
13 PVM methods We propose a class of finite volume methods defined by being P l (x) a polinomial of degree l, Q i+1/ = P l (A i+1/ ), P l (x) = l j=0 α i+1/ j x j, and A i+1/ a Roe matrix. That is, Q i+1/ can be seen as a Polynomial Viscosity Matrix (PVM). Some well-known solvers as Lax-Friedrichs, Rusanov, FORCE/GFORCE, HLL, Roe, Lax-Wendroff,... can be recovered as PVM methods
14 PVM-1U(S L, S R ) or HLL method P 1(x) = α 0 + α 1 x such as P 1(S L) = S L, P 1(S R) = S R. Q i+1/ = α 0Id + α 1A i+1/ PVM 1U(S L,S R ) S L λ 1 λ λ j... λ N S R
15 PVM-1U(S L, S R ) or HLL method The usual HLL scheme coincides with PVM-1U(S L, S R) in the case of conservative systems. Let us suppose that the system is conservative. Then, the conservative flux associated to PVM-1U(S L,S R) is F i+1/ = D i+1/ + F(wi). Taking into account that then F i+1/ = α 0 = SR SL SL SR SR SL, α 1 =, S R S L S R S L F(wi)(SR + SR SL SL ) + F(wi+1)(SR SR SL + SL ) S R S L (SR SL SL SR )(wi+1 wi) S R S L = S+ R F(wi) S L F(wi+1) + (S+ R S L )(wi+1 wi) S + R S L which is a compact definition of the HLL flux, being S + R = min(sl, 0). S L = max(sr, 0) and
16 PVM-U(S L, S R ) method such as where P (x) = α 0 + α 1x + α x, P (S m) = S m, P (S M) = S M, P (S M) = sgn(s M), { SL if S L S R, S M = if S L < S R. S R { SR if S L S R, S m = if S L < S R. S L PVM U(S L,S R ) SL λ1 λ λj... λn SR
17 WAF method Let us consider the following Riemann problem: w + F(w) = 0; t x w(x, 0) = { wi x < 0; w i+1 x > 0. We denote by S i for i = 1,, N some approximation of the characteristic velocities. S 1 t S! t "! t / " 1 " 3 0 #! x /! x / x F WAF i+1/ = 1 x x/ x/ F(w(x, t ))dx.
18 where F (k) i+1/ is the value of the flux function in the interval k, WAF method (Toro [1989]) S 1 t S! t "! t / " 1 " 3 0 #! x /! x / x ω k = 1 (c k c k 1 ), c 0 = 1, c N+1 = 1 and c l = t x S l, for 1 l N, N is the number of waves. So, we can write the numerical flux as follows: F WAF i+1/ = 1 (Fi + Fi+1) 1 N k=1 c k F (k) i+1/,
19 WAF method TVD WAF method (Toro [1989]) We denote χ(v) a flux limiter function, and Ψ(v, c) = 1 (1 c )χ(v). So the TVD-WAF flux function is defined by: where F WAF i+1/ = 1 (Fi + Fi+1) 1 N k=1 sign(s k )Ψ k F (k) i+1/, Ψ k = Ψ(v (k), c k ) = 1 (1 c k )χ(v (k) ). Some suitable choices for χ can be found in [Toro]. Let us consider here the Van Albada s limiter: p (k) 0 if v (k) i p (k) i 1 0 χ(v (k) ) = v (k) (1 + v (k) ) if v (k) 0, where v (k) p (k) i+1 = p(k) i p (k) 1 + v (k) i+ p(k) i+1 p (k) i+1 p(k) i being p (k) a scalar value. if S k > 0 if S k < 0,
20 HLL-WAF method We consider now N =, S 1 = S L, S = S R, F 1 i+1/ = F(w i), F 3 i+1/ = F(wi+1) and F () SRFi SLFi+1 + SRSL(wi+1 wi) i+1/ = S R S L F HLL-WAF i+1/ = 1 (Fi + Fi+1) 1 (ν1(χl, χr)(wi+1 wi) + ν(χl, χr)(fi+1 Fi)) 1 t (µ1(χl, χr)(wi+1 wi) + µ(χl, χr)(fi+1 Fi)), x where ν 1(χ L, χ R) = ν (χ L, χ R) = µ 1(χ L, χ R) = µ (χ L, χ R) = S Rχ R S Lχ L S R S L. SLSR((1 χl)sgn(sl) (1 χr)sgn(sr)) S R S L (1 χr) SR (1 χl) SL S R S L SLSR(SLχL SRχR) S R S L
21 HLL-WAF method as a PVM-type method Then, we can rewrite the HLL-WAF method as follows: F HLL-WAF i+1/ = 1 (Fi+1 + Fi) 1 QHLL WAF i+1/ (w i+1 w i), with Q HLL WAF i+1/ (χ L, χ R) = Q HLL WAF t o1,i+1/ (χl, χr) + x QHLL WAF o,i+1/ (χl, χr) Q HLL WAF o1,i+1/ (χl, χr) = ν1(χl, χr)i + ν(χl, χr)a i+1/ Q HLL WAF o,i+1/ (χl, χr) = µ1(χl, χr)i + µ(χl, χr)a i+1/. That is, the usual two-waves HLL-WAF method can be seen as a non-linear combination of two PVM schemes associated to the first order polynomials: P o1 1 (x) = ν 1(χ L, χ R) + ν (χ L, χ R)x and P o 1 (x) = µ 1(χ L, χ R) + µ (χ L, χ R)x.
22 HLL-WAF method as a PVM-type method y=p 1 HLL (x) y=p 1 LW (x) S1 λ1 λ λj... λn S
23 HLL-WAF method as a PVM-type method y=p 1 HLL (x) y=p 1 LW (x) S1 λ1 λ λj... λn S Remarks For systems with N=. If S 1 = λ 1,i+1/ and S = λ,i+1/ being λ j,i+1/ the eigenvalues of Roe matrix, then HLL-WAF method coincides with Lax-Wendroff method if (χ L = χ R = 1). For N > it is not true! HLL-WAF method coincides with HLL if (χ L = χ R = 0) (N ).
24 HLL-WAF method as a PVM-type method y=p 1 HLL (x) y=p 1 LW (x) S1 λ1 λ λj... λn S Objective We want to define a new two-wave WAF method so that coincides with Lax-Wendroff for N > if χ L = χ R = 1
25 Two-waves PVMU-WAF methods We consider F U-FL Fi + Fi+1 i+1/ = where Q U FL i+1/ (χ L, χ R ) is defined as follows: with 1 QU FL i+1/ (χl, χr)(wi+1 wi), Q U FL i+1/ (χ L, χ R ) = Q U FL o1,i+1/ (χ L, χ R ) + t x QU FL o,i+1/ (χ L, χ R ), Q U FL o1,i+1/ (χ L, χ R ) = sgn(s L)(1 χ L ) + sgn(s R )(1 χ R ) A i+1/ and + sgn(s R)(1 χ R ) P,αR (A i+1/ ) sgn(s L)(1 χ L ) P,αL (A i+1/ ), Q U FL o,i+1/ (χ L, χ R ) = S Lχ L + S R χ R A i+1/ + S Rχ R P,αR (A i+1/ ) S Lχ L P,α L (A i+1/ ), where α K = 1 (1 χ K )(1 α), K = L, R, α = (S R S L )sgn(s M ) (S R + S L ), 4S M (S L + S R )
26 Two-waves PVMU-WAF methods We consider F U-FL Fi + Fi+1 i+1/ = where Q U FL i+1/ (χ L, χ R ) is defined as follows: with 1 QU FL i+1/ (χl, χr)(wi+1 wi), Q U FL i+1/ (χ L, χ R ) = Q U FL o1,i+1/ (χ L, χ R ) + t x QU FL o,i+1/ (χ L, χ R ), Q U FL o1,i+1/ (χ L, χ R ) = sgn(s L)(1 χ L ) + sgn(s R )(1 χ R ) A i+1/ and + sgn(s R)(1 χ R ) P,αR (A i+1/ ) sgn(s L)(1 χ L ) P,αL (A i+1/ ), Q U FL o,i+1/ (χ L, χ R ) = S Lχ L + S R χ R A i+1/ + S Rχ R P,αR (A i+1/ ) S Lχ L P,α L (A i+1/ ), and P,α (x) = αp M (x) + (1 α)p 1M (x), P 1M (x) = S RS L S R S L + S R + S L S R S L x.
27 Two-wave PVM-U-WAF method Properties Only uses the information of the two fastest waves. If N = it coindices with the usual HLL-WAF scheme. If χ L = χ R = 0, then we recover the PVM-U first order scheme for N. If χ L = χ R = 1, then scheme reduces to Lax-Wendroff scheme for N.
28 Two-wave PVM-U-WAF method Remark A natural extension to balance laws and non-conservative system is straightforward: with w n+1 i = w n i t ( D + i 1/ x + i+1/) D, D ± i+1/ = 1 ( F(wi+1) F(w i) + B i+1/ (w i+1 w i) G i+1/ (H i+1 H i) ) ± Q i+1/ (w i+1 w i A 1 i+1/ G i+1/(h i+1 H i)), being But it is not second order. Q i+1/ = Q U FL i+1/ (χl, χr)
29 Two-wave PVM-U-WAF method To recover the second order, new terms appear in the Lax-Wendroff scheme due to the non-conservative products and source terms (see Castro, Pares & Toro Math. Comp. 010): w n+1 i = w n i t ( D + i 1/ x + ) t D i+1/ + 4 x (R(χL, χr)n i 1/ + R(χ L, χ R) n i+1/) with R(χ L, χ R ) n i+1/ = 1 ( χ L DA(W i )[A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i ), w i+1 w i ] + χ R DA(W i+1 )[A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i ), w i+1 w i ] χ L DA(W i )[w i+1 w i, A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i )] χ R DA(W i+1 )[w i+1 w i, A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i )] χ L G w(w i )(A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i ))(H i+1 H i ) ) χ R G w(w i+1 )(A i+1/ (w i+1 w i ) G i+1/ (H i+1 H i ))(H i+1 H i ) ( ) N being DA(W)[U, V] = l=1 u l wl A(W) V and wl A(W) is the N N matrix whose (i, j) element is wl a ij (W). G w(w) denotes the Jacobian matrix of G(w).
30 Multilayer shallow water system th j + xq j = 0, tq j + x( q j + 1 h j gh j ) + gh j x(z b + h k + ρ k h k ) = 0. ρ j k>j k<j j = 1,..., m, where m is the number of layers, h j, j = 1,..., m are the fluid depths, q j = h ju j are the discharges, u j are the velocites and z b (x) is the topography. g is the gravity constant and ρ j the densisites of the stratified fluid layers, with 0 < ρ 1 < < ρ m. We use the Van Albada s limiter with a smooth indicator of the fluid interfaces. Total energy of the system provides also good results.
31 1LSW: Stationary subcritical solution over bump I = [0, 0]. z b (x) = 0.e 0.16(x 10) Boundary conditions: q(0, t) = 4.4 and h(0, t) =.0 x = 1/0. cfl = 0.9. Nodes L 1 err h L 1 order h L 1 err q L 1 order q Table: Errors and order. Subcritical stationary solution.
32 1LSW: Order of accuracy I = [0, 1], T = 0.1. z b (x) = sen (πx). Initial condition: h(x, 0) = 5 + e cos(πx), q(x, 0) = sin(cos(πx)), CFL = 0.8. Reference solution computed with ROE scheme with x = 1/1800. Nodes L 1 err h L 1 order h L 1 err q L 1 order q Table: Errors and order.
33 1LSW stationary transcritical flow with a shock I = [0, 0] bottom topography: { (x 10) 8 < x < 1 z b (x) = 0 otherwise. Initial condition h(x, 0) = 0.33 z b, q(x, 0) = 0. Boundary conditions: q = 0.18 at x=0 and h = 0.33 at x = 0 x = 1/10 and CFL = 0.9. Concerning CPU time similar results are obtained for Roe solver and PVM-U WAF method.
34 1LSW stationary transcritical flow with a shock Exact solution 0.1 Roe PVM U WAF 0.05 HLL WAF PVM U Bottom (a) Free surface an bottom topography
35 LSW: Internal dam break problem. I = [0, 10]. zb(x) = 0 Initial condition: q 1(x, 0) = q { (x, 0) = 0, x [0, 10], 0.9 if x < 5, h 1(x, 0) = 0.1 if x 5, h (x, 0) = 1.0 h 1(x, 0) x [0, 10]. Open boundary conditions x = 1/0. A reference solution is computed Roe scheme with x = 1/00. r = 0.99, cfl = 0.9. Concerning CPU time PVM-U-WAF method is.8 times faster than Roe and similar to original HLL-WAF scheme.
36 LSW: Internal dam break problem Reference solution Roe PVM U WAF (b) PVM-U-WAF and Roe Reference solution HLL WAF PVM U (c) HLL-WAF and PVM-U Figure: Internal dam break: free surface and interface at t = 0 seg.
37 LSW: Stationary transcritical flow with an internal shock I = [0, 10] Bottom topography: z b (x) = 0.5e (x 5). Initial condition: q 1(x, 0) = q (x, 0) = 0. and h 1(x, 0) = { 0.48 if x < 5, 0.5 if x 5, h (x, 0) = 1 h 1(x, 0) z b (x), ρ 1/ρ = Free boundary conditions. CFL = 0.9, x = 1/0. Reference solution computed with x = 1/00.
38 LSW: Stationary transcritical flow with an internal shock Reference solution 0. Roe PVM U WAF HLL WAF PVM U bottom (a) Free surface, bottom topography and interface
39 4LSW: Internal dam breaks. I = [0, 10] Bottom topography: z b (x) = 0.0 Initial condition: q i(x, 0) = 0., i = 1,..., 4 and { 0.9 if x < 5, h 1(x, 0) = 0.1 if x 5, h (x, 0) = 1 h 1(x, 0), h 3(x, 0) = h 1(x, 0), h 4(x, 0) = h (x, 0) ρ 1/ρ 4 = 0.85, ρ /ρ 4 = 0.9, ρ 3/ρ 4 = Free boundary conditions. CFL = 0.9, x = 1/0. Reference solution computed with x = 1/00. Concerning CPU time PVM-U-WAF method is 9.8 times faster than Roe and similar to original HLL-WAF scheme.
40 4LSW: Internal dam breaks Reference solution Roe 0. Reference solution HLL WAF PVM U WAF PVM U (b) PVM-U-WAF and Roe (c) HLL-WAF and PVM-U Figure: Internal dam breaks: free surface and interfaces at t = 5 seg.
41 D 1LSW: Circular dam break D = [, ] [, ] Bottom topography: zb(x, y) = 0.8 e x y Initial condition: q x(x, y, 0) = q y(x, y, 0) = 0 and { 1 zb(x, y) if x + y h(x, y, 0) = < zb(x, y) otherwise Wall boundary conditions x = y = 1/100, CFL = 0.9. Three implementations are considered: first and second order HLL and PVM-U-WAF method. Algorithms implemented on GPUs: speedups of more than 00 for the three numerical schemes, The extension by the method of lines of 1D PVM-U-WAF method to multidimensional problems is NOT second order accurate.
42 D 1LSW: Circular dam break (a) PVM-U WAF t = 1.0 s (b) HLL t = 1.0 s Figure: D circular dam break: free surface at t = 1 seg.
43 D 1LSW: Circular dam break (a) PVM-U WAF t = 1.0 s (b) Second order HLL t = 1.0 s Figure: D circular dam break: free surface at t = 1 seg.
44 D 1LSW: Circular dam break (a) PVM-U WAF t =.0 s (b) HLL t =.0 s Figure: D circular dam break: free surface at t = 1 seg.
45 D 1LSW: Circular dam break (a) PVM-U WAF t =.0 s (b) Second order HLL t =.0 s Figure: D circular dam break: free surface at t = seg.
46 Conclusions Conclusions The original two-wave HLL-WAF method can be seen as a PVM-based flux-limiting scheme. A new two-wave WAF method that ensured second order of accuracy for N > is defined using PVM framework. It can be applied to conservative, balance laws and non-conservative systems. Its performance increases with the complexity of the system. It can be 10 times faster than Roe solver for the 1D 4LSW. Extension to D that preserves second order accuracy: comming soon, it is NOT straight forward.
HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF STATES FOR SOLVING HYPERBOLIC SYSTEMS WITH NONCONSERVATIVE PRODUCTS
MATHEMATICS OF COMPUTATION Volume 75, Number 255, July 26, Pages 113 1134 S 25-5718(6)1851-5 Article electronically published on March 21, 26 HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF
More informationWell-Balanced Finite Volume schemes for 2D non-homogeneous hyperbolic systems. Application to the Dam-break of Aznalcóllar.
Well-Balanced Finite Volume schemes for D non-homogeneous hyperbolic systems. Application to the Dam-break of Aznalcóllar. M. J. Castro Díaz, T. Chacón Rebollo, E. D. Fernández-Nieto, J.M. González Vida
More informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
More informationWELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL
WELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL VIVIEN DESVEAUX, MARKUS ZENK, CHRISTOPHE BERTHON, AND CHRISTIAN KLINGENBERG Abstract. The present paper concerns the derivation
More informationarxiv: v1 [physics.comp-ph] 26 Oct 2018
arxiv:1810.11285v1 [physics.comp-ph] 26 Oct 2018 Efficient analytical Roe solver for two-layer shallow-water equations with accurate treatment for loss of hyperbolicity Nino Krvavica a,, Miran Tuhtan a,
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 informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations 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,
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationCommissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation. MIT, November 3rd, 2008
ON THE NUMERICAL SIMULATION OF GENERAL MULTI-FIELD TWO-PHASE FLOW MODELS Anela KUMBARO Commissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation MIT, November
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 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 informationSistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationA Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids
Diese Arbeit wurde vorgelegt am Lehrstuhl für Mathematik (MathCCES) A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids CES Seminar Computational Engineering
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
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 informationDiscontinuous Galerkin Method for 1D Shallow Water Flow with Water Surface Slope Limiter
Clemson University TigerPrints All Theses Theses 1-010 Discontinuous Galerkin Method for 1D Shallow Water Flow with Water Surface Slope Limiter Wencong Lai Clemson University, wlai@clemson.edu Follow this
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationAdmissibility and asymptotic-preserving scheme
Admissibility and asymptotic-preserving scheme F. Blachère 1, R. Turpault 2 1 Laboratoire de Mathématiques Jean Leray (LMJL), Université de Nantes, 2 Institut de Mathématiques de Bordeaux (IMB), Bordeaux-INP
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationUn schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme
Un schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme Christophe Berthon, Anaïs Crestetto et Françoise Foucher LMJL, Université de Nantes ANR GEONUM Séminaire de l équipe
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 informationHigh Performance meshfree methods for fluid flows Computation
High Performance meshfree methods for fluid flows Computation F. Benkhaldoun 1 A. Halassi 2 D. Ouazar 3 M. Seaid 4 A. Taik 2 1 LAGA, Unstitut Galillé, Université Paris 13, Villetaneuse, France 2 LaboMAC
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
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 informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationA numerical study of a particular non-conservative hyperbolic problem
A numerical study of a particular non-conservative hyperbolic problem V. Dolejší a,b,, T. Gallouët b a Charles University Prague, Faculty of Mathematics and Physics, Soolovsá 83, 186 75 Prague, Czech Republic
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 informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationDEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 26 DEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES V.
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationA Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws
A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract
More informationManuel Castro 1,JorgeMacías 1 and Carlos Parés 1
Mathematical Modelling and Numerical Analysis ESAIM: MAN Modélisation Mathématique et Analyse Numérique Vol 35, N o,, pp 7 7 A Q-SCHEME FOR A CLASS OF SYSTEMS OF COUPLED CONSERVATION LAWS WITH SOURCE TERM
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationA kinetic scheme for transient mixed flows in non uniform closed pipes: a global manner to upwind all the source terms
A kinetic scheme for transient mixed flows in non uniform closed pipes: a global manner to upwind all the source terms C. BOURDARIAS M. ERSOY S. GERBI LAMA-Université de Savoie de Chambéry, France September
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationFlux Limiter Lattice Boltzmann for Compressible Flows
Commun. Theor. Phys. 56 (2011 333 338 Vol. 56, No. 2, August 15, 2011 Flux Limiter Lattice Boltzmann for Compressible Flows CHEN Feng (íô, 1 XU Ai-Guo (ÆÇÁ, 2, ZHANG Guang-Cai (¾, 2 and LI Ying-Jun (Ó
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 informationWELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES: RIPA MODEL
WELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES: RIPA MODEL VIVIEN DESVEAUX, MARKUS ZENK, CHRISTOPHE BERTHON, AND CHRISTIAN KLINGENBERG Abstract. The present paper concerns the derivation of
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 informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More informationComparison of Approximate Riemann Solvers
Comparison of Approximate Riemann Solvers Charlotte Kong May 0 Department of Mathematics University of Reading Supervisor: Dr P Sweby A dissertation submitted in partial fulfilment of the requirement for
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationOn the entropy stability of Roe-type finite volume methods
Proceedings of Symposia in Applied Mathematics On the entropy stability of Roe-type finite volume methods Mária Lukáčová - Medvid ová and Eitan Tadmor Abstract. We study the entropy stability of a class
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
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 informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
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 informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
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 informationA CLASS OF COMPUTATIONALLY FAST FIRST ORDER FINITE VOLUME SOLVERS: PVM METHODS.
A CLASS OF COMPUTATIONALLY FAST FIRST ORDER FINITE VOLUME SOLVERS: PVM METHODS. M. J. CASTRO DíAZ AND E. D. FERNÁNDEZ-NIETO Abstract. In this work, we present a cass of fast first order finite voume sovers,
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationEmpirical versus Direct sensitivity computation
Empirical versus Direct sensitivity computation Application to the shallow water equations C. Delenne P. Finaud Guyot V. Guinot B. Cappelaere HydroSciences Montpellier UMR 5569 CNRS IRD UM1 UM2 Simhydro
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationNumerical and mathematical analysis of a five-equation model for two-phase flow
SINTEF Energy Research Master s Thesis PFE-Master Numerical and mathematical analysis of a five-equation model for two-phase flow Author Pedro José Martínez Ferrer Supervisors Svend Tollak Munkejord Tore
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationOn the influence of the thickness of the sediment moving layer in the definition of the bedload transport formula in Exner systems
On the influence of the thickness of the sediment moving layer in the definition of the bedload transport formula in Exner systems E.D. Fernández-Nieto a, C. Lucas b, T. Morales de Luna c,, S. Cordier
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
More informationStabilized Spectral Element Approximation of the Saint Venant System Using the Entropy Viscosity Technique
Stabilized Spectral Element Approximation of the Saint Venant System Using the Entropy Viscosity Technique R. Pasquetti, J.L. Guermond, and B. Popov Abstract We consider the Saint Venant system (shallow
More informationSimulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes
J Sci Comput (2011) 47: 281 302 DOI 10.1007/s10915-010-9437-3 Simulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes Changna Lu Jianxian
More informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
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 informationScienceDirect. A pipe network simulation model with dynamic transition between free surface and pressurized flow
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 70 ( 2014 ) 641 650 12th International Conference on Computing and Control for the Water Industry, CCWI2013 A pipe network simulation
More informationJournal of Computational Physics
Journal of Computational Physics 8 (9) 77 798 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A central Rankine Hugoniot solver
More informationON THE COMPARISON OF EVOLUTION GALERKIN AND DISCONTINUOUS GALERKIN SCHEMES
1 ON THE COMPARISON OF EVOLUTION GALERKIN AND DISCONTINUOUS GALERKIN SCHEMES K. BAUMBACH and M. LUKÁČOVÁ-MEDVIĎOVÁ Institute of Numerical Simulation, Hamburg University of Technology, 21 079 Hamburg, Germany,
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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationA Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms
A Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms Randall J. LeVeque 1 Abstract. Systems of hyperbolic partial differential equations with source terms (balance laws)
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationApproximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD
Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD Andrea Mignone Collaborators: G. Bodo, M. Ugliano Dipartimento di Fisica Generale, Universita di Torino (Italy)
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationHyperconcentrated 1D Shallow Flows on Fixed Bed with Geometrical Source Term Due to a Bottom Step
J Sci Comput (2011) 48:319 332 DOI 10.1007/s10915-010-9457-z Hyperconcentrated 1D Shallow Flows on Fixed Bed with Geometrical Source Term Due to a Bottom Step L. Begnudelli G. Rosatti Received: 24 November
More informationFinite volume method for conservation laws V
Finite volume method for conservation laws V Schemes satisfying entropy condition Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationSome approximate Godunov schemes to compute. shallow-water equations with topography. Universite de Provence. 39 rue Joliot Curie.
Some approximate Godunov schemes to compute shallow-water equations with topography Thierry Gallouet, Jean-Marc Herard ;, Nicolas Seguin ; Laboratoire d'analyse Topologie et Probabilites - UMR 663 Centre
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationNotes: Outline. Shallow water equations. Notes: Shallow water equations. Notes:
Outline Nonlinear hyperbolic systems Shallow water equations Shock waves and Hugoniot loci Integral curves in phase plane Compression and rarefaction R.J. LeVeque, University of Washington IPDE 2011, July
More informationState of the Art MHD Methods for Astrophysical Applications p.1/32
State of the Art MHD Methods for Astrophysical Applications Scott C. Noble February 25, 2004 CTA, Physics Dept., UIUC State of the Art MHD Methods for Astrophysical Applications p.1/32 Plan of Attack Is
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationDiscontinuous Galerkin methods Lecture 2
y y RMMC 2008 Discontinuous Galerkin methods Lecture 2 1 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu y 1 0.75 0.5 0.25 0-0.25-0.5-0.75 y 0.75-0.0028-0.0072-0.0117 0.5-0.0162-0.0207-0.0252
More informationCentral-Upwind Schemes for Shallow Water Models. Alexander Kurganov. Tulane University Mathematics Department
Central-Upwind Schemes for Shallow Water Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov Supported by NSF and ONR Shallow Water Equations Shallow water
More informationFLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS
FLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS M.E.Hubbard The University of Cambridge, DAMTP, Silver Street, Cambridge, CB3 9EW, United Kingdom and P.Garcia-Navarro Area
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationThe Euler Equations! Advection! λ 1. λ 2. λ 3. ρ ρu. c 2 = γp/ρ. E = e + u 2 /2 H = h + u 2 /2; h = e + p/ρ. 0 u 1/ρ. u p. t + A f.
http://www.nd.edu/~gtryggva/cfd-course/! Advection! Grétar Tryggvason! Spring! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ E = e + u / H = h + u
More informationApplication of the relaxat ion met hod to model hydraulic jumps
Application of the relaxat ion met hod to model hydraulic jumps P. J. Montgomery Mathematics and Computer Science Program, University of Northern British Columbia, Prince George, Canada. Abstract A finite
More information