c 2003 Society for Industrial and Applied Mathematics
|
|
- Louisa Walton
- 5 years ago
- Views:
Transcription
1 SIAM J. NUMER. ANAL. Vol. 41, No. 6, pp c 3 Society for Industrial and Applied Mathematics CONVERGENCE ANALYSIS OF SPECTRAL COLLOCATION METHODS FOR A SINGULAR DIFFERENTIAL EQUATION WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN Abstract. Solutions of partial differential equations with coordinate singularities often have special behavior near the singularities, which forces them to be smooth. Special treatment for these coordinate singularities is necessary in spectral approximations in order to avoid degradation of accuracy and efficiency. It has been observed numerically in the past that, for a scheme to attain high accuracy, it is unnecessary to impose all the pole conditions, the constraints representing the special solution behavior near singularities. In this paper we provide a theoretical justification for this observation. Specifically, we consider an existingapproach, which uses a pole condition as the boundary condition at a singularity and solves the reformulated boundary value problem with a commonly used Gauss Lobatto collocation scheme. Spectral convergence of the Legene and Chebyshev collocation methods is obtained for a singular differential equation arising from polar and cylinical geometries. Key words. coordinate singularity, convergence, spectral collocation method AMS subject classifications. 65N35, 65N1, 65L1 DOI /S Introduction. Physical problems in polar, cylinical, or spherical geometries often give rise to mathematical models involving singular partial differential equations (PDEs) with smooth solutions. Acommon feature of these PDEs is that their solutions have special behavior near coordinate singularities, which forces the solutions to be smooth. For the spectral solution of this type of equation, special treatment for the coordinate singularities is needed, since a traditional spectral scheme either does not fully capture the special solution behavior or is ill-suited to fast transform techniques; e.g., see [6, 7, 11]. Anumber of spectral approaches have been developed in the past in attempts to capture the solution behavior near coordinate singularities. They include those expanding the solution in specially designed basis functions, such as spherical harmonics, parity-modified Fourier series, modified Robert functions, and eigenfunctions of singular Sturm Liouville problems [6, 9, 11,, 1,, 7]; approaches using inherent symmetries of the solution [9, 1]; and methods using pole conditions (i.e., compatibility conditions at the center of polar coordinates) as boundary conditions in the collocation context [1] and the Galerkin context [4, 5]. Many of these approaches have been successfully applied to steady state and time dependent problems including the Navier Stokes equations; e.g., see [13, 14,, 3, 6, 8]. Received by the editors August 1, ; accepted for publication (in revised form) June 11, 3; published electronically December 17, 3. This work was supported in part by the NSF (U.S.A.) through grant DMS 744, the Research Grants Council of the Hong Kong Special Administrative Region, China, through grant CityU-1141/1P, and the Special Funds for State Major Basic Research Projects of China through grant G Department of Mathematics, University of Kansas, Lawrence, KS 6645 (huang@math.ukans. edu). Department of Mathematics, Shanghai University, Shanghai 436, China (hpma@mail.shu. edu.cn). Department of Mathematics, City University of HongKong, Kowloon, HongKong, China (maweiw@math.cityu.edu.hk). 333
2 334 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN The solution behavior of PDEs near coordinate singularities can be described by an infinite number of pole conditions derived either from the underlying differential equation, by assuming some kind of smoothness of the solution, or more generally from the analyticity of the solution at singularities. Most of the existing methods are developed more or less to accommodate this behavior. However, it is observed numerically first by Orszag [1] and then by many other researchers (e.g., [6]) that it is unnecessary to impose all of the pole conditions in order for a numerical scheme to attain high accuracy. In fact, Huang and Sloan [1], using some of these pole conditions as the boundary conditions at the coordinate singularities, show numerically that the spectral collocation approximation of the Helmholtz equation on the unit disk has a spectral convergence rate. Three Legene-type pseudospectral schemes and their convergence analysis have been developed for axisymmetric domains by Bernardi, Dauge, and Maday in their recent book [1]. The basic idea behind these schemes is to incorporate the natural measure r of the coordinate singularity into the quaature formula defining the spectral approximation. In the radial direction the formula reads as (1.1) v(r)r = j= N v(r j ) j v P, where P is a polynomial space, r j =(1+ρ j )/, and the j s are the corresponding weights. Three sets of ρ j s and P are chosen, one for each scheme: Method A(Gauss Radau): P = P N, ρ N =1, ρ j ( j N 1) are the roots of L N+1 (ρ); Method B (Gauss Lobatto): P = P N 1, ρ = 1, ρ N =1, ρ j (1 j N 1) are the roots of ( ) LN (ρ)+l N+1 (ρ) ; 1+ρ Method C (Gauss Radau): P = P N 1, ρ N =1, ρ j ( j N 1) are the roots of L N+1 (ρ) L N (ρ), where L N is the Legene polynomial of degree N. The pseudospectral approximations are then defined through the boundary condition(s) and the Galerkin formulation of the underlying problem in the discrete inner product ((u, v)) N N j= u(r j)v(r j ) j induced from the quaature formula (1.1). It is noted that the nodes used in these schemes are different from those in a traditional (unweighted) spectral collocation method. Moreover, among these three schemes, only Method C is equivalent to a collocation system. Furthermore, the authors of the book suggest that two boundary conditions u() = (which is a pole condition, cf. (.6)) and u(1) = g be used for a reduced equation (see (.1) (.)) with n. Thus, only Method B, which uses the Gauss Lobatto nodes but cannot be interpreted as a collocation scheme, can be applied to the case n. The objective of this paper is to provide a theoretical justification for the method developed in [1], which uses a pole condition as the boundary condition at the coordinate singularity and solves the reformulated boundary value problem with a Gauss Lobatto collocation scheme. The method corresponds to the standard quaature formula (1.) v(r)(r) = j= N w j v(r j ) for v P N 1,
3 SPECTRAL COLLOCATION FOR A SINGULAR DE 335 where (r) is the weight function and r j =(1+ρ j )/. For example, (r) =1and ρ j ( j N) are the roots of (1 ρ )L N (ρ) for the Legene collocation method, and (r) =(r r ) 1/ and ρ j = cos π(n j) N ( j N) for the Chebyshev collocation method. We emphasize that the method of [1], which is no more than a traditional collocation method, is different from those considered by Bernardi, Dauge, and Maday [1]. The method shares with many existing methods the common feature of explicitly using pole conditions, and has been successfully applied to practical problems including the Navier Stokes equations; e.g., see [13, 14]. Our analysis is given for both the Legene and Chebyshev schemes. AChebyshev collocation scheme is often desirable in practical computation because the fast Fourier transformation (FFT) can be utilized. We find that for the current situation with coordinate singularities the corresponding bilinear form lacks the coercive property which is often crucial to the convergence analysis of a Chebyshev scheme. Because of this, the error estimate of the Chebyshev scheme is obtained in the weighted energy norm En, with (r) being the Chebyshev weight function for the reduced equation with n>, but in the unweighted norm En for the case n =. An outline of this paper is as follows. The method of [1] is briefly described in section. The convergence analysis of the Legene and Chebyshev methods is given in section 3. In section 4 we present numerical results to verify the theoretical findings. Finally, section 5 contains conclusions and further comments.. Pole conditions and spectral collocation approximation. In this section we briefly describe the spectral collocation method of [1] for a model problem (.1) d u 1 du r + n r u = f, (.) u(1) = g, <r<1, where n is a given integer. This problem is obtained using separation of variables for the Poisson equation on the unit disk..1. Pole conditions. Equation (.1) has a coordinate singularity at r =. Assume that both f and u are sufficiently smooth. ATaylor series expansion of u about r = yields the pole conditions ( ) 1 (.3) O r : n u()=, ( ) 1 (.4) O :(n 1) du r ()=, ( ) n d (.5) O (1) : u () = f(), These conditions contain full information about the solution behavior near r =. It was observed first by Orszag [1] and later by many other researchers (see [6]) that it is unnecessary to impose all of these pole conditions in order for a numerical scheme to obtain high accuracy. In fact, using one constraint (.6) { u() = for n, du () = for n =
4 336 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN as the boundary condition at r =, Huang and Sloan [1] obtain spectrally accurate solutions; also see [4, 5] for the spectral Galerkin approximation. Once a boundary condition has been defined at r =, it is straightforward to apply a traditional spectral collocation scheme to the singular problem (.1) and (.)... Legene and Chebyshev collocation approximations. Hereafter, the weight functions (r) =1and(r) =(r r ) 1/ will be associated with the Legene and Chebyshev methods, respectively. For simplicity, we use subscript in common notation for both methods and for those which apply only to the Chebyshev method, and suppress the subscript for the Legene method. For a given integer N>, let {ρ j, } N j= be a set of Gauss Lobatto points associated with the weight function (r). Define (.7) r j, = 1+ρ j,, j =, 1,...,N. The solution u(r) is approximated by (.8) N u N (r) = u j, l j, (r), j= where u j, denotes the approximation of u(r j, ) and l j, (r) is the Lagrangian interpolation polynomial (.9) l j, (r) = N i= i j r r i, r j, r i,. Acollocation approximation to (.1), (.), and (.6) is then defined by the collocation equations (.1) (.11) (.1) d u N (r j,) 1 du N r j, (r j,)+ n rj, u N (r j, )=f(r j, ), j =1,...,N 1, u N (1) = g, { u N () = for n, du N () = for n =. Recall that the transformed Gauss Lobatto quaature rule satisfies (.13) v(r)(r) = j= N w j, v(r j, ) for v P N 1, where the w j, s are the corresponding weights and P N 1 is the space of real polynomials (in r) of degree no more than N 1. The associated interpolation operator I N : C[, 1] P N is defined as (.14) I N v P N : (I N v)(r j, )=v(r j, ), j =, 1,...,N.
5 We use the notation SPECTRAL COLLOCATION FOR A SINGULAR DE 337 (.15) (.16) (.17) u, v = ( m u m, = d k u k u, v,n = uv, u = u, u 1/, k= N j= ) 1/ for u H m, w j, u(r j, )v(r j, ), u,n = u, u 1/,N, where H m (m ) is a (weighted) Sobolev space on [, 1]. For the Legene collocation scheme, the set of Legene Gauss Lobatto points is defined by ρ = 1, ρ N = 1, and ρ j (j =1,,...,N 1) being the roots of L N, the first derivative of the Legene polynomial L N of degree N. We have r(1 r) L N l j (r) = (r) N(N +1) (r r j ) L N (r j ), w j = 1 N(N +1) 1 L N (r j), with L N (r) being the transformed Legene polynomial L N (r 1). For the Chebyshev approximation, the set of Gauss Lobatto points is defined by ρ, = 1, ρ N, = 1, and ρ j, (j =1,,...,N 1) being the roots of T N, the derivative of the Chebyshev polynomial T N of degree N. We have l j, (r) =( 1) j+1 r(1 r) T N (r) c j N, w j, = π (r r j, ) c j N, where T k (r) is the transformed Chebyshev polynomial T k (r 1) and {, j =,N, c j = 1, j =1,...,N Convergence analysis Preliminary approximation results. To start with, we introduce some preliminary results. Hereafter, C is used to denote the generic constant. We shall assume that N m (the smoothness order of functions); otherwise, the estimates given below, especially those involving seminorms, will not be true. Lemma 3.1. Let P N denote the Legene (or Chebyshev) truncated operator; i.e., P N v is the truncated Legene (or Chebyshev) series of v. Then, for m and for (r) =1(the Legene case) or (r) =(r r ) 1/ (the Chebyshev case), v P N v CN m (r r ) m/ dm v (3.1) m v H m. Lemma 3.. For (r) =1(the Legene case) or (r) =(r r ) 1/ (the Chebyshev case) and for m 1, (r r ) 1/ (v I N v) + N 1 v I N v 1, CN m (r r ) m 1 d m v (3.) m v H m.
6 338 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN Lemma 3.3. Let (r) =1(for the Legene case) or (r) =(r r ) 1/ (for the Chebyshev case). There exists a positive constant C, independent of N and M, such that for all φ P M with M being any nonnegative integer, ( φ,n C 1+ M ) (3.3) φ, N (3.4) φ φ,n C φ. The interested reader is referred to [, 3, 4, 7, 16, 18, 19] for the proofs of these Lemmas. Lemmas 3.1 and 3. are the improvements of existing results in terms of the weight and can be obtained by the method in the aforementioned references. 3.. Convergence analysis of the Legene method. We now proceed to the convergence analysis for the Legene approximation (.1) (.1). Let φ be an arbitrary polynomial in P N satisfying { φ(1) = φ() = if n, (3.5) φ(1) = if n =. Multiplying (.1) by r j w j φ(r j ) and summing over the range of j from 1 to N 1, we have (3.6) N 1 j=1 d w j φ(r j ) [ r u N ] j (r j) dun (r j)+ n N 1 u N (r j ) = r j It is not difficult to see from (3.5) that d w N φ(r N ) [ r u N ] (3.7) N (r N) dun (r N)+ n u N (r N ) r N = r N w N φ(r N )f(r N ) =. j=1 r j w j φ(r j )f(r j ). Noticing that u N (r)/r is a polynomial of degree not greater than N 1, that φ(r )= (see (3.5)) when n, and that (du N /)(r ) = when n =,wehave d w φ(r ) [ r u N ] (3.8) (r ) dun (r )+ n u N (r ) r = r w φ(r )f(r ) =. Thus, (3.6) (3.8) imply that r d u N (3.9) dun + n r un,φ N = rf, φ N. Since rφ(d u N / ), φ(du N /), and n u N φ/r are in P N 1, (.13) and (3.9) lead to r d u N (3.1) dun + n r un,φ = rf, φ N
7 SPECTRAL COLLOCATION FOR A SINGULAR DE 339 or, taking integration by parts, du N u,rdφ + n N (3.11) r,φ = rf, φ N. Multiplying the continuous equation (.1) by rφ and integrating from r = to 1, we obtain du,rdφ + n u (3.1) r,φ = rf, φ. Then, subtracting (3.11) from (3.1) gives the error equation d(u u N ),r dφ u u + n N (3.13),φ = rf, φ rf, φ N, r which can be written in a simpler form as (3.14) where (3.15) a r,n (u u N,φ)=F r (φ) φ P N, du a r,n (u, v) =,rdv + n u, v, r v En = a r,n (v, v) 1/, F r (φ) = rf, φ rf, φ N. We first consider the case n. Recall that we have u N () = u() = (cf. (.6)). Let V g = {v H 1 (I) :v()=,v(1) = g}, V N g = V g P N. Lemma 3.4. Let u and u N be the solutions of the problem (.1) (.) (n,u()=)and the approximation (.1) (.11) (u N ()=), respectively. We have (3.16) u u N En sup ϕ V N F r (ϕ) ϕ En Proof. Equation (3.14) can be rewritten as a r,n (v u N,φ)=F r (φ)+a r,n (v u, φ) Taking φ = v u N V N results in Since v u N E n + inf v u En. v Vg N φ V N and v V N g. = a r,n (v u N F r (ϕ),φ) sup φ En + a r,n (v u, φ). ϕ V ϕ N En a r,n (v u, φ) r 1/ (v u) r r 1/ φ r + n r 1/ (v u) r 1/ φ ( r 1/ (v u) r + n r 1/ (v u) ) 1/ ( r 1/ φ r + n r 1/ φ ) 1/ Ca r,n (v u, v u) 1/ a r,n (φ, φ) 1/ = C v u En v u N En,
8 34 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN we have ( v u N E n C sup ϕ V N Then the desired result follows from ) F r (ϕ) + v u En v u N En. ϕ En u u N En u v En + v u N En. We now use Lemma 3.4 to obtain the estimate of u u N En. For the first term on the right-hand side of (3.16), we have from the Cauchy Schwarz inequality and Lemmas that, for any φ P N, F r (φ) = rf, φ rf, φ N = rf, φ P N 1 (rf),φ + P N 1 (rf),φ I N (rf),φ N rf P N 1 (rf),φ + P N 1 (rf) I N (rf),φ N (3.17) rf P N 1 (rf) φ + C P N 1 (rf) I N (rf) φ N rf P N 1 (rf) φ + C ( rf P N 1 (rf) + rf I N (rf) ) φ = C ( rf P N 1 (rf) + rf I N (rf) ) φ. By Lemmas 3.1 and 3. (and taking v = rf H m, m := max {m 1, 1}), we obtain F r (φ) CN 1 m (r r ) m 1 d m (rf) m φ CN 1 m (r r ) m 1 d m (rf) (3.18) m φ E n. For the second term on the right-hand side of (3.16), taking v = I N u and using the definition of the energy norm leads to (3.19) v u En v u 1 + [r(1 r)] 1/ (v u) CN 1 m (r r ) m 1 d m u m. Substituting (3.18) and (3.19) into (3.16), we obtain the estimate ( u u N En CN 1 m (r r ) m 1 d m u (3.) m + (r r ) m 1 We now consider the case n =. Recall again that we have dun Define W g = {v H 1 (I) :v(1) = g}, W N g = W g P N. d m (rf) m ). () = du ()=. Lemma 3.5. Let u and u N be the solutions of the problem (.1) (.) (n =, du dun ()=)and the approximation (.1) (.11) ( ()=), respectively. We have u u N E sup ϕ W N F r (ϕ) ϕ E + inf v W N g v u E.
9 SPECTRAL COLLOCATION FOR A SINGULAR DE 341 Proof. Equation (3.14) can be written as a r, (v u N,φ)=F r (φ)+a r, (v u, φ) φ W N and v W N g. Taking φ = v u N W N,wehave v u N E = a r, (v u N,φ) ( sup ϕ W N F r (ϕ) + v u E ϕ E ) v u N E. Then the conclusion follows. From the Hardy-type inequality v dv r = v E v W, Lemma 3.5 leads to the same result as in (3.). Hence, we have proved the following theorem. Theorem 3.1. For any integer m 1, the Legene-collocation approximation u N defined by the scheme (.1) (.1) for the problem (.1) and (.) satisfies ( u u N En CN 1 m (r r ) m 1 d m u m + (r r ) m 1 d m ) (rf) (3.1) m, where u is the exact solution of (.1) and (.) and m = max{m 1, 1}. This theorem shows that the Legene collocation approximation is convergent and the error decays faster than algebraically, provided that the right-hand-side term f and the solution u are infinitely differentiable. As shown in [1], the correct regularity requirement for u and f should be considered in a weighted Sobolev space { Hr s s = v dlv } (3.) r l < l= for some integer s. It is not difficult to see that the terms in the bracket on the right-hand side of (3.1) are bounded for u Hr m and f Hr m with m. In this sense, the result of Theorem 3.1 is optimal Convergence analysis of the Chebyshev method. We now consider the convergence of the Chebyshev method (.1) (.1). Let φ be the same as in (3.5). As for (3.1), we have r d u N (3.3) dun + n r un,φ = rf, φ,n or, taking integration by parts, du N u,rd(φ) + n N (3.4) r,φ = rf, φ,n. The error equation reads as d(u u N ),r d(φ) u u + n N,φ r (3.5) = rf, φ rf, φ,n.
10 34 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN We also write it in a simpler form (3.6) where (3.7) a r,n, (u u N,φ)=F r, (φ) φ P N, a r,n, (u, v) =b r, (u, v)+n u, v r du b r, (u, v) =,rd(v), F r, (φ) = rf, φ rf, φ,n., It is known that the nonsymmetric bilinear form b r, (, ), without the factor r, is coercive (see [3, 7, 8, 17, 15]). On the other hand, in the current situation, b r, (v, v) can become negative for some polynomials subject to the boundary conditions v( 1) = v(1) = or dv ( 1) = v(1) =. We have the following Gårding-type inequality. Lemma 3.6. For all u, v H, 1 we have (3.8) (3.9) 1 4 dv r + 3 v 8 1 v 1 r 8 b r, (v, v) dv r r b r, (u, v) 3 du r dv r. Proof. For notational simplicity, define ( ) dv (3.3) I 1 (v) = r. We have from integrating by parts dv b r, (v, v) =I 1 (v) v 1 r (r r r ) = I 1 (v) 1 8 v r(1 r +4r ) 5 (3.31) = I 1 (v) 3 v r Thus b r, (v, v) I 1 (v). On the other hand, v r(1 r) 5., (3.3) which gives (3.33) = I 1 (v)+ ( ) dv + vr r = I 1 (v) 1 b r, (v, v) 1 4 I 1(v)+ 3 8 v r v r 5, d(v ) r 3 v r (1 r) v r.
11 SPECTRAL COLLOCATION FOR A SINGULAR DE 343 The result (3.8) follows. To prove (3.9), we estimate b r, (u, v) by (3.34) b r, (u, v) du dv 1 r + du v 1 r (r r r ) [I 1 (u)] 1/ [I 1 (v)] 1/ + I (u, v), where (3.35) I (u, v) = du v 1 r (r r r ) [I 1(u)] 1/ [I 3 (v)] 1/ with (3.36) I 3 (v) = (1 r) v 4(r r ) r = 1 4 v r(1 r) 5. On the other hand, from integrating by parts, (3.37) I (v, v) = 1 8 v (r + r(1 r) ) 5 1 I 3(v). Thus we get from (3.37) and (3.35) (3.38) I 3 (v) I (v, v) [I 1 (v)] 1/ [I 3 (v)] 1/, which gives I 3 (v) 4I 1 (v) and (3.39) b r, (u, v) [I 1 (u)] 1/ ([I 1 (v)] 1/ +[I 3 (v)] 1/ ) 3[I 1 (u)] 1/ [I 1 (v)] 1/. We first consider the case n. Let V g and V N g ( dv v En, = r ) + n v 1/ r. be the same as before and Lemma 3.7. Let u and u N be the solutions of the problem (.1) (.) (n,u()=)and the Chebyshev approximation (u N ()=), respectively. We have (3.4) u u N En, C sup ϕ V N F r, (ϕ) ϕ En, + C inf v u En,. v Vg N Proof. Equation (3.6) can be rewritten as a r,n, (v u N,φ)=F r, (φ)+a r,n, (v u, φ) φ V N and v V N g. Taking φ = v u N V N and using the inequality (3.8) of Lemma 3.6 yields v u N E n, Ca r,n, (v u N F r, (ϕ),φ) C sup φ En, + Ca r,n, (v u, φ). ϕ V ϕ N En,
12 344 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN From (3.9) of Lemma 3.6 we have a r,n, (v u, φ) C d(v u) r dφ r + n v u φ r r ( ) d(v u) C r + n v u 1/ ( d(φ) r r = C v u En, v u N En,, and therefore ( v u N E n, C sup w V N + n φ r ) F r, (w) + v u En, v u N En,. w En, Then, the desired result follows from the triangular inequality: u u N En, u v En, + v u N En,. ) 1/ We now use Lemma 3.7 to obtain the estimate of u u N En,. For the first term on the right-hand side of (3.4), we have from the Cauchy Schwarz inequality and Lemmas that, for any φ P N, F r, (φ) = rf, φ rf, φ,n = rf, φ P N 1 (rf),φ + P N 1 (rf),φ I N (rf),φ,n = C ( rf P N 1 (rf) + rf I N ) (rf) φ. Then, by Lemmas 3.1 and 3., (3.41) F r, (φ) CN 1 m (r r ) m 1 CN 1 m (r r ) m 1 d m (rf) m d m (rf) m φ φ En,. For the second term on the right-hand side of (3.16), taking v = I N u leads to v u En, d(v u) r v u + r(1 r) CN 1 m (r r ) m 1 d m u (3.4) m. Substituting (3.41) and (3.4) into (3.16), we obtain the estimate (3.43) u u N En, ( CN 1 m (r r ) m 1 d m u m + (r r ) m 1 d m (rf) m ). We now consider the case n =. In this case, a r,, (, ) =b r, (, ) is not coercive. Thus it does not seem likely to us that an error bound can be obtained in the weighted
13 SPECTRAL COLLOCATION FOR A SINGULAR DE 345 energy norm E,. For this reason, we conduct the estimation in the energy norm E without the Chebyshev weight. We first note that the polar condition () = allows us to extend the collocation equations du N d u N (3.44) r j, (r j,) dun (r j,) =r j, f(r j, ), 1 j N 1, to the point r = r j, =. Since the left-hand side of (3.44) is a polynomial of degree N 1, (3.44) holds for all r [, 1], provided that f P N. We introduce an auxiliary interpolation operator ĨN : C([, 1]) P N defined by (3.45) Ĩ N v(r j )=v(r j ), 1 j N 1. Thus we are able to rewrite (3.44) as (3.46) r d u N dun (r) (r) =rĩn f(r), r 1. Let W g and Wg N be the same as before. Lemma 3.8. Let u and u N be the solutions of the problem (.1) (.) (n =, du dun ()=)and the Chebyshev collocation approximation ( ()=), respectively. We have u u N E sup ϕ W N where F r (ϕ) = rf rĩn f,ϕ. Proof. We have from (.1) and (3.46) F r (ϕ) ϕ E + inf v W N g v u E, a r, (v u N,φ)= F r (φ)+a r, (v u, φ) φ W N and v W N g. Taking φ = v u N W N,wehave v u N E = a r, (v u N,φ) ( ) F r (ϕ) + v u E v u N E, ϕ E sup ϕ W N which gives the desired result. We need further to estimate the term F r (ϕ). According to the definition, it is easy to see that (r r )ĨN f = I N ((r r )f). Therefore, we have from the Cauchy Schwarz inequality and Lemma 3. that, for any ϕ W N, (3.47) F r (ϕ) = (r r )f I N ((r r )f), (1 r) 1 ϕ (r r )f I N ((r r )f) (1 r) 1 ϕ CN 1 m (r r ) m 1 d m ((r r )f) m ϕ E, where we have used (3.48) (1 r) 1 ϕ dϕ r ϕ W N,
14 346 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN which can be derived from ( r dϕ 1 ) (1 r) 1 ϕ = dϕ r (1 r) 1 ϕ 1 d(ϕ) r1/ (1 r) 1 dϕ (3.49) r 1 4 (1 r) 1 ϕ. Hence, we have proved the following theorem. Theorem 3.. For any integer m 1, the Chebyshev-collocation approximation u N defined by the scheme (.1) (.1) for the problem (.1) and (.) satisfies ( u u N En CN 1 m (r r ) m 3 d m u 4 m + (r r ) m 3 d m (rf) 4 m + (r r ) m 3 d m ((r r ) )f) (3.5) 4 m, where u is the exact solution of (.1) and (.) and m = max{m 1, 1}. For the case n, (3.5) also holds in the stronger norm En,. Thus, we obtain a convergence result similar to that of the Legene collocation method. As in Theorem 3.1, when u Hr m and f Hr m with m 5/, the terms in the bracket on the right-hand side of (3.5) are bounded. 4. Numericalexperiments. In this section we present some numerical results to demonstrate the accuracy of the Legene and Chebyshev collocation methods (.1) (.1) for the model problem (.1) (.). Example 1. The function f(r) and the Dirichlet boundary condition at r = 1 are chosen such that the exact solution of (.1) and (.) is (4.1) u(r) =r cos(1πr), <r<1. We note that the energy norm En is stronger than the L norm for n>. (This is not true for n =.) For this reason, we use the maximum norm to measure the error for the case n>, but use the energy norm for the case n =. These norms are numerically approximated in the computations, viz., 1/ N du E,N = r j (r j) dun (r (4.) j) w j for n =, (4.3) j= E 1,N = max j N u(r j,) u N (r j, ) for n =1, where the Legene points {r j } are used in (4.) and both the Legene and Chebyshev points {r j, } are used in (4.3). The numerical results for the cases n = and n = 1 are listed in Table 1. The spectral convergence of the methods is clearly shown in the table. One may also notice that the Legene and Chebyshev collocation methods produce very comparable results. Example. The function f(r) and the boundary condition at r = 1 are chosen such that the problem has a less regular solution (4.4) u(r) =r 5/, <r<1.
15 SPECTRAL COLLOCATION FOR A SINGULAR DE 347 Table 1 Numerical results obtained with the Legene (LC) and Chebyshev (CC) collocation methods for Example 1. The convergence order is N Order. n = n =1 LC-method CC-method LC-method CC-method N E,N Order E,N Order E 1,N Order E 1,N Order 1 6.9e+1 1.5e+ 1.7e+1 4.1e+1.1e e e+. 5.e e e e e e e e e e e e e e e e e e 13 4.e 1 1.5e 14 1.e 13 Table Numerical results obtained with the Legene collocation method for Example. The convergence order is N Order. n = n =1 N E,N Order E 1,N Order 4 3.e 7 1.7e e e e e e e e e e e The computation is done with the Legene collocation method for n = and n =1. The solution error and the convergence order are listed in Table. One can easily see that E,N O ( N 5) and E 1,N O ( N 5). That is, the rate of convergence is nearly twice the exponent of r in (4.4), 5/. On the other hand, it is not difficult to show that the first term on the right-hand side of (3.1) is bounded for m<5, while the second term is bounded for m <3orm<4for the current example. Thus, the right-hand-side terms are bounded for m<4. From Theorem 3.1, we have u u N En O(N 3 ). This indicates that the convergence rate predicted by (3.1) is not sharp, although the estimate is optimal according to the regularity requirement ([1]; also cf. (3.)). Such an order loss seems typical in the convergence analysis of collocation schemes, especially for problems involving force terms; e.g., see [4] for comparison of typical estimates for the Legene Galerkin method ((8.7) on p. 74) and the Legene collocation method ((15.15) on p. 31). It is interesting to note that sharp estimates have been obtained for the p-version finite element method (which is of Galerkin type); e.g., see Babuska and Suri [5]. Finally, we mention that the Chebyshev collocation method leads to very comparable results. 5. Conclusions and comments. In the previous sections we have proved that the Legene and Chebyshev collocation approximations presented in [1] are convergent and that the error decays faster than algebraically when f and u are infinitely differentiable for the singular problem (.1) and (.). Our main results are given in Theorems 3.1 and 3.. The key feature of the spectral collocation approximation is that it uses a pole condition as the boundary condition at the singularity and employs a commonly used collocation scheme. Thus, the convergence result provides a theoretical justification for the well-known fact that it is unnecessary to impose all the pole conditions in order for numerical schemes to obtain high accuracy. Because most of the existing spectral
16 348 WEIZHANG HUANG, HEPING MA, AND WEIWEI SUN approaches for singular problems use more or less the pole conditions, we expect that our result can also be regarded as a theoretical justification for these methods. Finally we make a few comments on the method we analyzed. The method has been successfully applied to solving steady-state Navier Stokes equations in [13, 14]. However, since it uses the Chebyshev or Legene type of collocation methods in the r interval (,1), one may suspect that the clustering of grid points near r = leads to a very severe restriction on time steps for time dependent problems. To see this, let us consider the time dependent problem on the unit disk (5.1) u t = u + au x + bu y. In polar coordinates the equation becomes u t = u r + 1 u r r + 1 u r θ (5.) +(a cos θ + b sin θ) u r +( asin θ + b cos θ)1 u r θ. Assume that (5.) is approximated in r using the Legene or Chebyshev collocation, and in θ using Fourier collocation. Then it is not difficult to see that for the diffusion term the time step restriction for an explicit integration scheme is ( ) 1 (5.3) t max min{( r ),r 1 r } (r 1 ) = O N 4 at r and ( ) 1 (5.4) t max min{( r N 1 ),r N r N 1 } N 4 at r 1. Obviously these two time scales are the same. For the convection term we have at r ( ) 1 (5.5) t max min{ r,r 1 } O N and at r 1 ( ) 1 (5.6) t max min{ r N 1,r N } O N. Once again they are the same. Thus, the above simple analysis tells us that the clustering of grid points near the singularity does not result in a time restriction worse than that near the outer boundary. Of course, just like spectral methods applied to nonsingular problems, a restriction O(1/N 4 ) on time steps is too severe. Implicit or semi-implicit time integrators should be used. The resultant algebraic systems can be solved using either iterative methods with effective preconditioners [7, 1] or fast direct solvers [4]. For problems in spheric geometries the method can be applied straightforwardly. However, the severe restriction on time steps at the north and south poles could be a potential problem for the method (see discussion in [6, pp ]). This issue deserves further investigation. Acknowledgment. The authors would like to thank the anonymous referees for their valuable comments.
17 SPECTRAL COLLOCATION FOR A SINGULAR DE 349 REFERENCES [1] C. Bernardi, M. Dauge, and Y. Maday, Spectral Methods for Axisymmetric Domains, Gauthier-Villars, Éditions Scientifiques et Médicales Elsevier, Paris, [] C. Bernardi and Y. Maday, Properties of some weighted Sobolev spaces and application to spectral approximations, SIAM J. Numer. Anal., 6 (1989), pp [3] C. Bernardi and Y. Maday, Polynomial interpolation results in Sobolev space, J. Comput. Appl. Math., 43 (199), pp [4] C. Bernardi and Y. Maday, Spectral Methods, Techniques of Scientific Computing (Part ), Handb. Numer. Anal. 5, P. Ciarlet and J. L. Lions, eds., Elsevier, Amsterdam, 1997, pp [5] I. Babuška and M. Suri, The p and h-p versions of the finite element method, Basic principles and properties, SIAM Rev., 36 (1994), pp [6] J. P. Boyd, Chebyshev and Fourier Spectral Methods, Springer-Verlag, Berlin, [7] C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer-Verlag, Berlin, [8] C. Canuto and A. Quarteroni, Variational methods in the theoretical analysis of spectral approximations, in Spectral Methods for Partial Differential Equations (NASA Langley Research Center, 198), R. G. Voigt, D. Gottlieb, and M. Y. Hussaini, eds., SIAM, Philadelphia, PA, 1984, pp [9] H. Eisen, W. Heinrichs, and K. Witsch, Spectral collocation methods and polar coordinate singularities, J. Comput. Phys., 96 (1991), pp [1] B. Fornberg, A pseudospectral approach for polar and spherical geometries, SIAM J. Sci. Comput., 16 (1995), pp [11] D. Gottlieb and S. A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications, CBMS-NSF Regional Conf. Ser. in Appl. Math., 6, SIAM, Philadelphia, [1] W. Huang and D. M. Sloan, Pole condition for singular problems: The pseudospectral approximation, J. Comput. Phys., 17 (1993), pp [13] W. Huang and T. Tang, Pseudospectral solutions for steady motion of a viscous fluid inside a circular boundary, Appl. Numer. Math., 33 (), pp [14] A. Karageorghis and T. Tang, A spectral domain decomposition approach for steady Navier Stokes problems in circular geometries, Comput. Fluids, 5 (1996), pp [15] S. D. Kim and S. V. Parter, Preconditioning Chebyshev spectral collocation by finite-difference operators, SIAM J. Numer. Anal., 34 (1997), pp [16] J. Li, H.-P. Ma, and W.-W. Sun, Error analysis for solving the Korteweg de Vries equation by a Legene pseudo-spectral method, Numer. Methods Partial Differential Equations, 16 (), pp [17] H.-P. Ma and B.-Y. Guo, The Chebyshev spectral method for Burgers-like equations, J. Comput. Math., 6 (1988), pp [18] H.-P. Ma and W.-W. Sun, A Legene Petrov Galerkin and Chebyshev collocation method for third-order differential equations, SIAM J. Numer. Anal., 38 (), pp [19] H.-P. Ma and W.-W. Sun, Optimal error estimates of the Legene Petrov Galerkin method for the Korteweg de Vries equation, SIAM J. Numer. Anal., 39 (1), pp [] T. Matsushima and P. S. Marcus, A spectral method for polar coordinates, J. Comput. Phys., 1 (1995), pp [1] S. A. Orszag, Fourier series on spheres, Mon. Weather Rev., 1 (1974), pp [] V. G. Priymak, Pseudospectral algorithms for Navier Stokes simulation of turbulent flows in cylinical geometry with coordinate singularities, J. Comput. Phys., 118 (1995), pp [3] V. G. Priymak and T. Miyazaki, Accurate Navier Stokes investigation of transitional and turbulent flows in a circular pipe, J. Comput. Phys., 14 (1998), pp [4] J. Shen, Efficient spectral-galerkin methods III: Polar and cylinical geometries, SIAM J. Sci. Comput., 18 (1997), pp [5] J. Shen, Efficient spectral-galerkin methods IV: Spherical geometries, SIAM J. Sci. Comput., (1999), pp [6] D. J. Torres and E. A. Coutsias, Pseudospectral solution of the two-dimensional Navier Stokes equations in a disk, SIAM J. Sci. Comput., 1 (1999), pp [7] W. T. Verkley, A pseudo-spectral model for two-dimensional incompressible flow in a circular basin, I. Mathematical formulation, J. Comput. Phys., 136 (1997), pp [8] W. T. Verkley, A pseudo-spectral model for two-dimensional incompressible flow in a circular basin, II. Numerical examples, J. Comput. Phys., 136 (1997), pp
1462. 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 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 informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
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 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 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 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 informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
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 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 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 informationOptimal Spectral-Galerkin Methods Using Generalized Jacobi Polynomials
Journal of Scientific Computing ( 2006) DO: 10.1007/s10915-005-9055-7 Optimal Spectral-Galerkin Methods Using Generalized Jacobi Polynomials Ben-Yu Guo, 1 Jie Shen, 2 and Li-Lian Wang 2,3 Received October
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 informationCOMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE
COMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE BENGT FORNBERG and DAVID MERRILL Abstract. For modeling convective flows over a sphere or within a spherical
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationStability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation
J Sci Comput (01) 53:10 18 DOI 10.1007/s10915-01-961-8 Stability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation Sigal Gottlieb Cheng Wang
More informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationT u
WANG LI-LIAN Assistant Professor Division of Mathematical Sciences, SPMS Nanyang Technological University Singapore, 637616 T 65-6513-7669 u 65-6316-6984 k lilian@ntu.edu.sg http://www.ntu.edu.sg/home/lilian?
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
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 informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 informationCURRICULUM VITAE. Education: Teaching, Research, and Industrial Experience: List of Refereed Publications:
Last update: 12/12/2006 CURRICULUM VITAE Name : Weiming Cao Address : Department of Mathematics, The University of Texas at San Antonio San Antonio, TX 78249 Phone : (210)-458-7327 Fax: (210)-458-4439
More informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationInterpolation approximations based on Gauss Lobatto Legendre Birkhoff quadrature
Journal of Approximation Theory 161 009 14 173 www.elsevier.com/locate/jat Interpolation approximations based on Gauss Lobatto Legendre Birkhoff quadrature Li-Lian Wang a,, Ben-yu Guo b,c a Division of
More informationFast Direct Solver for Poisson Equation in a 2D Elliptical Domain
Fast Direct Solver for Poisson Equation in a 2D Elliptical Domain Ming-Chih Lai Department of Applied Mathematics National Chiao Tung University 1001, Ta Hsueh Road, Hsinchu 30050 Taiwan Received 14 October
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 informationDepartment of Mathematics
INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY Department of Mathematics MA 04 - Complex Analysis & PDE s Solutions to Tutorial No.13 Q. 1 (T) Assuming that term-wise differentiation is permissible, show that
More informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
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 informationPost-processing of solutions of incompressible Navier Stokes equations on rotating spheres
ANZIAM J. 50 (CTAC2008) pp.c90 C106, 2008 C90 Post-processing of solutions of incompressible Navier Stokes equations on rotating spheres M. Ganesh 1 Q. T. Le Gia 2 (Received 14 August 2008; revised 03
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
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 informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationAbsorbing boundary conditions for simulation of gravitational waves with spectral methods in spherical coordinates
Absorbing boundary conditions for simulation of gravitational waves with spectral methods in spherical coordinates Jérôme Novak, Silvano Bonazzola arxiv:gr-qc/0203102v2 29 Apr 2004 Laboratoire de l Univers
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
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 informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationAN ALGORITHM FOR COMPUTING FUNDAMENTAL SOLUTIONS OF DIFFERENCE OPERATORS
AN ALGORITHM FOR COMPUTING FUNDAMENTAL SOLUTIONS OF DIFFERENCE OPERATORS HENRIK BRANDÉN, AND PER SUNDQVIST Abstract We propose an FFT-based algorithm for computing fundamental solutions of difference operators
More informationRelationships between upper exhausters and the basic subdifferential in variational analysis
J. Math. Anal. Appl. 334 (2007) 261 272 www.elsevier.com/locate/jmaa Relationships between upper exhausters and the basic subdifferential in variational analysis Vera Roshchina City University of Hong
More informationOSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION
Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY
More informationA HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD
SIAM J. SCI. COMPUT. Vol. 26, No. 5, pp. 1504 1524 c 2005 Society for Industrial and Applied Mathematics A HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD E.
More information, /'/ NI'9 ; / Proceedings of the Third International. High Order Methods. Conference on Spectral and. Houston, Texas
; /, /'/ NI'9 - Edited by Andrew V. Ilin and L. Ridgway Scott Proceedings of the Third International Conference on Spectral and High Order Methods Houston, Texas 5-9 June 1995 Houston Journal of Mathematics
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationPSEUDOSPECTRAL SOLUTION OF THE TWO-DIMENSIONAL NAVIER STOKES EQUATIONS IN A DISK
SIAM J. SCI. COMPUT. Vol. 1, No. 1, pp. 378 43 c 1999 Society for Industrial and Applied Mathematics PSEUDOSPECTRAL SOLUTION OF THE TWO-DIMENSIONAL NAVIER STOKES EQUATIONS IN A DISK D. J. TORRES AND E.
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationThe need for de-aliasing in a Chebyshev pseudo-spectral method
The need for de-aliasing in a Chebyshev pseudo-spectral method Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (March 2) Abstract In the present report, we
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationON THE CONVERGENCE OF THE FOURIER APPROXIMATION FOR EIGENVALUES AND EIGENFUNCTIONS OF DISCONTINUOUS PROBLEMS
SIAM J. UMER. AAL. Vol. 40, o. 6, pp. 54 69 c 003 Society for Industrial and Applied Mathematics O THE COVERGECE OF THE FOURIER APPROXIMATIO FOR EIGEVALUES AD EIGEFUCTIOS OF DISCOTIUOUS PROBLEMS M. S.
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationComputing the effective diffusivity using a spectral method
Materials Science and Engineering A311 (2001) 135 141 www.elsevier.com/locate/msea Computing the effective diffusivity using a spectral method Jingzhi Zhu a, *, Long-Qing Chen a, Jie Shen b, Veena Tikare
More informationJACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL
Bull. Korean Math. Soc. 53 (016), No. 1, pp. 47 6 http://dx.doi.org/10.4134/bkms.016.53.1.47 JACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL Yin Yang Abstract.
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More informationSpherical Harmonics and Related Topics. David Randall 2 S = 0, r 2 r r S 2. S = r n Y n
! Revised April 10, 2017 1:32 PM! 1 Spherical Harmonics and Related Topics David Randall The spherical surface harmonics are convenient functions for representing the distribution of geophysical quantities
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
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 informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationLARGE TIME BEHAVIOR OF SOLUTIONS TO THE GENERALIZED BURGERS EQUATIONS
Kato, M. Osaka J. Math. 44 (27), 923 943 LAGE TIME BEHAVIO OF SOLUTIONS TO THE GENEALIZED BUGES EQUATIONS MASAKAZU KATO (eceived June 6, 26, revised December 1, 26) Abstract We study large time behavior
More informationA new class of highly accurate differentiation schemes based on the prolate spheroidal wave functions
We introduce a new class of numerical differentiation schemes constructed for the efficient solution of time-dependent PDEs that arise in wave phenomena. The schemes are constructed via the prolate spheroidal
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More information31. Successive Subspace Correction method for Singular System of Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 2003 DDM.org 31. Successive Subspace Correction method for Singular
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
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 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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationBernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of
More information