Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic analysis”

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 469 records · Page 26

Gaussian Process for Flight Delay Prediction: Learning a Stochastic Process

This paper presents a machine-learning approach to predict flight delays. Whereas neural networks are extensively studied for predictive capabilities, they involve non-intuitive design and extensive analysis, particularly in training and optimization processes. Instead, the proposed framework employs Gaussian Processes as a supervised learning technique for flight delay prediction. This data-driven approach trains the model using prior information, specifically the mean and covariance tied to existing data. The proposed Gaussian Process Regression (GPR) model employs the day of flight as a pivotal feature for delay forecasting. We analyze flights from various routes and gauge the accuracy of the presented learning technique by comparing the predicted delays with the actual ones. Given the inherent challenges in precisely forecasting delays, we predict the delays with a 95 % confidence interval. Also, an error propagation analysis in the prediction horizon is carried out to determine the optimal time frame for prediction. The proposed method for flight delay prediction is important as airlines can strategize flight operations and issue timely advisories.

stochastic↗

Modeling laser velocimeter signals as triply stochastic Poisson processes

Previous models of laser Doppler velocimeter (LDV) systems have not adequately described dual-scatter signals in a manner useful for analysis and simulation of low-level photon-limited signals. At low photon rates, an LDV signal at the output of a photomultiplier tube is a compound nonhomogeneous filtered Poisson process, whose intensity function is another (slower) Poisson process with the nonstationary rate and frequency parameters controlled by a random flow (slowest) process. In the present paper, generalized Poisson shot noise models are developed for low-level LDV signals. Theoretical results useful in detection error analysis and simulation are presented, along with measurements of burst amplitude statistics. Computer generated simulations illustrate the difference between Gaussian and Poisson models of low-level signals.

Mayo, W. T., Jr.↗

Analysis and Monte Carlo simulation of near-terminal aircraft flight paths

The flight paths of arriving and departing aircraft at an airport are stochastically represented. Radar data of the aircraft movements are used to decompose the flight paths into linear and curvilinear segments. Variables which describe the segments are derived, and the best fitting probability distributions of the variables, based on a sample of flight paths, are found. Conversely, given information on the probability distribution of the variables, generation of a random sample of flight paths in a Monte Carlo simulation is discussed. Actual flight paths at Dulles International Airport are analyzed and simulated.

Schiess, J. R.↗

Simulation of probabilistic wind loads and building analysis

Probabilistic wind loads likely to occur on a structure during its design life are predicted. Described here is a suitable multifactor interactive equation (MFIE) model and its use in the Composite Load Spectra (CLS) computer program to simulate the wind pressure cumulative distribution functions on four sides of a building. The simulated probabilistic wind pressure load was applied to a building frame, and cumulative distribution functions of sway displacements and reliability against overturning were obtained using NESSUS (Numerical Evaluation of Stochastic Structure Under Stress), a stochastic finite element computer code. The geometry of the building and the properties of building members were also considered as random in the NESSUS analysis. The uncertainties of wind pressure, building geometry, and member section property were qualified in terms of their respective sensitivities on the structural response.

Shah, Ashwin R.↗

Nth-plant supply: corn stover supplies and costs in a fleet of biorefineries

Feedstock cost and cost variability is expected to increase with the number of biorefineries. To quantify this effect, this spatial-economic analysis simulates feedstock cost and cost variability of an industry based on corn stover as a function of the number of biorefineries. Results are reported for nine scenarios (a base case and sensitivity analysis of four variables – harvest efficiency, sustainability constraints, opportunity cost, and corn grain yield) under deterministic and stochastic simulations, assuming biorefineries using 658 000 Mg (725 000 tons) year- 1 of corn stover in 2019. The resulting supply curves are highly elastic (i.e. little change in cost) for the first 50 of the 121 biorefineries, with price increases in subsequent biorefineries depending on scenario. In the base-case deterministic scenario, weighted-average stover costs are $66 Mg -1 ($60 ton- 1 ), $69 Mg -1 ($62 ton -1 ), and $156 Mg -1 ($142 ton -1 ), at the first, 60th, and 121st biorefineries, respectively. The stochastic simulations, subject to observed 30-year corn yield variability, follow a similar pattern, with price distributions that vary by scenario. The base-case stochastic simulations illustrate minimal cost variability for the first 60 biorefineries, but rapid increases in cost variability in the second half of potential biorefineries, with similar patterns observed in the other scenarios. Of the four variables explored, price was most sensitive to harvest efficiency, followed by sustainability constraints, corn yield, and opportunity cost. Results suggest that, under conventional logistics, about half of the US corn stover resource is reliably available with minimum cost increase and variability. Interactive visualization is available at https://doi.org/10.11578/1828779.

09 BIOMASS FUELS↗

A Machine Learning Assisted Development of a Model for the Populations of Convective and Stratiform Clouds

Abstract Traditional parameterizations of the interaction between convection and the environment have relied on an assumption that the slowly varying large‐scale environment is in statistical equilibrium with a large number of small and short‐lived convective clouds. They fail to capture nonequilibrium transitions such as the diurnal cycle and the formation of mesoscale convective systems as well as observed precipitation statistics and extremes. Informed by analysis of radar observations, cloud‐permitting model simulation, theory, and machine learning, this work presents a new stochastic cloud population dynamics model for characterizing the interactions between convective and stratiform clouds, with the goal of informing the representation of these interactions in global climate models. Fifteen wet seasons of precipitating cloud observations by a C‐band radar at Darwin, Australia are fed into a machine learning algorithm to obtain transition functions that close a set of coupled equations relating large‐scale forcing, mass flux, the convective cell size distribution, and the stratiform area. Under realistic large‐scale forcing, the derived transition functions show that, on the one hand, interactions with stratiform clouds act to dampen the variability in the size and number of convective cells and therefore in the convective mass flux. On the other, for a given convective area fraction, a larger number of smaller cells is more favorable for the growth of stratiform area than a smaller number of larger cells. The combination of these two factors gives rise to solutions with a few convective cells embedded in a large stratiform area, reminiscent of mesoscale convective systems.

54 ENVIRONMENTAL SCIENCES↗

Incorporating long-range dependence and fractal features in turbulence spectra

We introduce an advanced turbulence spectrum model developed from mathematical foundations from a covariance function class and empirically validated using extensive field data. This model captures the complex dynamics of long-range dependence, and fractal characteristics prevalent in riverine and atmospheric boundary layer (ABL) flows that are ignored by classical spectrum models, such as IEC (International Electrotechnical Commission) von Kármán and Kaimal model. The model delineates scaling behaviors across distinct frequency bands and offers substantial flexibility through five well-defined parameters each characterizing a distinct physical aspect of the velocity time series. A detailed procedure for obtaining each parameter from time series data is outlined. The comprehensive validations with field data from tidal currents and ABL flows substantiate the model’s fidelity in accurately replicating observed phenomena. This validation establishes the reliability of the proposed model and, when incorporated into stochastic full-field simulators such as TurbSim, demonstrates its potential to advance the predictive modeling and analysis of turbulent flows in environmental science and engineering contexts.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗

NO PLIF flow visualization and time-resolved temperature distributions in laser induced breakdown plumes

NO planar laser induced fluorescence (PLIF) is used to obtain images of laser-induced breakdown plasma plumes in NO-seeded nitrogen and dry air at near atmospheric pressure. Single-shot PLIF-images show that the plume development 5–50 μs after the breakdown pulse is fairly reproducible shot-to-shot, although the plume becomes increasingly stochastic on longer timescales, 100–500 μs. The stochastic behavior of the plume is quantified using probability distributions of the loci of the plume boundary. Analysis of the single-shot images indicates that the mixing of the plume with ambient gas on sub-ms time scale is insignificant. The induced flow velocity in the plume is fairly low, up to 30 m s –1 , suggesting that laser breakdowns are ineffective for mixing enhancement in high speed flows. The ensemble-averaged PLIF images indicate the evolution of the plume from an initially elongated shape to near-spherical to toroidal shape, with a subsequent radial expansion and formation of an axial jet in the center. Temperature distributions in the plume in air are obtained from the NO PLIF images, using two rotational transitions in the NO(X, v' = 0 → A, v'' = 0) band, J'' = 6.5 and 12.5 of the Q R 12 + Q 2 branch. The results indicate that the temperature in the plume remains high, above 1000 K, for approximately 100 μs, after which it decays gradually, to below 500 K at 500 μs. The residual NO fraction in the plume is ~0.1%, indicating that repetitive laser-assisted ignition may result in significant NO-generation. Furthermore, these measured temperature and velocity distributions can be used for detailed validation of kinetic models of laser-induced breakdown and assessment of their predictive capability.

42 ENGINEERING↗

Optimality of Gradient-MUSIC for Spectral Estimation

We introduce the Gradient-MUSIC algorithm for estimating the unknown frequencies and amplitudes of a nonharmonic signal from noisy time samples. While the classical MUSIC algorithm performs a computationally expensive search over a fine grid, Gradient-MUSIC is significantly more efficient and eliminates the need for discretization over a fine grid by using optimization techniques. It coarsely scans the 1D landscape to find initialization simultaneously for all frequencies followed by parallelizable local refinement via gradient descent. We also analyze its performance when the noise level is sufficiently small and the signal frequencies are separated by at least 8π/m, where π/m is the standard resolution of this problem. Even though the 1D landscape is nonconvex, we prove a global convergence result for Gradient-MUSIC: coarse scanning provably finds suitable initialization and gradient descent converges at a linear rate. In addition to convergence results, we also upper bound the error between the true signal frequencies and amplitudes with those found by Gradient-MUSIC. For example, if the noise has $\ell^\infty$ norm at most ϵ, then the frequencies and amplitudes are recovered up to error at most Cϵ/m and Cϵ respectively, which are minimax optimal in m and ϵ. Our theory can also handle stochastic noise with performance guarantees under nonstationary independent Gaussian noise. Our main approach is a comprehensive geometric analysis of the landscape, a perspective that has not been explored before.

97 MATHEMATICS AND COMPUTING↗

Energy-balance climate models

An introductory survey of the global energy balance climate models is presented with an emphasis on analytical results. A sequence of increasingly complicated models involving ice cap and radiative feedback processes are solved and the solutions and parameter sensitivities are studied. The model parameterizations are examined critically in light of many current uncertainties. A simple seasonal model is used to study the effects of changes in orbital elements on the temperature field. A linear stability theorem and a complete nonlinear stability analysis for the models are developed. Analytical solutions are also obtained for the linearized models driven by stochastic forcing elements. In this context the relation between natural fluctuation statistics and climate sensitivity is stressed.

North, G. R.↗

Energy balance climate models

An introductory survey of the global energy balance climate models is presented with an emphasis on analytical results. A sequence of increasingly complicated models involving ice cap and radiative feedback processes are solved, and the solutions and parameter sensitivities are studied. The model parameterizations are examined critically in light of many current uncertainties. A simple seasonal model is used to study the effects of changes in orbital elements on the temperature field. A linear stability theorem and a complete nonlinear stability analysis for the models are developed. Analytical solutions are also obtained for the linearized models driven by stochastic forcing elements. In this context the relation between natural fluctuation statistics and climate sensitivity is stressed.

North, G. R.↗

Probabilistic analysis of a materially nonlinear structure

A probabilistic finite element program is used to perform probabilistic analysis of a materially nonlinear structure. The program used in this study is NESSUS (Numerical Evaluation of Stochastic Structure Under Stress), under development at Southwest Research Institute. The cumulative distribution function (CDF) of the radial stress of a thick-walled cylinder under internal pressure is computed and compared with the analytical solution. In addition, sensitivity factors showing the relative importance of the input random variables are calculated. Significant plasticity is present in this problem and has a pronounced effect on the probabilistic results. The random input variables are the material yield stress and internal pressure with Weibull and normal distributions, respectively. The results verify the ability of NESSUS to compute the CDF and sensitivity factors of a materially nonlinear structure. In addition, the ability of the Advanced Mean Value (AMV) procedure to assess the probabilistic behavior of structures which exhibit a highly nonlinear response is shown. Thus, the AMV procedure can be applied with confidence to other structures which exhibit nonlinear behavior.

Millwater, H. R.↗

Generalized covariance analysis for partially autonomous deep space missions

A new covariance analysis method is presented that is suitable for the evaluation of multiple impulsive controllers acting on some stochastic process x. The method accommodates batch and sequential estimators with equal ease and accounts for time-delay effects in a natural manner. The formalism is developed in terms of a generalized state vector that is formed from the system state vector x, augmented by various fixed epoch estimates, and a data vector formed from discrete time observations of the system. Recursions are developed for time transition, measurement incorporation, and impulsive control updating of the generalized covariance matrix. Means of limiting the dimensional growth of the generalized state vector via the processes of estimator epoch adjustment and measurement vector deflation are described and the application of numerically stable matrix factorization methods to the generalized covariance recursions is outlined. The method is applied to the Magellan spacecraft to demonstrate the capability of ground-based optimal estimation and control of gyro/star scanner misalignment.

Boone, Jack N.↗

Assessing the Effect of Spaceflight on the Propensity for Astronauts to Develop Disk Herniation

BACKGROUND: A previous study [1] reported that the instantaneous risk of developing a Herniated Nucleus Pulposus (HNP) was higher in astronauts who had flown at least one mission, as compared with those in the corps who had not yet flown. However, the study only analyzed time to HNP after the first mission (if any) and did not account for the possible effects of multiple missions. While many HNP's occurred well into astronauts' careers or in some cases years after retirement, the higher incidence of HNPs relatively soon after completion of space missions appears to indicate that spaceflight may lead to an increased risk of HNP. The purpose of this study was to support the Human System Risk Board assessment of back pain, evaluate the risk of injury due to dynamic loads, and update the previous dataset which contained events up to December 31, 2006. METHODS: Data was queried from the electronic medical record and provided by the Lifetime Surveillance of Astronaut Health. The data included all 330 United States astronauts from 1959 through February 2014. Cases were confirmed by Magnetic Resonance Imaging, Computerized Tomography, Myelography, operative findings, or through clinical confirmation with a neurologist or neurosurgeon. In this analysis, astronauts who had an HNP at selection into the corps or had an HNP diagnosis prior to their first flight were excluded. The statistical challenges in using the available data to separate effects of spaceflight from those associated with general astronaut training and lifestyle on propensity to develop HNPs are many. The primary outcome is reported date of first HNP (if any), which at best is only an approximation to the actual time of occurrence. To properly analyze this data with a survival analysis model, one must also know the "exposure" time - i.e. how long each astronaut has been at risk for developing an HNP. If an HNP is reported soon after a mission, is it mission caused or general? If the former, exposure time should be counted from the time of landing (assuming the risk of HNP occurring during a mission is zero). If the latter, exposure time should be counted from the time of selection; however we can't directly know which one to use. In our analysis we take both of these possibilities into account with a competing risks model, wherein two distinct stochastic processes are going on: TG = time to HNP (general) and TS = time to HNP (spaceflight). Under this type of model, whichever of these occurs first is what we observe; in other words we don't observe TG or TS, only min(TG, TS). Here, we parameterized the model in terms of separate Weibull hazard functions for each process and estimated all parameters using maximum likelihood. In addition, we allowed for a "cured fraction" - i.e. the possibility that some astronauts may never develop an HNP. RESULTS: Results will include a depiction of the competing hazard functions as well as a probability curve for the relative likelihood that an HNP reported at a given time after a mission is actually mission caused. Other factors, such as dwell time in microgravity and vehicle landing environment will be explored. An overall assessment as to whether spaceflight truly exacerbates HNP risk will be made.

Feiveson, A. H.↗

Multifidelity, Multidisciplinary Design Under Uncertainty with Non-Intrusive Polynomial Chaos

The primary objective of this work is to develop an approach for multifidelity uncertainty quantification and to lay the framework for future design under uncertainty efforts. In this study, multifidelity is used to describe both the fidelity of the modeling of the physical systems, as well as the difference in the uncertainty in each of the models. For computational efficiency, a multifidelity surrogate modeling approach based on non-intrusive polynomial chaos using the point-collocation technique is developed for the treatment of both multifidelity modeling and multifidelity uncertainty modeling. Two stochastic model problems are used to demonstrate the developed methodologies: a transonic airfoil model and multidisciplinary aircraft analysis model. The results of both showed the multifidelity modeling approach was able to predict the output uncertainty predicted by the high-fidelity model as a significant reduction in computational cost.

West, Thomas K., IV↗

Modeling the Uncertainties of Solar System Ephemerides for Robust Gravitational-wave Searches with Pular-timing Arrays

The regularity of pulsar emissions becomes apparent once we reference the pulses’ times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earthʼs position with respect to the solar system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar-timing arrays. Indeed, recent array data sets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article, we describe the motivation, construction, and application of a physical model of solar system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiterʼs orbital elements) most relevant to gravitational-wave searches with pulsar-timing arrays. This model, BAYESEPHEM, was used to derive ephemeris-robust results in NANOGravʼs 11 yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11 yr data set and that this degeneracy will vanish with improved ephemerides and with pulsar-timing data sets that extend well beyond a single Jovian orbital period.

M Vallisneri↗

Stochastic Adaptive Droop Control in Frequency Regulation of Power Systems With Intermittent Generators

Modern power systems (MPSs), including microgrids (MGs), are increasingly incorporating multiple renewable energy sources (RESs) such as wind and solar power, as well as battery storage and controllable loads. While environmentally beneficial, these sources pose challenges for control and management due to their intermittent and stochastic nature, especially in maintaining frequency stability with multiple interconnected generators of varying capacities. Traditional droop control methods are effective in systems with generators that are dispatchable and have fixed generation capacities, but they fall short when applied to systems with RESs, where generation capacities are dynamic and affected by unpredictable environmental conditions. To address these challenges, this paper introduces a novel stochastic adaptive droop control (SADC) method for load frequency control (LFC). The proposed method adapts droop coefficients in real time, based on the measured stochastic data of power generation capacities, enabling more effective frequency regulation in systems with variable and intermittent power generation. Unlike traditional adaptive control methods, which assume constant or slowly-varying system parameters, this approach accounts for stochastic processes by modeling them as Markov chains, enabling robust performance under highly dynamic and unpredictable conditions. The key contributions of this work include the development of real-time droop coefficient adaptation algorithms, derivation of their stability and convergence properties, and the demonstration of the advantages of the method through simulations. Case studies highlight the improved performance of frequency regulation, particularly in addressing the impact of stochastic weather conditions and the benefits of reducing dependence on battery reserves in dealing with intermittency of RESs. Finally, this paper provides a comprehensive analysis of the theoretical foundations of the method, as well as practical implementation insights for future power systems with high penetration of RESs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Failure Analysis–Informed Risk Assessment Framework for Geological Carbon Storage Using Numerical Simulation and Machine Learning

Geological carbon storage (GCS) is recognized as a critical technology for achieving large-scale reductions in anthropogenic carbon dioxide (CO 2 ) emissions. Ensuring long-term containment and safety requires robust risk assessment frameworks that account for geological uncertainty and identify potential failure scenarios. Among various indicators, the area of review (AoR) serves as a key metric for evaluating storage performance, regulatory compliance, and monitoring design, as it delineates the spatial extent impacted by pressure buildup and plume migration. However, conventional AoR-based risk assessments typically perturb parameters within narrow uncertainty bounds, potentially overlooking rare but high-impact events arising from extreme geological conditions. In this study, we present a failure analysis–informed risk assessment framework for large-scale GCS projects to improve site prescreening and monitoring design. A suite of 300 numerical simulations was generated using stochastic geological models that vary five key parameters: net-to-gross ratio, anisotropy azimuth, porosity multiplier, permeability multiplier, and vertical-to-horizontal permeability ratio. Among these, 200 realizations represent normal geological uncertainty, while 100 additional cases explore extreme yet plausible conditions for failure-case analysis. The AoR was simulated and computed from pressure and CO 2 saturation fields, where the baseline AoR boundary, representing the extent predicted under typical geological uncertainty, was defined as the union of 200 normal-range simulations, and failure was identified when extreme-range cases exceeded this baseline. Results show that incorporating broader parameter uncertainty produces significantly larger AoR extents, underscoring the potential underestimation of risk under conventional uncertainty ranges. Furthermore, spatial probability maps derived from failure-induced AoR exceedance identify regions requiring enhanced monitoring attention. Various machine learning (ML)–based classifiers were developed to predict failure occurrence from geological parameters, with the random forest model achieving the highest performance (F1-score of 0.986). Consistent findings from correlation coefficient, feature importance, and Sobol sensitivity analyses reveal that low net-to-gross ratios and permeability multipliers are the dominant risk drivers, reflecting reduced reservoir connectivity and limited pressure dissipation. Altogether, these results provide a novel framework for risk-informed site prescreening and monitoring design that explicitly considers rare but high-impact geological scenarios in GCS projects.

25 ENERGY STORAGE↗