Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface coupling”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2026

Influence of Surface Treatment on the Interfacial and Mechanical Performance of Metal-Insert Overmolded Composites

Structural composites with metal inserts are gaining increasing interest enabling light weight design for energy efficiency, provide insulation, corrosion protection without compromising mechanical integrity. Such over molded metal-polymer composites show strong potential in many automotive and aerospace component parts, e.g. car bumper, chassis, door, brackets, fasteners, etc. Despite several advantages, metal-polymer dissimilar interface remains a weak spot in structural design, requiring further investigation to improve interfacial interlocking. The current work focuses on enhancing the metal-polymer adhesion between stainless steel inserts and polymer matrices (e.g. acrylonitrile butadiene styrene resin (ABS), and polyamide 66 (PA66) through surface treatment via chemical etching and mechanical roughening. In the chemical etching process, metal substrates were treated with aqua regia solution for different time intervals to evaluate surface roughness measured by water contact angle and morphology analysis. For mechanical roughening approach, metal inserts were subjected to blasting for surface abrasion. Treated samples were also functionalized with aminosilane coupling agent which showed successful grafting on the surface measured by water contact angle. Metal inserts with micro level surface roughness through etching exhibited significant 21% improvement in adhesion strength compared to untreated inserts; however, the effect of silane treatment on the interfacial strength was found to be highly matrix dependent.

Saha, Subhabrata [ORNL] (ORCID:000000033966064X)↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Behind the Meter Storage for Electric Vehicle Charging, Electrochemical and Thermal Energy Storage, and Solar Photovoltaic

In response to the potentially large and irregular demand from EVs, along with changing load profiles from buildings with on-site generation, utilities are evaluating multiple options for managing dynamic loads, including time-of-use pricing, demand charges, battery storage, and curtailment of variable generation. Buildings, as well as commercial, public, and workplace EV charging operations, can use a combination of electrochemical battery storage and thermal energy storage coupled with on-site generation to manage energy costs as well as provide resiliency and reliability for EV charging and building energy loads. We are completing a behind the meter storage analysis that focuses on determining the optimal system designs and energy flows for thermal and electrochemical behind the meter storage with on-site solar photovoltaic (PV) generation enabling electric vehicle charging in various climates, building types, and utility rate structures. In completing this analysis, we have developed a tool that combines existing battery models via the System Advisor Model (SAM) and building modeling software via EnergyPlus into a single interface. This tool allows us to simulate a building with a detailed battery model to properly size the battery, thermal energy storage, and solar PV systems to maximize profit for the system owner. This also allows us to assess how the battery degrades under various supervisory control dispatch algorithms to control charging/discharging; we can also see how thermal energy storage is created and used to complement the battery to reduce thermal loads in the building. With this project, we can analyze new batteries that are designed specifically for energy storage, rather than designed to be extremely energy dense for electric vehicle applications, using battery lifetime models from other national labs and the existing SAM battery model, which has detailed lifetime and degradation parameters. We can also assess novel thermal storage technologies by integrating them into the whole building energy simulation program EnergyPlus. Because the model calls both SAM and EnergyPlus, required inputs need to be compatible for both models. These inputs include, on a high-level, the following: weather files, building and electric vehicle load profiles, electricity rate tariff information, and system cost information for the stationary battery, solar PV, and thermal storage system. The various buildings we are studying for this analysis are retail big-box grocery store, commercial office building, fleet vehicle depot and operations facility, multi-family residential, and electric vehicle charging station. For these different applications, the battery and thermal storage will be dispatched differently, and the various technologies are sized differently to optimize cost.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Micro-tensile characteristics of As-fabricated and irradiated AGR-2 TRISO fuel particle buffer, IPyC, and buffer-IPyC interlayer regions

A recently developed micro-tensile sample preparation technique was implemented to evaluate the tensile strengths of the buffer, IPyC, and buffer-IPyC interlayer regions of the unirradiated and irradiated AGR-2 TRISO fuel particles. Understanding the mechanical properties of the buffer-IPyC interlayer is essential for developing thermomechanical models of buffer-IPyC separation, yet there is a lack of experimental data on its micro-tensile properties. TEM analysis was conducted on these regions to determine the microstructural changes relevant to the samples' tensile properties. In the unirradiated TRISO particle samples, the buffer layer demonstrated the weakest tensile strength, while the IPyC layer exhibited the highest. Conversely, in the irradiated TRISO particle samples, the buffer-IPyC interlayer region showed the lowest tensile strength, with the IPyC layer being the strongest. Fractures in the samples from the buffer-IPyC region predominantly occurred either in the buffer layer or at the buffer-IPyC interface. However, some buffer-IPyC interlayer samples displayed stress-strain and fracture behaviors more akin to the IPyC layer than the buffer layer. Analysis of diffraction patterns suggests that irradiation may have increased anisotropy in the three regions tested. Despite this suggested increase in anisotropy, there was no evidence that it affected the measured strengths. The irradiated TRISO particles demonstrated a considerable increase in void space and a decrease in ultimate tensile strength within the buffer-IPyC interlayer region due to the densification and contraction of the buffer layer. Minor variations in diffraction ring patterns were also observed. These changes, coupled with a significant reduction in the Weibull modulus/shape parameter, imply that irradiation-induced densification leads to tearing between the buffer and IPyC layers at locations of elevated porosity in the buffer-IPyC interlayer region.

Tristructural isotropic (TRISO)↗
DOE OSTI2020

Extraction and validation of a new set of CMS pythia8 tunes from underlying-event measurements

New sets of CMS underlying-event parameters ("tunes") are presented for the PYTHIA8 event generator. These tunes use the NNPDF3.1 parton distribution functions (PDFs) at leading (LO), next-to-leading (NLO), or next-to-next-to-leading (NNLO) orders in perturbative quantum chromodynamics, and the strong coupling evolution at LO or NLO. Measurements of charged-particle multiplicity and transverse momentum densities at various hadron collision energies are fit simultaneously to determine the parameters of the tunes. Comparisons of the predictions of the new tunes are provided for observables sensitive to the event shapes at LEP, global underlying event, soft multiparton interactions, and double-parton scattering contributions. In addition, comparisons are made for observables measured in various specific processes, such as multijet, Drell-Yan, and top quark-antiquark pair production including jet substructure observables. The simulation of the underlying event provided by the new tunes is interfaced to a higher-order matrix-element calculation. For the first time, predictions from PYTHIA8 obtained with tunes based on NLO or NNLO PDFs are shown to reliably describe minimum-bias and underlying-event data with a similar level of agreement to predictions from tunes using LO PDF sets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Total metals, sulfur and organic carbon data; Slate River floodplain, Crested Butte, CO; March 2021-October 2021

This data package includes processed and undiluted measurements for metal, sulfur and organic carbon concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between March and October of 2021. These measurements were all recorded at the Environmental Measurements Facility (EM-1) at Stanford University in Stanford, CA. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 2C until measured at EM-1.Analysis by ICP-MS:Measurements for total metals were performed on an inductively coupled plasma optical emission spectrometer (ICP-OES; iCAP 6300, Thermo Scientific, Cambridge, U.K.). Calibration standards were prepared from the mulit element stock solution (Sigma-Aldrich Multielement standard for ICP, St. Louis, MO) using matrix matched to sample solutions (2% HNO3). Calibration curves included 5 points with correlation coefficients of >0.99. The QC protocol includes a continuing calibration blank and quality control samples that are analyzed just after calibration and again every 20 samples and at the completion of the run. Acceptable QC responses must be between 90-110% of the certified value. Analysis by Total Organic Carbon:Dissolved organic carbon concentrations were quantified on a total organic carbon (TOC) analyzer (TOC-L, Shimadzu, Kyoto, Japan) running the NPOC method. Standard curves were developed using an Organic Carbon Standard from RICCA Chemical Company (Arlington,TX). A series of 4 to 6 standards were automatically diluted by the instrument in a concentration range that spans that of the samples. A blank sample was run just after the calibration curve and at the end of the run. A QC sample was run every 25 samples. Acceptable QC responses must be between 90-110% of the certified value.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measurement of Photons Emitted by High-Energy Charged Particles as Background in Single-Photon Resolving Image Sensors

This work introduces an advanced technique optimized for detecting photons generated by charged particles, leveraging Skipper charge coupled device (Skipper-CCD) image sensors. By analyzing background sources and detection efficiencies, the technique achieves strong agreement between experimental results and Cherenkov-based simulations. It also provides a robust framework for investigating secondary photon production in environments with high fluxes of ionizing particles, such as those anticipated in space-based astronomical instruments. These secondary photons present a critical challenge as background noise for next-generation single-photon resolving imagers used to study faint celestial objects. Furthermore, the method exhibits significant potential for broader applications, including exploring photon generation in various substrate materials and examining their transport through multiple interfaces.

Fernandez Moroni, Guillermo [Fermilab; Chicago U.,↗
DOE OSTI2022

Integrated Optimization and Control of a Hybrid Gas Turbine/sCO 2 Power System

During phase-I, the project team led by Echogen Power Systems (EPS) had two primary objectives based on investigating the application of gas turbines with supercritical carbon dioxide (sCO 2 ) power cycles. The first objective was to improve the overall efficiency and performance of a hybrid gas turbine/sCO 2 power system through a joint optimization of the two subsystems (gas turbine and sCO 2 power cycle) using non-linear optimization techniques that simultaneously evaluate thermal performance of the combined cycle. The hybrid power system included several points of interaction, including (but not limited to) gas turbine exhaust, fuel heating, inlet chilling and turbine cooling. The second objective was to establish a baseline transient response model of the hybrid power system and a notional microgrid and begin steps to integrate the control systems of the three major elements (gas turbine, sCO 2 cycle and grid controller). The project team established a baseline performance for a combined cycle power plant using a production gas turbine and scaled sCO 2 power cycle only utilizing exhaust heat recovery. Echogen’s non-linear techno-economic optimization code was extended by adding gas turbine component models derived from a in-house developed gas turbine design code. With the two cycles coupled by the gas turbine exhaust, design parameters of both cycles were allowed to vary simultaneously to determine performance opportunity versus isolated designs. Returning to the baseline gas turbine/sCO 2 power cycle transient models: Echogen had in-house developed sCO 2 cycle transient model in GT-Suite system simulation software, and had partnered with Siemens Finspång for gas turbine transient model, and Siemens PTI group to provide micro-grid load profile as well as hybrid power cycle generated load (power and frequency) analysis. The transient model for the SGT-750 Siemens gas turbine was a “black-box” functional mock-up interface (FMI) model developed by Siemens Industrial Turbomachinery in Finspång, Sweden. The SGT-750 is a twin-shaft gas turbine that produces 40 MW electricity with an efficiency of about 40% at ISO conditions. At 100% gas turbine throttle (load), the SGT-750 has average exhaust conditions of 114.6 kg/s and 469.8°C. The transient model for sCO 2 power cycle was developed by Echogen in GT-SUITE 1D system simulation software platform. The basic CO 2 flow circuit has single-shaft turbomachinery with net 11.5 MW electrical power output at design conditions. The power turbine has a double-ended shaft with one end connected to synchronous generator through a fixed-ratio gearbox. The other end of power turbine is connected to the compressor through a continuously variable transmission. The major components of the sCO 2 power cycle modeled include air cooled condenser/cooler, CO 2 compressor, recuperator, two waste heat exchanger coils, power turbine, continuous variable transmission, gearbox and generator. Integration of SGT-750 transient model and sCO 2 power cycle transient model was done in Matlab Simulink. In the integrated model, the gas turbine and sCO 2 power cycle interacted at two points, first one being the gas turbine exhaust gas flow rate and temperature, which were inputs to sCO 2 power cycle model. The second point was the distribution of micro-grid load demand signal between the SGT-750 generator and sCO 2 cycle generator. For a given combined-cycle load demand, the gas turbine load demand was equal to the total demand minus the sCO 2 cycle power generated. In the present study the integrated model was simulated for two cases of grid load demand: (i) for a step change, both positive-step and negative-step, in grid load demand (ii) for a micro-grid load demand curve provided by Siemens PTI group. Finally, the time series plots representing load demand versus integrated system response were presented including the sCO 2 power cycle control system performance plots. The actual generated power and frequency of both the generators, gas turbine and sCO 2 power cycle, was supplied to Siemens PTI group for dynamic grid assessment, results of which are provided in appendices.

03 NATURAL GAS↗
DOE OSTI2023

Fluorinated High-Voltage Electrolytes To Stabilize Nickel-Rich Lithium Batteries

As state-of-the-art (SOA) lithium-ion (Li-ion) batteries approach their specific energy limit (~250 Wh kg -1 ), layer-structured, nickel-rich (Ni-rich) lithium transition metal oxide-based cathode materials, e. g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), have attracted great interest owing to their practical high specific capacities (>200 mAhg -1 ). Coupled with their high average discharge voltages (~4 V vs Li/Li + ), Ni-rich cathode-based lithium batteries possess a great potential to achieve much higher specific energies (>350 Wh kg -1 at the cell level) than the SOA Li-ion counterparts. In addition, Ni-rich oxides are low-cost battery cathode materials due to their low cobalt contents. However, Ni-rich cathode-based lithium batteries suffer quick capacity degradations upon cycling, particularly at high upper cutoff voltages (e.g., ≥= 4.5 V vs Li/Li + ), due to crystal structure changes of the active cathode materials and parasitic side reactions at the electrolyte/electrode interfaces. In this study, a fluorinated-solvent-based, high-voltage stable electrolyte (HVE), i.e., 1 M Li bis(trifluoromethanesulfonyl)imide (LiTFSI) in fluoroethylene carbonate (FEC), bis(2,2,2-trifluoroethyl) carbonate (FDEC), and methyl (2,2,2-trifluoroethyl) carbonate (FEMC) with Li difluoro(oxalate)borate (LiDFOB) additive, was formulated and evaluated in Li/NMC811 battery cells. Further, to the best of our knowledge, this class of electrolyte has not been investigated for Ni-rich cathode-based lithium batteries. Li/NMC811 cells with HVE exhibited a superior long-term cycle performance stability, maintaining ~80% capacity after ~500 cycles at a high cutoff voltage of 4.5 V (vs Li/Li + ) than a baseline carbonate-solvent-based electrolyte (BE). The superior cycle stability of the Li/NMC811 cells is attributed to the inherently high-voltage stability of HVE, supported by the physical and electrochemical analyses. This conclusion is supported by our density functional theory (DFT) modeling where HVE shows a less tendency of deprotonation/oxidation than BE, leading to the observed cycle stability. The findings in this study are important to help tackle the technical challenges facing Ni-rich cathode-based lithium batteries to realize their high energy density potentials with a long cycle life.

25 ENERGY STORAGE↗
DOE OSTI2021

A New Modeling Framework for Multi-Scale Simulation of Hydraulic Fracturing and Production from Unconventional Reservoirs

This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.

58 GEOSCIENCES↗
DOE OSTI2022

Emergence of electric-field-tunable interfacial ferromagnetism in 2D antiferromagnet heterostructures

Van der Waals (vdW) magnet heterostructures have emerged as new platforms to explore exotic magnetic orders and quantum phenomena. Here, we study heterostructures of layered antiferromagnets, CrI 3 and CrCl 3 , with perpendicular and in-plane magnetic anisotropy, respectively. Using magneto-optical Kerr effect microscopy, we demonstrate out-of-plane magnetic order in the CrCl 3 layer proximal to CrI 3 , with ferromagnetic interfacial coupling between the two. Such an interlayer exchange field leads to higher critical temperature than that of either CrI 3 or CrCl 3 alone. We further demonstrate significant electric-field control of the coercivity, attributed to the naturally broken structural inversion symmetry of the heterostructure allowing unprecedented direct coupling between electric field and interfacial magnetism. These findings illustrate the opportunity to explore exotic magnetic phases and engineer spintronic devices in vdW heterostructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Characterization of high-flux cadmium zinc telluride detectors coupled to MM-PAD readout ASICs for high-flux X-ray applications

High-flux cadmium zinc telluride (HF-CZT) is a promising material for next-generation X-ray imaging due to its high stopping power, wide bandgap, and improved resistance to radiation-induced polarization [1]. We bonded 2 mm thick HF-CZT sensors with Cornell's Mixed-Mode Pixel Array ASICs, which combine analog charge integration and in-pixel digital counting to achieve single X-ray sensitivity, large full well capacity, and frame rates up to 10 kHz [2,3]. Detector characterization showed stable electron collection above 150 V bias, while hole collection remained limited by trapping even at 700 V. Temporal studies revealed non-negligible afterglow, with an unexpected increase in afterglow at higher biases, which was reduced by infrared illumination, suggesting that injected carriers reduce trapping near contact interfaces. Polarization measurements confirmed improved stability compared to conventional CZT, though polarization effects persisted, particularly with hole collection, with up to 5% signal loss and spatial distortion at high doses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unconventional spin-orbit torques from sputtered MoTe 2 films

Materials with strong spin-orbit coupling and low crystalline symmetry are promising for generating large unconventional spin-orbit torques (SOTs), such as in-plane fieldlike (FL) torques and out-of-plane dampinglike (DL) torques, which can effectively manipulate and deterministically switch an out-of-plane magnetization without the need for additional external in-plane magnetic fields. Here, we report SOTs generated by magnetron-sputtered 1⁢T' MoTe 2 /Permalloy (Py; Ni 80 ⁢Fe 20 )/MgO heterostructures using both spin-torque ferromagnetic resonance (ST-FMR) and second harmonic Hall measurements. We observed unconventional FL and DL torques in our samples due to spins polarized normal to the interface of MoTe 2 and Py layers, and studied the influence of crystallographic order and MoTe 2 layer thickness on the SOTs. By comparing the Raman spectra of 1⁢T' MoTe 2 samples prepared in different ways, we found a tensile strain in sputtered MoTe 2 films, which might further enhance the generation of unconventional torques by reducing the symmetry of 1⁢T' MoTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Controlling single rare earth ion emission in an electro-optical nanocavity

Abstract Rare earth emitters enable critical quantum resources including spin qubits, single photon sources, and quantum memories. Yet, probing of single ions remains challenging due to low emission rate of their intra-4 f optical transitions. One feasible approach is through Purcell-enhanced emission in optical cavities. The ability to modulate cavity-ion coupling in real-time will further elevate the capacity of such systems. Here, we demonstrate direct control of single ion emission by embedding erbium dopants in an electro-optically active photonic crystal cavity patterned from thin-film lithium niobate. Purcell factor over 170 enables single ion detection, which is verified by second-order autocorrelation measurement. Dynamic control of emission rate is realized by leveraging electro-optic tuning of resonance frequency. Using this feature, storage, and retrieval of single ion excitation is further demonstrated, without perturbing the emission characteristics. These results promise new opportunities for controllable single-photon sources and efficient spin-photon interfaces.

42 ENGINEERING↗
DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗