Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic process”

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,135 records · Page 63

Sensitivity analysis, surrogate modeling, and optimization of pebble-bed reactors considering normal and accident conditions

This research provides a valuable tool that streamlines the optimization process while significantly increasing its accuracy. This study creates a robust framework for reactor design optimization by incorporating comprehensive modeling using the Comprehensive Reactor Analysis Bundle, or BlueCRAB, within the Multiphysics Object-Oriented Simulation Environment (MOOSE). BlueCRAB is the United States Nuclear Regulatory Commission's code suite for non-light water reactor analysis and includes the Griffin, Pronghorn, and Bison applications. This not only improves the efficiency of the optimization process but also enhances the reliability of the results. Such a tool is essential for advancing the state-of-the-art in pebble-bed reactor technology and is critical for achieving the goals of Generation IV reactors, which aim for safe, sustainable, and economically viable nuclear energy solutions. This work presents and applies this workflow on pebble-bed reactors while considering both normal and off-normal conditions. A representative gas-cooled pebble-bed reactor at equilibrium core conditions serves as the nominal design specification for normal operation and is based on previous research. The depressurized loss-of-forced-cooling accident is deployed for off-normal conditions in this work. After defining design-related parameters and quantities of interest regarding reactor safety and performance, this multiphysics model is sampled using the MOOSE stochastic tools module. The result is a comprehensive dataset of configurations, enabling sensitivity analysis and the generation of surrogate models. Subsequently, the dataset and surrogate models are employed in two optimization studies aimed at maximizing fuel utilization and economic profit while adhering to safety and operational constraints. Performing the optimization process with fuel utilization as the metric leads to an improvement of approximately 10%, compared to engineering-judgment-based nominal conditions. The optimization on economic profit leads to an estimated increase of ~300 million USD over the lifetime of the reactor.

97 MATHEMATICS AND COMPUTING↗

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗

Quantifying Load Uncertainty Using Real Smart Meter Data

As we get closer to customers in distribution systems, load stochasticity increases. In the past, due to lack of real-time data, the comprehensive knowledge of load behavior was limited, and simplistic assumptions had to be made for distribution system modeling and analysis, especially in the processes of network design and expansion. With the deployment of Advanced Metering Infrastructure (AMI), ample real-time smart meter data has become available to utilities. In this paper, using real hourly smart meter data, we have quantified load uncertainty in terms of average, maximum and maximum noncoincident demands on a daily basis, as well as load factor and diversity factor. These uncertainty metrics are examined for individual residential, commercial and industrial customers, as well as distribution transformers serving residential customers. This paper provides a benchmark on load uncertainty quantification for practicing engineers and researchers.

Bu, Fankun↗

Machine Learning Assisted HPC Workload Trace Generation for Leadership Scale Storage Systems

Monitoring and analyzing a wide range of I/O activities in an HPC cluster is important in maintaining mission-critical performance in a large-scale, multi-user, parallel storage system. Center-wide I/O traces can provide high-level information and fine-grained activities per application or per user running in the system. Studying such large-scale traces can provide helpful insights into the system. It can be used to develop predictive methods for making predictive decisions, adjusting scheduling policies, or providing decisions for the design of next-generation systems. However, sharing real-world I/O traces to expedite such research efforts leaves a few concerns; i) the cost of sharing the large traces is expensive due to this large size, and ii) privacy concern is an issue.We address such issues by building an end-to-end machine learn- ing (ML) workflow that can generate I/O traces for large-scale HPC applications. We leverage ML based feature selection and gener- ative models for I/O trace generation. The generative models are trained on I/O traces collected by the darshan I/O characterization tool over a period of one year. We present a two-step generation process consisting of two deep-learning models, called the feature generator and the trace generator. The combination of two-step generative models provides robustness by reducing the bias of the model and accounting for the stochastic nature of the I/O traces across different runs of an application. We evaluate the performance of the generative models and show that the two-step model can generate time-series I/O traces with less than 20% root mean square error.

Paul, Arnab↗

Strong Red Noise Ocean Forcing on Atlantic Multidecadal Variability Assessed from Surface Heat Flux: Theory and Application

The role of ocean forcing on Atlantic multidecadal variability (AMV) is assessed from the (downward) heat flux–SST relation in the framework of a new stochastic climate theory forced by red noise ocean forcing. Previous studies suggested that atmospheric forcing drives SST variability from monthly to interannual time scales, with a positive heat flux–SST correlation, while heat flux induced by ocean processes can drive SST variability at decadal and longer time scales, with a negative heat flux–SST correlation. Here, first, we develop a theory to show how the sign of heat flux–SST correlation is affected by atmospheric and oceanic forcing with time scale. In particular, a red noise ocean forcing is necessary for the sign reversal of heat flux–SST correlation. Furthermore, this sign reversal can be detected equivalently in three approaches: the low-pass correlation at lag zero, the unfiltered correlation at long (heat flux) lead, and the real part of the heat flux–SST coherence. Second, we develop a new scheme in combination with the theory to assess the magnitude and time scale of the red noise ocean forcing for AMV in the GFDL SPEAR model (Seamless System for Prediction and Earth System Research) and observations. In both the model and observations, the ocean forcing on AMV is in general comparable with the atmospheric forcing, with a 90% probability greater than the atmospheric forcing in observations. In contrast to the white noise atmospheric forcing, the ocean forcing has a persistence time comparable or longer than a year, much longer than the SST persistence of ~3 months. This slow ocean forcing is associated implicitly with slow subsurface ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗

Energetic Electron Transport in Magnetized Plasma with Magnetic Islands (Final Technical Report)

The main scientific goal of this project was to investigate the interactions between energetic electrons and magnetic islands in magnetized plasmas from lab to space. The proposed plan of addressed the following specific objectives: (1) Determine if energetic electrons are trapped by magnetic islands. Specifically, test the hypothesis that energetic electrons are trapped near island O-points due to stable orbits and de-confined near island X-points due to X-point tangles. (2) Determine how electron transport changes with island width and location. Specifically, test the hypothesis that electrons are accelerated by contracting magnetic islands through a Fermi acceleration process. (3) Assess how electron transport is affected by island dynamics, including island rotation, bifurcation, and/or island overlap. Specifically, test the hypothesis that electrons are deconfined and possibly accelerated during island rotations, bifurcations, and overlap due to stochastization of the magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Persistent Elevated Soot Emissions Induced by Clustered Stochastic Preignition Events

Stochastic Preignition (SPI) is an abnormal combustion phenomenon that can occur in spark-ignition engines particularly under high-load operation. SPI is characterized by uncontrolled initiation of combustion prior to spark discharge, an abnormal combustion process that can lead to severe knock events and significant engine damage. SPI has been associated with fuel properties, lubricant composition, and engine design and operation. Here, in this work, a single-cylinder test engine with a dry-sump oil system was utilized to study the SPI response of E10 and E25 fuels with a range of Reid Vapor Pressure (RVP). An automated test procedure was employed, consisting of ten square-waved load profile segments, with each segment composed of 5 min of low-load operation followed by 25 min of sustained high-load operation. These tests were replicated across multiple days of testing including a lubricant triple flush between tests, and an online Fuel in Oil diagnostic measurement. Exhaust particulate emissions were continuously measured by an AVL microsoot sensor (MSS). Elevated particulate matter emissions were observed to occur concurrently with SPI events as blooms of soot. Particularly after clustered events (i.e., multiple SPI cycles occurring within 10 consecutive engine cycles), high soot emissions were observed to persist over several days of sequential operation despite daily lubricant changes, a complete warm-up procedure, and sustained low-load operation between test segments. This result implies that the particulate emissions trends may be dominated by deposit-based effects, where higher load operation is needed to alter deposition and formation processes. The observed soot blooms were also found to correspond to a reduction in the engine fueling and the fuel engine oil dilution rate despite the engine exhaust remaining at stoichiometric exhaust operation. These observations suggest that post-SPI events, pathways for lubricant migration and consumption into the combustion chamber may occur until these pathways are closed from deposit formation or ring dynamics during extended operation. These observed sooting propensity persisted with all fuels tests, but a linear correlation was observed between the summation of soot and particulate matter index (PMI) value for each fuel as well as SPI events, proving that PMI is a crucial fuel property for reducing SPI.

Splitter, Derek [Oak Ridge National Laboratory (OR↗

Navigation capability for an ion drive rendezvous with Halley's Comet

Navigation accuracies are presented for a 1985 rendezvous with Halley's Comet using an ion propulsion system. Individual error sources are examined to determine their relative contributions to final delivery errors. The sensitivity of delivery accuracy to stochastic thrust variations is demonstrated by considering a baseline thrust error model and a more benign alternative model. Also studied are the effects of increased comet ephemeris uncertainties and of operational time delays between the orbit determination process and the guidance implementation.

Thornton, C. L.↗

Discrete-time Markovian stochastic Petri nets

We revisit and extend the original definition of discrete-time stochastic Petri nets, by allowing the firing times to have a 'defective discrete phase distribution'. We show that this formalism still corresponds to an underlying discrete-time Markov chain. The structure of the state for this process describes both the marking of the Petri net and the phase of the firing time for each transition, resulting in a large state space. We then modify the well-known power method to perform a transient analysis even when the state space is infinite, subject to the condition that only a finite number of states can be reached in a finite amount of time. Since the memory requirements might still be excessive, we suggest a bounding technique based on truncation.

Ciardo, Gianfranco↗

A large-scale benchmarking of deterministic and stochastic derivative-free optimization algorithms

This presentation summarizes our work in the PrOMMiS project on benchmarking of data-driven optimization algorithms and their applications in self-driving laboratories. This work supports the broader project goal of accelerating the identification of promising separation methods and operating conditions for critical minerals separation processes. We present a systematic benchmarking study of 42 data-driven optimization algorithms on a broad collection of 502 test problems. The results identify BAM, GLCCLUSTER, and MULTIMIN as the most effective optimization solvers, with BAM showing the highest overall performance and solving more than 80% of the benchmark problems. The study also shows that no single solver consistently outperforms the others across all problem types, indicating that our future laboratory applications may benefit from using a small set of strong solvers rather than relying on a single method. The presentation also illustrates an in-silico chemical reactor case study showing that data-driven optimization methods can guide autonomous experimentation in a self-driving laboratory and identify optimal operating conditions within a small number of experiments. Overall, the results provide a basis for selecting efficient optimization methods and demonstrate the practical use of data-driven optimization in self-driving laboratory workflows.

36 MATERIALS SCIENCE↗

Machine Learning Guided Synthesis of Flash Graphene

Advances in nanoscience have enabled the synthesis of nanomaterials, such as graphene, from low-value or waste materials through flash Joule heating. Though this capability is promising, the complex and entangled variables that govern nanocrystal formation in the Joule heating process remain poorly understood. In this work, machine learning (ML) models are constructed to explore the factors that drive the transformation of amorphous carbon into graphene nanocrystals during flash Joule heating. An XGBoost regression model of crystallinity achieves an r 2 score of 0.8051 ± 0.054. Feature importance assays and decision trees extracted from these models reveal key considerations in the selection of starting materials and the role of stochastic current fluctuations in flash Joule heating synthesis. Furthermore, partial dependence analyses demonstrate the importance of charge and current density as predictors of crystallinity, implying a progression from reaction-limited to diffusion-limited kinetics as flash Joule heating parameters change. Finally, a practical application of the ML models is shown by using Bayesian meta-learning algorithms to automatically improve bulk crystallinity over many Joule heating reactions. Furthermore, these results illustrate the power of ML as a tool to analyze complex nanomanufacturing processes and enable the synthesis of 2D crystals with desirable properties by flash Joule heating.

01 COAL, LIGNITE, AND PEAT↗

Attack-Resilient Weighted $\ell_{1}$ Observer with Prior Pruning

Security related questions for Cyber Physical Systems (CPS) have attracted much research attention in searching for novel methods for attack-resilient control and/or estimation. Specifically, false data injection attacks (FDIAs) have been shown to be capable of bypassing bad data detection (BDD), while arbitrarily compromising the integrity of state estimators and robust controller even with very sparse measurements corruption. Moreover, based on the inherent sparsity of pragmatic attack signals, ℓ1 -minimization scheme has been used extensively to improve the design of attack-resilient estimators. For this, the theoretical maximum for the percentage of compromised nodes that can be accommodated has been shown to be 50%. In order to guarantee correct state recoveries for larger percentage of attacked nodes, researchers have begun to incorporate prior information into the underlying resilient observer design framework. For the most pragmatic cases, this prior information is often obtained through some data-driven machine learning process. Existing results have shown strong positive correlation between the tolerated attack percentages and the precision of the prior information. In this paper, we present a pruning method to improve the precision of the prior information, given corresponding stochastic uncertainty characteristics of the underlying machine learning model. Then a weighted ℓ1 -minimization is proposed based on the pruned prior. The theoretical and simulation results show that the pruning method significantly improves the observer performance for much larger attack percentages, even when moderately accurate machine learning model used.

Resilient observer, Cyber-physical systems, prunin↗

A theoretical and experimental study of turbulent nonevaporating sprays

Measurements and analysis limited to the dilute portions of turbulent nonevaporating sprays injected into a still air environment were completed. Mean and fluctuating velocities and Reynolds stress were measured in the continuous phase. Liquid phase measurements included liquid mass fluxes, drop sizes and drop size and velocity correlation. Initial conditions needed for model evaluation were measured at a location as close to the injector exit as possible. The test sprays showed significant effects of slip and turbulent dispersion of the discrete phase. The measurements were used to evaluate three typical models of these processes: (1) a locally homogenous flow (LHF) model, where slip between the phases were neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for drop motion. The LHF and DSF models did not provide very satisfactory predictions for the present measurements. In contrast, the SSF model performed reasonably well with no modifications in the prescription of eddy properties from its original calibration. Some effects of drops on turbulence properties were observed near the dense regions of the sprays.

Solomon, A. S. P.↗

A theoretical and experimental study of turbulent evaporating sprays

Measurements and analysis limited to the dilute portions of turbulent evaporating sprays, injected into a still air environment were completed. Mean and fluctuating velocities and Reynolds stress were measured in the continuous phase. Liquid phase measurements included liquid mass fluxes, drop sizes and drop size and velocity correlation. Initial conditions needed for model evaluation were measured at a location as close to the injector exit as possible. The test sprays showed significant effects of slip and turbulent dispersion of the discrete phase. The measurements were used to evaluate three typical models of these processes: (1) a locally homogeneous flow (LHF) model, where slip between the phases were neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for drop motion. For all three models, a k-epsilon model as used to find the properties of the continuous phase. The LHF and DSF models did not provide very satisfactory predictions for the present measurements. In contrast, the SSF model performed reasonably well--with no modifications in the prescription of eddy properties from its original calibration.

Solomon, A. S. P.↗

A coalescence/dispersion model for turbulent flame stability

A stochastic mixing model is employed that couples finite rate mixing with finite rate combustion chemistry to examine the effects of combustor inlet conditions on the ignition limits of lean premixed turbulent flames. The kinetics of the fuel-oxidation process is described by means of a two-step reaction mechanism. The model predictions of the lean ignition limits of premixed turbulent propane-air flames are verified through comparisons with data in the literature. The variations of the lean limits are found to compare favorably with experimental magnitudes and trends.

Radhakrishnan, K.↗

Statistics and thermodynamics of fracture

A probabilistic model of the fracture processes unifying the phenomenological study of long term strength of materials, fracture mechanics and statistical approaches to fracture is briefly outlined. The general framework of irreversible thermodynamics is employed to model the deterministic side of the failure phenomenon. The stochastic calculus is used to account for thg failure mechanisms controlled by chance; particularly, the random roughness of fracture surfaces.

Chudnovsky, A.↗