Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometric modelling”

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

Attribute-Guided Adversarial Training for Robustness to Natural Perturbations

We report while existing work in robust deep learning has focused on small pixel-level norm-based perturbations, this may not account for perturbations encountered in several real world settings. In many such cases although test data might not be available, broad specifications about the types of perturbations (such as an unknown degree of rotation) may be known. We consider a setup where robustness is expected over an unseen test domain that is not i.i.d. but deviates from the training domain. While this deviation may not be exactly known, its broad characterization is specified a priori, in terms of attributes. We propose an adversarial training approach which learns to generate new samples so as to maximize exposure of the classifier to the attributes-space, without having access to the data from the test domain. Our adversarial training solves a min-max optimization problem, with the inner maximization generating adversarial perturbations, and the outer minimization finding model parameters by optimizing the loss on adversarial perturbations generated from the inner maximization. We demonstrate the applicability of our approach on three types of naturally occurring perturbations --- object-related shifts, geometric transformations, and common image corruptions. Our approach enables deep neural networks to be robust against a wide range of naturally occurring perturbations. We demonstrate the usefulness of the proposed approach by showing the robustness gains of deep neural networks trained using our adversarial training on MNIST, CIFAR-10, and a new variant of the CLEVR dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Investigation of Point-Contact Strategies for CFD Simulations of Pebble-Bed Reactor Cores

This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y + wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗
DOE OSTI2025

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗
DOE OSTI2020

Long-Range Exciton Diffusion in Two-Dimensional Assemblies of Cesium Lead Bromide Perovskite Nanocrystals

Förster resonant energy transfer (FRET)-mediated exciton diffusion through artificial nanoscale building block assemblies could be used as an optoelectronic design element to transport energy. However, so far, nanocrystal (NC) systems supported only diffusion lengths of 30 nm, which are too small to be useful in devices. In this work, we demonstrate a FRET-mediated exciton diffusion length of 200 nm with 0.5 cm 2 /s diffusivity through an ordered, two-dimensional assembly of cesium lead bromide perovskite nanocrystals (CsPbBr 3 PNCs). Exciton diffusion was directly measured via steady-state and time-resolved photoluminescence (PL) microscopy, with physical modeling providing deeper insight into the transport process. This exceptionally efficient exciton transport is facilitated by PNCs' high PL quantum yield, large absorption cross section, and high polarizability, together with minimal energetic and geometric disorder of the assembly. This FRET-mediated exciton diffusion length matches perovskites' optical absorption depth, thus enabling the design of device architectures with improved performances and providing insight into the high conversion efficiencies of PNC-based optoelectronic devices.

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

Geometry-complete diffusion for 3D molecule generation and optimization

Abstract Generative deep learning methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a denoising diffusion framework. However, such methods are unable to learn important geometric properties of 3D molecules, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which notably hinders their ability to generate valid large 3D molecules. In this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset and the larger GEOM-Drugs dataset, respectively. Importantly, we demonstrate that GCDM’s generative denoising process enables the model to generate a significant proportion of valid and energetically-stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but can be repurposed to consistently optimize the geometry and chemical composition of existing 3D molecules for molecular stability and property specificity, demonstrating new versatility of molecular diffusion models. Code and data are freely available on GitHub .

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

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The age–chemical abundance structure of the Galactic disc – II. α -dichotomy and thick disc formation

ABSTRACT We extend our previous work on the age–chemical abundance structure of the Galactic outer disc to the inner disc (4 < r < 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]–[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multiphase inner disc formation scenario that matches the observed bimodal [Mg/Fe]–[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes $6\,$Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times ($t\sim 1\,$ Gyr) and the second one $6\,$ Gyr later at a cosmic time of $t\sim 7\,$ Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-α sequence, and the second starburst leads to the formation of the metal-poor low-α sequence. The metal-rich low-α stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the α-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z >1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.

Lian, Jianhui↗
DOE OSTI2020

The TESS Phase Curve of KELT-1b Suggests a High Dayside Albedo

We measured the optical phase curve of the transiting brown dwarf KELT-1b (TOI 1476) using data from the TESS spacecraft. We found that KELT-1b shows significant phase variation in the TESS bandpass, with a relatively large phase amplitude of 234{sub −44}{sup +43} ppm and a secondary eclipse depth of 371{sub −49}{sup +47} ppm. We also measured a marginal eastward offset in the dayside hot spot of 18.°3 ± 7.°4 relative to the substellar point. We detected a strong phase-curve signal attributed to ellipsoidal distortion of the host star with an amplitude of 399 ± 19 ppm. Our results are roughly consistent with the Spitzer phase curves of KELT-1b, but the TESS eclipse depth is deeper than expected. Our cloud-free 1D models of KELT-1b’s dayside emission are unable to fit the full combined eclipse spectrum. Instead, the large TESS eclipse depth suggests that KELT-1b may have a significant dayside geometric albedo of A {sub g} ∼ 0.5 in the TESS bandpass, which would agree with the tentative trend between equilibrium temperature and geometric albedo recently suggested by Wong et al. We posit that if KELT-1b has a high dayside albedo, it is likely due to silicate clouds that form on KELT-1b’s nightside and are subsequently transported onto the western side of KELT-1b’s dayside hemisphere before breaking up.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermophysical Modeling of 20 Themis Family Asteroids with WISE/NEOWISE Observations

Themis family is one of the largest and oldest asteroid populations in the main belt. Water ice may widely exist on the parent body (24) Themis. In this work, we employ the Advanced Thermophysical Model as well as midinfrared measurements from NASA’s Wide-field Infrared Survey Explorer to explore thermal parameters of 20 Themis family members. Here we show that the average thermal inertia and geometric albedo are 39.5 ± 26.0 J m{sup −2} s{sup −1/2} K{sup −1} and 0.067 ± 0.018, respectively. The family members have a relatively moderate roughness fraction on their surfaces. We find that the relatively low albedos of Themis members are consistent with the typical values of B-type and C-type asteroids. As mentioned, the Themis family bears a very low thermal inertia, which indicates a fine and mature regolith on their surfaces. The resemblance of thermal inertia and geometric albedo of Themis members may reveal their close connection in origin and evolution. In addition, we present the compared results of thermal parameters for several prominent families.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Aerial Data and Processed Models of Port Arthur Coastal Neighborhood and Pleasure Island Golf Course, June 2024

Our Co-design team is from the University of Texas, working on a Department of Energy-funded project focused on the Beaumont-Port Arthur area. As part of this project, we will be developing climate-resilient design solutions for areas of the region. More on www.caee.utexas.edu.We captured aerial photos in the Port Arthur Coastal Neighborhood Community and the Golf Course on Pleasure Island, Texas, in June 2024.Aerial photos taken were through DroneDeploy autonomous flight, and models were processed through the DroneDeploy engine as well. All aerial photos are in .JPG format and contained in zipped files for each area.The processed data package includes 3D models, geospatial data, mappings, and point clouds. Please be aware that DTM, Elevation toolbox, Point cloud, and Orthomosaic use EPSG: 6588. And 3D Model uses EPSG: 3857.For using these data:- The Adobe Suite gives you great software to open .Tif files.- You can use LASUtility (Windows), ESRI ArcGIS Pro (Windows), or Blaze3D (Windows, Linux) to open a LAS file and view the data it contains.- Open an .OBJ file with a large number of free and commercial applications. Some examples include Microsoft 3D Builder, Apple Preview, Blender, and Autodesk.- You may use ArcGIS, Merkaartor, Blender (with the Google Earth Importer plug-in), Global Mapper, and Marble to open .KML files.- The .tfw world file is a text file used to georeference the GeoTIFF raster images, like the orthomosaic and the DSM. You need suitable software like ArcView to open a .TFW file.This dataset provides researchers with sufficient geometric data and the status quo of the land surface at the locations mentioned above. This dataset could streamline researchers' decision-making processes and enhance the design as well.

2D mapping↗
DOE OSTI2020

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

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

Data-driven Whitney forms for structure-preserving control volume analysis

Control volume analysis models physics via the exchange of generalized fluxes between subdomains. Here, we introduce a scientific machine learning framework adopting a partition of unity architecture to identify physically-relevant control volumes, with generalized fluxes between subdomains encoded via Whitney forms. The approach provides a differentiable parameterization of geometry which may be trained in an end-to-end fashion to extract reduced models from full field data while exactly preserving physics. The architecture admits a data-driven finite element exterior calculus allowing discovery of mixed finite element spaces with closed form quadrature rules. An equivalence between Whitney forms and graph networks reveals that the geometric problem of control volume learning is equivalent to an unsupervised graph discovery problem. The framework is developed for manifolds in arbitrary dimension, with examples provided for H(div) problems in $\mathbb{R}$ establishing convergence and structure preservation properties. Finally, we consider a lithium-ion battery problem where we discover a reduced finite element space encoding transport pathways from high-fidelity microstructure resolved simulations. The approach reduces the 5.89M finite element simulation to 136 elements while reproducing pressure to under 0.1% error and preserving conservation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Renormalizing two-fermion operators in the SMEFT via supergeometry

We extend the geometric framework of field-space covariance for loop computations, thereby unifying the treatment of scalars, fermions, and gauge bosons in effective field theories. This allows us to derive a manifestly covariant formula for one-loop UV divergences that includes contributions from mixed boson-fermion graphs. The result is expressed in terms of geometric invariants of the field-space supermanifold. As a demonstration of this formula, we compute the renormalization group equations for two-fermion operators at the dimension-eight level in the Standard Model Effective Field Theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗