Non-linear Methods for Scalar Equations
|
|
- Lorraine Lawson
- 5 years ago
- Views:
Transcription
1 Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it October 3, 04 / 56
2 Abstract We construct non-linear second order methods for scalar equations, including source terms Two non-linear constraints are enforced: TVD and ENO Approaches: flux limiter methods and reconstruction methods Schemes derived: MUSCL-Hancock (one upwind and one centred), WAF and ADER / 56
3 Flux Limiter Methods We study a class of TVD schemes called flux limiter methods, as applied to The schemes have conservative form qi n+ = qi n t x with numerical flux f LO i+ ψ i+ t q + x f(q) = 0, f(q) = λq. () f i+ : a monotone flux; f HO i+ f LO i+ f HO i+ = f GodU i+ = f LW i+ = f LO + ψ i+ i+ ( f i+ f i ) () (f HO f LO ). (3) i+ i+ : flux for second-order linear scheme. Here { λq n = i if λ > 0, λqi+ n if λ < 0, = ( + c)λqn i + ( c)λqn i+. is called flux limiter and is constructed on TVD considerations. (4) 3 / 56
4 We construct the TVD region, which is the range for the limiter function. First we consider the case of positive characteristic speed, λ > 0. ψ i+ Substitution of (4) into (3) gives f i+ = λq n i + ( c)λ(qn i+ q n i )ψ i+ Inserting this flux into the conservative formula () gives qi n+ = qi n c[ q i + ( c)(ψ i+ q i+ which if written in incremental form ψ i. (5) q i )] (6) has coefficients q n+ i = q n i C i q i + D i+ q i+ (7) C i = c[ ( c)ψ i ] ; D i+ = ( c)cψ i+ (8) 4 / 56
5 We note that (7) can be simplified by introducing the definition qi n qn i r i+ = q qi+ n = q i if λ > 0, qn i q i+ upw = q loc Then for λ > 0 we have q n i+ qn i+ q n i+ qn i = q i+ 3 q i+ if λ < 0. (9) which if inserted in (7) gives with Ĉ i = c ( q i+ = q i r i+ q n+ i = q n i Ĉi q i + ( c) ( ψi+ r i+, (0) + ˆD i+ q i+ ψ i )), ˆDi+, () = 0. () 5 / 56
6 Recall that Harten s theorem requires C i+ 0, D i+ 0, 0 C i+ + D i+. (3) For this to apply we simply need to enforce the following inequalities ( ( )) 0 Ĉi = c + ψi+ ( c) ψ r i. (4) i+ From now on we assume the following condition on the Courant number 0 c. (5) Then, enforcing the left-hand side and the right-hand side of (4) leads to c ψ i+ r i+ ψ i c. (6) 6 / 56
7 We impose the following conditions on the limiter functions 0 ψ i, i. (7) c Then from (6) and (7) we can deduce the following pair of inequalities c ψ i+ r i+, i (8) c( c) Inequalities (7) give lower and upper bounds for the limiter functions. Inequalities (8) determine two inclined straight lines of slopes c and ( c)c. The analysis for λ < 0 is analogous, so that c is replaced by c in all equations. The resulting general TVD region in the r-ψ plane is depicted in Fig.. 7 / 56
8 ψ(r; c ) ( ) ( ψ L = c r ψr = ( c ) c ) r ψ = c TVD TVD ψ = 0 r = r = 0 r = c Fig. : General TVD region for the flux limiter method based on the Godunov upwind and Lax-Wendroff methods. r 8 / 56
9 The complete general TVD region is given by 0 ψ i+ ψ max (r), i, (9) with ψ max (r) = Negative r Positive r c if r, r if r 0, ( c ) r if 0 r c, ( c ) c c if r c. (0) 9 / 56
10 The Sweby TVD region [0] is a subset of the general TVD region derived above and is obtained by defining 0 if r 0, ψ Sweby max (r) = r if 0 r, () if r in (9). The Sweby TVD region is depicted by the shaded zone of Fig., to the right of ψ = r, above ψ = 0 and below ψ = ; also, ψ = 0 for r 0 is part of the TVD region. Any function ψ(r) within this region, called a flux limiter, produces a TVD scheme Not every choice produces a good TVD scheme. 0 / 56
11 ψ(r; c ) ψ R = r ψ = ψ = TVD ψ = 0 r = 0 r = r = r Fig. : Sweby s TVD region for the flux limiter method based on the Godunov upwind and the Lax-Wendroff methods. / 56
12 Limiter functions Popular choices for flux limiters are the SUPERBEE limiter 0 if r 0, ψ sb (r) = r if 0 r, if r, r if r, if r, () the VANLEER limiter ψ vl (r) = 0 if r 0, r + r if r 0, (3) / 56
13 the VANALBADA limiter ψ va (r) = 0 if r 0, r( + r) + r if r 0, (4) and the MINBEE (or MINMOD) limiter 0 if r 0, ψ mb (r) = r if 0 r, if r. (5) 3 / 56
14 The WAF method The WAF method (Weighted Average Flux) is a conservative method qi n+ = qi n t ( ) f x i+ f i, (6) with numerical flux f waf = x f(q i+ x i+ (x, t))dx. (7) x See Toro (989) []. The integration range in (7) goes from the middle of cell I i (left of the interface) to the middle of cell I i+ (right of the interface) The integrand is the physical flux function f(q) in the conservation law q i+ (x, t) is the solution of the classical Riemann problem with (left) and qi+ n ) (right) data (qi n The flux formula can be generalised in a number of ways, to include for example, irregular meshes 4 / 56
15 WAF applied to the linear advection equation For the linear advection equation the classical Riemann problem is PDE: t q(x, t) + x f(q(x, t)) = 0, f(q(x, t)) = λq(x, t), IC: q(x, 0) = The solution is q i+ (x, t) = { q n i if x < 0, qi+ n if x > 0. q n i if x t < λ, qi+ n x if t > λ. Simple calculations show that the WAF flux is (8) (9) f waf i+ = ( + c)(λqn i ) + ( c)(λqn i+). (30) See Fig / 56
16 x/t = a β x β x t A B C t x 0 Fig. 3: The waf flux for the linear advection equation, for λ > 0. The WAF flux is identical to the Lax-Wendroff flux. This is only true for the linear advection equation. For non-linear equations WAF is a Riemann-problem based generalisation of the Lax-Wendroff method The flux is a weighted average of the upwind weight ( + c) and the downwind weight ( c) The upwind weight controls monotonicity and the downwind weight controls accuracy A TVD version will reduce the downwind contribution, when necessary, in favour of the upwind weight x 6 / 56
17 TVD version of WAF First note that the WAF flux can be written as f waf i+ = (λqn i + λq n i+) c(qn i+ q n i ). (3) The TVD version is f waf i+ = (λqn i + λq n i+) sign (c)ψwaflim(λq i+ i+ n λqi n ). (3) Here ψ waflim i+ = ( c )ψ fluxlim i+ (33) where ψ fluxlim is a flux limiter function. This version of the scheme i+ circumvents Godunov s theorem: the resulting scheme is non-linear. For complete details on the WAF method and the construction of its TVD version see Chapter 3 of Toro (009) []. 7 / 56
18 Reconstruction-based Methods Recall that cell averages are q n i x xi+ within the cell [x i, x i+ ], at time t = t n. x i q(x, t n )dx, (34) Given a set of cell averages {qi n } at time level n we construct polynomials p i (x), see Fig. 4, such that they satisfy: The conservation property. One expects the integral averages of the reconstructed functions to coincide with the original averages, that is q n i = x xi+ x i p i (x)dx. (35) The non-oscillatory property. To this end one imposes a TVD condition or an Essentially Non-Oscillatory (ENO) property, for example. This part of the scheme circumvents Godunov s theorem by constructing a non-linear scheme. 8 / 56
19 p i (x) q n i q n i q n i+ x i x i x i+ Fig. 4. Reconstruction polynomial function p i (x) from cell averages {q n j }. We study the simplest case, first-degree polynomials of the form within the cell [x i, x i+ p i (x) = q n i + (x x i ) i (36) ], at time t = t n. Here x i = (x i + x i+ ) is the cell centre and i is the slope of p i (x), still to be determined. Fig. 5 depicts two polynomials of the form (36) and are determined by respectively choosing the slopes x i i = qn i qn i x, i i+ = qn i+ qn i x. (37) 9 / 56
20 TVD slopes i / i+/ p L (x) p R (x) x i x i / x i x i+/ x i+ Fig. 5. Reconstruction functions p L (x), p R (x) from q n i, qn i, qn i+. The minmod slope. This has the following form i = minmod( i, i+ ), (38) where the minmod function is defined as minmod(x, y) = x if x y and xy > 0, y if x > y and xy > 0, 0 if xy < 0. x (39) 0 / 56
21 The superbee slope. This has the following form where () and i = minmod( i+ i = maxmod( () i, () i ), (40), i ), () i = minmod( i+ x if x y and xy > 0, maxmod(x, y) = y if x y and xy > 0, 0 if xy < 0. The MC slope. This has the following form where is the centred slope., i ) (4) (4) i = minmod( (c) i, i, i+ ), (43) (c) i = ( i + i+ ). (44) / 56
22 ENO reconstruction ENO (Essentially Non-Oscillatory) approach due to Harten et al. (987) [5] is used for reconstructions that satisfy the conservation and the non-oscillatory properties. In ENO the choice of slope i is made by analysing two potential candidate stencils S L and S R and their corresponding polynomials p L (x) and p R (x), namely Left stencil: S L = {i, i} p L (x) = a L + b L x, Right stencil: S R = {i, i + } p R (x) = a R + b R x. The coefficients a L, b L, a R and b R are to be found. See Fig. 5. Each candidate polynomial for each stencil must satisfy the conservation property, for each cell in the stencil. For the left stencil S L the conservation property requires x xi x i 3 p L (x)dx = q n i, x xi+ } (45) x i p L (x)dx = q n i. (46) / 56
23 On exact integration one obtains i i = qn i qn i. (47) x Similarly, by imposing conservation on the polynomial p R (x) based on the right stencil S R we have x xi+ x i p R (x)dx = q n i, After exact integration we obtain xi+ 3 p R (x)dx = qi+ n. (48) x x i+ i i+ = qn i+ qn i. (49) x Then the resulting two possible choices for the reconstruction polynomial are ( q p L (x) = q n n i + (x x i ) i q n ) i, x p i (x) = ( q p R (x) = qi n n + (x x i ) i+ qi n ) (50). x 3 / 56
24 The ENO method selects the polynomial with the smallest slope in absolute value, namely ENO i = i i+ Thus the resulting ENO polynomial is p L (x) = qi n + (x x i ) p ENO i (x) = if i i+, if i > i+. ( q n i qi n ) x ( q p R (x) = qi n n + (x x i ) i+ qi n ) x (5) if i i+ ; Compare the ENO slope (5) with the minmod TVD slope (38). if i > i+ (5) 4 / 56
25 Remarks on reconstruction schemes A modification of the ENO approach, called WENO (for weighted ENO), uses all candidate stencils in the ENO approach to produce the WENO polynomial, which is a linear combination of all candidate ENO polynomials. See the work of Jiang and Shu (996) [6]. An ENO reconstruction is not TVD an the resulting ENO numerical scheme admits spurious oscillations with new extrema. However, on mesh refinement the amplitude of the spurious oscillations tends to zero as the mesh width tends to zero. This property is not satisfied by linear schemes. 3 The ENO and WENO reconstruction approaches can be extended to any order of accuracy in one, two and three space dimensions. See [3], [4] on recent advances on reconstruction procedures in multiple space dimensions. 5 / 56
26 The MUSCL-Hancock scheme To compute the numerical flux f i+ one performs the following four steps: Step (I): Reconstruction. Use polynomial (36). See Fig. 6. q(x, t n ) q R i p i (x i+/ ) p i (x i ) p i (x i / ) q L i q L i q n i p i (x) q R i x i / x i x i+/ x Fig. 6. MUSCL-Hancock scheme. The reconstructed st-degree polynomial p i (x) has boundary extrapolated values qi L (left) and qi R (right). These are then evolved in time to become q L i and q R i. 6 / 56
27 Step (II): Cell-boundary values, or boundary extrapolated values qi L = p i (x i ) = qi n x i ; qi R = p i (x i+ ) = qi n + x i. (53) Step (III): Evolution of cell-boundary values for a time t q L i = qi L t x [f(qr i ) f(qi L )] = qi n ( + c) x i, (54) q R i = qi R t x [f(qr i ) f(qi L )] = qi n + ( c) x i. See Fig. 6. Step (IV): The Riemann problem. To compute f i+ one solves the classical Riemann problem. For the linear advection equation this is: PDE: t q + x (λq) = 0, q R i = qi n + ( c) x i if x < 0, IC: q(x, 0) = q L i+ = qn i+ ( + c) x i+ if x > 0. (55) 7 / 56
28 See Fig. 7 for al illustration of the Riemann problem with evolved initial conditions. The similarity solution q i+ (x/t) is q R i = qi n + ( c) x i if x/t < λ, q i+ (x/t) = q L i+ = qi+ n ( + c) x i+ if x/t > λ. and the flux is f i+ = f(q i+ (0)) = λq R i = λ ( qi n + ( c) x ) i λq L i+ = λ ( qi+ n ( + c) x ) i+ if λ > 0, if λ < 0. (56) 8 / 56
29 t Godunov state q i+/ (0) (b) x = 0 x q(x, 0) q R i (a) q L i+ x = 0 x Fig. 7. Step IV of the MUSCL-Hancock method. (a) shows the the piece-wise constant data for the classical Riemann problem obtained from the evolved extrapolated values. (b) illustrates the similarity solution of the Riemann problem and the Godunov state right at the interface, along the t axis. 9 / 56
30 Well-known, linear second-order methods can be reproduced, in particular: i = i i+ = qn i qn i x = qn i+ qn i x ( Warming-Beam ) λ > 0, ( Lax-Wendroff ) λ > 0, (57) c = qn i+ qn i x ( Fromm). These schemes are linear, of second-order accuracy and produce spurious oscillations with new extrema In order to avoid or control the spurious oscillations in a MUSCL-Hancock type scheme we must select the slopes in a nonlinear manner, as discussed previously 30 / 56
31 FORCE-based variant of MUSCL-Hancock Instead of solving the Riemann problem in step (IV), see Fig. 7, for the MUSCL-Hancock scheme and computing the Godunov flux, one can use a simple, monotone numerical flux, such as the force flux: f force i+ = f force (r, s), r = q R i+ i, s = q L i+. (58) For the linear advection equation the force flux is given as f force (r, s) = i+ ( + c) (λr) 4c ( c) (λs). (59) 4c This approach is applicable to general non-linear systems and the resulting schemes are simple and effective, though not as accurate as Riemann-problem based schemes, particularly for certain types or problems involving slowly moving intermediate waves. 3 / 56
32 Recalling the finite volume framework Consider the general balance law t q + x f(q) = s(q). (60) Integrate (60) on finite volume [x i, x i+ ] [t n, t n+ ] to obtain = qi n t where with q n+ i t n f i ] + ts i, (6) x [f i+ qi n = xi+ q(x, t n )dx, x x i f i+ = tn+ f(q(x t i+, t))dt, t n s i = tn+ xi+ s(q(x, t))dxdt, t x x i x = x i+ x i (6), t = t n+ t n, x i = (x i + x i+ ). (63) 3 / 56
33 Formula (6) is the starting point for constructing finite volume methods This is achieved by constructing approximations to the integrals in (6), resulting in the, approximate, finite volume formula: ˆq n+ i = ˆq i n t x [ ˆf i+ ˆf i ] + tŝ i. (64) This formula is used to evolve, in time, approximations to the cells averages in (6) In order to fully determine a finite volume formula (64) one requires: ˆfi+ : the numerical flux, an approximation to f i+ in (6) ŝ i : the numerical source, an approximation to s i in (6) Often one refers directly to (6) as the finite volume scheme, but (64) is the correct interpretation. The accuracy of the resulting numerical scheme depends crucially on the accuracy used to approximate the integrals in (6) 33 / 56
34 The ADER Approach ADER, first suggested by Toro et al. (00) [3], is an approach for constructing non-linear numerical schemes of arbitrary order of accuracy in both space and time. Here we restrict the presentation to second order as applied to the linear advection equation. ADER has the following stages: Step (A): Reconstruction. The simplest case is p i (x) = q n i + i(x x i ), p i+ (x) = q n i+ + i+(x x i+ ). (65) Step (B): Generalized Riemann problem. For the linear advection equation with linear source term, the simplest case, denoted by GRP, is PDE: t q + λ x { q = βq, pi (x) if x < 0, (66) IC: q(x, 0) = p i+ (x) if x > 0, whose solution is used to find the numerical flux ˆf i+ according to (6). Step C: Numerical source. One requires an approximation ŝ i to the volume integral in (6) to find the numerical source. 34 / 56
35 Step (B): generalised Riemann problem and the numerical flux We express the solution q i+ (τ) as a truncated power series expansion (τ) = q(0, 0 + ) + τ t q(0, 0 + ). (67) q i+ B: The leading term q(0, 0 + ). One solves the classical Riemann problem PDE: t q + λ x q = 0, IC: q(x, 0) = The similarity d (0) i+ d (0) (x/t) = i+ { pi (0) q n i + x i if x < 0, p i+ (0) q n i+ x i+ if x > 0, (x/t) of this problem is q n i + x i if x/t < λ, q n i+ x i+ if x/t > λ. (68) (69) 35 / 56
36 The leading term in (67) is called the Godunov state; it is the solution (69) along the interface x/t = 0, that is q q(0, 0 + ) = d (0) i n + x i if λ > 0, (0) = (70) i+ qi+ n x i+ if λ < 0. Remark: Upto this point we have the Kolgan scheme (Kolgan, 97) [7], whose linear version is however unstable. The reader is encouraged to derive the Kolgan scheme and analyse its linear stability. Note: Kolgan is credited to have been the first researcher to have constructed a non-linear scheme (97), a scheme that circumvents Godunov s theorem. See [8] for a English translation of the original Russian article. 36 / 56
37 B: The higher order term. To compute the second term in (67) we do the following: Time derivative in terms of the spatial derivative. We apply the Cauchy-Kowalewski procedure whereby the PDE in (68) is used to express the time derivative in terms of the space derivative and the source term, that is t q(x, t) = λ x q(x, t) + βq(x, t). (7) The sought term t q(x, t) will be known once the source function βq(x, t) and the spatial gradient x q(x, t) are known. These are required at the leading term, so that the expansion (67) reads q i+ (τ) = q(0, 0 + ) + [ λ x q(0, 0 + ) + βq(0, 0 + )]τ. (7) The gradient x q(0, 0 + ) remains unknown. 37 / 56
38 Evolution equation for the space derivative. t ( x q(x, t)) + λ x ( x q(x, t)) = β x q(x, t). (73) Riemann problem for the spatial derivative. We simplify (73) by neglecting the source term for the gradient. For the model equation this can be justified. Then we pose the classical homogeneous Riemann problem for the spatial derivative, namely PDE: t ( x q(x, t)) + λ x ( x q(x, t)) = 0, i if x < 0, (74) IC: x q(x, 0) = i+ if x > / 56
39 We solve the classical Riemann problem (74) i if x/t < λ, d () (x/t) = i+ i+ if x/t > λ. We need the Godunov state at x/t = 0 at the interface, which is i if λ > 0, x q(0, 0 + ) = d () (0) = i+ i+ if λ < 0 (75) (76) The sought solution (67) is q i+ (τ) = qi n + x i + τ [ λ i + β(qi n + x i) ], λ > 0, q n i+ x i+) + τ [ λ i+ + β(q n i+ x i+) ], λ < 0. (77) 39 / 56
40 The numerical flux in (6) is obtained according to (6), namely f i+ tn+ t t n f(q i+ (τ))dτ. (78) By performing exact integration one obtains λ ( qi n + ( c) x i + r(qn i + x i) ), λ > 0 f i+ = λ ( qi+ n ( + c) x i+) + r(qn i+ x i+) ), λ < 0 (79) where c = λ t/ x is the Courant number and r = tβ is a dimensionless reaction number. The ADER numerical flux (79) includes a contribution due to the source term. In other words, the flux knows of the source term, as one would expect. But note that there is still another contribution to the scheme resulting from the numerical source, still to be specified. 40 / 56
41 Step (C): numerical source The numerical source in (6) is obtained according to (6), namely s i = t x tn+ xi+ t n x i s(q i (x, t))dxdt. (80) q i (x, t) is an approximation in the control volume. To second-order, the mid-point integration rule suffices, and one requires q i (x i, t). Taylor expanding plus Cauchy-Kowalewski give q i (x i, τ) = q i (x i, 0 + ) + τ t q(x i, 0 + ) = q i (x i, 0 + ) + τ( λ x q i (x i, 0 + ) + βq i (x i, 0 + )). In the case considered we have q i (x i, 0 + ) = q n i = p i (x i ), x q i (x i, 0 + ) = i. } (8) The numerical source from (80) is ( s i = β qi n + ) t( λ i + βqi n ). (8) 4 / 56
42 The ADER scheme for linear advection-reaction The final form of the ADER scheme is qi n+ = qi n t x [f i+ f i ] + ts i, (83) with numerical flux λ ( qi n + ( c) x i + r(qn i + x i) ), λ > 0 f i+ = λ ( qi+ n ( + c) x i+ + r(qn i+ x i+) ), λ < 0 (84) and numerical source ( s i = β qi n + ) t( λ i + βqi n ). (85) In full, for the case λ > 0, the scheme reads qi n+ = qi n c [ (qi n qn i ) + ( c) x( i i ) ] cr [ 4 (qn i qn i ) + x( i i ) ] +r [ ( + r)qn i c x i] (86). 4 / 56
43 Remarks on the ADER approach Note that if the slopes are all set to zero and we also set to zero the contribution of the source to the numerical flux, then we produce a first-order numerical scheme for the model advection-reaction equation that consists of the first-order upwind method of Godunov for the advection part and an upwind discretisation of the source term, namely q n+ i = q n i c(q n i q n i ) + r [ cqn i + ( ] c)qn i. (87) This corresponds to one of the schemes analysed in (Roe, 986) [9] and (Bermudez and Vazquez, 994) [] in which the concept of upwinding the source term was employed By choosing the slopes i appropriately, we can reproduce the well-known schemes (57), when β = 0 (no source term) The second-order ADER scheme with ENO reconstruction is effectively a blend of the Warming-Beam and the Lax-Wendroff schemes (verify). 43 / 56
44 The ADER approach may look complex for constructing a second order method. However for higher order, the steps remain essentially the same The ADER approach is capable of producing non-linear schemes of arbitrary order of accuracy in both space and time Here we have presented the finite volume version. A discontinuous Galerkin finite element version is also available. See Dumbser et al. (008) [] An introduction to ADER schemes is found in chapters 9 and 0 of Toro (009) [] In the last 0 years there have been many developments of the ADER approach due to many researchers. See chapters 9 and 0 of Toro (009) [] for references. 44 / 56
45 Sample numerical results Godunov upwind Lax--Wendroff Position Position Limiter mina Limiter supera Position Position Fig. 8: LAE at T out =.Top left: Godunov. Top right: Lax-Wendroff Bottom: WAF method with limiters MINBEE (left) and SUPERBEE (right). 45 / 56
46 Sample numerical results.5 Exact MINBEE limiter, 00 cells.5 Exact SUPERBEE limiter, 50 cells Solution Solution X X.5 Exact SUPERBEE limiter, 00 cells.5 Exact ADER5, 50 cells Solution Solution X X Fig. 9: LAE at T out = 000π. Bottom right: ADER 5th order. Rest: Second-order TVD schemes. 46 / 56
47 Second-Order Schemes for Equation with Source Consider the model advection-reaction equation t q + x f(q) = s(q), f(q) = λq, s(q) = βq, (88) where λ > 0 and β 0 are two constants. Recall that for the ADER method, for λ > 0 the numerical flux is ( = λ qi n + ( c) x i + r(qn i + ) x i), (89) f i+ where c = λ t, r = tβ, (90) x are the CFL number and a reaction number. The numerical source is ( s i = β ( + r)qn i ) c x i (9) and is independent of the sing of λ, positive or negative. 47 / 56
48 MUSCL-Hancock for source term To account for the source term, the evolution step to compete the flux reads q i L = qi L t x [f(qr i ) f(qi L )] + t(βql i ), (9) q i R = qi R t x [f(qr i ) f(qi L )] + t(βqr i ). Compare this to (54). For λ > 0, simple manipulations give the numerical flux for the MUSCL-Hancock method accounting for the presence of the source as f i+ [ = λ qi n + ( c) x i + r(qn i + ] x i). (93) The numerical source for the MUSCL-Hancock method is taken, by construction, to be identical to that of the ADER scheme, namely (9). 48 / 56
49 The WAF method The Weighted Average Flux (WAF) method has numerical flux, for the homogeneous linear advection equation, given as f i+ = f i+ (r, s) = ( + φ i+ )(λr) + ( φ i+ )(λs), (94) where φ i+ is a WAF limiter. To apply WAF to the inhomogeneous case we perform the following steps Evolution of data: Evolve cell averages q n i, qn i+ to by solving the ODE q n i = q n i + tβqn i = ( + r)qn i, q n i+ = q n i+ + tβqn i+ = ( + r)qn i+, (95) d q = s(q), (96) dt by half a time step and initial condition q n i and q n i+. 49 / 56
50 Riemann problem and flux: Solve the classical Riemann problem with evolved data q i n, qn i+ to compute the usual WAF flux, which now reads f i+ = (+φ i+ )(+ r)(λqn i )+ ( φ i+ )(+ r)(λqn i+) (97) Numerical source: Compute the numerical source using the ADER procedure, namely [ s i = β ( + r)qn i ] c x i. (98) with WAF slope i = q i+ q i. (99) x Here the intercell states are taken as the WAF averaged states q i q i+ = ( + φ i )qi n + ( φ i )qi n, = ( + φ i+ )qi n + ( φ i+ )qi+ n. (00) 50 / 56
51 Test problem We solve the initial-value problem for the linear advection-reaction equation t q + λ x q = βq x [, ] q(x, 0) = h(x) = sin πx with periodic boundary conditions, whose exact solution is (0) q(x, t) = h(x λt)e βt. (0) M L -error L -order L -error L -order 0 3.8E xx E E E 3..33E E E 4.0 Convergence rates for LAE with source term. WAF method (non limited). Error measured in two norms, for a sequence of four meshes (M = 0, 0, 40, 80). 5 / 56
52 Exercises Problem : TVD schemes. Show that the following schemes for the linear advection equation (LAE) are TVD: Lax-Friedrichs, FORCE and Godunov upwind (for λ > 0 and for λ < 0). Problem : Monotone schemes are TVD. Show that a three-point, monotone scheme for LAE is TVD. Problem 3: Viscosity form of a scheme. Consider the three-point scheme written in viscous form q n+ i = q n i c(qn i+ q n i ) + d(qn i+ q n i + q n i ) (03) to solve LAE. By performing a truncation error analysis show that the coefficient of numerical viscosity in the sense of local truncation error is α visc = ( ) d c xλ. (04) c 5 / 56
53 Problem 4: Coefficient of viscosity form of a scheme. Show that by choosing appropriate values of d one can reproduce some well known schemes, that is Lax-Friedrichs, d(c) = ( + c ) FORCE, (05) c Godunov upwind, c Lax-Wendroff. Problem 5: Coefficient of viscosity for known schemes. Consider the functions d(c) above. Plot all the function d(c) above in the rectangle [, ] [0, ]. Identify the regions of monotone schemes and of non-monotone schemes. Discuss accuracy of the considered schemes in terms of the plot above. 53 / 56
54 Problem 6: WAF method. Consider the WAF method for LAE with numerical flux f i+ where φ i+ = f i+ (r, s) = ( + φ i+ )(λr) + ( φ i+ )(λs), (06) is a WAF limiter function. Applying the methodology used for deriving flux limiters ψ(r), find the TVD region for WAF and plot it. Show that ψ(r) is related to φ i+ (r) via where c is the Courant number. φ i+ (r) = ( c )ψ(r), (07) 54 / 56
55 Problem 7: The Kolgan scheme (97) Kolgan [7] is now recognized to have been the first to propose non-linear schemes. The linear version of Kolgan s scheme proceeds as for MUSCL-Hancock but neglecting the data evolution step. For LAE show that the numerical flux for the Kolgan scheme is λ(qi n + x i) if λ > 0, f i+ = (08) λ(qi+ n x i+) if λ < 0. For λ > 0 consider the linear version of Kolgan s scheme by taking linear (fixed) slopes i = i = qn i qn i x and write the resulting three-point scheme as q n+ i = b q n i + b 0 q n i + b q n i+, (09) explicitly identifying the coefficients b, b 0 and b. 55 / 56
56 Derive the local truncation error and state the order of accuracy of the scheme. Analyse the stability of the scheme using the (a) the von Neumann method, (b) Billett s method and (c) the viscosity form of the scheme. Repeat the above with slopes chosen as i = i+ = qn i+ qn i x Remark: non-linear Kolgan scheme. The linear version of Kolgan s scheme is unconditionally unstable, as shown above. However the non-linear version proposed originally by Kolgan uses a non-linear reconstruction whereby the choice of the slopes is non-linear. In particular Kolgan proposed the use of minimum slopes, which is equivalent to choosing the MINMOD slope limiter.. 56 / 56
57 L. Bermúdez and M. E. Vázquez. Upwind Methods for Hyperbolic Conservation Laws with Source Terms. Computers and Fluids, 3:049 07, 994. M. Dumbser, D. Balsara, E. F. Toro, and C. D. Munz. A Unified Framework for the Construction of One Step Finite Volume and Discontinuous Galerkin Schemes. J. Comput. Phys., 7: , 008. M. Dumbser and M. Käser. Arbitrary High Order Non-Oscillatory Finite Volume Schemes on Unstructured Meshes for Linear Hyperbolic Systems. J. Comput. Phys., ():693 73, 007. M. Dumbser, M. Käser, V. A. Titarev, and E. F. Toro. Quadrature-Free Non-Oscillatory Finite Volume Schemes on Unstructured Meshes for Nonlinear Hyperbolic Systems. J. Comput. Phys., 6(8):04 43, 007. A. Harten and S. Osher. 56 / 56
58 Uniformly High Order Accurate Nonoscillatory Schemes I. SIAM J. Numer. Anal., 4():79 309, 987. G. S. Jiang and C. W. Shu. Efficient Implementation of Weighted ENO Schemes. J. Comp. Phys., 6:0 8, 996. V. P. Kolgan. Application of the Principle of Minimum Derivatives to the Construction of Difference Schemes for Computing Discontinuous Solutions of Gas Dynamics (in Russian). Uch. Zap. TsaGI, Russia, 3(6):68 77, 97. V. P. Kolgan. Application of the principle of minimising the derivative to the construction of finite-difference schemes for computing discontinuous solutions of gas dynamics. J. Comp. Phys., 30: , 0. P. L. Roe. 56 / 56
59 Upwind Differencing Schemes for Hyperbolic Conservation Laws with Source Terms. In Proc. First International Conference on Hyperbolic Problems, Carasso, Raviart and Serre (Editors), pages 4 5. Springer, 986. P. K. Sweby. High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws. SIAM J. Numer. Anal., :995 0, 984. E. F. Toro. A Weighted Average Flux Method for Hyperbolic Conservation Laws. Proc. Roy. Soc. London, A43:40 48, 989. E. F. Toro. Riemann Solvers and Numerical Methods for Fluid Dynamics. A Practical Introduction. Third Edition. Springer Verlag, 009. E. F. Toro, R. C. Millington, and L. A. M. Nejad. Towards Very High Order Godunov Schemes. 56 / 56
60 In Godunov Methods: Theory and Applications. Edited Review, E. F. Toro (Editor), pages Kluwer Academic/Plenum Publishers, / 56
Basics 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 informationSimplified ADER schemes based on a time-reconstruction GRP solver
Simplified ADER schemes based on a time-reconstruction GRP solver Eleuterio F. TORO Professor Emeritus Laboratory of Applied Mathematics University of Trento, Italy Email: eleuterio.toro@unitn.it Webpage:
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 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 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 informationDepartment of Computing and Mathematics, Manchester Metropolitan University, Chester Street, Manchester, M1 5GD, U.K. s:
When citing this article please use: E F Toro, R C Millington and L A M Nejad. Towards very high-order Godunov schemes. In Godunov Methods: Theory and Applications. Edited Review. E F Toro (Editor). Kluwer
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 informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
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 informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
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 informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
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 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 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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
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 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 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 informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
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 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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
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 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 informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
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 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 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 informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationSpatial Discretization
Spatial Discretization Shengtai Li and Hui Li Los Alamos National Laboratory 1 Overview Finite Difference Method Finite Difference for Linear Advection Equation 3 Conservation Laws Modern Finite Difference
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationarxiv: v3 [math.na] 10 Mar 2015
A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws Michael Dumbser a, Olindo Zanotti a, Raphaël Loubère b, Steven Diot c a Department of
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationDiscontinuous Galerkin methods for Hyperbolic PDEs. Lecture 1
Discontinuous Galerkin methods for Hyperbolic PDEs. Lecture 1 Olindo Zanotti olindo.zanotti@unitn.it University of Trento Laboratory of Applied Mathematics PiTP 2016 Computational Plasma Astrophysics 18
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationHigh Order ADER FV/DG Numerical Methods for Hyperbolic Equations
Monographs of the School of Doctoral Studies in Environmental Engineering 10 High Order ADER FV/DG Numerical Methods for Hyperbolic Equations Cristóbal E. Castro 2007 Based on the Doctoral Thesis in Environmental
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
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 informationHIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 HIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
More informationDefinition and Construction of Entropy Satisfying Multiresolution Analysis (MRA)
Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 6 Definition and Construction of Entropy Satisfying Multiresolution Analysis (MRA Ju Y. Yi Utah State University
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationMAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS
MATHEMATICS OF COMPUTATION S 005-5780054-7 Article electronically published on May 6, 0 MAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS ORHAN MEHMETOGLU AND BOJAN POPOV Abstract.
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationCompact Central WENO Schemes for Multidimensional Conservation Laws
arxiv:math/9989v [math.na] 3 Nov 999 Compact Central WENO Schemes for Multidimensional Conservation Laws Doron Levy Gabriella Puppo Giovanni Russo Abstract We present a new third-order central scheme for
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 informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
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 informationThe multi-stage centred-scheme approach applied to a drift-flux two-phase flow model
The multi-stage centred-scheme approach applied to a drift-flux two-phase flow model Svend Tollak Munkejord 1 Steinar Evje 2 Tore Flåtten 2 1 Norwegian University of Science and Technology (NTNU), Department
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationLecture 6) Multidimensional Schemes for Hyperbolic Conservation Laws: Bringing it Together
Lecture 6) Multidimensional Schemes for Hyperbolic Conservation Laws: Bringing it Together By Prof. Dinshaw S. Balsara (dbalsara@nd.edu) Les Houches Summer School in Computational Astrophysics http://www.nd.edu/~dbalsara/numerical-pde-course
More informationSecond Order Positive Schemes by means of Flux Limiters for the Advection Equation
Second Order Positive Schemes by means of Flux Limiters for the Advection Equation Riccardo Fazio and Alessandra Jannelli Abstract In this paper, we study first and second order positive numerical methods
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No., pp. 57 7 c 00 Society for Industrial and Applied Mathematics ON THE ARTIFICIAL COMPRESSION METHOD FOR SECOND-ORDER NONOSCILLATORY CENTRAL DIFFERENCE SCHEMES FOR SYSTEMS
More informationOrdinary Differential Equations
Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation
More informationOn Convergence of Minmod-Type Schemes
On Convergence of Minmod-Type Schemes Sergei Konyagin, Boan Popov, and Ognian Trifonov April 16, 2004 Abstract A class of non-oscillatory numerical methods for solving nonlinear scalar conservation laws
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationComparison of (Some) Algorithms for Edge Gyrokinetics
Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
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 informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods 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 informationRunge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu 4 Abstract In this paper, we propose a new type of weighted
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 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 informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
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 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 informationRECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes
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 informationNon-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics
Non-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics Alexander Kurganov and Anthony Polizzi Abstract We develop non-oscillatory central schemes for a traffic flow
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationA unified flow theory for viscous fluids
Laboratoire Jacques-Louis Lions, Paris 27/10/2015 A unified flow theory for viscous fluids ILYA PESHKOV CHLOE, University of Pau, France joint work with EVGENIY ROMENSKI Sobolev Institute of Mathematics,
More informationA genuinely high order total variation diminishing scheme for one-dimensional. scalar conservation laws 1. Abstract
A genuinely high order total variation diminishing scheme for one-dimensional scalar conservation laws Xiangxiong Zhang and Chi-Wang Shu 3 Abstract It is well known that finite difference or finite volume
More informationTHIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS
Palestine Journal of Mathematics Vol. 5(Special Issue: ) (6), 8 37 Palestine Polytechnic University-PPU 6 THIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS A.
More informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
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 informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More information4.30 Thursday, Session 3 (afternoon): Godunov and ENO schemes
160 HYP 00 4.30 Thursday, Session 3 (afternoon): Godunov and ENO schemes The generalised Riemann problem: the basis of ADER schemes E. F. Toro Λ and V. A. Titarev Department of Civil and Environmental
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationA new type of Hermite WENO schemes for Hamilton-Jacobi equations
A new type of Hermite WENO schemes for Hamilton-Jacobi equations Jun Zhu and Jianxian Qiu Abstract In this paper, a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes is designed
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
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 information