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 757 records · Page 42

Higher-derivative corrections to entropy and the weak gravity conjecture in Anti-de Sitter space

We compute the four-derivative corrections to the geometry, extremality bound, and thermodynamic quantities of AdS-Reissner-Nordström black holes for general dimensions and horizon geometries. We confirm the universal relationship between the extremality shift at fixed charge and the shift of the microcanonical entropy, and discuss the consequences of this relation for the Weak Gravity Conjecture in AdS. The thermodynamic corrections are calculated using two different methods: first by explicitly solving the higher-derivative equations of motion and second, by evaluating the higher-derivative Euclidean on-shell action on the leading-order solution. In both cases we find agreement, up to the addition of a Casimir energy in odd dimensions. We derive the bounds on the four-derivative Wilson coefficients implied by the conjectured positivity of the leading corrections to the microcanonical entropy of thermodynamically stable black holes. These include the requirement that the coefficient of Riemann-squared is positive, meaning that the positivity of the entropy shift is related to the condition that c – a is positive in the dual CFT. We discuss implications for the deviation of η/s from its universal value and a potential lower bound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fundamental Issues Identified for Thermodynamic Description of Molten Salt Systems

Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.

42 ENGINEERING↗

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: I. Fully vibrationally-resolved electronic excitation of H 2 $X^1Σ^+_g)$

Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of molecular hydrogen suitable for implementation in collisional-radiative models. The adiabatic-nuclei molecular convergent close-coupling method is used to calculate cross sections for excitation of all bound vibrational levels and dissociative excitation of the B 1 Σ u + , C 1 Π u , E F 1 Σ g + , B ′ 1 Σ u + , G K 1 Σ g + , I 1 Π g , J 1 Δ g , D 1 Π u , H 1 Σ g + , b 3 Σ u + , c 3 Π u , a 3 Σ g + , e 3 Σ u + , d 3 Π u , h 3 Σ g + , g 3 Σ g + , i 3 Π g , and j 3 Δ g electronic states from all –14 bound vibrational levels of the ground electronic ($\text{X}$ $^1&#x3A3^+_g$) state. The data set consists of cross sections from threshold to 500 eV for over 5000 transitions, representing all possible electronic and vibrational transitions between the state and the –3 singlet and triplet states (where refers to the united-atoms-limit principle quantum number). The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interaction of extended dislocations with nanovoid clusters

Voids of nanoscale dimensions in irradiated metals can act as obstacles to dislocation motion and cause strengthening. In this work, nanovoid strengthening and the influences of void size, void spacing and material properties, such as stacking fault energies, on dislocation bypass mechanisms are investigated using Phase Field Dislocation Dynamics, a three-dimensional mesoscale model that predicts the minimum energy pathway taken by discrete dislocations. A broad range of face centered cubic metals (copper, nickel, silver, rhodium, and platinum) and nanovoid sizes and spacings are treated, altogether spanning void size–to–dislocation stacking fault width ratios from less than unity to ten. Material γ-surfaces, calculated from ab initio methods, are input directly into the formulation. The analysis reveals that the critical bypass stress scales linearly with the linear void fraction, effective isotropic shear modulus, and ratio of the intrinsic to unstable stacking fault energies. With only a few exceptions, the critical stress is controlled by the stress required for the leading partial to impinge the voids (to move within range of the attractive image stress field of the void). When the void diameter is nearly an order of magnitude greater than the stacking fault width, the mechanism determining critical strength shifts to the stress for the dislocation to breakaway after partially cutting the void. Furthermore, this situation corresponds to that treated by line tension models and is realized here for Pt, with a sub-nanometer stacking fault width.

36 MATERIALS SCIENCE↗

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗

Radiation hot spot formation on UF 6 30B cylinders

The accurate material accountancy of a country's enriched uranium hexafluoride (UF 6 ) stockpiles is imperative to the International Atomic Energy Agency (IAEA) for determining appropriate material safeguards. Uranium hexafluoride is typically stored in large cylinders with enrichment verification typically performed with a NaI gamma detector using the enrichment meter method. This approach for calculating enrichment assumes that the UF 6 is homogeneous inside the cylinder. Here, the measurements and tests presented in this manuscript show that the UF 6 in cylinders left outside in the elements can experience fractionation, leading to inhomogeneity, and subsequent development of radiation hot spots. These hot spot locations can produce errant enrichment measurements. In dark-colored cylinders that experience significant surface heating from the sun, UF 6 sublimates off the cylinder walls, leaving only the uranium daughter products in a horizontal patch (stripe) along the cylinder's midline. This daughter patch is hypothesized to be the source of the radiation hot spots. The patch forms during the summer and tends to decay away during the winter, when solar intensity and temperatures are reduced.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Blob velocities and sizes in the Alcator C-Mod scrape-off layer for ohmic and high confinement mode plasmas

An improved time delay estimation method is used to calculate the velocity of cross-field blob motion in the scrape-off layer of Alcator C-Mod for an ohmic and two high confinement (H-mode) plasmas; an edge localized mode free and an enhanced D-alpha H-mode. The gas puff imaging data analysis results are interpreted in the framework of a stochastic model that describes the fluctuations as a super-position of uncorrelated blob-like structures. In all confinement modes investigated, the scrape-off layer is dominated by large amplitude, blob-like filaments moving radially outwards with velocities in the range from 400 to 1000 m/s. Blobs in high confinement mode plasmas have similar velocities and sizes as in ohmic plasma, which is consistent with the close similarity of conditionally averaged burst shapes and frequency spectra for the confinement modes investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First measurement of the $Λ–Ξ$ interaction in proton–proton collisions at the LHC

The first experimental information on the strong interaction between $Λ$ and $Ξ$ - strange baryons is presented in this Letter. The correlation function of $Λ–Ξ$ - and $\overline{Λ}$ –$\overline{Ξ}$ + pairs produced in high-multiplicity proton–proton (pp) collisions at $\sqrt{s}$=13 TeV at the LHC is measured as a function of the relative momentum of the pair. The femtoscopy method is used to calculate the correlation function, which is then compared with theoretical expectations obtained using a meson exchange model, chiral effective field theory, and Lattice QCD calculations close to the physical point. Data support predictions of small scattering parameters while discarding versions with large ones, thus suggesting a weak $Λ–Ξ$ - interaction. The limited statistical significance of the data does not yet allow one to constrain the effects of coupled channels like $Σ–Ξ$ and $N–Ω$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗

Deep Learning enabled spectral energy conversion for in situ exposure measurements

A detector-specific deep learning (DL) approach is presented for spectra-to-exposure conversion using large-format sodium iodide (NaI(Tl)) detectors deployed for in situ environmental radiation measurements in emergency response scenarios. Accurate determination of exposure from NaI spectra is challenging due to poor energy resolution, partial energy absorption, and the strong sensitivity of traditionally deployed analytical conversion methods to calibrated source geometry and pre-deployment assumptions. Here, to address these limitations, a multi-layer perceptron model was trained on a hybrid in situ /Monte Carlo dataset constructed to span a broad range of photon energies, spatial extents, and realistic deployment variability, representative of general in situ emergency response conditions. The DL model was evaluated against commonly fielded analytical approaches under matched simulation conditions, including a single-factor method, a G-function method, and a modeled pressurized ion chamber (PIC) baseline. This study was intentionally computational in scope to enable controlled, like-for-like comparisons between conversion techniques while minimizing confounding real-world variability. Comparison to the modeled PIC provides contextual benchmarking and is not intended as a field inter-comparison with deployed instruments. Across the evaluated 20 keV to 3 MeV energy range, the DL approach consistently exhibited higher accuracy and reduced variance relative to the analytical methods against a deterministically calculated exposure. This may indicate improved robustness to spectral complexity without reliance on source-, geometric-, or spectral region-specific optimization. While results do not represent real-world validation, the presented work demonstrates that deep learning may effectively learn the nonlinear detector response-to-exposure relationship for asymmetric NaI(Tl) detectors and offers a promising pathway for improving in situ exposure estimation using spectroscopic systems already integrated into initial real-time emergency response operations.

61 RADIATION PROTECTION AND DOSIMETRY↗

A first principles approach for determining solute effect on electrical resistivity of Aluminum solid solution

There are no general approaches for predicting the effect of solutes on electrical resistivity changes in aluminum (Al) solid solutions. Here, we introduce a first principles-based approach for predicting temperature-dependent electrical resistivity of Al solid solutions. A combination of electronic structure methods was applied to calculate electronic scattering effects on chemical disorder and temperature induced atomic vibrations or phonons. Our calculations are validated by contrasting experimental resistivity measurements for two elements in Al solid solution: at 0.3 at. % concentration, Sn (an sp-valence solute) increases the residual resistivity by 0.2 μΩ·m, whereas Zr (a d-valence solute) causes a 2.5 μΩ·m increase. This trend may deviate for solute elements in Al that form semiconductors in pure state, like Si. Our findings provide an understanding of the role of valence electron character of solute atoms in determining electrical resistivity of Al alloys.

Samolyuk, German [ORNL] (ORCID:0000000168778255)↗

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗

Site Suitability and Air Pollution Impacts of Composting Infrastructure for California’s Organic Waste Diversion Law

The emerging field of polaritonic chemistry explores the behavior of molecules under strong coupling with cavity modes. Despite recent developments in ab initio polaritonic methods for simulating polaritonic chemistry under electronic strong coupling, their capabilities are limited, especially in cases where the molecule also features strong electronic correlation. To bridge this gap, we have developed a novel method for cavity QED calculations utilizing the Density Matrix Renormalization Group (DMRG) algorithm in conjunction with the Pauli–Fierz Hamiltonian. Our approach is applied to investigate the effect of the cavity on the S0–S1 transition of n-oligoacenes, with n ranging from 2 to 5, encompassing 22 fully correlated π orbitals in the largest pentacene molecule. Our findings indicate that the influence of the cavity intensifies with larger acenes. Additionally, we demonstrate that, unlike the full determinantal representation, DMRG efficiently optimizes and eliminates excess photonic degrees of freedom, resulting in an asymptotically constant computational cost as the photonic basis increases.

54 ENVIRONMENTAL SCIENCES↗

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal Reversible Hydrogen Potential for Electrocatalytic Ammonia Splitting Reactions in Nonaqueous Solvents from Unified pH Measurements

In this work, we introduce a new approach of using differential potentiometric measurements in four nonaqueous solvents─MeCN, THF, DMF, and PC─to determine the universal pH abs H 2 O values aligned to the aqueous pH scale for dilute NH 4 + /NH 3 solutions. Knowledge of the pH abs H 2 O values allows simple determination of the reversible hydrogen potential in any given solvent relative to the aqueous standard hydrogen electrode (SHE) and, most importantly, ensures comparability across different solvents. As an independent method, Open Circuit Potenial measurements were carried out in the same solvents titrated with NH 4 + /NH 3 to obtain alternative values for the reversible hydrogen potential in these solvents. The close agreement of these two methods, as well as calculated potentials from literature values when available, substantiates the new, simpler, and more robust approach to determine the reversible hydrogen potential introduced here. We further use the reversible hydrogen potential values established here to report, for the first time, the overpotential for ammonia oxidation as a function of solvent, with a recently discovered ruthenium catalyst.

ammonia↗

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗