Efficient Simulation of Stiff Stochastic Differential Equations.
Abstract not provided.
Han, Xiaoying↗
SEARCH · Engineering Papers
Results for “Stochastic simulations”
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.
Validity of path thermodynamic description of reactive systems: Microscopic simulations
Traditional stochastic modeling of reactive systems limits the domain of applicability of the associated path thermodynamics to systems involving a single elementary reaction at the origin of each observed change in composition. An alternative stochastic modeling has recently been proposed to overcome this limitation. These two ways of modeling reactive systems are in principle incompatible. Here, the question thus arises about choosing the appropriate type of modeling to be used in practical situations. In the absence of sufficiently accurate experimental results, one way to address this issue is through the microscopic simulation of reactive fluids, usually based on hard-sphere dynamics in the Boltzmann limit. In this paper, we show that results obtained through such simulations unambiguously confirm the predictions of traditional stochastic modeling, invalidating a recently proposed alternative.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Machine Learning-Assisted High-Temperature Reservoir Thermal Energy Storage Optimization: Numerical Modeling and Machine Learning Input and Output Files
This data set includes the numerical modeling input files and output files used to synthesize data, and the reduced-order machine learning models trained from the synthesized data for reservoir thermal energy storage site identification. In this study, a machine-learning-assisted computational framework is presented to identify High-Temperature Reservoir Thermal Energy Storage (HT-RTES) site with optimal performance metrics by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems. All the simulations and the neural network model were done by Idaho National Laboratory. A detailed description of the work was reported in publication linked below.
15 GEOTHERMAL ENERGY↗
Simulation-based assessment on stochastic load scheduling for building cooling systems
Here, to fill knowledge gaps related to stochastic load scheduling, we performed a comprehensive evaluation of the stochastic load scheduling for building cooling systems. Specifically, we studied the common uncertain variables in the load scheduling process for building cooling systems and categorized those variables based on their dynamic patterns. We then developed a generic stochastic load scheduling framework and applied it to building cooling systems that served a simulated community. This community consists of 100 heterogeneous houses and serves as a virtual testbed for evaluating the performance of stochastic load scheduling. In this evaluation, we considered representatives of uncertain variables with different dynamic patterns and included 100 realizations of the considered uncertainty in the evaluation to better catch the probability distribution of the control performance. The evaluation results suggest that deterministic load scheduling can reduce the operating energy cost by 18% but its performance can be affected by uncertainty. Stochastic load scheduling can further decrease the operating energy cost under uncertainty compared to deterministic load scheduling. We also found that the effectiveness of stochastic load scheduling in handling uncertainty is not directly associated with the number of uncertainty scenarios that are considered in its formulation.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Particle Swarm Optimisation for group structure optimization for radiotherapy shielding
Neutron transport simulations are ubiquitous in nuclear engineering because they allow one to model experimental systems and render a model platform for easy perturbation of experimental designs. In addition, simulations allow one to gain experimental insight without actually having to go through the trouble of building a physical experiment. Neutron transport simulations can be stochastic or deterministic based. Stochastic neutron transport simulations are typically simulated using the Monte Carlo method and yield very accurate solutions but are computationally expensive, while deterministic methods are typically faster but can be less accurate. Here we focus on optimizing the accuracy of deterministic neutron transport simulations for radiotherapy simulations. Deterministic neutron transport requires discretization of angle, energy, and space to appropriately analyze the system one is trying to model. Discretization of energy is challenging because of the highly variable neutron flux at certain neutron energies. Improper discretization of energy in the transport model can lead to erroneous results and therefore inaccurate interpretations of the solution. In this study, we evaluate Particle Swarm Optimization (PSO) as a mechanism for selecting optimal group structures for radiotherapy shielding. We tested the particle swarm optimization algorithm on radiotherapy shielding problems using Los Alamos National Laboratory's (LANL) main deterministic transport code PARTISN. Results show that the optimized energy group structures generated from the optimization algorithm outperformed LANL's standard energy group structures, and therefore demonstrate utility in using PSO to expedite computation times due to the increased accuracy obtained with a smaller but optimized group structure. (authors)
43 PARTICLE ACCELERATORS↗
Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics
Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NRAP-open-IAM: A flexible open-source integrated-assessment-model for geologic carbon storage risk assessment and management
Large-scale implementation of geologic carbon storage (GCS) to help reduce atmospheric greenhouse gas emissions requires stakeholder confidence that injected CO2 will remain contained and that potential subsurface environmental risks are acceptably small and manageable. The U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) has developed an open-source integrated assessment model (NRAP-Open-IAM) to help address questions about a potential GCS site’s ability to effectively contain injected CO 2 and protect groundwater and other overlying environmentally sensitive receptors. NRAP-Open-IAM allows a user to: (1) incorporate relevant site geologic and injection scenario data; (2) characterize important site features and events;(3) couple fast prediction models of various system components of the engineered geologic system; and (4) execute stochastic, dynamic simulation of whole GCS system performance, leakage risk assessment, and uncertainty quantification. NRAP-Open-IAM is available on GitLab (https://gitlab.com/NRAP/OpenIAM), and is accompanied by multiple application examples and detailed user and developer guides.
54 ENVIRONMENTAL SCIENCES↗
Crystallization of trans-stilbene derivatives via bound-sphere kinetic Monte Carlo in SPPARKS (Stochastic Parallel PARticle Kinetic Simulator)
Explore the source record for details and available documents.
Janicki, Tesia Danielle↗
Thermally induced mimicry of quantum cluster excitations and implications for the magnetic transition in FePSe 3
In two dimensional magnets, the interplay of thermal fluctuations and spin anisotropy control the existence of long-range magnetic order. In the van der Waals antiferromagnets FePX 3 , orbital degeneracy in the 𝑡 2𝑔 levels of the Fe 2+ ions in octahedral coordination yields strong uniaxial anisotropy, which stabilizes magnetic order up to 𝑇 ≈ 100 K. Recent inelastic neutron scattering measurements around the magnetic ordering transition have shown the existence of a broad spectrum of magnetic fluctuations with nontrivial momentum dependence, which has been interpreted as evidence for localized entangled cluster excitations. In this paper, we offer an alternative interpretation using classical nonlinear spin dynamics simulations. We present stochastic Landau Lifshitz dynamics simulations that reproduce the neutron scattering measurements of Chen et al. [npj Quantum Mater. 9, 40 (2024)] on FePSe 3 . These calculations faithfully explain the dynamical structure factor's momentum and energy dependence and point to a classical origin for the excitations observed in neutron spectroscopy and that the order-disorder transition can be understood in terms of thermal fluctuations overcoming the anisotropy energy.
Landau-Lifschitz-Gilbert equation↗
"Hybrid fracture/matrix modeling for well completion options evaluation"
Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Here we present the development of a methodology to rapidly and efficiently numerically simulate mixed fracture-matrix flow systems for evaluation of well design and completion options. The methodology is based on a loose coupling framework, allowing the fracture and matrix systems to be meshed separately. The fracture system includes the integration of fracture growth and aperture data from well stimulation simulations of stochastically generated fracture networks. Automatic mesh refinement is used in the matrix simulation to resolve heat transfer near the fracture network. This simulation framework is used to efficiently determine optimal production and injection well placement using adaptive sampling.
15 GEOTHERMAL ENERGY↗
Exploiting stochastic locality in lattice QCD: hadronic observables and their uncertainties
Abstract Because of the mass gap, lattice QCD simulations exhibit stochastic locality: distant regions of the lattice fluctuate independently. There is a long history of exploiting this to increase statistics by obtaining multiple spatially-separated samples from each gauge field; in the extreme case, we arrive at the master-field approach in which a single gauge field is used. Here we develop techniques for studying hadronic observables using position-space correlators, which are more localized, and compare with the standard time-momentum representation. We also adapt methods for estimating the variance of an observable from autocorrelated Monte Carlo samples to the case of correlated spatially-separated samples.
Physics↗
Machine-learning-assisted high-temperature reservoir thermal energy storage optimization
High-temperature reservoir thermal energy storage (HT-RTES) has the potential to become an indispensable component in achieving the goal of the net-zero carbon economy, given its capability to balance the intermittent nature of renewable energy generation. In this study, a machine-learning-assisted computational framework is presented to co-optimize the performance metrics of HT-RTES by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Further, neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems.
15 GEOTHERMAL ENERGY↗
Hydroxyl radical yields in the heavy ion radiolysis of water
The yields of hydroxyl radicals in the radiolysis of water with protons and carbon ions have been examined using experimental scavenging techniques coupled with Monte Carlo track simulations. Combined with previous results using helium ions, this set of data gives valuable information on the potential for radiation induced damage to biological systems for light ions with very different track structures. Carbon dioxide production from aerated formic acid solutions in concentrations ranging from 10 -3 to 1 M was used as a probe of hydroxyl radical yields from about 8 ns to 8 µs. Numerical interpolation of the results at slightly different ion energies in combination with data from gamma radiolysis allows for a systematic analysis of both track average and track segment yields. As expected, considerable track chemistry is found to occur on the nanosecond to microsecond time scales. Monte Carlo track simulations employing stochastic diffusion-kinetic calculations of product yields are found to reproduce experimental observations satisfactorily. The track simulations are used to extract hydroxyl radical kinetics in pure water at neutral conditions.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Insight on electrolyte infiltration of lithium ion battery electrodes by means of a new three-dimensional-resolved lattice Boltzmann model
Electrolyte filling takes place between sealing and formation in Lithium Ion Battery (LIB) manufacturing process. This step is crucial as it is directly linked to LIB quality and affects the subsequent time consuming electrolyte wetting process. Although having fast, homogeneous and complete wetting is of paramount importance, this process has not been sufficiently examined and fully understood. For instance, experimentally available data is insufficient to fully capture the complex interplay upon filling between electrolyte and air inside the porous electrode. We report here for the first time a 3D-resolved Lattice Boltzmann Method (LBM) model able to simulate electrolyte filling upon applied pressure of LIB porous electrodes obtained both from experiments (micro X-ray tomography) and computations (stochastic generation, simulation of the manufacturing process using Coarse Grained Molecular Dynamics and Discrete Element Method). The model allows obtaining advanced insights about the impact of the electrode mesostructures on the speed of electrolyte impregnation and wetting, highlighting the importance of porosity, pore size distribution and pores interconnectivity on the filling dynamics. Furthermore, we identify scenarios where volumes with trapped air (dead zones) appear and evaluate the impact of those on the electrochemical behavior of the electrodes.
25 ENERGY STORAGE↗
Effect of diffusion constant on the morphology of dendrite growth in lithium metal batteries
Lithium dendrites can lead to a short circuit and battery failure, and developing strategies for their suppression is of considerable importance. In this work, we study the growth of dendrites in a simple model system where the solvent is a continuum and the lithium ions are hard spheres that can deposit by sticking to existing spheres or the electrode surface. Using stochastic dynamics simulations, we investigate the effect of applied voltage and diffusion constant on the growth of dendrites. We find that the diffusion constant is the most significant factor, and the inhomogeneity of the electric field does not play a significant role. Additionally, the growth is most pronounced when the applied voltage and diffusion constant are both low. We observe a structural change from broccoli to cauliflower shape as the diffusion constant is increased. The simulations suggest that a control of electrolyte parameters that impact lithium diffusion might be an attractive route to controlling dendrite growth.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Selective capture and recovery of uranium oxide colloids from aqueous soil suspensions using high gradient magnetic filtration
High Gradient Magnetic Filtration (HGMF) is a promising method for the selective capture and recovery of uranium oxide from surface soils. To date, however, magnetic filtration of uranium oxide has only been demonstrated at a proof-of-principle scale using relatively small filters (<5 cm 3 ) at low flowrates (<60 mL/min). Here, to explore the efficacy of magnetic filtration of uranium oxide at a larger scale, a newly designed HGMF apparatus that is more than an order of magnitude larger than our earlier filters (106 cm 3 ) was designed, fabricated, and tested at relatively high flowrates. Filtration experiments were performed using aqueous uranium oxide particle suspensions with Arizona Road Dust (ARD) as a soil simulant. At a flowrate of 125 mL/min, the apparatus’ uranium capture rate was exceptionally high (96 %), but selectivity was poor due to the high rate of capture for diamagnetic soil constituents (e.g., 77 % for silicon). All particles were captured at a lower rate when the flowrate was increased to 250 mL/min, but uranium selectivity was significantly increased due to the more substantial reduction in diamagnetic particle capture (i.e., capture rate of 77 % and 15 % for uranium and silicon, respectively). When backwashing the apparatus at the same flowrates used during filtration experiments, the rate of uranium recovery tended to be fairly low. Nevertheless, higher flowrates (1 L/min) and sonication were both shown to be highly effective methods of increasing uranium recovery. Magnetic field simulations were also performed to investigate potential optimizations to the design of the apparatus. These simulations showed that the intensity of the applied magnetic field could be increased by increasing the thickness of the steel magnetic housing. Additionally, stochastic trajectory simulations were performed to investigate the potential mechanisms of particle capture.
HGMF↗
Surrogate optimization of variational quantum circuits
Variational quantum eigensolvers are touted as a near-term algorithm capable of impacting many applications. However, the potential has not yet been realized, with few claims of quantum advantage and high resource estimates, especially due to the need for optimization in the presence of noise. Finding algorithms and methods to improve convergence is important to accelerate the capabilities of near-term hardware for VQE or more broad applications of hybrid methods in which optimization is required. To this goal, we look to use modern approaches developed in circuit simulations and stochastic classical optimization, which can be combined to form a surrogate optimization approach to quantum circuits. Using an approximate (classical CPU/GPU) state vector simulator as a surrogate model, we efficiently calculate an approximate Hessian, passed as an input for a quantum processing unit or exact circuit simulator. This method will lend itself well to parallelization across quantum processing units. We demonstrate the capabilities of such an approach with and without sampling noise and a proof-of-principle demonstration on a quantum processing unit utilizing 40 qubits.
Gustafson, Erik J. [RIACS, Mtn. View] (ORCID:00000↗
Effects of coupling a stochastic convective parameterization with the Zhang–McFarlane scheme on precipitation simulation in the DOE E3SMv1.0 atmosphere model
Abstract. A stochastic deep convection parameterization is implemented into the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1.0 (EAMv1). This study evaluates its performance in simulating precipitation. Compared to the default model, the probability distribution function (PDF) of rainfall intensity in the new simulation is greatly improved. The well-known problem of “too much light rain and too little heavy rain” is alleviated, especially over the tropics. As a result, the contribution from different rain rates to the total precipitation amount is shifted toward heavier rain. The less frequent occurrence of convection contributes to suppressed light rain, while more intense large-scale and convective precipitation contributes to enhanced heavy total rain. The synoptic and intraseasonal variabilities of precipitation are enhanced as well to be closer to observations. The sensitivity of the rainfall intensity PDF to the model vertical resolution is examined. The relationship between precipitation and dilute convective available potential energy in the stochastic simulation agrees better with that in the Atmospheric Radiation Measurement (ARM) observations compared with the standard model simulation. The annual mean precipitation is largely unchanged with the use of the stochastic scheme except over the tropical western Pacific, where a moderate increase in precipitation represents a slight improvement. The responses of precipitation and its extremes to climate warming are similar with or without the stochastic deep convection scheme.
54 ENVIRONMENTAL SCIENCES↗