Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computing”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗
DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of the temperature-dependent failure processes in PVD Cr-coated ZIRLO nuclear fuel cladding using in situ X-ray micro-tomography imaging

In this work, an accident tolerant fuel cladding system with Cr coating deposited using physical vapour deposition (PVD) method onto commercial Optimized ZIRLO™ was studied. The cladding tubes were machined into C-rings and loaded to failure under compression with real-time synchrotron micro-X-ray computed tomography (XCT) imaging at room temperature (RT), 345 °C, 650 °C and 950 °C in argon (Ar) atmosphere. The mechanical behaviour and failure processes were found strongly temperature-dependent where the Cr coating showed brittle fracture at RT and 345 °C, ductile fracture at 650 °C and a reversion to brittle fracture at 950 °C. Nanoindentation measurements and scanning electron microscopy (SEM) imaging were conducted on the materials after high temperature testing. It was found that recrystallisation of the Cr coating occurred at elevated temperatures, which significantly affected its local properties hence the failure behaviour at different temperatures. This work represents the first in situ 3D XCT observation of progressive failure processes in PVD Cr-coated ZIRLO claddings up to 950 °C providing critical insights into its brittle-to-ductile transition (BDTT) behaviour and subsequent ductile-to-brittle reversion with the increase in temperature. Moreover, results are compared with other PVD Cr-coated Zircaloy materials from open literature; the influences of the coating microstructure and local properties on the failure stress/strain and fracture processes are discussed.

Yuan, Guanjie↗
DOE OSTI2023

Novel Structural Motif To Promote Mg-Ion Mobility: Investigating ABO 4 Zircons as Magnesium Intercalation Cathodes

There is an increasing need for sustainable energy storage solutions as fossil fuels are replaced by renewable energy sources. Multivalent batteries, specifically Mg batteries, are one energy storage technology that researchers continue to develop with hopes to surpass the performance of Li-ion batteries. However, the limited energy density and transport properties of Mg cathodes remain critical challenges preventing the realization of high-performance multivalent batteries. In this work, ABO 4 zircon materials (A = Y, Eu and B = V, Cr) are computationally and experimentally evaluated as Mg intercalation cathodes. Remarkably good Mg-ion transport properties were predicted and Mg-ion intercalation was experimentally verified in sol-gel synthesized zircon YVO 4 , EuVO 4 , and EuCrO 4 . Among them, EuVO 4 exhibited the best electrochemical performance and demonstrated repeated reversible cycling. While we believe that the one-dimensional diffusion channels and redox-active species tetragonal coordination limit the value of many zircons as high-performance cathodes, their unique structural motif of overlapping polyhedra along the diffusion pathway appears instrumental for promoting good Mg-ion mobility. The motif results in a favorable "6-5-4" change in coordination that avoids unfavorable sites with lower coordination along the diffusion pathway and a structural design metric for future Mg cathode development.

25 ENERGY STORAGE↗
DOE OSTI2026

Hot-Carrier Injection and Millisecond Charge Separation from a Robust Heteroleptic Iron(II) Chromophore Immobilized on TiO 2

The synthesis, spectroscopic characterization, computational analysis, and photoelectrochemical behavior of a new iron-based chromophore, [(Cpy) 2 Fe- (deeb)](PF 6 ) 2 (Fe(Cpy) 2 (deeb)), where Cpy is 1-methyl-3-(2-pyridyl)imidazole and deeb is 4,4’-(CO 2 CH 2 CH 3 ) 2 -2,2’-bipyridine, is reported. Electrochemically reversible waves assigned to a metal-centered E o (Fe III/II ) = +0.48 and a ligand-centered E o (Fe 2+/+ ) = −1.47 V vs Fc +/0 reduction were evident in cyclic voltammetry measurements. The combination of a strong σ-donor and a π-acceptor lowered the energy of the metal-toligand charge-transfer (MLCT) excited state relative to the metal-centered state. Two MLCT transitions appear in the visible region at 424 and 580 nm. TDDFT calculations revealed that the lower-energy band was well formulated as Fe(II)→deeb, and the higher-energy transition was charge transfer to both the deeb and Cpy ligands. Resonance Raman spectroscopy supports these findings showing enhanced deeb vibrational modes with 532 nm excitation, both deeb and Cpy modes with 473 nm excitation, and exclusively Cpy with 405 nm excitation. Ultrafast spectroscopy reveals a short-lived (∼2 ps) MLCT excited state and a longer-lived (∼20 ps) metal-centered state. Efficient methods to deprotect the ester groups and anchor the complex to mesoporous TiO 2 (anatase) thin films in high surface coverages, Fe(Cpy) 2 (dcb)|TiO 2 σ = 3 × 10 −8 mol/cm 2 , were established. Pulsed light excitation of Fe(Cpy) 2 (dcb)|TiO 2 resulted in rapid excited state injection (k inj > 10 8 s −1 ) and formation of a charge-separated state, Fe III (Cpy) 2 (dcb)|TiO 2 (e), which persists on the millisecond time scale before returning cleanly to the ground state with secondorder kinetics. Injection yields measured 50 ns after light excitation were found to double from Φ = 0.15 with green (532 nm) light to 0.30 with blue (457 nm) light excitation. Incident photon-to-current efficiency (% IPCE) measurements as a function of excitation wavelength in a 0.5 M LiI/I 2 /CH 3 CN electrolyte provide clear evidence for band-selective “hot carrier” injection from the remote Cpy-localized excited state. Collectively, the spectroscopic and photoelectrochemical data indicate that a semiconductor can intercept hot electrons from iron chromophores even when the excited-state dipole is oriented away from the surface-anchoring ligand.

Charge transfer↗
DOE OSTI2020

Highly Reversible Plating/Stripping of Porous Zinc Anodes for Multivalent Zinc Batteries

Zinc continues to garner immense interest due to its versatility as an anode material in several configurations utilizing either alkaline or mild-pH electrolytes. Current research on using mild-pH electrolytes has improved the rechargeability aspect of Zn-based batteries since Zn2+ is solely utilized for plating/stripping of Zn. Several studies have incorporated Zn metal foils, yet, dramatic improvements can be achieved by expressing Zn as a porous structure. Herein, we use a quasi-pulsed electrodeposition process to prepare a conformal Zn coating onto 3D porous copper foam. By tuning the electrodeposition parameters, we achieved an optimal Zn coating that undergoes reversible plating/stripping when tested in symmetric Zn cells, which supported a low overpotential of ~60 mV for up to 100 cycles. We further investigated changes in the surface morphology by studying the Zn surface of both foil and 3D structure using scanning electron microscopy and X-ray micro-computed tomography. Both techniques showed that the Zn foil undergoes dramatic alterations at the surface, which results in inhomogeneous deposition of Zn, whereas the 3D form exhibited minimal changes. Lastly, we paired both Zn foil and 3D Zn with vanadium oxide and demonstrated that the porous structure supports high rate capability and high specific capacity.

25 ENERGY STORAGE↗
DOE OSTI2021

Resilient NdFeB magnet recycling under the impacts of COVID-19 pandemic: Stochastic programming and Benders decomposition

Neodymium-iron-boron (NdFeB) magnets are the most powerful magnets per unit volume sold in the commercial market. Despite the increasing demand for clean energy applications such as electric vehicles and wind turbines, disruptive events including the COVID-19 pandemic have caused significant uncertainties in the supply and demand for NdFeB magnets. Therefore, this study aims to alleviate the risk of supply shortage for NdFeB magnets and the containing critical materials, rare-earth elements (REEs), through the development of a resilient reverse supply chain and logistics network design. We develop scenarios to model the unique impact of the COVID-19 pandemic on the proposed business, incorporating both disruption intensity and recovery rate. We formulate a chance-constrained two-stage stochastic programming model to maximize the profit while guaranteeing the network resiliency against disruption risks. To solve the problem in large-scale instances, we develop an efficient Benders decomposition algorithm that reduces the computational time by 98.5% on average compared to the default CPLEX algorithm. When applied to the United States, the model suggests the optimal facility locations, processing capacities, inventory levels, and material flows for NdFeB magnet recyclers that could meet 99.7% of the demand. To the best of our knowledge, this study is the first to incorporate the impacts of the COVID-19 pandemic to design a resilient NdFeB magnet recycling supply chain and logistics network, leveraging risk-averse stochastic programming.

42 ENGINEERING↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Two-Stage Approach for PV Inverter Engagement in Power Factor Correction and Voltage Regulation

The rapid integration of distributed energy resources, like solar photovoltaics (PVs), can lead to overvolt-age challenges due to reverse power flow and a noticeable decrease in power factor at the substation interface. While existing literature extensively explores utilizing smart inverter capabilities for reactive power flexibility using a volt-var curve (VVC), obtaining time-varying operating points of such curves in real-time is challenging due to computational demands and communication requirements. Similarly, employing optimization-based approaches for reactive power control and active voltage regulation in large-scale distribution feeders is difficult due to the complexity of the problem and the challenges in effectively engaging customer-owned resources. This paper proposes a two-stage strategy to harness smart inverters for reactive power support. The first stage formulates short-term planning by optimally designing VVCs (on a daily or hourly basis) for large-scale solar PVs based on projected system needs and communicating optimal curves to smart inverters in advance. Subsequently, the second stage employs a transactive-based method to involve customer-owned PVs for reactive power support, effectively enhancing overall system performance and addressing real-time demands. In conclusion, the efficacy of this approach will be demonstrated using real-world distribution circuits provided by Vermont Electric Power Company (VELCO) and Vermont Electric Cooperative (VEC).

Poudel, Shiva [Pacific Northwest National Laborato↗
DOE OSTI2022

Self-assembly of nanocrystals into strongly electronically coupled all-inorganic supercrystals

Colloidal nanocrystals of metals, semiconductors, and other functional materials can self-assemble into long-range ordered crystalline and quasicrystalline phases, but insulating organic surface ligands prevent the development of collective electronic states in ordered nanocrystal assemblies. We reversibly self-assembled colloidal nanocrystals of gold, platinum, nickel, lead sulfide, and lead selenide with conductive inorganic ligands into supercrystals exhibiting optical and electronic properties consistent with strong electronic coupling between the constituent nanocrystals. The phase behavior of charge-stabilized nanocrystals can be rationalized and navigated with phase diagrams computed for particles interacting through short-range attractive potentials. By finely tuning interparticle interactions, the assembly was directed either through one-step nucleation or nonclassical two-step nucleation pathways. In the latter case, the nucleation was preceded by the formation of two metastable colloidal fluids.

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

A Combined Computational and Mathematical Analysis of Interconnect Fatigue Potential in Photovoltaic Modules

A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.

14 SOLAR ENERGY↗
DOE OSTI2023

MHD simulation on magnetic compression of field reversed configurations with NIMROD

Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

E × B flow driven electron temperature bifurcation in a closed slot divertor with ion B × ∇ B away from the X-point in the DIII-D tokamak

An electron temperature bifurcation is observed in the small angle slot divertor, which has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Experiments in the DIII-D tokamak and associated SOLPS-ITER modeling with full drifts find a strong interplay between drifts and divertor geometry on divertor dissipation. The coupling of divertor geometry and drift flows can strongly affect the path towards divertor detachment onset as the plasma density is raised. With the strike point on the inner slanted surface and ion B × ∇B away from the magnetic X-point, bifurcative transitions were observed with sharp decrease of T e towards detachment onset both experimentally and computationally. This differs from the situation for the open divertor where the T e cliff was only observed for ion B × ∇B towards the X-point. SOLPS-ITER modeling with full drifts demonstrates that the magnitude of the E × B drift flow is comparable with the main plasma flow. The reversal of both the poloidal and radial E × B flows near the strike point leads to rapid density accumulation right near the separatrix, which results in bifurcative step transition of divertor conditions with cold plasma across the entire divertor target plate. Furthermore, these results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mapping Composition Evolution through Synthesis, Purification, and Depolymerization of Random Heteropolymers

Random heteropolymers (RHPs) consisting of three or more comonomers have been routinely used to synthesize functional materials. While increasing the monomer variety diversifies the side-chain chemistry, this substantially expands the sequence space and leads to ensemble-level sequence heterogeneity. Most studies have relied on monomer composition and simulated sequences to design RHPs, but the questions remain unanswered regarding heterogeneities within each RHP ensemble and how closely these simulated sequences reflect the experimental outcomes. Here, we quantitatively mapped out the evolution of monomer compositions in four-monomer-based RHPs throughout a design-synthesis-purification-depolymerization process. By adopting a Jaacks method, we first determined 12 reactivity ratios directly from quaternary methacrylate RAFT copolymerization experiments to account for the influences of competitive monomer addition and the reversible activation/deactivation equilibria. The reliability of in silico analysis was affirmed by a quantitative agreement (<4% difference) between the simulated RHP compositions and the experimental results. Furthermore, we mapped out the conformation distribution within each ensemble in different solvents as a function of monomer chemistry, composition, and segmental characteristics via high-throughput computation based on self-consistent field theory (SCFT). These comprehensive studies confirmed monomer composition as a viable design parameter to engineer RHP-based functional materials as long as the reactivity ratios are accurately determined and the livingness of RHP synthesis is ensured.

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

Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

We formulate the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that both the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we show that the adjoint wave propagator preserves the dispersion and compensates the amplitude, while the time-reversed adjoint wave propagator behaves identically as the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit Q parameterization, which avoids the implicit Q in the conventional viscoacoustic/viscoelastic full waveform inversion (Q-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. Here, we find that the lossless propagation kernel dominates the velocity kernel, while the dissipation kernel dominates the attenuation kernel over the dispersion kernel.

58 GEOSCIENCES↗
DOE OSTI2026

Assessing the Accuracy of Property Model Predictions for Cost Optimization of Desalination Technologies

Accurate modeling of seawater thermophysical and thermodynamic properties is critical for optimizing desalination processes. This study compares three seawater property models, a Reaktoro multicomponent model, the thermophysical seawater properties library from the Massachusetts Institute of Technology, and a simplified sodium chloride model, in the context of levelized cost of water (LCOW) minimization for reverse osmosis (RO) and mechanical vapor compression systems. Process simulations and cost optimizations reveal that although all three models yield comparable LCOW and specific energy consumption (SEC) estimates under baseline conditions, deviations among their predictions increase with salinity. Relative differences in LCOW and SEC reach up to 6% and 8%, respectively. RO results show greater variability due to differences in osmotic pressure predictions, which affect pressure constraints at high recoveries. Computational performance varies substantially; specifically, Reaktoro simulations are up to 28 times slower than empirical models due to their detailed equilibrium calculations. These results suggest that empirical models offer acceptable accuracy for routine desalination process design, while Reaktoro provides advantages in scenarios requiring detailed speciation, such as scaling or pH adjustment studies. These findings underscore the importance of selecting appropriate property models based on the modeling objective of desalination applications and motivate future work integrating thermodynamic rigor with empirical efficiency.

Physical and chemical properties↗