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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2015

Vapor–Liquid Equilibrium and Polarization Behavior of the GCP Water Model: Gaussian Charge-on-Spring versus Dipole Self-Consistent Field Approaches to Induced Polarization

Here we implemented the Gaussian charge-on-spring (GCOS) version of the original self-consistent field implementation of the Gaussian Charge Polarizable water model and test its accuracy to represent the polarization behavior of the original model involving smeared charges and induced dipole moments. Moreover, for that purpose we adapted the recently developed multiple-particle-move (MPM) within the Gibbs and isochoric-isothermal ensembles Monte Carlo methods for the efficient simulation of polarizable fluids. We also assessed the accuracy of the GCOS representation by a direct comparison of the resulting vapor-liquid phase envelope, microstructure, and relevant microscopic descriptors of water polarization along the orthobaric curve against the corresponding quantities from the actual GCP water model.

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

Elucidation of Peptide-Directed Palladium Surface Structure for Biologically Tunable Nanocatalysts

Peptide-enabled synthesis of inorganic nanostructures represents an avenue to access catalytic materials with tunable and optimized properties. This is achieved via peptide complexity and programmability that is missing in traditional ligands for catalytic nanomaterials. Unfortunately, there is limited information available to correlate peptide sequence to particle structure and catalytic activity to date. As such, the application of peptide-enabled nanocatalysts remains limited to trial and error approaches. In this paper, a hybrid experimental and computational approach is introduced to systematically elucidate biomolecule-dependent structure/function relationships for peptide-capped Pd nanocatalysts. Synchrotron X-ray techniques were used to uncover substantial particle surface structural disorder, which was dependent upon the amino acid sequence of the peptide capping ligand. Nanocatalyst configurations were then determined directly from experimental data using reverse Monte Carlo methods and further refined using molecular dynamics simulation, obtaining thermodynamically stable peptide-Pd nanoparticle configurations. Sequence-dependent catalytic property differences for C-C coupling and olefin hydrogenation were then eluddated by identification of the catalytic active sites at the atomic level and quantitative prediction of relative reaction rates. This hybrid methodology provides a clear route to determine peptide-dependent structure/function relationships, enabling the generation of guidelines for catalyst design through rational tailoring of peptide sequences

atomic pair distribution function↗
DOE OSTI2023

Local Chemical Clustering Enabled Ultrahigh Capacitive Energy Storage in Pb-Free Relaxors

Designing Pb-free relaxors with both a high capacitive energy density (W rec ) and high storage efficiency (η) remains a remarkable challenge for cutting-edge pulsed power technologies. Local compositional heterogeneity is crucial for achieving complex polar structure in solid solution relaxors, but its role in optimizing energy storage properties is often overlooked. Here, in this paper, we report that an exceptionally high W rec of 15.2 J cm –3 along with an ultrahigh η of 91% can be achieved through designing local chemical clustering in Bi 0.5 Na 0.5 TiO 3 –BaTiO 3 -based relaxors. A three-dimensional atomistic model derived from neutron/X-ray total scattering combined with reverse Monte Carlo method reveals the presence of subnanometer scale clustering of Bi, Na, and Ba, which host differentiated polar displacements, and confirming the prediction by density functional theory calculations. This leads to a polar state with small polar clusters and strong length and direction fluctuations in unit-cell polar vectors, thus manifesting improved high-field polarizability, steadily reduced hysteresis, and high breakdown strength macroscopically. The favorable polar structure features also result in a unique field-increased η, excellent stability, and superior discharge capacity. Our work demonstrates that the hidden local chemical order exerts a significant impact on the polarization characteristic of relaxors, and can be exploited for accessing superior energy storage performance.

25 ENERGY STORAGE↗
DOE OSTI2017

Nanoscale π–π stacked molecules are bound by collective charge fluctuations

Non-covalent π-π interactions are central to chemical and biological processes, yet the full understanding of their origin that would unite the simplicity of empirical approaches with the accuracy of quantum calculations is still missing. Here we employ a quantum-mechanical Hamiltonian model for van der Waals interactions, to demonstrate that intermolecular electron correlation in large supramolecular complexes at equilibrium distances is appropriately described by collective charge fluctuations. We visualize these fluctuations and provide connections both to orbital-based approaches to electron correlation, as well as to the simple London pairwise picture. The reported binding energies of ten supramolecular complexes obtained from the quantum-mechanical fluctuation model joined with density functional calculations are within 5% of the reference energies calculated with the diffusion quantum Monte-Carlo method. Our analysis suggests that π-π stacking in supramolecular complexes can be characterized by strong contributions to the binding energy from delocalized, collective charge fluctuations—in contrast to complexes with other types of bonding.

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

Interaction-driven spontaneous ferromagnetic insulating states with odd Chern numbers

Motivated by recent experimental work on moiré systems in a strong magnetic field, we compute the compressibility as well as the spin correlations and Hofstadter spectrum of spinful electrons on a honeycomb lattice with Hubbard interactions using the determinantal quantum Monte Carlo method. While the interactions in general preserve quantum and anomalous Hall states, emergent features arise corresponding to an antiferromagnetic insulator at half-filling and other incompressible states following the Chern sequence ± (2N + 1). These odd integer Chern states exhibit strong ferromagnetic correlations and arise spontaneously without any external mechanism for breaking the spin-rotation symmetry. Analogs of these magnetic states should be observable in general interacting quantum Hall systems. In addition, the interacting Hofstadter spectrum is qualitatively similar to the experimental data at intermediate values of the on-site interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Low-iridium stabilized ruthenium oxide anode catalyst for durable proton-exchange membrane water electrolysis

While mixing iridium (Ir) with ruthenium oxide (RuO 2 ) has proven to be an effective strategy for reducing Ir loading in anode catalysts for proton-exchange membrane (PEM) water electrolysers, achieving industrially relevant long-term stability typically requires an Ir-rich, Ru-lean combination. Here, in this study, by combining density functional theory with Metropolis Monte Carlo methods, we discovered that sufficient stabilization in the RuO 2 lattice could be achieved with less than 50 at.% of Ir, and that Ir in the first subsurface layer plays a critical role. By effectively dispersing Ir dopants within the RuO 2 lattice, we demonstrated an Ir:Ru atomic ratio of only 1:6 that exhibited exceptional stability for over 1,500 h of continuous water electrolysis at 2 A cm −2 . Our Ru 6 IrO x catalyst has the potential to reduce Ir loading by 80% compared with current commercial PEM water electrolysers, and its stability was further validated under industrial testing conditions in a 25-cm 2 PEM electrolyser.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Atomic-level structural correlations across the morphotropic phase boundary of a ferroelectric solid solution: xBiMg 1/2 Ti 1/2 O 3 -(1$-$x)PbTiO 3

Revelation of unequivocal structural information at the atomic level for complex systems is uniquely important for deeper and generic understanding of the structure property connections and a key challenge in materials science. Here in this paper we report an experimental study of the local structure by applying total elastic scattering and Raman scattering analyses to an important non-relaxor ferroelectric solid solution exhibiting the so-called composition-induced morphotropic phase boundary (MPB), where concomitant enhancement of physical properties have been detected. The powerful combination of static and dynamic structural probes enabled us to derive direct correspondence between the atomic-level structural correlations and reported properties. The atomic pair distribution functions obtained from the neutron total scattering experiments were analysed through big-box atom-modelling implementing reverse Monte Carlo method, from which distributions of magnitudes and directions of off-centred cationic displacements were extracted. We found that an enhanced randomness of the displacement-directions for all ferroelectrically active cations combined with a strong dynamical coupling between the A- and B-site cations of the perovskite structure, can explain the abrupt amplification of piezoelectric response of the system near MPB. Finally, altogether this provides a more fundamental basis in inferring structure-property connections in similar systems including important implications in designing novel and bespoke materials.

36 MATERIALS SCIENCE↗