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

Title: A simple finite element method for non-divergence form elliptic equation

Abstract

Here, we develop a simple finite element method for solving second order elliptic equations in non-divergence form by combining least squares concept with discontinuous approximations. This simple method has a symmetric and positive definite system and can be easily analyzed and implemented. We could have also used general meshes with polytopal element and hanging node in the method. We prove that our finite element solution approaches to the true solution when the mesh size approaches to zero. Numerical examples are tested that demonstrate the robustness and flexibility of the method.

Authors:
 [1];  [2]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Computer Science and Mathematics Division
  2. Univ. of Arkansas, Little Rock, AR (United States). Dept. of Mathematics
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)
OSTI Identifier:
1346680
Grant/Contract Number:
AC05-00OR22725
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
International Journal of Numerical Analysis and Modeling
Additional Journal Information:
Journal Volume: 14; Journal Issue: 2; Journal ID: ISSN 1705-5105
Publisher:
Institute for Scientific Computing and Information
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; finite element methods; non-divergence form elliptic equations; polyhedral meshes

Citation Formats

Mu, Lin, and Ye, Xiu. A simple finite element method for non-divergence form elliptic equation. United States: N. p., 2017. Web.
Mu, Lin, & Ye, Xiu. A simple finite element method for non-divergence form elliptic equation. United States.
Mu, Lin, and Ye, Xiu. Wed . "A simple finite element method for non-divergence form elliptic equation". United States. doi:. https://www.osti.gov/servlets/purl/1346680.
@article{osti_1346680,
title = {A simple finite element method for non-divergence form elliptic equation},
author = {Mu, Lin and Ye, Xiu},
abstractNote = {Here, we develop a simple finite element method for solving second order elliptic equations in non-divergence form by combining least squares concept with discontinuous approximations. This simple method has a symmetric and positive definite system and can be easily analyzed and implemented. We could have also used general meshes with polytopal element and hanging node in the method. We prove that our finite element solution approaches to the true solution when the mesh size approaches to zero. Numerical examples are tested that demonstrate the robustness and flexibility of the method.},
doi = {},
journal = {International Journal of Numerical Analysis and Modeling},
number = 2,
volume = 14,
place = {United States},
year = {Wed Mar 01 00:00:00 EST 2017},
month = {Wed Mar 01 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
The DOI is not currently available

Save / Share:
  • We study a general class of quasilinear non-uniformly elliptic pdes in divergence form with linear growth in the gradient. We examine the notions of BV and viscosity solutions and derive for such generalized solutions various a priori pointwise and integral estimates, including a Harnack inequality. In particular we prove that viscosity solutions are unique (on strictly convex domains), are contained in the space BV{sub loc} and are C{sup 1,{alpha}} almost everywhere. 10 refs.
  • Let L be a linear uniformly elliptic operator of the second order in R{sup n}, n>=2, with bounded measurable real coefficients, that satisfies the weak uniqueness property. The removability of compact subsets of a domain D subset of R{sup n} is studied for weak solutions of the equation Lf=0 (in the sense of Krylov and Safonov) in some classes of continuous functions in D. In particular, a metric criterion for removability in Hoelder classes with small exponent of smoothness is obtained. Bibliography: 20 titles.
  • In this paper, we study a coupling of spectral methods and the [rho]-version finite element methods for elliptic boundary value problems containing singularities. The method of auxiliary mapping, which is a recent development to deal with domain singularities in the [rho]-version of the finite element method, is employed to remove the pollution effect caused by singularities. An iterative interfacial coupling between spectral methods and the [rho]-version of the finite element method is used and investigated numerically. The advantages of such an approach are demonstrated by the high accuracy of spectral methods for the smooth part of solutions and the flexibilitymore » of the [rho]-version of the finite element method for dealing with singularities and irregular domains. The efficiency of the coupling method is also evaluated by comparing results obtained by this method with those obtained by the full finite element algorithm. 26 refs., 9 figs., 6 tabs.« less
  • Babuska and Oh have introduced a new approach called the method of auxiliary mapping (MAM), to deal with elliptic boundary value problems with singularities, in the context of the [rho]-version of the finite element method, MAM yielded an exponential rate of convergence at virtually no extra cost. In this paper those results are extended by showing the exponential convergence of MAM for homogeneous Laplace equations with boundary singularities and Helmholtz equations with both corner and boundary singularities. In addition a convergence result is developed for MAM as applied to the h-[rho]-version of the finite element method. To clarify the powermore » of MAM a series of benchmark runs are made for three examples using our implementation, MAPFEM. Comparisons are made with two h-version finite element codes (PLTMG6 and FESOP), both of which use adaptive meshes, the finite difference code ELLPACK, and the recent singular element code ISBFM. The examples include the well-known Motz problem and both homogeneous and nonhomogeneous Helmholtz equations. 41 refs., 12 figs., 8 tabs.« less
  • In this paper, we study a multiscale finite element method for solving a class of elliptic problems arising from composite materials and flows in porous media, which contain many spatial scales. The method is designed to efficiently capture the large scale behavior of the solution without resolving all the small scale features. This is accomplished by constructing the multiscale finite element base functions that are adaptive to the local property of the differential operator. Our method is applicable to general multiple-scale problems without restrictive assumptions. The construction of the base functions is fully decoupled from element to element; thus, themore » method is perfectly parallel and is naturally adapted to massively parallel computers. For the same reason, the method has the ability to handle extremely large degrees of freedom due to highly heterogeneous media, which are intractable by conventional finite element (difference) methods. In contrast to some empirical numerical upscaling methods, the multiscale method is systematic and self-consistent, which makes it easier to analyze. We give a brief analysis of the method, with emphasis on the {open_quotes}resonant sampling{close_quotes} effect. Then, we propose an oversampling technique to remove the resonance effect. We demonstrate the accuracy and efficiency of our method through extensive numerical experiments, which include problems with random coefficients and problems with continuous scales. Parallel implementation and performance of the method are also addressed. 27 refs., 5 figs., 4 tabs.« less