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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microscopic study of the fusion reactions 40,48 Ca+ 78 Ni and the effect of the tensor force

In this work, we provide a microscopic description of the fusion reactions between 40,48 Ca and 78 Ni. The internuclear potentials are obtained using the density-constrained (DC) time-dependent Hartree-Fock (TDHF) approach and fusion cross sections are calculated via the incoming wave boundary condition method. By performing DC-TDHF calculations at several selected incident energies, the internuclear potentials for both systems are obtained and the energy-dependence of fusion barrier are revealed. The influence of tensor force on internuclear potentials of 48 Ca+ 78 Ni is more obvious than those of 40 Ca+ 78 Ni. By comparing the calculated fusion cross sections between 40 Ca+ 78 Ni and 48 Ca+ 78 Ni, an interesting enhancement of subbarrier fusion cross sections for the former system is found, which can be explained by the narrow width of internuclear potential for 40 Ca+ 78 Ni while the barrier heights and positions are very close to each other. The tensor force suppresses the subbarrier fusion cross sections of both two systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Parametric estimation of Poisson's ratio for thin hinged-hinged plates

Cornu's method is an elegant calculation besieged by an impractical approach to obtain accurate estimates for Poisson's ratio. Conventionally, Cornu's method requires several components, and each component can adversely affect the accuracy of the measurement. Furthermore, Cornu's conventional method requires a long beam because beams with short length-to-width ratios cause the estimate of Poisson's ratio to diverge from the true value of Poisson's ratio. We believe that, with the right modifications, Cornu's method can become an attractive approach to obtaining precise estimates for Poisson's ratio from mode shapes. Here we use finite element simulations to show how to use Cornu's method to estimate Poisson's ratio from a mode shape. Our modified Cornu's method removes knife-edges and loading components for a hinged-hinged plate under steady-state excitation. Given the true value of Poisson's ratio, as a performance specification, we show that simple parametric expressions can fit estimates for Poisson's ratio for different length-to-width ratios of thin hinged-hinged plates. Additionally, we show that estimates for Poisson's ratio from higher modes align with the results from the first mode and explain our expectation for this outcome. Furthermore, our results challenge the idea that anticlastic, monoclastic, and synclastic deformation uniquely correspond to positive, zero, and negative estimates of Poisson's ratio, respectively. With the emergence of materials-by-design, we expect that this parametric technique will be able to assist in experimental qualification of thin beam and plate structures with respect to the desired value of Poisson's ratio.

42 ENGINEERING↗
DOE OSTI2022

Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics

Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.

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

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
DOE OSTI2025

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗
DOE OSTI2021

A nodal transport method for coupled fast-thermal reactor analysis

We report that for efficient core analyses of coupled fast-thermal reactors, a new deterministic method has been developed based on the variational nodal transport method of VARIANT. A new multigroup cross section generation procedure was devised by combining Monte Carlo lattice calculations for thermal assemblies and the two-step procedure of the MC 2 -3 code for fast assemblies. To reduce assembly homogenization errors in nodal transport calculations, a new nodal equivalence method was developed based on the partial current discontinuity factor (PCDF) and incorporated in VARIANT. A practical procedure to calculate PCDFs with fixed source supercell calculations was also devised. The performance of the proposed method was investigated using a test problem derived from the versatile coupled test reactor (VCTR) design. The new procedure produced multigroup cross sections accurately for all assemblies. VARIANT transport calculations with PCDFs produced accurate multiplication factor and power distribution compared to reference Serpent-2 Monte Carlo solutions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Correction: Enhanced 5f-δ bonding in [U(C 7 H 7 ) 2 ] - : C K-edge XAS, magnetism, and ab initio calculations

The authors regret that two references to the multi-configurational pair-density functional theory methods used in the wavefunction calculations did not appear in the main article. The references were originally included in the Methods section of the Electronic Supplementary Information as ref. 17a and b, but should also have appeared in the second paragraph on page 9564. The references are listed below as ref. 1 and 2. The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers.

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

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗
DOE OSTI2023

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

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

Performing k eff Validation of As-Loaded Criticality Safety Calculations Using UNF-ST&DARDS: Applicable Experiment Selection

The general method for performing validation of as loaded criticality safety calculations using UNF ST&DARDS is presented in a paper by Clarity, which includes a description of the UNF-ST&DARDS system. Proof-of-principle analyses were performed in the summer of 2019 for MPC-32 dual purpose canisters (DPCs) containing pressurized water reactor (PWR) fuel assemblies. Summaries of these results are presented in this and a companion paper for this conference. The current paper describes the TSUNAMI-IP calculations performed to select applicable experiments for validation of 11 MPC-32 DPCs. The companion paper discusses the TSUNAMI-3D calculations used to generate sensitivity data to support the experiment selections discussed here. Experiment selection is based on the sensitivity/uncertainty (S/U) methods used to validate criticality safety calculations of as-loaded DPCs containing pressurized water reactor (PWR) spent nuclear fuel (SNF). This process has been demonstrated and is summarized in this paper. The approach is similar to that used in NUREG/CR-7109, which provides an approach for validation of PWR burnup credit (BUC), including major and minor actinides and major fission products. The premise of S/U-based validation is that applicable experiments—those having a similar bias to a given application system—will have similar sensitivities for each isotope and reaction in the two systems. It is assumed that cross sections with larger uncertainties are more likely to contain data errors which contribute to the bias. The integral index c k thus propagates the system sensitivities with the nuclear covariance data to calculate a correlation coefficient representing the similarity of the two systems. In this work, a c k value of 0.8 or higher is interpreted as identifying an experiment with sufficient similarity for use in validation. This paper presents a brief summary of the characteristics of the 11 MPC-32 DPCs used in the proof of-principle analysis for as-loaded criticality safety calculation validation and an overview of the critical experiment suite with which each of these DPC models was compared. A summary and discussion of c k results is also presented, followed by conclusions and a discussion of future work to be performed for validation of UNF ST&DARDS as-loaded criticality safety calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

ARDP Natrium Neutronic Methodology: Argonne Neutronic Assessment of ABR-1000

The Natrium Sodium-cooled Fast Reactor (SFR) concept developed by TerraPower, LLC was selected as one of the advanced reactor designs for demonstration under the Advanced Reactor Demonstration Program Natrium Demo project. In collaboration with TerraPower, the Argonne National Laboratory (ANL) team provide independent reviews of the neutronic methodology being used for the Natrium design. As part of this activity, both the TerraPower and ANL teams agreed to perform independent neutronics analyses of the ABR-1000 reactor problem as a demonstration on the usage of the neutronics methodologies. Comparison of the results from TerraPower and ANL would expose any methodological differences in the modeling approaches. This report describes and presents the Argonne analysis work on the ABR-1000, and will be provided to TerraPower, along with data files, to complete the comparison study. The comparison will cover the methodology to evaluate the reactivity coefficients for a safety analysis, the shutdown margin, and the reconstructed pin power distribution. The methods and models used to calculate the kinetic parameters and the reactivity coefficients (density, Doppler, axial/radial expansion, etc.) are described in detail with recommendations of specific method options and modelling techniques. The control system reactivity worth and the excess reactivity due to temperature defect are calculated for the shutdown margin analysis. The method to evaluate the excess reactivity is discussed in detail. The pin power reconstruction methodology used in SE2ANL, SE2RCT, and DASSH are discussed in this report. The pin power distribution evaluated by different methods are compared and discussed in the report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science↗
DOE OSTI2022

A Novel Data Segmentation Method for Data-driven Phase Identification

This paper presents a smart meter phase identification algorithm for two cases: meter-phase-label-known and meter-phase-label-unknown. To improve the identification accuracy, a data segmentation method is proposed to exclude data segments that are collected when the voltage correlation between smart meters on the same phase is weakened. Then, using the selected data segments, a hierarchical clustering method is used to calculate the correlation distances and cluster the smart meters. If the phase labels are unknown, a Connected-Triple-based Similarity (CTS) method is adapted to further improve the phase identification accuracy of the ensemble clustering method. The methods are developed and tested on both synthetic and real feeder data sets. Here, simulation results show that the proposed phase identification algorithm outperforms the state-of-the-art methods in both accuracy and robustness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗