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 1,153 records · Page 64

File Specification for GEOS-5 FP-IT (Forward Processing for Instrument Teams)

The GEOS-5 FP-IT Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP-IT incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP-IT uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The analysis is performed at a horizontal resolution of 0.625-degree longitude by 0.5-degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational GEOS-5 FP-IT processing in support of the EOS instrument teams. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary.

3D Var↗

Dust Abrasion Damage on Martian Solar Arrays: Experimental Investigation and Opportunity Rover Performance Analysis

Here we investigate the effects of erosion and weathering that occur on III-V cover-glass interconnected cells (CICs) after exposure to Mars dust storm conditions. The durability of these materials in a Martian environment is not well characterized so we perform analogous experimentation. To replicate the dust impingement, test coupons were placed in an enclosure and sandblasted with Mars dust simulant. We show the J-V response dependency on both incident angle and exposure times. We find that the simulated Martian dust storm often results in damage to the anti-reflective coating and subsequent reduced short circuit current. Reduction in the open circuit voltage is observed, likely caused by structural damage to the crystal lattice after CIC fracture. We employ data-driven modeling to determine a performance degradation rate that is consistent with zero within uncertainty. We also quantify the soiling contribution and power degradation of the photovoltaic cells on Mars through analysis of 4.95 Martian years of report-out power conditions from the Opportunity rover. We find that atmospheric dust suspended due to a weather event does not result in instantaneous settled dust on the PV cells. We calculate via autocorrelation function that the dust settling rate is approximately 21 Sols from atmospheric dust suspension. The findings presented here deliver a realistic approximation for the insolation values and subsequent PV power expected over time on the Martian surface thus informing future dust abatement systems.

McMillon-Brown, Lyndsey↗

Dust Abrasion Damage on Martian Solar Arrays: Experimental Investigation and Opportunity Rover Performance Analysis

Here we investigate the effects of erosion and weathering that occur on III-V cover-glass interconnected cells (CICs) after exposure to Mars dust storm conditions. The durability of these materials in a Martian environment is not well characterized so we perform analogous experimentation. To replicate the dust impingement, test coupons were placed in an enclosure and sandblasted with Mars dust simulant. We show the J-V response dependency on both incident angle and exposure times. We find that the simulated Martian dust storm often results in damage to the anti-reflective coating and subsequent reduced short circuit current. Reduction in the open circuit voltage is observed, likely caused by structural damage to the crystal lattice after CIC fracture. We employ data-driven modeling to determine a performance degradation rate that is consistent with zero within uncertainty. We also quantify the soiling contribution and power degradation of the photovoltaic cells on Mars through analysis of 4.95 Martian years of report-out power conditions from the Opportunity rover. We find that atmospheric dust suspended due to a weather event does not result in instantaneous settled dust on the PV cells. We calculate via autocorrelation function that the dust settling rate is approximately 21 Sols from atmospheric dust suspension. The findings presented here deliver a realistic approximation for the insolation values and subsequent PV power expected over time on the Martian surface thus informing future dust abatement systems.

McMillon-Brown, Lyndsey↗

USM3D Simulations for the Third AIAA Sonic Boom Prediction Workshop

The NASA USM3D flow solver was used to compute test cases for the Third AIAA Sonic Boom Prediction Workshop (SBPW3). The test cases included a nearfield biconvex shock-plume interaction wind tunnel model and the C608 Low Boom Flight Demonstrator. Numerical simulations were conducted on the mixed-element and tetrahedral grids provided by the workshop committee, as well as a family of grids generated by an in-house approach known as BoomGrid. The nearfield pressure signatures were extracted and propagated to the ground, and the perceived loudness levels on the ground were computed. The USM3D nearfield pressure signatures, corresponding ground signatures, and loudness levels on the ground were compared with that of mean values from other workshop participants. The effect of three flux-limiters on the accuracy of nearfield pressure signature prediction was investigated and results were compared with that of mean values from other workshop participants. The effect of using wall-function models for predicting nearfield pressure signatures was also evaluated. Results showed that the overpressure signatures extracted from the wall-modeled simulations and the wall-resolved simulations are in good agreement. The use of wall functions (wall-modeled simulations) allowed for approximately 23% savings in the time to solution and yielded comparable perceived loudness levels on the ground to the wall-resolved Reynolds-averaged Navier Stokes simulations.

Supersonics↗

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗

Thermal diffuse X-ray scattering in simple metals

Calculations are reported for the ionic structure factor and X-ray scattering cross section of sodium (at T=0 K and 90 K) and lithium (both isotopes at T=0 K) within the harmonic approximation. An evaluation of the appropriate displacement- displacement correlation function by the special point method circumvents the need for a multiphonon expansion. In the case of sodium, the structure in the one-phonon scattering was straightforwardly accounted for, and an approximate expansion was obtained for all multiphonon scattering. By treating core and conduction electrons on an equal footing, it is shown that information on the conduction electron system is present in the forward scattering component. In lithium the one-phonon cross section at small angles aids in the determination of the effective electron-ion interaction.

Straus, D. M.↗

Modal propagation angles in a cylindrical duct with flow and their relation to sound radiation

The main emphasis is upon the propagation angle with respect to the duct axis and its relation to the far-field acoustic radiation pattern. When the steady flow Mach number is accounted for in the duct, the propagation angle in the duct is shown to be coincident with the angle of the principal lobe of far-field radiation obtained using the Wiener-Hopf technique. Different Mach numbers are allowed within the duct and in the external field. For static tests with a steady flow in an inlet but with no external Mach number the far-field radiation pattern is shifted considerably toward the inlet axis when compared to zero Mach number radiation theory. As the external Mach number is increased the noise radiation pattern is shifted away from the inlet axis. The theory is developed using approximations for sound propagation in circular ducts. An exact analysis using Hankel function solutions for the zero Mach number case is given to provide a check of the simpler approximate theory.

Rice, E. J.↗

Modified Coulomb-Dipole Theory for 2e Photoionization

In the light of recent experiment on 2e photoionization of Li near threshold, we have considered a modification of the Coulomb-dipole theory, retaining the basic assumption that the threshold is dominated by asymmetric events in phase space [implies r(sub 1), k(sub 1)) greater than or equal to 2(r(sub 2), k(sub )]. In this region [in a collinear model, 2/r(sub 12) approached + 2/(r(sub 1)+r(sub 2)] the interaction reduces to V(rIsub 1) is greater than or equal to 2r(sub 2) is identically equal to [(-Z/r(sub 2)-(A-1)/r(sub 1)] + [(-2r(sub 2)/r(sub 1 exp 2)] is identically equal to V(sub c)+[V(sub d)]. For two electron emission Z = 2, thus both electrons see a Coulomb potential (V(sub c)) asymptotically, albeit each seeing a different charge. The residual potential (V(sub d)) is dipole in character. Writing the total psi = psi (sub c) + psi(sub d) = delta psi, and noting that. (T+V(sub c)-E)psy(sub c) = 0 and (T+V(sub c))psi(sub d) = 0 can be solved exactly, we find, substituting psi into the complete Schrod. Eq., that delta psi = -(H-E)(exp -1)(V(sub d) psi(sub 0)+V(sub c psi (sub 1). Using the fact that the absolute value of V(sub c) is much more than the absolute value of V(sub d) in almost all of configuration space, we can replace H by H(sub 0) in 9H-E)(exp -1) to obtain an improved approximation psi (improved) = psi(sub c) + psi(sub d) -(H(sub 0)-E)(exp -1) (V(sub c) psi (sub 0) + V(sub c) psi(sub 1). Here's the Green's function (H(sub 0)-E)(exp -1), can be exhibited explicitly, but the last term in psi (improved) is small, compared to the first two terms. Inserting them into the transition matrix element, which one handles in the usual way, we obtain in the limit E approaches 0, the threshold law: Q(E) alpha E + M(E)E(exp 5/4) + higher order (Eq. 1a). The modulation function, M(E), is a well-defined (but very non-trivial integral, but it is expected to be well approximated by a sinusoidal function containing a dipole phase term (M(E) = c sin[alpha log (E) + micron] (Eq. 1b). Experimental results show definite modulations, and are well fitted by Eqs (1).

Source record↗

Vapor ingestion in a cylindrical tank with a concave elliptical bottom

An approximate analytical technique is presented for estimating the liquid residual in a tank of arbitrary geometry due to vapor ingestion at any drain rate and acceleration level. The bulk liquid depth at incipient pull-through is defined in terms of the Weber and Bond numbers and two functions that describe the fluid velocity field and free surface shape appropriate to the tank geometry. Numerical results are obtained for the Centaur LH2 tank using limiting approximations to these functions.

Klavins, A.↗

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗

NanoSail-D: The First Flight Demonstration of Solar Sails for Nanosatellites

The NanoSail-D mission is currently scheduled for launch onboard a Falcon Launch Vehicle in the late June 2008 timeframe. The NanoSail-D, a CubeSat-class satellite, will consist of a sail subsystem stowed in a Cubesat 2U volume integrated with a CubeSat 1U volume bus provided by the NASA Ames Research Center (ARC). Shortly after deployment of the NanoSail-D from a Poly Picosatellite Orbital Deployer (P-POD) ejection system, the solar sail will deploy and mission operations will commence. This demonstration flight has two primary mission objectives: 1) to successfully stow and deploy the sail and 2) to demonstrate de-orbit functionality. Given a nearterm opportunity for launch, the project was met with the challenge of delivering the flight hardware in approximately six months, which required a significant constraint on flight system functionality. As a consequence, passive attitude stabilization will be achieved using permanent magnets to de-tumble and orient the body with the magnetic field lines and then rely on atmospheric drag to passively stabilize the sailcraft in an essentially maximum drag attitude. This paper will present an introduction to solar sail propulsion systems, overview the NanoSail-D spacecraft, describe the performance analysis for the passive attitude stabilization, and present a prediction of flight data results from the mission.

Whorton, Mark↗

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z = 4 to 8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1- 2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 ×10(exp 6) Mpc(exp 3) over this epoch, allowing us to perform a robust search for bright (M(sub UV) less than −21) and faint (M(sub UV) = −18) galaxies. We select galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5 less than z less than 8.5, with more than 1000 galaxies at z of approximately 6 - 8. We measure both a stepwise luminosity function for galaxies in our redshift samples, as well as a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z of greater than or equal to 6. Our bestfit value of the characteristic magnitude M* is consistent with −21 at z of greater than or equal to 5, different than that inferred based on previous trends at lower redshift. At z = 8, a single power-law provides an equally good fit to the UV luminosity function, while at z = 6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift, though a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to M(sub UV) = −17, correcting for dust attenuation, and find that the SFR density declines proportionally to (1 + z)((exp −4.3)(+/-)(0.5)) at z greater than 4, consistent with observations at z greater than or equal to 9. Our observed luminosity functions are consistent with a reionization history that starts at redshift of approximately greater than 10, completes at z greater than 6, and reaches a midpoint (x(sub HII) = 0.5) at 6.7 less than z less than 9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z similar to 10 galaxies.

Galaxy↗

Optimal thermionic energy conversion with established electrodes for high-temperature topping and process heating

Applied research-and-technology (ART) work reveals that optimal thermionic energy conversion (TEC) with approximately 1000 K to approximately 1100 K collectors is possible using well established tungsten electrodes. Such TEC with 1800 K emitters could approach 26.6% efficiency at 27.4 W/sq cm with approximately 1000 K collectors and 21.7% at 22.6 W/sq cm with approximately 1100 K collectors. These performances require 1.5 and 1.7 eV collector work functions (not the 1 eV ultimate) with nearly negligible interelectrode losses. Such collectors correspond to tungsten electrode systems in approximately 0.9 to approximately 6 torr cesium pressures with 1600 K to 1900 K emitters. Because higher heat-rejection temperatures for TEC allow greater collector work functions, interelectrode loss reduction becomes an increasingly important target for applications aimed at elevated temperatures. Studies of intragap modifications and new electrodes that will allow better electron emission and collection with lower cesium pressures are among the TEC-ART approaches to reduced interelectrode losses. These solutions will provide very effective TEC to serve directly in coal-combustion products for high-temperature topping and process heating. In turn this will help to use coal and to use it well.

Morris, J. F.↗

A minimum entropy principle in the gas dynamics equations

Let u(x bar,t) be a weak solution of the Euler equations, governing the inviscid polytropic gas dynamics; in addition, u(x bar, t) is assumed to respect the usual entropy conditions connected with the conservative Euler equations. We show that such entropy solutions of the gas dynamics equations satisfy a minimum entropy principle, namely, that the spatial minimum of their specific entropy, (Ess inf s(u(x,t)))/x, is an increasing function of time. This principle equally applies to discrete approximations of the Euler equations such as the Godunov-type and Lax-Friedrichs schemes. Our derivation of this minimum principle makes use of the fact that there is a family of generalized entrophy functions connected with the conservative Euler equations.

Tadmor, E.↗

Cosmological parameters and evolution of the galaxy luminosity function

The relationship between the observed distribution of discrete sources of a flux limited sample, the luminosity function of these sources, and the cosmological model is discussed. It is stressed that some assumptions about the form and evolution of the luminosity function must be made in order to determine the cosmological parameters from the observed distribution of sources. Presented is a method to test the validity of these assumptions using the observations. It is shown how, using higher moments of the observed distribution, one can determine, independently of the cosmological model, all parameters of the luminosity function except those describing evolution of the density and the luminosity of the luminosity function. These methods are applied to the sample of approximately 1000 galaxies recently used by Loh and Spillar to determine a value of the cosmological density parameter Omega approx = 1. It is shown that the assumptions made by Loh and Spillar about the luminosity function are inconsistent with the data, and that a self-consistent treatment of the data indicates a lower value of Omega approx = 0.2 and a flatter luminosity function. It should be noted, however, that incompleteness in the sample could cause a flattening of the luminosity function and lower the calculated value of Omega and that uncertainty in the values of these parameters due to random fluctuations is large.

Caditz, David↗

Cosmological parameters and evolution of the galaxy luminosity function

The relationship between the observed distribution of discrete sources of a flux limited sample, the luminosity function of these sources, and the cosmological model is discussed. It is stressed that some assumptions about the form and evolution of the luminosity function must be made in order to determine the cosmological parameters from the observed distribution of sources. Presented is a method to test the validity of these assumptions using the observations. It is shown how, using higher moments of the observed distribution, one can determine, independently of the cosmological model, all parameters of the luminosity function except those describing evolution of the density and the luminosity of the luminosity function. These methods are applied to the sample of approximately 1000 galaxies recently used by Loh and Spillar to determine a value of the cosmological density parameter Omega approx = 1. It is shown that the assumptions made by Loh and Spillar about the luminosity function are inconsistent with the data, and that a self-consistent treatment of the data indicates a lower value of Omega approx = 0.2 and a flatter luminosity function. It should be noted, however, that incompleteness in the sample could cause a flattening of the luminosity function and lower the calculated value of Omega and that uncertainty in the values of these parameters due to random fluctuations is large.

Caditz, David↗

Metal interfaces - Adhesive energies and electronic barriers

We report a fully self-consistent calculation of the electron number density, barrier height and adhesive energy as a function of separation in an aluminum-aluminum (100) contact. The local density approximation is used for exchange and correlation. The electron number density and barrier heights are strong functions of the separation. The range of strong chemical bonding is about 0.2 nm.

Ferrante, J.↗

Avalanche-diode oscillator circuit with tuning at multiple frequencies

Detailed theoretical analysis of three different modes or types of high efficiency oscillation in a PIN diode are presented. For the TRAPATT mode in a PIN diode, it is shown that a traveling avalanche zone is not necessary to generate a dense trapped plasma. An economical computer program for TRAPATT oscillations in a PIN diode is described. Typical results of diode power, dc-to-RF conversion efficiency, and required circuit impedances are presented for several different current waveforms. A semianalytical solution for a second type of high efficiency mode in a PIN diode is derived assuming a rectangular current waveform. A quasi-static approximation is employed to derive a semianalytical solution for the voltage across a PIN diode in a third mode, where avalanching occurs during a major portion of a half cycle. Calculations for this mode indicate that the power increases proportionally to the magnitude of the drive current with a small decrease in efficiency relative to the ordinary TRAPATT mode. An analytical solution is also given for a PIN diode, where it is assumed that the ionization coefficient is a step function. It is shown that the step-ionization approximation permits one to draw possible patterns of avalanche region in the depletion layer as a function of time. A rule governing admissible patterns is derived and an example solution given for one admissible pattern. Preliminary experimental results on the high-efficiency oscillations are presented and discussed. Two different experimental circuits, which used channel-dropping filters to provide independent harmonic tuning, are described. Simpler circuits used to produce high-efficiency oscillations are discussed. Results of experiments using inexpensive Fairchild FD300 diodes are given.

Parker, D.↗