SEARCH · Search NASA
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.
Exnihilo is a transport code system that contains both deterministic and stochastic (Monte Carlo) solvers. The deterministic solvers are contained in the Denovo package while the Monte Carlo solvers are implmented in the Shift package. This document explains the mathematical and numerical methods employed within Exnihilo.
recision timing distribution systems require determin istic and repeatable high speed serial data link latency. The EIC Timing Data Link will employ a specialized link layer protocol for deterministic multigigabit communication us ing AMD Ultrascale+ GTH and GTY transceivers. Link la tency must also be measurable to accurately compensate in real time for latency variations of the physical medium. The point-to-point link is full duplex and fully synchronous with a latency control algorithm to align the internal clocks of each system with picosecond resolution. Deterministic clock domain crossing between systems is ensured using static timing analysis of the elastic buffer control signals.
Simulations of quantum systems with Hamiltonian classical stochastic noise can be challenging when the noise exhibits temporal correlations over a multitude of time scales, such as for 1/f noise in solid-state quantum information processors. Here we present an approach for simulating Hamiltonian classical stochastic noise that performs temporal coarse-graining by effectively integrating out the high-frequency components of the noise. We focus on the case where the stochastic noise can be expressed as a sum of Ornstein-Uhlenbeck processes. Temporal coarse-graining is then achieved by conditioning the stochastic process on a coarse realization of the noise, expressing the conditioned stochastic process in terms of a sum of smooth, deterministic functions and bridge processes with boundaries fixed at zero, and performing the ensemble average over the bridge processes. For Ornstein-Uhlenbeck processes, the deterministic components capture all dependence on the coarse realization, and the stochastic bridge processes are not only independent but taken from the same distribution with correlators that can be expressed analytically, allowing the associated noise propagators to be precomputed once for all simulations. This combination of noise trajectories on a coarse time grid and ensemble averaging over bridge processes has practical advantages, such as a simple concatenation rule, that we highlight with numerical examples.
This is an investigation on two experimental datasets of laminar hypersonic flows, over a double-cone geometry, acquired in Calspan—University at Buffalo Research Center’s Large Energy National Shock (LENS)-XX expansion tunnel. These datasets have yet to be modeled accurately. A previous paper suggested that this could partly be due to mis-specified inlet conditions. The authors of this paper solved a Bayesian inverse problem to infer the inlet conditions of the LENS-XX test section and found that in one case they lay outside the uncertainty bounds specified in the experimental dataset. However, the inference was performed using approximate surrogate models. Here in this paper, the experimental datasets are revisited and inversions for the tunnel test-section inlet conditions are performed with a Navier–Stokes simulator. The inversion is deterministic and can provide uncertainty bounds on the inlet conditions under a Gaussian assumption. It was found that deterministic inversion yields inlet conditions that do not agree with what was stated in the experiments. An a posteriori method is also presented to check the validity of the Gaussian assumption for the posterior distribution. This paper contributes to ongoing work on the assessment of datasets from challenging experiments conducted in extreme environments, where the experimental apparatus is pushed to the margins of its design and performance envelopes.
A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].
Ankle sprains are the most common musculoskeletal injury in sport and military activity, despite existing prophylactic strategies. The purpose of this report was to develop a probabilistic simulation of lateral ankle sprains during single-limb drop landing, towards accelerating innovation in ankle sprain prevention. A deterministic, subject-specific musculoskeletal model was extended with automation and probabilistic distributions on sprain-related biomechanical factors. Probabilistic simulations were generated using traditional Monte Carlo techniques and the advanced mean value method, a more computationally-efficient approach. Predicted distributions of peak ankle joint rotations, velocities, and moments borne by supporting passive structures agreed favorably with the deterministic model and with reports of real sprain biomechanics. Parameter sensitivities identified that predictions were most strongly influenced by drop height, subtalar joint posture at contact, invertor/evertor co-activation, and passive ankle stiffness. The advanced mean value method predicted confidence bounds comparable to a 1000-trial Monte Carlo simulation, and required only 14 model evaluations and 4-min processing time. The extended probabilistic simulation may be useful to virtually test new prophylactic strategies for ankle sprains, and is made available for open-source use (https://simtk.org/projects/sprain-sim).
We present a probabilistic framework tailored for solar energy applications referred to as the Weather Research and Forecasting-Solar ensemble prediction system (WRF-Solar EPS). WRF-Solar EPS has been developed by introducing stochastic perturbations into the most relevant physical variables for solar irradiance predictions. In this study, we comprehensively discuss the impact of the stochastic perturbations of WRF-Solar EPS on solar irradiance forecasting compared to a deterministic WRF-Solar prediction (WRF-Solar DET), a stochastic ensemble using the stochastic kinetic energy backscatter scheme (SKEBS), and a WRF-Solar multi-physics ensemble (WRF-Solar PHYS). The performances of the four forecasts are evaluated using irradiance retrievals from the National Solar Radiation Database (NSRDB) over the contiguous United States. We focus on the predictability of the day-ahead solar irradiance forecasts during the year of 2018. The results show that the ensemble forecasts improve the quality of the forecasts, compared to the deterministic prediction system, by accounting for the uncertainty derived by the ensemble members. However, the three ensemble systems are under-dispersive, producing unreliable and overconfident forecasts due to a lack of calibration. In particular, WRF-Solar EPS produces less optically thick clouds than the other forecasts, which explains the larger positive bias in WRF-Solar EPS (31.7 W/m 2 ) than in the other models (22.7–23.6 W/m 2 ). This study confirms that the WRF-Solar EPS reduced the forecast error by 7.5% in terms of the mean absolute error (MAE) compared to WRF-Solar DET, and provides in-depth comparisons of forecast abilities with the conventional scientific probabilistic approaches (i.e., SKEBS and a multi-physics ensemble). Guidelines for improving the performance of WRF-Solar EPS in the future are provided.
In the Columbia and Snake River basins, several species of Pacific salmon were listed under the Endangered Species Act of 1973 due to significant declines of fish population. Dam operators and design engineers are thus faced with the task of making those hydroelectric facilities more ecologically friendly through changes in hydro-turbine design and operation. Public Utility District No. 2 of Grant County, Washington, applied for re-licensing from the U.S. Federal Energy Regulatory Commission to replace the 10 turbines at Wanapum Dam with advanced hydropower turbines that were designed to increase power generation and improve fish passage conditions. We applied both deterministic and stochastic blade-strike models to the newly installed turbine and an existing turbine. Modeled probabilities were compared to the results of a large-scale live fish survival study and a sensor fish study under the same operational parameters. Overall, injury rates predicted by the deterministic model were higher than experimental rates of injury while those predicted by the stochastic model were in close agreement with experiment results. Fish orientation at the time of entry into the plane of the leading edges of the turbine runner blades was an important factor contributing to uncertainty in modeled results. The advanced design turbine had slightly higher modeled injury rates than the existing turbine design; however, there was no statistical evidence that suggested significant differences in blade-strike injuries between the two turbines and the hypothesis that direct fish survival rate through the advanced hydropower turbine is equal or better than that through the conventional turbine could not be rejected.
Ecological mechanisms influence relationships among microbial communities, which in turn impact biogeochemistry. In particular, microbial communities are assembled by deterministic (e.g., selection) and stochastic (e.g., dispersal) processes, and the relative balance of these two process types is hypothesized to alter the influence of microbial communities over biogeochemical function. We used an ecological simulation model to evaluate this hypothesis, defining biogeochemical function generically to represent any biogeochemical reaction of interest. We assembled receiving communities under different levels of dispersal from a source community that was assembled purely by selection. The dispersal scenarios ranged from no dispersal (i.e., selection-only) to dispersal rates high enough to overwhelm selection (i.e., homogenizing dispersal). We used an aggregate measure of community fitness to infer a given community’s biogeochemical function relative to other communities. We also used ecological null models to further link the relative influence of deterministic assembly to function. We found that increasing rates of dispersal decrease biogeochemical function by increasing the proportion of maladapted taxa in a local community. Niche breadth was also a key determinant of biogeochemical function, suggesting a tradeoff between the function of generalist and specialist species. Finally, we show that microbial assembly processes exert greater influence over biogeochemical function when there is variation in the relative contributions of dispersal and selection among communities. Taken together, our results highlight the influence of spatial processes on biogeochemical function and indicate the need to account for such effects in models that aim to predict biogeochemical function under future environmental scenarios.
The relations between observable stellar parameters are usually assumed to be deterministic. That is, given an infinitely precise measurement of independent variable, “x”, and some model, the value of dependent variable, “y” can be known exactly. In practice this assumption is rarely valid and intrinsic stochasticity means that two stars with exactly the same “x” will have slightly different “y”s. The relation between short-timescale brightness fluctuations (flicker) of stars and both surface gravity and stellar density are two such stochastic relations that have until now been treated as deterministic ones. We recalibrate these relations in a probabilistic framework, using hierarchical Bayesian modeling to constrain the instrinsic scatter in the relations. We find evidence for additional scatter in the relationship, that cannot be accounted for by the observational uncertainties alone. The scatter in surface gravity and stellar density does not depend on flicker, suggesting that using flicker as a proxy for logg and ρ {sub ⋆} is equally valid for dwarf and giant stars, despite the fact that the observational uncertainties tend to be larger for dwarfs.
Dilute combustion with exhaust gas recirculation (EGR) in spark-ignition engines presents a cost-effective method for achieving higher levels of engine efficiency. At high levels of EGR, however, cycle-to-cycle variability (CCV) of the combustion process is exacerbated by sporadic occurrences of misfires and partial burns. Previous studies have shown that temporal deterministic patterns emerge at such conditions and certain combustion cycles have a significant influence over future events. Due to the complexity of the combustion process and the nature of CCV, harnessing all the deterministic information for control purposes has remained challenging even with physics based 0-D, 1-D, and high-fidelity computational fluid dynamics (CFD) models. In this study, we present a data-driven approach to optimize the combustion process by controlling CCV adjusting the cycle-to-cycle fuel injection quantity. Readily available data from in-cylinder pressure was used to train a spiking neural network (SNN) which learns the optimal way to manage fuel injection in order to reduce CCV while maintaining acceptable levels of fuel consumption. SNNs are particularly well suited for powertrain control applications due to their ability to be deployed on FPGA-based neuromorphic hardware which are small, inexpensive, and have a low power demand. The high-performance computing (HPC) resources of Oak Ridge National Laboratory were used to run an evolutionary-based training approach for choosing the best SNN configuration that minimizes the size of the network while achieving the desired goal. The neuromorphic hardware with the optimized SNN deployed was connected to the rapid prototyping engine control system for real-time control implementation and tested on a single cylinder version of a GM LNF 4-cylinder engine. The results show a significant reduction of CCV with a small percentage of additional fuel used to stabilize the charge.
Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.
Abstract. Conceptual frameworks linking microbial community membership, properties, and processes with the environment and emergent function have been proposed but remain untested. Here we refine and test a recent conceptual framework using hyporheic zone sediments exposed to wetting–drying transitions. Our refined framework includes relationships between cumulative properties of a microbial community (e.g., microbial membership, community assembly properties, and biogeochemical rates), environmental features (e.g., organic matter thermodynamics), and emergent ecosystem function. Our primary aim was to evaluate the hypothesized relationships that comprise the conceptual framework and contrast outcomes from the whole and putatively active bacterial and archaeal communities. Throughout the system we found threshold-like responses to the duration of desiccation. Membership of the putatively active community – but not the whole bacterial and archaeal community – responded due to enhanced deterministic selection (an emergent community property). Concurrently, the thermodynamic properties of organic matter (OM) became less favorable for oxidation (an environmental component), and respiration decreased (a microbial process). While these responses were step functions of desiccation, we found that in deterministically assembled active communities, respiration was lower and thermodynamic properties of OM were less favorable. Placing the results in context of our conceptual framework points to previously unrecognized internal feedbacks that are initiated by disturbance and mediated by thermodynamics and that cause the impacts of disturbance to be dependent on the history of disturbance.
Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.
Quantitative characterization of soil organic carbon (OC) content is essential due to its significant impacts on surface–subsurface hydrological–thermal processes and microbial decomposition of OC, which both in turn are important for predicting carbon–climate feedbacks. While such quantification is particularly important in the vulnerable organic-rich Arctic region, it is challenging to achieve due to the general limitations of conventional core sampling and analysis methods, and to the extremely dynamic nature of hydrological–thermal processes associated with annual freeze–thaw events. In this study, we develop and test an inversion scheme that can flexibly use single or multiple datasets – including soil liquid water content, temperature and electrical resistivity tomography (ERT) data – to estimate the vertical distribution of OC content. Our approach relies on the fact that OC content strongly influences soil hydrological–thermal parameters and, therefore, indirectly controls the spatiotemporal dynamics of soil liquid water content, temperature and their correlated electrical resistivity. We employ the Community Land Model to simulate nonisothermal surface–subsurface hydrological dynamics from the bedrock to the top of canopy, with consideration of land surface processes (e.g., solar radiation balance, evapotranspiration, snow accumulation and melting) and ice–liquid water phase transitions. For inversion, we combine a deterministic and an adaptive Markov chain Monte Carlo (MCMC) optimization algorithm to estimate a posteriori distributions of desired model parameters. For hydrological–thermal-to-geophysical variable transformation, the simulated subsurface temperature, liquid water content and ice content are explicitly linked to soil electrical resistivity via petrophysical and geophysical models. We validate the developed scheme using different numerical experiments and evaluate the influence of measurement errors and benefit of joint inversion on the estimation of OC and other parameters. We also quantify the propagation of uncertainty from the estimated parameters to prediction of hydrological–thermal responses. We find that, compared to inversion of single dataset (temperature, liquid water content or apparent resistivity), joint inversion of these datasets significantly reduces parameter uncertainty. We find that the joint inversion approach is able to estimate OC and sand content within the shallow active layer (top 0.3 m of soil) with high reliability. Due to the small variations of temperature and moisture within the shallow permafrost (here at about 0.6 m depth), the approach is unable to estimate OC with confidence. However, if the soil porosity is functionally related to the OC and mineral content, which is often observed in organic-rich Arctic soil, the uncertainty of OC estimate at this depth remarkably decreases. Our study documents the value of the new surface–subsurface, deterministic–stochastic inversion approach, as well as the benefit of including multiple types of data to estimate OC and associated hydrological–thermal dynamics.
Optimized static and dynamic linear system compensator designs with a priori fixed dimensions are discussed. Expressions are derived to determine the parameters of these designs, including the fixed-gain output feedback, the time-varying output feedback, and time-invariant compensation. A relevant deterministic off-line optimization problem is formulated, and the properties of a solution to the problem are analyzed. A unified design philosophy for limited state-variable control problems of essentially deterministic systems is developed. The philosophy incorporates design constraints, a careful selection of performance criteria, and interpretation of the results.
A unified design philosophy is presented for limited state feedback control problems with parameter uncertainty for both deterministic and stochastic problems. Two approaches are considered: linear compensator for the deterministic problem with parameter uncertainty, and for the single input-single output system with parameter uncertainty, a model on order equal to that of the system less the number of zeroes. The limitations of these approaches are discussed along with suggestions for further research.
A generalized problem of stochastic control is discussed in which multiple controllers with different data bases are present. The vehicle for the investigation is the finite state, finite memory (FSFM) stochastic control problem. Optimality conditions are obtained by deriving an equivalent deterministic optimal control problem. A FSFM minimum principle is obtained via the equivalent deterministic problem. The minimum principle suggests the development of a numerical optimization algorithm, the min-H algorithm. The relationship between the sufficiency of the minimum principle and the informational properties of the problem are investigated. A problem of hypothesis testing with 1-bit memory is investigated to illustrate the application of control theoretic techniques to information processing problems.