Studia Geophysica et Geodaetica. Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids
|
|
- Alice Aubrey Marsh
- 5 years ago
- Views:
Transcription
1 Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids Journal: Studia Geophysica et Geodaetica Manuscript ID: SGEG-0-00 Manuscript Type: Original Article Date Submitted by the Author: -Jun-0 Complete List of Authors: Costea, Adrian; Leibniz University Hanover, Institute for Applied Mathematics Stephan, Ernst; Leibniz University Hanover, Institute for Applied Mathematics LeGia, Quoc; University of New South Wales, School of Mathematics and Statistics Tran, Thanh; University of New South Wales, School of Mathematics and Statistics Keywords: Boundary integral equations, meshless methods, spheroids, overlapping Schwarz, domain decomposition
2 Page of Studia Geophysica et Geodaetica Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids A. Costea Q.T. Le Gia E.P. Stephan T. Tran Abstract We consider the exterior Neumann problem of the Laplacian with boundary condition on a prolate spheroid. We propose to use spherical radial basis functions in the solution of the boundary integral equation arising from the Dirichlet-to-Neumann map. Our meshless approach with radial basis functions is particularly suitable for handling scattered satellite data. We also propose a preconditioning technique based on an overlapping domain decomposition method to deal with ill-conditioned matrices arising from the approximation problem. Here we report from Le Gia et al. (0), Costea et al. (Solution to the Neumann problem exterior to an oblate spheroid by radial basis functions, in preparation, 0) on a meshless method with radial basis functions for the Neumann problem for the Laplacian exterior to an oblate or a prolate spheroid. In geophysical applications, see Freeden et al. () and Krumm et al. (00), one is interested in such exterior Neumann problems where the orbits of satellites are located on spheroids. The satellites create data which amount to boundary conditions given in scattered points. A key tool of our approach is the use of the Dirichlet-to-Neumann map which directly converts the bounday value problem into a pseudodifferential equation on the spheroid. This integral equation is then handled with Fourier techniques by expansion into appropriate spherical harmonics. In Huang and Yu (00) this pseudodifferential equation was solved numerically with Institut für Angewandte Mathematik and QUEST (Centre for Quantum Engineering and Space-Time Research), Leibniz Universität Hannover, Welfengarten, 0 Hannover, Germany. costea@ifam.uni-hannover.de, stephan@ifam.uni-hannover.de School of Mathematics and Statistics, University of New South Wales, Sydney, NSW 0, Australia. qlegia@unsw.edu.au, Thanh.Tran@unsw.edu.au,
3 Page of standard boundary elements on a regular grid on the angular domain of the spherical coordinates. Our approach uses spherical radial basis functions (SRBF s) instead, allowing for better handling of scattered data. The error analysis for the meshless method with radial basis functions on the sphere as done in Tran et al. (00) has been extended to spheroids in Le Gia et al. (0) and Costea et al. (Solution to the Neumann problem exterior to an oblate spheroid by radial basis functions, in preparation, 0). Again the smooth solution of the pseudodifferential equation can be approximated with high convergence rates by the Galerkin solution consisting of radial basis functions. The numerical solution is obtained by an appropriate implementation of the prolate and oblate spheroidal harmonics and by truncating the resulting infinite dimensional discrete Galerkin system. We present numerical results of our meshless boundary element method for points from an equal area partitioning algorithm, so-called Saff points (Fig.), and for scattered data points from satellite observations (Fig.). Furthermore we present iteration numbers for the conjugate gradient method applied to solve the discrete systems in case of scattered data. We list the iteration numbers for the preconditioned conjugate gradient method when an overlapping additive Schwarz method is used as preconditioner. For both oblate and prolate spheroids we obtain only mildly growing iteration numbers for the overlapping additive Schwarz method. This technique was firstly analysed for pseudodifferential equations on the sphere in Tran et al. (0). Let Γ 0 = {(x,x,x ) R : x +x + x b =, a > b > 0} be a prolate a spheroid and Ω c be the unbounded domain outside the boundary Γ 0. We consider the exterior Neumann problem: given g defined on Γ 0, find U defined on Ω c satisfying U = 0 in Ω c ν U = g on Γ 0 (0.) U(x) = O( x ) as x where x denotes the Euclidean norm of x and ν denotes the unit outward normal vector on Γ 0. The solution is given by the series U(µ,θ,ϕ) = n=0 m= n Q m n (cosh µ) Q m n (cosh µ 0 )ûnmy nm (θ,ϕ), µ µ 0 > 0. (0.) where (µ,θ,ϕ) denote the prolate spheroidal coordinates of x Γ 0 (Huang and Yu (00)). Here Q m n are the associated Legendre functions of second kind,
4 Page of Studia Geophysica et Geodaetica Y nm the spherical harmonics of degree n, and û nm the expansion coefficients. In Huang and Yu (00) it is shown that (0.) is equivalent to Ku = g on Γ 0 (0.) with the Dirichlet-to-Neumann map (Steklov-Poincaré operator) K. Its weak formulation is: Find u H / (Γ 0 ) satisfying D(u,v) := (Ku)v ds = gv ds v H / (Γ 0 ) (0.) Γ 0 Γ 0 with D(u,v) = f 0 n=0 m= n H m n (cosh µ 0 )û nm v nm, (0.) where H m n (x) := (x )dq m n (x)/dx Q m n (x). (0.) Here we recall that the Sobolev space H s (Γ 0 ) is defined by u H s (Γ 0 ),s R u H s (Γ 0 ) := n=0 m= n ( + n ) s û nm < (0.) Proposition: There exists a unique solution for the variational problem (0.). This is due to the Lax-Milgram theorem since D(, ) is continuous and coercive on the Sobolev space H / (Γ 0 ) and the right hand side in(0.) defines a continuous linear functional on H / (Γ 0 ). The approximate solution to (0.) is sought in a finite dimensional subspace of H / (Γ 0 ). In order to use SRBFs we take the bijection ω : Γ 0 S (where S is the unit sphere in R ) ω(x) = (sin θ cos ϕ,sin θ sin ϕ,cos θ), (0.) and we define a reproducing kernel on Γ 0 as Ψ(x,x ) = Φ(ω(x),ω(x )), x,x Γ 0, (0.) where Φ is defined via a univariate function φ : [,] R by Φ(y,z) = φ(y z) y,z S (see Schoenberg () and Xu and Cheney ()). Here φ has a series expansion in terms of Legendre polynomials P n of degree n, as φ(t) = π + (n + )ˆφ(n)P n (t) with ˆφ(n) = π φ(t)p n (t)dt. (0.) n=0
5 Page of The kernel Ψ can be expanded into a series of spherical harmonics as Ψ(x,x ) = n=0 m= n where we choose φ such that ˆφ(n)Y nm (ω(x))y nm(ω(x )) (0.) ˆφ(n) ( + n ) τ for some τ > 0. (0.) Given a set of scattered data points X = {x,...,x M } Γ 0, we take V τ := span{ψ j := Ψ(x j, ) : x j X}. Now, the solution of (0.) is approximated by u X V τ satisfying the Galerkin equation D(u X,v) = gv ds v V τ. (0.) Γ 0 To this end we have to solve a linear system Ac = g (0.) where A is a matrix with entries A i,j = D(Ψ i,ψ j ), i,j =,...,M, and g is a vector with entries g j = Γ 0 gψ j ds, j =,..,M. Using (0.) and (0.) we can write A i,j = f 0 n=0 m= n [ˆφ(n)] H m n (cosh µ 0 )Y nm (ω(x i ))Y nm(ω(x j )). (0.) Further let N denote the number of the series terms of the truncated matrix element A N i,j. We compute the actual Galerkin approximation u X by solving (0.) with the Galerkin matrix A N obtained via the truncated entries A N i,j. Let Y = {y,...,y M } be the image of X under the map ω, i.e. y j = ω(x j ) for j =,...,M. As Y is a set of scattered points on S, we define the mesh norm h Y of Y as usual as h Y = sup min y S y j Y cos (y y j ). Proposition (Le Gia et al. (0)): For a prolate spheroid and truncation number N chosen sufficiently large there holds the quasioptimal error estimate for the difference between the exact solution u H s (Γ 0 ) of (0.) and the Galerkin solution u X V τ, τ > s + of (0.) u u X H / (Γ 0 ) Chs / Y, (0.)
6 Page of Studia Geophysica et Geodaetica where the constant C is independent of the set X used to define Ψ j, but may depend on N. (For the corresponding result for the oblate spheroid see Costea et al. (Solution to the Neumann problem exterior to an oblate spheroid by radial basis functions, in preparation, 0).) Our numerical experiments plotted in (Fig. ) for the mesh of Saff points show these predicted convergence rates α, namely α =. for P m = 0(τ =.),α =. for Pm = (τ =.),α =. for Pm = (τ =.), cf. Table. To compute the entries A ij of the stiffness matrix given in (0.), we need to compute the spherical harmonics Y nm and the functions Hn m. The term Hn m(cosh µ 0) in the entry A i,j of the stiffness matrix (see (0.) and (0.)) can be computed by using the relation H m n (x)= (n m+)q m n+(x)/q m n(x)+(n+)x,m=0,...,n; n=0,,... (0.) For negative values of m we use the relation Hn m (x) = H m n (x). The right hand side terms g j are computed by using the Fourier coefficients of g and Ψ j and Parseval s identity. We use different kernels Ψ(x,x ) = Φ(ω(x),ω(x )) in our numerical experiments, where the functions Φ are restrictions to the sphere of two different classes of positive definite RBFs defined by Matérn and Wendland functions. The Matérn functions (or Sobolev splines) were introduced for statistical applications in Matérn (). They are defined by φ ν (r) = ν Γ(ν) rν / K ν / (r), ν > /, (0.) where K ν is the K-Bessel function of order ν. In R, the Fourier transform of φ decays like φ(ξ) ( + ξ ) ν. The Matérn kernels used in our experiments are listed in Table. When restricting to the sphere, the native space associated with the kernel Φ(y,z) := φ ν ( y z) is the Sobolev space H ν / (S ), see Narcowich et al. (00). The Wendland functions discussed in Wendland (00) are positive definite functions with compact support. For arbitrary y,z S, Φ(y,z) = ρ m ( y z) with ρ m being locally supported radial basis functions defined by Wendland; see Table. In this case (0.) holds for τ = m + / see Narcowich and Ward (00). Example : Problem (0.) with prolate spheroid Γ 0 (f 0 =, µ 0 = ), Neumann condition and exact solution sin θ cos ϕ( cosh µ + cosh µ cos θ) g(µ,θ,ϕ) = (0.) cosh µ cos θ(cosh µ + cos θ) / f 0
7 Page of sinhµ sin θ cos ϕ U(µ,θ,ϕ) = f0 (0.0) (cosh µ + cos θ)/. Let e := u u X where u(θ,ϕ) = U(µ 0,θ,ϕ), solves (0.) and u X solves (0.) with truncation number N trunc = 0. We compute e L (Γ 0 ) and e H / (Γ 0 ) approximately by N max = and e L (Γ 0 ) ( Nmax n=0 m= n ( Nmax e H / (Γ 0 ) ( + n ) / n=0 û Xnm û nm ) / (0.) m= n û Xnm û nm ) / (0.) with û Xnm = M ˆφ(n)c i= i Y nm (ω(x i )) and û nm is computed by an appropriate quadrature formula as shown in Driscoll and Healy (). Table gives the errors in the L (Γ 0 ) and H / (Γ 0 ) norms for scattered points. The matrix is ill conditioned and a preconditioner is required. The symbol of the operator K behaves like (n + ) /, see Le Gia et al. (0). Hence K is a pseudodifferential operator of order. This allows us to extend our analysis discussed in Tran et al. (0) for the overlapping Schwarz preconditioner on the sphere to the prolate spheroid. The preconditioner is defined by the additive Schwarz operator, using a subspace decomposition of V τ as V τ = V V J. These subspaces V j, j = 0,...,J, are defined from a decomposition of the set Y into overlapping subsets Y j, j = 0,...,J. These subsets are generated by the following algorithm: Select α (0,π/), β (0,π]; p = y Y ; Y 0 := {p }; Y := {y Y : cos (y p ) α}; J =, k = ; while Y Y k Y do k = k + ; p k is chosen from Y \ Y 0 such that cos (p k p k ) β; Y 0 := Y 0 {p k }; Y k := {y Y : cos (y p k ) α} end J = k. Each set Y k is a collection of points inside a spherical cap of radius α centered at p k. These sets define subdomains on the prolate spheroid.
8 Page of Studia Geophysica et Geodaetica The centers p k are chosen so that the geodesic distance between two successive centers is not less than β, see Le Gia et al. (0) This is to ensure that finally all points are considered, and thus the algorithm terminates. The set Y 0 serves as the coarse mesh in domain decomposition methods for finite element approximations. This set defines the space V 0 which provides global communication of data. Table shows the corresponding numbers of iteration of the preconditioned conjugate gradient method (PCG) with an overlapping additive Schwarz preconditioner; stopping criteria in both cases is relative tolerance. Errors of the same order as in the non-preconditioned case are obtained but PCG needs much smaller iteration numbers.(cpu: computational times in seconds, iter: numbers of iterations) In Costea et al. (Solution to the Neumann problem exterior to an oblate spheroid by radial basis functions, in preparation, 0) we have analysed the foregoing meshless boundary element method for an oblate spheroid. We obtain similar results both with respect to convergence of the meshless Galerkin solution as well as the behaviour of the iteration number of PCG. For brevity we present only some numerical results. Example : Problem (0.) with oblate spheroid Γ 0 (see also Costea et al. (Solution to the Neumann problem exterior to an oblate spheroid by radial basis functions, in preparation, 0)), with Neumann condition and exact solution g = sinhµsin θ cos ϕ(cosh µ+cos θ) f 0 (cosh µ cos θ) cosh µ sin θ cos ϕ, U =. (0.) f0 (cosh µ cos θ) Figure shows the experimental orders of convergence (EOC) for the errors in the H / (Γ 0 )-norm (energy norm) when Saff points and Wendland s kernels are used, for the prolate (Pm) and oblate (Om) spheroid. Figure shows the corresponding results when Matérn kernels are used for the prolate (Example ) and oblate (Example ) spheroid. Tables and give the errors in the L (Γ 0 ) and H / (Γ 0 ) norms and the iteration numbers for CG and PCG for scattered points for Example (stopping criteria in both cases is relative tolerance ). We observe for the oblate spheroid that the errors of the meshless boundary element method and the performance of the PCG with overlapping Schwarz preconditioner behave as in the case of the prolate spheroid.
9 Page of Figure : Image of Saff points on the prolate spheroid Figure : Image of satellite points on the prolate spheroid
10 Page of Studia Geophysica et Geodaetica ν φ ν (r) τ e r. e r ( + r). Table : Matérn s rbfs m ρ m (r) τ 0 ( r) +. ( r) +(r + ). ( r) +(r + r + ). Table : Wendland s RBFs m M q Y e L (Γ 0 ) e H / (Γ 0 ) cpu iter 0 π/0.0e-.e-. π/00.e-.0e-. π/0.0e-.e- π/00.0e-.e Table : Errors with scattered points from MAGSAT, using CG, ρ m (r), Ex. m M cos α cos β J cpu iter e L (Γ 0 ) E E E E E E E E E E E E E E E E- Table : Errors with scattered points from MAGSAT, for PCG, ρ m (r), Ex.
11 Page of H / -Error H / -Error H / -Error e-0 e-0 e-0 e-0 e-0 s=. α=. s=. α=.0 Pm=0 Pm= Pm= u-u X H / (Γ0 ) C h Y s-/ s=. α=. e Mesh norm h Y e-0 e-0 e-0 e-0 e-0 Om=0 Om= Om= Pm=0 Pm= Pm= e Mesh norm h Y e-0 e-0 e-0 e-0 Example φ Example φ Example φ Example φ e Mesh norm h Y Figure : Log-log plot for H / (Γ 0 ) errors using Wendland RBFs ρ m (r) for the prolate spheroid Figure : Log-log plot for H / (Γ 0 ) errors using Wendland RBFs ρ m (r) (Saff points) for the prolate (Pm) and oblate spheroid (Om) Figure : Log-log plot for H / (Γ 0 ) errors using Matérn RBFs φ ν (r) (Saff points) for the prolate (Ex.) and oblate spheroid (Ex.)
12 Page of Studia Geophysica et Geodaetica M q Y e L (Γ 0 ) e H / (Γ 0 ) ITER CPU 0 π/.0e-00.0e π/00.e-00.e π/0.e-00.0e Table : Errors with scattered points from MAGSAT, using CG, ρ 0 (r), Ex. M cos α cos β e L (Γ 0 ) e H / (Γ 0 ) ITER CPU E-00.0E E-00.0E E-00.0E E-00.0E E-00.E E-00.E E-00.E E-00.E E-00.0E E-00.0E E-00.0E E-00.0E-00. Table : Errors with scattered points from MAGSAT, for PCG, ρ 0 (r), Ex.
13 Page of References Driscoll, James R. and Healy, Jr., Dennis M.. Computing Fourier transforms and convolutions on the -sphere. Adv. in Appl. Math., ():0 0. Freeden, W., Gervens, T., and Schreiner, M.. Constructive approximation on the sphere. Numerical Mathematics and Scientific Computation. The Clarendon Press Oxford University Press, New York. With applications to geomathematics. Grafarend, E. W., Krumm, F. W., and Schwarze, V. 00. Geodesy: the Challenge of the rd Millennium. Springer, Berlin. Huang, H. and Yu, D. 00. The ellipsoid artificial boundary method for threedimensional unbounded domains. J. Comput. Math., (-):. Le Gia, Q. T., Tran, T., and Stephan, E. P. 0. Solution to the Neumann problem exterior to a prolate spheroid by radial basis functions. Adv. Comput. Math. DOI.0/s-0--. Matérn, B.. Spatial variation, volume of Lecture Notes in Statistics. Springer-Verlag, Berlin, second edition. Narcowich, F. J., Sun, X., and Ward, J. D. 00. Approximation power of RBFs and their associated SBFs: a connection. Adv. Comput. Math., ():. Narcowich, F. J. and Ward, J. D. 00. Scattered data interpolation on spheres: error estimates and locally supported basis functions. SIAM J. Math. Anal., (): 0 (electronic). Schoenberg, I. J.. Positive definite functions on spheres. Duke Math. J., :. Tran, T., Le Gia, Q. T., Sloan, I. H., and Stephan, E. P. 00. Boundary integral equations on the sphere with radial basis functions: error analysis. Appl. Numer. Math., ():. Tran, T., Le Gia, Q. T., Sloan, I. H., and Stephan, E. P. 0. Preconditioners for pseudodifferential equations on the sphere with radial basis functions. Numer. Mathematik, ():. Wendland, H. 00. Scattered data approximation, volume of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge. Xu, Y. and Cheney, E. W.. Strictly positive definite functions on spheres. Proc. Amer. Math. Soc., ():.
Solving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method
Solving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method A. Costea a, Q.T. Le Gia b,, D. Pham c, E. P. Stephan a a Institut für Angewandte Mathematik and QUEST (Centre for Quantum
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationOverlapping additive Schwarz preconditioners for interpolation on the unit sphere by spherical radial basis functions
Overlapping additive Schwarz preconditioners for interpolation on the unit sphere by spherical radial basis functions Q. T. Le Gia T. Tran I. H. Sloan October, 008 Abstract We present an overlapping domain
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationApplied and Computational Harmonic Analysis
Appl Comput Harmon Anal 32 (2012) 401 412 Contents lists available at ScienceDirect Applied and Computational Harmonic Analysis wwwelseviercom/locate/acha Multiscale approximation for functions in arbitrary
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 informationMultiscale RBF collocation for solving PDEs on spheres
Multiscale RBF collocation for solving PDEs on spheres Q. T. Le Gia I. H. Sloan H. Wendland Received: date / Accepted: date Abstract In this paper, we discuss multiscale radial basis function collocation
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationNEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE
NEEDLET APPROXIMATION FOR ISOTROPIC RANDOM FIELDS ON THE SPHERE Quoc T. Le Gia Ian H. Sloan Yu Guang Wang Robert S. Womersley School of Mathematics and Statistics University of New South Wales, Australia
More informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationINTEGRATION BY RBF OVER THE SPHERE
INTEGRATION BY RBF OER THE SPHERE ALISE SOMMARIA AND ROBERT S. WOMERSLEY Abstract. In this paper we consider numerical integration over the sphere by radial basis functions (RBF). After a brief introduction
More informationSOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS
SOLVING HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS IN SPHERICAL GEOMETRY WITH RADIAL BASIS FUNCTIONS NATASHA FLYER 1. Introduction Mathematical modeling of space and climate phenomena generally requires
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationLocal radial basis function approximation on the sphere
Local radial basis function approximation on the sphere Kerstin Hesse and Q. T. Le Gia School of Mathematics and Statistics, The University of New South Wales, Sydney NSW 05, Australia April 30, 007 Abstract
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
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 informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationLecture 2: Reconstruction and decomposition of vector fields on the sphere with applications
2013 Dolomites Research Week on Approximation : Reconstruction and decomposition of vector fields on the sphere with applications Grady B. Wright Boise State University What's the problem with vector fields
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 informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationAPPROXIMATION OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS ON SPHERES. A Dissertation QUOC THONG LE GIA
APPROIMATION OF LINEAR PARTIAL DIFFERENTIAL EQUATIONS ON SPHERES A Dissertation by QUOC THONG LE GIA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2 Part 3: Native Space for Positive Definite Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationResearch Institute of Geodesy, Topography and Cartography Zdiby, Prague-East, Czech Republic
Leibniz Society of Science at Berlin, Scientific Colloquium Geodesy - Mathematic - Physics - Geophysics in honour of Erik W. Grafarend on the occasion of his 75th birthday Berlin, Germany, 3 February 05
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
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 informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
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 informationMultiscale Modeling of Ocean Circulation
Multiscale Modeling of Ocean Circulation Willi Freeden (1), Dominik Michel (1 ), and Volker Michel (1) (1) Geomathematics Group, University of Kaiserslautern, P.O.Box 3049, 67653 Kaiserslautern, Germany
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 5: Completely Monotone and Multiply Monotone Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationConsistency Estimates for gfd Methods and Selection of Sets of Influence
Consistency Estimates for gfd Methods and Selection of Sets of Influence Oleg Davydov University of Giessen, Germany Localized Kernel-Based Meshless Methods for PDEs ICERM / Brown University 7 11 August
More informationThe Coupling Method with the Natural Boundary Reduction on an Ellipse for Exterior Anisotropic Problems 1
Copyright 211 Tech Science Press CMES, vol.72, no.2, pp.13-113, 211 The Coupling Method with the Natural Boundary Reduction on an Ellipse for Exterior Anisotropic Problems 1 Quan Zheng 2, Jing Wang 2 and
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 7: Conditionally Positive Definite Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationSolving approximate cloaking problems using finite element methods
Solving approximate cloaking problems using finite element methods Q. T. Le Gia 1 H. Gimperlein M. Maischak 3 E. P. Stephan 4 June, 017 Abstract Motivating from the approximate cloaking problem, we consider
More informationSpherical designs with close to the minimal number of points 1
Spherical designs with close to the minimal number of points Robert S. Womersley School of Mathematics and Statistics, University of New South Wales, Sydney, NSW, 05, Australia Abstract In this paper we
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationCubic spline Numerov type approach for solution of Helmholtz equation
Journal of Linear and Topological Algebra Vol. 03, No. 01, 2014, 47-54 Cubic spline Numerov type approach for solution of Helmholtz equation J. Rashidinia a, H. S. Shekarabi a and M. Aghamohamadi a a Department
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationCompression on the digital unit sphere
16th Conference on Applied Mathematics, Univ. of Central Oklahoma, Electronic Journal of Differential Equations, Conf. 07, 001, pp. 1 4. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationGeneralized Lax-Milgram theorem in Banach spaces and its application to the mathematical fluid mechanics.
Second Italian-Japanese Workshop GEOMETRIC PROPERTIES FOR PARABOLIC AND ELLIPTIC PDE s Cortona, Italy, June 2-24, 211 Generalized Lax-Milgram theorem in Banach spaces and its application to the mathematical
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationNumerical Minimization of Potential Energies on Specific Manifolds
Numerical Minimization of Potential Energies on Specific Manifolds SIAM Conference on Applied Linear Algebra 2012 22 June 2012, Valencia Manuel Gra f 1 1 Chemnitz University of Technology, Germany, supported
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 informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationOn a Generalization of the Busemann Petty Problem
Convex Geometric Analysis MSRI Publications Volume 34, 1998 On a Generalization of the Busemann Petty Problem JEAN BOURGAIN AND GAOYONG ZHANG Abstract. The generalized Busemann Petty problem asks: If K
More informationPositive Definite Functions on Spheres
Positive Definite Functions on Spheres Tilmann Gneiting Institute for Applied Mathematics Heidelberg University, Germany t.gneiting@uni-heidelberg.de www.math.uni-heidelberg.de/spatial/tilmann/ International
More informationOn the interplay between discretization and preconditioning of Krylov subspace methods
On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
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 informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes Fachrichtung 6.1 Mathematik Preprint Nr. 362 An iterative method for EIT involving only solutions of Poisson equations. I: Mesh-free forward solver Thorsten Hohage and Sergej
More informationQMC designs and covering of spheres by spherical caps
QMC designs and covering of spheres by spherical caps Ian H. Sloan University of New South Wales, Sydney, Australia Joint work with Johann Brauchart, Josef Dick, Ed Saff (Vanderbilt), Rob Womersley and
More informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationNumerical Integration in Meshfree Methods
Numerical Integration in Meshfree Methods Pravin Madhavan New College University of Oxford A thesis submitted for the degree of Master of Science in Mathematical Modelling and Scientific Computing Trinity
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationInvertibility of the biharmonic single layer potential operator
Invertibility of the biharmonic single layer potential operator Martin COSTABEL & Monique DAUGE Abstract. The 2 2 system of integral equations corresponding to the biharmonic single layer potential in
More informationHarmonic Polynomials and Dirichlet-Type Problems. 1. Derivatives of x 2 n
Harmonic Polynomials and Dirichlet-Type Problems Sheldon Axler and Wade Ramey 30 May 1995 Abstract. We take a new approach to harmonic polynomials via differentiation. Surprisingly powerful results about
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 information