skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems

Authors:
ORCiD logo [1];  [1];  [2]
  1. Department of Applied Mathematics, University of Colorado at Boulder, USA
  2. Department of Applied Mathematics, University of Colorado at Boulder, USA, Fakultät für Mathematik, Universität Duisburg-Essen, Germany
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1417505
Grant/Contract Number:
FC02-03ER25574; NA0002376; B614452
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Numerical Linear Algebra with Applications
Additional Journal Information:
Journal Volume: 25; Journal Issue: 3; Related Information: CHORUS Timestamp: 2018-04-09 16:06:04; Journal ID: ISSN 1070-5325
Publisher:
Wiley Blackwell (John Wiley & Sons)
Country of Publication:
United Kingdom
Language:
English

Citation Formats

Kalchev, Delyan Z., Manteuffel, Thomas A., and Münzenmaier, Steffen. Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems. United Kingdom: N. p., 2018. Web. doi:10.1002/nla.2150.
Kalchev, Delyan Z., Manteuffel, Thomas A., & Münzenmaier, Steffen. Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems. United Kingdom. doi:10.1002/nla.2150.
Kalchev, Delyan Z., Manteuffel, Thomas A., and Münzenmaier, Steffen. Fri . "Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems". United Kingdom. doi:10.1002/nla.2150.
@article{osti_1417505,
title = {Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems: Mixed (LL∗)−1 and LL∗ least-squares finite element methods with application to linear hyperbolic problems},
author = {Kalchev, Delyan Z. and Manteuffel, Thomas A. and Münzenmaier, Steffen},
abstractNote = {},
doi = {10.1002/nla.2150},
journal = {Numerical Linear Algebra with Applications},
number = 3,
volume = 25,
place = {United Kingdom},
year = {Fri Jan 19 00:00:00 EST 2018},
month = {Fri Jan 19 00:00:00 EST 2018}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on January 19, 2019
Publisher's Accepted Manuscript

Save / Share:
  • No abstract prepared.
  • In this paper we first introduce the least-squares finite element method (LSFEM) for two-dimensional steady-state pure convection problems with smooth solutions and compare the LSFEM with other finite element methods. We prove that the LSFEM has the same stability estimate as the original equation; i.e., the LSFEM has better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method. Numerical convergence rates are given to show that the LSFEM is almost optimal. Then we use this LSFEM as a framework to develop an iteratively reweighted least-squares finite element method (IRLSFEM) to obtain non-oscillatory and non-diffusive solutions for problems withmore » contact discontinuities. This new method produces a highly accurate numerical solution that has a sharp discontinuity in one element. A number of examples solved by using triangular and bilinear elements are presented to show that the method can convect contact discontinuities without error.« less
  • Least-squares finite-element methods for Darcy flow offer several advantages relative to the mixed-Galerkin method: the avoidance of stability conditions between finite-element spaces, the efficiency of solving symmetric and positive definite systems, and the convenience of using standard, continuous nodal elements for all variables. However, conventional C{sup o} implementations conserve mass only approximately and for this reason they have found limited acceptance in applications where locally conservative velocity fields are of primary interest. In this paper, we show that a properly formulated compatible least-squares method offers the same level of local conservation as a mixed method. The price paid for gainingmore » favourable conservation properties is that one has to give up what is arguably the least important advantage attributed to least-squares finite-element methods: one can no longer use continuous nodal elements for all variables. As an added benefit, compatible least-squares methods inherit the best computational properties of both Galerkin and mixed-Galerkin methods and, in some cases, yield identical results, while offering the advantages of not having to deal with stability conditions and yielding positive definite discrete problems. Numerical results that illustrate our findings are provided.« less
  • A significant amount of the computational time in large Monte Carlo simulations of lattice quantum chromodynamics (QCD) is spent inverting the discrete Dirac operator. Unfortunately, traditional covariant finite difference discretizations of the Dirac operator present serious challenges for standard iterative methods. For interesting physical parameters, the discretized operator is large and ill-conditioned, and has random coefficients. More recently, adaptive algebraic multigrid (AMG) methods have been shown to be effective preconditioners for Wilson's discretization of the Dirac equation. This paper presents an alternate discretization of the Dirac operator based on least-squares finite elements. The discretization is systematically developed and physical propertiesmore » of the resulting matrix system are discussed. Finally, numerical experiments are presented that demonstrate the effectiveness of adaptive smoothed aggregation ({alpha}SA ) multigrid as a preconditioner for the discrete field equations resulting from applying the proposed least-squares FE formulation to a simplified test problem, the 2d Schwinger model of quantum electrodynamics.« less
  • Abstract not provided.