DOCSLIB.ORG

Asymptotic Bounds to Outputs of the Exact Weak Solution of the Three-Dimensional Helmholtz Equation

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Asymptotic Bounds to Outputs of the Exact Weak Solution of the Three-Dimensional Helmholtz Equation Asymptotic Bounds to Outputs of the Exact Weak Solution of the Three-Dimensional Helmholtz Equation Shahin Ghomeshi A Thesis in The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy (Mechanical Engineering) at Concordia University Montreal, Quebec, Canada March 2010 ©Shahin Ghomeshi, 2010 Library and Archives Bibliotheque et 1*1 Canada Archives Canada Published Heritage Direction du Branch Patrimoine de I'edition 395 Wellington Street 395, rue Wellington Ottawa ON K1A 0N4 OttawaONK1A0N4 Canada Canada Your file Votre reference ISBN: 978-0-494-80163-5 Our file Notre reference ISBN: 978-0-494-80163-5 NOTICE: AVIS: The author has granted a non­ L'auteur a accorde une licence non exclusive exclusive license allowing Library and permettant a la Bibliotheque et Archives Archives Canada to reproduce, Canada de reproduire, publier, archiver, publish, archive, preserve, conserve, sauvegarder, conserver, transmettre au public communicate to the public by par telecommunication ou par I'lnternet, preter, telecommunication or on the Internet, distribuer et vendre des theses partout dans le loan, distribute and sell theses monde, a des fins commerciales ou autres, sur worldwide, for commercial or non­ support microforme, papier, electronique et/ou commercial purposes, in microform, autres formats. paper, electronic and/or any other formats. The author retains copyright L'auteur conserve la propriete du droit d'auteur ownership and moral rights in this et des droits moraux qui protege cette these. Ni thesis. Neither the thesis nor la these ni des extraits substantiels de celle-ci substantial extracts from it may be ne doivent etre imprimes ou autrement printed or otherwise reproduced reproduits sans son autorisation. without the author's permission. In compliance with the Canadian Conformement a la loi canadienne sur la Privacy Act some supporting forms protection de la vie privee, quelques may have been removed from this formulaires secondaires ont ete enleves de thesis. cette these. While these forms may be included Bien que ces formulaires aient inclus dans in the document page count, their la pagination, il n'y aura aucun contenu removal does not represent any loss manquant. of content from the thesis. 1*1 Canada ABSTRACT Asymptotic Bounds to Outputs of the Exact Weak Solution of the Three-Dimensional Helmholtz Equation Shahin Ghomeshi, Ph.D. Concordia University, 2010 In engineering practice, the design is based on certain design quantities or "outputs of interest" which are functionals of field variables such as displacement, velocity field, or pressure. In order to gain confidence in the numerical approximation of "outputs," a method of obtaining sharp, rigorous upper and lower bounds to outputs of the exact solution have been developed for symmetric and coercive problems (the Poisson equation and the elasticity equation), for non-symmetric coercive problems (advection-diffusion-reaction equation), and more recently for certain constrained problems (Stokes equation). In this thesis we develop the method for the Helmholtz equation. The common approach relies on decomposing the global problem into independent local elemental sub-problems by relaxing the continuity along the edges of a partition­ ing of the entire domain, using approximate Lagrange multipliers. The method exploits the Lagrangian saddle point property by recasting the output problem as a constrained minimization problem. The upper and lower computed bounds then hold for all levels of refinement and are shown to approach the exact solution at the same rate as its underlying finite element approach. The certificate of precision can then determine the best as well as the worst case scenario in an engineering design problem. This thesis addresses bounds to outputs of interest for the complex Helmholtz equation. The Helmholtz equation is in gen- iii IV eral non-coercive for high wave numbers and therefore, the previous approaches that relied on duality theory of convex minimization do not directly apply. Only in the asymptotic regime does the Helmholtz equation become coercive, and reliable (guaranteed) bounds can thus be obtained. Therefore, in order to achieve good bounds, several new ingredients have been introduced. The bounds procedure is firstly formulated with appropriate extension to complex-valued equations. Secondly, in the computation of the inter-subdomain continuity multipliers we follow the FETI-H approach and regularize the system matrix with a complex term to make the system non-singular. Finally, in order to obtain sharper output bounds in the presence of pollution errors, higher order nodal spectral element method is employed which has several computational advantages over the traditional finite element approach. We performed verification of our results and demonstrate the bounding properties for the Helmholtz problem. Acknowledgements First and foremost, I would like to thank my advisor, Professor Marius Paraschivoiu, for his continual supply of ideas, as well as for his encouragement and patience during the course of this research. Secondly, I am grateful to Professor Frederic Magoules from Ecole Centrale Paris for his feedback and insight into my work, which was not only very helpful, but also very encouraging. I would also like thank my committee members at Concordia university for their useful feedback and advice on the thesis. These include Professors Georgios Vatistas, and Ramin Sedaghati from the mechanical engineering department and Professor Tien Bui from the computer science department. I am also extremely thankful for the warm atmosphere of the CFD lab where many of my colleagues there provided me with very fruitful discussions. Finally, I could not have overcome the many frustrations and disappointments experienced during the course of my studies without the support of my family and friends. I would especially like to thank my parents Eshrat and Behrang, whose constant love and encouragement gave me the strength and encouragement when times seemed most difficult. v Table of Contents Chapter 1 Introduction 1 1.1 Acoustic Wave Propagation Problem 2 1.1.1 Model Problems in Acoustics 3 1.1.1 Interor Problems 4 1.1.1 Exterior Problems 5 1.1.2 Numerical Challenges 7 1.2 Error Estimation 7 1.3 Objectives and Scope of Thesis 13 Chapter 2 Numerical Methods for the Helmholtz Equation 15 2.1 Model Problem 15 2.1.1 Strong Formulation 15 2.1.2 Preliminaries from Functional Analysis 16 2.1.3 Notations and Weak Formulation 19 2.1.4 Linear Functional Outputs 20 2.2 Discretization Methods 21 vi Table of Contents vn 2.2.1 Nodal Points 24 2.2.2 Basis Functions 25 2.3 Some Key Concepts 28 2.3.1 Positive Definite Forms 29 2.3.2 inf — sup Condition 30 2.3.3 V-Coercive Forms 32 2.3.4 Variational Methods 33 2.3.4 Convergence Properties 35 2.3.5 An Observation 37 2.4 Pollution Effect 38 2.4.1 Piecewise linear case 38 2.4.2 Higher-Order Elements-Zip FEM 41 Chapter 3 Domain Decomposition Methods 43 3.1 The FETI Procedure - Poisson Example 45 3.1.1 Discretization 47 3.1.2 The FETI PCPG Iterative Procedure 50 3.2 Domain Decomposition for Helmholtz 51 3.2.1 Discretization 53 3.2.2 Iterative Method Used 56 3.2.2 Some Comments on Preconditioning 60 Table of Contents viii Chapter 4 Exact Bounds Method for the Poisson Equation 61 4.1 Bounds on Energy 62 4.2 Bounds on Quantitative Outputs 69 4.2.2 Lagrange Multiplier Approximation 71 4.2.3 Local Dual Sub-Problems 71 4.2.3 Sub-problem Approximation 73 4.2.4 Output Bounds 75 4.3 Implementation 76 4.3.1 Interpolation of Hybrid Fluxes 76 4.4 Discrete forms and Sub-problems Computation 78 4.4.1 Numerical Examples 80 4.4.1 Constructed Exact Solution 80 4.4.1 Uniformly Forced Domain 82 4.5 Some Results on Stokes' Problem 85 4.5.1 Model Problem 85 4.5.2 Output Functional 87 4.5.3 Numerical Example of Bound Calculation for Stokes 87 Chapter 5 Bounds for the Helmholtz Equation 91 5.1 Problem Statement 92 5.1.1 Governing Equations 92 Table of Contents IX 5.1.2 Approximation Spaces 92 5.2 Error Bound Formulation 94 5.2.1 Output Functional 94 5.2.2 Error Formulation 94 5.2.3 Lagrange Multiplier Approximation 95 5.2.3 Modified Lagrangian 95 5.2.4 Local Problems 99 5.3 Localized Lagrangian 102 5.4 An Equivalent Lagrangian 103 5.5 Elemental Sub-problems 105 5.6 Bounding Property 107 5.7 Sub-problem Computations 109 5.7.1 Discrete Forms and Approximations 110 5.8 Computation of Bounds 112 5.9 Convergence Properties 113 Chapter 6 Numerical Validation and Results 117 6.1 Discussion of Results 118 6.1.1 For p = 2, q = 3 118 6.1.1 Case I: k = 1 118 6.1.1 Case II: k = 3 118 6.1.1 Case III: k = 5 119 Table of Contents x 6.1.2 For p = 4, q = 5 120 6.1.2 Case I: k = 3 120 6.1.2 Case II: k = 5 121 Chapter 7 Conclusions and Future Work 129 7.1 Conclusions 129 7.2 Future Work 131 Bibliography 134 List of Figures Figure 1.1 Computational domain £7 for a vibro-acoustic problem 5 Figure 1.2 Radiation problem for a vibrating body D 6 Figure 3.1 Representation of the checkerboard partitioning of the mesh 55 Figure 4.1 Conforming nature of the upper bound 63 Figure 4.2 Lower bounding property 68 Figure 4.3 Hybrid flux interpolation on the faces of the subdomains 79 Figure 4.4 Output bounds obtained for constructed solution 83 Figure 4.5 Output bounds obtained for constant forcing, p = 1, q = 2 84 Figure 4.6 Output bounds obtained for constant forcing, p = 2, q = 3 84 Figure 4.7 Geometry for the Stokes problem 86 Figure 4.8 Bounds for the Stokes output 88 Figure 6.1 Convergence study for k = 1 123 xi List of Figures xn Figure 6.2 Bounding property of the bounds for k = 1 123 Figure 6.3 Convergence study for k = 3 124 Figure 6.4 Upper and lower bounds for k = 3 124 Figure 6.5 Convergence study for k = 5 125 Figure 6.6 Upper and lower bounds for k — 5 125 Figure 6.7 Convergence study for k — 3 using higher-order elements 126 Figure 6.8 Upper and lower bounds for k = 3 using higher-order elements.
Load more

Recommended publications
  • Efficient Spectral-Galerkin Method I. Direct Solvers for the Second And
    SIAM J. SCI. COMPUT. °c 1994 Society for Industrial and Applied Mathematics Vol. 15, No. 6, pp. 1489{1505, November 1994 013 E±cient Spectral-Galerkin Method I. Direct Solvers for the Second and Fourth Order Equations Using Legendre Polynomials¤ Jie Sheny Abstract. We present some e±cient algorithms based on the Legendre-Galerkin approxima- tions for the direct solution of the second and fourth order elliptic equations. The key to the e±ciency of our algorithms is to construct appropriate base functions, which lead to systems with sparse matri- ces for the discrete variational formulations. The complexities of the algorithms are a small multiple of N d+1 operations for a d dimensional domain with (N ¡ 1)d unknowns, while the convergence rates of the algorithms are exponential for problems with smooth solutions. In addition, the algorithms can be e®ectively parallelized since the bottlenecks of the algorithms are matrix-matrix multiplications. Key words. spectral-Galerkin method, Legendre polynomial, Helmholtz equation, biharmonic equation, direct solver AMS subject classi¯cations. 65N35, 65N22, 65F05, 35J05 1. Introduction. This article is the ¯rst in a series for developing e±cient spectral Galerkin algorithms for elliptic problems. The spectral method employs global polynomials as the trial functions for the discretization of partial di®erential equations. It provides very accurate approximations with a relatively small number of unknowns. Consequently it has gained increasing popularity in the last two decades, especially in the ¯eld of computational fluid dynamics (see [11], [8] and the references therein). The use of di®erent test functions in a variational formulation leads to three most commonly used spectral schemes, namely, the Galerkin, tau and collocation versions.
    [Show full text]
  • Family Name Given Name Presentation Title Session Code
    Family Name Given Name Presentation Title Session Code Abdoulaev Gassan Solving Optical Tomography Problem Using PDE-Constrained Optimization Method Poster P Acebron Juan Domain Decomposition Solution of Elliptic Boundary Value Problems via Monte Carlo and Quasi-Monte Carlo Methods Formulations2 C10 Adams Mark Ultrascalable Algebraic Multigrid Methods with Applications to Whole Bone Micro-Mechanics Problems Multigrid C7 Aitbayev Rakhim Convergence Analysis and Multilevel Preconditioners for a Quadrature Galerkin Approximation of a Biharmonic Problem Fourth-order & ElasticityC8 Anthonissen Martijn Convergence Analysis of the Local Defect Correction Method for 2D Convection-diffusion Equations Flows C3 Bacuta Constantin Partition of Unity Method on Nonmatching Grids for the Stokes Equations Applications1 C9 Bal Guillaume Some Convergence Results for the Parareal Algorithm Space-Time ParallelM5 Bank Randolph A Domain Decomposition Solver for a Parallel Adaptive Meshing Paradigm Plenary I6 Barbateu Mikael Construction of the Balancing Domain Decomposition Preconditioner for Nonlinear Elastodynamic Problems Balancing & FETIC4 Bavestrello Henri On Two Extensions of the FETI-DP Method to Constrained Linear Problems FETI & Neumann-NeumannM7 Berninger Heiko On Nonlinear Domain Decomposition Methods for Jumping Nonlinearities Heterogeneities C2 Bertoluzza Silvia The Fully Discrete Fat Boundary Method: Optimal Error Estimates Formulations2 C10 Biros George A Survey of Multilevel and Domain Decomposition Preconditioners for Inverse Problems in Time-dependent
    [Show full text]
  • A Space-Time Discontinuous Galerkin Spectral Element Method for the Stefan Problem
    Noname manuscript No. (will be inserted by the editor) A Space-Time Discontinuous Galerkin Spectral Element Method for the Stefan problem Chaoxu Pei · Mark Sussman · M. Yousuff Hussaini Received: / Accepted: Abstract A space-time Discontinuous Galerkin (DG) spectral element method is presented to solve the Stefan problem in an Eulerian coordinate system. In this scheme, the level set method describes the time evolving interface. To deal with the prior unknown interface, two global transformations, a backward transformation and a forward transformation, are introduced in the space-time mesh. The idea of the transformations is to combine an Eulerian description, i.e. a fixed frame of reference, with a Lagrangian description, i.e. a moving frame of reference. The backward transformation maps the unknown time- varying interface in the fixed frame of reference to a known stationary inter- face in the moving frame of reference. In the moving frame of reference, the transformed governing equations, written in the space-time framework, are discretized by a DG spectral element method in each space-time slab. The forward transformation is used to update the level set function and then to project the solution in each phase onto the new corresponding time-dependent domain. Two options for calculating the interface velocity are presented, and both methods exhibit spectral accuracy. Benchmark tests in one spatial di- mension indicate that the method converges with spectral accuracy in both space and time for both the temperature distribution and the interface ve- locity. The interrelation between the interface position and the temperature makes the Stefan problem a nonlinear problem; a Picard iteration algorithm is introduced in order to solve the nonlinear algebraic system of equations and it is found that only a few iterations lead to convergence.
    [Show full text]
  • Spectral Element Simulations of Laminar and Turbulent Flows in Complex Geometries
    Applied Numerical Mathematics 6 (1989/90) 85-105 85 North-Holland SPECTRAL ELEMENT SIMULATIONS OF LAMINAR AND TURBULENT FLOWS IN COMPLEX GEOMETRIES George Em KARNIADAKIS* Center for Turbulence Research, Stanford University/NASA Ames Research Center, Stanford, CA 94305, USA Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes and transport equations . Their success in the recent past in simulating flows of industrial complexity derives from the flexibility of the method to accurately represent nontrivial geometries while preserving the good resolution properties of the spectral methods . In this paper, we review some of the main ideas of the method with emphasis placed on implementation and data management . These issues need special attention in order to make the method efficient in practice, especially in view of the fact that high computing cost as well as strenuous storage requirements have been a major drawback of high-order methods in the past . Several unsteady, laminar complex flows are simulated, and a direct simulation of turbulent channel flow is presented, for the first time, using spectral element techniques . 1 . Introduction Spectral methods have proven in recent years a very powerful tool for analyzing fluid flows and have been used almost exclusively in direct simulations of transitional and turbulent flows [8]. To date, however, only simple geometry turbulent flows have been studied accurately [14], chiefly because of the poor performance of global spectral methods in representing more complex geometries. The development of the spectral element method [26] and its success to accurately simulate highly unsteady, two-dimensional laminar flows [6,13] suggests that a three-dimensional implementation may be the right approach to direct simulations of turbulent flows in truly complex geometries.
    [Show full text]
  • Hybrid Multigrid/Schwarz Algorithms for the Spectral Element Method
    Hybrid Multigrid/Schwarz Algorithms for the Spectral Element Method James W. Lottes¤ and Paul F. Fischery February 4, 2004 Abstract We study the performance of the multigrid method applied to spectral element (SE) discretizations of the Poisson and Helmholtz equations. Smoothers based on finite element (FE) discretizations, overlapping Schwarz methods, and point-Jacobi are con- sidered in conjunction with conjugate gradient and GMRES acceleration techniques. It is found that Schwarz methods based on restrictions of the originating SE matrices converge faster than FE-based methods and that weighting the Schwarz matrices by the inverse of the diagonal counting matrix is essential to effective Schwarz smoothing. Sev- eral of the methods considered achieve convergence rates comparable to those attained by classic multigrid on regular grids. 1 Introduction The availability of fast elliptic solvers is essential to many areas of scientific computing. For unstructured discretizations in three dimensions, iterative solvers are generally optimal from both work and storage standpoints. Ideally, one would like to have computational complexity that scales as O(n) for an n-point grid problem in lRd, implying that the it- eration count should be bounded as the mesh is refined. Modern iterative methods such as multigrid and Schwarz-based domain decomposition achieve bounded iteration counts through the introduction of multiple representations of the solution (or the residual) that allow efficient elimination of the error at each scale. The theory for these methods is well established for classical finite difference (FD) and finite element (FE) discretizations, and order-independent convergence rates are often attained in practice. For spectral element (SE) methods, there has been significant work on the development of Schwarz-based methods that employ a combination of local subdomain solves and sparse global solves to precondition conjugate gradient iteration.
    [Show full text]
  • SCHWARZ Rr'"YPE DOMAIN DECOMPOSITION ]\1ETHODS FOR
    SCHWARZ rr'"YPE DOMAIN DECOMPOSITION ]\1ETHODS FOR St<EC'l'RAL ELEMENT DISCRETIZATIONS Shannon S. Pahl University of the Witwatersrand December 1993 .() Degree awarded with distinction on 30 June 1994. , o A research report submitted to the Faculty of Science in partial fulfillment of the requirements 'L for the degree of Mastel' of Science at the University of the Witwatersrand, Johannesburg. () II " Abstract Most cf the theory for domain decomposition methods of Schwarz tyne has been set in the framework of the h ~',ndp~version finite element method. In this study, Schwarz methods are formulated for both the conforming and nonconforming spectral element discretizations applied to linear scalar self adjoint second order elliptic problems in two dimensions. An overlapping additive Schwarz method for the conforming spectral element discretization is formulated. However, unlike p~nnite element discretizations, a minimum overlap strategy com- mon to h-J11niteelement discret)zations can be accommodated. Computational results indicate if that the convergence rate of the overlapping Schwarz method is similar for the corresponding method defined for the h-finite element discretization. It is also shown that the minimum over- lap additive Schwarz method results in significantly fewer floating point operations than finite element preconditioning for the spectral element method. Iterative substructuring and Neumann-Neumann methods carealso formulated for the conforming spectral element discretization. These methods demonstrate greater robustness when increasing the degree p of the elements than the minimum overlap additive Schwarz method. Computa- tionally, convergencerates are also similar for the corresponding methods defined for the h-finite element discretization. In addition, an efficient interface preconditioner for iterative substruc- turing methods is employed which does not degrade the performance of the method for tY..e problems considered.
    [Show full text]
  • NONCONFORMING FORMULATIONS with SPECTRAL ELEMENT METHODS a Dissertation by CUNEYT SERT Submitted to the Office of Graduate Studi
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Texas A&amp;M Repository NONCONFORMING FORMULATIONS WITH SPECTRAL ELEMENT METHODS A Dissertation by CUNEYT SERT Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY August 2003 Major Subject: Mechanical Engineering NONCONFORMING FORMULATIONS WITH SPECTRAL ELEMENT METHODS A Dissertation by CUNEYT SERT Submitted to Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Approved as to style and content by: Ali Beskok (Chair of Committee) J. N. Reddy N. K. Anand (Member) (Member) Paul Cizmas Dennis O’Neal (Member) (Head of Department) August 2003 Major Subject: Mechanical Engineering iii ABSTRACT Nonconforming Formulations with Spectral Element Methods. (August 2003) Cuneyt Sert, B.S., Middle East Technical University; M.S., Middle East Technical University Chair of Advisory Committee: Dr. Ali Beskok A spectral element algorithm for solution of the incompressible Navier-Stokes and heat transfer equations is developed, with an emphasis on extending the classical conforming Galerkin formulations to nonconforming spectral elements. The new algorithm employs both the Constrained Approximation Method (CAM), and the Mortar Element Method (MEM) for p-and h-type nonconforming elements. Detailed descriptions, and formulation steps for both methods, as well as the performance com- parisons between CAM and MEM, are presented. This study fills an important gap in the literature by providing a detailed explanation for treatment of p-and h-type nonconforming interfaces. A comparative eigenvalue spectrum analysis of diffusion and convection operators is provided for CAM and MEM.
    [Show full text]
  • Fekete-Gauss Spectral Elements for Incompressible Navier-Stokes Flows: the Two-Dimensional Case 1 1, 1 Laura Lazar , Richard Pasquetti ∗ and Francesca Rapetti 1 Lab
    Commun. Comput. Phys. Vol. 13, No. 5, pp. 1309-1329 doi: 10.4208/cicp.180112.110612a May 2013 Fekete-Gauss Spectral Elements for Incompressible Navier-Stokes Flows: The Two-Dimensional Case 1 1, 1 Laura Lazar , Richard Pasquetti ∗ and Francesca Rapetti 1 Lab. J.A. Dieudonn´e, UMR 7351 CNRS UNS, Universit´ede Nice - Sophia Antipolis, 06108 Nice Cedex 02, France. Received 18 January 2012; Accepted (in revised version) 11 June 2012 Available online 8 October 2012 Abstract. Spectral element methods on simplicial meshes, say TSEM, show both the advantages of spectral and finite element methods, i.e., spectral accuracy and geomet- rical flexibility. We present a TSEM solver of the two-dimensional (2D) incompressible Navier-Stokes equations, with possible extension to the 3D case. It uses a projection method in time and piecewise polynomial basis functions of arbitrary degree in space. The so-called Fekete-Gauss TSEM is employed, i.e., Fekete (resp. Gauss) points of the triangle are used as interpolation (resp. quadrature) points. For the sake of consistency, isoparametric elements are used to approximate curved geometries. The resolution al- gorithm is based on an efficient Schur complement method, so that one only solves for the element boundary nodes. Moreover, the algebraic system is never assembled, therefore the number of degrees of freedom is not limiting. An accuracy study is car- ried out and results are provided for classical benchmarks: the driven cavity flow, the flow between eccentric cylinders and the flow past a cylinder. AMS subject classifications: 65M60, 65M70, 76D05 Key words: Spectral elements, simplicial meshes, Fekete-Gauss approximation, Navier-Stokes equations, projection methods, domain decomposition.
    [Show full text]
  • Additive Schwarz Preconditioners with Minimal Overlap For
    Commun. Comput. Phys. Vol. 13, No. 2, pp. 411-427 doi: 10.4208/cicp.080711.160212a February 2013 Additive Schwarz Preconditioners with Minimal Overlap for Triangular Spectral Elements Yuen-Yick Kwan∗ Center for Computational Science, Tulane University, New Orleans, LA 70118, USA. Received 8 July 2011; Accepted (in revised version) 16 February 2012 Communicated by Jie Shen Available online 28 June 2012 Abstract. The additive Schwarz preconditioner with minimal overlap is extended to triangular spectral elements (TSEM). The method is a generalization of the correspond- ing method in tensorial quadrilateral spectral elements (QSEM). The proposed precon- ditioners are based on partitioning the domain into overlapping subdomains, solving local problems on these subdomains and solving an additional coarse problem asso- ciated with the subdomain mesh. The results of numerical experiments show that the proposed preconditioner are robust with respect to the number of elements and are more efficient than the preconditioners with generous overlaps. AMS subject classifications: 65F08, 65N35 Key words: Spectral elements, additive Schwarz. 1 Introduction Spectral element method, which combines the flexibility of the low-order finite element methods and the high accuracy of the spectral methods, are popular in computational science and engineering. The efficiency of the spectral element method depends on the solution method employed to solve the resultant linear system. Since direct methods are infeasible when the number of elements get large, preconditioned iterative methods are usually used. Preconditioners can be based on an overlapping or non-overlapping domain decomposition [24, 28]. The former include the multiplicative Schwarz [18], ad- ditive Schwarz [9] and the restricted Schwarz preconditioners [6].
    [Show full text]
  • The Spectral Element Method (Sem): Formulation and Implementation of the Method for Engineering Seismology Problems
    th The 14 World Conference on Earthquake Engineering October 12-17, 2008, Beijing, China THE SPECTRAL ELEMENT METHOD (SEM): FORMULATION AND IMPLEMENTATION OF THE METHOD FOR ENGINEERING SEISMOLOGY PROBLEMS 1 2 Kristel C. Meza-Fajardo and Apostolos S. Papageorgiou 1Professor, Depto. de Ingenieria Civil,Universidad Nacional Autónoma de Honduras,Tegucigalpa, Honduras 2 Professor, Department of Civil Engineering, University of Patras, GR-26500 Patras, Greece Email: [email protected], [email protected] ABSTRACT Various numerical methods have been proposed and used to investigate wave propagation in realistic earth media. Recently an innovative numerical method, known as the Spectral Element Method (SEM), has been developed and used in connection with wave propagation problems in 3D elastic media (Komatitsch and Tromp, 1999). The SEM combines the flexibility of a finite element method with the accuracy of a spectral method and it is not cumbersome in dealing with non-flat free surface and spatially variable anelastic attenuation. The SEM is a highly accurate numerical method that has its origins in computational fluid dynamics. One uses a weak formulation of the equations of motion, which are solved on a mesh of hexahedral elements that is adapted to the free surface and to the main internal discontinuities of the model. The wavefield on the elements is discretized using high-degree Lagrange interpolants, and integration over an element is accomplished based upon the Gauss-Lobatto-Legendre integration rule. This combination of discretization and integration results in a diagonal mass matrix, which greatly simplifies the algorithm. The most important property of the SEM is that the mass matrix is exactly diagonal by construction, which drastically simplifies the implementation and reduces the computational cost because one can use an explicit time integration scheme without having to invert a linear system.
    [Show full text]
  • Fast Iterative Solver for Convection-Diffusion Systems with Spectral Elements P
    Fast Iterative Solver for Convection-Diffusion Systems with Spectral Elements P. Aaron Lott,1 Howard Elman2 1Mathematical and Computational Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 2Department of Computer Science and Institute for Advanced Computer Studies, University of Maryland, College Park, Maryland 20742 Received 9 March 2009; accepted 12 July 2009 Published online 16 October 2009 in Wiley Online Library (wileyonlinelibrary.com). DOI 10.1002/num.20518 We introduce a solver and preconditioning technique based on Domain Decomposition and the Fast Diagonal- ization Method that can be applied to tensor product based discretizations of the steady convection–diffusion equation. The method is based on a Robin–Robin interface preconditioner coupled to a fast diagonalization solver which is used to efficiently eliminate the interior degrees of freedom and perform subsidiary subdo- main solves. Using a spectral element discretization, we first apply our technique to constant wind problems, and then propose a means for applying the technique as a preconditioner for variable wind problems. We demonstrate that iteration counts are mildly dependent on changes in mesh size and convection strength. © 2009 Wiley Periodicals, Inc.∗ Numer Methods Partial Differential Eq 27: 231–254, 2011 Keywords: convection-diffusion; domain decomposition; preconditioning; spectral element method I. INTRODUCTION Numerical simulation of fluid flow allows for improved prediction and design of natural and engi- neered systems such as those involving water, oil, or blood. The interplay between inertial and viscous forces in a fluid flow dictates the length scale where energy is transferred, thus determin- ing the resolution required to capture flow information accurately.
    [Show full text]
  • Hybrid Schwarz-Multigrid Methods for the Spectral Element Method: Extensions to Navier-Stokes
    Hybrid Schwarz-Multigrid Methods for the Spectral Element Method: Extensions to Navier-Stokes Paul F. Fischer1 and James W. Lottes2 1 Argonne National Laboratory, Mathematics and Computer Science Division (http://www.mcs.anl.gov/~fischer/) 2 University of Illinois, Dept. of Theoretical and Applied Mechanics Summary. We study the performance of the multigrid method applied to spectral element (SE) discretizations of the Poisson and Helmholtz equations. Smoothers based on finite element (FE) discretizations, overlapping Schwarz methods, and point-Jacobi are considered in conjunction with conjugate gradient and GMRES acceleration techniques. It is found that Schwarz methods based on restrictions of the originating SE matrices converge faster than FE-based methods and that weighting the Schwarz matrices by the inverse of the diagonal counting matrix is essential to effective Schwarz smoothing. Several of the methods considered achieve convergence rates comparable to those attained by classic multigrid on regular grids. 1 Introduction The spectral element method (SEM) is a high-order weighted residual tech- nique that combines the geometric flexibility of finite elements with the rapid convergence properties and tensor-product efficiencies of global spectral meth- ods. Globally, elements are coupled in an unstructured framework with in- terelement coupling enforced through standard matching of nodal interface values. Locally, functions are represented as tensor products of stable Nth- order Lagrangian interpolants based on Gauss-Lobatto (GL) or Gauss (G) quadrature points. For problems having smooth solutions, such as the incom- pressible Navier-Stokes equations, the SEM converges exponentially fast with the local approximation order N. Because of its minimal numerical dissipa- tion and dispersion, the SEM is particularly well suited for the simulation of flows at transitional Reynolds numbers, where physical dissipation is small and turbulence-model dissipation is absent.
    [Show full text]