Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
DOE OSTI2017

Topology optimization under uncertainty via non-intrusive polynomial chaos expansion

This paper presents a systematic approach for topology optimization under uncertainty that integrates non-intrusive polynomial chaos expansion with design sensitivity analysis for reliability-based and robust topology optimization. Uncertainty is introduced in loading and in geometry to address the manufacturing variability. The manufacturing variability is modeled via a thresholding technique in which the threshold field is represented by a reduced dimensional random field. Response metrics such as compliance and volume are characterized as polynomial chaos expansions of the underlying uncertain parameters thus allowing accurate and efficient estimation of statistical moments, failure probabilities and their sensitivities. The number of simulations is reduced for linear structures under loading uncertainty by means of superposition. Efficiency of the non-intrusive polynomial chaos approach is highlighted by comparison with the Monte Carlo method in terms of the number of simulations. To demonstrate the effect of uncertainty, optimized designs that consider uncertainty are compared to those that do not. Finally, comparisons of polynomial chaos expansion to existing analytical methods on a benchmark numerical example are also provided for reliability-based and worst case designs.

42 ENGINEERING↗
DOE OSTI2016

Multi-objective optimization techniques to design the Pareto front of organic dielectric polymers

Here, we present two Monte Carlo algorithms to find the Pareto front of the chemical space of a class of dielectric polymers that is most interesting with respect to optimizing both the bandgap and dielectric constant. Starting with a dataset generated from density functional theory calculations, we used machine learning to construct surrogate models for the bandgaps and dielectric constants of all physically meaningful 4-block polymers (that is, polymer systems with a 4-block repeat unit). We parameterized these machine learning models in such a way that the surrogates built for the 4-block polymers were readily extendable to polymers beyond a 4-block repeat unit. By using translational invariance, chemical intuition, and domain knowledge, we were able to enumerate all possible 4, 6, and 8 block polymers and benchmark our Monte Carlo sampling of the chemical space against the exact enumeration of the surrogate predictions. We obtained exact agreement for the fronts of 4-block polymers and at least a 90% agreement for those of 6 and 8-block polymers. We present fronts for 10-block polymer that are not possible to obtain by direct enumeration. We note that our Monte Carlo methods also return polymers close to the predicted front and a measure of the closeness. Both quantities are useful information for the design and discovery of new polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First operation and validation of simulations for the divertor cryo-vacuum pump in Wendelstein 7-X

Ten cryo-vacuum pumps (CVPs) were installed in the subdivertor region of each island divertor in the stellarator Wendelstein 7-X (W7-X) and operated for the first time during the recently completed plasma campaign OP2.1. A pumping speed of 70 ± 1 $\frac{m^3}{s}$ was measured during dedicated tests with known hydrogen gas injection. Based on a conductance model, the estimated pumping speed ranges from 86-93 $\frac{m^3}{s}$ for different sticking coefficients between 0.6 and 0.8. After completion of the initial tests the CVPs were operated successfully throughout the campaign, with regeneration performed once a week. Neutral gas pressures in the subdivertor in the range of 10 −4 mbar are well within the molecular flow regime and limit the particle exhaust capabilities of the CVPs. Simulations of the neutral gas pressure in the three-dimensional complex geometry of the subdivertor were performed using the DIVGAS code based on the direct simulation Monte Carlo method and a model implemented in the steady-state thermal package in ANSYS, which are in agreement with the measured values during plasma operation.

Cryo-vacuum pumping↗
DOE OSTI2017

Equilibrium clumped-isotope effects in doubly substituted isotopologues of ethane

We combine path-integral Monte Carlo methods with a new intramolecular potential energy surface to quantify the equilibrium enrichment of doubly substituted ethane isotopologues due to clumped-isotope effects. Ethane represents the simplest molecule to simultaneously exhibit 13C–13C, 13C–D, and D–D clumped-isotope effects, and the analysis of corresponding signatures may provide useful geochemical and biogeochemical proxies of formation temperatures or reaction pathways. Utilizing path-integral statistical mechanics, we predict equilibrium fractionation factors that fully incorporate nuclear quantum effects, such as anharmonicity and rotational-vibrational coupling which are typically neglected by the widely used Urey model. The magnitude of the calculated fractionation factors for the doubly substituted ethane isotopologues indicates that isotopic clumping can be observed if rare-isotope substitutions are separated by up to three chemical bonds, but the diminishing strength of these effects suggests that enrichment at further separations will be negligible. The Urey model systematically underestimates enrichment due to 13C–D and D–D clumped-isotope effects in ethane, leading to small relative errors in the apparent equilibrium temperature, ranging from 5 K at 273.15 K to 30 K at 873.15 K. We additionally note that the rotameric dependence of isotopologue enrichment must be carefully considered when using the Urey model, whereas the path-integral calculations automatically account for such effects due to configurational sampling. These findings are of direct relevance to future clumped-isotope studies of ethane, as well as studies of 13C–13C, 13C–D, and D–D clumped-isotope effects in other hydrocarbons.

Geochemistry & Geophysics↗
DOE OSTI2012

A Computational Model of Cellular Response to Modulated Radiation Fields

Purpose: To develop a model to describe the response of cell populations to spatially modulated radiation exposures of relevance to advanced radiotherapies. Materials and Methods: A Monte Carlo model of cellular radiation response was developed. This model incorporated damage from both direct radiation and intercellular communication including bystander signaling. The predictions of this model were compared to previously measured survival curves for a normal human fibroblast line (AGO1522) and prostate tumor cells (DU145) exposed to spatially modulated fields. Results: The model was found to be able to accurately reproduce cell survival both in populations which were directly exposed to radiation and those which were outside the primary treatment field. The model predicts that the bystander effect makes a significant contribution to cell killing even in uniformly irradiated cells. The bystander effect contribution varies strongly with dose, falling from a high of 80% at low doses to 25% and 50% at 4 Gy for AGO1522 and DU145 cells, respectively. This was verified using the inducible nitric oxide synthase inhibitor aminoguanidine to inhibit the bystander effect in cells exposed to different doses, which showed significantly larger reductions in cell killing at lower doses. Conclusions: The model presented in this work accurately reproduces cell survival following modulated radiation exposures, both in and out of the primary treatment field, by incorporating a bystander component. In addition, the model suggests that the bystander effect is responsible for a significant portion of cell killing in uniformly irradiated cells, 50% and 70% at doses of 2 Gy in AGO1522 and DU145 cells, respectively. This description is a significant departure from accepted radiobiological models and may have a significant impact on optimization of treatment planning approaches if proven to be applicable in vivo.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Efficient particle control in systems with large density gradients

Simulations of large density gradients present a number of challenges for direct Monte Carlo methods, since they lead to too few particles in dilute regions and prohibitively many in the dense regions. Here, we propose a particle control methodology that gives the user more control of the number of particles per cell by introducing a variable weight for each particle. The proposed scheme is based on the stochastic weighted particle method, requires only small modifications to DSMC, and exactly conserves mass, momentum, and energy. In validation tests of systems with density ratios of order 102-104, we observe 101-102 times less variance in the dilute region compared to a DSMC solution with the same number of system particles, while introducing a moderate additional computational cost.

DSMC↗
DOE OSTI2024

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A likelihood ratio anomaly detector for identifying within-perimeter computer network attacks

The rapid detection of attackers within firewalls of enterprise computer networks is of paramount importance. Anomaly detectors address this problem by quantifying deviations from baseline statistical models of normal network behavior and signaling an intrusion when the observed data deviates significantly from the baseline model. But, many anomaly detectors do not take into account plausible attacker behavior. As a result, anomaly detectors are prone to a large number of false positives due to unusual but benign activity. Our paper first introduces a stochastic model of attacker behavior which is motivated by real world attacker traversal. Then, we develop a likelihood ratio detector that compares the probability of observed network behavior under normal conditions against the case when an attacker has possibly compromised a subset of hosts within the network. Since the likelihood ratio detector requires integrating over the time each host becomes compromised, we illustrate how to use Monte Carlo methods to compute the requisite integral. We then present Receiver Operating Characteristic (ROC) curves for various network parameterizations that show for any rate of true positives, the rate of false positives for the likelihood ratio detector is no higher than that of a simple anomaly detector and is often lower. Finally, we demonstrate the superiority of the proposed likelihood ratio detector when the network topologies and parameterizations are extracted from real-world networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Application of the Hartmann–Tran profile to precise experimental data sets of 12 C 2 H 2

Self- and nitrogen-broadened line shape data for the P e (11) line of the ν₁ + ν₃ band of acetylene, recorded using a frequency comb-stabilized laser spectrometer, have been analyzed using the Hartmann–Tran profile (HTP) line shape model in a multispectrum fitting. In total, the data included measurements recorded at temperatures between 125 K and 296 K and at pressures between 4 and 760 Torr. New, sub-Doppler, frequency comb-referenced measurements of the positions of multiple underlying hot band lines have also been made. These underlying lines significantly affect the P e (11) line profile at temperatures above 240 K and poorly known frequencies previously introduced errors into the line shape analyses. Thus, the behavior of the HTP model was compared to the quadratic speed dependent Voigt profile (QSDVP) expressed in the frequency and time domains. A parameter uncertainty analysis was carried out using a Monte Carlo method based on the estimated pressure, transmittance and frequency measurement errors. From the analyses, the P e (11) line strength was estimated to be 1.2014(50) × 10 -20 in cm.molecules⁻¹ units at 296 K with the standard deviation in parenthesis. For analyzing these data, we found that a reduced form of the HTP, equivalent to the QSDVP, was most appropriate because the additional parameters included in the full HTP were not well determined. As a supplement to this work, expressions for analytic derivatives and a lineshape fitting code written in Matlab for the HTP are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Measurements of scattering and absorption properties of submillimeter bauxite and silica particles

Submillimeter solid particles have been considered as thermal storage media for concentrated solar power applications. Knowledge of the scattering and absorption properties of individual particles is crucial for modeling the radiative heat transfer of the particle bed. In this work, a laser scatterometer is used to measure the single-particle scattering properties at a wavelength of 635 nm by using two configurations: (1) a falling particle curtain and (2) a taped particle layer. Because a one-particle nominal thickness is formed with area fractions of 5-55% depending on the configuration, multiple scattering is minimized and hence the single scattering phase function, averaged over all illuminated particles, is directly measured. Bauxite-based ceramic particles that are strongly absorbing in the solar spectrum and silica particles that are nonabsorbing in the visible and near-infrared are investigated. The directional-hemispherical reflectance and transmittance of the taped particles are also measured to deduce the forward and backward scattering efficiency factors and the absorption efficiency factors. Only weak wavelength dependence is observed in the measured region from 380 nm to 1020 nm. Furthermore, the scattering phase functions of all bauxite-based particles with varying sizes and compositions are very similar and can be fitted to a Henyey-Greenstein phase function with an asymmetry factor g = -0.20. For the silica particles, forward scattering dominates and g = 0.45 yields the best fit. A Monte Carlo method is developed to model the particle scattering characteristics, and reasonable agreements between the modeling and experimental results are observed by introducing a specularity parameter.

42 ENGINEERING↗
DOE OSTI2015

Atomic charges for modeling metal–organic frameworks: Why and how

Atomic partial charges are parameters of key importance in the simulation of Metal–Organic Frameworks (MOFs), since Coulombic interactions decrease with the distance more slowly than van der Waals interactions. But despite its relevance, there is no method to unambiguously assign charges to each atom, since atomic charges are not quantum observables. There are several methods that allow the calculation of atomic charges, most of them starting from the electronic wavefunction or the electronic density or the system, as obtained with quantum mechanics calculations. In this work, we describe the most common methods employed to calculate atomic charges in MOFs. In order to show the influence that even small variations of structure have on atomic charges, we present the results that we obtained for DMOF-1. We also discuss the effect that small variations of atomic charges have on the predicted structural properties of IRMOF-1. - Graphical abstract: We review the different method with which to calculate atomic partial charges that can be used in force field-based calculations. We also present two examples that illustrate the influence of the geometry on the calculated charges and the influence of the charges on structural properties. - Highlights: • The choice of atomic charges is crucial in modeling adsorption and diffusion in MOFs. • Methods for calculating atomic charges in MOFs are reviewed. • We discuss the influence of the framework geometry on the calculated charges. • We discuss the influence of the framework charges on structural the properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Simulation of energy loss of fractionally charged particles using Geant4

Over the years several rare-event search experiments have shown increasing interest to search for Fractionally Charged Particles (FCPs) within the data acquired in those experiments. This, in turn, has required the need for a framework to calculate accurately the energy loss expected from FCPs as they pass through the detector material. We have developed Monte Carlo methods within Geant4 to simulate for the first time energy loss of FCPs through different physics processes. The program can be used to calculate energy deposition distributions, a major component of any FCP search analysis, as they interact with detector material. We discuss the implementation and validation of the energy loss of FCPs using a typical detector material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Prognostic model and failure mechanisms of steam generators in Sodium-Cooled fast reactors

This paper presents a prognostic model for sodium-cooled fast reactor (SFR) steam generators (SGs). Here, the purpose of the model is to estimate the remaining useful life of SFR SGs and thus to support the decision-making of autonomous control. SFR SGs are of great interest for plant integrity due to their harsh operation environment. They operate at higher temperatures and higher coolant-to-steam pressure differences than those of current light water reactors (LWR). The severity of the SFR SG failure consequences, which include water-sodium contact, is another reason. Understanding its failure mechanisms is important for the development of its prognostic model. Based on our literature review and physics-based analysis, we concluded that creep would be a dominant degradation mode of SFR SGs due to SFR’s high-temperature and high-pressure environment. Thus, creep is the focus of the prognostic model. Various other failure modes were also investigated in this study. The mechanical fatigue due to flow-induced vibration is usually observed in early developmental phases and would not be an important issue during normal operation. The thermal fatigue due to thermal stripping is occasionally observed in other components in SFRs but does not affect SG integrity. Pure water stress corrosive crack, fretting, and so on are commonly observed in LWR SGs but are significantly less important in SFR SGs because of the high temperatures, high pressure differences, and chemical properties of liquid sodium. Based on these investigations, a prognostic model focusing on creep failures was developed. It estimates the failure probability profile by sampling the Larson-Miller parameter (LMP) and temperature and associated uncertainties through Monte Carlo methods. Two case studies were presented. The first one demonstrated the model’s capability to calculate the failure probability for a new specimen within a given time. If the tube is working under temperatures of 500C ± 3 and pressures of 15.2 MPa, which leads to LMP of 20,100 ± 50, the failure probabilities within 50, 70, and 90 years are approximately 0.1 %, 3.2 %, and 17.5 %. The second demonstrated how adjusting the workload can help to protect the integrity of the component. For an old tube reaching 99 % of its lifespan, continuing to run at temperatures of 500C ± 3 and LMP of 20,100 ± 50 leads to a failure probability of about 29.0 % within a year. If the temperature is reduced by 5C, the failure probability can be reduced to 2.5 %, and, if the pressure is also reduced such that the LMP is increased by 100, the failure probability can further be reduced to 0.12 %.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗
DOE OSTI2017

Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations

Small-angle X-ray scattering (SAXS) has become an increasingly popular technique for characterizing the solution ensemble of flexible biomolecules. However, data resulting from SAXS is typically low-dimensional and is therefore difficult to interpret without additional structural knowledge. In theory, molecular dynamics (MD) trajectories can provide this information, but conventional simulations rarely sample the complete ensemble. Here, we demonstrate that accelerated MD simulations can be used to produce higher quality models in shorter time scales than standard simulations, and we present an iterative Bayesian Monte Carlo method that is able to identify multistate ensembles without overfitting. This methodology is applied to several ubiquitin trimers to demonstrate the effect of linkage type on the solution states of the signaling protein. Here, we observe that the linkage site directly affects the solution flexibility of the trimer and theorize that this difference in plasticity contributes to their disparate roles in vivo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗