Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH 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.

936 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

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

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2023

Progress Report of Alloy 709 Performance in Liquid Sodium

With improved high-temperature properties, Alloy 709 shows great potential for advanced reactor applications. To understand the effects of sodium on A709, small tensile specimens were machined from two commercial heats and subject to sodium exposure in two forced convection loops operated by Argonne. Thermal aging experiments were also conducted in parallel to isolate any sodium-specific effects. An increasing trend in yield strength and a decrease trend in ultimate tensile strength can be seen with the increasing exposure time, suggesting microstructural changes at these temperatures. Nonetheless, the thermally aged and sodium exposed samples behaved similarly, implying no sodium-specific effect at these temperatures. In contrast, 316H SS exposed to sodium at 650°C showed a strong sodium effect. The different responses of A709 and 316H SS to sodium exposure is of interest and requires additional effort to understand. The cold traps of the sodium loop were replaced early in FY23. During a maintenance operation, an unexpected sodium leak incident occurred at one of the sodium loops, triggering an emergency response. An investigation into the incident revealed that the leak happened at a VCR fitting on the recently replaced cold trap. The metal gasket of the leaked VCR fitting was not compressed sufficiently. A dimensional mismatch between the new component and old piping configuration was responsible for a significant misalignment. A new cold trap with improved design is in production. The new design provides additional flexibility and an access point for post-installation helium leak check and for venting of trapped gas. These design features will make the system more error-tolerant and the installation quality verifiable. A thermodynamic analysis was performed to understand the carburization-decarburization behavior of A709 in sodium environment. The equilibrium carbon activity in the alloy was evaluated with M 23 C 6 , TiC, and NbC. The carburization-decarburization boundary was calculated with different methods developed previously. The results suggests that A709 will be carburized in SMT-1 and SMT-2, and can be decarburized in low-carbon-activity environments such as EBR-II at temperatures above 620-650°C. Nonetheless, A709 should be less likely to be decarburized than 316H SS when exposed to the same low-carbon-activity environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Hierarchical excitations from correlated spin tetrahedra on the breathing pyrochlore lattice

The hierarchy of the coupling strengths in a physical system often engenders an effective model at low energies where the decoupled high-energy modes are integrated out. Here, using neutron scattering, we show that the spin excitations in the breathing pyrochlore lattice compound CuInCr 4 S 8 are hierarchical and can be approximated by an effective model of correlated tetrahedra at low energies. At higher energies, intratetrahedron excitations together with strong magnon-phonon couplings are observed, which suggests the possible role of the lattice degree of freedom in stabilizing the spin tetrahedra. Finally, our work illustrates the spin dynamics in CuInCr 4 S 8 and demonstrates a general effective-cluster approach to understand the dynamics on the breathing-type lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗
DOE OSTI2025

Tuning the magnetic properties of the spin-split antiferromagnet MnTe through pressure

The hexagonal antiferromagnet MnTe has attracted enormous interest as a prototypical example of a spin-compensated magnet in which the combination of crystal and spin symmetries lifts the spin degeneracy of the electron bands without the need for spin-orbit coupling, a phenomenon called nonrelativistic spin splitting (NRSS). Subgroups of NRSS are determined by the specific spin-interconverting symmetry that connects the two opposite-spin sublattices. In MnTe, this symmetry is rotation, leading to the subgroup with spin splitting away from the Brillouin zone center, often called altermagnetism. MnTe also has the largest spontaneous magnetovolume effect of any known antiferromagnet, implying strong coupling between the magnetic moment and volume. This magnetostructural coupling offers a potential knob for tuning the spin-splitting properties of MnTe. Here, we use neutron diffraction with in situ applied pressure to determine the effects of pressure on the magnetic properties of MnTe and further explore this magnetostructural coupling. We find that applying pressure significantly increases the Néel temperature, but decreases the ordered magnetic moment. We explain this as a consequence of strengthened magnetic exchange interactions under pressure, resulting in higher 𝑇 N , with a simultaneous reduction of the local moment of individual Mn atoms, described here via density functional theory. This reflects the increased orbital hybridization and electron delocalization with pressure. In conclusion, these results shed light on the competition between magnetic exchange interactions and the strength of individual magnetic moments and show that the magnetic properties of MnTe can be controlled by pressure, opening the door to improved properties for spintronic applications through tuning via physical or chemical pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

First study of the $^{139}\text{Ba}(n,y)^{140}\text{Ba}$ reaction to constrain the conditions for the astrophysical $\mathcal{i}$ process

New astronomical observations point to a nucleosynthesis picture that goes beyond what was accepted until recently. The intermediate “i” process was proposed as a plausible scenario to explain some of the unusual abundance patterns observed in metal-poor stars. The most important nuclear physics properties entering i-process calculations are the neutron-capture cross sections and they are almost exclusively not known experimentally. In this report we provide the first experimental constraints on the 139 Ba(n,γ) 140 Ba reaction rate, which is the dominant source of uncertainty for the production of lanthanum, a key indicator of i-process conditions. This is an important step towards identifying the exact astrophysical site of stars carrying the i-process signature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Summary of Graphite Data Stored within NDMAS

The Graphite Technology Development Project provides data to support the design of graphite core components within specific reactor service conditions of the next generation of high-temperature, gas-cooled nuclear reactors. Physical, mechanical, and thermal properties of nuclear grade graphite were characterized for specimens that were unirradiated, irradiated, and irradiated under various stress conditions. The material properties include diameter, length, mass, density, compressive strength, tensile strength, flexural strength, modulus, resistivity, thermal diffusivity, and thermal expansion coefficient. Baseline graphite specimens are unirradiated from different grades (2114, IG-110, NBG-17, NBG-18, and PCEA) and different types used in different characterization tests (compressive, flexural, tensile, one-inch cylinder, and quarter-inch cylinder). The Advanced Graphite Creep (AGC) irradiation specimens are from a much larger number of grades and cylinder types (creep, piggyback, and pencil). For baseline graphite, characterization data for 7,756 specimens extracted from thirteen graphite billets were captured to the NDMAS database. For AGC experiments, four irradiation campaigns have been completed: AGC-1, AGC-2, AGC-3, and AGC-4. The ongoing HDG-1 (High-Dose Graphite) experiment, which began irradiation with Cycle 168B on August 26, 2020, includes specimens previously irradiated in AGC-2 in addition to the specimens originally destined for AGC-5.. Besides the characterization data, the AGC data includes irradiation monitoring and physics data representing the irradiation conditions of AGC specimens. Currently, all data for AGC-1, AGC-2, and AGC-3 have been captured to the NDMAS database. Only AGC-4 pre-irradiation and irradiation monitoring data have been captured, and HDG-1 pre-irradiation data are in process of being captured for unirradiated specimens. To date, a total of 38,149 material property records have been captured into NDMAS database for baseline and AGC specimens. All characterization data are qualified for use according to their perspective data verification reports. For the AGC irradiation campaigns, a total of 41,502 qualified physics calculation records were added for AGC-1, AGC-2, and AGC-3. Finally, a total of 173,698,505 AGC irradiation monitoring data records (thermocouple temperature, gas flow rate, gas pressure, gas moisture, applied load, and specimen displacement) have been captured to NDMAS; the majority of those records (~94%) are qualified data records.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Constraining capture cross sections using proton inelastic scattering as a surrogate reaction

The surrogate reaction method is an alternative to direct measurements of compound nuclear reaction cross sections. We introduce theory tools for extracting capture cross sections from experiments that use proton inelastic scattering as a surrogate reaction mechanism. This makes it possible to constrain compound nucleus decay models which are typically the largest source of uncertainty in capture cross section calculations. This letter describes the theory developments that were used to simultaneously infer 89 Y(p, γ) and 89 Zr(n, γ) cross sections from 90 Zr(p, p′γ) surrogate measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Computational Assessment of Actinide Dioxide Cations AnO 2 2+ for An = U to Lr: The Limited Stability Range of the Hexavalent Actinyl Moiety, [O=An=O] 2+

The isolated gas-phase actinide dioxide dications, AnO 2 2+ , were evaluated by DFT and coupled cluster CCSD(T) calculations for 12 actinides, An = U–Lr. CASSCF calculations were used to define the orbitals for the CCSD(T) calculations. The characteristic linear [O$=$An$=$O] 2+ hexavalent actinyl(VI) was found to be the lowest energy structure for An = U, Np, and Pu, which also form stable actinyl(VI) species in solution and possibly for Am when spin–orbit effects are included. For Am, there is a divalent [An II (O 2 )] 2+ structure where the dioxygen is an end-on physisorbed η 1 - 3 O 2 2 kcal/mol above the actinyl when spin–orbit effects are included which lower the energy of the actinyl structure. For An = Cm, Bk, and Lr, the lowest energy structure is trivalent [An III (O 2 – )] 2+ where the dioxygen is a side-on superoxide, η 2 -O 2 – . For Cm, the actinyl is close in energy to the ground state when spin–orbit effects are included. For An = Cf, Es, Fm, Md, and No, the lowest energy structure is divalent [An II (O 2 )] 2+ where the dioxygen is an end-on physisorbed η 1 - 3 O 2 . The relative energies suggest that curyl(VI) and berkelyl(VI), like well-known americyl(VI), might be stabilized by coordinating ligands in condensed phases. The results further indicate that for californyl and beyond, the actinyl(VI) moieties will probably be elusive even using strong donor ligands. The prevalence of low oxidation states (OSs) An(II) and An(III) for transplutonium actinides reflects stabilization of the 5f orbitals and validates established trends, including the remarkably high stability of divalent No. Finally, bond distances and other parameters suggest maximum bond covalency around plutonyl(VI), with a particularly substantial decrease in bond strength between americyl(VI) and curyl(VI).

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

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
DOE OSTI2020

Gating of Hydrophobic Nanopores with Large Anions

In this report, understanding ion transport in nanoporous materials is critical to a wide variety of energy and environmental technologies, ranging from ion-selective membranes, drug delivery, and biosensing, to ion batteries and supercapacitors. While nanoscale transport is often described by continuum models that rely on a point charge description for ions and a homogeneous dielectric medium for the solvent, here, we show that transport of aqueous solutions at a hydrophobic interface can be highly dependent on the size and hydration strength of the solvated ions. Specifically, measurements of ion current through single silicon nitride nanopores that contain a hydrophobic-hydrophilic junction show that transport properties are dependent not only on applied voltage but also the type of anion. We find that in Cl - -contained solutions the nanopores only conducted ionic current above a negative voltage threshold. On the other hand, introduction of large polarizable anions, such as Br - and I - , facilitated the pore wetting, making the pore conductive at all examined voltages. Molecular dynamics simulations revealed that the large anions, Br - and I - , have a weaker solvation shell compared to that of Cl - , and consequently were prone to migrate from the aqueous solution to the hydrophobic surface, leading to the anion accumulation responsible for pore wetting. The results are essential for designing nanoporous systems that are selective to ions of the same charge, for realization of ion-induced wetting in hydrophobic pores, as well as for fundamental understanding on the role of ion hydration shell on the properties of solid/liquid interfaces.

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

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Reactivity descriptors in acid catalysis: acid strength, proton affinity and host–guest interactions

Brønsted acids mediate chemical transformations via proton transfer to bound species and interactions between the conjugate anion and bound cationic intermediates and transition states that are also stabilized by van der Waals forces within voids of molecular dimensions in inorganic hosts. This Feature Article describes the relevant descriptors of reactivity in terms of the properties of acids and molecules that determine their ability to donate and accept protons and to reorganize their respective charges to optimize their interactions at bound states. The deprotonation energy (DPE) of the acids and the protonation energy (Eprot) of the gaseous analogs of bound intermediates and transition states reflect their respective properties as species present at non-interacting distances. These properties accurately describe the reactivity of acids of a given type, such as polyoxometalates (POM) with a given type of addenda atom but different central atoms and heterosilicates, for different families of reactions. They do not fully capture, however, differences among acid types (e.g., Mo and W POM, heterosilicates, mineral acids) for diverse types of chemical transformations (e.g., elimination, isomerization, dimerization, condensation). The incompleteness of such descriptors reflects their inability to describe how protonated molecular species and conjugate anions restructure their respective charges when present as a binding pair at interacting distances. Such interaction energies represent electrostatic forces that depend on charge distributions in the cations and anions and the ability to reorganize the distributions to maximize the interactions. In the case of deprotonation, the electrostatic and charge reorganization components of DPE for various acids solely reflect the ability of the conjugate anion to accept and distribute the negative charge, a characteristic unique of each type of solid acid and specifically of the composition of its extended conjugate anion framework. The energy required to accept and rearrange the positive charge in bound intermediates and transition states reflects, in turn, their respective ability to recover the ionic and covalent components of DPE, the energy required to detach proton from conjugate anions. The DPE components and the recovery fractions together lead to a modified DPE, which captures only the part of DPE that remains unrecovered by the ion-pair interactions at bound intermediates and transition states, as the unifying descriptor for broad families of acids and reactions. The electrostatic and charge reorganization energies involved in these general descriptors are placed in historical context by assessing their connections to the heuristics of hard–soft acid–base displacements. Further development of these concepts requires benchmarking and extension of electrostatic and reorganization components of energies for a more diverse set of reaction types and acid families and advancement of methods for more efficient calculations of electrostatic interactions. Reactivity descriptors must also account for dispersive interactions between host cavities and guest molecules, requiring a framework analogous to the one described here for ion-pair interactions; these dispersive interactions depend on the fit between their shapes and sizes as well as their ‘‘structural stiffness’’ that determines the ability to modify the shapes of molecules and voids to minimize free energy. Entropy considerations and estimates of their dependence on properties of catalysts and molecules are also required for accurately determining Gibbs free energies that ultimately determine reaction rates.

Deshlahra, Prashant↗
DOE OSTI2025

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗