Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional approximation”

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 577 records · Page 32

High temperature ablative foam

An ablative foam composition is formed of approximately 150 to 250 parts by weight polymeric isocyanate having an isocyanate functionality of 2.6 to 3.2; approximately 15 to 30 parts by weight reactive flame retardant having a hydroxyl number range from 200-260; approximately 10 to 40 parts by weight non-reactive flame retardant; approximately 10 to 40 parts by weight nonhydrolyzable silicone copolymer having a hydroxyl number range from 75-205; and approximately 3 to 16 parts by weight amine initiated polyether resin having an isocyanate functionality greater than or equal to 3.0 and a hydroxyl number range from 400-800.

Liu, Matthew T.↗

Simplified computational approach for dual-probe heat-pulse method

Two equations are currently available for estimating soil volumetric heat capacity (pc) with the dual-probe heat-pulse (DPHP) method. One is simple but gives only approximate results because it assumes that the DPHP sensor releases an impulse of heat instantaneously. The other explicitly accounts for the finite duration of heating and gives exact results. Unfortunately, the equation that gives exact results involves the exponential integral function, which is not available in most computer spreadsheet software packages or data logger function libraries. In this note we introduce an approximation of the exact equation that contains only simple algebraic functions. The approximation consists of the first five terms of a Taylor series, which are written as a telescoped polynomial for computational purposes. For most applications of the DPHP method, the polynomial approximation gives estimates of pc that are at least an order of magnitude more accurate than estimates obtained from the simple equation based on instantaneous heating.

NASA Center JSC↗

Control law synthesis and sensor design for active flutter suppression.

Methods are presented for representing unsteady aerodynamic loadings, valid for arbitrary motion, for the two-dimensional typical section with a trailing-edge control surface in incompressible flow and a three-dimensional lifting surface with leading- or trailing-edge control surfaces in subsonic compressible flow. Loads for the two-dimensional incompressible case are obtained via analytic continuation of the Theodorsen function into the complex plane, with guidance from a time-domain approximation of the Wagner indicial function. For the three-dimensional case, it is proposed that oscillatory generalized aerodynamic forces be approximated by Pade fractions, thereby permitting an obvious continuation into the complex plane. A theoretical justification for this procedure is briefly outlined, and some practical aspects of its implementation are discussed.

Lyons, M. G.↗

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗

Direct Calculation of the Scattering Amplitude Without Partial Wave Decomposition: Inclusion of Correlation Effects - III

In the first two papers in this series, we developed a method for studying electron-hydrogen scattering that does not use partial wave analysis. We constructed an ansatz for the wave function in both the static and static exchange approximations and calculated the full scattering amplitude. Here we go beyond the static exchange approximation, and include correlation in the wave function via a modified polarized orbital. This correlation function provides a significant improvement over the static exchange approximation: the resultant elastic scattering amplitudes are in very good agreement with fully converged partial wave calculations for electron-hydrogen scattering. A fully variational modification of this approach is discussed in the conclusion of the article Popular summary of Direct calculation of the scattering amplitude without partial wave expansion. III ....." by J. Shertzer and A. Temkin. In this paper we continue the development of In this paper we continue the development of a new approach to the way in which researchers have traditionally used to calculate the scattering cross section of (low-energy) electrons from atoms. The basic mathematical problem is to solve the Schroedinger Equation (SE) corresponding the above physical process. Traditionally it was always the case that the SE was reduced to a sequence of one-dimensional (ordinary) differential equations - called partial waves which were solved and from the solutions "phase shifts" were extracted, from which the scattering cross section was calculated.

Shertzer, Janine↗

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Application

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Applications

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

MAFT: The Multicomputer Architecture for Fault-Tolerance

Multicomputer Architecture for Fault-Tolerance (MAFT) is a loosely coupled multiprocessor system designed to achieve an unreliability of less than 10(exp -10)/hr in flight-critical real time applications. The MAFT design objectives and architecture are presented. The fault-tolerance implementation of major functions in MAFT is also presented, including communication; task scheduling; reconfiguration; clock synchronization; and data handling and voting. The need for Byzantine agreement or approximate agreement in various functions is discussed. Different methods were selected to achieve agreement in various subsystems. These methods are illustrated by a more detailed description of the task scheduling and error handling subsystems.

Kieckhafer, Roger M.↗

Practical Probabilistic Programming

Recent advances in probabilistic programming languages (PPLs) have provided the capability for exact inference: computing a closed-form probability distribution for a given probabilistic program. In particular, the new language Roulette uses a language oriented programming (LOP) approach, wherein analysts build new programming languages on top of a set of primitives provided by Roulette, which then translates these structures into a weighted model counting problem which can be solved by automated reasoning tools. However, because Roulette provides few convenience features, developing these new languages is challenging even for expert users. We developed a standard library of common probability functions for Roulette with the goal of improved usability. This included approximation of continuous probability density functions using discrete probability mass functions. We demonstrated this approach by modeling a cosmic ray striking a RAM controller. We found that Roulette provides a powerful interface for highly expressive probabilistic programs to be generated. In collaboration with the NNSA Advanced Simulation and Computing program, which resulted in development of a tool called Circulette, we were able to model complex circuits expressed in Verilog using probabilistic programs with an expressivity not previously possible. Our research question that motivated the development of a Roulette standard library was to determine whether non-experts could use a PPL to model relevant problems regarding radiation effects on microelectronics. This standard library improved the expressivity of Roulette by implementing common probability density functions, mathematical operators on distributions, and support for empirical distributions. While Roulette is a powerful modeling language, the untyped, LOP approach makes error messages difficult to understand and requires expert aid. We recommend further research on Roulette, especially with its error messages, to enable improved usability. At the same time, this project demonstrated that for users familiar with Roulette and the LOP approach, Roulette provides powerful new capabilities that can be integrated with other Sandia modeling capabilities.

97 MATHEMATICS AND COMPUTING↗

Comparison of the Radiative Two-Flux and Diffusion Approximations

Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.

Spuckler, Charles M.↗

Thermodynamics and phase separation of dense fully-ionized hydrogen-helium fluid mixtures

The free energy of a hydrogen-helium fluid mixture is evaluated for the temperatures and densities appropriate to the deep interior of a giant planet such as Jupiter. The electrons are assumed to be fully pressure-ionized and degenerate. In this regime, an appropriate first approximation to the ionic distribution functions can be found by assuming hard sphere interactions. Corrections to this approximation are incorporated by means of the perturbation theory of Anderson and Chandler. Approximations for the three-body interactions and the nonlinear response of the electron gas to the ions are included. It is predicted that a hydrogen-helium mixture, containing 10% by number of helium ions, separates into hydrogen-rich and helium-rich phases below about 8000 K, at the pressures relevant to Jupiter (4-40 Megabars). It is also predicted that the alloy occupies less volume per ion than the separated phases. The equations of state and other thermodynamic derivatives are tabulated. Implications of these results are discussed.

Stevenson, D. J.↗

Quadratic Liapunov functions

Criteria for estimating errors in quadratic approximations to asymptotic stability involving Popov condition based on existence of quadratic Liapunov functions

Weissenberger, S.↗

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Variance of phase fluctuations of waves propagating through a random medium

As an electromagnetic wave propagates through a random scattering medium, such as a forest, its energy is attenuated and random phase fluctuations are induced. The magnitude of the random phase fluctuations induced is important in estimating how well a Synthetic Aperture Radar (SAR) can image objects within the scattering medium. The two-layer random medium model, consisting of a scattering layer between free space and ground, is used to calculate the variance of the phase fluctuations induced between a transmitter located above the random medium and a receiver located below the random medium. The scattering properties of the random medium are characterized by a correlation function of the random permittivity fluctuations. The effective permittivity of the random medium is first calculated using the strong fluctuation theory, which accounts for large permittivity fluctuations of the scatterers. The distorted Born approximation is used to calculate the first-order scattered field. A perturbation series for the phase of the received field in the Rytov approximation is then introduced and the variance of the phase fluctuations is also calculated assuming that the transmitter and receiver are in the paraxial limit of the random medium, which allows an analytic solution to be obtained. Results are compared using the paraxial approximation, scalar Green's function formulation, and dyadic Green's function formulation. The effects studied are the dependence of the variance of the phase fluctuations on receiver location in lossy and lossless regions, medium thickness, correlation length and fractional volume of scatterers, depolarization of the incident wave, ground layer permittivity, angle of incidence, and polarization.

Chu, Nelson C.↗