Search NASA⌕ Search

SEARCH · Search NASA

Results for “elliptic problems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Annotations on the virtual element method for second-order elliptic problems

This document contains working annotations on the Virtual Element Method (VEM) for the approximate solution of diffusion problems with variable coefficients. To read this document you are assumed to have familiarity with concepts from the numerical discretization of Partial Differential Equations (PDEs) and, in particular, the Finite Element Method (FEM). This document is not an introduction to the FEM, for which many textbooks (also free on the internet) are available. Eventually, this document is intended to evolve into a tutorial introduction to the VEM (but this is really a long-term goal).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Programming substructure computations for elliptic problems on the CHiP system

A number of studies have been conducted with the aim to apply parallel computation to problems associated with solving finite element equations arising in structural mechanics and fluid dynamics. These studies have provided many important results. The present investigation is concerned with a set of experiments designed to test two ideas, including configurability and substructuring. The considered algorithms and tests are intended for implementation on the Configurable, Highly Parallel (CHiP) family of architecture described by Snyder (1982). The ChiP computer is composed of homogeneous processing elements (PEs) placed at regular intervals in a lattice of programmable switches. Two examples of the role of configurability and substructuring for simple iterative algorithms are considered, giving attention to conjugate gradient iterations, and tridiagonal systems of equations.

Gannon, D.↗
NASA NTRS1996

Multiple Coarse Grid Multigrid Methods for Solving Elliptic Problems

In this paper we describe some classes of multigrid methods for solving large linear systems arising in the solution by finite difference methods of certain boundary value problems involving Poisson's equation on rectangular regions. If parallel computing systems are used, then with standard multigrid methods many of the processors will be idle when one is working at the coarsest grid levels. We describe the use of Multiple Coarse Grid MultiGrid (MCGMG) methods. Here one first constructs a periodic set of equations corresponding to the given system. One then constructs a set of coarse grids such that for each grid corresponding to the grid size h there are four grids corresponding to the grid size 2*h. Multigrid operations such as restriction of residuals and interpolation of corrections are done in parallel at each grid level. For suitable choices of the multigrid operators the MCGMG method is equivalent to the Parallel Superconvergent MultiGrid (PSMG) method of Frederickson and McBryan. The convergence properties of MCGMG methods can be accurately analyzed using spectral methods.

Xiao, Shengyou↗
NASA NTRS1985

Absorbing boundary conditions for exterior problems

Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary condition at infinity and be accurate with respect to the interior numerical scheme. Such boundary conditions are commonly referred to as absorbing boundary conditions. A treatment is given of these boundary conditions for wave-like equations.

Hariharan, S. I.↗
NASA NTRS1985

Nonlinear acoustic wave propagation in atmosphere. Absorbing boundary conditions for exterior problems

Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary conditions at infinity and be accurate with respect to the interior numerical scheme. The treatment of these boundary conditions for wave-like equations is discussed.

Hariharan, S. I.↗
NASA NTRS1982

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies is introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The lineaar stabilities of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.

Kwok, J. H.↗
DOE OSTI2018

Operator-Norm Convergence Estimates for Elliptic Homogenization Problems on Periodic Singular Structures

For an arbitrary periodic Borel measure μ we prove order O(ε) operator-norm resolvent estimates for the solutions to scalar elliptic problems in L{sup 2}(ℝ{sup d}, dμ{sup ε}) with ε-periodic coefficients, ε > 0. Here, μ{sup ε} is the measure obtained by ε-scaling of μ. Our analysis includes the case of a measure absolutely continuous with respect to the standard Lebesgue measure, as well as the case of “singular” periodic structures (or “multistructures”), when μ is supported by lower-dimensional manifolds.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

A uniformly valid solution for motion about the interior libration point of the perturbed elliptic-restricted problem

Bounded motion about Lagrange collinear libration points is considered for a perturbed elliptic-restricted problem. A practical application is the motion of a satellite near a libration point collinear with the sun and the earth-moon barycenter. Such a system is treated here as an earth-sun-satellite elliptic restricted problem with lunar perturbations. The method of dual time scales is utilized to develop a uniformly valid three-dimensional analytical solution to the satellite's equations of motion. The analytical development applies somewhat generally to that class of four-body problems where the second primary mass is much greater than the first, and the third primary mass is much greater than the second.

Richardson, D. L.↗
DOE OSTI2026

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A new weak Galerkin finite element method for elliptic interface problems

We introduce and analyze a new weak Galerkin (WG) finite element method in this paper for solving second order elliptic equations with discontinuous coefficients and interfaces. Comparing with the existing WG algorithm for solving the same type problems, the present WG method has a simpler variational formulation and fewer unknowns. Moreover, the new WG algorithm allows the use of finite element partitions consisting of general polytopal meshes and can be easily generalized to high orders. Optimal order error estimates in both H1 and L2 norms are established for the present WG finite element solutions. We conducted extensive numerical experiments in order to examine the accuracy, flexibility, and robustness of the proposed WG interface approach. In solving regular elliptic interface problems, high order convergences are numerically confirmed by using piecewise polynomial basis functions of high degrees. Moreover, the WG method is shown to be able to accommodate very complicated interfaces, due to its flexibility in choosing finite element partitions. Finally, in dealing with challenging problems with low regularities, the piecewise linear WG method is capable of delivering a second order of accuracy in L∞ norm for both C1 and H2 continuous solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗