Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential evolution”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗
DOE OSTI2024

Effect of Zn Addition on Phase Evolution in AlCrFeCoNiZn High–Entropy Alloy

The addition of Zn to AlCrFeCoNi high-entropy alloy (HEA) poses intriguing questions as to how it would affect phase evolution. Herein, the phase evolution in AlCrFeCoNiZn is studied using a combination of experimental techniques (X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, and differential scanning calorimetry) and computational (density-functional theory [DFT], calculation of phase diagrams, and machine-learning) methods. Mechanically alloyed and spark-plasma-sintered AlCrFeCoNiZn assumes a metastable single-phase, body-centered-cubic (BCC) structure that undergoes diffusion-controlled phase separation upon subsequent heat treatment to form separate (Al, Cr)-rich, (Fe, Co)-rich, and (Zn, Ni)-rich phases. The formation of (Al, Cr)-rich phase, not reported previously in AlCrFeCoNi-based HEAs, is attributed to strong clustering tendency of Cr–Zn and Cr–Ni pairs, combined with the strong ordering of Zn–Ni pair, driving out Cr that in turn combines with Al to form a (Al, Cr)-rich phase. In the DFT results, the formation of thermodynamically stable L1 2 phase is shown wherein Cr–Fe–Zn [Al–Ni-Co] preferably occupy1a (000) [3c (0 ½ ½)] positions. Furthermore, the sluggish diffusional transformation to L1 2 phase from BCC precursors is attributed to the small stacking-fault energy of AlCrFeCoNiZn. The equilibrated HEA exhibits a high microhardness of 8.24 GPa with an elastic modulus of 184 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonintrusive projection-based reduced order modeling using stable learned differential operators

Nonintrusive projection-based reduced order models (ROMs) are essential for dynamics prediction in multi-query applications where underlying governing equations are known but the access to the source of the underlying full order model (FOM) is unavailable; that is, FOM is a glass-box. This article proposes a learn-then-project approach for nonintrusive model reduction. In the first step of this approach, high-dimensional stable sparse learned differential operators (S-LDOs) are determined using the generated data. In the second step, the ordinary differential equations, comprising these S-LDOs, are used with suitable dimensionality reduction and low-dimensional subspace projection methods to provide equations for the evolution of reduced states. This approach allows easy integration into the existing intrusive ROM framework to enable nonintrusive model reduction while allowing the use of Petrov–Galerkin projections. The applicability of the proposed approach is demonstrated for Galerkin and LSPG projection-based ROMs through four numerical experiments: 1-D scalar advection, 1-D Burgers, 2-D scalar advection and 1-D scalar advection–diffusion–reaction equations. In conclusion, the results indicate that the proposed nonintrusive ROM strategy provides accurate and stable dynamics prediction.

42 ENGINEERING↗
DOE OSTI2025

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Formation and interior evolution of Mars

The formation and earliest evolution of Mars involved core formation followed by magma ocean solidification and primordial crust formation. Here, in this study, we review the formation and early interior evolution of Mars as understood through meteorite studies. Martian meteorites exhibit ages that span the entire history of the Solar System, but most derive from compositionally distinct sources that formed early in the history of Mars. Geochemical and isotopic studies of Martian meteorites have provided critical information on Mars’ building blocks, bulk composition, and its differentiation into a core, mantle and crust. Radiogenic isotope systems provide vital insights into the petrogenesis of Martian meteorite source regions and into the nature and timescales of primordial differentiation on Mars. Application of short-lived isotope chronometers to Martian meteorites indicates that the Martian core formed early within <10 million years (Ma), while magma ocean differentiation and crust formation occurred soon thereafter, within the first tens of Ma of Solar System history. Bulk silicate Mars is variably depleted in moderately volatile elements relative to chondrites, but to a lesser extent than the Earth. The nature and evolution of more volatile elements and volatile species on Mars (e.g., H 2 O, C, S, halogens) remains an active area of investigation, but overall suggest that Mars is more enriched in volatile elements than the Moon but less than the Earth. While geochemical exploration of Mars coupled with new discoveries of Martian meteorites has provided a wealth of information, the compositional and lithologic diversity on Mars significantly exceeds that seen in the suites of Martian meteorites currently studied. This implies that our current perspective on the planet's evolution is very likely biased, highlighting the need for future sample return missions to Mars.

Kruijer, Thomas S.↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
DOE OSTI2026

Thermal Decomposition Kinetics of 4,6‐Diamino‐5,7‐dinitro‐benzo‐furazan

This experimental study investigated the thermal decomposition kinetics of 4,6-diamino-5,7-dinitro-benzo-furazan (referred to as F1 hereafter)—an important decomposition product of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB—a prototypical insensitive high explosive). Simultaneous differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), and mass spectrometry (MS) measurements were employed to determine the decomposition kinetics of F1 and to track the evolution of product gases. The DSC profiles were measured at 10 different heating rates between 0.025°C/min and 10°C/min. The measured exotherms were influenced by F1 melting at heating rates above 0.25°C/min, and corresponding changes in decomposition enthalpy and TGA mass-loss-rate profiles indicated a transition from solid-to-gas decomposition to an increasing contribution from liquid-to-gas decomposition. Analysis of low-heating-rate DSC data between 0.025°C/min and 0.17°C/min with the extended Prout–Tompkins model yielded an activation energy of 305 kJ/mol for solid-to-gas F1 decomposition, higher than previous values inferred from TATB decomposition models involving F1. This study provides the first direct experimental determination of the energy barrier for F1 decomposition. MS measurements showed that the major gaseous products matched species previously reported for TATB decomposition (e.g., CO 2 , HCN, C 2 N 2 , etc.), with water identified as the dominant product. Furthermore, these results provide important experimental constraints for improving chemical kinetics models of TATB decomposition and for predicting the reactivity, stability, and safety of TATB-based high explosives under long-term aging conditions and abnormal thermal environments.

4,6-Diamino-5,7-dinitro-benzo-furazan↗
DOE OSTI2025

Kinetics of Pyrolysis and Thermal Evolution of Negev Desert Lithologies

The Negev desert in Israel is home to large quantities of organic-rich, shallow marine sedimentary lithologies that could potentially accommodate the disposal of spent nuclear fuel. Previous thermal analyses of Negev carbonates have focused on industrially relevant considerations such as natural gas and oil extraction or pyrolysis for recovering hydrocarbon fuels. Here, this study addresses thermal evolution of the Negev organic-rich carbonate, siliceous, and phosphorite rocks and associated chemical, mineralogical, and microstructural changes that may occur under prolonged thermal loading in the vicinity of spent nuclear fuel disposal systems. Our employed methods include high-temperature X-ray diffraction, high-temperature infrared spectroscopy, and thermal analysis integrating thermogravimetry, differential scanning calorimetry, and mass spectrometry. Further, we apply iterative iso-conversional model-free methods to derive kinetic parameters for thermal decomposition of the Negev organic-rich carbonate rocks from 200 to 550 °C. Our results have provided mechanistic insights into the thermal evolution encompassing water desorption, decomposition of organic matter, and decarbonation of carbonate phases.

58 GEOSCIENCES↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗
DOE OSTI2024

Capturing O2 Desorption Through Isoconversional Kinetics for CFD Application

The present work examines desorption mechanisms and kinetics of the calcium-doped strontium perovskite at 25% calcium content (i.e., Sr_0.75 Ca_0.25 FeO_(3-δ)). Laboratory-scale, low inventory, fixed bed redox cycle experiments were developed and conducted in both dry and steam-based environments for three isothermal temperatures, 450, 500 and 550C. The temporal evolution of the conversion extent of the material sample over the redox cycles were obtained through the analysis of the gaseous products of the redox reactions. Using subsequently isoconversional differential methods, it was found that the activation energy exhibits a strong dependence on the extent of conversion during the desorption of this perovskite. In addition, the reconstruction of the reaction model showed that the desorption kinetics are controlled by a three-step mechanism. A successful a priori verification of the isoconversional desorption kinetics against the dry and steam-based environment desorption experimental data are shown.

Konan, Ndri A.↗
DOE OSTI2025

Virulence and Genetic Diversity of Puccinia spp., Causal Agents of Rust on Switchgrass (Panicum virgatum L.) in the USA

Switchgrass (Panicum virgatum L.) is an important cellulosic biofuel grass native to North America. Rust, caused by Puccinia spp. is the most predominant disease of switchgrass and has the potential to impact biomass conversion. In this study, virulence patterns were determined on a set of 38 switchgrass genotypes for 14 single-spore rust isolates from 14 field samples collected in seven states. Single nucleotide polymorphism (SNP) variation was also assessed in 720 sequenced cloned amplicons representing 654 base pairs of the elongation factor 1-α gene from the field samples. Five major haplotypes were identified differing by 11 out of the 39 SNP positions identified. STRUCTURE, Principal Coordinate Analysis, and phylogenetic analyses divided the rust population into two genetic clusters. Virginia and Georgia had the highest and lowest rust genetic diversity, respectively. Only nine accessions showed a differential disease response between the 14 isolates, allowing the identification of eight races, differing by 1–3 virulence factors. Overall, the results suggested clonal reproduction of the pathogen and a North–South differentiation via local adaptation. However, similar haplotypes and races were also recovered from several states, suggesting migration events, and highlighting the need to further investigate the switchgrass rust population structure and evolution in the USA.

Bahri, Bochra A. (ORCID:0000000159055880)↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗
DOE OSTI2025

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A combinatorially complete epistatic fitness landscape in an enzyme active site

Protein engineering often targets binding pockets or active sites which are enriched in epistasis—nonadditive interactions between amino acid substitutions—and where the combined effects of multiple single substitutions are difficult to predict. Few existing sequence-fitness datasets capture epistasis at large scale, especially for enzyme catalysis, limiting the development and assessment of model-guided enzyme engineering approaches. We present here a combinatorially complete, 160,000-variant fitness landscape across four residues in the active site of an enzyme. Assaying the native reaction of a thermostable β-subunit of tryptophan synthase (TrpB) in a nonnative environment yielded a landscape characterized by significant epistasis and many local optima. These effects prevent simulated directed evolution approaches from efficiently reaching the global optimum. There is nonetheless wide variability in the effectiveness of different directed evolution approaches, which together provide experimental benchmarks for computational and machine learning workflows. The most-fit TrpB variants contain a substitution that is nearly absent in natural TrpB sequences—a result that conservation-based predictions would not capture. Thus, although fitness prediction using evolutionary data can enrich in more-active variants, these approaches struggle to identify and differentiate among the most-active variants, even for this near-native function. Overall, this work presents a large-scale testing ground for model-guided enzyme engineering and suggests that efficient navigation of epistatic fitness landscapes can be improved by advances in both machine learning and physical modeling.

biocatalysis↗