Show simple item record

dc.contributor.advisorCastillo, Paul E.
dc.contributor.authorVelazquez-Suarez, Esov S.
dc.description.abstractIn this work, we consider second order elliptic problems arising in the modeling of single phase flows in porous media in 2D and in the analysis of transverse electromagnetic modes in wave guides using a discontinuous Galerkin (DG) method, the so-called Local Discontinuous Galerkin (LDG) method. We designed and developed an object oriented framework for performing DG computations on unstructured meshes that allows the use of arbitrary degree in the polynomial approximations and non conformal meshes with an arbitrary number of hanging nodes per edge. We present numerical studies of an automatic mesh adaptation technique and a semi-algebraic multilevel preconditioner for the LDG method. DG methods may be viewed as high-order extensions of the classical finite volume method. Since no inter-element continuity is imposed, they can be defined on very general meshes, including non-conforming meshes, making these methods suitable for h-adaptivity. Our adaptive algorithm starts with an initial, conformal spatial discretization of the domain where the numerical solution of the partial differential equation is obtained using the LDG method. In each step, the error of the solution is estimated and the mesh is modified successively by performing two local operations: refining a fraction of the cells where the estimated error is greater and agglomerating a fraction of the cells where the estimated error is smallest. This procedure is repeated until the solution reaches a desired accuracy. It has been recently shown that the spectral condition number of the stiffness matrix exhibits an asymptotic behavior of O(h?2) on structured and unstructured meshes, where h is the mesh size, making the use of effective preconditioners a practical requirement. We present a semi-algebraic multilevel preconditioner for the LDG method and show through several numerical experiments that its performance does not degrade as the number of unknowns augments. The performance of these techniques is explored on problems with high jumps in the coefficients, which is the typical scenario of problems arising in practical applications.en_US
dc.subjectGalerkin methodsen_US
dc.titleAn object-oriented framework for hp-adaptive discontinuous galerkin methodsen_US
dc.rights.licenseAll rights reserveden_US
dc.rights.holder(c) 2006 Esov S. Velazquez Suarezen_US
dc.contributor.committeeIerkic, Henrick Mario
dc.contributor.committeeRamos, Rafael A.
dc.contributor.committeeSteinberg, Lev
dc.contributor.representativeAlemar, Jose D. and Information Sciences and Engineeringen_US
dc.contributor.collegeCollege of Engineeringen_US
dc.contributor.departmentDepartment of Electrical and Computer Engineeringen_US

Files in this item


This item appears in the following Collection(s)

  • Theses & Dissertations
    Items included under this collection are theses, dissertations, and project reports submitted as a requirement for completing a graduate degree at UPR-Mayagüez.

Show simple item record

All rights reserved
Except where otherwise noted, this item's license is described as All Rights Reserved