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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Methods and devices for preventing computationally explosive calculations in a computer for model parameters distributed on a hierarchy of geometric simplices

A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.

Source record↗
DOE OSTI2020

Characterization of the β-Phrase Region of the SAS Compressor Bed 1 Material (LANA10)

Isotherm data have been acquired and analyzed on the SAS Compressor Bed 1 material in order to define a method to set the initial bed Q/M value when freshly filled after inventory. β-phase data from eleven protium and two deuterium isotherms have been combined to determine a method for calculating the bed Q/M value given a β-phase equilibrium pressure at temperature. This method is a two-step method; first, the plateau pressure is calculated based on the isotope used and the temperature of the bed, and second, the measured equilibrium pressure is divided by the computed plateau pressure to obtain the reduced pressure, which is then inserted into the β-phase equation to calculate the bed’s Q/M value. Error sources and effects are discussed, with the conclusion that the ‘true’ reduced isotherm has been determined in the β-phase which will allow virgin bed inventory to be determined to ~3-4% accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Calculating elements of matrix functions using divided differences

In this work, we introduce a method for calculating individual elements of matrix functions. Our technique makes use of a novel series expansion for the action of matrix functions on basis vectors that is memory efficient even for very large matrices. We showcase our approach by calculating the matrix elements of the exponential of a transverse-field Ising model and evaluating quantum transition amplitudes for large many-body Hamiltonians of sizes up to 2 64 x 2 64 on a single workstation. We also discuss the application of the method to matrix inverses. We relate and compare our method to the state-of-the-art and demonstrate its advantages. We also discuss practical applications of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

Here we report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

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

Examination of Digital-Delay Rossi-$\alpha$ for 252 Cf-Driven Highly Enriched Uranium Using Organic Scintillators

Neutron noise techniques constitute several analysis methods applicable to non-destructive assay. One technique is the Rossi-α method to calculate the prompt neutron decay constant (α) or its inverse, the prompt neutron period (1/α), for assemblies of fissionable material. This work evaluates a high-data-throughput measurement at the National Criticality Experiments Research Center (NCERC) with organic scintillators measuring a subcritical assembly of highly enriched uranium (HEU) metal (93% 235 U). The assembly comprises hemi shells stacked together to form fully closed shells and is driven by a Cf source at the center. The assembly is a total of 59.85 kg HEU and a k eff of 0.99, according to MCNP6.2 KCODE simulations. Measurements were acquired with a three-by-four array of 5.08-cm-diameter by 5.08-cm-length trans-stilbene crystals 166 cm from the assembly center. Two types of coincident binning methods are used to build the Rossi-α distribution of coincident neutron detections: A) type 1 binning, also known as any-and-all forward time differences, and B) type 1 binning with a digital-delay technique that is analogous to Orndoff’s use of delay cabling in early Rossi-alpha measurements. Method B) disregards same-detector coincidences and once all nearest time coincidences are collected between all detectors after a trigger, a time delay of 0.75 µs is implemented. The measured prompt neutron decay constants for both techniques are calculated from single exponential fits and the two methods are compared. Method A) shows an apparent timing discontinuity near 500 ns time differences. This timing discontinuity interferes with the fitting method used to calculate the prompt neutron decay constant. Method B) mitigates the artificial timing discontinuity and removes the disagreement of the fit. Future work will model the time-dependent detector response to discern why this timing discontinuity occurs, discern how Method B) reduces the amplitude of this timing discontinuity, and apply this method to critical assembly measurements to ultimately confirm the recommendation to use Method B) for high-data-throughput measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Rovibrationally resolved Rayleigh and Raman scattering cross sections for molecular hydrogen

Accurate Rayleigh and Raman scattering cross sections, tensor components, depolarization ratios, and reversal coefficients for all rovibrational transitions within the X1Σg+ ground electronic state of H2 have been calculated. Raman spectra have been generated using these data. A method for calculating Raman scattering cross sections is formulated that is valid below the ionization threshold and in the region containing resonances, which explicitly accounts for all bound and dissociative vibrational levels of the bound intermediate electronic states and approximately accounts for the ionization continuum. A representative set of cross sections is presented for incident photon energies below 15 eV and compared with existing results in the literature where possible. Convergence of our results with an increasing number of bound intermediate electronic states is demonstrated. The accuracy of the Placzek–Teller approximation is discussed. The effect of accounting for the intermediate ionization continuum is investigated. Local thermal equilibrium cross sections are calculated for Rayleigh and Raman scattering. This work represents the most accurate and complete treatment of Raman scattering for molecular hydrogen to date. A total of 9582 Rayleigh and Raman scattering cross sections have been generated and are openly available on Zenodo under an open-source Creative Commons Attribution license at https://zenodo.org/doi/10.5281/zenodo.13441471.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

The ABACUS cosmological N -body code

We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct 1/r 2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10 -5 . We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dispersion interactions in proposed covalent superhydride superconductors

Recent developments in high temperature superconducting materials under high pressure have made numerical evaluation of the superconducting transition temperature (T c ) of predicted materials critically important as a means to identification. Existing methods of calculating T c often do not agree with each other or experiment, often due to the large number of complex factors that contribute to this property; among them is the neglect of dispersion interactions in commonly used density functionals. We evaluate the effect of including dispersion interactions on the predicted superconducting properties of two examples of the covalent superhydride class of very high-T c superconducting materials. In both cases, dispersion is found to have sizable effects, increasing the electron-phonon coupling as compared to the reference case of elemental niobium. A detailed investigation traces the origin of this effect in a 270 GPa $\mathrm {R3m CSH_7}$ structure to structural distortions driven by long-range electron-phonon interactions rather than novel bonding networks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale, mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, the effective diffusion coefficient of the fission product Ag through the grain boundary (GB) network is calculated using a combination of atomistic and phase-field methods. Atomistic calculations of Ag diffusivity in SiC bulk and GBs are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Since GBs serve as pathways for Ag diffusion, Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Higher spin partition functions via the quasinormal mode method in de Sitter quantum gravity

In this note we compute the 1-loop partition function of spin- s s fields on Euclidean de Sitter space S^{2n+1} S 2 n + 1 using the quasinormal mode method. Instead of computing the quasinormal mode frequencies from scratch, we use the analytic continuation prescription L_{\text{AdS}}\to iL_{\text{dS}} L AdS → i L dS , appearing in the dS/CFT correspondence, and Wick rotate the normal mode frequencies of fields on thermal \text{AdS}_{2n+1} AdS 2 n + 1 into the quasinormal mode frequencies of fields on de Sitter space. We compare the quasinormal mode and heat kernel methods of calculating 1-loop determinants, finding exact agreement, and furthermore explicitly relate these methods via a sum over the conformal dimension. We discuss how the Wick rotation of normal modes on thermal \text{AdS}_{2n+1} AdS 2 n + 1 can be generalized to calculating 1-loop partition functions on the thermal spherical quotients S^{2n+1}/\mathbb{Z}_{p} S 2 n + 1 / ℤ p . We further show that the quasinormal mode frequencies encode the group theoretic structure of the spherical spacetimes in question, analogous to the analysis made for thermal AdS in [1-3] .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗
DOE OSTI2022

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

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

Hamiltonian Truncation Effective Theory

Hamiltonian truncation is a non-perturbative numerical method for calculating observables of a quantum field theory. The starting point for this method is to truncate the interacting Hamiltonian to a finite-dimensional space of states spanned by the eigenvectors of the free Hamiltonian H_0 H 0 with eigenvalues below some energy cutoff E_\text{max} E max . In this work, we show how to treat Hamiltonian truncation systematically using effective field theory methodology. We define the finite-dimensional effective Hamiltonian by integrating out the states above E_\text{max} E max . The effective Hamiltonian can be computed by matching a transition amplitude to the full theory, and gives corrections order by order as an expansion in powers of 1/E_\text{max} 1 / E max . The effective Hamiltonian is non-local, with the non-locality controlled in an expansion in powers of H_0/E_\text{max} H 0 / E max . The effective Hamiltonian is also non-Hermitian, and we discuss whether this is a necessary feature or an artifact of our definition. We apply our formalism to 2D \lambda\phi^4 λ ϕ 4 theory, and compute the the leading 1/E_\text{max}^2 1 / E max 2 corrections to the effective Hamiltonian. We show that these corrections nontrivially satisfy the crucial property of separation of scales. Numerical diagonalization of the effective Hamiltonian gives residual errors of order 1/E_\text{max}^3 1 / E max 3 , as expected by our power counting. We also present the power counting for 3D \lambda \phi^4 λ ϕ 4 theory and perform calculations that demonstrate the separation of scales in this theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Calculations of electronic excitation by protons and alpha particles in silicon

This paper presents an ab initio method of calculating the valence electron contribution to the stopping and interaction cross sections of swift ions in matter using linear response calculations of the dielectric function in crystals. We estimate the contribution to the stopping and cross sections by repeated electron capture and stripping of the projectile ion by the medium. We also predict the electronic stopping of protons and alpha particles in silicon and compare to measured values. The predictions generally agree with the data at ion speeds greater than 1 atomic unit. We find that electron capture and stripping is at most a small contribution to the stopping for both ions studied.

electron excitation, silicon, nuclear interactions↗