Numerische Mathematik
|
|
- Octavia Carter
- 5 years ago
- Views:
Transcription
1 Numer. Math. ( 88: Digital Object Identifier (DOI.7/s5 Numerische Mathematik A third-order semi-discrete genuinely multidimensional central scheme for hyperbolic conservation laws and related problems Alexander Kurganov, Guergana Petrova Department of Mathematics, University of Michigan, Ann Arbor, MI 489-9, USA; {kurganov,petrova}@math.lsa.umich.edu Received February 7, / Published online December 9, c Springer-Verlag Summary. We construct a new third-order semi-discrete genuinely multidimensional central scheme for systems of conservation laws and related convection-diffusion equations. This construction is based on a multidimensional extension of the idea, introduced in [7] the use of more precise information about the local speeds of propagation, and integration over nonuniform control volumes, which contain Riemann fans. As in the one-dimensional case, the small numerical dissipation, which is independent of O( t, allows us to pass to a limit as t. This results in a particularly simple genuinely multidimensional semi-discrete scheme. The high resolution of the proposed scheme is ensured by the new twodimensional piecewise quadratic non-oscillatory reconstruction. First, we introduce a less dissipative modification of the reconstruction, proposed in [9]. Then, we generalize it for the computation of the two-dimensional numerical fluxes. Our scheme enjoys the main advantage of the Godunov-type central schemes simplicity, namely it does not employ Riemann solvers and characteristic decomposition. This makes it a universal method, which can be easily implemented to a wide variety of problems. In this paper, the developed scheme is applied to the Euler equations of gas dynamics, a convectiondiffusion equation with strongly degenerate diffusion, the incompressible Euler and Navier-Stokes equations. These numerical experiments demonstrate the desired accuracy and high resolution of our scheme. Mathematics Subject Classification (99: 65M, 65M5 Correspondence to: A. Kurganov
2 684 A. Kurganov, G. Petrova Introduction We are interested in Godunov-type schemes for solving the multidimensional system of conservation laws (. u t + x f(u =, x IR d, and the closely related convection-diffusion equations (. subject to the initial data (.3 u t + x f(u = x Q(u, x u, x IR d, u(x, = u (x. Here Q satisfies the weak parabolicity condition s Q(u, s for all (u, s. Godunov-type schemes are projection-evolution methods, which at every time step consist of three consecutive stages reconstruction, evolution and projection. First, a piecewise polynomial interpolant is reconstructed from the piecewise constant data (the cell-averages, computed at the previous time step. Then, the interpolant is evolved in time according to (.. Finally, the solution is projected onto a space of piecewise constants. Depending on the projection step, we distinguish two kinds of Godunov-type schemes central and upwind. The Godunov-type central schemes are based on exact evolution and averaging over Riemann fans. Therefore, in contrast to the Godunov-type upwind schemes, no characteristic decomposition or Riemann solvers are involved. This makes the central schemes simple, and they can be implemented as a black box solver for the general multidimensional systems (. and (.. In the one-dimensional case, examples of such schemes are the first-order (staggered Lax-Friedrichs scheme [9, 7], the second-order Nessyahu- Tadmor scheme [3], and their higher-order generalizations [9, 4, ]. Second-order multidimensional central schemes were introduced in [, ], and their higher-order extensions were developed in [, 3]. The advantage of the higher-order schemes is that they reduce the excessive numerical viscosity, typical for the Lax-Friedrichs scheme, and give much sharper resolution of the shocks and rarefactions. We would also like to mention the central schemes for the Hamilton- Jacobi equations, recently developed in [5, 6, 8], the central schemes for incompressible flows in [ 4, 6, 4], and their applications to various systems, for example, [6,,33]. Unfortunately, the aforementioned staggered central schemes may not provide a satisfactory resolution when small time steps are enforced by stability restrictions (e.g., this may occur when one applies these schemes to
3 Genuinely multidimensional central scheme 685 the (degenerate parabolic systems (.. Another disadvantage of these schemes is that they cannot be used for steady-state computations, or equivalently, they do not admit semi-discrete form. The reason is in the accumulation of numerical dissipation, which is of order O( ( xr t, where r is the formal order of the scheme. These problems have been recently resolved in [7], where a new secondorder Godunov-type central scheme is introduced. The proposed scheme has smaller dissipation ( ( t r, and thus it can be efficiently used with time steps as small as required. In the new construction, the evolution step is performed by integrating over non-equal control volumes, whose sizes are proportional to the local speeds of propagation. The evolved solution is then projected back onto the original grid in a more complicated manner, which requires an additional piecewise polynomial reconstruction. In this way, a non-staggered fully-discrete central scheme is derived. It can be naturally reduced to a particularly simple semi-discrete form (for details see [7]. The same idea was used in [6] to develop a third-order semi-discrete central scheme. Both in [7] and [6], the multidimensional semi-discrete schemes were obtained by using a dimension-by-dimension approach the one-dimensional numerical flux was straightforwardly applied in all coordinate directions. In this paper, we present new high-order genuinely multidimensional semi-discrete central schemes. They are based on a multidimensional generalization of the one-dimensional construction, used in [7, 6], namely we first integrate over non-equal volumes to derive a fully-discrete scheme, and then we pass to a semi-discrete limit (we let t, keeping x fixed. Semi-discrete central schemes consist of three independent building blocks a piecewise polynomial reconstruction, a spatial flux discretization and an ODE solver. Our genuinely multidimensional recipe provides the numerical fluxes. For the other two blocks one can pick any non-oscillatory reconstruction and an efficient ODE solver. For second-order schemes one uses a piecewise linear reconstruction a library of such (essentially non-oscillatory reconstructions is available (see e.g. [8, 9, 5,, 3, 3]. To achieve a third-order accuracy, a piecewise quadratic approximation is needed. However, to build a third-order nonoscillatory interpolant is a highly nontrivial problem. In the one-dimensional case, one of the possibilities is to use the essentially non-oscillatory (ENO reconstruction [9, 35], or its weighted generalizations [8,,, 3]. Multidimensional ENO-type reconstructions were introduced in [, 3]. A disadvantage of this approach is that it employs smoothness indicators, which require a certain a-priori information about the solution. Alternative one-dimensional non-oscillatory piecewise quadratic reconstructions were proposed in [7, 9]. They both satisfy the number of extrema
4 686 A. Kurganov, G. Petrova diminishing (NED property and do not require the use of smoothness indicators. In this work, we present a less dissipative one-dimensional modification of the reconstruction from [9]. We also generalize the one-dimensional reconstruction for the computation of the two-dimensional numerical fluxes in our scheme. The paper is organized as follows. In Sect., we briefly describe the recent development of Godunov-type central schemes. Next, in Sect. 3., we present the one-dimensional modification of the non-oscillatory piecewise quadratic reconstruction from [9]. In Sects. 3. and 3.3, we develop our new genuinely multidimensional central scheme, which employs the multidimensional third-order reconstruction, introduced in Sect The results from our numerical experiments are shown in Sect. 4. We apply our thirdorder scheme to a variety of one- and two-dimensional problems, including the compressible Euler equations of gas dynamics, a convection-diffusion equation with a nonlinear degenerate diffusion, the incompressible Euler and Navier-Stokes equations. These numerical examples demonstrate the expected third-order accuracy and the high resolution of our scheme. Godunov-type central schemes an overview In this section, we discuss Godunov-type central schemes. In Sect.., we recall the construction of staggered central schemes. Then, in Sect.., we review the recently proposed [7, 6] semi-discrete central schemes. Finally, in Sect..3, we briefly describe the one-dimensional non-oscillatory piecewise parabolic reconstruction from [9]. For the sake of simplicity, we will consider only uniform grids. Throughout this paper, we will use the following notation. Let x j := j x, x j± (j ± / x, t n := n t, u n j := u(x j,t n, where x and t are small spatial and time scales, respectively. :=. Central schemes for conservation laws We begin with the definition of the sliding average of u(,t, ū(x, t := (. u(ξ,t dξ, I(x ={ξ : ξ x < x x }. I(x We then integrate (. over I(x [t, t + t], and arrive at the equivalent formulation of the system (., ū(x, t + t =ū(x, t x
5 Genuinely multidimensional central scheme 687 (. [ t+ t τ=t f(u(x + x,τ dτ t+ t τ=t ] f(u(x x,τ dτ. This is the starting point for the construction of Godunov-type schemes, which consists of three steps reconstruction,evolution and projection. Reconstruction. At time level t = t n a piecewise polynomial function ũ(x, t n, (.3 ũ(x, t n =p n j (x, x j <x<x j+, j, is reconstructed from the data ū n j := ū(x j,t n, computed at the previous time step. Here, {p n j } are algebraic polynomials, and therefore ũ(,tn may be discontinuous only at x j±. The reconstruction is required to satisfy two basic properties conservation of the given cell averges : x j+ x j ũ(ξ,t n dξ =ū n j, j; accuracy (for smooth u(x, t: ũ(x, t n =u(x, t n +O(( x r, where r is the (formal order of the scheme to be constructed. Evolution. The piecewise polynomial ũ(,t n is evolved exactly according to the integral equation (., and the solution at time t = t n+ is obtained in terms of its sliding averages. Projection. An evaluation of these sliding averages at particular grid points provides the approximate cell averages of the solution at the next time level, {ū n+ }. Choosing x = x j in (. results in an upwind Godunov-type scheme. In this case, the solution of the initial value problem (.,(.3 with initial data u (x =ũ(x, t n, prescribed at t = t n, may be non-smooth in the neighborhood of the points {x j+ }. Therefore, the evaluation of the flux integrals in (. requires the use of a computationally expensive (approximate Riemann solver and chraracteristic decomposition. An alternative staggered grid approach, namely the choice of x = x j+ in (., was proposed by Nessyahu and Tadmor [3]. In this case, (.
6 688 A. Kurganov, G. Petrova leads to the Godunov-type central scheme ū n+ j+ (.4 [ x j+ = p n j (x dx + x x j λ t [ t n+ t n x j+ x j+ f(u(x j+,t dt ] p n j+(x dx t n+ t n ] f(u(x j,t dt, λ := t x. In contrast to the upwind schemes, the RHS of (.4 can be easily evaluated. The spatial integrals can be computed exactly for a given reconstruction {p n j }. Also, due to the finite speed of propagation, the solution of (.,(.3 with u (x =ũ(x, t n is smooth in (x j ε, x j + ε [t n,t n+, under an appropriate CFL restriction on t. Hence, the flux integrals in (.4 can be approximated by a quadrature formula of order r. Notice that the function values, needed in the quadrature, may be computed either by Taylor expansion, or by a Runge-Kutta method ([9,4]. For example, if we use the second-order reconstruction p n j (x =ū n j + s n j (x x j, and the midpoint quadrature rule ū n+ j+ t n+ f(u(x j,t dt = tf(u(x j,t n+/ + O(( t, t n we arrive at the Nessyahu-Tadmor (NT staggered central scheme = (ūn j +ū n j+ [ ] 8 (sn j s n j+ λ f(u n+ j+ f(un+ j, where u n+ j are evaluated, for example, by the Taylor expansion u n+ j =ū n j t (f x n j. To ensure the non-oscillatory nature of the NT scheme, one should apply a nonlinear limiter while computing the slopes {s n j }. In the original NT scheme [3] the so-called minmod limiter [, 8, 3] is used, namely, the slopes are determined by (ūn s n j ū n j (.5 j = minmod x where minmod(a, b := sgn(a + sgn(b, ūn j+ ūn j x, min( a, b.
7 Genuinely multidimensional central scheme 689. Semi-discrete central schemes As mentioned in the Introduction, the staggered central schemes do not admit semi-discrete formulation. This makes them inappropriate for steady state computations or for solving the diffusion dominated convection-diffusion systems (.. The encountered difficulties were recently resolved by the non-staggered central schemes, proposed in [7]. Next, we give a brief overview of these schemes. We begin with a piecewise polynomial reconstruction {p n j }, which may be discontinuous at the interface points {x j+ }. This discontinuities propagate with different speeds, whose upper bounds are ( f a n := max j+ u (ω. Here, C(u n,u n+ j+ j+ p n j (x j+ and u n+ j+ Jacobian f ω C(u n j+,u n+ j+ ρ is the curve in the phase space that connects u n := j+ := p n j+ (x j+, and ρ is the spectral radius of the ( f u u. For example, in the genuinely nonlinear or linearly degenerate case, we have { ( f ( f a n = max ρ j+ u (un,ρ } j+ u (un+. j+ The main idea in the constrtuction of the schemes in [7] is based on utilizing these CFL related local speeds of propagation. We consider the domains (.6[x n j,r,xn j+,l] [tn,t n+ ] and [x n j+,l,xn j+,r] [tn,t n+ ], with x n j+,l := x j+ ta n, and x n j+ j+,r := x j+ + ta n, where j+ the solution of (.,(.3 with u (x =ũ(x, t n is smooth and nonsmooth, respectively. Given the reconstruction {p n j (x}, we integrate over these domains and obtain the cell averages w n+ j = x n j+,l xn j,r [ xn j+,l x n j,r p n j (x dx (.7 t n+ t n ] (f(u(x nj+,l,t f(u(xnj,r,t dt,
8 69 A. Kurganov, G. Petrova n+ u j n+ w j+ n+ w j n+ w j+/ n+ wj-/ n+ wj- n u j+ n u j n u j- n x x n n n n j-3/,r j-/,l xj-/,r xj+/,l xj+/,r n x j+3/,l x j- xj-/ x j xj+/ xj+ xj+3/ Fig... Modified central differencing x and (.8 w n+ j+ = x n j+,r xn j+,l t n+ t n [ x j+ x n j+,l p n j (x dx + x n j+,r x j+ ] (f(u(x nj+,r,t f(u(xnj+,l,t dt p n j+(x dx over the corresponding non-equal spatial cells (see Fig... The spatial integrals in (.7 and (.8 can be computed explicitly. For the approximation of the flux integrals an appropriate quadrature rule can be used, since the solution is smooth along the vertical lines x = x n j+,l t < t n+. Further, for the construction of a and x = x n for j+ tn,r fully-discrete scheme, we need to project the computed intermediate data { w n+ } onto the original grid. To keep the high order of accuracy, we use another non-oscillatory piecewise polynomial reconstruction, this time over the nonuniform grid (for details see [7,6]. We would like to emphasize that the resulting fully-discrete scheme is still Godunov-type central scheme, which does not employ any Riemann solver or characteristic decomposition. Moreover, it admits an extremely simple semi-discrete formulation, which can be achieved in the following way. We substitute ū n+ j from the fully-discrete scheme in
9 Genuinely multidimensional central scheme 69 d ū n+ j ū n j dtūj(t = lim, t t and take a limit as t, with x fixed. This was done in [7] and [6] for the second- and the third-order schemes, respectively. Actually, in a similar way, we derive the rth-order semi-discrete scheme d dtūj(t = H j+ (.9 with the numerical flux (t H j x (t, (. f(u + (t + f(u (t j+ H j+ (t := j+ a j+ (t[ ] u + (t u (t. j+ j+ The intermediate values u ± j+ are given by (. u + j+ := p j+ (x j+, u j+ := p j (x j+, where {p j } is the rth-order piecewise polynomial reconstruction at time t. We leave for the reader the details of the calculations. Remark. The semi-discrete scheme (.9 (. is a system of time dependent ODEs, which can be solved by an ODE solver of an appropriate order (see the discussion in [7]. In contrast to the staggered central schemes, described in the previous section, the semi-discrete central schemes can be efficiently applied to the convection-diffusion equations (.. The resulting scheme is (. d dtūj(t = H j+ (t H j (t + Q j+ x (t Q j x (t, where Q j+ is a reasonable approximation to the diffusion flux at x = x j+ (for example, one may use a standard central differencing, see [6]. This technique can also be directly applied to multidimensional problems if one uses the so-called dimension-by-dimension approach (see [7, 6]. For instance, in the two-dimensional case, the corresponding semi-discrete scheme for the system (.3 u t + f(u x + g(u y =
10 69 A. Kurganov, G. Petrova is (.4 d H x j+ dtūj,k(t =,k(t Hx j,k(t x H y (t H y (t j,k+ j,k. y Here, the numerical fluxes are a straightforward generalization of the onedimensional numerical flux (., (.5 f(u + H x j+,k(t := j+,k(t + f(u j+,k(t a x j+,k(t ] [u +j+,k(t u j+,k(t, and (.6 g(u + (t + g(u (t H y j,k+ (t := j,k+ j,k+ a y (t[ ] j,k+ u + (t u (t. j,k+ j,k+ For more details, we refer the reader to [7] and [6]. In this paper (Sect. 3., we present an alternative genuienly multidimensional extension of the scheme (.9 (., which is different from (.4 (.6 for third- and higher-order schemes. Our numerical experiments demonstrate the advantage of our genuienly multidimensional construction over the straightforward dimension-by-dimension approach..3 A piecewise quadratic non-oscillatory reconstruction In this section, we give a brief description of the third-order non-oscillatory reconstruction of Liu and Tadmor [9]. Later on, in Sect. 3., we propose a new, less dissipative modification of this reconstriction. Let us denote by D ± and D the one-sided and the central divided differences v(x ± x v(x D ± v(x :=±, x v(x + x v(x x D v(x :=. x
11 Genuinely multidimensional central scheme 693 We start with the basic piecewise quadratic reconstruction (.7 qj n (x =(ū n j ( x 4 D +D ū n j +D ū n j (x x j + D +D ū n j (x x j, which obeys the requirements of conservation and accuracy. Moreover, it satisfies a shape-preserving property, namely qj n(x is monotone on (x j ū n j+ is monotone, qj n(x has extremum in the interior of (x j local extremum.,x j+ if and only if the sequence ū n j, ūn j,,x j+ if and only if ū n j is a This property implies that in the process of reconstruction no new extrema are created at the interior of the intervals (x j,x j+. The only places where any new extrema may appear are the interface points {x j± }.Toavoid this, one should use new quadratic polynomials {p n j }, which are convex combination of the basic parabolas qj n and the cell averages ūn j, this is (.8 p n j (x =( θ n j ū n j + θ n j q n j (x, <θ n j <. Here, θj n are determined from the formula { M n ū n j+ min j m n ū n } j Mj n, j ūn j m n j,, if ū n ūn j < ūn j < ūn j+, j θj n := { M n ū n j min j m n ū n } j+ Mj n, j ūn j m n j,, if ū n ūn j > ūn j > ūn j+, j, otherwise, (.9 where and M n j± m n j± M n j { = max { m n j = min q n j (x j+ q n j (x j+ { = max ū ( n j +ū n j± { = min ū ( n j +ū n j±,q n j (x j,qj n (x j }, },,q n j±(x j±,q n j±(x j± }, }.
12 694 A. Kurganov, G. Petrova As it was proved in [9], such a choice of the limiter θj n ensures that in smooth regions θj n ( x3. Thus, the reconstruction {p n j } is a thirdorder accurate. The polynomials {p n j } also inherit the conservation and the shape-preserving properties of {qj n}. Moreover, the limiters θn j are designed to prevent oscillations at the interface points in the sense that new extrema are created neither at the interior of the intervals (x j,x j+, nor at their endpoints. We refer the reader to [9] for details. 3 New genuine multidimensional semi-discrete central scheme In this section, we introduce a new semi-discrete central scheme, which is a genuine multidimensional extension of the scheme (.9 (.. It employs a multidimensional generalization (Sect. 3.4 of the new piecewise quadratic reconstruction (Sect A new one-dimensional reconstruction As in Sect..3 (see also [9], we use the basic piecewise quadratic function {qj n}, defined in (.7. The new reconstruction pn j is obtained as a convex combination of qj n and the piecewise linear interpolant Ln j ( vs. the convex combination of qj n and the piecewise constant function ūn j in (.8, (3. p n j (x =( θj n L n j (x+θj n qj n (x, <θj n <, where (3. L n j (x =ū n j + s n j (x x j. The non-oscillatory property of the new piecewise quadratic reconstruction {p n j } is ensured by an appropriate choice of {θn j } and by a nonoscillatory property of {L n j } (which depends on the choice of the slopes {s n j }. We take {θn j } to be { M n L n j+ min j (x j+ mn L n Mj n Ln j (x j+, j j (x } j m n j Ln j (x j,, if ū n j < ūn j < ūn j+, θj n { := M n L n j min j (x j mn L n Mj n Ln j (x j, j+ j (x } (3.3 j+ m n j Ln j (x j+,, if ū n j > ūn j > ūn j+,, otherwise,
13 Genuinely multidimensional central scheme 695 where M n j { = max { m n j = min q n j (x j+ q n j (x j+,q n j (x j,qj n (x j }, }, and M n j± m n j± { = max ( { ( = min L n j (x j± L n j (x j± +L n j±(x j± +L n j±(x j±,qj±(x n j±,q n j±(x j± }, }. In this paper, we use the non-oscillatory TVD minmod interpolant {L n j }, defined in (3.,(.5. This guarantees that the number of local extrema of {p n j } does not exceed the number of extrema of {ūn j } (NED property. The proof of the NED property is the same as in [9], and we omit the details. One can also easily check that our new reconstruction (3. (3.3 possesses all the properties of the basic parabolas {qj n } conservation, accuracy and shape-preserving. Remark 3. The new reconstruction (3. (3.3 is a less dissipative generalization of (.7 (.9. This fact allows us to achieve a better resolution, as illustrated in the numerical examples in Sect. 4. At the same time, the modified reconstruction (3. (3.3 is only slightly more computationally expensive than the original one. Notice, that if one sets all slopes s n j to be zero, then (3. (3.3 is reduced to (.7 (.9. In summary, the resulting piecewise parabolic reconstruction has the form (3.4 p n j (x = where θ n j and sn j ( ū n j θj n ( x ( 4 D +D ū n j + θj n D ū n j +( θj n s n j (x x j + θn j D +D ū n j (x x j, are given by (3.3 and (.5, respectively. 3.A genuinely multidimensional approach fully-discrete set up In this section, we generalize the idea of the one-dimensional construction, described in Sect.. (see Fig.. for the multidimensional case. For simplicity, we will discuss only the two-dimensional system of hyperbolic conservation laws (.3.
14 696 A. Kurganov, G. Petrova Let x j := j x, y k := k y, x j± := x j ± x, y k± := y k ± y. Assume that at time t = t n we have already constructed a conservative piecewise polynomial interpolant of an appropriate order ũ n (x, y := j,k p n j,k (x, yχ j,k(x, y, with possible discontinuities along the interface lines x = x j± and y = y k±. Later on, in Sect. 3.4, we present a new piecewise quadratic reconstruction, especially designed for the purposes of this scheme. We denote the corresponding point values by (3.5 u j,k := p n j,k (x j,y k, u N j,k := pn j,k (x j,y k+, u S j,k := pn j,k (x j,y k, u E j,k := pn j,k (x j+,y k, u W j,k := pn j,k (x j,y k, u NE j,k := pn j,k (x j+,y k+, u NW j,k := pn j,k (x j u SW j,k := pn j,k (x j,y k+, u SE j,k := pn j,k (x j+,,y k,y k, and the cell average by ū j,k := x y x j+ y k+ y k p n j,k (x, y dxdy. x j As in the one-dimensional case, our construction employs the CFL related maximal local speeds of propagation of the discontinuities, and we denote these speeds by {a x j±,k} and {ay }. They are not easy to compute, but j,k± in practice one may use, for example, the values (3.6 ( f a {ρ x j+,k := max { ( g := max ρ a y j,k+ u (uw j+,k u (us j,k+ Further, we consider the nonuniform domains ( f },ρ u (ue j,k, ( g },ρ u (un j,k. D j,k [t n,t n+ ], D j±,k [tn,t n+ ], D j,k± [t n,t n+ ], D j±,k± [t n,t n+ ],
15 Genuinely multidimensional central scheme 697 y k+ y k+/ y k D jk y k-/ y k- x j- x x x x j-/ j j+/ j+ Fig. 3.. Two-dimensional modified central differencing outlined in Fig. 3.. They are determined by the corresponding local speeds, where D j,k+ := [x j + A x j,k+ t, x j+ A x j+,k+ t] t], [y k+ a y j,k+ D j+,k := [x j+ a x j+,k t, x j+ D j+,k+ t, y k+ + a y j,k+ + a x j+,k t] [y k + A y t, y j+,k k+ A y t], j+,k+ := [x j+ A x j+,k+ t, x j+ + A x j+,k+ t] [y k+ A y t, y j+,k+ k+ + A y j+,k+ A x j+,k+ A y j+,k+ } := max {a xj+,k,axj+,k+, := max { a y,a y j,k+ j+,k+ }. t], Under an appropriate CFL restriction, the solution of the system (.3, subject to the initial data u(x, y, = ũ n (x, y, is smooth in D j,k (which is, in general, a nonrectangular domain inside the (j, k-cell, see Fig. 3., and may be nonsmooth in the other domains. Given the reconstruction ũ n, we integrate the system (.3 over these domains, and obtain the corresponding cell averages w n+ j,k+ = ( a y t x (A x + A x t j,k+ j+,k+ j,k+ ( ũ n (x, y dx dy D j,k+
16 698A. Kurganov, G. Petrova (3.7 t n+ t n t n+ t n y k+ +a y j,k+ t y k+ a y j,k+ t x j+ A x j+,k+ t x j +A x j,k+ t f(u(x, y, t g(u(x, y, t x j+ A x j+,k+ t x=x j +A x j,k+ t y k+ +a y j,k+ t y=y k+ a y j,k+ t dydt dx dt, = ( j+,k a x j+,k t y (A y + A y t j+,k+ j+,k ( ũ n (x, y dx dy w n+ D j+,k (3.8 t n+ t n t n+ t n y k+ A y j+,k+ t y k +A y j+,k t x j+ +a x j+,k t g(u(x, y, t x j+ a x j+,k t x j+ +a x j+,k t f(u(x, y, t dydt x=x j+ a x j+,k t y k+ A y j+,k+ t y=y k +A y j+,k t dx dt, w n+ j+,k+ = 4A x j+,k+ ( A y ( t j+,k+ ũ n (x, y dx dy (3.9 t n+ t n t n+ t n y k+ +A y j+,k+ t y k+ A y j+,k+ t x j+ +A x j+,k+ t x j+ A x j+,k+ t D j+,k+ f(u(x, y, t g(u(x, y, t x j+ +A x j+,k+ t x=x j+ A x j+,k+ t y k+ +A y j+,k+ t y=y k+ A y j+,k+ t dydt dx dt.
17 Genuinely multidimensional central scheme 699 The computation of w n+ j,k is more complicated, but it can be handled in a similar way. The spatial integrals in (3.7 (3.9 can be computed exactly and for the flux integrals one may use an appropriate quadrature formula (following the approach in [,]. As in the one-dimensional case, we complete the construction of the fully-discrete scheme in two steps. We first build a new piecewise polynomial reconstruction w n+ (x, y := [ w n+ j,k χ j,k + w n+ χ j+,k j+,k + wn+ χ j,k+ j,k+ j,k ] + w n+ χ j+,k+ j+,k+, where the χ s are the characteristic functions of the corresponding D s. Then, we project w n+ back onto the original grid and obtain the new cell averages (3. ū n+ j,k = x y x j+ y k+ x j y k w n+ (x, y dxdy. The explicit form of the resulting fully-discrete scheme is extremely complicated, and is of no practical use. Fortunately, all the computations, which have been omitted in this section, are drastically simplified if they are made in a semi-discrete context (as t. 3.3 A genuinely multidimensional third-order semi-discrete scheme In this section, we derive our new genuinely multidimensional semi-discrete scheme. We proceed as in Sect.. (see also [7,6]. Using (3., we obtain d ū n+ j,k dtūj,k(t = lim ūn j,k t t [ (3. = lim t t x y x j+ y k+ x j y k ] w n+ (x, y dxdy ū n j,k.,x j+ ] for the four corners, The following notation for the intersections of the cell [x j [y k,y k+ S j±,k, S j,k± ] with the D-domains is used: C j±,k± for the four side domains and D j,k for the center. Notice,
18 7 A. Kurganov, G. Petrova that since the size of the C-domains is proportional to ( t, and the size of the S-domains t,wehave (3. w n+ (x, y = w n+ + O(( t, j±,k± j±,k± w n+ j±,k(x, y = wn+ + O( t,,k j± (3.3 w n+ (x, y = w n+ + O( t, j,k± j,k± and due to the conservation property of the reconstruction (3.4 D j,k w n+ j,k D j,k (x, y dxdy = wn+ The substitution of (3. (3.4 in (3. yields ( d dtūj,k(t = lim t t x y + ± + ± ± S j,k± w n+ j,k± C j±,k± dxdy w n+ dxdy j±,k j,k. w n+ dxdy j±,k± (3.5 S j±,k [ ] D j,k + lim t t x y wn+ j,k ū n j,k. The first sum in (3.5 tends to zero as t, since C j±,k± = O(( t. To pass to the limit in the second sum, we use the fact that S j,k+ (3.6 = a y j,k+ t x + O(( t, and obtain lim t t x y S j,k+ w n+ j,k+ dxdy = a y j,k+ y It follows from (3.7 that the value of lim t wn+ is j,k+ lim t wn+. j,k+ x [ x j+ x j { p n j,k+ (x, y k+ +p n j,k (x, y k+
19 Genuinely multidimensional central scheme 7 g(pn j,k+ (x, y k+ g(p n j,k (x, y ] k+ } dx. a y j,k+ We then apply the Simpson s quadrature formula to the above integral, and using (3.6, we arrive at lim dxdy t t x y j,k+ (3.7 S j,k+ w n+ a y [ ] j,k+ u SW j,k+ y + unw j,k +4(uS j,k+ + un j,k +use j,k+ + une j,k y [ g(u SW j,k+ g(unw j,k +4(g(uS j,k+ g(un j,k +g(u SE j,k+ g(une j,k ]. A similar treatment of the third sum in (3.5 leads to lim w n+ dxdy t t x y j+,k (3.8 a x j+,k [ x x S j+,k u NW j+,k + une j,k +4(uW j+,k + ue j,k +usw [ f(u NW j+,k f(une j,k +4(f(uW j+,k f(ue j,k +f(u SW j+,k f(use j,k ]. j+,k + use j,k Finally, we consider the last term in (3.5. Note, that in the limit as t, the domain D j,k becomes rectangular, up to small corners of a negligible size O(( t (see Fig. 3., and thus one can easily integrate equation (.3 over D j,k [t n,t n + t]. This observation together with Simpson s rule give lim t t [ D j,k x y wn+ j,k ū n j,k a x j+,k [ ] u NE j,k 6 x +4uE j,k + use j,k a y j,k+ 6 y 6 x ] a x j,k [ u NW j,k 6 x [ ] a y [ u NW j,k +4uN j,k + une j,k j,k u SW j,k 6 y [ f(u NE j,k f(unw j,k +4(f(uE j,k f(uw j,k ] ] +4uW j,k + usw j,k ] +4uS j,k + use j,k
20 7 A. Kurganov, G. Petrova ] +f(u SE j,k f(usw j,k [ g(u NW j,k 6 y g(usw j,k +4(g(uN j,k (3.9 g(u S j,k + g(une j,k g(use j,k ]. At the end, plugging (3.7 (3.9 into (3.5, we arrive at our new semi-discrete scheme, which can be presented in the conservative form d H x j+ dtūj,k(t =,k(t Hx j,k(t x H y (t H y (t j,k+ j,k (3., y with the numerical fluxes { H x j+,k(t := f(u NW j+,k (t + f(une j,k (t+4(f(uw j+,k (t + f(ue j,k (t } a +f(u SW j+,k (t + f(use j,k }/{ x (t j+,k(t [ u NW j+,k (t une j,k (t+4(uw j,k+ (t ue j,k (t (3. +u SW j+,k (t use j,k ], (t and { H y (t := g(u SW j,k+ j,k+ (t + g(unw j,k (t+4(g(us j,k+ (t + g(un j,k (t (3. } a +g(u SE j,k+ (t + g(une j,k }/{ y (t (t j,k+ [ u SW j,k+ (t unw j,k (t+4(us j,k+ (t un j,k (t +u SE j,k+ (t une j,k ]. (t Here, a x j+,k(t,ay (t are given in (3.6, and the u s are determined by j,k+ the corresponding values of the piecewise quadratic reconstruction {p j,k } at time t, see (3.5. Remark 3.. We would like to point out that the scheme (3. (3. has been constructed as a genuinely two-dimensional Godunov-type central scheme. Our numerical experiments demonstrate that it has sharper resolution, in comparison to the dimension-by-dimension third-order central semi-discrete scheme in [6].
21 Genuinely multidimensional central scheme 73 At the same time, our new scheme enjoys all the main advantages of central schemes, the most important of which is simplicity. Namely, it does not involve Riemann solvers and does not require any information about the eigenstructure of the Jacobians f u speeds {a x j+,k} and {ay }. j,k+ and g u beyond the CFL-related local. If in the derivation of (3. (3., we use piecewise linear reconstruction and the second-order midpoint rule (instead of Simpson s rule, we arrive exactly at the scheme (.4 (.6, obtained in [7], where the dimension-by-dimension approach was used. 3. The scheme (3. can be generalized for the convection-diffusion equation as follows u t + f(u x + g(u y = Q x (u, u x,u y x + Q y (u, u x,u y y, d H x j+ dtūj,k(t =,k Hx j,k x H y j,k+ H y j,k y (3.3 Q x j+ +,k Qx Q y Q y j,k j,k+ + j,k. x y Here, the numerical convection fluxes H x, H y are given by (3.,(3. and, as in the one-dimensional case, Q x j+,k and Qy are appropriate j,k+ approximations to the diffusion fluxes. 4. Our semi-discrete approach can be easily extended to the multidimensional case, d A new two-dimensional reconstruction via dimension-by-dimension approach We generalize the idea from Sect. 3. for the two-dimensional case we construct a piecewise quadratic reconstruction {p n j,k } as a convex combination of the basic parabolas {qj,k n } and the linear functions {Ln j,k }. The basic parabolas are given by ( qj,k n (x, y = ū n j,k ( x 4 Dx +D ū x n j,k ( y 4 Dy + Dy ūn j,k +Dū x n j,k (x x j+d y ūn j,k (y y k + Dx +D x ū n j,k (x x j + D x D y ūn j,k (x x j(y y k (3.4 + Dy + Dy ūn j,k (y y k,
22 74 A. Kurganov, G. Petrova where the corresponding divided differences are D±v(x, x v(x ± x, y v(x, y y :=±, x D y y ± y v(x, y ± v(x, y :=±v(x,, y Dv(x, x v(x + x, y v(x x, y y :=, x D y y + y v(x, y y v(x, y :=v(x,. y One can easily check that, as in the one-dimensional case, these parabolas satisfy the requirements of conservation and accuracy. The piecewise linear functions have the form (3.5 L n j,k (x, y =ūn j,k + sx j,k (x x j+s y j,k (y y k, where the slopes {s x j,k,sy j,k } may be computed, for example, using the minmod limiter (3.6 (ūn s x j,k = minmod j,k ū n j,k, ūn j+,k ūn j,k x x (ūn s y j,k = minmod j,k ū n j,k, ūn j,k+ ūn j,k y y Then the reconstruction is p n j,k (x, y =( θn j,k Ln j,k (x, y+θn j,k qn j,k (x, y, (3.7 <θj,k n <. In general, the purpose of the limiters {θj,k n } is to guarantee a non-oscillatory nature of the reconstruction (3.7. Unfortunately, in the two-dimensional case, we do not know how such θ s can be chosen. Moreover, it is not even clear what the definition of a non-oscillatory propery should be. Notice, that we do not need to recover the whole reconstruction only the eight point values u N j,k, ue j,k, us j,k, uw j,k,une j,k, unw j,k, use j,k, and usw j,k in each (j, k-cell are required to compute the numerical fluxes in (3. (3.. These values may be computed using the dimension-by-dimension approach, applied in two steps. First, we use the reconstruction (3.4 (3.7 with (3.8 θj,k n = min{θx j,k,θy j,k }, to compute the point values in the coordinate directions u N j,k, ue j,k, us j,k and u W j,k (see Fig. 3.. The limiters {θx j,k } and {θy j,k } are determined by (3.3 with L n j ( =L n j,k (,y k, q j ( =qj,k n (,y k,,.
23 Genuinely multidimensional central scheme 75 y k W N S E x j Fig. 3.. Reconstruction in x- and y-directions - d + d NW NE y k SW SE x j Fig Reconstruction in the diagonal directions and L n j ( =L n j,k (x k,, q n j ( =L n j,k (x k,. The choice of {θj,k n } is based on the one-dimensional non-oscillatory reconstruction, described in Sect. 3.. This guarantees the third-order accuracy and the lack of oscillations in the x- and y-directions. To ensure that there are no oscillations in the diagonal directions, the point values u NE j,k, unw j,k, use j,k and usw j,k should be computed by another reconstruction (see Fig This additional reconstruction is determined as a convex combination of different basic parabolas {ˆq j,k n } and linear functions {ˆL n j,k }, adjusted to the diagonal directions d and d +, (3.9 ˆp n j,k (x, y =( ˆθ j,k n ˆL n j,k (x, y+ˆθ j,k n ˆqn j,k (x, y, < ˆθ j,k n <.
24 76 A. Kurganov, G. Petrova The adjusted basic parabolas are ( ˆq j,k n (x, y = ū n j,k (3.3 +D d+ ū n j,k +D d ū n j,k 48 Dd+ + D d+ +D d+ D d ū n j,k + Dd+ + D d+ ū n j,k + Dd + D d ū n j,k ū n j,k 48 Dd + D d ū n j,k [ y (y y k+ ] x (x x j [ y (y y k ] x (x x j [ 4( y (y y k 4( x (x x j ] [ y (y y k+ ] x (x x j [ y (y y k ], x (x x j where := ( x +( y, and the divided differences in the diagonal directions are v(x ± x, y ± y v(x, y D± d+ v(x, y :=±, v(x + x, y + y v(x x, y y D d+ v(x, y :=, v(x x, y ± y v(x, y D± d v(x, y :=±, D d v(x, y := v(x x, y + y v(x + x, y y. The corresponding linear functions are given by [ ˆL n j,k (x, y =ūn j,k +ŝ+ j,k y (y y k+ ] x (x x j [ +ŝ j,k y (y y k ] (3.3 x (x x j, where the slopes are computed using the minmod limiter, applied in the diagonal directions, that is, (3.3, (3.33 (ūn ŝ + j,k = minmod j,k ū n j,k, ūn j+,k+ ūn j,k (ūn ŝ j,k = minmod j,k ū n j+,k, ūn j,k+ ūn j,k The parameter ˆθ j,k n in (3.9 is determined by. (3.34 ˆθ n j,k = min{ˆθ + j,k, ˆθ j,k },
25 Genuinely multidimensional central scheme 77 where ˆθ + j,k and ˆθ j,k are designed to prevent oscillations in the diagonal directions d + and d, respectively. The ˆθ s are computed similarly to (3.3, namely { M ± ˆL n j± min,k+ j,k (x j±,y k+ M ± j,k ˆL n j,k (x, j±,y k+ m ± ˆL n j,k j,k (x } j,y k m ± j,k ˆL n j,k (x,, j,y k if ū n j,k < ūn j,k < ūn j±,k+, (3.35 θ ± j,k := { M ± ˆL n j min,k j,k (x j,y k M ± j,k ˆL n j,k (x, j,y k m ± ˆL n j±,k+ j,k (x } j±,y k+ m ± j,k ˆL n j,k (x,, j±,y k+ if ū n j,k > ūn j,k > ūn j±,k+, where and M + j±,k± m + j±,k± M j,k± m j,k±, otherwise, { M ± j,k = max ˆq j,k n (x j±,y k+, ˆq j,k n (x j { m ± j,k = min ˆq j,k n (x j±,y k+, ˆq j,k n (x j { = max (ˆLn j,k (x j± ˆq n j±,k± (x j±,y k± { = min (ˆLn j,k (x j± ˆq n j±,k± (x j±,y k± { = max (ˆLn j,k (x j ˆq n j,k± (x j { = min,y k±,y k±,y k,y k }, } +ˆL n j±,k± (x j± },,y k± +ˆL n j±,k± (x j± },,y k± +ˆL n j,k± (x j }, (ˆLn j,k (x j,y k± +ˆL n j,k± (x j,,y k±,,y k±,,y k±,,y k±,
26 78A. Kurganov, G. Petrova } ˆq j,k± n (x j,y k±. Remark 3.3. Direct calculations show that the reconstruction ˆp n j,k (x, y conserves the original cell averages, that is x y x j+ x j y k+ y k ˆp n j,k (x, y dx dy =ūn j,k.. We would like to emphasize that we have actually used two different piecewise quadratic reconstructions, (3.4 (3.8 and (3.9 (3.35, in order to compute the eight point values along the inner boundary of the (j, k-cell. In the derivation of our genuinely two-dimensional scheme (3. (3. we use Simpson s rule, which is applied to a single smooth function, defined on the (j, k-cell. For example, this smooth function can be viewed as the following convex combination, ( ϕ(x, yp n j,k (x, y+ϕ(x, yˆpn j,k (x, y, where the weight function ϕ can be any smooth function satisfying ϕ(x j± ϕ(x j±,y k =ϕ(x j,y k± =, =, ϕ(x, y.,y k± In summary, the values u N j,k, us j,k, ue j,k and uw j,k are computed as the corresponding point values of p n j,k (x, y ( = ū n j,k ( x θn j,k 4 Dx +D ū x n j,k ( y θn j,k ( + θj,k n Dx ū n j,k +( θn j,k sx j,k (x x j ( + θj,k n Dy ūn j,k +( θn j,k sy j,k (y y k 4 Dy + Dy ūn j,k + θn j,k Dx +D ū x n j,k (x x j + θj,k n Dx D y ūn j,k (x x j(y y k + θn j,k (3.36 Dy + Dy ūn j,k (y y k, where θj,k n, sx j,k and sy j,k are given by (3.8 and (3.6, respectively. To evaluate the corner values u NE j,k, unw j,k, use j,k and usw j,k,weuse ˆp n j,k (x, y
27 Genuinely multidimensional central scheme 79 (3.37 ( = ū n j,k ˆθ j,k n 48 (Dd+ + D d+ + D+ d D d ū n j,k + (ˆθn x j,k (D d+ D d ū n j,k +( ˆθ j,k n (ŝ+ j,k ŝ j,k (x x j + y (ˆθn j,k (D d+ + D d ū n j,k +( ˆθ n j,k (ŝ+ j,k +ŝ j,k (y y k +ˆθ j,k n ( 8( x D+ d+ D d+ +ˆθ n j,k +ˆθ n j,k 4 x y ( 8( y D d+ D d + D d + D d ū n j,k (x x j ( D+ d+ D d+ D+ d D d ū n j,k (x x j(y y k D+ d+ D d+ +D d+ D d + D d + D d ū n j,k (y y k, with ˆθ j,k n, ŝ+ j,k and ŝ j,k, given by (3.34, (3.3 and (3.33, respectively. 4 Numerical examples In this section, we present a number of numerical experiments that have been performed using the one-dimensional scheme (.9 (. together with the new reconstruction (3.4, and our genuinely two-dimensional scheme (3. (3., coupled with the new piecewise quadratic reconstruction (3.36, (3.37. The third-order semi-discrete methods, presented in this paper, require at least a third-order time discretization scheme. In the numerical examples, shown below, we have used the third-order TVD Runge-Kutta method, proposed in [34,36]. Our choice is based on the stability properties of this method, each time step of which can be viewed as a convex combination of small forward Euler steps. To solve problem (., we may use either an implicit or an explicit time discretization for the parabolic term. To ensure stability of the explicit method, an additional restriction on the time step is imposed by the parabolic term, namely t must be proportional to ( x. One may accelerate the computations by using an implicit, or explicit-implicit ODE solver. These methods are unconditionally stable, but they require inverting of nonlinear operators (in the case of nonlinear diffusion, which is a computationally expensive and analytically complicated procedure. Alternatively, the acceleration may be achieved by employing an efficient explicit Runge-Kutta type ODE solver with larger stability domain (in comparison with the standard Runge-Kutta methods, see [3] and the references therein. In all the numerical experiments below, the CFL number is equal to.475.
28 7 A. Kurganov, G. Petrova 4. Linear accuracy tests Example One-dimensional transport equation. First, we test the accuracy of the scheme (.9 (. with the reconstruction (3.4 on the initial boundary value problem (IBVP for the linear transport equation { ut + u x =, x [, π], (4. u(x, = sin x, with periodic boundary conditions. This problem has a global classical solution, which is computed at time T =. We consider N grid points, N =4, 8,...,8, and measure the L - and L -errors, respectively. To calculate the errors, we use the computed values {ū j } at the final time, and the sliding averages {ū(x j, } of the exact solution of (4., namely ū ū(, L := x j ū j ū(x j,, ū ū(, L := max ū j ū(x j,. j The results of these computations are presented in Table 4.. They clearly demonstrate that the scheme is third-order. Example Two-dimensional transport equation. Second, we apply our two-dimensional scheme (3. (3. with the reconstruction (3.36 (3.37 to the IBVP for the two-dimensional linear transport equation { ut + u x + u y =, (x, y [, ] [, ], (4. u(x, y, = sin[π(x +y], subject to periodic boundary conditions. We again calculate the L - and L - errors at time T =, using the cell averages of the computed and the exact solutions of (4.. As in the one-dimensional case, the results, presented in Table 4., indicate the third-order convergence rate. Table 4.. Accuracy test for the linear advection problem (4., T= N L -error rate L -error rate 4.355e e e e e e e e e e e e-8 3.
29 Genuinely multidimensional central scheme 7 Table 4.. Accuracy test for the linear advection problem (4., T= Grid L -error rate L -error rate 6 6.6e- 9.4e e-.89.7e e e e e e e One-dimensional problems Example 3 Burgers equation. In this example, we solve the IBVP for the one-dimensional Burgers equation ( u u t + =, x [, π], (4.3 x u(x, =.5 + sin x, with periodic boundary conditions. It is known that the unique entropy solution of (4.3 develops a shock discontinuity at time t =. In Fig. 4. we present the approximate solution at the post-shock time T =, computed by the scheme (.9 (. with the reconstruction (3.4. Notice, that even though the method provides a high resolution of the shock, one can observe the over- and undershootings near the discontinuity. This happens because the limiter θj n, defined in (3.3, is switched off (i.e., θj n = at the local extrema. It is possible to reduce these oscillations with the following recipe we choose θj n to be { ( θj n = max min qn j (x j+ ū n j α+3 ( x α, qn j (x j ū n } j α+3 (4.4 ( x α,, α, when ū n j > ūn j±,orūn j < ūn j±. The modified reconstruction (3.4,(4.4 is still third-order in smooth regions, and at the same time, it reduces the oscillations near the discontinuities. In general, any α can be used in (4.4. Our numerical experiments have not indicated which value of α is optimal. Our experience shows that larger α s lead to smaller oscillations, but increase the numerical dissipation. In this example (see Fig. 4. and in the examples below, we have used α =.
30 7 A. Kurganov, G. Petrova.5 EXACT N= Fig. 4.. Burgers equation (4.3; using reconstruction (3.4.5 EXACT N= Fig. 4.. Burgers equation (4.3; using reconstruction (3.4,(4.4 Example 4 One-dimensional Euler equations of gas dynamics. Here we consider the one-dimensional Euler system, (4.5 t ρ m E + x m ρu + p =, p =(γ u(e + p ( E ρ u, where ρ, u, m = ρu, p and E are the density, velocity, momentum, pressure and the total energy, respectively, and γ =.4. We solve this system with the initial data (4.6 u L =(,, 5 T, u(x, = u M =(,,.5 T, u R =(,, 5 T, x<.,. x<.9,.9 x<,
31 Genuinely multidimensional central scheme 73 7 N=6 N= Fig Problem (4.5 (4.6, density at T=.; using reconstruction (3.4 7 N=6 N= Fig Problem (4.5 (4.6, density at T=.; using reconstruction (3.4, (4.4 and solid boundary conditions, applied to both ends. The example describes the interaction of blast waves and was proposed by Woodward and Colella in [37]. To compute the approximate solution of (4.5 (4.6, we use the scheme (.9 (. with the reconstruction (3.4. The computations are done, using N = 4 grid points, and the solution is plotted together with a reference solution, obtained by the same method with N = 6. Figures 4.3, 4.5, and 4.7 show the density, the velocity, and the pressure at time T =.. Notice, that if we use 4 grid points, the second density spike has a height of 5.75, which is better in comparison with the heights of 5. obtained by the third-order staggered central scheme in [9], and 3.7, obtained by the second-order Nessyahu-Tadmor scheme in [3]. This illustrates the higher resolution and smaller numerical dissipation of our method.
32 74 A. Kurganov, G. Petrova N=6 N= Fig Problem (4.5 (4.6, velocity at T=.; using reconstruction (3.4 N=6 N= Fig Problem (4.5 (4.6, velocity at T=.; using reconstruction (3.4, (4.4 Figures 4.4, 4.6, and 4.8show the solution of the same problem, computed again by the scheme (.9 (., but this time coupled with the modified reconstruction (3.4,(4.4. The height of the second density spike here is 5.5 which, as we have mentioned above, is due to the more dissipative nature of the modified reconstruction. The advantage of this approach is that it reduces the oscillations, as one can see on Figs. 4.5 and 4.6. We also perform the computations at time T =.3. The results are presented in Figs We would like to point out that for 4 grid points the maximum value of the density is 3 (see Fig. 4.9, and it is (see Fig. 4.9, if we apply the limiter (4.4 at local extrema. Notice, that for the more dissipative third-order staggered central scheme in [9] this height is. Finally, we compute the solution of (4.5 (4.6 at time T =.38, shown in Figs Here, the value of the second density spike is 5.5 (see
33 Genuinely multidimensional central scheme N=6 N= Fig Problem (4.5 (4.6, pressure at T=.; using reconstruction ( N=6 N= Fig Problem (4.5 (4.6, pressure at T=.; using reconstruction (3.4, (4.4. Fig. 4.7, and 5., if we apply (4.4 (see Fig The corresponding number from [9] is Two-dimensional problems Example 5 Two-dimensional convection-diffusion problem. We consider the two-dimensional Burgers-type equation (4.7 u t +(u x +(u y = ε(ν(uu x x + ε(ν(uu y y, (x, y [.5,.5] [.5,.5], with a strongly degenerate diffusion coefficient {, u.5, ν(u =, u >.5.
34 76 A. Kurganov, G. Petrova 5 N=6 N= Fig Problem (4.5 (4.6, density at T=.3; using reconstruction (3.4 5 N=6 N= Fig. 4.. Problem (4.5 (4.6, density at T=.3; using reconstruction (3.4, (4.4 4 N=6 N= Fig. 4.. Problem (4.5 (4.6, velocity at T=.3; using reconstruction (3.4
35 Genuinely multidimensional central scheme 77 4 N=6 N= Fig. 4.. Problem (4.5 (4.6, velocity at T=.3; using reconstruction (3.4, ( N=6 N= Fig Problem (4.5 (4.6, pressure at T=.3; using reconstruction ( N=6 N= Fig Problem (4.5 (4.6, pressure at T=.3; using reconstruction (3.4, (4.4
36 78A. Kurganov, G. Petrova 9 8 N=6 N= Fig Problem (4.5 (4.6, pressure at T=.3; using reconstruction (3.4; zoom at [.65,.75] 9 8 N=6 N= Fig Problem (4.5 (4.6, pressure at T=.3; using reconstruction (3.4, (4.4; zoom at [.65,.75] 7 N=6 N= Fig Problem (4.5 (4.6, density at T=.38; using reconstruction (3.4
37 Genuinely multidimensional central scheme 79 7 N=6 N= Fig Problem (4.5 (4.6, density at T=.38; using reconstruction (3.4, ( N=6 N= Fig Problem (4.5 (4.6, velocity at T=.38; using reconstruction ( N=6 N= Fig. 4.. Problem (4.5 (4.6, velocity at T=.38; using reconstruction (3.4, (4.4
38 7 A. Kurganov, G. Petrova 45 4 N=6 N= Fig. 4.. Problem (4.5 (4.6, pressure at T=.38; using reconstruction ( N=6 N= Fig. 4.. Problem (4.5 (4.6, pressure at T=.38; using reconstruction (3.4, (4.4 This equation is of hyperbolic nature when u [.5,.5], and is parabolic elsewhere. We solve (4.7, subject to the initial data u(x, y, =, if (x.5 +(y.5.6,, if (x +.5 +(y +.5.6,, otherwise. The numerical experiments are performed for ε =., and for ε =.In the first case we use the scheme (3.3, and in the second, pure hyperbolic case, we apply the scheme (3. (3., coupled with the reconstruction (3.36 (3.37. Figures 4.3 and 4.4 show the computed solutions at time T =.5 in the hyperbolic and the hyperbolic-parabolic case, respectively. We would like to emphasize, that in both cases the resolution of the shock discontinuities is very high, and the transition between the hyperbolic and parabolic regions
39 Genuinely multidimensional central scheme 7 y Fig Pure hyperbolic problem (4.7 with ε =, T=.5; 6 6 grid y x Fig Degenerate parabolic problem (4.7 with ε =., T=.5; 6 6 grid x.8.6 PURE HYPERBOLIC PROBLEM CONVECTION-DIFFUSION PROBLEM Fig Cross-section along the line y = x
CS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationOn the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes
Journal of Computational Physics 8, 87 09 (00) doi:0.006/jcph.00.79 On the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes Jianxian Qiu, and Chi-Wang
More informationCentral-Upwind Schemes on Triangular Grids for Hyperbolic Systems of Conservation Laws
Central-Upwind Schemes on Triangular Grids for Hyperbolic Systems of Conservation Laws Alexander Kurganov, 1 Guergana Petrova 2 1 Department of Mathematics, Tulane University, New Orleans, Louisiana 70118
More informationNumerical Methods for Hyperbolic and Kinetic Equations
Numerical Methods for Hyperbolic and Kinetic Equations Organizer: G. Puppo Phenomena characterized by conservation (or balance laws) of physical quantities are modelled by hyperbolic and kinetic equations.
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
More informationNon-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. 933-963 Commun. Comput. Phys. October 7 Non-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes Yingjie Liu,, Chi-Wang Shu,
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More information1. Introduction. Numerical methods for approximating solutions of hyperbolic conservation laws,
SIAM J. SCI. COMPUT. Vol., No. 4, pp. 1461 1488 c 000 Society for Industrial and Applied Mathematics A THIRD-ORDER SEMIDISCRETE CENTRAL SCHEME FOR CONSERVATION LAWS AND CONVECTION-DIFFUSION EQUATIONS ALEXANDER
More informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More informationHierarchical Reconstruction for Spectral Volume Method on Unstructured Grids
Hierarchical Reconstruction for Spectral Volume Method on Unstructured Grids Zhiliang Xu, Yingjie Liu and Chi-Wang Shu April 4, 2009 Abstract The hierarchical reconstruction (HR) [, 24] is applied to a
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
More informationObservations on the fifth-order WENO method with non-uniform meshes
Observations on the fifth-order WENO method with non-uniform meshes Rong Wang, a, Hui Feng, b, Raymond J. Spiteri a,, a Department of Computer Science, University of Saskatchewan, Saskatoon, SK, S7N 5C9,
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
AIAA SciTech Forum 8 January 08, Kissimmee, Florida 08 AIAA Aerospace Sciences Meeting 0.54/6.08-0067 A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley and Dr.
More informationNON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES
NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. This is the continuation of the paper central discontinuous
More informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationA SIMPLE EULERIAN FINITE-VOLUME METHOD FOR COMPRESSIBLE FLUIDS IN DOMAINS WITH MOVING BOUNDARIES
A SIMPLE EULERIAN FINITE-VOLUME METHOD FOR COMPRESSIBLE FLUIDS IN DOMAINS WITH MOVING BOUNDARIES ALINA CHERTOCK AND ALEXANDER KURGANOV Abstract. We introduce a new simple Eulerian method for treatment
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationA Jiang Tadmor Scheme on Unstructured Triangulations
A Jiang Tadmor Scheme on Unstructured Triangulations Ivan Christov, Bojan Popov Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA Abstract Nonoscillatory central schemes
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationCentral Runge-Kutta Schemes for Conservation Laws
Central Runge-Kutta Schemes for Conservation Laws Lorenzo Pareschi Gabriella Puppo Giovanni Russo Abstract In this work, a new formulation for central schemes based on staggered grids is proposed. It is
More informationHigh Order Schemes for CFD: A Review. Juan Cheng 1. Institute of Applied Physics and Computational Mathematics, Beijing , China.
High Order Schemes for CFD: A Review Juan Cheng 1 Institute of Applied Physics and Computational Mathematics, Beijing 100088, China and Chi-Wang Shu 2 Division of Applied Mathematics, Brown University,
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 13: The Lecture deals with:
The Lecture deals with: Some more Suggestions for Improvement of Discretization Schemes Some Non-Trivial Problems with Discretized Equations file:///d /chitra/nptel_phase2/mechanical/cfd/lecture13/13_1.htm[6/20/2012
More informationOn the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes
Commun. Comput. Phys. doi: 8/cicp.99.84s Vol. 9, No. 3, pp. 87-87 March On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes Rui Zhang, Mengping Zhang and Chi-Wang
More informationDiscontinuous Fluctuation Distribution
Discontinuous Fluctuation Distribution Matthew Hubbard School of Computing, University of Leeds, Leeds, LS2 9JT, UK. Abstract This paper describes a new numerical scheme for the approximation of steady
More informationCentPack version 1.0 User s Guide
CentPack version 1.0 User s Guide Jorge Balbás University of Michigan, Ann Arbor Eitan Tadmor University of Maryland, College Park July 19, 2006 http://www.cscamm.umd.edu/centpack/ c 2006 Jorge Balbás
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationThird Order WENO Scheme on Three Dimensional Tetrahedral Meshes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu
More informationGENUINELY MULTIDIMENSIONAL NON-DISSIPATIVE FINITE VOLUME SCHEMES FOR TRANSPORT
GENUINELY MULTIDIMENSIONAL NON-DISSIPATIVE FINITE VOLUME SCHEMES FOR TRANSPORT BRUNO DESPRÉS, FRÉDÉRIC LAGOUTIÈRE Commissariat à l Énergie Atomique, Bruyères-le-Châtel. despres@cmpax.polytechnique.fr e-mail:
More informationAPPROXIMATING PDE s IN L 1
APPROXIMATING PDE s IN L 1 Veselin Dobrev Jean-Luc Guermond Bojan Popov Department of Mathematics Texas A&M University NONLINEAR APPROXIMATION TECHNIQUES USING L 1 Texas A&M May 16-18, 2008 Outline 1 Outline
More informationSpectral(Finite) Volume Method for Conservation Laws on Unstructured Grids
Journal of Computational Physics 179, 665 697 (2002) doi:10.1006/jcph.2002.7082 Spectral(Finite) Volume Method for Conservation Laws on Unstructured Grids II. Extension to Two-Dimensional Scalar Equation
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More information= f (a, b) + (hf x + kf y ) (a,b) +
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
More informationWeno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach
Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 6, No. 3, pp. 785 89 c 5 Society for Industrial and Applied Mathematics A HIGH-RESOLUTION ROTATED GRID METHOD FOR CONSERVATION LAWS WITH EMBEDDED GEOMETRIES CHRISTIANE HELZEL,
More informationA Review on the Numerical Solution of the 1D Euler Equations. Hudson, Justin. MIMS EPrint:
A Review on the Numerical Solution of the D Euler Equations Hudson, Justin 6 MIMS EPrint: 6.9 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available
More informationHigh Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study
Commun. Comput. Phys. doi:.48/cicp.375.6a Vol., No. 4, pp. 835-869 October 6 High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study Liang
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationHigh-Order Finite Difference Schemes for computational MHD
High-Order Finite Difference Schemes for computational MHD A. Mignone 1, P. Tzeferacos 1 and G. Bodo 2 [1] Dipartimento di Fisica Generale, Turin University, ITALY [2] INAF Astronomic Observatory of Turin,,
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationA third order conservative Lagrangian type scheme on curvilinear meshes for the compressible Euler equations
A third order conservative Lagrangian type scheme on curvilinear meshes for the compressible Euler equations Juan Cheng and Chi-Wang Shu 2 Abstract Based on the high order essentially non-oscillatory (ENO)
More information1. Introduction. We consider Cauchy problems for Hamilton Jacobi (HJ) equations
SIAM J. SCI. COMPUT. Vol. 8, No. 6, pp. 9 47 c 6 Society for Industrial and Applied Mathematics CENTRAL WENO SCHEMES FOR HAMILTON JACOBI EQUATIONS ON TRIANGULAR MESHES DORON LEV, SUHAS NAAK, CHI-WANG SHU,
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationA CONSERVATIVE FRONT TRACKING ALGORITHM
A CONSERVATIVE FRONT TRACKING ALGORITHM Vinh Tan Nguyen, Khoo Boo Cheong and Jaime Peraire Singapore-MIT Alliance Department of Mechanical Engineering, National University of Singapore Department of Aeronautics
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationRunge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh 1
Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh Jun Zhu, inghui Zhong 3, Chi-Wang Shu 4 and Jianxian Qiu 5 Abstract In this paper we generalize a new type
More informationNUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD
18th Engineering Mechanics Division Conference (EMD007) NUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD Abstract S. Fu University of Texas at Austin, Austin,
More informationLax-Wendroff and McCormack Schemes for Numerical Simulation of Unsteady Gradually and Rapidly Varied Open Channel Flow
Archives of Hydro-Engineering and Environmental Mechanics Vol. 60 (2013), No. 1 4, pp. 51 62 DOI: 10.2478/heem-2013-0008 IBW PAN, ISSN 1231 3726 Lax-Wendroff and McCormack Schemes for Numerical Simulation
More informationA Technique of Treating Negative Weights in WENO Schemes
NASA/CR--63 ICASE Report No. -49 A Technique of Treating Negative Weights in WENO Schemes Jing Shi, Changqing Hu, and Chi-Wang Shu Brown University, Providence, Rhode Island ICASE NASA Langley Research
More informationHigh-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh
High-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh by Al-Amin Aziz A thesis submitted in conformity with the requirements for the degree of Masters of Applied Science Graduate
More informationA high order moving boundary treatment for compressible inviscid flows 1. Abstract
A high order moving boundary treatment for compressible inviscid flows Sirui Tan and Chi-Wang Shu Abstract We develop a high order numerical boundary condition for compressible inviscid flows involving
More informationNon-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws
Non-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws Matthew Hubbard School of Computing, University of Leeds, Leeds, LS 9JT, UK. Abstract This paper addresses
More informationComputing Nearly Singular Solutions Using Pseudo-Spectral Methods
Computing Nearly Singular Solutions Using Pseudo-Spectral Methods Thomas Y. Hou Ruo Li January 9, 2007 Abstract In this paper, we investigate the performance of pseudo-spectral methods in computing nearly
More informationDDFV Schemes for the Euler Equations
DDFV Schemes for the Euler Equations Christophe Berthon, Yves Coudière, Vivien Desveaux Journées du GDR Calcul, 5 juillet 2011 Introduction Hyperbolic system of conservation laws in 2D t W + x f(w)+ y
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationChapter 6. Petrov-Galerkin Formulations for Advection Diffusion Equation
Chapter 6 Petrov-Galerkin Formulations for Advection Diffusion Equation In this chapter we ll demonstrate the difficulties that arise when GFEM is used for advection (convection) dominated problems. Several
More informationA limiting strategy for the back and forth error compensation and correction method for solving advection equations
A limiting strategy for the back and forth error compensation and correction method for solving advection equations Lili Hu, Yao Li, Yingjie Liu April 22, 2014 Abstract We further study the properties
More informationAN ADAPTIVE MESH REDISTRIBUTION ALGORITHM FOR CONVECTION-DOMINATED PROBLEMS. Zheng-Ru Zhang and Tao Tang
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume, Number 3, september pp. 34 357 AN ADAPTIVE MESH REDISTRIBUTION ALGORITHM FOR CONVECTION-DOMINATED PROBLEMS Zheng-Ru Zhang
More informationAn adaptive discretization of compressible flow using a multitude of moving Cartesian grids
An adaptive discretization of compressible flow using a multitude of moving Cartesian grids Linhai Qiu, Wenlong Lu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationSemi-Conservative Schemes for Conservation Laws
Semi-Conservative Schemes for Conservation Laws Rosa Maria Pidatella 1 Gabriella Puppo, 2 Giovanni Russo, 1 Pietro Santagati 1 1 Università degli Studi di Catania, Catania, Italy 2 Università dell Insubria,
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationA New Trouble-Cell Indicator for Discontinuous Galerkin Methods for. Hyperbolic Conservation Laws ABSTRACT
A New Trouble-Cell Indicator for Discontinuous Galerkin Methods for Hyperbolic Conservation Laws Guosheng Fu and Chi-Wang Shu ABSTRACT We introduce a new troubled-cell indicator for the discontinuous Galerkin
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationOn the simulation of wave propagation with a higher-order finite volume scheme based on Reproducing Kernel Methods
On the simulation of wave propagation with a higher-order finite volume scheme based on Reproducing Kernel Methods X. Nogueira a, I. Colominas,a, L. Cueto-Felgueroso b, S. Khelladi c a Group of Numerical
More informationAn Introduction to Viscosity Solutions: theory, numerics and applications
An Introduction to Viscosity Solutions: theory, numerics and applications M. Falcone Dipartimento di Matematica OPTPDE-BCAM Summer School Challenges in Applied Control and Optimal Design July 4-8, 2011,
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationComputational Astrophysics 5 Higher-order and AMR schemes
Computational Astrophysics 5 Higher-order and AMR schemes Oscar Agertz Outline - The Godunov Method - Second-order scheme with MUSCL - Slope limiters and TVD schemes - Characteristics tracing and 2D slopes.
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationA-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain
2th Annual CFD Symposium, August 9-1, 218, Bangalore A-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain S. M. Joshi & A. Chatterjee Department of Aerospace Engineering Indian Institute
More informationA new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes. Abstract
A new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes Xuliang Liu 1, Shuhai Zhang, Hanxin Zhang 3 and Chi-Wang Shu 4 Abstract In this paper, we develop
More informationStudy on the Numerical Accuracy for the CFD
Study on the Numerical Accuracy for the CFD T.Yamanashi 1, H.Uchida, and M.Morita 1 Department of Mathematics, Master s Research Course of Faculty of Science, Tokyo University of Science,1-3 Kagurazaka,
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationMath 690N - Final Report
Math 690N - Final Report Yuanhong Li May 05, 008 Accurate tracking of a discontinuous, thin and evolving turbulent flame front has been a challenging subject in modelling a premixed turbulent combustion.
More informationATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, Paul A. Ullrich (HH 251)
ATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, 2013 Paul A. Ullrich (HH 251) paullrich@ucdavis.edu Outline 1. Introduction / Motivation 2. Finite Difference Methods
More informationAn Adaptive Stencil Linear Deviation Method for Wave Equations
211 An Adaptive Stencil Linear Deviation Method for Wave Equations Kelly Hasler Faculty Sponsor: Robert H. Hoar, Department of Mathematics ABSTRACT Wave Equations are partial differential equations (PDEs)
More informationThis is the main idea of the evolution Galerkin (EG) methods, which evolve the initial data using the bicharacteristic cone and then project them onto
Finite volume evolution Galerkin methods for multidimensional hyperbolic problems M. Lukacova { Medvid'ova 1 3, K. W. Morton 2,G.Warnecke 1 1 Institut fur Analysis und Numerik, Otto-von-Guericke-Universitat
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationHigh order regularization of Dirac-delta sources in two space dimensions.
INTERNSHIP REPORT High order regularization of Dirac-delta sources in two space dimensions. Author: Wouter DE VRIES s1010239 Supervisors: Prof. Guustaaf JACOBS Jean-Piero SUAREZ HOST INSTITUTION: San Diego
More informationDownloaded 10/27/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 37, No. 3, pp. C354 C383 c 2015 Society for Industrial and Applied Mathematics EFFICIENT IMPLEMENTATION OF NONLINEAR COMPACT SCHEMES ON MASSIVELY PARALLEL PLATFORMS DEBOJYOTI
More informationAccurate Monotonicity-Preserving Schemes with Runge Kutta Time Stepping
JOURNAL OF COMPUTATIONAL PHYSICS 136, 83 99 (1997) ARTICLE NO. CP975745 Accurate Monotonicity-Preserving Schemes with Runge Kutta Time Stepping A. Suresh*, and H. T. Huynh *NYMA, Inc., Brook Park, Ohio
More informationRobust and Efficient Adaptive Moving Mesh Solution of the 2-D Euler equations
Contemporary Mathematics Robust and Efficient Adaptive Moving Mesh Solution of the -D Euler equations P. A. Zegeling, W. D. de Boer, and H. Z. Tang Abstract. In this paper we describe an adaptive moving
More informationA HYBRID SEMI-PRIMITIVE SHOCK CAPTURING SCHEME FOR CONSERVATION LAWS
Eighth Mississippi State - UAB Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conf. 9 (), pp. 65 73. ISSN: 7-669. URL: http://ejde.math.tstate.edu
More informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationRunge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes
Runge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes Jun Zhu, Jianxian Qiu,Chi-WangShu 3 and Michael Dumbser 4 Abstract In [], Qiu and Shu investigated using weighted essentially
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationA new type of multi-resolution WENO schemes with increasingly higher order of accuracy
A new type of multi-resolution WENO schemes with increasingly higher order of accuracy Jun Zhu and Chi-Wang Shu Abstract In this paper, a new type of high-order finite difference and finite volume multi-resolution
More information