Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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,189 records · Page 66

Solid earth as a recycling systems and the lateral growth of Precambrian North America

If plotted on mass vs time diagrams, geologic entities (for example, continental and oceanic crust, sediments, and mineral resources) display an exponential (power law) relationship, with entity per unit time increasing toward the present. This relationship is consistent with the concept of recycling and can be simulated mathematically. The approach is based on the plate tectonic theory and considers area-age or mass-age distributions of crystalline basement and sediments for major global tectonic realms. Each tectonic realm is characterized by a specific lifespan, which is an inverse function of its recycling rate. The estimated average half-area of half-mass ages are given. The corresponding parameters for continental crust are 690 Ma for K/Ar, and approximately 1200 Ma for Rb/St and U-Th/Pb dating pairs. Tectonic diversity preserved in the geologic record is therefore a function of time, with oceanic tectonic realms, because of their rapid recycling, underrepresented in the rocks older than approximately 300 Ma. The Sm/Nd isotopic systematic of sediments suggest that, for a near steady-state post-Archean sedimentary mass, recycling is approximately 90 + or - 5 percent cannibalistic. This yields an estimated upper limit on crust-mantle exchange via sediment subduction of approximately 1.1 + or - 0.5 x 10 g a(sup -1) considerably less than demanded by isotopic constraints. The discrepancy may indicate the existence of additional loci, such as orogenic belts, for significant crust-mantle interaction.

Veizer, Jan↗

Metal-insulator transition and magnetism of SU(3) fermions in the square lattice

We study the SU(3) symmetric Fermi-Hubbard model (FHM) in the square lattice at 1/3-filling using numerically exact determinant quantum Monte Carlo and numerical linked-cluster expansion techniques. We present the different regimes of the model in the T–U plane, which are characterized by local and short-range correlations, and capture signatures of the metal-insulator transition and magnetic crossovers. These signatures are detected as the temperature scales characterizing the rise of the compressibility, and an interaction-dependent change in the sign of the diagonal spin-spin correlation function. The analysis of the compressibility estimates the location of the metal-insulator quantum critical point at U c /t ~ 6, and provides a temperature scale for observing Mott physics at finite T. Furthermore, from the analysis of the spin-spin correlation function we observe that for U/t ≳ 6 and T ~ J = 4⁢t 2 /U there is a development of a short-range two-sublattice (2SL) antiferromagnetic structure, as well as an emerging three-sublattice (3SL) antiferromagnetic structure as the temperature is lowered below T/J ≲ 0.57. This crossover from 2SL to 3SL magnetic ordering agrees with Heisenberg limit predictions, and has observable effects on the density of on-site pairs. Finally, we describe how the features of the regimes in the T–U plane can be explored with alkaline-earth-like atoms in optical lattices with currently achieved experimental techniques and temperatures. Furthermore, the results discussed in this paper provide a starting point for the exploration of the SU(3) FHM upon doping.

74 ATOMIC AND MOLECULAR PHYSICS↗

Boundedness regions of discrete-time dynamic systems

Techniques for obtaining quantitative information about boundedness properties are developed and applied to the sampled-data control of satellite attitude with quantization. Relevant stability concepts are introduced as a series of definitions, and interrelationships between various definitions are discussed. The boundedness regions are estimated by means of quadratic Liapunov functions, and a sufficient condition for the existence of a boundedness region is given for a certain class of systems. A quadratic Liapunov function is applied to the Lur'e-Postinkov class of systems, where the linear part of the system is not asymptotically stable and the quantizer represents the nonlinear characteristic. A numerical calculation of the region of boundedness estimates is performed for satellite attitude control and is compared with simulation results. It is tentatively concluded that the Liapunov results may be good and that simulation results may be difficult to interpret and time-consuming to generate. The Lur'e-based technique yields estimates of regions of absolute boundedness, but at the cost of greater analytical complexity.

Siljak, D.↗

Satellite thermal modeling using the illuminographic camera

The application is discussed of illuminographic camera techniques to the measurement of radiative form factors for satellite modeling applications. The equation for the illuminographic mirror is derived, and the advantages of the method are discussed. Proximity error for the illuminographic mirror is defined, and estimates of this error as a function of the scale factor of the model are given, together with a method for halving this error. Designs for illuminographic mirrors are presented: the first provides a further threefold reduction in proximity error; the second permits the estimation of thermal radiation transfer coefficients to interior regions of otherwise inaccessible volumes. Finally, reflex Polaroid camera techniques appropriate to the relatively small size of the illuminographic mirror are described.

Wiesner, S. J.↗

Maximizing TDRS Command Load Lifetime

The GNC software onboard ISS utilizes TORS command loads, and a simplistic model of TORS orbital motion to generate onboard TORS state vectors. Each TORS command load contains five "invariant" orbital elements which serve as inputs to the onboard propagation algorithm. These elements include semi-major axis, inclination, time of last ascending node crossing, right ascension of ascending node, and mean motion. Running parallel to the onboard software is the TORS Command Builder Tool application, located in the JSC Mission Control Center. The TORS Command Builder Tool is responsible for building the TORS command loads using a ground TORS state vector, mirroring the onboard propagation algorithm, and assessing the fidelity of current TORS command loads onboard ISS. The tool works by extracting a ground state vector at a given time from a current TORS ephemeris, and then calculating the corresponding "onboard" TORS state vector at the same time using the current onboard TORS command load. The tool then performs a comparison between these two vectors and displays the relative differences in the command builder tool GUI. If the RSS position difference between these two vectors exceeds the tolerable lim its, a new command load is built using the ground state vector and uplinked to ISS. A command load's lifetime is therefore defined as the time from when a command load is built to the time the RSS position difference exceeds the tolerable limit. From the outset of TORS command load operations (STS-98), command load lifetime was limited to approximately one week due to the simplicity of both the onboard propagation algorithm, and the algorithm used by the command builder tool to generate the invariant orbital elements. It was soon desired to extend command load lifetime in order to minimize potential risk due to frequent ISS commanding. Initial studies indicated that command load lifetime was most sensitive to changes in mean motion. Finding a suitable value for mean motion was therefore the key to achieving this goal. This goal was eventually realized through development of an Excel spreadsheet tool called EMMIE (Excel Mean Motion Interactive Estimation). EMMIE utilizes ground ephemeris nodal data to perform a least-squares fit to inferred mean anomaly as a function of time, thus generating an initial estimate for mean motion. This mean motion in turn drives a plot of estimated downtrack position difference versus time. The user can then manually iterate the mean motion, and determine an optimal value that will maximize command load lifetime. Once this optimal value is determined, the mean motion initially calculated by the command builder tool is overwritten with the new optimal value, and the command load is built for uplink to ISS. EMMIE also provides the capability for command load lifetime to be tracked through multiple TORS ephemeris updates. Using EMMIE, TORS command load lifetimes of approximately 30 days have been achieved.

Brown, Aaron J.↗

Estimation, Navigation and Control of Multi-Rotor Drones in an Urban Wind Field

The paper presents an on-board estimation, navigation and control architecture for multi-rotor drones flying in urban environment. It consists of adaptive algorithms to estimate vehicle's aerodynamic drag coefficients with respect to still air and the urban wind components along the flight trajectory, with guaranteed fast and reliable convergence to the true values; navigation algorithms to generate feasible trajectories between given way-points that take into account the estimated wind; and of control algorithms to track the generated trajectories as long as the vehicle retains sufficient number of functioning rotors capable of compensating for the estimated wind. All components of this on-board system are computationally effective and are intended for a real time implementation. The algorithms were tested in simulations.

Stepanyan, Vahram↗

Quantifying radiation quality for space relevant radiation types: Fitting excess risk models to three combined HZE-irradiated mouse datasets

Radiation health risks are predominantly derived from low linear energy transfer (LET) terrestrial exposures; however, space radiation includes exposure to high-LET and high-charge, high-energy (HZE) particles. Accurately quantifying the differences in radiation quality between the space and terrestrial radiation environments is important for assessing and predicting health risks for astronauts. Weil et al. 2009 and 2014 used two different inbred mouse strains to study differences in hepatocellular carcinoma (HCC) tumorigenesis after exposures to low- and high- LET radiation. More recently, Edmundson et al. 2020 provided valuable new tumor data in outbred mice that were exposed to low- and high-LET radiation. The present study aims to rigorously investigate a relative biological effectiveness (RBE) factor by leveraging the HCC tumor data from Weil et al. 2009, Weil et al. 2014, and Edmundson et al. 2020. The three experiments were similarly designed, allowing the raw data to be combined into a pooled dataset to estimate excess relative risk (ERR) and excess absolute risk (EAR) models using Bayesian Poisson regression. These effect estimates from the pooled data provide greater power to calculate a data driven RBE. Extensive sensitivity analyses test the robustness of RBE estimates to various model assumptions. The following questions will be explored through the sensitivity analyses: • Is the shape of the dose response different for low-LET radiation and HZE radiation, indicating that RBE is a function of dose? • Does attained age modify the effect estimates differently for low-LET radiation and HZE radiation, indicating RBE is a function of attained age? • Are the effect estimates and RBE estimates different for inbred mouse strains and outbred mouse strains? • Do assumptions about differences in ERR models and EAR models change the estimated RBE? Additional studies would be needed to validate the findings from these exploratory analyses.

Lori J. Chappell↗

Photon Radiation Scatter from Heterogeneous Shields using Green's Functions

The effect of shielding on ionizing photon radiation can be estimated using radiation transport simulations. This report covers the methodology and implementation of using Green's Functions to pre-compute this effect, which allows the radiation field exiting a variety of shielding configurations to be quickly computed. It also covers a weighting function that makes a relatively small pre-computed library applicable to a large variety of heterogeneous shields. The method enables rapid computation of the intensity versus energy for scattered radiation exiting a variety of shield materials and thicknesses without running a full transport simulation.

61 RADIATION PROTECTION AND DOSIMETRY↗

On estimating the curvature attributes and strain invariants of deformed surface through radial basis functions

Numerical methods in approximation of derivatives on regular or irregular gridded and scattered data points are of significance in the study of numerous problems in engineering and geoscience. A mathematically accurate and stable method is assessed to compute the differential quantities, while being flexible and easy to implement. Strain-invariant parameters in deformation analysis and curvature attributes in geometrical analysis of approximated surfaces are independent quantities which could be derived numerically through radial basis functions (RBFs) for scattered data points. Approximation of a function or its derivatives through RBFs is highly correlated to its shape parameter which depends on the number and distribution of the particles on the support domain. A procedure to find the optimal shape parameter for RBFs in the influence domain of each nodal point together with the presentation of a manner of calculating the curvature attributes and strain parameters through Gaussian RBFs is proposed here. Results indicate a significant improvement in the approximation accuracy of strain parameters and curvature attributes, while this approximation is more stable when RBFs with optimal shape parameter are implemented rather than the traditional moving least squares.

97 MATHEMATICS AND COMPUTING↗

Estimating Basic Preliminary Design Performances of Aerospace Vehicles

Aerodynamics and Performance Estimation Toolset is a collection of four software programs for rapidly estimating the preliminary design performance of aerospace vehicles represented by doing simplified calculations based on ballistic trajectories, the ideal rocket equation, and supersonic wedges through standard atmosphere. The program consists of a set of Microsoft Excel worksheet subprograms. The input and output data are presented in a user-friendly format, and calculations are performed rapidly enough that the user can iterate among different trajectories and/or shapes to perform "what-if" studies. Estimates that can be computed by these programs include: 1. Ballistic trajectories as a function of departure angles, initial velocities, initial positions, and target altitudes; assuming point masses and no atmosphere. The program plots the trajectory in two-dimensions and outputs the position, pitch, and velocity along the trajectory. 2. The "Rocket Equation" program calculates and plots the trade space for a vehicle s propellant mass fraction over a range of specific impulse and mission velocity values, propellant mass fractions as functions of specific impulses and velocities. 3. "Standard Atmosphere" will estimate the temperature, speed of sound, pressure, and air density as a function of altitude in a standard atmosphere, properties of a standard atmosphere as functions of altitude. 4. "Supersonic Wedges" will calculate the free-stream, normal-shock, oblique-shock, and isentropic flow properties for a wedge-shaped body flying supersonically through a standard atmosphere. It will also calculate the maximum angle for which a shock remains attached, and the minimum Mach number for which a shock becomes attached, all as functions of the wedge angle, altitude, and Mach number.

Luz, Paul L.↗

Robust and scalable uncertainty estimation with conformal prediction for machine-learned interatomic potentials

Abstract Uncertainty quantification (UQ) is important to machine learning (ML) force fields to assess the level of confidence during prediction, as ML models are not inherently physical and can therefore yield catastrophically incorrect predictions. Established a-posteriori UQ methods, including ensemble methods, the dropout method, the delta method, and various heuristic distance metrics, have limitations such as being computationally challenging for large models due to model re-training. In addition, the uncertainty estimates are often not rigorously calibrated. In this work, we propose combining the distribution-free UQ method, known as conformal prediction (CP), with the distances in the neural network’s latent space to estimate the uncertainty of energies predicted by neural network force fields. We evaluate this method (CP+latent) along with other UQ methods on two essential aspects, calibration, and sharpness, and find this method to be both calibrated and sharp under the assumption of independent and identically-distributed (i.i.d.) data. We show that the method is relatively insensitive to hyperparameters selected, and test the limitations of the method when the i.i.d. assumption is violated. Finally, we demonstrate that this method can be readily applied to trained neural network force fields with traditional and graph neural network architectures to obtain estimates of uncertainty with low computational costs on a training dataset of 1 million images to showcase its scalability and portability. Incorporating the CP method with latent distances offers a calibrated, sharp and efficient strategy to estimate the uncertainty of neural network force fields. In addition, the CP approach can also function as a promising strategy for calibrating uncertainty estimated by other approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗

A statistical analysis of elevated temperature gravimetric cyclic oxidation data of 36 Ni- and Co-base superalloys based on an oxidation attack parameter

A large body of high temperature cyclic oxidation data generated from tests at NASA Lewis Research Center involving gravimetric/time values for 36 Ni- and Co-base superalloys was reduced to a single attack parameter, K(sub a), for each run. This K(sub a) value was used to rank the cyclic oxidation resistance of each alloy at 1000, 1100, and 1150 C. These K(sub a) values were also used to derive an estimating equation using multiple linear regression involving log(sub 10)K(sub a) as a function of alloy chemistry and test temperature. This estimating equation has a high degree of fit and could be used to predict cyclic oxidation behavior for similar alloys and to design an optimum high strength Ni-base superalloy with maximum high temperature cyclic oxidation resistance. The critical alloy elements found to be beneficial were Al, Cr, and Ta.

Barrett, Charles A.↗

Constrained optimization using design of experiment surfaces

An algorithm for solving constrained optimization problems is presented. First, design of experiment techniques are used to survey the design space. After evaluating the objective and constraint functions, as specified by Taguchi orthogonal arrays, analytical models of these functions are generated using a least-squares regression analysis. Next, a nonlinear programming package is used to optimize the analytical model. Based on the optimization information, the design space is reduced so as to close in around the minimum, and the entire procedure is repeated until convergence. An important feature of the algorithm is that function gradients are not required; therefore, for problems in which gradients would have to be estimated using finite-differences the number of function evaluations required for the optimization is significantly reduced, when compared with traditional nonlinear programming techniques. In addition, there is no requirement that the gradients must be smooth and continuous.

Bolt, Marvin Vance↗

Initial mass functions from ultraviolet stellar photometry: A comparison of Lucke and Hodge OB associations near 30 Doradus with the nearby field

UV stellar photometry is presented for 1563 stars within a 40 minutes circular field in the Large Magellanic Cloud (LMC), excluding the 10 min x 10 min field centered on R136 investigated earlier by Hill et al. (1993). Magnitudes are computed from images obtained by the Ultraviolet Imaging Telescope (UIT) in bands centered at 1615 A and 2558 A. Stellar masses and extinctions are estimated for the stars in associations using the evolutionary models of Schaerer et al. (1993), assuming the age is 4 Myr and that the local LMC extinction follows the Fitzpatrick (1985) 30 Dor extinction curve. The estimated slope of the initial mass function (IMF) for massive stars (greater than 15 solar mass) within the Lucke and Hodge (LH) associations is Gamma = -1.08 +/- 0.2. Initial masses and extinctions for stars not within LH associations are estimated assuming that the stellar age is either 4 Myr or half the stellar lifetime, whichever is larger. The estimated slope of the IMF for massive stars not within LH associations is Gamma = -1.74 +/- 0.3 (assuming continuous star formation), compared with Gamma = -1.35, and Gamma = -1.7 +/- 0.5, obtained for the Galaxy by Salpeter (1955) and Scalo (1986), respectively, and Gamma = -1.6 obtained for massive stars in the Galaxy by Garmany, Conti, & Chiosi (1982). The shallower slope of the association IMF suggests that not only is the star formation rate higher in associations, but that the local conditions favor the formation of higher mass stars there. We make no corrections for binaries or incompleteness.

Hill, Jesse K.↗

Intra-day variability of three Seyfert galaxies measured with XMM-Newton

We analyze the variability of three Seyfert galaxies on intra-day timescales. We have analyzed in a uniform manner the 38 longest (>30 ks) observations made for NGC 4051, MCG-06–30–15 and NGC 4151 by XMM-Newton between 2000 and 2015. The nuclei were quite active during essentially all of these observations and the overall X-ray fluxes (0.3–10 keV) varied by an order of magnitude. Most of the observations do appear to show characteristic timescales, estimated through their auto-correlation functions, ranging between ~2.9 ks and ~45.3 ks. The hard (2–10 keV) and soft (0.3–2.0 keV) bands are very well correlated but consideration of their hardness ratios shows that the sources typically soften during flares. We also provide new estimates of the central black hole masses for these three AGNs that support the hypothesis that Narrow Line Seyfert 1 galaxies have relatively small central black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Shape Servoing of Deformable Objects using Adaptive Deformation Model Estimation

In this paper, we propose an adaptive shape servoing method to deform a soft object into a desired 3-D shape. The high dimensional representation and the unknown deformation properties of the soft object pose a challenge to actively manipulate its shape. To address this issue, we develop a method to compute the deformation Jacobian matrix in real-time. The Jacobian is estimated using a set of basis functions and its corresponding parameters to capture the dynamics of the system and relate the applied input motion to changes in the soft object's shape. An integral concurrent learning (ICL) based adaptive update law is derived using Lyapunov analysis to estimate the deformation parameters and prove its convergence. A physics-based simulation is used to validate the proposed method and controller by performing manipulation tasks with different desired configurations. The performance is compared with a standard gradient update law to demonstrate the accuracy and robustness of our approach.

Vrithik Raj Guthikonda↗

Integration, Validation, and Application of a PV Snow Coverage Model in SAM

Due to the increasing deployment of PV systems in snowy climates, there is significant interest in a method capable of estimating PV losses resulting from snow coverage that has been verified for a variety of system designs and locations. Many independent snow coverage models have been developed over the last 15 years; however, there has been very little effort verifying these models beyond the system designs and locations on which they were based. Moreover, major PV modeling software products have not yet incorporated any of these models into their workflows. In response to this deficiency, we have integrated the methodology of the snow model developed in the paper by Marion et al. (2013) into the National Renewable Energy Laboratory's (NREL) System Advisor Model (SAM). In this work, we describe how the snow model is implemented in SAM and we discuss our demonstration of the model's effectiveness at reducing error in annual estimations for three PV arrays. Next, we use this new functionality in conjunction with a long term historical data set to estimate average snow losses across the United States for two typical PV system designs. The open availability of the snow loss estimation capability in SAM to the PV modeling community, coupled with our results of the nationwide study, will better equip the industry to accurately estimate PV energy production in areas affected by snowfall.

14 SOLAR ENERGY↗

A Drift‐Kinetic Method for Obtaining Gradients in Plasma Properties From Single‐Point Distribution Function Data

Abstract In this paper, we derive a new drift‐kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle‐in‐cell simulations.

Wetherton, B. A.↗