Extremum-Preserving Limiters for MUSCL and PPM
|
|
- Ashley Ann Richards
- 5 years ago
- Views:
Transcription
1 arxiv: v [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 Phillip Colella Applied Numerical Algorithms Group, Lawrence Berkeley National Laboratory, Cyclotron Road, Berkeley, CA 9470, USA 5 March 009 Limiters are nonlinear hybridization techniques that are used to preserve positivity and monotonicity when numerically solving hyperbolic conservation laws. Unfortunately, the original methods suffer from the truncation-error being st order accurate at all extrema despite the accuracy of the higher-order method [,, 3, 4]. To remedy this problem, higherorder extensions were proposed that relied on elaborate analytic and geometric constructions [5, 6, 7, 8]. Since extremum-preserving limiters are applied only at extrema, additional computational cost is negligible. Therefore, extremum-preserving limiters ensure higher-order spatial accuracy while maintaining simplicity. This report presents higher-order limiting for (i) computing van Leer slopes and (ii) adjusting parabolic profiles. This limiting preserves monotonicity and accuracy at smooth extrema, maintains stability in the presence of discontinuities and under-resolved gradients, and is based on constraining the interpolated values at extrema (and only at extrema) by using nonlinear combinations of nd derivatives. The van Leer limiting can be done separately and implemented in MUSCL (Monotone Upstreamcentered Schemes for Conservation Laws) [] or done in concert with the parabolic profile limiting and implemented in PPM (Piecewise Parabolic Method) [9, 0]. The extremumpreserving limiters elegantly fit into any algorithm which uses conventional limiting techniques. Limiters are outlined for scalar advection and nonlinear systems of conservation laws. This report also discusses the 4 th order correction to the point-valued, cell-centered initial conditions that is necessary for implementing higher-order limiting. The material herein complements Colella and Sekora []. Lastly, there is no guarantee that extremumpreserving limiters preserve positivity. To ensure this property, one should combine the limiting with FCT (Flux-Corrected Transport) [3].
2 M Sekora, P Colella Algorithms for Scalar Advection Consider the following scalar equation in one spatial dimension: a t + λ a x = 0, λ = constant. () At time-step n, the average-valued, cell-centered quantity a over a finite volume of length h = x is: a n i h (i+ )h (i )h a(x, n t)dx. () MUSCL/PPM are conservative finite volume methods that are used to compute a n+ a n+ i = a n i λ t x ( a n+ i+ i : ) a n+, (3) i where a n+ is the average of a linear/parabolic interpolant over the interval swept out by i+ the characteristics crossing the cell face at (i + )h and is given by: a n+ i+ = a + i+ a i+ = σh = σh (i+ )h (i+ σ)h ai i (x)dx λ 0 (i+, (4) +σ)h a I (i+ )h i+ (x)dx λ 0 where σ = λ t [0, ] is the CFL number and a I x i (x) is the linear/parabolic interpolant, such that x [(i )h, (i + )h]. There are three variations of Godunov-type methods for which limiters can be implemented:. van Leer Limiter in MUSCL. van Leer Limiter + Parabolic Profile Limiter in PPM 3. Parabolic Profile Limiter in PPM Each of these algorithm variations are discussed below.. MUSCL. van Leer limit the differences ( a i ), giving nd order results. Apply the corresponding boundary conditions.. Use the van Leer limited differences to compute 4 th order differences: 4 a i = ( (a i+ 3 4 a i+) (a i + ) 4 a i ). (5)
3 Extremum-Preserving Limiters for MUSCL and PPM 3 3. Employ piecewise linear reconstruction by computing spatially extrapolated face-centered values at the low and high (left and right) edges of cells: a i,± = a i + (± σ) 4a i. (6). PPM. van Leer limit the differences ( a i ), giving nd order results. Apply the corresponding boundary conditions.. Employ either 4 th or 6 th order piecewise parabolic reconstruction: Method : use the van Leer limited differences and compute spatially extrapolated face-centered values at the low and high (left and right) edges of cells: 4 a i,+ = (a i+ + a i ) 6 ( a i+ a i ), (7) 4 a i, = (a i + a i ) 6 ( a i a i ), (8) 6 a i,+ = 4 a i,+ 30 (3( a i+ a i ) ( a i+ a i )), (9) 6 a i, = 4 a i, 30 (3( a i a i ) ( a i+ a i )), (0) α ± = a ± a i. () Method : employ either 4 th or 6 th order piecewise parabolic reconstruction without using the van Leer limited differences in Step : 4 a i+ 6 a i+ = 7 (a i+ + a i ) (a i+ + a i ), () = (a i+ + a i ) 8 60 (a i+ + a i ) + 60 (a i+3 + a i ), (3) α ± = a i± a i. (4) It is important to note that when a i = c a i (centered difference) the two formulations for the piecewise parabolic reconstruction are identical. 3. Limit the parabolic profile (α ± ). 4. Use the PPM predictor values to reconstruct the parabolic profile: a ± = a i + α ± + σ (±(α α + )) (α + α + )(3 σ). (5)
4 4 M Sekora, P Colella.3 Update Solution for MUSCL or PPM. Compute fluxes: F i+ = { λa + λ 0, λa λ 0.. (6). Use the divergence of the fluxes to update the solution: a n+ i.4 Conventional Limiters = a n i t x ( a n i = σ a n i σ ( F i+ a + i+ ) F i, (7) ) a + λ 0, i ). (8) λ 0. ( a a i+ i To complete the specification of the MUSCL and PPM schemes, one defines the conventional limiters used to constrain the interpolated profiles within each cell..4. Conventional van Leer Limiter Given a sequence of average-valued, cell-centered quantities, the conventional van Leer limiter proceeds with the following steps [, 9, 0]:. Compute one-sided and centered differences: a i = a i a i, (9) c a i = (a i+ a i ), (0) + a i = a i+ a i. (). Apply the conventional van Leer limiter: lim a i = min( a i, + a i ), () S = sign( c a i ), (3) { min( c a i, lim a i )S a i + a i > 0 a i = 0 a i + a i 0. (4) One significant defect of this method is the clipping of extremum when a i + a i 0. This clipping sets a i 0 as a precautionary measure for suppressing spurious oscillations.
5 Extremum-Preserving Limiters for MUSCL and PPM 5 How one arrives at the formula for the conventional van Leer limiter can be understood by considering the following example. Assume that one is not at an extremum, da dx ih < 0, and a i > + a i. The value of a((i + )h) at face-center (i + )h is approximated by: a((i + )h) = a((i + )h) h da dx (i+)h + a((i + )h) = a(ih) + h da dx ih + Clearly, a((i + )h) a((i + )h) and a((i + )h) a(ih) + h da dx ih. Therefore: (5) h da dx ih (a((i + )h) a(ih)), (6) and one arrives at Eq.. However, at extrema da 0, which leads one to Eq. 4. dx.4. Conventional Parabolic Profile Limiter Given the higher-order reconstruction of a ± such that α ± = a ± a i, the conventional parabolic profile limiter proceeds with the following steps [9, 0]:. Adjust α ± according to the following cases:. Reconstruct a ± given the adjusted values for α ±. α ± 0 α + α 0, α + α α+ > 4α, (7) α α + α > 4α +. One significant defect of this method is the adjustments to α ± make the reconstruction a monotone profile as a precautionary measure for suppressing spurious oscillations. However, this constraint is more restrictive than is required to preserve monotonicity [3]. The adjustments to α ± are derived from the interpolation polynomial described in the original Piecewise Parabolic Method: a i = a + σ(δ ± a + a 6 ( σ)), (8) δ ± a = a + a = α + α, (9) ( a 6 = 6 a i ) (a + + a ) = 3(α + + α ). (30) where σ [0, ] is the CFL number, which also corresponds to a dimensionless length scale. This scale is associated with each grid cell such that the left side of a cell is designated σ = 0 and the right side of a cell is designated σ =. Differentiating a i with respect to σ gives: da i dσ = δ ±a + a 6 ( σ) = (α + α ) (α + + α )( σ). (3)
6 6 M Sekora, P Colella By evaluating da i /dσ at the left and right sides of the cell: da i dσ da i dσ σ 0 + σ = δ ± a + a 6 = (α + + α ), (3) = δ ± a a 6 = (α + + α ). (33) Maximize a i with respect to σ by setting the derivatives equal to zero and solving the resulting equations. One arrives at the following result: α + α σ 0, α α + σ. (34) Extremum-Preserving Limiters For smooth solutions away from extrema, the MUSCL scheme is nd order accurate for linear advection whereas the PPM is 3 rd order accurate for linear advection and 4 th order accurate in the limit of vanishing CFL number. However, the monotonicity constraints at extrema reduce the truncation error to O(h) even at smooth extrema. This reduction in the overall accuracy of the method also introduces a non-smooth component to the error. To eliminate this problem, one constructs a new limiting scheme at extrema. The defect of the standard approach to limiting is most easily seen in the MUSCL limiter. Away from extrema, the magnitude of the slope is computed as the minimum of three undivided differences: the centered difference and twice the one-sided differences. In smooth regions away from extrema, the centered difference and the one-sided differences all approximate h da such that the minimum is always defined by the centered difference. At dx discontinuities, one of the one-sided differences is typically much smaller than the other two differences. Therefore, this one-sided difference is chosen because it leads to a reduction in the slope and suppresses oscillations. However, the idea behind this method fails at extrema because the derivative vanishes. Furthermore, the one-sided differences have opposite signs and this non-constant multiple bounds the centered difference. In the original van Leer and PPM limiters, the solution is to simply drop the order of the method to st order. In the approach used in this report, one changes the limiters at extrema, and only at extrema, by using comparisons of different estimates of the nd derivatives as a basis for whether to limit the underlying linear scheme. If the solution is smooth at the extremum, then all of the estimates of the nd derivative are comparable and the limiter leaves the underlying linear scheme unchanged. Discontinuities, underresolved gradients, and high-wavenumber oscillations are detected either by one of the estimates of the nd derivative being much
7 Extremum-Preserving Limiters for MUSCL and PPM 7 smaller than the others or by the various estimates of the nd derivatives changing sign. Either effect triggers a nontrivial limiting of the interpolating function in the cell and a resulting suppression of oscillations.. Extremum-Preserving van Leer Limiter The extremum-preserving van Leer limiter parallels the conventional van Leer limiter:. Compute one-sided and centered differences as well as one-sided differences that are an additional spatial step-size away: a i = a i a i, (35) a i = a i a i, (36) c a i = (a i+ a i ), (37) + a i = a i+ a i, (38) ++ a i = a i+ a i+. (39). Test for extrema. An extremum is defined if the following condition is satisfied: min( a i + a i, a i ++ a i ) < 0. (40) 3. If the above extremum condition is not satisfied, then limiting follows the conventional van Leer method: 4. If the above extremum condition is satisfied, then: lim a i = min( a i, + a i ), (4) Compute one-sided and centered nd derivatives: S = sign( c a i ), (4) a i = min( c a i, lim a i )S. (43) D a i = (a(ih) a((i )h) + a((i )h)), (44) h Dca i = (a((i + )h) a(ih) + a((i )h)), (45) h D+ a i = (a((i + )h) a((i + )h) + a(ih)). (46) h Find the minimum nd derivative over the five-cells in question: S = sign(d c a i), (47) D lima i = min( D ca i, max(s D a i, 0), max(s D +a i, 0)). (48)
8 8 M Sekora, P Colella Apply the modified van Leer limiter at the extremum: ( ) 3h min C V L lim a i = D lim a i, a i ( ) 3h min C V L D lim a i, + a i S c a i < 0 else (49) S = sign( c a i ), (50) a i = min( c a i, lim a i )S. (5) D designates a derivative while designates a difference. Derivatives and differences are similar operators that can be transformed back-and-forth when one considers the relevant Taylor expansion and multiplies/divides each operator by factors of h. C V L is a constant that is independent of the mesh spacing and as C V L 0, the extremum-preserving van Leer limiter reduces to the conventional van Leer limiter. For most calculations, C V L =.5. How one arrives at the tighter bound of 3h D lim a i for peak height in the extremum-preserving van Leer limiter can be understood by considering the following example. Assume that one is at a local maximum such that d a dx ih < 0 and da dx ih < 0. The van Leer limiting condition bounds the derivative on the high (right) edge of the cell: a i (a i+ a i ). (5) When looking for a bound on a i from the low (left) edge of the cell, an extremum near i implies that: 0 a((i )h) a((i )h) = h da dx (i 3 )h + O(h3 ). (53) Therefore, one arrives at the following bound on a i : a i = h da dx ih + O(h 3 ) = h da dx 3h d a (i 3 )h + dx + O(h3 ) 3h d a dx + O(h3 ). (54) These bounds are summarized in the following inequality: 0 a i C V L 3h d a dx, C V L >, (55) where the 3/ comes from the condition that (i 3 )h is the nearest face-centered point to the cell-centered point ih for which one can unequivocally assert that h da dx (i 3 )h O(h).. Extremum-Preserving Parabolic Profile Limiter If the piecewise parabolic reconstruction is done without using van Leer limited differences, then one has to include an additional step that limits extrema at cell faces such that a i+ lies between adjacent cell averages. Van Leer limiting automatically enforces this constraint.
9 Extremum-Preserving Limiters for MUSCL and PPM 9 Test for extremum at cell faces: (a i+ a i )(a i+ a i+ ) < 0. (56) If no extremum is found, then proceed without any adjustment to a i+. If an extremum is found, compute one-sided and centered nd derivatives: D a i+ D ca i+ D + a i+ = (a((i + )h) a(ih) + a((i )h)), (57) h = 3 ( a((i + )h) a((i + ) h )h) + a(ih), (58) = (a((i + )h) a((i + )h) + a(ih)). (59) h Find the minimum difference with respect to a i+ that is greater than zero: S i+ D lim a i+ = sign(dca i+ ), (60) = max(min(c PPM S D i+ a i+, S D i+ c a i+, C PPM S D i+ + a i+ 0). (6) ), Adjust a i+ according to: a i+ = (a i+ + a i ) 6 D lim a i+. (6) The above formulas are arrived at by considering the centered nd derivative: D c a i+ = (h/) (a((i + )h) a((i + )h) + a(ih) ). (63) The spatial resolution of h/ enters because one is considering differences between a(ih), a((i+ )h), a((i + )h). By substituting in the following expressions: a i a i+ = a(ih) h 4 D c a i + O(h 4 ), (64) = a((i + )h) h 4 D c a i+ + O(h 4 ), (65) which convert point-valued to average-valued quantities, one arrives at: h 4 D c a i+ = a i+ a((i + )h) + a i h D c a i+ + O(h 4 ). (66) Rearranging the above equation gives the expression: Dc a i+ = 3 (a h i+ a((i + ) )h) + a i. (67)
10 0 M Sekora, P Colella After finding the minimum nd derivative Dlim a i+, one can again use the above equation to assign a value to a i+ : Lastly, a ± = a i± a i+ = (a i+ + a i ) 6 D lim a i+. (68) and one proceeds with limiting α ±. Now given the higher-order reconstruction of a ± such that α ± = a ± a i, adjust α ± according to the following cases:. α + α 0 (ai+ a i )(a i a i ) 0 Compute the nd derivative with respect to σ from the the interpolation polynomial described in the original Piecewise Parabolic Method: D PPM a i = d a i d σ = a 6 = 6(α + + α ). (69) It is important to note that this nd derivative with respect to σ is also a nd order difference for a i. Furthermore, this difference represents the maximum difference that can occur across a cell given a parabolic profile. Compute one-sided and centered nd order derivatives: D a i = (a(ih) a((i )h) + a((i )h)), (70) h Dca i = (a((i + )h) a(ih) + a((i )h)), (7) h D+a i = (a((i + )h) a((i + )h) + a(ih)). (7) h Find the minimum difference that is greater than zero, given the nd derivatives over the five-cells in question as well as the above difference that was derived from the interpolation polynomial: S PPM = sign(d PPMa i ), (73) D lim a i = max(min(s PPM D PPM a i, C PPM S PPM D a i, Adjust α ± according to D lim a i and D PPM a i: C PPM S PPM D c a i, C PPM S PPM D + a i), 0). (74) α ± α ± D lim a i D PPM a i. (75). α + > 4α
11 Extremum-Preserving Limiters for MUSCL and PPM Compute the maximum value of α over a given cell: α max = α + 4(α + + α ). (76) This formula is arrived at by averaging the interpolation polynomial from the original Piecewise Parabolic Method on the left side of cell i: σ f = a(ξ)ξ (77) σ 0 = a + σ ( ( δ ± a + σ )) 3 a 6 (78) = a + σ (α + α ) 3σ ( σ ) (α + + α ). (79) 3 Maximize f with respect to σ: Therefore: df dσ = 0 σ max = α + + α (α + + α ). (80) α max = f (σ max ) a i = α + 4(α + + α ). (8) Preserve monotonicity by ensuring that the following condition is just satisfied: Solve the above equation for α + : I a i a i I = a i a i = α max. (8) S = sign(α ), (83) α + = I S ( I α I ) /. (84) 3. α > 4α + Compute the maximum value of α + over a given cell: α max + = α 4(α + + α ). (85) This formula is arrived at by averaging the interpolation polynomial from the original Piecewise Parabolic Method on the right side of cell i: f + = a(ξ)ξ (86) σ σ = a + σ ( ( δ ± a σ )) 3 a 6 (87) = a + σ (α + α ) 3σ ( σ ) (α + + α ). (88) 3
12 M Sekora, P Colella Maximize f + with respect to σ: Therefore: df + dσ = 0 σ max = α + + α (α + + α ). (89) α max + = f +(σ max ) a i = α 4(α + + α ). (90) Preserve monotonicity by ensuring that the following condition is just satisfied: Solve the above equation for α : 4. Reconstruct a ± given the adjusted values for α ±. I + a i+ a i I + = a i+ a i = α max +. (9) S + = sign(α + ), (9) α = I + S + ( I + α +I + ) /. (93) 3 Numerical Results Results are presented for D scalar advection. The standard test problems were employed to demonstrate improvements in accuracy [4]. The following parameters were used in all test problems: σ = 0., x [0, ], t = 0, λ =. (94) Periodic boundary conditions were used with the following number of ghost cells: 4 th Order Reconstruction 3 ghost cells per side (95) 6 th Order Reconstruction 4 ghost cells per side. (96) The test problems were defined by the following initial conditions: Gaussian Wave: a = exp ( 56(x 0.5) ) { (97) Semi-Circle Wave: a = (0.5 (x 0.5) ) 0.5 < x < else { (98) Square Wave: a = 0.5 < x < else (99) The conventional initialization for nd order accurate methods is obtained by approximating the average over a cell by the value at the center of the cell, i.e., the midpoint rule for
13 Extremum-Preserving Limiters for MUSCL and PPM 3 integrals. However, that initialization is only nd order accurate, and less accurate than what one expects for PPM when applied to the advection equation. For that reason, one uses a 4 th order accurate approximation to the cell average, following []: a i = a(ih) + h 4 D c a i, Dc a i = (a(i + )h) a(ih) + a((i )h)). (00) h This 4 th order correction is only administered when defining the initial conditions. The corresponding boundary conditions are applied and one begins advancing the temporal loop. For smooth solutions away from extrema, PPM is 3 rd order accurate for linear advection and 4 th order accurate in the limit of vanishing CFL number σ 0. The following definitions for the n-norms and convergence rates are used throughout this note. Given the numerical solution a r at resolution r and the analytic solution u, the error at a given point i is: The -norm and -norm of the error are: ǫ r i = a r i u i. (0) L = ǫ = i ǫ i x, L = ǫ = max i ǫ i. (0) The convergence rate is measured using Richardson extrapolation: R n = ln (L n(ǫ r )/L n (ǫ r+ )). (03) ln ( x r / x r+ )
14 4 M Sekora, P Colella 3. 4 th Order Reconstruction 4 th Order, Gaussian Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3 8.0E E E-.6.5E E E E E E E E-.5.6E E E E E E E E-.8.4E E-3.3.8E E E E E E-.8.5E E-3.3.8E E E-3 3.
15 Extremum-Preserving Limiters for MUSCL and PPM 5 4 th Order, Semi-Circle Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3.E- - 5.E E E E-3. 3.E E-3..5E E-3-4.E E E E-3..7E E-4..E E-3-4.E E E E-3..7E E-4..E E-3-4.E E E E-3..7E E-4..E- 0.3
16 6 M Sekora, P Colella 4 th Order, Square Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3.E E E E E E E E E- - 4.E E E E E E E E- - 4.E E E E E E E E E E E E E E E- 0.0
17 Extremum-Preserving Limiters for MUSCL and PPM th Order Reconstruction 6 th Order, Gaussian Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3 5.0E- -.6E E E E-3.7.8E E-4.9.5E E E E-.5.3E E E E-4.5.7E E- -.9E E E E-3.5.8E E-4.9.5E E- - 3.E E-.9.E E-3.6.8E E-3.9.6E-3.9
18 8 M Sekora, P Colella 6 th Order, Semi-Circle Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3 8.4E-3-4.E E E E-3..7E E-4..E E-3-3.6E E E E-3..4E E-4..9E E-3-3.7E E E E-3..4E E-4..9E E-3-3.6E E-3.3.9E E-3.3.3E E-4..9E- 0.3
19 Extremum-Preserving Limiters for MUSCL and PPM 9 6 th Order, Square Wave No Limiter PPM Limiting C PPM = 0 PPM Limiting C PPM =.5 VL+PPM Limiting C V L, C PPM =.5 N cell L R L R 3 9.8E- - 4.E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E- -0.
20 0 M Sekora, P Colella 3.3 Comparing Conventional and Extremum-Preserving Limiters This section compares the conventional limiter and extremum-preserving limiter when used in PPM for the one-dimensional scalar advection tests. The conventional limiter uses nd order van Leer limited differences to compute spatially 4 th order accurate extrapolated face-centered values [9, 0]. The extremum-preserving limiter employs 6 th order accurate piecewise parabolic reconstruction without using the van Leer limited differences []. As expected, the extremum-preserving limiters significantly reduce the error for the Gaussian wave (G), reduce the error only slightly for the semi-circle wave (SC), and does not reduce the error for the square wave (S). This comparison is shown in Figures and. Conventional Limiter N cell G L R G L R SC L R S L R 3 7.6E E E-3-8.4E E-.5.7E E E E-3.8.0E-.4.9E-3..8E E E E-4..6E- 0.8 Extremum-Preserving Limiter (C =.5) N cell G L R G L R SC L R S L R 3 4.E- -.9E E-3-7.7E E E-.6 3.E E E-3.5.8E-.4.4E-3..6E E-4.9.5E E-4..5E Sensitivity to C V L, C PPM Using a Gaussian wave as a test problem, the sensitivity in the extremum-preserving limiter was analyzed by plotting ln( x) versus ln(l ) for multiple values of C V L = C PPM between.5 and 5. The extremum-preserving limiter is insensitive to the limiting constants C V L, C PPM, where L changed by 4% when comparing C V L, C PPM =.5 and 5.
21 Extremum-Preserving Limiters for MUSCL and PPM Analytic Ext Pres Conven a x Figure : Extremum-Preserving vs Conventional Limiter. Gaussian Pulse, N cell = 8, t = 0 periods, CFL = 0., λ = a 0.4 Analytic Ext Pres Conven x Figure : Extremum-Preserving vs Conventional Limiter. Square Wave, N cell = 8, t = 0 periods, CFL = 0., λ =.
22 M Sekora, P Colella 4 Algorithms for Nonlinear Systems of Hyperbolic Conservation Laws Higher-order limiting can be extended to nonlinear systems of conservation laws. In one spatial dimension, these systems having the following form: U t + F x = 0, U(x, t) RN, F = F(U), F : R N R N, (04) where the quasilinear form is: U t + A U x = 0, A = A(U) = F U. (05) A has N real eigenvalues corresponding to N linearly independent eigenvectors: Ar k = λ k r k, l k A = λ k l k, (06) r k = r k (U k i ), l k = l k (U k i ). (07) Like the limiting done for scalar advection, each of the variations (VL, PPM, VL+PPM) can be implemented for nonlinear systems of hyperbolic conservation laws. The limiters are either applied componentwise to the primitive variables or applied to one characteristic field at a time. The former approach was implemented in Chombo (adaptive mesh refinement infrastructure for solving a wide variety partial differential equations) [5] and Athena (magnetohydrodynamics code for astrophysical applications) [6]. Results for extremum-preserving limiters applied to nonlinear systems of hyperbolic conservation laws are shown in [6]. In particular, [6] presents the Shu-Osher shock entropy wave interaction [7], where the numerical solution obtained using an extremum-preserving limiter is comparable to WENO3 and WENO5 [8]. These results are replotted below in Figures 3 and 4 with the permission of the authors. Figures 5 and 6 show results from the Woodward-Colella ramp problem [9] that were computed by Chombo using the unsplit PPM method of Miller-Colella [0] and the extremum-preserving limiters of Colella-Sekora []. This calculation employed a 64 6 base grid with two levels of adaptive mesh refinement such that there is a factor of four between each level. Thus, there is an effective resolution of The two figures show material density for the entire domain and the double Mach region, respectively. From these plots, it is clear that extremum-preserving limiters are robust enough to handle multidimensional shocks. Acknowledgment MS acknowledges support from the DOE CSGF Program which is provided under grant DE-FG0-97ER5308.
23 Extremum-Preserving Limiters for MUSCL and PPM density 3 N cell = 800 Conven N cell = 00 Conven x Figure 3: Conventional Limiter for the Shu-Osher Shock Entropy Wave Interaction. Computed using Athena. 5 4 density 3 N cell = 800 Ext Pres N cell = 00 Ext Pres x Figure 4: Extremum-Preserving Limiter for the Shu-Osher Shock Entropy Wave Interaction. Computed using Athena.
24 4 M Sekora, P Colella Figure 5: Extremum-Preserving Limiter for the Woodward-Colella Ramp Problem. Material density for the entire domain. Computed using Chombo. Figure 6: Extremum-Preserving Limiter for the Woodward-Colella Ramp Problem. Material density for the double Mach region. Computed using Chombo.
25 Extremum-Preserving Limiters for MUSCL and PPM 5 References [] J P Boris and D L Book. Flux-corrected transport. III. Minimal-error FCT algorithms. Journal of Computational Physics, 0:397-43, 976. [] B van Leer. Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection. Journal of Computational Physics, 3:76-99, 977. [3] S T Zalesak. Fully multidimensional flux-corrected transport algorithms for fluids. Journal of Computational Physics, 3:335-36, 979. [4] A Harten. High resolution schemes for hyperbolic conservation laws. Journal of Computational Physics, 49: , 983. [5] A Harten, B Engquist, S Osher, and S R Chakravarthy. Uniformly high order accurate essentially non-oscillatory schemes, III. Journal of Computational Physics, 77:3-303, 987. [6] G S Jiang and C W Shu. Efficient implementation of weighted ENO schemes. Journal of Computational Physics, 6:0-8, 996. [7] H T Huynh. Accurate upwind methods for the Euler equations. SIAM Journal on Numerical Analysis, 3:565-69, 995. [8] W J Rider, J A Greenough, and J R Kamm. Accurate monotonicity- and extremapreserving methods through adaptive nonlinear hybridizations. Journal of Computational Physics, 5:87-848, 007. [9] P Colella and P R Woodward. The Piecewise Parabolic Method (PPM) for gasdynamical simulations. Journal of Computational Physics, 54:74-0, 984. [0] G H Miller and P Colella. A conservative three-dimensional eulerian method for coupled solid-fluid shock capturing. Journal of Computational Physics, 83:6-8, 00. [] P Colella and M D Sekora. A limiter for PPM that preserves accuracy at smooth extrema. Journal of Computational Physics, 7: , 008. [] M Barad and P Colella. A fourth-order accurate local refinement method for Poisson s equation. Journal of Computational Physics, 09:-8, 005. [3] P R Woodward. Piecewise-parabolic methods for astrophysical fluid dynamics. K. H. A. Winkler and M. L. Norman, editors. Astrophysical Radiation Hydrodynamics, pg 45-36, 986.
26 6 M Sekora, P Colella [4] S Zalesak. R Vichnevetsky and R Stepleman, editors. Advances in Computer Methods for Partial Differential Equations V, 984. [5] Applied Numerical Algorithms Group, Lawrence Berkeley National Laboratory. [6] J M Stone, T A Gardiner, P Teuben, J F Hawley, and J B Simon. Athena: a new code for astrophysical MHD. Astrophysical Journal Supplement Series, 78:37-77, 008 (pending publication). [7] C W Shu and S Osher. Efficient implementation of essentially non-oscillatory shockcapturing schemes, II. Journal of Computational Physics, 83:3-78, 989. [8] D S Balsara and C W Shu. Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. Journal of Computational Physics, 60:405-45, 000. [9] P Woodward and P Colella. The numerical simulation of two-dimensional fluid flow with strong shocks. Journal of Computational Physics, 54:5-73, 984.
Numerical 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 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 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 informationHigh-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory!
High-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory! Why Higher Order?! Locally-refined grids, mapped-multiblock grids smooth except
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 informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
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 informationPart 1: Numerical Modeling for Compressible Plasma Flows
Part 1: Numerical Modeling for Compressible Plasma Flows Dongwook Lee Applied Mathematics & Statistics University of California, Santa Cruz AMS 280C Seminar October 17, 2014 MIRA, BG/Q, Argonne National
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
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 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 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 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 informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
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 informationA Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes )
A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) Takahiro MIYOSHI and Kanya KUSANO 1) Hiroshima University, Higashi-Hiroshima 739-856, Japan 1) Nagoya University,
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 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 informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
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 informationAn Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations. Thesis by Cheng Chi
An Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations Thesis by Cheng Chi In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah
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 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 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 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 informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
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 high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More 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 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 informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA 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 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 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 informationTotal Variation Theory and Its Applications
Total Variation Theory and Its Applications 2nd UCC Annual Research Conference, Kingston, Jamaica Peter Ndajah University of the Commonwealth Caribbean, Kingston, Jamaica September 27, 2018 Peter Ndajah
More informationComputational Fluid Dynamics. PHY 688: Numerical Methods for (Astro)Physics
Computational Fluid Dynamics Hydrodynamics When we discussed PDEs, we focused so far on scalar PDEs Often we wish to study systems of PDEs. Here we'll look at the equations of hydrodynamics Nonlinear system
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 informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 informationANALYSIS AND IMPLEMENTATION OF THE GAS-KINETIC BGK SCHEME FOR COMPUTATIONAL GAS DYNAMICS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 25, 21±49 (1997) ANALYSIS AND IMPLEMENTATION OF THE GAS-KINETIC BGK SCHEME FOR COMPUTATIONAL GAS DYNAMICS CHONGAM KIM*, KUN XU, LUIGI MARTINELLI
More informationThe PLUTO code for astrophysical gasdynamics
Mem. S.A.It. Suppl. Vol. 13, 67 c SAIt 009 Memorie della Supplementi The PLUTO code for astrophysical gasdynamics A. Mignone 1, 1 Dipartimento di Fisica Generale Amedeo Avogadro, Università degli Studi
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 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 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 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 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 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 informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall
COMPUTATIONAL METHODS FOR COMPRESSIBLE FLOW PROBLEMS Rémi Abgrall Université Bordeaux I, Talence Cedex, France Keywords: computational methods, numerical schemes, 1- D problems, multidimensional problems,
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
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 informationThe Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1
The Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 94305
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 information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
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 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 informationVarious Hydro Solvers in FLASH3
The Center for Astrophysical Thermonuclear Flashes Various Hydro Solvers in FLASH3 Dongwook Lee FLASH3 Tutorial June 22-23, 2009 An Advanced Simulation and Computing (ASC) Academic Strategic Alliances
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 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 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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
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 informationGas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations
Gas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations Iizln Xu* and Antony ~arneson~ Department of Mechanical and Aerospace Engineering Princeton University Princeton, New Jersey
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
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 informationEssentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws
NASA/CR-97-0653 ICASE Report No. 97-65 Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws Chi-Wang Shu Brown University Institute for Computer
More informationA New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods
A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods Raymond J. Spiteri Steven J. Ruuth Technical Report CS-- May 6, Faculty of Computer Science 65 University Ave.,
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
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 informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
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 informationConjugate lter approach for shock capturing
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2003; 19:99 110 (DOI: 10.1002/cnm.573) Conjugate lter approach for shock capturing Yun Gu 1 1; 2; ; and G. W. Wei 1 Department
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
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 informationGas Dynamics Equations: Computation
Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/
More information