Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION 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 937 records · Page 52

Calibration verification for stochastic agent-based disease spread models

Accurate disease spread modeling is crucial for identifying the severity of outbreaks and planning effective mitigation efforts. To be reliable when applied to new outbreaks, model calibration techniques must be robust. However, current methods frequently forgo calibration verification (a stand-alone process evaluating the calibration procedure) and instead use overall model validation (a process comparing calibrated model results to data) to check calibration processes, which may conceal errors in calibration. In this work, we develop a stochastic agent-based disease spread model to act as a testing environment as we test two calibration methods using simulation-based calibration, which is a synthetic data calibration verification method. The first calibration method is a Bayesian inference approach using an empirically-constructed likelihood and Markov chain Monte Carlo (MCMC) sampling, while the second method is a likelihood-free approach using approximate Bayesian computation (ABC). Simulation-based calibration suggests that there are challenges with the empirical likelihood calculation used in the first calibration method in this context. These issues are alleviated in the ABC approach. Despite these challenges, we note that the first calibration method performs well in a synthetic data model validation test similar to those common in disease spread modeling literature. We conclude that stand-alone calibration verification using synthetic data may benefit epidemiological researchers in identifying model calibration challenges that may be difficult to identify with other commonly used model validation techniques.

60 APPLIED LIFE SCIENCES↗

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Peierls transition, ferroelectricity, and spin-singlet formation in monolayer VOI 2

Using ab initio density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI 2 with a 3d 1 electronic configuration. Our phonon calculations indicate that the orthorhombic Pmm2 FE-II phase is the most likely ground state, involving a ferroelectric (FE) distortion along the a axis and V-V dimerization along the b axis. Specifically, the “pseudo Jahn-Teller” effect caused by the coupling between empty V (d xz/yz and d 3z 2 -r 2 ) and O 2p states is proposed as the mechanism that stabilizes the FE distortion from the paraelectric phase. Moreover, the half-filled metallic d xy band displays a Peierls instability along the b axis, inducing a V-V dimerization. Additionally, we found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly decoupled spin singlets in the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculated fission-fragment mass yields and average total kinetic energies of heavy and superheavy nuclei

Fission-fragment mass and total-kinetic-energy (TKE) distributions following fission of even-even nuclides in the region \(74 \le Z \le 126\) and \(92 \le N \le 230\) , comprising 896 nuclides have been calculated using the Brownian shape-motion method. The emphasis is the region of superheavy nuclei. To show compatibility with earlier results the calculations are extended to include earlier studied regions. An island of asymmetric fission is obtained in the superheavy region, \(106\le Z\le 114\) and \(162\le N\le 176\) , where the heavy fragment is found to be close to \(^{208}\) Pb and the light fragment adjusts accordingly. Most experimentally observed \(\alpha \) -decay chains of superheavy nuclei with \(Z > 113 \) terminate by spontaneous fission in our predicted region of asymmetric fission. In these cases, the pronounced large asymmetry is accompanied by a low TKE value compatible with measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SLICE

For most of their lifetime, pebble-bed reactors (PBRs) operate at an equilibrium state in which the core is filled with fuel pebbles at various levels of burnup. A fuel pebble travels multiple times in so-called passes through the reactor before it reaches its target discharge burnup and is replaced with a fresh fuel pebble. Given the stochastic nature of the fuel pebble travel paths and consequently the individual fuel pebble histories, it is not possible with standard methods developed for traditional reactor concepts to calculate the fuel inventory in the reactor core. An iterative approach, the SCALE Leap-In method for Cores at Equilibrium (SLICE), was developed to generate region-average fuel inventory for a PBR. The SLICE code enables automatic generation of input files for the SCALE code system (https://www.ornl.gov/scale), management of the SCALE result files, and analysis of results.

Bostelmann, Friederike [Oak Ridge National Laborat↗

Two-dimensional isometric tensor networks on an infinite strip

The exact contraction of a generic two-dimensional (2D) tensor network state (TNS) is known to be exponentially hard, making simulation of 2D systems difficult. The recently introduced class of isometric TNS (isoTNS) represents a subset of TNS that allows for efficient simulation of such systems on finite square lattices. The isoTNS ansatz requires the identification of an “orthogonality column” of tensors, within which one-dimensional matrix product state (MPS) methods can be used for calculation of observables and optimization of tensors. Here we extend isoTNS to infinitely long strip geometries and introduce an infinite version of the Moses Move algorithm for moving the orthogonality column around the network. Using this algorithm, we iteratively transform an infinite MPS representation of a 2D quantum state into a strip isoTNS and investigate the entanglement properties of the resulting state. In addition, we demonstrate that the local observables can be evaluated efficiently. Lastly, we introduce an infinite time-evolving block decimation algorithm (iTEBD 2 ) and use it to approximate the ground state of the 2D transverse field Ising model on lattices of infinite strip geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Level scheme of Nb 92 and observation of an oblate collective rotational band

Excited states of 92 Nb were investigated using the reactions 82 Se( 14 N, 4$n$) and 89 Y( 6 Li, 1 $p$2$n$). Based on these experimental results, the level scheme was extended abundantly. The structures of both the positive and negative parity states have been compared with shell model calculations using the SNET interaction. Potential energy surface calculation using the configuration-constrained method was performed to study a rotational-like band in 92 Nb, and possible configuration of the band head has been discussed. In conclusion, the calculation indicates that it may be an oblate collective rotational band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Comparison of As-loaded Dose Calculations to Measured Dose Rates

Demonstrating that the radiation dose fields surrounding spent nuclear fuel (SNF) storage and transportation systems meet the applicable limits set forth in 10 CFR 72 for storage and 10 CFR 71 for transportation is essential for the safe handling of radioactive material. The Used Nuclear Fuel-Storage, Transportation, and Disposal Analysis Resource and Data System (UNF-ST&DARDS) [2] is used to provide realistic estimates of SNF-related safety margins. The UNF-ST&DARDS dose rate analysis approach differs from that used in typical licensing approaches, which use design-basis assemblies with bounding source term characteristics to demonstrate the packaging design complies with the regulations. These bounding licensing approaches can increase the time needed to qualify fuel for loading into dry storage and the time a loaded system must wait prior to transportation. UNF-ST&DARDS dose rate assessments allow quantification of realistic, uncredited safety margins associated with actual fuel loading compared with the regulatory limits. While realistic estimation of the dose field surrounding SNF systems may allow for additional flexibility in operations, it is essential to understand how these predictions compare to measured doses. The U.S. Department of Energy Office of Integrated Waste Management and the Prairie Island Indian Community conducted a transportation dose assessment to estimate the site-specific incident-free radiation doses from shipping SNF by rail from the Prairie Island Nuclear Generating Plant (PINGP) through the Prairie Island Indian Community Reservation and Trust Land [5,6]. For that effort, the dose rates were obtained for 50 TN-40 and TN-40HT systems in storage configurations. This work compares the predicted dose rates from UNF-ST&DARDS as-loaded calculations with dose rates measured from 50 SNF storage systems at PINGP. The remainder of this paper discusses the data obtained for the evaluation, the modeling methods, and the results of the calculations.

spent nuclear fuel (SNF), UNF-ST&DARDS, Validation↗

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Density functional based methods for light harvesting molecules

The goal of the proposal was to computationally study the photoinduced charge transfer excited states of organic molecules. The systems of interest are the large organic supramolecules that have potential for photovoltaic applications. The simulation of the charge transfer excited states in large supramolecular assemblies require quantum mechanical methods that can be applied to systems with hundreds of atoms. Through this project we have implemented a perturbative delta-self-consistent field method that can easily describe the charge transfer excited states with regular density functionals without range separation. The cost of the calculations of an excited state is similar to that of the ground state. The method also allows for calculation of the energy gradients of the excited state. The developed method was applied to study a number of different supramolecular systems containing hundreds of atoms with excellent agreement with experiment. The absorption characteristics and the influence of density functional approximations were also examined. These led to several studies on the variety of organic electron donor and acceptor molecules. Apart from these developments through this project we also improved the code capabilities by implementing new meta-GGA functionals. This was a major step since the code is purely density functional and earlier the functional capabilities were limited only to generalized gradient approximations. A total of 23 peer-reviewed publications in top chemical physics journals resulted from this grant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT modelling of explicit solid–solid interfaces in batteries: methods and challenges

Density Functional Theory (DFT) calculations of electrode material properties in high energy density storage devices like lithium batteries have been standard practice for decades. By contrast, DFT modelling of explicit interfaces in batteries arguably lacks universally adopted methodology and needs further conceptual development. In this paper, we focus on solid–solid interfaces, which are ubiquitous not just in all-solid state batteries; liquid-electrolyte-based batteries often rely on thin, solid passivating films on electrode surfaces to function. We use metal anode calculations to illustrate that explicit interface models are critical for elucidating contact potentials, electric fields at interfaces, and kinetic stability with respect to parasitic reactions. The examples emphasize three key challenges: (1) the “dirty” nature of most battery electrode surfaces; (2) voltage calibration and control; and (3) the fact that interfacial structures are governed by kinetics, not thermodynamics. To meet these challenges, developing new computational techniques and importing insights from other electrochemical disciplines will be beneficial.

25 ENERGY STORAGE↗

Magnetic structure of A ≤ 10 nuclei using the Norfolk nuclear models with quantum Monte Carlo methods

Here we present quantum Monte Carlo calculations of magnetic moments, form factors, and densities of A ≤ 10 nuclei within a chiral effective field theory approach. We use the Norfolk two- and three-body chiral potentials and their consistent electromagnetic one- and two-nucleon current operators. We find that two-body contributions to the magnetic moment can be large (up to ≈ 33% in A = 9 systems). We study the model dependence of these observables and place particular emphasis on investigating their sensitivity to using different cutoffs to regulate the many-nucleon operators. Calculations of elastic magnetic form factors for A ≤ 10 nuclei show excellent agreement with the data out to momentum transfers q ≈ 3 fm -1 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influence of substituents on the vectorial difference static dipole upon excitation in synthetic bacteriochlorins

Organic dye aggregates have been shown to exhibit exciton delocalization in natural and synthetic systems. Such dye aggregates show promise in the emerging area of quantum information science (QIS). We believe that the difference static dipole (Δd) is an essential dye parameter in the development of molecular QIS systems. However, a foundational understanding of the structural factors influencing Δd remains elusive. Bacteriochlorins play a vital role in photosynthesis due to their exceptional photophysical properties. Therefore, bacteriochlorins are particularly suitable dyes for the construction of aggregate systems for QIS. Synthetic bacteriochlorins further offer stability and tunability via chemical modifications. Here, the influence of substituents on the Δd of monomeric (non-aggregated) dyes was investigated via density functional theory (DFT) and time-dependent (TD-)DFT in a set of 5-methoxybacteriochlorins progressively substituted with ethynyl, phenyl, and phenylethynyl substituents at the 3,13- and 3,13,15-positions of the macrocycle. Symmetrically substituted 5-methoxybacteriochlorins were shown to have the largest Δd. The increase in Δd in the series of dyes was largely due to changes in orientation of the static dipole upon excitation rather than large changes in magnitude. In addition, the transition dipole (μ) and the angle between Δd and μ (ζ) were calculated. Three 5-methoxybacteriochlorins with high predicted Δd and μ were synthesized and characterized spectroscopically. The trend in Δd values empirically determined using the solvatochromic Stokes shift method was comparable with the DFT calculations.

36 MATERIALS SCIENCE↗

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗

Fragmentation dynamics through geometrical distortion in low-energy electron attachment to carbon disulfide

The dissociation dynamics of dissociative electron attachment to carbon disulfide is studied using the velocity map imaging technique. Ion yields of the S− and CS− fragments are measured as a function of incident electron energy between 2 and 12 eV. Energy and angle-differential yields of the S− fragments are analyzed to investigate the dynamics of the CS2− anion around the 3.8-eV peak. The resonances involved in this fragmentation are identified and characterized in calculations using the well-known R-matrix method. Based on our measurements and the geometry dependence of the calculated resonances, we propose that a geometrical distortion of the transient negative-ion state of CS2 is necessary for S-CS bond dissociation to occur.

Paul, Anirban↗

The development and application of the stirred‐reactor coupon analysis (SRCA) test method

A new technique, termed the stirred‐reactor coupon analysis (SRCA) method, has been developed to measure the rate of glass dissolution in forward‐rate conditions. Monolithic glass coupons are partially masked with an inert material before placement in a large volume of well‐mixed solution with known chemistry and temperature for a predetermined duration. After the test, the mask is removed, and the difference in step height between the protected area and the exposed corroded portions of the sample coupon is measured to determine the extent of glass dissolution. The step height is converted to a rate measurement using the test duration and glass density. Test parameters such as sample surface preparation and test duration were evaluated to determine their effects on the measured rates. Additionally, results from an interlaboratory study (ILS) consisting of 12 laboratories from 11 different institutions are presented, where each laboratory performed 12 independent tests. When removing experimental outlier data, the 95% reproducibility limits for the SRCA method has no statistical difference with previously published standardized test methods used to determine the forward rate of glass dissolution. Overall, this paper describes steps necessary to perform the test method and provides the statistical calculations to evaluate test accuracy.

chemical durability↗