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

Coaxial Cables for Martian Extreme Temperature Environments

Work was conducted to validate the use of the rover external flexible coaxial cabling for space under the extreme environments to be encountered during the Mars Science Laboratory (MSL) mission. The antennas must survive all ground operations plus the nominal 670-Martian-day mission that includes summer and winter seasons of the Mars environment. Successful development of processes established coaxial cable hardware fatigue limits, which were well beyond the expected in-flight exposures. In keeping with traditional qualification philosophy, this was accomplished by subjecting flight-representative coaxial cables to temperature cycling of the same depth as expected in-flight, but for three times the expected number of in-flight thermal cycles. Insertion loss and return loss tests were performed on the coaxial cables during the thermal chamber breaks. A vector network analyzer was calibrated and operated over the operational frequency range 7.145 to 8.450 GHz. Even though some of the exposed cables function only at UHF frequencies (approximately 400 MHz), the testing was more sensitive, and extending the test range down to 400 MHz would have cost frequency resolution. The Gore flexible coaxial cables, which were the subject of these tests, proved to be robust and displayed no sign of degradation due to the 3X exposure to the punishing Mars surface operations cycles.

Ramesham, Rajeshuni↗

Large Amplitude Whistlers in the Magnetosphere Observed with Wind-Waves

We describe the results of a statistical survey of Wind-Waves data motivated by the recent STEREO/Waves discovery of large-amplitude whistlers in the inner magnetosphere. Although Wind was primarily intended to monitor the solar wind, the spacecraft spent 47 h inside 5 R(sub E) and 431 h inside 10 R(sub E) during the 8 years (1994-2002) that it orbited the Earth. Five episodes were found when whistlers had amplitudes comparable to those of Cattell et al. (2008), i.e., electric fields of 100 m V/m or greater. The whistlers usually occurred near the plasmapause. The observations are generally consistent with the whistlers observed by STEREO. In contrast with STEREO, Wind-Waves had a search coil, so magnetic measurements are available, enabling determination of the wave vector without a model. Eleven whistler events with useable magnetic measurements were found. The wave vectors of these are distributed around the magnetic field direction with angles from 4 to 48deg. Approximations to observed electron distribution functions show a Kennel-Petschek instability which, however, does not seem to produce the observed whistlers. One Wind episode was sampled at 120,000 samples/s, and these events showed a signature that is interpreted as trapping of electrons in the electrostatic potential of an oblique whistler. Similar waveforms are found in the STEREO data. In addition to the whistler waves, large amplitude, short duration solitary waves (up to 100 mV/m), presumed to be electron holes, occur in these passes, primarily on plasma sheet field lines mapping to the auroral zone.

Kellogg, P. J.↗

The Net Decay Time of Anomalies in Concentrations of Atmospheric Pollutants

This paper introduces a new parameter to characterize the random component in temporal variability of atmospheric pollutants and proposes a simple statistical technique for its evaluation. That parameter is the net decay time (or the time scale) of the local anomalies in concentrations of atmospheric pollutants, rather than the traditionally used chemical lifetimes of total amounts of the species. Using widely available data of hourly multi-year surface trace gas pollutant concentrations we demonstrate a simplified way to estimate the net decay time with an exponential approximation of lag-correlation functions. We assessed the decay times of fluctuations in observations of eight atmospheric pollutants (SO2, NO, NO2, NOy, O3, CO, NH3, and HNO3) at two urban sites and one cleaner rural site in the Eastern US. The time scales of temporal fluctuations (net decay times) vary from about one hour to slightly more than one day. These scales are generally much shorter in urban environments than in remote regions. We also compared day- and night-time observations in warm and cold seasons. At night in the cold season, time scales of fluctuations in atmospheric pollutants are usually the longest. Such estimates should be useful to air quality prediction, public health, and satellite remote sensing research communities.

air quality prediction↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Downlink Probability Density Functions for EOS-McMurdo Sound

The visibility times and communication link dynamics for the Earth Observations Satellite (EOS)-McMurdo Sound direct downlinks have been studied. The 16 day EOS periodicity may be shown with the Goddard Trajectory Determination System (GTDS) and the entire 16 day period should be simulated for representative link statistics. We desire many attributes of the downlink, however, and a faster orbital determination method is desirable. We use the method of osculating elements for speed and accuracy in simulating the EOS orbit. The accuracy of the method of osculating elements is demonstrated by closely reproducing the observed 16 day Landsat periodicity. An autocorrelation function method is used to show the correlation spike at 16 days. The entire 16 day record of passes over McMurdo Sound is then used to generate statistics for innage time, outage time, elevation angle, antenna angle rates, and propagation loss. The levation angle probability density function is compared with 1967 analytic approximation which has been used for medium to high altitude satellites. One practical result of this comparison is seen to be the rare occurrence of zenith passes. The new result is functionally different than the earlier result, with a heavy emphasis on low elevation angles. EOS is one of a large class of sun synchronous satellites which may be downlinked to McMurdo Sound. We examine delay statistics for an entire group of sun synchronous satellites ranging from 400 km to 1000 km altitude. Outage probability density function results are presented three dimensionally.

Christopher, P.↗

Generalized Sternheimer potential.

Generalization of Sternheimer potential to include wave functions involving spin in nonfactorable way, exhibiting approximate Hamiltonian

QUANTUM THEORY↗

The delta-four-stream approximation for radiative flux transfer

An approximate technique for solving the radiative transfer problem for fluxes reflected, transmitted, and absorbed by a homogeneous scattering layer is presented. The technique is a straightforward delta-function modification of a solution to the transfer equation in closed form via the discrete-ordinates method with four streams. By the use of four streams, significant improvement in accuracy is obtained for anisotropic phase functions over that obtained by similarly employed two-stream approximations with delta-functions. However, the computational effort is increased only slightly. Equations are developed for the general case of a layer underlaid by a reflecting Lambert surface. Graphical comparisons are given of fractional error resulting from the use of this method with that resulting from the use of typical four-stream and delta-two-stream techniques.

Cuzzi, J. N.↗

A Test of Thick-Target Nonuniform Ionization as an Explanation for Breaks in Solar Flare Hard X-Ray Spectra

Solar nonthermal hard X-ray (HXR) flare spectra often cannot be fitted by a single power law, but rather require a downward break in the photon spectrum. A possible explanation for this spectral break is nonuniform ionization in the emission region. We have developed a computer code to calculate the photon spectrum from electrons with a power-law distribution injected into a thick-target in which the ionization decreases linearly from 100% to zero. We use the bremsstrahlung cross-section from Haug (1997), which closely approximates the full relativistic Bethe-Heitler cross-section, and compare photon spectra computed from this model with those obtained by Kontar, Brown and McArthur (2002), who used a step-function ionization model and the Kramers approximation to the cross-section. We find that for HXR spectra from a target with nonuniform ionization, the difference (Delta-gamma) between the power-law indexes above and below the break has an upper limit between approx.0.2 and 0.7 that depends on the power-law index delta of the injected electron distribution. A broken power-law spectrum with a. higher value of Delta-gamma cannot result from nonuniform ionization alone. The model is applied to spectra obtained around the peak times of 20 flares observed by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI from 2002 to 2004 to determine whether thick-target nonuniform ionization can explain the measured spectral breaks. A Monte Carlo method is used to determine the uncertainties of the best-fit parameters, especially on Delta-gamma. We find that 15 of the 20 flare spectra require a downward spectral break and that at least 6 of these could not be explained by nonuniform ionization alone because they had values of Delta-gamma with less than a 2.5% probability of being consistent with the computed upper limits from the model. The remaining 9 flare spectra, based on this criterion, are consistent with the nonuniform ionization model.

Holman, gordon↗

Modeling the Variability of the BL Lacertae Object PKS 2155-304

The bright X-ray-selected BL Lacertae object PKS 2155-304 has been the target of two intense multi-wavelength campaigns in 1991 November and in 1994 May. Although the spectral energy distributions at both epochs were quite similar, the source exhibited two very distinct variability patterns that cannot be easily reconciled with homogeneous one-zone C7 jet models. During the first epoch the variability was almost achromatic in amplitude, with a time lag between X-rays and ultraviolet radiation (UV) of approximately three hours, while during the second epoch the variability amplitude increased as a function of wavelength with the EUV flare peaking approximately one day after the X-ray flare. We model the source using a time-dependent inhomogeneous accelerating jet model. We reproduce the general characteristics of the different variability signatures by assuming that plasma disturbances with different physical properties propagate downstream in an underlying jet characterized by the same set of physical parameters at both epochs. A time delay of approximately one day between the hardening of the UV spectral index and the UV flux present at both epochs is modeled with stochastic fluctuations in the particle acceleration manifested through small variations of the maximum energy of the injected electrons. We predict that similar time delays will be present in future observations even in the absence of strong variability event. We stress the importance of observations at neighboring frequencies as a diagnostic tool for the structure of the quiescent jet in blazars especially in the seemingly dull case in which strong variability is absent.

Georganopuolos, Markos↗

Reflectance Modeling

The main thrust of the work conducted is in the modeling of spectral responses of discontinuous conifer forest canopies, considered as collections of conical forms that are illuminated at an angle and cast shadows on contracting background. The three first-year objectives are: (1) modification of the variance-driven model to incorporate overlapping of shadows and crowns; (2) collection and analysis of photographic data to calibrate tree-spacing functions; and (3) collection of field data to determine height, spacing, and shape of conifers in open and closed stands in the Goosenest (California) test area. Modification of the invertible model uses a linear approximation to correct for the reduction in shadow and canopy area that occurs when shadows fall on crowns or tree crowns intersect. Comparison of the results of the approximation with a series of Monte Carlo simulations shows that the mean area covered by crowns or shadows is approximated quite accurately, but the variance calculated departs significantly from simulated values. The analysis of tree spacing patterns using air photos showed that spacing tended to be somewhat regular at a scale approximately equal to the average distance between trees, but at larger distance scales counts of trees in cells could be approximated acceptably by a Poisson (random) function.

Strahler, A. H.↗

A frequency domain analysis for damped space structures

A method is presented for the analysis of damped structural systems in which the structural components are represented by impedance models and analyzed in the frequency domain. Methods are presented to assemble and condense system impedance matrices, and then to identify approximate mass, stiffness, and damping matrices for systems whose impedances are complicated functions of frequency. Formulas are derived for determination of approximate values for system natural frequencies and damping using frequency domain quantities. The sensitivities of these approximate values to system parameter changes are analyzed. The implementation of these analysis tools is discussed and applied to a simple mechanical system.

Hagood, Nesbitt W.↗

A theoretical derivation of the dependence of the remotely sensed reflectance of the ocean on the inherent optical properties

An expression for the ratio of the upwelling nadir radiance L(pi, z) and the downwelling scalar irradiance E(sub od)(Z) is derived from the following equation of radiative transfer. This expression is given by RSR(z) = (L(pi, z))/E(sub od) = (f(sub b)(z)b(sub b)(z))/2 pi(k(pi, z) + c(z) - F(sub L)(z)b(sub f)(z)), where b(sub b)(z) is the backscattering coefficient, k(pi, z) is the vertical attenuation coefficient of the nadir radiance, c(z) is the beam attenuation coefficient, and f(sub b)(z) and f(sub L)(z) are shape parameters that depend on the shape of the volume scattering function and the radiance distribution. Successive approximations are subsequently applied to the above exact equation. These are f(sub b)(z) = (2 pi beta(pi - theta(sub m), z)/(b(sub b)(z))), where beta(pi - theta(sub m), z) is the volume scattering function at 180 deg minus the zenith angle of the maximum radiance, and k(pi, z) = am = c(1 - 0.52 b/c - 0.44 (b/c)(exp 2)), where m is a parameter that is numerically equal to the inverse of the average cosine of the asymptotic light field for a medium with the same inherent optical properties, a is the absorption coefficient, and b/c is the single scattering albedo. Together with f(sub L)(z) = 1.05 and application of Gershun's equation, it is shown that for nearly all oceanic cases RSR(z) identical to L(pi, z)/E(sub od)(z) = (Beta(pi - theta(sub m), z))/(a(z)(1 + m(z))).

Zaneveld, J. Ronald V.↗

Extension of the quasistatic far-wing line shape theory to multicomponent anisotropic potentials

The formalism developed previously for the calculation of the far-wing line shape function and the corresponding absorption coefficient using a single-component anisotropic interaction term and the binary collision and quasistatic approximations is generalized to multicomponent anisotropic potential functions. Explicit expressions are presented for several common cases, including the long-range dipole-dipole plus dipole-quadrupole interaction and a linear molecule interacting with a perturber atom. After determining the multicomponent functional representation for the interaction between the CO2 and Ar from previously published data, we calculate the theoretical line shape function and the corresponding absorption due to the nu(sub 3) band of CO2 in the frequency range 2400-2580 cm(exp -1) and compare our results with previous calculations carried out using a single-component anisotropic interaction, and with the results obtained assuming Lorentzian line shapes. The principal uncertainties in the present results, possible refinements of the theoretical formalism, and the applicability to other systems are discussed briefly.

Ma, Q.↗

Modal density function and number of propagating modes in ducts

Often raised questions in duct sound propagation studies involve the total number of propagating modes, the number of propagating radial modes for a particular spinning lobe number, and the number of modes possible between two given values of cutoff ratio or eigenvalue. These questions can be answered approximately by using the modal distribution function which is the integral of the modal density function for ducts in a manner similar to that previously published for architectural acoustics. The modal density functions are derived for rectangular and circular ducts with a uniform steady flow. Results from this continuous theory are compared to the actual (discrete) modal distributions.

Rice, E. J.↗

Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids

This paper is concerned with two important elements in the high-order accurate spatial discretization of finite volume equations over arbitrary grids. One element is the integration of basis functions over arbitrary domains, which is used in expressing various spatial integrals in terms of discrete unknowns. The other consists of quadrature approximations to those integrals. Only polynomial basis functions applied to polyhedral and polygonal grids are treated here. Non-triangular polygonal faces are subdivided into a union of planar triangular facets, and the resulting triangulated polyhedron is subdivided into a union of tetrahedra. The straight line segment, triangle, and tetrahedron are thus the fundamental shapes that are the building blocks for all integrations and quadrature approximations. Integrals of products up to the fifth order are derived in a unified manner for the three fundamental shapes in terms of the position vectors of vertices. Results are given both in terms of tensor products and products of Cartesian coordinates. The exact polynomial integrals are used to obtain symmetric quadrature approximations of any degree of precision up to five for arbitrary integrals over the three fundamental domains. Using a coordinate-free formulation, simple and rational procedures are developed to derive virtually all quadrature formulas, including some previously unpublished. Four symmetry groups of quadrature points are introduced to derive Gauss formulas, while their limiting forms are used to derive Lobatto formulas. Representative Gauss and Lobatto formulas are tabulated. The relative efficiency of their application to polyhedral and polygonal grids is detailed. The extension to higher degrees of precision is discussed.

Liu, Yen↗

Optimal filtering in the presence of unmodeled time correlated driving disturbances.

In many realistic data filtering problems, the cross correlation of the state estimation error and the state forcing function is unknown due to the poor knowledge of the time history of the forcing function. In this paper, the conservative and minimal approximation to the cross correlation terms is presented. It requires only the knowledge of the estimation error covariance and the forcing function covariance, with the choice of an associated free parameter left to the user. If the estimation error covariance and/or the forcing function covariance are bounded from above but not known exactly, the cross correlation approximation using those upper bounds remains conservative. This cross correlation approximation leads to a conservative approximation to the estimation error covariance matrix differential equation between measurement times. The free parameter is determined as the analytic solution to an associated optimal control problem. The procedure is expanded to include discrete linear measurement incorporation.

Fraser, D. C.↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗