Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function 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 883 records · Page 49

On the detectability of infrared echo arcs around supernova 1987A

The ring-like interstellar visual echoes of radii 33 and 54 arcsec detected around SN 1987A should coincide with infrared echoes (thermal reradiation) from dust at T approximately equal to 15 to 30 K. Whether these infrared echoes are detectable at present is considered. They will be brightest at approximately 100 microns, the range of the Texas infrared photometer. Detectability depends on the ratio zeta congruent to tan(sub a)/tan(sub s)P(theta), where tan(sub a) and tan(sub s) are the visual absorption and scattering optical thicknesses of the echo layer, and P is the phase function function for small-angle scattering (theta approximately equal to 2 to 4 degrees). Zeta approximately greater than 1 is needed for a detectable signal (approximately 0.3 Jy), but zeta cannot be much less than 1; otherwise the visual echoes could not be as bright as they are. Typical dust mixtures of Mathis-Rumpl-Nordsieck type have zeta much less than 1. Zeta remains small even if a population of very small grains with power-law index as steep as approximately 5.5 is added. A population with even more small grains and/or fewer large grains could have a zeta similar to 1 and be detectable at present, but this seems unlikely. The echoes will move, but should remain accessible for many years and should be detected eventually.

Felten, James E.↗

Satellite Attitude Control Utilizing the Earth's Magnetic Field

A study was conducted to determine the feasibility of a satellite attitude fine-control system using the interaction of the earth's magnetic field with current-carrying coils to produce torque. The approximate intensity of the earth's magnetic field was determined as a function of the satellite coordinates. Components of the magnetic field were found to vary essentially sinusoidally at approximately twice orbital frequency. Amplitude and distortion of the sinusoidal components were a function of satellite orbit. Two systems for two-axis attitude control evolved from this study, one using three coils and the other using two coils. The torques developed by the two systems differ only when the component of magnetic field along the tracking line is zero. For this case the two-coil system develops no torque whereas the three-coil system develops some effective torque which allows partial control. The equations which describe the three-coil system are complex in comparison to those of the two-coil system and require the measurement of all three components of the magnetic field as compared with only one for the two-coil case. Intermittent three-axis torquing can also be achieved. This torquing can be used for coarse attitude control, or for dumping the stored momentum of inertia reaction wheels. Such a system has the advantage of requiring no fuel aboard the satellite. For any of these magnetic torquing schemes the power required to produce the magnetic moment and the weight of the coil seem reasonable.

White, John S.↗

Empirical Modeling of the Storm Time Innermost Magnetosphere Using Van Allen Probes and THEMIS Data: Eastward and Banana Currents

The structure of storm time currents in the inner magnetosphere, including its innermost region inside 4R(sub E), is studied for the first time using a modification of the empirical geomagnetic field model TS07D and new data from Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms missions. It is shown that the model, which uses basis-function expansions instead of ad hoc current modules to approximate the magnetic field, consistently improves its resolution and magnetic field reconstruction with the increase of the number of basis functions and resolves the spatial structure and evolution of the innermost eastward current. This includes a connection between the westward ring current flowing largely at R > or approx. 3R(sub E) and the eastward ring current concentrated at R < or approx. 3R(sub E) resulting in a vortex current pattern. A similar pattern coined 'banana current' was previously inferred from the pressure distributions based on the energetic neutral atom imaging and first-principles ring current simulations. The morphology of the equatorial currents is dependent on storm phase. During the main phase, it is complex, with several asymmetries forming banana currents. Near SYM-H minimum, the banana current is strongest, is localized in the evening-midnight sector, and is more structured compared to the main phase. It then weakens during the recovery phase resulting in the equatorial currents to become mostly azimuthally symmetric.

Stephens, G. K.↗

Incomplete fuzzy data processing systems using artificial neural network

In this paper, the implementation of a fuzzy data processing system using an artificial neural network (ANN) is discussed. The binary representation of fuzzy data is assumed, where the universe of discourse is decartelized into n equal intervals. The value of a membership function is represented by a binary number. It is proposed that incomplete fuzzy data processing be performed in two stages. The first stage performs the 'retrieval' of incomplete fuzzy data, and the second stage performs the desired operation on the retrieval data. The method of incomplete fuzzy data retrieval is proposed based on the linear approximation of missing values of the membership function. The ANN implementation of the proposed system is presented. The system was computationally verified and showed a relatively small total error.

Patyra, Marek J.↗

Two-photon excitation cross section in light and intermediate atoms in frozen-core LS-coupling approximation

Using the method of explicit summation over the intermediate states two-photon absorption cross sections in light and intermediate atoms based on the simplistic frozen-core approximation and LS coupling have been formulated. Formulas for the cross section in terms of integrals over radial wave functions are given. Two selection rules, one exact and one approximate, valid within the stated approximations are derived. The formulas are applied to two-photon absorptions in nitrogen, oxygen, and chlorine. In evaluating the radial integrals, for low-lying levels, the Hartree-Fock wave functions, and for high-lying levels, hydrogenic wave functions obtained by the quantum-defect method have been used. A relationship between the cross section and the oscillator strengths is derived.

Omidvar, K.↗

Galerkin approximations for dissipative magnetohydrodynamics

A Galerkin approximation scheme is proposed for voltage-driven, dissipative magnetohydrodynamics. The trial functions are exact eigenfunctions of the linearized continuum equations and represent helical deformations of the axisymmetric, zero-flow, driven steady state. The lowest nontrivial truncation is explored: one axisymmetric trial function and one helical trial function each for the magnetic and velocity fields. The system resembles the Lorenz approximation to Benard convection, but in the region of believed applicability, its dynamical behavior is rather different, including relaxation to a helically deformed state similar to those that have emerged in the much higher resolution computations of Dahlburg et al.

Chen, Hudong↗

Comparison of measurements by different instruments of the far-infrared reflectance of rough, optically black coatings

Two comparisons of available specular infrared measurements using identical or very similar coatings reveal that specular BRDF (bidirectional reflectance distribution function) measurements of rough coatings are comparable but that simple measured reflectances are not. This is seen as indicating that the diffuse component of the measured specular reflectance is significantly larger than that of the specular component. It is noted that the diffuse component of the measurement is a convolution of the intrinsic directional reflectivity of the surface with the instrument function and that because the value of that function increases at the wavelengths under consideration, the diffuse component is larger than would otherwise be expected. The BRDF normalization is found to cause a primitive deconvolution of the instrument function and to produce a quantity that is approximately equal to the mean intrinsic directional reflectivity very near the specular direction. The analysis and comparison of data reveal that nonspecular BRDF measurements at angles well outside the instrument function can be calibrated as well by an image of the source upon the detector as by reference to a known surface.

Smith, S. M.↗

Spectral filtering for plant production

Both plants and animals have one general commonality in their perception of light. They both are sensitive primarily to the 400 to 700 nm wavelength portion of the electromagnetic spectrum. This is referred to as the visible spectrum for animals and as the photosynthetically active radiation (PAR) spectrum for plants. Within this portion of the spectrum, animals perceive colors. Relatively recently it has been learned that within this same spectral range plants also demonstrate varying responses at different wavelengths, somewhat analogous to the definition of various colors at specific wavelengths. Although invisible to the human eye, portions of the electromagnetic spectrum on either side of the visible range are relatively inactive photosynthetically but have been found to influence important biological functions. These portions include the ultraviolet (UV approximately equal to 280-400 nm) and the far-red (FR approximately equal to 700-800 nm). The basic photoreceptor of plants for photosynthesis is chlorophyll. It serves to capture radiant energy which combined with carbon dioxide and water produces oxygen and assimulated carbon, used for the synthesis of cell wall polysaccarides, proteins, membrane lipids and other cellular constituents. The energy and carbon building blocks of photosynthesis sustain growth of plants. On the other hand, however, there are other photoreceptors, or pigments, that function as signal transducers to provide information that controls many physiological and morphological responses of how a plant grows. Known photomorphogenic receptors include phytochrome (the red/far-red sensor in the narrow bands of 655-665 nm and 725-735 nm ranges, respectively) and 'cryptochrome' (the hypothetical UV-B sensor in the 280-320 nm range). Since the USDA team of W. L. Butler, S. B. Hendricks, H. A. Borthwick, H. A. Siegleman and K. Norris in Beltsville, MD detected by spectroscopy, extracted and identified phytochrome as a protein in the 1950's, many other investigators have found evidence of its control functions in plants. Considerably less, however, is known about the yet non-isolated cryptochrome. The information-transferring roles of photoreceptors in plants at specific spectral ranges quite naturally stimulated plant scientists and engineers to consider physically manipulating light to achieve desired physiological and morphological characteristics. One way to manipulate light is to filter it through materials that selectively transmit portions of the sun's spectrum in and near the PAR range.

Young, Roy E.↗

Feature selection for best mean square approximation of class densities

A criterion for linear feature selection is proposed which is based on mean square apporximation of class density functions. It is shown that for the widest possible class of approximants, the criterion reduces to Devijver's Bayesian distance. For linear approximants the criterion is equivalent to well known generalized Fisher criteria.

Peters, C.↗

Stability of dithered non-linear systems with backlash or hysteresis

A study is conducted of the effect of dither on the nonlinear element of a single-input single-outout feedback system. Nonlinearities are considered with memory (backlash, hysteresis), in the feedforward loop; a dither of a given amplitude is injected at the input of the nonlinearity. The nonlinearity is followed by a linear element with low-pass characteristic. The stability of the dithered system and an approximate equivalent system (in which the nonlinearity is a smooth function) are compared. Conditions on the input and on the dither frequency are obtained so that the approximate-system stability guarantees that of the given hysteretic system.

Desoer, C. A.↗

Performance test of a laboratory pump for liquid transfer based on the fountain effect

A laboratory scale pump which determines the flow characteristics of a heater-activated all-metal thermomechanical fountain-effect pump (FEP) is evaluated. Results are presented on the functional dependence of the fountain pressure difference versus mass throughput. A modified Vote et al. (1971) power law approximation is used to obtain the flow rate as a function of the driving force. Despite a large nominal pore size of the porous plug used for the FEP, flow rates of up to 10 liter/hr sq cm have been obtained.

Chen, W. E. W.↗

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

VEXT: A Virtual Observatory Exploration Toolkit

This final report consists of two main parts. The first is taken from a paper by the PiCA (Pittsburgh Computational Astrostatistics) Group which describes our ongoing work in fast computation of n-point correlation functions. We present here a new algorithm for the fast computation of N-point correlation functions in large astronomical data sets. The algorithm is based on kd-trees which are decorated with cached sufficient statistics thus allowing for orders of magnitude speed-ups over the naive non-tree-based implementation of correlation functions. We further discuss the use of controlled approximations within the computation which allows for further acceleration. In summary, our algorithm now makes it possible to compute exact, all-pairs, measurements of the two, three and four-point correlation functions for cosmological data sets like the Sloan Digital Sky Survey and the next generation of Cosmic Microwave Background experiments. The second part summarizes the progress made by the PiCA Group in this area through the AISR grant.

Schneider, Jeff↗

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

Corbetta, M.↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

stochastic calculus↗