Spectral collocation and waveform relaxation methods with Gegenbauer reconstruction for nonlinear conservation laws
|
|
- Stanley Jenkins
- 5 years ago
- Views:
Transcription
1 Spectral collocation and waveform relaxation methods with Gegenbauer reconstruction for nonlinear conservation laws Z. Jackiewicz and B. Zubik Kowal November 2, 2004 Abstract. We investigate Chebyshev spectral collocation and waveform relaxation methods for nonlinear conservation laws. Waveform relaxation methods allow to replace the system of nonlinear differential equations resulting from the application of spectral collocation methods by a sequence of linear problems which can be effectively integrated by highly stable implicit methods. The obtained numerical solution is then enhanced on the intervals of smoothness by Gegenbauer reconstruction. The effectiveness of this approach is illustrated by numerical experiments. Keywords: Nonlinear conservation law, pseudospectral methods, waveform relaxation iterations, Gegenbauer reconstruction. Department of Mathematics, Arizona State University, Tempe, Arizona 85287, jackiewi@math.la.asu.edu. The work of this author was partially supported by the National Science Foundation under grant NSF DMS Department of Mathematics, Boise State University, 1910 University Drive, Boise, Idaho 83725, zubik@diamond.boisestate.edu. 1
2 1 Introduction It is the purpose of this paper to investigate spectral collocation and waveform relaxation methods for the nonlinear conservation law u(x, t) + t x f ( u(x, t) ) = 0, L x L, t, u(x, ) = g(x), L x L, (1.1) u( L, t) = α(t), u(l, t) = β(t), t. The solution u(x, t) to (1.1) is a limit of u ν (x, t) as ν 0, where u ν (x, t) satisfies the boundary value problem with added viscosity ν t uν (x, t) + x f ( u ν (x, t) ) = ν 2 x (x, t), L x L, t, u ν (x, ) = g(x), L x L, u ν ( L, t) = α(t), u ν (L, t) = β(t), t, (1.2) and it may be beneficial to consider (1.2) with small viscosity ν instead of (1.2) to stabilize the resulting schemes. Spectral and Legendre pseudospectral viscosity methods for nonlinear conservation laws have been considered by Gottlieb, Lustman and Orszag [17], Maday, Ould Kaber and Tadmor [32], and Guo, Ma and Tadmor [24]. In the next section we will describe the Chebyshev pseudospectral methods for (1.1) and (1.2). This leads to the nonlinear system of ordinary differential equations (ODEs) which will be integrated in time by explicit Runge Kutta methods of order four. The resulting numerical schemes are spectrally accurate if the solution to (1.1) or (1.2) is smooth. For discontinuous problems we can improve the convergence away from the discontinuity by applying a filter of order p. High order filters allow for high resolution away from the discontinuities with strong oscillations (Gibbs phenomenon) in the neighborhoods of the discontinuities. The application of low order filters will reduce these oscillations, but at the expense of excessive blurring. In this paper we will use the former approach (high order filters) with subsequent application of Gegenbauer reconstruction, as in [24], to resolve the Gibbs phenomenon and enhance the accuracy of the resulting numerical approximations. This will be described in Section 3. In Section 4 the 2
3 results of numerical experiments are presented which illustrate the effectiveness of pseudospectral Chebyshev methods with Gegenbauer reconstruction for discontinuous problems. In Section 5 we will discuss waveform relaxation methods for (1.1) and (1.2) in conjunction with Chebyshev pseudospectral methods. These methods can be viewed as a way to replace the nonlinear systems of differential equations obtained by semidiscretization in space by the method of lines of (1.1) or (1.2) by a sequence of linear problems which are easier to solve. As a result, we can employ implicit numerical schemes for integration in time. These methods have much better stability properties than the explicit formulas employed in Section 2 and, as a consequence, allow the use of much larger stepsizes for time integration which fulfill stability restrictions, as compared with the explicit Runge Kutta schemes. This is confirmed by numerical experiments presented in Section 6. The additional advantage of waveform relaxation is the fact that depending on the choice of splitting we can decouple the resulting system of differential equations which can then be efficiently integrated in a parallel computing environment. Finally, in Section 7 some concluding remarks are given and plans for future research are briefly outlined. 2 Pseudospectral Chebyshev method We will describe the pseudospectral Chebyshev method for (1.2), where to simplify the notation, we will skip the dependence of the solution to this problem on the parameter ν. Let M be a nonnegative integer and denote by x i = L cos(πi/m), i = 0, 1,..., M, the Chebyshev Gauss Lobatto points in the interval [ L, L]. We discretize (1.2) in space by the method of lines replacing 2 f(u(x, t)) and u(x, t) by pseudospectral approximations given x x 2 by and i = 0, 1,..., M. Here, x f ( u(x i, t) ) f ( u(x i, t) ) M 2 x 2 u(x i, t) M j=0 j=0 D (k) = [ d (k) ] M ij, k = 1, 2, i,j=0 3 d (1) ij u(x j, t), (2.1) d (2) ij u(x j, t), (2.2)
4 are differentiation matrices of order k. The explicit expressions for these matrices are given in [7], [10] and [38]. Put u i (t) = u(x i, t). Substituting (2.1) and (2.2) into (1.2) and taking into account that u 0 (t) = α(t), u M (t) = β(t), we obtain M 1 u i(t) + f (u i (t)) d (1) M 1 ij u j (t) = ν d (2) ij u j (t) j=1 j=1 f (u i (t)) ( d (1) i0 α(t) + d (1) imβ(t) ) + ν ( d (2) i0 α(t) + d (2) imβ(t) ), (2.3) u i ( ) = g(x i ), i = 1, 2,..., M 1. Put and D (k) = [ d (k) ] M 1 ij, i,j=1 a(k) = u(t) = [ u 1 (t) u M 1 (t) d (k) 1,0. d (k) M 1,0, b(k) = ] T, g = [ f (u(t)) = [ f (u 1 (t)) f (u M 1 (t)) d (k) 1,M. d (k) M 1,M g(x 1 ) g(x M 1 ) ] T., k = 1, 2, ] T, Then the system (2.3) can be rewritten in the following vector form u (t) = f (u(t)) D (1) u(t) + ν D (2) u(t) f (u(t)) (a (1) α(t) + b (1) β(t) ) + ν ( a (2) α(t) + b (2) β(t) ), (2.4) u( ) = g, t, where stands for componentwise multiplication. This system is then integrated in time by the classical Runge Kutta method of order four, compare [6]. If the initial condition g(x) or the solution to (1.1) or (1.2) is discontinuous it may be necessary to use filters to stabilize the resulting scheme and to enhance the accuracy of numerical approximations. The use of filters for discontinuous problems is discussed in [39], [17], [23], [1], [15], [36], and 4
5 [16]. If the initial function g(x) is discontinuous we can obtain the filtered approximation M gm(x σ i ) = fijg(x σ j ), i = 0, 1,..., M, where F σ = [ f σ ij f σ ij = w j ] M M n=0 i,j=0 j=0 is the filter matrix defined by σ ( n ) T n (x i )T n (x j ). M γ n Here, T n (x) = cos(n arccos(x)), n = 0, 1,..., are Chebyshev polynomials, π, j = 0, M, 2M π, n = 0, M, w j = γ π n = M, j = 1, 2,..., M 1, π, n = 1, 2,..., M 1, 2 are the weights of the Chebyshev Gauss Lobatto quadrature formula and the normalization constants, respectively, and σ(η) is a filter of order p, i.e., the real function σ C p 1 (, ) which has the following properties: σ(η) = 0 for η > 1; σ(0) = 1; σ(1) = 0; and σ (m) (0) = σ (m) (1) = 0 for m = 1, 2,..., p 1, compare [39], [15], [16]. The frequently used filters are the Cesáro, raised cosine, sharpened raised cosine, Lanczos, and exponential cut off filter, see [39], [7], [15], [16]. The family of exponential filters of order p is defined by 1, η η c, σ(η) = ( exp α ( η η ) ) c p, η > η c, 1 η c where α is a measure of how strong the various modes should be filtered. Usually α is chosen so that σ(1) is of the order of the machine accuracy eps of the computer. This leads to α = ln(c eps) where C is a constant of moderate size, and we have chosen C = 1 in our implementation of this filter. We have found this family of filters to be quite effective in our numerical experiments reported in Section 4 and 6. We can also obtain matrices which combine the effect of filtering and differentiation of order k. These matrices D (k),σ = [ d (k),σ ] M ij are defined by d (k),σ ij = w j M n=0 σ ( n ) T n (x j )T n (k) (x i ), M γ n 5 i,j=0
6 compare again [15], [16]. The general system of differential equations resulting from the semidiscretization of (1.2) in space using filtering and differentiation matrices D (k),σ takes the form u (t) = f (u(t)) D (1),σ u(t) + ν D (2),σ u(t) f (u(t)) (a (1),σ α(t) + b (1),σ β(t) ) + ν ( a (2),σ α(t) + b (2),σ β(t) ), (2.5) u( ) = F σ g, ] M 1 t, where F σ = [ fij σ, and D (k),σ, a (k),σ, and b (k),σ are defined similarly i,j=1 as D (k), a (k), and b (k). Observe that for σ(η) 1 the system (2.5) reduces to (2.4). It is also possible to use different functions σ(η) for D (k),σ, k = 1, 2, and F σ, but it seems to be difficult to decide which strategy leads to the optimal results. We experimented with and compared many choices and results of some numerical experiments are presented in Sections 4 and 6. We have introduced small viscosity ν in (1.2) and (2.5) in the physical space. Guo, Ma and Tadmor [24] consider alternative approach, the so called spectral vanishing viscosity method, which is implemented directly on high modes of the computed solution in the spectral domain. 3 Location of discontinuities in the solution and Gegenbauer reconstruction After the approximation u h (x i, t end ) to the solution u(x i, t end ) of (1.1) or (1.2) is computed we can compute spectral coefficients a n (t end ), i.e., the coefficients of the expansion of u h (x i, t end ) in terms of Chebyshev polynomials T n (x) u h (x i, t end ) = M n=0 i = 0, 1,..., M, from the formula listed in [17] a n (t end ) = 2 M c n M j=0 a n (t end )T n ( x i ), (3.1) L 1 c j u h (x M j, t end ) cos πjn M, (3.2) where c 0 = c M = 2 and c n = 1 for n = 1, 2,..., M 1. By examination of these coefficients we can then find the number S and the location of discontinuities in the solution u to (2.5). This can be done in many different 6
7 ways and we have adopted the approach by Gottlieb, Lustman and Orszag described in [17]. In this approach the Chebyshev coefficients a n (t end ) (3.2) are fitted by an expression of the form S a n (t end ) = B s A n (X s ) s=1 for M/3 n 2M/3, which corresponds to the spectral representation with coefficients A n (X s ) of the sum of the Heaviside functions with discontinuities at X s and discontinuity jumps B s. As explained in [17] the values of X s are related to the eigenvalues of some summation operator and can be determined independently of B s. This leads to the exact determination of subintervals on which discontinuities are located and they are then assigned to the midpoints of these subintervals. As reported in [17] the number of shocks does not have to be specified in advance and this procedure works quite well for as many as seven shocks. Once X s are computed the values of B s can be found, if needed, by a least squares procedure. We refer to [17] for a complete description of this process. Another approach advocated by Gelb and Tadmor [13] in the context where Fourier spectral coefficients are given, is based on generalized conjugate partial sums whose convergence to the jump function is then accelerated by the so called concentration factors. This approach was generalized in [14] to general Jacobi polynomials which include Chebyshev polynomials as special case. This results in effective edge detectors, where both the location of discontinuities and discontinuity jumps are recovered. We refer again to [13], [14] for details. After the locations of discontinuities are detected we will use next Gegenbauer reconstruction as described in [19], [20], [21] to improve the accuracy of the numerical solution on the subintervals of [ L, L] where the solution to (1.1) or (1.2) is smooth. This reconstruction is based on ultraspherical or Genenbauer polynomials Cn(x), λ i.e., the polynomials which are orthonormal with respect to the inner product C λ k, C λ n = 1 h λ n 1 1 where, for λ > 0, h λ n is defined by (1 x 2 ) λ 1 2 C λ k (x)c λ n(x) dx h λ n = π 1 2 C λ Γ(λ n (1) ) Γ(λ)(n + λ), 7
8 and normalized so that C λ n(1) = Γ(n + 2λ) n! Γ(2λ). They can be computed from the three term recurrence relation (n + 1)C λ n+1(x) = 2(n + λ) x C λ n(x) (n + 2λ 1)C λ n 1(x), C λ 0 (x) = 1, C λ 1 (x) = 2λ x, n = 1, 2,..., compare [2], [8]. To illustrate this procedure let us assume that the solution to (1.1) or (1.2) is smooth on the subinterval [a, b] [ L, L], and, as in [19], define the local variable ξ by x = x(ξ) = ɛ ξ+δ, where ɛ = (b a)/2, δ = (b+a)/2, which transfoms the interval [ 1, 1] onto [a, b]. Denote by x i, i = 0, 1,..., M, the Chebyshev Gauss Lobatto points corresponding to [a, b] and compute the Chebyshev partial sum by ũ h ( x i, t end ) = M n=0 ( xi ) a n (t end )T n, (3.3) L i = 0, 1,..., M. We compute next the Gegenbauer series on a subinterval [a, b] by the formula m gm(x, λ t end ) = ĝɛ λ (l)c ( l λ x δ ), (3.4) l=0 ɛ where ĝ λ ɛ (l) are approximations to the Gegenbauer expansion coefficients g λ ɛ (l) based on a subinterval [a, b]. These coefficients are defined by g λ ɛ (l) = 1 h λ l 1 1 (1 ξ 2 ) λ 1 2 C λ l (ξ)ũ h (ɛ ξ + δ, t end ) dξ. The approximations to these coefficients are computed by evaluating the above integral by the Chebyshev Gauss Lobatto quadrature formula. This leads to ĝɛ λ (l) = 1 M w h λ j (1 x 2 j) λ Cl λ ( x j )ũ(ɛ x j + δ, t end ), (3.5) l j=0 where x j are the Chebyshev Gauss Lobatto points corresponding to the interval [ 1, 1] and w j are the corresponding weights. We expect that g λ m(x, t end ) 8
9 will be a better approximation to u(x, t end ) than u h (x, t end ) on the intervals of smoothness of u. This expectation is supported by the result of Gottlieb and Shu [19] who proved that Gegenbauer reconstruction m ĝɛ λ (l)cl λ (ξ) l=0 applied to the piecewise smooth function u(x) defined on the interval [ 1, 1] with λ = m = β ɛ M, β < 2/27, is uniformly convergent to u(x) on the interval of smoothness [a, b] of u as M. To be more precise, the resulting error can be bounded by max 1 ξ 1 u(ɛ ξ + δ) m l=0 ĝ λ ɛ (l)c λ l (ξ) A ( q ɛ M T ) + qr ɛ M, where ( ) 27β β q T = < 1, qr = 2 ( ) 27ɛ β < 1, 32ρ and A grows at most as M 1+2λ. Here, ρ is the distance from [a, b] to the nearest singularity of u(x) in the complex plane. The constants q T and q R correspond to truncation and regularization errors, see [19], [20], [21]. The implementation of Gegenbauer reconstruction is quite sensitive to round off errors and it is not clear how to choose m and λ to obtain optimal results. Gottlied and Shu [19] recommend the values and Guo, Ma, and Tadmor [24] choose m = 0.1 ɛ M, λ = 0.2 ɛ M, (3.6) m = λ = 0.05 M, but they did not attempt to optimize these parameters. To our knowledge the first attempt to optimize these parameters was made by Gelb [11] for Gegenbauer reconstruction where the Fourier coefficients are given for a smooth but non periodic function. This approach was further refined in [12] by taking into account the smoothness characteristics of the function. The determination of optimal parameters based on the Chebyshev spectral coefficients is discussed in [25]. 9
10 M=20 M=20 M=40 M= M=80 M=80 error M= x Figure 1: Errors of the Genenbauer reconstruction with m and λ defined by (3.6) for M = 20, 40, 80, and 160. To illustrate the effectiveness of Gegenbauer reconstruction we have plotted in Fig. 1 the error of this procedure for the function f(x) given by exp(x), x [ 1, 0], f(x) = sin(cos(x)) + 1, x (0, 1], with the parameters m and λ defined by (3.6) for M = 20, 40, 80, and 160. The performance of this procedure could be further improved by appropriate tuning of the parameters m and λ, compare again [11], [12], [25]. The graphs of errors versus M in double logarithmic scale which illustrate spectral convergence of Gegenbauer reconstruction method for functions with varying degrees of smoothness are presented in [25] and [12]. 10
11 4 Numerical experiments with pseudospectral methods and Gegenbauer reconstruction In this section we present the results of some numerical experiments for the Burgers equation u(x, t) + u(x, t) t x u(x, t) = 0, L x L, t [, t end ], u(x, ) = x, x < 2A, u(x, ) = 0, x > 2A, (4.1) u(±l, t) = ± L t, t [, t end ], where A > 0 is a parameter. The solution to this problem is the N wave x u(x, t) = t, x < 2At, 0, x > 2At, compare [40]. This solution is the limit as ν 0 of the solutions to the problem u(x, t) + u(x, t) t x 2 u(x, t) = ν u(x, t), (4.2) x2 L x L, t [, t end ], with appropriate initial and boundary conditions. The solution to (4.2) is given by u(x, t) = x ( t e x2/(4νt) ) 1 1 +, t e A/(2ν) 1 compare again [40]. The pseudospectral Chebyshev method applied to (4.2) leads to the system of ODEs u (t) = u(t) D (1),σ u(t) + ν D (2),σ u(t) u(t) (a (1),σ α(t) + b (1),σ β(t) ) + ν ( a (2),σ α(t) + b (2),σ β(t) ), (4.3) u( ) = F σ g, 11
12 x x M Figure 2: Numerical results for M = 64 (left graphs) and error versus M (right graph). t, where g correspond to initial function at the grid points x i and F σ, D (k),σ, a (k),σ, and b (k),σ are defined in Section 2. We have solved the problem (4.3) with = 1, t end = 2 and A = 1 by pseudospectral Chebyshev method (2.5). This system was integrated by the classical Runge Kutta method of order four. The efficiency of numerical computations could be increased somewhat by the use of total variation diminishing (TVD) or strong stability preserving (SSP) Runge Kutta methods (compare [22], [35]). However, time integration is not the main focus of this paper and we do not pursue these issues in more detail here. After integrating (2.5) in time, we then computed Chebyshev spectral coefficients (3.2) which were used to locate the discontinuities in the solution 12
13 and for Gegenbauer reconstruction as described in Section 3. In Fig. 2 (left graphs) we present the results for M = 64, ν = 1/M, time step t = 10 3, D (1),σ = D (1), D (2),σ = D (2), where we filter the initial condition using F σ corresponding to the exponential filter with p = 4 and η c = 2/M. The parameters used for Gegenbauer reconstruction were chosen by trial and error and they are m = 1 and λ = 4 for left and right subintervals and m = 2 and λ = 5 for the middle subinterval. In the left upper graph of this figure the exact solution is plotted by thick solid line, the numerical solution before postprocessing and the not a knot spline fit of it by black square and thin solid line and the numerical solution obtained by Gegenbauer reconstruction by white circle. In the left lower graph we have plotted the not a knot spline fit of the error before postprocessing by dashed line and the error after postprocessing by solid line. These global errors are defined as norms of the differences between the exact solution to (4.2) and the corresponding numerical solutions. We can see that Gegenbauer postprocessing is very effective in reconstructing the solution to (4.1) on the subintervals of smoothness. In Fig. 2 (right graph) we have plotted in double logarithmic scale error versus M (black circles connected by thin solid line) together with the spline fit (thick solid line) generated using program csaps.m from Matlab Spline Toolbox (Version 3.2.1). This graph illustrates spectral convergence of the overall numerical scheme. The parameters m and λ of Gegenbauer reconstruction were chosen again by trial and error for different values of M and different subintervals. The systematic approach to choosing these parameters is the subject of current research, compare [12]. 5 Waveform relaxation methods Waveform relaxation method is an iterative technique to compute successive approximations u (k) (t), k = 0, 1,..., to the solution u(t) of (2.5). To be more specific, we use previous iterate u (k) (t) in the argument of the function f and the next iterate u (k+1) (t) elsewhere. This has the effect of replacing the 13
14 nonlinear system (2.5) by a sequence of linear problems of the form d dt u(k+1) (t) = f (u (k) (t)) D (1),σ u (k+1) (t) + ν D (2),σ u (k+1) (t) f (u (k) (t)) (a (1),σ α(t) + b (1),σ β(t) ) + ν ( a (2),σ α(t) + b (2),σ β(t) ), (5.1) u (k+1) ( ) = F σ g, t, k = 0, 1,.... We can also consider more general schemes corresponding to appropriate splittings of matrices D (l),σ, i.e., D (l),σ = (l),σ (l),σ D 1 + D 2, l = 1, 2. For example, choosing Gauss Jacobi or block Gauss Jacobi splittings D (l),σ ) (l),σ (l),σ 1 = diag ( D 1,1,..., D 1,s, l = 1, 2, (l),σ (l),σ where D 1,j, j = 1, 2,..., s, are diagonal blocks of D 1 of dimension m j, m m s = M 1, leads to the sequence of linear problems d dt u(k+1) (t) = f (u (k) (1),σ (t)) ( D 1 u (k+1) (1),σ (t) + D 2 u (k) (t) ) (2),σ + ν ( D 1 u (k+1) (2),σ (t) + D 2 u (k) (t) ) (5.2) f (u (k) (t)) (a (1),σ α(t) + b (1),σ β(t) ) + ν ( a (2),σ α(t) + b (2),σ β(t) ), u (k+1) ( ) = F σ g, t, k = 0, 1,..., where u (0) (t) is an arbitrary starting function usually chosen as u (0) (t) = F σ (l),σ g, t. Observe that the form of the matrices D effectively decouples the resulting system (5.2) into s subsystems which can be assigned to different processors for efficient solution in a parallel computing environment. For any function u : [, T ] R M 1 define the supremum norm on the interval [, T ] by u [t0,t ] = sup{ u(t) : t [, T ]}, where is any norm in R M 1. To study boundedness of the sequence u (k) defined by (5.2) we assume that f (u) F (Q) for u Q, (5.3) 14
15 where F (Q) is a constant depending on Q. We also introduce the following notation: µ i (T ) = a (i),σ α + b (i),σ β [t0,t ], i = 1, 2, and C i = C i (Q, ν, σ) = F (Q) D (1),σ i + ν D (2),σ i, i = 1, 2, C 3 = C 3 (Q, ν, T ) = F (Q)µ 1 (T ) + νµ 2 (T ). We can demonstrate that the sequence u (k) defined by (5.2) is uniformly bounded for sufficiently small values of T. To be more precise we have the following theorem. Theorem 5.1 Assume that the function f satisfies (5.3) and the starting function u (0) satisfies u (0) [t0,t ] Q, for some constant Q 2 F σ g + 1. Then k = 0, 1,..., if T is chosen so that u (k) [t0,t ] Q, (5.4) (T )C 1 < 1 2 and (T )(C 2 Q + C 3 ) < 1 2. (5.5) Proof: The result will follow if we can show that (5.4) is satisfied for k + 1 assuming that it is satisfied for k. Integrating (5.2) from to T and taking norms in the resulting equation and using (5.3) we obtain or u (k+1) [t0,t ] F σ g + (T ) ( C 1 u (k+1) [t0,t ] + C 2 u (k) [t0,t ] + C 3 ), u (k+1) [t0,t ] F σ g + (T )(C 2 u (k) [t0,t ] + C 3 ) 1 (T )C 1. Since, by assumption, (5.4) is satisfied for k if T satisfies (5.5) it follows that u (k+1) [t0,t ] 2 ( F σ g + (T )(C 2 Q + C 3 ) ) 2 F σ g + 1 Q which proves (5.4) for k
16 To study the convergence of the sequence u (k) defined by (5.2) in the uniform norm we assume that f (u 1 ) f (u 2 ) L u 1 u 2, (5.6) for some constant L 0. We also introduce the notation ( C 4 = C 4 (L, Q, σ, T ) = L Q ( D (1),σ 1 + D (1),σ 2 ) ) + µ 1 (T ). Put e (k) (t) = u(t) u (k) (t), k = 0, 1,..., where u(t) is the solution to (2.5) and u (k) (t) are defined by (5.2). We have the following theorem. Theorem 5.2 Assume that the function f satisfies (5.3) and (5.6) and that T satisfies (5.5). Then e (k) (t) e C 1(t ) (C 5(t )) k E 0, (5.7) k! k = 0, 1,..., t [, T ], where C 5 = C 2 + C 4 and E 0 is some constant such that e (0) [t0,t ] E 0. Proof: Subtracting (5.2) from (2.5), then integrating from to t and taking norms on both sides of the resulting equation and using (5.3), (5.4) and (5.6) after some computations we obtain e (k+1) (t) C 1 k = 0, 1,..., t [, T ]. recurrence relation ɛ (0) (t) = e (0) (t), k = 0, 1,..., t [, T ]. inequalities that e (k+1) (s) ds + C 5 e (k) (s) ds, Define the sequence of functions ɛ (k) (t) by the ɛ (k+1) (t) = C 1 The equation (5.8) is equivalent to ɛ (k+1) (s)ds + C 5 ɛ (k) (s)ds, (5.8) Then it follows from the theory of differential e (k) (t) ɛ (k) (t). d dt ɛ(k+1) (t) C 1 ɛ (k+1) (t) = C 5 ɛ (k) (t), t [, T ], ɛ (k+1) ( ) = 0, 16
17 which has the solution e C 1(t ) ɛ (k+1) (t) = C 5 Define the sequence ɛ (k) (t) by e C 1(s ) ɛ (k) (s)ds. ɛ (k) (t) = e C 1(t ) ɛ (k) (t). Then ɛ (k) (t) satisfies the recurrence equation ɛ (k+1) (t) = C 5 k = 0, 1,..., which has the solution ɛ (k) (t) = C k 5 ɛ (k) (s)ds, (t s) k 1 (k 1)! ɛ(0) (s)ds. (5.9) To see that this is the case observe that this relation is satisfied for k = 1. Assuming that (5.9) is satisfied for k we obtain ɛ (k+1) (t) = C k+1 5 = C k+1 5 = C k+1 5 s ξ (t ξ) k k! (s ξ) k 1 ɛ (0) (ξ) dξ ds (k 1)! (s ξ) k 1 (k 1)! ɛ (0) (ξ) dξ ɛ (0) (ξ) ds dξ which proves (5.9) for k + 1. Since ɛ (0) (t) = e C 1(t ) ɛ (0) (t) E 0, (5.9) implies that ɛ (k) (t) C k 5 (t s) k 1 (k 1)! ds E 0 = (C 5(t )) k k! E 0, (5.10) k = 0, 1,..., t [, T ], which proves (5.7). The estimate (5.10) proves the superlinear convergence of the sequence u (k) to the solution u of (4.2). This estimate can be improved by using the relation ɛ (k) (t) = C5 k (t s) k 1 E 0 (k 1)! ec 1(t s) ds 17
18 instead of (5.9). This leads to the error bound e (k) (t) Ck 5 E 0 (k 1)! t0 0 s k 1 e C 1s ds, k = 0, 1,..., expressed in terms of the integral of the function s k 1 exp(c 1 s), which is a sharper estimate than (5.7). We can also show that the sequence u (k) is uniformly bounded on arbitrarily long intervals [, T ]. To do so we assume that Q > F σ g + 1 and use the exponential norm where α is chosen so that u α [,T ] := sup { u(t) e α(t ) : t [, T ] }, (5.11) α > max { C 1, C 3 /( F σ g + 1), (C 1 + C 2 )Q/(Q F σ g 1) }. We have the following theorem. Theorem 5.3 Assume that the function f satisfies (5.3) and that the starting function u (0) satisfies with Q > F σ g. Then for k = 0, 1,.... u (0) (t) Q, t [, T ] (5.12) u (k) α [,T ] Q (5.13) Proof: Note that (5.12) implies (5.13) with k = 0. Assume (5.13) for k. Integrating (5.2) from to t, then taking norms on both sides of the resulting equation and using (5.3) we obtain u (k+1) (t) F σ g + C 1 + C 2 u (k+1) (s) e α(s ) e α(s ) ds u (k) (s) e α(s ) e α(s ) ds + C 3 F σ g + C 1 u (k+1) α [,T ] + C 2 u (k) α [,T ] e α(s ) ds e α(s ) ds + C 3 (t ) F σ g + C 1 α u(k+1) α [,T ] e α(t ) + C 2 α u(k) α [,T ] e α(t ) + C 3 (t ). 18 ds
19 Hence, u (k+1) (t) e α(t ) ( F σ g + C 3 (t ) ) e α(t ) Observe that since α > C 3 ( F σ g + 1) we have + C 1 α u(k+1) α [,T ] + C 2 α u(k) α [,T ]. ( F σ g + C 3 (t ) ) e α (t ) F σ g + 1 and it follows that This leads to α u (k+1) α [,T ] α ( F σ g + 1) + C 1 u (k+1) α,t ] + C 2 Q. u (k+1) α [,T ] α ( F σ g + 1) + C 2 Q α C 1 Q. The last inequality follows from the definition of α. This completes the proof. Finally, we study the convergence of u (k) to the solution u of (2.5) in the exponential norm (5.11) on arbitrary interval [, T ]. Define the constants Q = Qe α(t ) C i = C i ( Q, (1),σ ν, σ) = F ( Q) D i + ν D (2),σ i, i = 1, 2, C 3 = C 3 ( Q, ν, T ) = F ( Q)µ 1 (T ) + νµ 2 (T ), C 4 = C 4 (L, Q, ( σ, T ) = L Q ( D (1),σ 1 + D (1),σ 2 ) ) + µ 1 (T ). We have the following theorem. Theorem 5.4 Assume that the function f satisfies (5.3) and (5.6). Then e (k) α [,t] e C 1 (t ) ( C 5 (t )) k Ẽ 0, (5.14) k! k = 0, 1,..., t [, T ], where C 5 = C 2 + C 4 and Ẽ0 is some constant such that e (0) α [,T ] Ẽ0. 19
20 Proof: Observe that u k (t) α [,T ] Q implies that uk (t) Q. Subtracting (5.2) from (2.5) and using (5.13) and (5.3) with Q replaced by Q after some computations we obtain e (k+1) (t) C 1 + C 5 k = 0, 1,..., t [, T ]. Hence, e (k+1) α [,t] C 1 C 1 e α(t ) + C 5 e α(t ) e (k+1) (s) e α(s ) e α(s ) ds e (k) (s) e α(s ) e α(s ) ds e (k+1) α [,s] ds e (k) α [,s] ds, e (k+1) α [,s] ds + C 5 e (k) α [,s] ds. The theorem now follows using exactly the same arguments as in the proof of Theorem 5.2. Similarly as before the error bound in Theorem 5.4 can be replaced by e (k) α [,t] C 5k Ẽ0 t0 s k 1 e C 1 s ds, (k 1)! 0 k = 0, 1,..., expressed in terms of the integral of the function s k 1 exp( C 1 s), which is a sharper estimate than (5.14). 6 Numerical experiments with waveform relaxation methods and Gegenbauer reconstruction In this section we will use again Burgers equations (4.1) and (4.2) as test problems for waveform relaxation. The sequence of linear systems of ODEs 20
21 corresponding to (5.2) takes now the form d dt u(k+1) (t) = u (k) (1),σ (t) ( D 1 u (k+1) (1),σ (t) + D 2 u (k) (t) ) (2),σ + ν ( D 1 u (k+1) (2),σ (t) + D 2 u (k) (t) ) u (k) (t) (a (1),σ α(t) + b (1),σ β(t) ) + ν ( a (2),σ α(t) + b (2),σ β(t) ), (6.1) u (k+1) ( ) = F σ g, t, k = 1, 2,.... To integrate these systems in time we apply the implicit x x M Figure 3: Numerical results after three iterations of Gauss Jacobi waveform relaxation method with M = 64. Runge Kutta methods of order four based on the Gauss Legendre quadra- 21
22 ture formula [6]. This method can be represented by Butcher tableaux (6.2) We have solved problem (5.2) with = 1, t end = 2 and A = 1 by waveform relaxation method corresponding to the Gauss Jacobi splitting of the matrices D (l),σ, l = 1, 2. We then computed Chebyshev spectral coefficients (3.2) which were used to locate the discontinuities in the solution and for Gegenbauer reconstruction as described in Section 3. In Fig. 3 (left graphs) we present the results of numerical experiments after three iterations, i.e., solving the system (6.1) for k = 0, 1 and 2, for M = 64, ν = 1/M and time step t = 10 2 which is an order of magnitude larger than that used in Section 4 (l),σ (l),σ for explicit Runge Kutta of order 4. The matrices D 1, D 2 and F σ correspond to the exponential filter with p = 4 and η c = 2/M. The parameters used for Gegenbauer reconstruction are m = 1 and λ = 3. As in Section 4 on the upper graph of this figure the exact solution is plotted by thick solid line, the numerical solution before postprocessing and the not a knot spline fit of it by black square and thin solid line and the numerical solution obtained by Gegenbauer reconstruction by white circle. In the lower graph we have plotted the not a knot spline fit of the error before postprocessing by dashed line and the error after postprocessing by solid line. These errors are defined as in Section 4. In Fig. 3 (right graph) we have plotted in double logarithmic scale error versus M (black circles connected by thin solid line) together with the spline fit (thick solid line) generated using program csaps.m from Matlab Spline Toolbox (Version 3.2.1). This graph illustrates spectral convergence of the overall numerical scheme. As in Section 4 the parameters m and λ of Gegenbauer reconstruction were chosen again by trial and error for different values of M and different subintervals. Similar results are obtained if BDF3, the backward differentiation method of order three, is used with the same time step for integrating the system (6.1). 22
23 7 Concluding remarks We investigate Chebyshev spectral collocation and waveform relaxation methods for nonlinear conservation laws. The location of discontinuities in the solution is determined by fitting the spectral coefficients corresponding to the numerical solution by spectral representation of a sum of Heaviside functions, and the numerical solution is then enhanced on the intervals of smoothness by Gegenbauer reconstruction. The systems of differential equations resulting from application of spectral collocation methods are solved by explicit Runge Kutta method of order four which require rather small time steps for stable integration. The main advantage of waveform relaxation consists of the fact that we can replace these nonlinear systems of ODEs by a sequence of linear problems which can then be effectively integrated by A stable implicit Runge Kutta methods or A(α) stable backward differentiation methods. This allows for much larger time steps than those used for explicit methods. This is confirmed by numerical experiments presented in Sections 4 and 6. Another advantage of waveform relaxation is that it can be applied in a parallel computing environment. Future work will address the numerical solution of nonlinear conservation laws and Gegenbauer reconstruction in two or three space dimensions, and the implementation of waveform relaxation methods with adaptive window control strategy. Some progress for Gegenbauer reconstruction in two dimension for rectangular domains using tensor product approach has been reported in [30]. The theoretical background for adaptive window strategy was formulated in [5]. Acknowledgements. The authors wish to express their gratidude to the anonymous referees for their useful comments. References [1] S. Abarbanel, D. Gottlieb and E. Tadmor, Spectral methods for discontinuous problems. In: Numerical Methods for Fluid Dynamics II. Proceedings of the 1985 Conference on Numerical Methods in Fluid Dynamics (K.W. Morton and M.J. Baines, eds.), pp , Clarendon Press, Oxford
24 [2] M. Abramowitz and I.S. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Washington [3] K. Burrage, Parallel and Sequential Methods for Ordinary Differential Equations, Oxford University Press, Oxford [4] K. Burrage, Z. Jackiewicz, S.P. Nørsett and R. Renaut, Preconditioning waveform relaxation iterations for differential systems, BIT 36(1996), [5] K. Burrage, Z. Jackiewicz and B.D. Welfert, Spectral approximation of time windows in the solution of dissipative linear differential equations, submitted to IMA J. Numer. Anal. [6] J.C. Butcher, The Numerical Analysis of Ordinary Differential Equations. Runge Kutta and General Linear Methods. John Wiley & Sons, Chichester, New York [7] C. Canuto, M.Y. Hussaini, A. Quarteroni and T.A. Zang, Spectral Methods in Fluid Mechanics, Springer Verlag, New York, Berlin, Heidelberg [8] P.J. Davis, Methods of Numerical Integration, Academic Press, Ontario, San Diego, New York [9] K. Dekker and J.G. Verwer, Stability of Runge Kutta Methods for Stiff Nonlinear Differential Equations, North Holland, Amsterdam, New York, Oxford [10] B. Fornberg, A Practical Guide to Pseudospectral Methods, Cambridge University Press, Cambridge [11] A. Gelb, On the reduction of round off error for the Gegenbauer reconstruction method, J. Sci. Comput. 20(2004), [12] A. Gelb and Z. Jackiewicz, Determining analyticity for parameter optimization of the Gegenbauer reconstruction method, submitted to SIAM J. Sci. Comput. [13] A. Gelb and E. Tadmor, Detection of edges in spectral data, Appl. Comput. Harmon. Anal. 7(1999),
25 [14] A. Gelb and E. Tadmor, Detection of edges in spectral data. II. Nonlinear enhancement, SIAM J. Numer. Anal. 38(2000), [15] D. Gottlieb and J.S. Hesthaven, Spectral methods for hyperbolic problems, J. Comput. Appl. Math. 128(2001), [16] D. Gottlieb and J.S. Hesthaven, Spectral Approximation of Partial Differential Equations, manuscript. [17] D. Gottlieb, L. Lustman and S.A. Orszag, Spectral calculations of one dimensional invisid compressible flows, SIAM J. Sci. Stat. Comput. 2(1981), [18] D. Gottlieb and S.A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications. Society for Industrial and Applied Mathematics, Philadelphia [19] D. Gottlieb and C.W. Shu, On the Gibbs phenomenon IV: Recovering exponential accuracy in a subinterval from a Gegenbauer partial sum of a piecewise analytic function. Math. Comp. 64(1995), [20] D. Gottlieb and C.W. Shu, On the Gibbs phenomenon and its resolution, SIAM Rev. 39(1997), [21] D. Gottlieb and C.W. Shu, A general theory for the resolution of the Gibbs phenomenon, Atti dei Convegni Lincei 147(1998), [22] D. Gottlieb and C.W. Shu, Total variation diminishing Runge Kutta schemes, Math. Comp. 67(1998), [23] D. Gottlieb and E. Tadmor, Recovering pointwise values of discontinuous data within spectral accuracy. In: Progress and Supercomputing in Fluid Dynamics. Proceedings of a 1984 U.S. Israel Workshop (E.M. Murman and S.S. Abarbanel, eds.), pp , Birkhauser, Boston [24] B.Y. Guo, H.P. Ma and E. Tadmor, Spectral vanishing viscosity method for nonlinear conservation laws, SIAM J. Numer. Anal. 39(2001),
26 [25] Z. Jackiewicz, Determination of optimal parameters for the Chebyshev Gegenbauer reconstruction method, SIAM J. Sci. Comput. 25(2003/2004), [26] Z. Jackiewicz, M. Kwapisz and E. Lo, Waveform relaxation methods for functional differential systems of neutral type, J. Math. Anal. Appl. 207(1997), [27] Z. Jackiewicz, B. Owren and B.D. Welfert, Pseudospectra of waveform relaxation operators, Comp. & Math. Appls 36(1998), [28] Z. Jackiewicz and B.D. Welfert, Stability of Gauss Radau pseudospectral approximations of the one dimensional wave equation, J. Sci. Comput 18(2003), [29] J. Janssen, Acceleration of waveform relaxation methods for linear ordinary and partial differential equations, Ph.D. Thesis, Department of Computer Science, Katholieke Universiteit Leuven, [30] J.H. Jung and B.D. Shizgal, Inverse polynomial reconstruction of two dimensional Fourier images, submitted to J. Sci. Comput. [31] B. Leimkuhler, Estimating waveform relaxation convergence, SIAM J. Sci. Comput. 14(1993), [32] Y. Maday, S.M. Ould Kaber and E. Tadmor, Legendre pseudospectral viscosity method for nonlinear conservation laws, SIAM J. Numer. Anal. 30(1993), [33] U. Miekkala and O. Nevanlinna, Convergence of dynamic iteration methods for initial value problems, SIAM J. Sci. Stat. Comput. 8(1987), [34] O. Nevanlinna, Remarks on Picard Lindelöf iteration, Part I, BIT 29(1989), [35] R. Spiteri and S. Ruuth, A new class of optimal high order strong stability preserving time discretization methods, SIAM J. Numer. Anal. 40(2002),
27 [36] E. Tadmor and J. Tanner, Adaptive mollifiers for high resolution recovery of piecewise smooth data from its spectral information, Found. Comput. Math. 2(2002), [37] L.N. Trefethen, Pseudospectra of linear operators, SIAM Rev. 39(1997), [38] L.N. Trefethen, Spectral Methods in Matlab. Society for Industrial and Applied Mathematics, Philadelphia [39] H. Vandeven, Family of spectral filters for discontinuous problems, J. Sci. Comput. 6(1991), [40] G.B. Whitman, Linear and Nonlinear Waves, John Wley & Sons, New York [41] B. Zubik Kowal, Chebyshev pseudospectral method and waveform relaxation for differential and differential functional parabolic equations, Appl. Numer. Math. 34(2000),
A Wavelet-based Method for Overcoming the Gibbs Phenomenon
AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 A Wavelet-based Method for Overcoming the Gibbs Phenomenon NATANIEL GREENE Department of Mathematics Computer
More information(4) k = 2N + I. A Numerical Study on the Accuracy of Fourier Spectral Methods Applied to the Nonlinear Burgers Equation.
A Numerical Study on the Accuracy of Fourier Spectral Methods Applied to the Nonlinear Burgers Equation Chi-Wang Shu* Peter S. Wong* Abstract It is well known that spectral method can achieve exponential
More informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
More informationEnhanced spectral viscosity approximations for conservation laws
Applied Numerical Mathematics 33 (000) 3 1 Enhanced spectral viscosity approximations for conservation laws Anne Gelb a,1, Eitan Tadmor b, a Department of Mathematics, Arizona State University, Tempe,
More informationStrong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods
Strong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods L. Ferracina and M. N. Spijker 2007, June 4 Abstract. This paper deals with the numerical solution of initial value problems, for systems
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationFOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE BOUSSINESQ CONVECTION PROBLEM
MATHEMATICS OF COMPUTATION Volume 7, Number 9, July 7, Pages 7 9 S -78(7)8- Article electronically published on February, 7 FOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE
More informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
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 informationSuper Viscosity and Spectral Approximations of Nonlinear Conservation Laws
umerical Methods for Fluid Dynamics. IV M. J. Baines and K. W. Morton, eds. Clarendon Press, Oxford, 1993, pp. 69-82 Super Viscosity and Spectral Approximations of onlinear Conservation Laws EITA TADMOR
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationStepsize Restrictions for Boundedness and Monotonicity of Multistep Methods
Stepsize Restrictions for Boundedness and Monotonicity of Multistep Methods W. Hundsdorfer, A. Mozartova, M.N. Spijker Abstract In this paper nonlinear monotonicity and boundedness properties are analyzed
More information1. Numerical Analysis of Spectral Methods/Theory and Applications, with S. Orszag, CBMS-SIAM No. 26, 1977, 170 pages.
DAVID GOTTLIEB Ford Foundation Professor Professor of Applied Mathematics Brown University Publications Updated February 2007 Books 1. Numerical Analysis of Spectral Methods/Theory and Applications, with
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationSpectral Methods for Hyperbolic Problems
Spectral Methods for Hyperbolic Problems J. S. Hesthaven EPFL-SB-MATHICSE-MCSS, Ecole Polytechnique Fédérale de Lausanne CH-1015 Lausanne, Switzerland E-mail: Jan.Hesthaven@epfl.ch We review spectral methods
More informationThe Resolution of the Gibbs Phenomenon for Fourier Spectral Methods
The Resolution of the Gibbs Phenomenon for Fourier Spectral Methods Anne Gelb Sigal Gottlieb December 5, 6 Introduction Fourier spectral methods have emerged as powerful computational techniques for the
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More information* ICASE G 11t w. t )t Operated by Universities Space Research Association
[ NASA Contractor Report 194959 0 ICASE Report No. 94-67 AD-A284 063 * ICASE A NOTE ON THE ACCURACY OF SPECTRAL METHOD APPLIED TO NONLINEAR CONSERVATION LAWS Chi-Wang Shu Peter S. Wong ( ELECTED biic SEP
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More informationTHE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands
THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay
More informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationPAPER Resolution of the Gibbs Phenomenon for Fractional Fourier Series
57 PAPER Resolution of the Gibbs Phenomenon for Fractional Fourier Series Hongqing ZHU a), Member, Meiyu DING, and Daqi GAO, Nonmembers SUMMARY The nth partial sums of a classical Fourier series have large
More informationA Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms
Universal Journal of Integral Equations 4 (2016), 54-62 www.papersciences.com A Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms N. H. S. Haidar CRAMS: Center for
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationEvaluation of Chebyshev pseudospectral methods for third order differential equations
Numerical Algorithms 16 (1997) 255 281 255 Evaluation of Chebyshev pseudospectral methods for third order differential equations Rosemary Renaut a and Yi Su b a Department of Mathematics, Arizona State
More information(Received November , accepted May )
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.10(010) No.,pp.164-178 Spectral Solutions for Some Hyperbolic Partial Differential Equations by the Ultraspherical
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
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 informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG
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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationSingular Padé-Chebyshev Reconstruction: The Case When Jump Locations are Unknown
Singular Padé-Chebyshev Reconstruction: The Case When Jump Locations are Unknown Arnel L. Tampos Jose Ernie C. Lope Jan S. Hesthaven Abstract We present a singularity-based approach to resolve the Gibbs
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationSolution of Parabolic Partial Differential Equations
Solution of Parabolic Partial Differential Equations in Complex Geometries by a Modified Fourier Collocation Method Knut S. Eckhoff* Carl Erik Wasberg* Abstract The heat equation with Dirichlet boundary
More informationApplication of the relaxat ion met hod to model hydraulic jumps
Application of the relaxat ion met hod to model hydraulic jumps P. J. Montgomery Mathematics and Computer Science Program, University of Northern British Columbia, Prince George, Canada. Abstract A finite
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
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 informationApplied and Computational Harmonic Analysis 7, (1999) Article ID acha , available online at
Applied and Computational Harmonic Analysis 7, 0 35 (999) Article ID acha.999.06, available online at http://www.idealibrary.com on Detection of Edges in Spectral Data Anne Gelb Department of Mathematics,
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationOptimal Runge Kutta Methods for First Order Pseudospectral Operators
Journal of Computational Physics 5, 404 49 (999) Article ID jcph.999.660, available online at http://www.idealibrary.com on Optimal Runge Kutta Methods for First Order Pseudospectral Operators J. L. Mead
More informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationSpectral methods for compressible reactive flows
C. R. Mecanique 333 (2005) 3 16 http://france.elsevier.com/direct/cras2b/ High-Order Methods for the Numerical Simulation of Vortical and Turbulent Flows Spectral methods for compressible reactive flows
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationOptimal Implicit Strong Stability Preserving Runge Kutta Methods
Optimal Implicit Strong Stability Preserving Runge Kutta Methods David I. Ketcheson, Colin B. Macdonald, Sigal Gottlieb. August 3, 2007 Abstract Strong stability preserving (SSP) time discretizations were
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
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 informationParallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationStrong stability preserving high order time discretizations.
SIGAL GOTTLIEB Mathematics Department 285 Old Westport Road North Dartmouth, MA 02747 sgottlieb@umassd.edu 322 Cole Avenue Providence, RI 02906 Phone: (401) 751-9416 sigalgottlieb@yahoo.com Current Research
More informationWojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS
Opuscula Mathematica Vol. 30 No. 2 2010 http://dx.doi.org/10.7494/opmath.2010.30.2.133 Wojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS Abstract. We present
More informationarxiv: v1 [math.na] 4 Nov 2015
ON WELL-CONDITIONED SPECTRAL COLLOCATION AND SPECTRAL METHODS BY THE INTEGRAL REFORMULATION KUI DU arxiv:529v [mathna 4 Nov 25 Abstract Well-conditioned spectral collocation and spectral methods have recently
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationSpectral Methods and Inverse Problems
Spectral Methods and Inverse Problems Omid Khanmohamadi Department of Mathematics Florida State University Outline Outline 1 Fourier Spectral Methods Fourier Transforms Trigonometric Polynomial Interpolants
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationA Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationSpectral method for the unsteady incompressible Navier Stokes equations in gauge formulation
Report no. 04/09 Spectral method for the unsteady incompressible Navier Stokes equations in gauge formulation T. W. Tee I. J. Sobey A spectral method which uses a gauge method, as opposed to a projection
More informationWe consider in this paper Legendre and Chebyshev approximations of the linear hyperbolic
LEGEDRE AD CHEBYSHEV DUAL-PETROV-GALERKI METHODS FOR HYPERBOLIC EQUATIOS JIE SHE & LI-LIA WAG 2 Abstract. A Legendre and Chebyshev dual-petrov-galerkin method for hyperbolic equations is introduced and
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationRATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages 72 83 c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
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 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 informationAdaptive Edge Detectors for Piecewise Smooth Data Based on the minmod Limiter
Journal of Scientific Computing, Vol. 8, os. /3, September 6 ( 6) DOI:.7/s95-988 Adaptive Edge Detectors for Piecewise Smooth Data Based on the minmod Limiter A. Gelb and E. Tadmor Received April 9, 5;
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationStarting Methods for Two-Step Runge Kutta Methods of Stage-Order 3 and Order 6
Cambridge International Science Publishing Cambridge CB1 6AZ Great Britain Journal of Computational Methods in Sciences and Engineering vol. 2, no. 3, 2, pp. 1 3 ISSN 1472 7978 Starting Methods for Two-Step
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More information