SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Theoretical analysis of rough-surface shadowing from point-source radiation.

The shadowing function, defined as the probability that a point on a rough surface will be illuminated by radiation from a source at fixed height above the surface, is considered. An exact expression for this probability is derived and evaluated approximately when the surface elevation is a normal random process. The effect of surface correlation is considered and shown to be appreciable. Several numerical computations of the shadowing function are presented.

Hardin, J. C.↗

NASA NTRS1976

Incremental analysis of large elastic deformation of a rotating cylinder

The effect of finite deformation upon a rotating, orthotropic cylinder was investigated using a general incremental theory. The incremental equations of motion are developed using the variational principle. The governing equations are derived using the principle of virtual work for a body with initial stress. The governing equations are reduced to those for the title problem and a numerical solution is obtained using finite difference approximations. Since the problem is defined in terms of one independent space coordinate, the finite difference grid can be modified as the incremental deformation occurs without serious numerical difficulties. The nonlinear problem is solved incrementally by totaling a series of linear solutions.

Buchanan, G. R.↗

NASA NTRS1985

Approximation methods for the solution of inverse problems in lake and sea sediment core analysis

A theoretical model employing one-dimensional (depth) transport equations to describe vertical redistribution of ocean-floor and lake-floor sediment (particulates, volcanic ash, microtektites, or radioactive tracers) by episodic and nonepisodic events including bioturbation is developed analytically and demonstrated. The principles underlying the model are explained; the model equations are derived; the inverse problem of identifying the depth-dependent bioturbation coefficient is addressed; two approximation theorems are presented; and numerical results for two sample problems are presented graphically. It is suggested that compatification, porosity effects, and depth-dependent sedimentation be taken into account when formulating future models.

Banks, H. T.↗

DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Infinite quantum signal processing

Quantum signal processing (QSP) represents a real scalar polynomial of degree d using a product of unitary matrices of size 2 × 2 , parameterized by ( d + 1 ) real numbers called the phase factors. This innovative representation of polynomials has a wide range of applications in quantum computation. When the polynomial of interest is obtained by truncating an infinite polynomial series, a natural question is whether the phase factors have a well defined limit as the degree d → ∞ . While the phase factors are generally not unique, we find that there exists a consistent choice of parameterization so that the limit is well defined in the ℓ 1 space. This generalization of QSP, called the infinite quantum signal processing, can be used to represent a large class of non-polynomial functions. Our analysis reveals a surprising connection between the regularity of the target function and the decay properties of the phase factors. Our analysis also inspires a very simple and efficient algorithm to approximately compute the phase factors in the ℓ 1 space. The algorithm uses only double precision arithmetic operations, and provably converges when the ℓ 1 norm of the Chebyshev coefficients of the target function is upper bounded by a constant that is independent of d . This is also the first numerically stable algorithm for finding phase factors with provable performance guarantees in the limit d → ∞ .

Dong, Yulong [Department of Mathematics, Universit↗

DOE OSTI2020

Uncertainty quantification of an empirical shell-model interaction using principal component analysis

Recent investigations have emphasized the importance of uncertainty quantification (UQ) in nuclear theory. Here, we carry out UQ for configuration-interaction shell-model calculations in the 1$\textit{s}$–0$\textit{d}$ valence space, investigating the sensitivity of observables to perturbations in the 66 parameters (matrix elements) of a high-quality empirical interaction. The large parameter space makes computing the corresponding Hessian numerically costly, so we compare a cost-effective approximation, using the Feynman-Hellmann theorem, to the full Hessian and find it works well. Diagonalizing the Hessian yields the principal components of the interaction: linear combinations of parameters ordered by sensitivity. This approximately decoupled distribution of parameters facilitates theoretical uncertainty propagation onto structure observables: electromagnetic transitions, Gamow-Teller decays, and dark matter-nucleus scattering matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1984

An investigation of the transient thermal analysis of spur gears

A finite element computer program is developed for evaluating the transient behavior of surface temperature in high performance spur gears. The time dimension is implemented using two and three point finite difference schemes. The different schemes are provided for the purpose of numerical stability and convergence studies. A detailed explanation of the gear cooling process leading to the establishment of a modified Blok model is also included. Other conventional models for approximating the heat transfer coefficients are available for comparison. Preliminary results are given showing snap shots of gear temperature contours at the initial stages of tooth engagement.

El-Bayoumy, L. E.↗

NASA NTRS1985

An investigation of the transient thermal analysis of spur gears

El-Bayoumy, L. E.↗

NASA NTRS1983

Selection of optimum frequencies for atmospheric electric path length measurement by satellite-borne microwave radiometers

Numerical experiments using regressions by leaps and bounds have been performed to determine the optimum frequencies for satellite-borne microwave radiometers to estimate atmospheric electrical path length over the sea. The frequency range 5-40 GHz was searched. The effect of surface wind speed, sea surface temperature, and clouds was considered in the optimum frequency selection. The analysis indicates that approximately 0.6-cm rms accuracy is possible for one-way path length measurement using a proper pair of frequencies. The best two-channel subset selected by the leaps and bounds techniques is (16.0, 21.0) GHz.

Pandey, P. C.↗

NASA NTRS1988

A nonlinear analysis of an equilibrium craze. I - Problem formulation and solution. II - Simulations of craze and crack growth

A numerical method is developed for the determinations of the equilibrium shape of a craze in an infinite elastic plane whose fibrils exhibit very general nonlinear force-displacement behavior. The problem formulation is based on the superposition of the relevant elasticity Green's function; the solution of the resulting nonlinear problem is obtained by using Picard's successive approximation scheme. The model is used to investigate the effect of nonlinear fibril behavior on the mechanics of craze and crack growth, and the results are compared with the Dugdale model.

Ungsuwarungsri, T.↗

NASA NTRS2008

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗

NASA NTRS1984

Engineering analysis of drooped leading-edge wings near stall

Numerical results obtained with a three-dimensional vortex panel computer program for the calculation of inviscid, incompressible (potential) flow over infinitely thin finite wings with camber are presented. The program, which is designed for wings with drooped leading edge discontinuities, is essentially a numerical representation of lifting surface theory, involving both spanwise and chordwise distributions of vorticity. The results also include two engineering approximations. First, the effect of the leading-edge discontinuities is modeled by assuming that the vortices emanating from the discontinuities aerodynamically divide the wing into three distinct sections of lower aspect ratio. Secondly, the effect of the separated flow at high angle of attack is modeled by applying rectangular vortex panels with a varying vortex strength over the portions of the wing with attached flow.

Anderson, J. D., Jr.↗

NASA NTRS1992

Experimental results and numerical modeling of solidification during aircraft high-g arcs

The question of how the Coriolis and gravity gradient forces during high-gravity maneuvers compare to that for common centrifuges is addressed. Microstructural examination of samples solidified during high-gravity arcs reveals no evidence of convection dampening. As a first approximation, the high-gravity arc is modeled as a centrifuge with a radius of 20,480 ft and angular speed of 0.318 rpm. Scaling analysis indicates that the Coriolis and gravity gradient expected on the aircraft high gravity arc are less than that for the centrifuges by a factor of 100. Detailed Navier-Stokes analysis of the fluid flow and thermal fields during solidification of aluminum and Cd-Te during KC-135 high gravity show that convective flows of about 1 mm/s are induced. The thermal field is only slightly modified by the convection. Coriolis and gravity gradient during solidification in KC-135 high-gravity arcs, even at accelerations that have been shown to produce significant convective flow dampening in the centrifuge systems, are found to have no significant influence on the melt thermal and flowfields.

Curreri, P. A.↗

NASA NTRS1952

Analysis of temperature distribution in liquid-cooled turbine blades

The temperature distribution in liquid-cooled turbine blades determines the amount of cooling required to reduce the blade temperature to permissible values at specified locations. This report presents analytical methods for computing temperature distributions in liquid-cooled turbine blades, or in simplified shapes used to approximate sections of the blade. The individual analyses are first presented in terms of their mathematical development. By means of numerical examples, comparisons are made between simplified and more complete solutions and the effects of several variables are examined. Nondimensional charts to simplify some temperature-distribution calculations are also given.

Livingood, John N B↗

NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗

DOE OSTI2021

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1973

The parallel diffusion of cosmic rays in a random magnetic field.

Within the quasi-linear approximation, the existence of the parallel diffusion coefficient for cosmic rays in a random magnetic field (homogeneous, isotropic), despite the slow decay of the interaction between particles and random field, is demonstrated. As an example, the results of a numerical calculation of the parallel diffusion coefficient for a Gaussian random-field correlation function are presented. The numerical results are corroborated by asymptotic analysis and are compared to those of other theories.

Klimas, A.↗

NASA NTRS1982

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗