SEARCH · Engineering Papers
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.
The author discussed the theoretical and computational results for spline based approximation schemes used in parameter estimation algorithms for distributed systems. Specific applications include beam-like structures described by the Euler-Bernoulli and Timoshenko theories and antenna surfaces such as that in the deployable Maypole Hoop/Column model.
ABSTRACT In this work, we demonstrate how differentiable stochastic sampling techniques developed in the context of deep reinforcement learning can be used to perform efficient parameter inference over stochastic, simulation-based, forward models. As a particular example, we focus on the problem of estimating parameters of halo occupation distribution (HOD) models that are used to connect galaxies with their dark matter haloes. Using a combination of continuous relaxation and gradient re-parametrization techniques, we can obtain well-defined gradients with respect to HOD parameters through discrete galaxy catalogue realizations. Having access to these gradients allows us to leverage efficient sampling schemes, such as Hamiltonian Monte Carlo, and greatly speed up parameter inference. We demonstrate our technique on a mock galaxy catalogue generated from the Bolshoi simulation using a standard HOD model and find near-identical posteriors as standard Markov chain Monte Carlo techniques with an increase of ∼8× in convergence efficiency. Our differentiable HOD model also has broad applications in full forward model approaches to cosmic structure and cosmological analysis.
We calculate field isotope shifts for nobelium atoms using nuclear charge distributions which come from different nuclear models.We demonstrate that comparing calculated isotope shifts with experiment can serve as a testing ground for nuclear theories. It also provides a way of extracting parameters of nuclear charge distribution beyond nuclear root mean square (rms) radius, e.g., parameter of quadrupole deformation β. The measurements of at least two atomic transitions is needed to disentangle the contributions of the changes in deformation and nuclear rms radius into field isotopic shift. We argue that a previous interpretation of the isotope measurements in terms of δ < r 2 > between 252,254 No isotopes should be amended when nuclear deformation is taken into account. As a result, we calculate isotope shifts for other known isotopes and for hypothetically metastable isotope 286 No for which the predictions of nuclear models differ substantially.
Treated jointly, two recent results imply that the distribution of winds in the polar ionosphere should change as a function of the direction of the interplanetary magnetic field. (1) From the motions of chemically released ion and neutral clouds, it is apparent that neutral winds in the high latitude ionosphere are driven principally by ion drag forces. (2) OGO-6 electric field measurements have demonstrated that there are definite relationships between the time-latitude distribution of ionospheric plasma convection and interplanetary magnetic field parameters, and also that the distribution is most sensitive to the azimuthal angle of the interplanetary field. Although direct neutral wind to interplanetary magnetic field comparisons are not available, logic clearly implies a close relationship. The lower altitude, meteorological effects of these externally driven ionospheric winds are not known. However, observations of infrasonic waves following sudden ionization enhancements indicate the existence of momentum transfer.
Two recent results imply that the distribution of winds in the polar ionosphere should change as a function of the direction of the interplanetary magnetic field. From the motions of chemically released ion and neutral clouds, it is apparent that neutral winds in the high latitude ionosphere are driven principally by ion drag forces. OGO-6 electric field measurements have demonstrated that there are definite relationships between the time latitude distribution of ionospheric plasma convection and interplanetary magnetic field parameters, and also that the distribution is most sensitive to the azimuthal angle of the interplanetary field. The lower altitude, meteorological effects of these externally driven ionospheric winds are not known. However, observations of infrasonic waves following sudden ionization enhancements indicate the existence of momentum transfer.
Treated jointly, two recent results imply that the distribution of winds in the polar ionosphere should change as a function of the direction of the interplanetary magnetic field. From the motions of chemically released ion and neutral clouds, it is apparent that neutral winds in the high-latitude ionosphere are driven principally by ion drag forces. OGO 6 electric field measurements demonstrate that there are definite relationships between the time/latitude distribution of ionospheric plasma convection and interplanetary magnetic field parameters, and also that the distribution is most sensitive to the azimuthal angle of the inter planetary field. Although direct neutral wind to interplanetary magnetic field comparisons are not available, logic clearly implies a close relationship. Observations of infrasonic waves following sudden ionization enhancements indicate the existence of momentum transfer.
The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).
Frequency distributions prepared from synoptic meteorological parameters associated with turbulent and nonturbulent conditions encountered by the XB-70 aircraft at altitudes generally between 12 and 20 km in the stratosphere show that turbulence is associated with certain ranges or critical values of the parameters. Certain combinations of the variables are associated with turbulence for at least one-fourth to one-third of the observations.
High latitude (>50N) natural wetlands emit a substantial amount of methane to the atmosphere, and are located in a region of amplified warming. Northern hemisphere high latitudes are characterized by cold climates, extensive permafrost, poor drainage, short growing seasons, and slow decay rates. Under these conditions, organic carbon accumulates in the soil, sequestering CO2 from the atmosphere. Methanogens produce methane from this carbon reservoir, converting stored carbon into a powerful greenhouse gas. Methane emission from wetland ecosystems depends on vegetation type, climate characteristics (e.g, precipitation amount and seasonality, temperature, snow cover, etc.), and geophysical variables (e.g., permafrost, soil type, and landscape slope). To understand how wetland methane dynamics in this critical region will respond to climate change, we have to first understand how wetlands themselves will change and therefore, what the primary controllers of wetland distribution and type are. Understanding these relationships permits data-anchored, physically-based modeling of wetland distribution and type in other climate scenarios, such as paleoclimates or future climates, a necessary first step toward modeling wetland methane emissions in these scenarios. We investigate techniques and datasets for predicting the distribution and type of high latitude (>50N) natural wetlands from a suite of geophysical and climate predictors. Hierarchical clustering is used to derive an empirical methane-centric wetland model. The model is applied in a multistep process first to predict the distribution of wetlands from relevant geophysical parameters, and then, given the predicted wetland distribution, to classify the wetlands into methane-relevant types using an expanded suite of climate and biogeophysical variables. As the optimum set of predictor variables is not known a priori, the model is applied iteratively, and each simulation is evaluated with respect to observed high-latitude wetlands.
Determination of the atomic-scale structures of certain fine-grained minerals using single-crystal X-ray diffraction (XRD) has been challenging because they commonly occur as submicron and nanocrystals in the geological environment. Synchrotron powder diffraction and scattering techniques are useful complementary methods for studying this type of minerals. In this review, we discussed three example studies investigated by combined methods of synchrotron radiation XRD and pair distribution function (PDF) techniques: (1) low-temperature cristobalite; (2) kaolinite; and (3) vernadite. Powder XRD is useful to determine the average structure including unit-cell parameters, fractional atomic coordinates, occupancies and isotropic atomic displacement parameters. X-ray/Neutron PDF methods are sensitive to study the local structure with anisotropic atomic displacement parameters (ADP). The results and case studies suggest that the crystal structure and high-quality ADP values can be obtained using the combined methods. The method can be useful to characterize crystals and minerals that are not suitable for single-crystal XRD.
This report describes the development and testing of a new method for extrapolating short-term creep rupture test data to predict long-term rupture strength. The goal of this work is to reduce the time required to qualify new materials for nuclear service by reducing the lead time required for dedicated, long-term material testing to establish key long-term material properties. The new approach described here uses a physics-based model to predict the long-term creep rupture strength of 316H stainless steel using only short-term test data. The key idea is to use Bayesian inference to find the statistical distribution of the model parameters that best explain the short-term rupture data. Because the model is physics-based these parameters are all microstructural quantities that can be measured through detailed material characterization experiments. The Bayesian prior distributions provide a means for incorporating this characterization data into the final model to improve the accuracy of the long-term model predictions. However, where such data is not available the process still produces an accurate model based on an uniformed prior. Our hypothesis is that this approach more accurately extrapolates the short-term test data when compared to current, empirical methods. The report proves this hypothesis using actual long-term rupture data available for 316H, including tests with rupture times greater than 200,000 hours. The general approach developed here could be applied to other materials and other time-dependent material properties. Applying this new technique to develop long-term qualified material properties, potentially in conjunction with other accelerated qualification approaches like staggered qualification test programs, could greatly reduce the time required to qualify new materials for nuclear service.
We report the discovery of a gas-giant planet orbiting a low-mass host star in the microlensing event MOA-bin-29 that occurred in 2006. We find five degenerate solutions with the planet/host-star mass ratio of q ∼ 10(exp−2). The Einstein radius crossing time of all models are relatively short (∼4–7 days), which indicates that the mass of host star is likely low. The measured lens-source proper motion is 5–9 mas/yr depending on the models. Since only finite source effects are detected, we conduct a Bayesian analysis in order to obtain the posterior probability distribution of the lens physical properties. As a result, we find the lens system is likely to be a gas-giant orbiting a brown dwarf or a very late M-dwarf in the Galactic bulge. The probability distributions of the physical parameters for the five degenerate models are consistent within the range of error. By combining these probability distributions, we conclude that the lens system is a gas giant with a mass of M(p)= 0.63(sub +1.13, sub 0.39) M(Jup) orbiting a brown dwarf with a mass of M(h)= 0.06 (sup +0.11, sub -0.04) M(ʘ) at a projected star–planet separation of r⊥ = 0.53(sup +0.89, sub -0.18)au . The lens distance is D(L)= 6.89(sup +1.19, sub -1.19) kpc, i.e., likely within the Galactic bulge.
A smooth empirical Bayes estimator is derived for the intensity parameter (hazard rate) in the Poisson distribution as used in life testing. The reliability function is also estimated either by using the empirical Bayes estimate of the parameter, or by obtaining the expectation of the reliability function. The behavior of the empirical Bayes procedure is studied through Monte Carlo simulation in which estimates of mean-squared errors of the empirical Bayes estimators are compared with those of conventional estimators such as minimum variance unbiased or maximum likelihood. Results indicate a significant reduction in mean-squared error of the empirical Bayes estimators over the conventional variety.
In order to quantitatively predict nano- as well as other particle-size distributions, one needs to have both a mathematical model and estimates of the parameters that appear in these models. Here, we show how one can use Bayesian inversion to obtain statistical estimates for the parameters that appear in recently derived mechanism-enabled population balance models (ME-PBM) of nanoparticle growth. The Bayesian approach addresses the question of “how well do we know our parameters, along with their uncertainties?.” The results reveal that Bayesian inversion statistical analysis on an example, prototype $\mathrm{lr(0)_n}$ nanoparticle formation system allows one to estimate not just the most likely rate constants and other parameter values, but also their SDs, confidence intervals, and other statistical information. Moreover, knowing the reliability of the mechanistic model's parameters in turn helps inform one about the reliability of the proposed mechanism, as well as the reliability of its predictions. Importantly, the paper can also be seen as a tutorial with the additional goal of achieving a “Gold Standard” Bayesian inversion ME-PBM benchmark that others can use as a control to check their own use of this methodology for other systems of interest throughout nature. Overall, the results provide strong support for the hypothesis that there is substantial value in using a Bayesian inversion methodology for parameter estimation in particle formation systems.
Aerosols are particles suspended in the atmosphere, ranging in size from nanometers to micrometers. Their size distribution affects key atmospheric processes, including nucleation, coagulation, scavenging, activation, and radiative properties (Seinfeld and Pandis 2016). Aerosol size distribution is a critical parameter in atmospheric science, influencing processes such as cloud formation and radiative forcing. Accurate representation of aerosol size distributions is essential for understanding their impact on climate, air quality, and human health. However, aerosol size and composition vary significantly across time and space due to meteorological conditions and natural or anthropogenic sources. (Wu and Boor 2021). Various instruments are used to measure aerosol size distributions, each with distinct principles, advantages, and limitations. This report begins with an in-depth overview of aerosol size-distribution comparison studies, focusing on the passive cavity aerosol spectrometer probe (PCASP), portable optical particle spectrometer (POPS), ultra-high-sensitivity aerosol spectrometer (UHSAS), aerodynamic particle sizer (APS), and scanning mobility particle sizer (SMPS). All of these instruments are used by the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Aerial Facility (AAF), which commissioned this comparison and report. The report evaluates the strengths and weaknesses of these instruments, highlights their applications, and discusses efforts to merge data from multiple instruments for comprehensive analysis.
Croplands play a critical role in regulating the energy and moisture exchanges between the land surface and atmosphere. However, the interactions between cropland and climate are usually poorly represented due to a lack of detailed representation in crop types and field management. Here, we coupled the Noah-MP-Crop model with the state-of-the-art Weather Research and Forecasting (WRF) model to explore and evaluate the crop growth dynamics in response to climate variations across Northeast China. The default parameters of the crop model were not exactly suitable for the agricultural ecosystems in Northeast China. The detailed cropland distribution, and crop phenology parameters including growing degree days (GDD) and planting (harvesting) date were first created using multi-source remote sensing products and reanalysis data, and was then successfully used to simulate the growth and yield for corn and soybean and associated energy exchanges. We also optimized and calibrated other crop parameters using the time-series of the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface products. The modified crop model substantially improved the simulation of crop growth, plant physiology, and biomass accumulation for both corn and soybean. Coupling the localized dynamic crop model into the WRF led to considerable decreases in the simulated mean-absolute-errors (MAEs) and biases of the leaf area index, evapotranspiration, and gross primary production compared with the MODIS observed values. Compared with the statistical yield from each province, the modified crop model underestimated the corn yield from 11.1% to 48.6%, whereas overestimated the soybean yield from 16.5% to 162.6%.
High temperature ceramic matrix composites (CMC) are being explored as viable candidate materials for hot section gas turbine components. These advanced composites can potentially lead to reduced weight, enable higher operating temperatures requiring less cooling and thus leading to increased engine efficiencies. However, these materials are brittle and show degradation with time at high operating temperatures due to creep as well as cyclic mechanical and thermal loads. In addition, these materials are heterogeneous in their make-up and various factors affect their properties in a specific design environment. Most of these advanced composites involve two- and three-dimensional fiber architectures and require a complex multi-step high temperature processing. Since there are uncertainties associated with each of these in addition to the variability in the constituent material properties, the observed behavior of composite materials exhibits scatter. Traditional material failure analyses employing a deterministic approach, where failure is assumed to occur when some allowable stress level or equivalent stress is exceeded, are not adequate for brittle material component design. Such phenomenological failure theories are reasonably successful when applied to ductile materials such as metals. Analysis of failure in structural components is governed by the observed scatter in strength, stiffness and loading conditions. In such situations, statistical design approaches must be used. Accounting for these phenomena requires a change in philosophy on the design engineer s part that leads to a reduced focus on the use of safety factors in favor of reliability analyses. The reliability approach demands that the design engineer must tolerate a finite risk of unacceptable performance. This risk of unacceptable performance is identified as a component's probability of failure (or alternatively, component reliability). The primary concern of the engineer is minimizing this risk in an economical manner. The methods to accurately determine the service life of an engine component with associated variability have become increasingly difficult. This results, in part, from the complex missions which are now routinely considered during the design process. These missions include large variations of multi-axial stresses and temperatures experienced by critical engine parts. There is a need for a convenient design tool that can accommodate various loading conditions induced by engine operating environments, and material data with their associated uncertainties to estimate the minimum predicted life of a structural component. A probabilistic composite micromechanics technique in combination with woven composite micromechanics, structural analysis and Fast Probability Integration (FPI) techniques has been used to evaluate the maximum stress and its probabilistic distribution in a CMC turbine stator vane. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Since the measured data for the ceramic matrix composite properties is very limited, obtaining a probabilistic distribution with their corresponding parameters is difficult. In case of limited data, confidence bounds are essential to quantify the uncertainty associated with the distribution. Usually 90 and 95% confidence intervals are computed for material properties. Failure properties are then computed with the confidence bounds. Best estimates and the confidence bounds on the best estimate of the cumulative probability function for R-S (strength - stress) are plotted. The methodologies and the results from these analyses will be discussed in the presentation.
Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year−1 (1Tg = 1012 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH4 year−1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of −34 to −46 Tg CH4 year−1, we estimated that the global net land methane emissions are 149–176 Tg CH4 year−1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.