Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculations”

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

Method for calculating radon activity and radon rejection using a beta-gamma detector

Radon interferes with concentration measurements used by atmospheric radioxenon systems. We demonstrate a method to quantify the amount of radon that is present in the detectors, the impact of radon activity on the minimum-detectable-concentrations, and how to determine the needed radon rejection levels. An example calculation shows a radon rejection level of 10 5 is sufficient to limit impact on the detector sensitivity. We anticipate this method will give analysist a better understanding of radon present in their measurements and allow system designers to tailor their systems’ radon rejection better for its location.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Interlaboratory Reproducibility of Contour Method Data Analysis and Residual Stress Calculation

While the contour method for residual stress assessment has developed rapidly, no published study documents its interlaboratory reproducibility. Objective: Here we report an initial reproducibility experiment focused on contour method data analysis and residual stress calculation. The experiment uses surface topography data from a physical process simulation of elastic-plastic beam bending. The simulation provides surface topography, for input to the contour method data analysis, as well as a known residual stress field with 130 MPa peak magnitude. Additionally, to increase realism, noise and specific artifacts are added to the topography data. A group of participants received the topography data (without the known residual stress), independently analyzed the data, and submitted results as a two-dimensional residual stress field. Analysis of submissions provides a group average residual stress field and the spatial distribution of reproducibility standard deviation. The group average residual stress agrees with the known stress in magnitude and spatial trend. The reproducibility standard deviation ranges from 2 to 54 MPa over the measurement plane, with an average of 5.4 MPa. Reproducibility standard deviation is smaller in the cross-section interior (≤ 5 MPa), modest near local extrema in the stress field (5 to 10 MPa), and larger near the cross-section boundaries (10 to 30 MPa). Overall, the largest values of reproducibility standard deviation (up to 54 MPa) occur in limited areas where artifacts had been added to the topography data; while some participants identified and removed these artifacts, some did not, leading to systematic differences that elevated the standard deviation.

42 ENGINEERING↗
DOE OSTI2020

Anharmonic multi-phonon nonradiative transition: An ab initio calculation approach

Nonradiative carrier recombinations at deep centers in semiconductors are of great importance for both fundamental physics and device engineering. Here, we provide a revised analysis of Huang’s original nonradiative multi-phonon (NMP) theory with ab initio calculations. First, we confirmed at the first-principles level that Huang’s concise formula gives the same results as the matrix-based formula, and that Huang’s high-temperature formula provides an analytical expression for the coupling constant in Marcus theory. Secondly, we correct for anharmonic effects by taking into account local phonon-mode variations for different charge states of a defect. The corrected capture rates for defects in GaN and SiC agree well with experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Accurate Density Calculation for Molten Slags in SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO-‘FeO’-‘Fe{sub 2}O{sub 3}’ Systems

A simple and accurate density calculation model for molten slags containing iron oxides in SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO-‘FeO’-‘Fe{sub 2}O{sub 3}’ systems is established in the present study, by extending a previous model developed for SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO systems. To quantitate the iron oxides in slags, a new method is proposed to recalculate all the iron oxides as ‘FeO’ and ‘Fe{sub 2}O{sub 3},’ which represent the iron oxides in equilibrium with metallic iron and the air atmosphere, respectively. This method is essentially different from most of the other models that dealt with ferrous (FeO) and ferric (Fe{sub 2}O{sub 3}) iron oxides. A total of 860 experimental data points for iron oxide-containing melts are critically collected for the optimization of the model parameters. The model achieves an excellent agreement with literature values with an average error of 1.96 pct and an overall absolute error of 0.092 g/cm{sup 3}. The compositional dependences from unary to multicomponent systems are discussed using the present model. The results show that the densities always increase with the increasing ‘FeO’ or ‘Fe{sub 2}O{sub 3},’ and the substitution of ‘FeO’ by ‘Fe{sub 2}O{sub 3}’ will cause a decrease in density. In the multicomponent systems, SiO{sub 2} tends to have a larger impact on the density of ‘Fe{sub 2}O{sub 3}’-containing systems compared with ‘FeO’-containing systems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Optimization of CCT Equations Using Calculated Grain Boundary Soluble Compositions for the Simulation of Austenite Decomposition of Steels

New CCT equations have been developed and optimized to simulate the start temperatures of the austenite decomposition process in low-alloyed steels using experimental CCT data published in the literature. Exceptionally, this optimization does not apply the nominal compositions of the steels, but the corresponding soluble compositions of the grain boundaries calculated using IDS software, depending on the reported austenitization treatments of the steels. These compositions, rather than the nominal ones, are expected to control the start of the austenite decomposition, which usually initiates at the grain boundaries. The new optimization treatment takes into account the solute microsegregation and the possible precipitate formation. Using IDS software, the new equations were validated with new experimental CCT data. Agreement was good not only for the austenite decomposition start temperatures, but also for the final phase fractions, indicating fairly reasonable predictions of phase transformation kinetics by the IDS. In addition, IDS simulations were compared with the experimental CCT data of five high-carbon steels, applying both the new equations based on grain boundary soluble compositions as well as the equations based on the nominal compositions. With the same experimental CCT data used in optimization, better agreement was obtained with the new equations, indicating the importance of determining the soluble compositions at the grain boundaries where the austenite decomposition process is likely to begin.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Temperature-Dependent Configurational Entropy Calculations for Refractory High-Entropy Alloys

The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N , the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N , and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.

36 MATERIALS SCIENCE↗
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

Status on lattice calculations of the proton spin decomposition

Abstract Lattice calculations of the proton spin components is reviewed. The lattice results of the quark spin from the axial-vector current matrix element at ∼ 0.3−0.4 is smaller than those from the constituent quark models. This is largely due to the fact that the vacuum polarization contribution from the disconnected insertion is negative. Its connection with the anomalous Ward identity is clarified and verified numerically. This resolves the contentious issue in the “proton spin crisis.” The glue spin and angular momentum are found to be large and there is notable contribution from the quark orbital angular momentum. Renormalization, mixing, and normalization of the quark and glue angular momenta are discussed. With sufficient precision, they can be compared with more precise experimental measurements when the electron-ion collider facility is available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Local order of orthorhombic weberite-type Y 3 TaO 7 as determined by neutron total scattering and density functional theory calculations

Weberite-type oxides are a family of oxides with A 3 BO 7 stoichiometry that can adopt different space groups depending on chemical composition. There is a discrepancy in previous studies as to whether Y 3 TaO 7 is the orthorhombic C222 1 or the Ccmm space group. Here, we describe the short- and long-range structural properties of weberite-type Y 3 TaO 7 using neutron total scattering data, which has a high sensitivity to the oxygen sublattice. Simultaneous analysis of both short- and long-range structural data via conventional Rietveld and small box modelling demonstrates that Y 3 TaO 7 is best modeled as C222 1 . While the Ccmm describes equally well the long-range structure, pair distribution function analysis revealed that the local atomic configuration can be best modelled as C222 1 . This is corroborated by first-principles calculations that confirm the energetic preference of the C222 1 over Ccmm. Further, neutron total scattering data are reported for the first time for the Y 3 TaO 7 weberite-type ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A thermochemical database from high-throughput first-principles calculations and its application to analyzing phase evolution in AM-fabricated IN718

A comprehensive thermochemical database is constructed based on high–throughput first-principles phonon calculations of over 3000 atomic structures in limited concentrations in Ni, Fe, and Co alloys involving a total of 26 elements including Al, B, C, Cr, Cu, Hf, La, Mn, Mo, N, Nb, O, P, Re, Ru, S, Si, Ta, Ti, V, W, Y, and Zr, providing thermochemical data largely unavailable from existing experiments. Here, the database can be employed to predict the equilibrium phase compositions and fractions directly from first-principles by minimizing the chemical potential of a multicomponent system with a fixed overall chemical composition and a fixed temperature. It is applied to the additively manufactured nickel-based IN718 superalloy to analyze the phase evolution with temperature. IN718 is known for its great performance in tensile, fatigue, creep, and rupture strength, combined with easy fabrication and corrosion resistance. In particular, we successfully predicted the formation of L1 0 -FeNi, γ’-Ni 3 (Fe,Al), α-Cr, δ-Ni 3 (Nb,Mo), γ”-Ni 3 Nb, and η-Ni 3 Ti at low temperatures (below 680 K), γ’-Ni 3 Al, δ-Ni 3 Nb, γ”-Ni 3 Nb, α-Cr, and γ-Ni(Fe,Cr,Mo) at intermediate temperatures (between 680 and 1140 K), and δ-Ni 3 Nb and γ-Ni(Fe,Cr,Mo) at high temperatures (above 1140 K) in IN718. These predictions are validated by EDS mapping of compositional distributions and corresponding identifications of phase distributions. The database is expected to be a valuable source for future thermodynamic analysis and microstructure prediction of alloys involving the 26 elements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗
DOE OSTI2020

Lattice physics calculations using the embedded self-shielding method in polaris, Part II: Benchmark assessment

Polaris is a 2D multigroup lattice physics capability in the SCALE code system used for the analysis of light water reactor fuel designs. In this work, verification and validation results of Polaris in SCALE 6.2 are presented. The SCALE/KENO continuous-energy Monte Carlo code is used as a reference to compare the results of more than 2,000 combinations of pressurized water reactor and boiling water reactor lattice designs and state points. Radiochemical assay measurements were also used to assess the accuracy of depletion calculations. Results show that Polaris, with ENDF/B-VII.1 cross section library, produces accurate lattice physics data for core simulators.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Data reduction in deterministic neutron transport calculations using machine learning

Neutron cross section matrices for fission and scattering data are required for each material, temperature, and enrichment level to calculate the neutron transport equation accurately. Here, this information can be a limiting factor when using the multigroup discrete ordinates (S N ) method when the number of energy groups is large. Machine Learning (ML) can be used to replace the need for the cross section matrices by reproducing the function that maps the scalar flux to the scattering and fission sources. Through the use of autoencoders and Deep Jointly-Informed Neural Networks (DJINN), the data storage requirements are reduced by 94% of the original data for a 618 group problem. This is accomplished while preserving the scalar flux, maintaining generality, and decreasing wall clock times.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The importance of delayed neutron precursors in gamma dose calculations for activated primary heat exchanger components in the Molten Salt Breeder Reactor

The Molten Salt Reactor (MSR) Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation program has supported the development of 3D Monte Carlo models of the Molten Salt Breeder Reactor (MSBR) over the last couple of years. This MSBR model was previously run with the Shift Monte Carlo code to perform radiation shielding calculations in the reactor cell area. The MSBR is a 2250 MWth (1000 MWe) liquid-fueled molten salt reactor design developed at Oak Ridge National Laboratory in the 1970s. Determining the source terms from activated primary heat exchanger (HX) components is important because delayed and prompt neutron fluxes incident on these components affect the dose rate in the primary HX maintenance areas. This information can be used in the development of remote handling procedures required during shutdown for maintenance. A methodology has been developed and is proposed in this paper to quantify the activated source term from the primary HX components as a result of the movement of the delayed neutron precursors in flowing primary fuel salt through the primary HXs in the MSBR. The goal of this research is to evaluate the gamma dose rates in the maintenance hatches above the primary HX using the activated HX source terms. The study showed that the gamma dose rates are approximately two orders of magnitudes higher when accounting for the neutron activation from the movement of delayed neutron precursors through the HXs than when flowing fuel is not considered. Thus, the movement of delayed neutron precursors must be taken into account for accurately predicting the neutron activation of primary loop components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗
DOE OSTI2021

A theoretical calculation of stacking fault energy of Ni alloys: The effects of temperature and composition

Combining cluster expansion (CE) method with one dimensional axial Ising model, this work investigated the effects of alloying elements and configurational variations due to temperature on the stacking fault energy (SFE) of FCC Ni binary alloys. Ensembles of large numbers of atomistic structures, each with more than ~400 atoms, were generated to consider sufficient long-range chemical disorder and temperature effects due to configurational entropy. Here, a Monte Carlo Metropolis algorithm was used to generate these structures, whose energies were then evaluated based on the effective cluster interactions obtained from CE. As a baseline, this work had shown the SFE of pure Ni and Al to be 127 mJ/m 2 and 137 mJ/m 2 , respectively, which agreed with the experimental values of 125 mJ/m 2 and 150 mJ/m 2 reported in the literature. Additions of Al, Ti, Cr and Co to pure Ni were found to decrease the SFE to different extents. Although temperature does not strongly influence the SFE of the FCC Ni-Al and Ni-Cr binary alloys, it can lead to significant changes to the SFE of the FCC Ni-Ti and Ni-Co alloys. While effects of temperature and composition on SFE observed in this work were calculated from binary Ni alloys, the general trends are nonetheless expected to be valid in the γ phase of multicomponent Ni superalloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗