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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development of a Model Correlating Axial HEPA Filter Performance Between Elevated and Standard Conditions

The ASME AG-1 Code specifies minimum rated airflow based on standard cubic feet per minute (SCFM) for axial high efficiency particulate air (HEPA) filters and actual cubic feet per minute (ACFM) for radial HEPA filters with corresponding maximum pressure drop. A quantitative analysis can be performed to demonstrate the variation in filter efficiency and develop a correlation of performance between axial filters subjected to ACFM and SCFM. This study illustrates the variation in filter efficiency of axial HEPA filters subjected to ACFM and SCFM flow rate conditions at varying atmospheric pressures. It also presents the methods used to calculate correction factors to convert filter efficiency and pressure drop from elevated to standard conditions. The Axial Large Scale Test Stand (ALSTS) in the Institute for Clean Energy Technology (ICET) at Mississippi State University has been utilized to control aerosol concentration, flowrate, atmospheric pressure, relative humidity, and temperature. The filter test matrix consists of axial HEPA filters rated at 1,000 and 1,500 CFM with varying filter pack types. Each filter is subjected to varying pressures of 14.70 psi, 13.55 psi, 12.40 psi, and 11.25 psi, temperatures of 70 deg. F, 110 deg. F, and 150 deg. F, and relative humidities of 50% and 80%, resulting in twenty-four unique testing points. The corresponding SCFM flow rates are computed for each point in order to form a comparison. Evaluation of the effects of designating ACFM or SCFM will be based upon the filter efficiencies of generated Dioctyl Phthalate (DOP) particles at a diameter of 0.3 micrometers. The results show filter efficiency decreasing as atmospheric pressure decreases due to increased filter face velocity. A correlation model is developed to convert filter efficiency and pressure drop at elevated conditions to those of standard pressure, relative humidity, and temperature.

42 ENGINEERING↗
DOE OSTI2023

3D full core BEAVRS using OpenMOC with transport equivalence

Using an optimized implementation of the three-dimension method of characteristics for neutron transport, as well as a novel equivalence method for transport calculations de- signed to recover self-shielding errors from neglecting the angular dependence of reso- nant group absorption, a 3D full core light water reactor hybrid stochastic-deterministic eigenvalue calculation was achieved. This paper presents the optimizations developed and compares the transport solutions obtained. It realizes a runtime near 6,000 CPU hours for the statepoint, an improvement of an order of magnitude on previous works, with less than 1% error on pin fission/U-238 capture ratios and a few dozen pcms on the eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗
DOE OSTI2021

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

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

Microscopic calculations of nuclear level densities with the Lanczos method

A new method for computing the density of states in nuclei making use of an extrapolated form of the tridiagonal matrix obtained from the Lanczos method is presented. It will be shown that the global, average properties of the entire Lanczos matrix can be predicted from just four Lanczos iterations. The extrapolated Lanczos matrix (ELM) approach provides for an accurate computation of the density of states described within the configuration space, which, in some cases, is sufficient to accurately calculate the density of states at, or near, the neutron separation energy. Comparisons between theory and experiment are shown for $^{57}Fe$, $^{74}Ge$, and $^6Ge$. In addition, we show results for the J -dependence of moments and the level density for these three nuclei.

Shell model↗
DOE OSTI2025

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
DOE OSTI2020

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

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

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

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

PI3NN: Out-of-distribution-aware Prediction Intervals from Three Neural Networks

We propose a novel prediction interval (PI) method for uncertainty quantification, which addresses three major issues with the state-of-the-art PI methods. First, existing PI methods require retraining of neural networks (NNs) for every given confidence level and suffer from the crossing issue in calculating multiple PIs. Second, they usually rely on customized loss functions with extra sensitive hyperparameters for which fine tuning is required to achieve a well-calibrated PI. Third, they usually underestimate uncertainties of out-of-distribution (OOD) samples leading to over-confident PIs. Our PI3NN method calculates PIs from linear combinations of three NNs, each of which is independently trained using the standard mean squared error loss. The coefficients of the linear combinations are computed using root-finding algorithms to ensure tight PIs for a given confidence level. We theoretically prove that PI3NN can calculate PIs for a series of confidence levels without retraining NNs and it completely avoids the crossing issue. Additionally, PI3NN does not introduce any unusual hyperparameters resulting in a stable performance. Furthermore, we address OOD identification challenge by introducing an initialization scheme which provides reasonably larger PIs of the OOD samples than those of the in-distribution samples. Benchmark and real-world experiments show that our method outperforms several state-of-the-art approaches with respect to predictive uncertainty quality, robustness, and OOD samples identification.

Liu, Siyan↗
DOE OSTI2023

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ZrCoBi based half Heuslers with high thermoelectric conversion efficiency

A ZrCoBi-based p-type half-Heusler material can have a formula: ZrCoBi$_{1-x-y}$Sn$_{x}$Sb$_{y}$, where x can vary between 0.01 and 0.25, and y can vary between 0 and 0.2. An average dimensionless figure-of-merit (ZT) for the material can be greater than or equal to about 0.80 as calculated by an integration method for temperatures between 300 and 973 K. A ZrCoBi-based n-type half-Heusler material can have a formula: ZrCo$_{1-x}$Ni$_{x}$Bi$_{1-y}$Sb$_{y}$, where x can vary between 0.01 and 0.25, and y can vary between 0 and 0.3. The material has an average dimensionless figure-of-merit (ZT) is greater than or equal to about 0.65 as calculated by an integration method for temperatures between 300 and 973 K.

Ren, Zhifeng↗
DOE OSTI2025

Structural Changes in Metal Chalcogenide Nanoclusters Associated with Single Heteroatom Incorporation

Atomically precise nanoclusters (NCs) are promising building blocks for designing materials and interfaces with unique properties. By incorporating heteroatoms into the core, the electronic and magnetic properties of NCs can be precisely tuned. To accurately predict these properties, density functional theory (DFT) is often employed, making the rigorous benchmarking of DFT results particularly important. In this study, we present a benchmarking approach based on metal chalcogenide NCs as a model system. We synthesized a series of bimetallic, iron-cobalt chalcogenide NCs [Co 6-x Fe x S 8 (PEt 3 ) 6 ] + (x = 0-6) (PEt = triethyl phosphine) and investigated the effect of heteroatoms in the octahedral metal chalcogenide core on their size and electronic properties. Using ion mobility-mass spectrometry (IM-MS), we observed a gradual increase in the collision cross section (CCS) with an increase in the number of Fe atoms in the core. DFT calculations combined with trajectory method CCS simulations successfully reproduced this trend, revealing that the increase in cluster size is primarily due to changes in metal-ligand bond lengths, while the electronic properties of the core remain largely unchanged. Moreover, this method allowed us to exclude certain multiplicity states of the NCs, as their CCS values were significantly different from those predicted for the lowest-energy structures. Here, this study demonstrates that gas-phase IM-MS is a powerful technique for detecting subtle size differences in atomically precise NCs, which are often challenging to observe using conventional NC characterization methods. Accurate CCS measurements are established as a benchmark for comparison with theoretical calculations. The excellent correspondence between experimental data and theoretical predictions establishes a robust foundation for investigating structural changes of transition metal NCs of interest to a broad range of applications.

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

Thermal expansion and phonon anharmonicity of cuprite studied by inelastic neutron scattering and ab initio calculations

Inelastic neutron scattering measurements were performed with a time-of-flight chopper spectrometer to observe phonons in all parts of the Brillouin zone of a single crystal of cuprite Cu 2 O. We reduced the experimental data to phonon dispersions in the high-symmetry directions, and changes between 10 and 300 K are reported. In this paper, we show ab initio quasiharmonic (QH) and anharmonic (AH) calculations of phonon dispersions. We performed all AH calculations with a temperature-dependent effective potential method. Both QH and AH calculations account for the small negative thermal expansion of cuprite at low temperatures. However, the measured temperature-dependent phonon behavior was predicted more accurately with the AH calculations than the QH ones. Nevertheless, at 300 K, the cubic AH used in this paper did not entirely account for the experimental phonon dispersions in cuprite.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Hydrological control of chemical weathering and rock-carbon fluxes: East River, Colorado

This data package is used in the manuscript entitled “Hydrological control of chemical weathering and rock-carbon fluxes”. The field study was conducted in a lower montane hillslope of the East River watershed, underlain by Mancos Shale, within the lower 140 m section of a transect that extends nearly 1 km to its local peak. The data in this package, in CSV file format, were collected from Fall 2016 to fall 2021, including depth-resolved dynamic water table depths, subsurface water fluxes, solid phase (soil to bedrock) chemical compositions, geochemical properties of porewater and pore-gas, including radio carbon concentrations. The detailed methods of field studies, laboratory chemical analyses, and calculations are described in the Methods section and in the dataset file: Wan_et_al_Methods.pdf. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the on-the-fly slowing down method for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Among various on-the-fly slowing down methods, the equivalent Dancoff factor cell method was employed. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for the very high temperature reactor (VHTR), high temperature test reactor (HTTR), and Empire benchmark problems with various resonance self-shielding conditions, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also show that the methodology works well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗
DOE OSTI2023

MOOSE Optimization Module Overview: Application to Residual Stress Inversion in Nuclear Fuel Plates

Manufacturing methods, material selections and unique plate geometries used in low enriched uranium plate fuels for the High Performance Research Reactor (HPRR) have led to concerns about the residual stresses effect on bond strength of the fuel/cladding interfaces and overall fuel performance under irradiation. The current HPRR plate fuel test program performed residual stress calculations using the contour method with experimental displacement data obtained by incrementally cutting the plate fuel. In this work, we use a gradient based inverse optimization module developed within the MOOSE framework to estimate normal and shear residual stresses within the plate fuels using displacement data taken from a cut placed in the HPRR fuel plate. These results are verified against residual stress calculations obtained from the contour method. The inversion protocol used in this work utilizes two meshes to separately discretize the physics and parameter space. This allows a fine mesh to be used to fully resolve the physics and geometry of the forward problem. A coarser mesh is used to resolve the parameter space with spatially varying mesh density determined from a sensitivity analysis. The adaptive resolution of the parameter space essentially has the effect of regularization as it is designed to achieve the bias-variance tradeoff.

42 ENGINEERING↗