Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “utility scale”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Data-Driven Discovery of Bimetallic Nanoparticles Catalysts for the Hydrogenolysis of Polyethylene

Supported platinum nanoparticles are known to convert polyolefins to high-quality liquid hydrocarbons with hydrogen under relatively mild conditions. However, no systematic study has been undertaken using bimetallic catalysts for polyethylene upcycling. Specifically, a total of 98 monometallic and bimetallic combinations (Ag, Cr, Co, Cu, Fe, Ga, In, Mn, Ni, Pd, Pt, Rh, Ru, Zr) on alumina were synthesized utilizing surface organometallic chemistry (SOMC) technique via robotic platform. These were investigated at a small scale (10 mg of catalyst and 50 mg of polyethylene) for their activity for the hydrogenolysis of polyethylene in a high-throughput batch reactor. Combinations of Ni and Co were selected as candidates with high activity toward conversion into paraffin oils. Reaction conditions were optimized with Ni/Co/Al 2 O 3 catalyst at a larger scale (300 mg catalyst and 3 g polyethylene) to obtain a high yield (93.1%) of paraffin wax with desired properties (M n = 380 Da) and low polydispersity (Đ = 1.2). Ni/Co/Al 2 O 3 was compared against Co/Ni/Al 2 O 3 to understand the role of the deposition sequence. When Co is deposited before Ni, a layer of cobalt aluminate is formed upon reduction, stabilizing the deposition of 5 nm metallic Ni particles. When nickel is deposited before Co, particles are larger (average >20 nm) and more oxidized (Ni δ+ in NiAl 2 O 4 ), decreasing the availability of the catalytically active metallic Ni. In conclusion, the difference in electronic environments was also described by DFT calculations, which revealed that smaller 3D clusters of Ni are preferred on CoAl2O4 over the 3D clusters on NiAl 2 O 4 and that these smaller clusters are more reducible, as confirmed experimentally.

Polymer↗
DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Towards a RAG-based summarization for the Electron Ion Collider

Abstract The complexity and sheer volume of information — encompassing documents, papers, data, and other resources — from large-scale experiments demand significant time and effort to navigate, making the task of accessing and utilizing these varied forms of information daunting, particularly for new collaborators and early-career scientists.To tackle this issue, a Retrieval Augmented Generation (RAG)-based Summarization AI for EIC (RAGS4EIC) is under development. This AI-Agent not only condenses information but also effectively references relevant responses, offering substantial advantages for collaborators. Our project involves a two-step approach: first, querying a comprehensive vector database containing all pertinent experiment information; second, utilizing a Large Language Model (LLM) to generate concise summaries enriched with citations based on user queries and retrieved data. We describe the evaluation methods that use RAG assessments (RAGAs) scoring mechanisms to assess the effectiveness of responses. Furthermore, we describe the concept of prompt template based instruction-tuning which provides flexibility and accuracy in summarization. Importantly, the implementation relies on LangChain [1], which serves as the foundation of our entire workflow. This integration ensures efficiency and scalability, facilitating smooth deployment and accessibility for various user groups within the Electron Ion Collider (EIC) community. This innovative AI-driven framework not only simplifies the understanding of vast datasets but also encourages collaborative participation, thereby empowering researchers. As a demonstration, a web application has been developed to explain each stage of the RAG Agent development in detail. The application can be accessed athttps://rags4eic-ai4eic.streamlit.app.[A tagged version of the source code can be found inhttps://github.com/ai4eic/EIC-RAG-Project/releases/tag/AI4EIC2023_PROCEEDING.]

Instruments & Instrumentation↗
DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗
DOE OSTI2026

DP-TwoLevel: two-stage gradient subspace learning for differentially private federated learning

Federated learning (FL) enables collaborative model training across distributed data sources without sharing raw data, but faces fundamental challenges in communication efficiency and privacy. Differentially private (DP) training mitigates information leakage but introduces noise that degrades model performance, especially in high-dimensional settings. We propose DP-TwoLevel, a hierarchical gradient projection method that improves utility under fixed DP constraints by exploiting low-dimensional structure in model updates. Our approach learns a two-level PCA-based representation of gradients and applies DP noise in a reduced-dimensional subspace, thereby lowering the effective noise magnitude while preserving dominant signal components. We evaluate the method across three datasets (MNIST, Fashion-MNIST, CIFAR-10) and three privacy regimes (ϵ∈0.5, 1.0, 2.0). Across nine experimental settings, DP-TwoLevel consistently outperforms DP-FedAvg, achieving an average accuracy improvement of 9.44%, with larger gains observed in lower ϵ(higher-noise) regimes (up to +22.31%). We further analyze scalability across models ranging from 100K to 1.49M parameters and identify a variance-based success criterion: performance remains strong when the projection preserves more than 75% of gradient variance, degrades in a marginal regime (65–75%), and fails below this threshold. Our results demonstrate that structure-aware dimensionality reduction can significantly improve the privacy–utility tradeoff in FL without modifying formal privacy guarantees. We also provide empirical evidence of scaling limitations for global projections and motivate per-layer extensions for larger models.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗
DOE OSTI2022

Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers

In recent years, the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) has been redesigned to better enable the use of application-specific and third-party algebraic solvers and data structures. Throughout this work, we have adhered to specific guiding principles that minimized the impact to current users while providing maximum flexibility for later evolution of solvers and data structures. The redesign was done through the addition of new linear and nonlinear solvers classes, enhancements to the vector class, and the creation of modern Fortran interfaces. The vast majority of this work has been performed “behind-the-scenes,” with minimal changes to the user interface and no reduction in solver capabilities or performance. These changes allow SUNDIALS users to more easily utilize external solver libraries and create highly customized solvers, enabling greater flexibility on extreme-scale, heterogeneous computational architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data for High-Level β-Carotene Production from Xylose by Engineered Saccharomyces cerevisiae without Overexpression of a Truncated HMG1 (t HMG1 )

β‐Carotene is a natural pigment and health‐promoting metabolite, and has been widely used in the nutraceutical, feed, and cosmetic industries. Here, we engineered a GRAS yeast Saccharomyces cerevisiae to produce β‐carotene from xylose, the second most abundant and inedible sugar component of lignocellulose biomass. Specifically, a β‐carotene biosynthetic pathway containing crtYB, crtI, and crtE from Xanthophyllomyces dendrorhous was introduced into a xylose‐fermenting S. cerevisiae . The resulting strain produced β‐carotene from xylose at a titer threefold higher than from glucose. Interestingly, overexpression of t HMG1 , which has been reported as a critical genetic perturbation to enhance metabolic fluxes in the mevalonate pathway and β‐carotene production in yeast when glucose is used, did not further improve the production of β‐carotene from xylose. Through fermentation profiling, metabolites analysis, and transcriptional studies, we found the advantages of using xylose as a carbon source, instead of glucose, for β‐carotene production to be a more respiratory feature of xylose consumption, a larger cytosolic acetyl‐CoA pool, and an upregulated expression level of rate‐limiting genes in the β‐carotene‐producing pathway, including ACS1 and HMG1 . As a result, 772.8 mg/L of β‐carotene was obtained in a fed‐batch bioreactor culture with xylose feeding. Considering the inevitable large scale production of xylose when cellulosic biomass‐based bioeconomy is implemented, our results suggest xylose utilization is a promising strategy for overproduction of carotenoids and other isoprenoids in engineered S. cerevisiae .

Conversion↗
DOE OSTI2022

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Developing nanometer scale, atomically precise metallo-catalysts with Molecular Lego

We propose to create atomically precise, highly robust, nanometer scale macromolecules that organize biomimetic metal binding groups in designed three-dimensional pockets to accelerate specific polymerization reactions. These atomically precise catalysts will create atomically precise polymers with improved mechanical and environmental properties from renewable building blocks. We seek to create and study nanometer scale catalysts that mimic metallo-enzymes in their activity and selectivity because these larger catalysts will be able to make more extensive non-covalent contacts with the transition state of the growing polymer. These catalysts will have extended lifetimes because they will better protect the metals and higher activity because additional reactive groups can be brought in close contact with the active site to enhance reactivity. They will be far more tolerant to extremes of temperature and will not denature in non-aqueous solvents because they are held together by multiple, strong, covalent bonds. Our approach to catalysts could be used to develop catalysts to replace noble metal-based catalysts when targeting other reactions, because the geometry of the active site enforced by the Molecular Lego scaffolding will access new reactivity, mimicking how nature uses earth abundant metals for much of its catalysis. We will initially target Lewis acid-based catalysts that assemble aliphatic polyesters with improved stereoselectivity, tacticity and complex alternation of monomers. Aliphatic polyesters have received growing attention as attractive, environmentally benign, and sustainable alternatives to polymers developed from petroleum feedstocks. Aliphatic polyesters are environmentally friendly because they undergo facile hydrolytic degradation to benign products, and they have excellent properties and high biocompatibility. We will assemble these catalysts using unique “Spiroligomers” (aka Molecular Lego) chemistry developed in the Schafmeister group combined with the organometallic chemistry expertise of the Dobereiner group. Polymers are an attractive synthetic target of nanoscale atomically precise catalysts because their resulting structure, stereochemistry and tacticity has profound impact on their properties (crystalline vs amorphous) and provides a readout of the catalytic mechanism. Polymer synthesis also serves as an example of “atomically precise manufacturing” where atomically precise nanoscale catalysts construct polymers with atomic precision at increasing length scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unsupervised Process Anomaly Detection and Identification Using the Leave-One-Variable-Out Approach

Automated anomaly detection and identification can signal equipment issues and pinpoint causes in large-scale industrial systems. For systems with limited failure history, unsupervised machine learning methods can be utilized as they do not require past failures. This study introduces the leave-one-variable-out (LOVO) model, which masks one variable at a time to predict the others, learning underlying process correlations. Detection performance was assessed with synthetic and experimental data, while identification performance used only synthetic data due to its ability to generate labeled anomaly types. For detection using synthetic data, the LOVO model generally outperformed comparative models; while using experimental data, the comparative methods outperformed the LOVO model. However, the comparative methods required selecting a latent size, and these conclusions pertain to using the optimal size. In practice, it would not be feasible to always select the optimal value, and incorrect selections impacted performance. In contrast, the LOVO model does not require a latent space. For identification using synthetic data, the LOVO model was slightly outperformed in interpretability and repeatability but still demonstrated impressive results. These outcomes suggest that the LOVO model is an effective model and may be more easily implemented without the challenging tuning process of selecting a latent size.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Impact of Self-shielding Minihalos on the Lyα Forest at High Redshift

Abstract Dense gas in minihalos with masses of 10 6 −10 8 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly α forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Ly α opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10 −12 s −1 , as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Ly α absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Ly α flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc −1 (or 10 −3 km −1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10 −12 s −1 ) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly α forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Agent-based modeling for multimodal transportation of CO 2 for carbon capture, utilization, and storage: CCUS-agent

Here, to understand the system-level interactions between the entities in Carbon Capture, Utilization, and Storage (CCUS), an agent-based foundational modeling tool, CCUS-Agent, is developed for a large-scale study of transportation flows and infrastructure in the United States. Key features of the tool include (i) modular design, (ii) multiple transportation modes, (iii) capabilities for extension, and (iv) testing against various system components and networks of small and large sizes. Five matching algorithms for CO 2 supply agents (e.g., powerplants and industrial facilities) and demand agents (e.g., storage and utilization sites) are explored: Most Profitable First Year (MPFY), Most Profitable All Years (MPAY), Shortest Total Distance First Year (SDFY), Shortest Total Distance All Years (SDAY), and Shortest distance to long-haul transport All Years (ACAY). Before matching, the supply agent, demand agent, and route must be available, and the connection must be profitable. A profitable connection means the supply agent portion of revenue from the 45Q tax credit must cover the supply agent costs and all transportation costs, while the demand agent revenue portion must cover all demand agent costs. A case study employing over 5500 supply and demand agents and multimodal CCUS transportation infrastructure in the contiguous United States is conducted. The results suggest that it is possible to capture over 9 billion tonnes (GT) of CO 2 from 2025 to 2043, which will increase significantly to 22 GT if the capture costs are reduced by 40 %. The MPFY and SDFY algorithms capture more CO 2 earlier in the time horizon, while the MPAY and SDAY algorithms capture more later in the time horizon.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Vibrational metrics for the evaluation of internal conditions in a scaled nuclear dry storage cask

In many countries, permanent repositories for high burnup spent nuclear fuel assemblies have yet to be established. As a result, much of the spent fuel from nuclear reactors is contained within dry storage casks, with many being beyond their designed service life. The assessment of the internal structural integrity of these casks and the fuel assemblies contained therein is of critical importance for both extended storage and transport to permanent repositories. The large size, structural complexity, and inaccessibility of the interior of the casks make this task challenging. To address these difficulties, a 1:6 scaled model based closely on the design of a Transnuclear (TN-32) dry storage cask was fabricated to facilitate controlled studies of these structures in the laboratory. Vibrational spectroscopy was used to evaluate the state of the cargo and internal structures within the cask utilizing only measurements on the outer surface. Using modes identified through Finite Element modeling corresponding to those previously measured on a full-scale TN-32 cask, we report on the development of amplitude- and phase-based metrics that are sensitive to internal conditions in the lab cask. Steel rod bundles and steel shot were used as surrogates for intact and damaged fuel assemblies, and various internal configurations of these materials were investigated. Finally, the metrics were based on acquired spectra involving the (1, 2) global bending mode and the (2, 1) radial-with-shearing mode. The results show that the metrics are sensitive to the condition of a single assembly and have some ability to determine the locations of damaged and empty slots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗
DOE OSTI2026

Fluorescence Signatures of Rare Earth Metals during Precipitation in Various Conditions

Fluorescence spectroscopy is a widely used sensor methodology that analyzes light emitted from a compound or element as it decays from an excited state. This technique is very sensitive and selective, which is ideal to characterize analytes at lower limits of detection. Key example targets of significant industry and research interest include rare earth elements (REEs) such as dysprosium (Dy) and europium (Eu). These are widely used in advanced technologies including semiconductors, electric vehicle motors, lasers, and permanent magnets. Identifying new sources and responsible reutilization of REEs is essential, and new approaches to extract and recycle REEs could be notably enhanced through the integration of on-line sensors. The sensors can support faster process design, informed scale-up, and cost-effective deployment. This study covers the initial exploration of applying fluorescence-based on-line monitoring to REEs within a precipitation process. This study demonstrates the successful scale-up of a fluorescence -based sensing approach, from stationary cuvettes and small-volume microfluidic devices to continuous flow systems operating at the bench scale (10-25mL). This work also provides initial insight into the challenges of signal’s effects and utility within a turbid environment. Using a modular design for monitoring flowing solutions in a flow tube, fluorescence can be characterized for a variety of analytical targets. In this study, detection performance parameters between the cuvette and flow tube system were compared. Additionally, the response of Dy during precipitation by sodium bicarbonate in the two measurement designs was explored. This letter represents a starting point to bridge the gap between traditional fluorescence sensor measurements in a cuvette to future developments that explore the ability to integrate fluorescence sensors into extraction and separation processes at industrially relevant scales.

fluorescence↗
DOE OSTI2020

Voluntary Renewable Energy Procurement Programs in Regulated Utility Markets

Multinational corporations are increasingly purchasing renewable energy (RE) to fulfill international commitments, reduce supply chain emissions, limit environmental impact, and secure stable and affordable electricity. The potential scale of global corporate RE purchasing has been estimated to be nearly 100 GW and growing; however, many markets still lack supportive enabling environments for corporations to access RE through on-site project development or private sector transactions. Even where allowed by law, on-site generation may be insufficient due to space and technical constraints, introducing additional challenges. Utility green pricing and utility green tariff programs backed by renewable energy certificates (RECs), referred to as utility green procurement programs (GPPs) in this report, offer powerful market-based solutions to utilities, regulators, and policymakers to provide corporate and other consumers with RE product options while generating revenue to support RE development. As RE markets expand around the globe, GPPs have proven to be effective mechanisms in market regimes, ranging from fully integrated state-owned utilities to broadly liberalized power markets. GPPs utilize RECs, which are a type of energy attribute certificate and closely resemble guarantees of origin, to track and ultimately monetize RE attributes that corporations and other buyers must procure to demonstrate progress against their RE commitments and make public claims of RE use. Time-tested and transparent REC accounting mechanisms provide market confidence, while at the same time offering flexibility to utilities, regulators, and customers for a range of applications. RECs are used in all types of electricity market structures, as indicated in Figure 1, but REC-based program designs and supporting components differ depending on the type of market. Liberalized markets allow for customers to contract directly with generators for electricity and RECs, while traditionally regulated markets with vertically integrated utility structures may have greater restrictions on generation asset ownership and electricity sales. RE markets can be further defined as being either mandatory or voluntary. Mandatory markets require suppliers to deliver specified amounts of RE to grid customers, such as under a renewable portfolio standard (RPS), while voluntary markets involve no legal mandates, but demand is driven by self-imposed customer goals.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗