Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for both the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.

nite element methods↗
NASA NTRS2007

Optimal Low-Thrust Limited-Power Transfers between Arbitrary Elliptic Coplanar Orbits

In this work, a complete first order analytical solution, which includes the short periodic terms, for the problem of optimal low-thrust limited-power transfers between arbitrary elliptic coplanar orbits in a Newtonian central gravity field is obtained through Hamilton-Jacobi theory and a perturbation method based on Lie series.

daSilvaFernandes, Sandro↗
NASA NTRS1984

Development of a dynamically adaptive grid method for multidimensional problems

An approach to solution adaptive grid generation for use with finite difference techniques, previously demonstrated on model problems in one space dimension, has been extended to multidimensional problems. The method is based on the popular elliptic steady grid generators, but is 'dynamically' adaptive in the sense that a grid is maintained at all times satisfying the steady grid law driven by a solution-dependent source term. Testing has been carried out on Burgers' equation in one and two space dimensions. Results appear encouraging both for inviscid wave propagation cases and viscous boundary layer cases, suggesting that application to practical flow problems is now possible. In the course of the work, obstacles relating to grid correction, smoothing of the solution, and elliptic equation solvers have been largely overcome. Concern remains, however, about grid skewness, boundary layer resolution and the need for implicit integration methods. Also, the method in 3-D is expected to be very demanding of computer resources.

Holcomb, J. E.↗
DOE OSTI2015

Fully-Implicit Orthogonal Reconstructed Discontinuous Petrov-Galerkin Method for Multiphysics Problems

This document is focussed on the development of reconstructed discontinuous Galerkin (RDG P n P m ) method, based on orthogonal basis functions. Orthogonality of basis functions is essential for enabling robust and efficient fully-implicit Newton-Krylov based time discretization. The method is designed for generic partial differential equations, including transient hyperbolic, parabolic and/or elliptic operators, which are attributed to many multiphysics problems. Our driving application of interest is additive manufacturing, which requires simulations of fluid-solid-gas systems, with compressibility effects and phase change (melting/ solidification). Thus, we demonstrate the method capabilities for solving nonlinear heat conduction equation coupled with multi-material Navier-Stokes based compressible fluid (in low Mach number limit) systems. We focus our attention to the method accuracy (in both time and space), as well as solvability of linear algebra involved in linear steps of the Newton-based solver. We demonstrate the accuracy and efficacy of the method, emphasizing the advantages from the orthogonality of basis functions, which makes better conditioning of underlying (approximate) Jacobian matrices, and rapid convergence of Krylov methods. iii

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Likelihood-informed dimension reduction for nonlinear inverse problems

The intrinsic dimensionality of an inverse problem is affected by prior information, the accuracy and number of observations, and the smoothing properties of the forward operator. From a Bayesian perspective, changes from the prior to the posterior may, in many problems, be confined to a relatively low-dimensional subspace of the parameter space. We present a dimension reduction approach that defines and identifies such a subspace, called the 'likelihood-informed subspace' (LIS), by characterizing the relative influences of the prior and the likelihood over the support of the posterior distribution. This identification enables new and more efficient computational methods for Bayesian inference with nonlinear forward models and Gaussian priors. In particular, we approximate the posterior distribution as the product of a lower-dimensional posterior defined on the LIS and the prior distribution marginalized onto the complementary subspace. Markov chain Monte Carlo sampling can then proceed in lower dimensions, with significant gains in computational efficiency. We also introduce a Rao–Blackwellization strategy that de-randomizes Monte Carlo estimates of posterior expectations for additional variance reduction. We demonstrate the efficiency of our methods using two numerical examples: inference of permeability in a groundwater system governed by an elliptic PDE, and an atmospheric remote sensing problem based on Global Ozone Monitoring System (GOMOS) observations.

42 ENGINEERING↗
DOE OSTI2017

Quantum geometry of resurgent perturbative/nonperturbative relations

For a wide variety of quantum potentials, including the textbook ‘instanton’ examples of the periodic cosine and symmetric double-well potentials, the perturbative data coming from fluctuations about the vacuum saddle encodes all non-perturbative data in all higher non-perturbative sectors. Here we unify these examples in geometric terms, arguing that the all-orders quantum action determines the all-orders quantum dual action for quantum spectral problems associated with a classical genus one elliptic curve. Furthermore, for a special class of genus one potentials this relation is particularly simple: this class includes the cubic oscillator, symmetric double-well, symmetric degenerate triple-well, and periodic cosine potential. These are related to the Chebyshev potentials, which are in turn related to certain N = 2 supersymmetric quantum field theories, to mirror maps for hypersurfaces in projective spaces, and also to topological c = 3 Landau-Ginzburg models and ‘special geometry’. These systems inherit a natural modular structure corresponding to Ramanujan’s theory of elliptic functions in alternative bases, which is especially important for the quantization. Insights from supersymmetric quantum field theory suggest similar structures for more complicated potentials, corresponding to higher genus. Lastly, our approach is very elementary, using basic classical geometry combined with all-orders WKB.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1966

Elliptic motion

Newtonian theory, Keplerian law, and two-body problem applied to orbital motion

NEWTON THEORY↗
NASA NTRS1981

Transonic airfoil calculations using solution-adaptive grids

A new algorithm for generating solution-adaptive grids (SAG) about airfoil configurations embedded in transonic flow is presented. The present SAG approach uses only the airfoil surface solution to recluster grid points on the airfoil surface, i.e., the reclustering problem is one dimension smaller than the flow-field calculation problem. Special controls automatically built into the elliptic grid generation procedure are then used to obtain grids with suitable interior behavior. This concept of redistributing grid points greatly simplifies the idea of solution-adaptive grids. Numerical results indicate significant improvements in accuracy for SAG grids relative to standard grids using the same number of points.

Holst, T. L.↗
NASA NTRS1991

Compressible flows with periodic vortical disturbances around lifting airfoils

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method that is presented fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. The numerical scheme is validated through extensive comparisons with known solutions to unsteady vortical flow problems. In general, it is seen that the agreement between the numerical and analytical results is very good for reduced frequencies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. Numerical results are also presented for a wide variety of flow configurations for the purpose of determining the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. It is seen that each of these parameters can have a significant effect on the unsteady airfoil response to the incident disturbances, and that the effect depends strongly upon the reduced frequency and the dimensionality of the gust. For a one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, the results indicate that airfoil thickness increases the unsteady lift and moment at the low reduced frequencies but decreases it at the high reduced frequencies. The results show that an increase in airfoil Mach number leads to a significant increase in the unsteady lift and moment for the low reduced frequencies, but a significant decrease for the high reduced frequencies.

Scott, James R.↗
NASA NTRS1991

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗
NASA NTRS1990

On the three-dimensional instability of strained vortices

The three-dimensional (3-D) instability of a two-dimensional (2-D) flow with elliptical streamlines has been proposed as a generic mechanism for the breakdown of many 2-D flows. A physical interpretation for the mechanism is presented together with an analytical treatment of the problem. It is shown that the stability of an elliptical flow is governed by an Ince equation. An analytical representation for a localized solution is given and establishes a direct link with previous computations and experiments.

Waleffe, Fabian↗
NASA NTRS1987

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.

Shankar, Vijaya↗
NASA NTRS2005

How Does Abundance Affect the Strength of UV Emission in Elliptical Galaxies?

This program used the Far Ultraviolet Spectroscopic Explorer (FUSE) to observe elliptical galaxies with the intention of measuring the chemical abundances in their hot stellar populations. It was designed to complement an earlier FUSE program that observed elliptical galaxies with strong UV emission. The current program originally planned observations of two ellipticals with weak UV emission (M32 and M49). Once FUSE encountered pointing control problems in certain regions of the sky (particularly Virgo, which is very unfortunate for the study of ellipticals in general), M49 was replaced with the bulge of M31, which has a similar UV-to-optical flux ratio as the center of M49. As the closest elliptical galaxy and the one with the weakest UV-to-optical flux ratio, M32 was an obvious choice of target, but M49 was the ideal complementary target, because it has a very low reddening (unlike M32). With the inability of FUSE to point at Virgo, nearly all of the best elliptical galaxies (bright galaxies with low foreground extinction) were also lost, and this severely hampered three FUSE programs of the PI, all focused on the hot stellar populations of ellipticals. M31 was the best replacement for M49, but like M32, it suffers from significant foreground reddening. Strong Galactic ISM lines heavily contaminate the FUSE spectra of M31 and M32. These ISM lines are coincident with the photospheric lines from the stellar populations (whereas M49, with little foreground ISM and significant redshift, would not have suffered from this problem). We have reduced the faint (and thus difficult) data for M31 and M32, producing final co-added spectra representing all of the exposures, but we have not yet finished our analysis, due to the complication of the contaminating ISM. The silver lining here is the set of CHI lines at 1175 Angstroms, which are not significantly contaminated by the ISM. A comparison of the M31 spectrum with other galaxies observed by FEE showed a surprising result: the hot stars in M31 seem to have a similar carbon abundance to those stars in galaxies with much brighter UV emission. The fraction of these hot stars in a population should be a strong function of chemical abundances, so this finding warrants further exploration, and we are proceeding with our analysis. Because the UV emission in these galaxies comes from a population of extreme horizontal branch stars, the PI (Brown) presented this result at a June 2003 conference on such stars.

Sonneborn, George↗
NASA NTRS1986

Survey of the status of finite element methods for partial differential equations

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Temam, Roger↗
DOE OSTI2018

Contact Problem of Wear of the Elastic Half Plane with Winkler’s Coating Caused by Punches of Canonical Shapes

The problem of contact interaction of punches of canonical shapes (cylindrical, elliptic, and hyperbolic) with elastic half planes covered with protective coatings is solved with regard for the wear of the material. We propose a method for the solution of the singular integrodifferential equation of the problem. The numerical analysis of contact pressure is performed for different times, and some specific features of wear of the material by the punches of different shapes are discovered. The time for which the coating is totally worn out at a given point of the contact zone is determined.

42 ENGINEERING↗