SEARCH · Engineering Papers

Results for “Variational”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Comparison of Pre‐Construction and Operational Wake Loss Estimates for Land‐Based Wind Plants

The overall bias between pre‐construction energy yield assessment (EYA) estimates of wind plant energy production and the achieved operational production is improving in the wind industry, but uncertainty remains high for individual wind plants. Wake effects within wind plants are one of the largest sources of energy loss considered in the EYA process, and previous work shows wake loss estimates to be a major source of disagreement among wind energy consultants who perform EYAs. To better understand the accuracy of wake loss predictions, we compare overall operational wake loss estimates based on supervisory control and data acquisition data to pre‐construction estimates provided by six wind energy consultants for five land‐based wind plants in North America. By augmenting existing approaches for quantifying operational wake losses, we estimate wake losses during the period of record for which operational data are available as well as the expected long‐term wake losses, based on historical reanalysis weather data, to which the EYA estimates are compared. To account for power variations at different turbine locations caused by terrain‐induced wind resource heterogeneity, we correct the operational wake loss estimates using predicted freestream wind speed variations from the Wind Systems Engineering Reynolds‐averaged Navier–Stokes (RANS) tool. We identify long‐term corrected operational wake losses between 1.9% and 6.4% for the five plants, with a mean loss of 4%. For the project deemed most acceptable for operational wake loss assessment, which is located in the simplest terrain and isolated from neighboring plants, the mean EYA wake loss estimate is within 0.7 percentage points of the operational value of 6.4%. For most of the remaining plants, results suggest that wake losses are generally overpredicted by 2.6–6.3 percentage points. However, operational wake losses may be underestimated for many of these projects because of spatial wind resource variations not captured by the RANS model, external wake effects that are unaccounted for in the estimation process, and wind plant blockage effects. To better understand factors that contribute to the observed wake losses, we investigate operational wake losses as a function of wind direction and wind speed. As expected, wake losses are generally concentrated near wind directions that are aligned with rows of closely spaced turbines and at below‐rated wind speeds; however, for some projects, the energy produced by the wind plant exceeds the estimated potential energy of the plant without wake interactions for certain wind directions and wind speeds, suggesting inaccurate assumptions in the wake loss estimation method for those plants. Lastly, we compare predicted and operational wake losses for individual wind turbines, finding that even when overall wake losses are predicted accurately, large uncertainty exists at the turbine level.

17 WIND ENERGY↗

DOE OSTI2020

Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey

Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a sampling-based approach that has widely used for quantification and propagation of uncertainties. However, the standard MC method is often time-consuming if the simulation-based model is computationally intensive. This article gives an overview of modern MC methods to address the existing challenges of the standard MC in the context of UQ. Specifically, multilevel Monte Carlo (MLMC) extending the concept of control variates achieves a significant reduction of the computational cost by performing most evaluations with low accuracy and corresponding low cost, and relatively few evaluations at high accuracy and corresponding high cost. Multifidelity Monte Carlo (MFMC) accelerates the convergence of standard Monte Carlo by generalizing the control variates with different models having varying fidelities and varying computational costs. Multimodel Monte Carlo method (MMMC), having a different setting of MLMC and MFMC, aims to address the issue of UQ and propagation when data for characterizing probability distributions are limited. Multimodel inference combined with importance sampling is proposed for quantifying and efficiently propagating the uncertainties resulting from small data sets. All of these three modern MC methods achieve a significant improvement of computational efficiency for probabilistic UQ, particularly uncertainty propagation. Finally, an algorithm summary and the corresponding code implementation are provided for each of the modern MC methods. The extension and application of these methods are discussed in detail.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Alpha-alpha scattering in the Multiverse

We investigate the phase shifts of low-energy α-α scattering under variations of the fundamental parameters of the Standard Model, namely the light quark mass, the electromagnetic fine-structure constant as well as the QCD θ-angle. As a first step, we recalculate α-α scattering in our Universe utilizing various improvements in the adiabatic projection method, which leads to an improved, parameter-free prediction of the S- and D-wave phase shifts for laboratory energies below 10 MeV. We find that positive shifts in the pion mass have a small effect on the S-wave phase shift, whereas lowering the pion mass adds some repulsion in the two-alpha system. The effect on the D-wave phase shift turns out to be more pronounced as signaled by the D-wave resonance parameters. Variations of the fine-structure constant have almost no effect on the low-energy α-α phase shifts. We further show that up-to-and-including next-to-leading order in the chiral expansion, variations of these phase shifts with respect to the QCD θ-angle can be expressed in terms of the θ-dependent pion mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Random Tensor Networks with Non-trivial Links

Abstract Random tensor networks are a powerful toy model for understanding the entanglement structure of holographic quantum gravity. However, unlike holographic quantum gravity, their entanglement spectra are flat. It has therefore been argued that a better model consists of random tensor networks with link states that are not maximally entangled, i.e., have non-trivial spectra. In this work, we initiate a systematic study of the entanglement properties of these networks. We employ tools from free probability, random matrix theory, and one-shot quantum information theory to study random tensor networks with bounded and unbounded variation in link spectra, and in cases where a subsystem has one or multiple minimal cuts. If the link states have bounded spectral variation, the limiting entanglement spectrum of a subsystem with two minimal cuts can be expressed as a free product of the entanglement spectra of each cut, along with a Marchenko–Pastur distribution. For a class of states with unbounded spectral variation, analogous to semiclassical states in quantum gravity, we relate the limiting entanglement spectrum of a subsystem with two minimal cuts to the distribution of the minimal entanglement across the two cuts. In doing so, we draw connections to previous work on split transfer protocols, entanglement negativity in random tensor networks, and Euclidean path integrals in quantum gravity.

Physics↗

DOE OSTI2023

Biomechanical and Cortical Control of Tongue Movements During Chewing and Swallowing

Tongue function is vital for chewing and swallowing and lingual dysfunction is often associated with dysphagia. Better treatment of dysphagia depends on a better understanding of hyolingual morphology, biomechanics, and neural control in humans and animal models. Recent research has revealed significant variation among animal models in morphology of the hyoid chain and suprahyoid muscles which may be associated with variation in swallowing mechanisms. The recent deployment of XROMM (X-ray Reconstruction of Moving Morphology) to quantify 3D hyolingual kinematics has revealed new details on flexion and roll of the tongue during chewing in animal models, movements similar to those used by humans. XROMM-based studies of swallowing in macaques have falsified traditional hypotheses of mechanisms of tongue base retraction during swallowing, and literature review suggests that other animal models may employ a diversity of mechanisms of tongue base retraction. There is variation among animal models in distribution of hyolingual proprioceptors but how that might be related to lingual mechanics is unknown. In macaque monkeys, tongue kinematics—shape and movement—are strongly encoded in neural activity in orofacial primary motor cortex, giving optimism for development of brain–machine interfaces for assisting recovery of lingual function after stroke. However, more research on hyolingual biomechanics and control is needed for technologies interfacing the nervous system with the hyolingual apparatus to become a reality.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Seasonal non-structural carbohydrate dynamics differ between twig bark and xylem tissues

Despite extensive research non-structural carbohydrate (NSC), the distribution of total NSC (TNC) and its primary components (sugar and starch) to the bark and xylem of twigs remains poorly understood. We determined seasonal NSC dynamics in twig bark and xylem in seven temperate tree species exhibiting different xylem anatomies and leaf habits. Seasonal trends in sugar were similar across species with concentrations peaking at least 15 days earlier in the spring for xylem than bark. Furthermore, evergreens exhibited maximum bark starch concentrations in early spring, followed by declines throughout the growing season, whereas deciduous species exhibited early spring declines in bark starch, followed by late growing season increases. Evergreens exhibited limited seasonal variation in xylem starch concentrations, whereas deciduous species exhibited variation in xylem starch concentrations that was similar to variation in bark starch. With a few exceptions, concentrations and seasonal amplitudes for sugar, starch, and TNC were generally higher in bark than xylem. Sugar concentrations were generally higher than starch, especially in the bark, which resulted in variability of TNC concentrations in bark or xylem. NSC concentrations varied significantly between xylem and the entire twig when bark was not explicitly considered. Averaged across species, sugar, starch, and TNC content in bark accounted for 66%, 54%, and 61% of total twig content, respectively, even though bark only accounted for 53% of the total twig mass. We recommend quantifying bark sugar and starch separately from xylem when determining twig NSC concentration and content.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Multiscale investigation of thermomechanical and compositional developments in Ni alloy 718 under laser processing

Laser processing has been widely employed in various applications due to its exceptional spatial resolution. However, the rapid temperature gradients generated in localized areas present significant challenges for experimental characterization using conventional instruments. To characterize Ni alloy 718 during laser processing, we employed in-situ synchrotron X-ray diffraction with a high-speed detector, a method particularly well-suited for probing processes with high temporal and spatial resolution. Through a series of in-situ experiments, we investigated the local variations in the evolution of microstructures and thermomechanical behaviors within a keyhole mode melt pool. The in situ macroscopic thermomechanical behaviors were quantified using an empirical model derived from diffraction patterns, with experimental results showing reasonable agreement with finite element analysis. Various laser parameters were tested to assess their influences on the residual strains in the melt pools. The results revealed that the residual strain in the keyhole mode melt pool is relatively insensitive to variations in the parameters and is smaller than that in the melt pool created under conduction mode laser scanning. Additionally, we analyzed the shapes of individual diffraction spots, providing insights into the plastic behaviors and compositional developments in the resolidified alloy. The analysis confirmed that compositional variations in a dendritic microstructure manifest as asymmetric broadening of the diffraction spots.

Ni alloy 718↗

DOE OSTI2023

Leveraging the digital thread for physics-based prediction of microstructure heterogeneity in additively manufactured parts

A major limitation of additive manufacturing (AM) processes is that local conditions of material deposition frequently lead to unintentional heterogeneities in microstructure and properties within a single component, despite nominally uniform process conditions. Up to now, there has been no way to a priori determine the distribution of these heterogeneities, requiring expensive trial-and-error approaches to fabrication, testing, and characterization. Here, a physics-based framework for creating a digital representation of the laser powder bed fusion (PBF) process is proposed to predict the variation in solidification behavior that leads to heterogeneous microstructures in an as-built part. By leveraging in situ process data stored in the part’s digital thread, the scan path and process parameters were input into a heat transfer model which predicted solidification data at the melt pool scale. A two-step unsupervised clustering algorithm was used to first cluster the local solidification conditions (12.5µm 3 voxels) and then to cluster the regional behavior on the scale of multiple scan passes and print layers (250µm 3 super-voxels). This process was used to identify regions with similar solidification characteristics for multiple locations in a Stainless Steel 316-L component. The corresponding as-built part was sectioned and characterized using electron backscatter diffraction (EBSD). Quantitative analysis of the pole figures confirmed that the predicted regions of heterogeneity in the solidification conditions corresponded with differences in the observed microstructure. In conclusion, this work shows a viable path for estimating the microstructural heterogeneity for additively manufactured parts to either limit microstructural variation throughout a part or to enable functionality-based variation of the microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗

DOE OSTI2025

Impact of U-10Mo HALEU Fuel Element Tolerances on the Massachusetts Institute of Technology Reactor safety and operational performance – Neutronics

The U.S. is coordinating efforts for the conversion of six U.S. high performance research reactors (USHPRR), including one critical assembly from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts Institute of Technology Reactor (MITR), and achieve similar performance, high assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being considered. Following the preliminary design of the proposed MITR LEU fuel elements using the U-10Mo monolithic alloy, the impact of the fabrication specification was assessed. This work focuses on the analysis of select neutronics characteristics of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., U-10Mo composition, fuel plate thickness, etc.). A separate article submitted to this journal addresses the impact on the thermal hydraulic performance. The analyses in these works are performed based on an all-LEU conversion management plan identified in previous work, in which only the proposed elements are utilized for achieving the conversion of MITR. The variations of two main neutronics characteristics are assessed as a function of the variability of the specification parameters resulting from the fabrication process: the MITR LEU core reactivity and the fuel cycle length. The main findings of this work show that the MITR core can meet the operational requirements during the LEU transition plan under the limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core neutronics characteristics on the specification parameters is highly linear within the specification tolerances. The rates of variation are reported in detail for each parameter and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

DOE OSTI2020

Impact of climate change on water availability and its propagation through the Western U.S. power grid

Climate impact assessment for power systems has tended to focus on the effects of water availability on generation without explicitly considering how these impacts propagate through the interconnected grid to affect generation in other locations. We investigate the impact of projected water availability on the operations of the contemporary Western U.S. power grid, focusing specifically on regional interdependencies. Hydrologic simulations derived from three GCMs (CCSM4, INMCM4 and GFDL-CM3) and two radiative scenarios (RCP4.5 and RCP8.5) are used to force a large-scale, distributed water management model (MOSART-WM), which translates water availability into power generation constraints at hydropower plants and water-dependent thermo-electric plants. Power system dynamics are evaluated using the production cost model PLEXOS. We find that the inter-regional connections across the Western U.S. electricity infrastructure are key in managing drastic variations in regional water availability. Net generation in the Northwest and Northern California is substantially sensitive to local water availability. Generation from the Desert Southwest emerges as a critical resource to compensate for variations in water availability in those regions. Climate change impacts on Northwestern water availability therefore drive future changes in other regions’ generation and regional power flows. The large inter-annual hydrological variability of Northern California exacerbates grid stress arising from long-term changes in water availability in the Northwest. While the regional power flow directions seem insensitive to long-term variations in water availability, the analysis highlights the need to consider other compounding regional factors such as changes in Southern California net load or changes in regional fuel prices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Enhancing flexibility for climate change using seasonal energy storage (aquifer thermal energy storage) in distributed energy systems

Long-term energy storage is expected to play a vital role in the deep decarbonization of building energy sectors, while enhancing the flexibility of buildings to withstand future climate variations. However, it is challenging to design distributed multi-energy systems (DMES) while taking into account the uncertainties introduced by climate change, since stochastic optimization of such systems is difficult. The present study introduces a stochastic optimization model to address this bottleneck, taking into account DMES including aquifer thermal energy storage (ATES) as the long-term thermal storage. For the first time, a novel optimization algorithm links ATES with the DMES optimization model with the support of a simplified geotechnical model. Subsequently, a case study was conducted, focusing on a residential district in Chicago where the impact of future climate condition, energy demand, and solar and wind energy potentials were evaluated using Weather Research and Forecasting (WRF) data (up to 2080) and the EnergyPlus model. The study revealed that ATES is an attractive way to improve the renewable energy penetration level and minimize the dependence on fossil fuels with reasonable support from the grid to assist the fluctuations in both demand and generation. Furthermore, ATES notably reduces fuel consumption and dependence while greatly enhancing the flexibility of the energy system to withstand fluctuations in demand and renewable energy generation brought by future climate variations. These qualities will make ATES an important part of distributed energy systems, even though it is not currently the lowest cost alternative due to lack of technology maturity. Furthermore, the design platform introduced in the present study can be used to design DMES enhancing flexibility to accommodate future climate variations.

25 ENERGY STORAGE↗

DOE OSTI2021

Behavior of the Mo, Tl, and U isotope systems during differentiation in the Kilauea Iki lava lake

Stable molybdenum (Mo), thallium (Tl), and uranium (U) isotope ratios were determined in a suite of samples from the 1959 Kilauea eruption and from Kilauea Iki lava lake with the aim of understanding the effects of igneous differentiation on these isotope systems. The samples range from olivine cumulate with MgO up to 27% to internal differentiates with MgO less than 3%, representing a tholeiitic differentiation series. Molybdenum, Tl, and U behave incompatibly during differentiation, and Mo and U isotope ratios do not systematically vary amongst the different samples. δ 98 Mo values range from -0.17 to -0.31‰ and δ 238 U values range from -0.20 to -0.38‰. Most individual analyses for both isotope systems overlap within measurement uncertainty (± ~0.7 and ~ 0.6, respectively). Mean δ 98 Mo and δ238U values are -0.22 ± 0.08‰ (2σ) and - 0.29 ± 0.09‰ (2σ), respectively, which overlap with Pacific mid ocean ridge basalt (MORB). In contrast, Tl isotopes show small but resolvable variations, with ε 205 Tl ranging from +1.20 to -1.38. The most negative ε 205 Tl values are confined to some of the lowest [Tl] samples, but the ε 205 Tl values do not otherwise vary smoothly with MgO or [Tl]. Possible mechanisms for thallium isotope fractionation are considered (e.g., degassing, water leaching, sulfide fractionation) but none are found to be satisfactory. Overall, the lack of resolvable variation in the Mo and U isotope systems and the small magnitude of heterogeneity in the Tl isotope system indicate that differentiation in tholeiitic systems is unlikely to be a major contributor to global variation in these isotope systems.

58 GEOSCIENCES↗

DOE OSTI2025

A tri-level distribution locational marginal price-based demand response framework

Here, in this paper, we propose a tri-level, nested, two-stage price-based demand response (PBDR) framework that considers distribution locational marginal price (DLMP) as DR enabler between load-serving entities (LSE), demand response providers (DRPs), and customers in the day-ahead distribution market. It enables LSE and customer interactions by using multiple DRPs, positioned in-between, and independently optimizes their objectives. The problem is formulated using linear power flow with approximated power losses and its application in DLMP as DR pricing. The tri-level problem is solved using a nested reformulation & decomposition (R&D) method and tested on the real Indian-108 bus distribution system under various dynamic pricings. Further, the temporal–spatial variations in DLMPs are assessed using fairness criteria. Numerical analyses demonstrate that DLMP applications can effectively improve economic efficiency, and transparency in DR programs valuation with a favorable fairness margin. The results show that DLMP as DR pricing signal induces (0-2) % variation in DLMP for DR participation up to 10 %. Further, it gives over 90 % fairness over temporal–spatial variation for all the customers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Off-design operation and performance of pumped thermal energy storage

In this article, we describe off-design models and control strategies for a Pumped Thermal Energy Storage (PTES) system that uses liquid thermal energy storage: specifically molten salt for hot storage and methanol for cold storage. Off-design conditions arise when load-following, or due to variations in storage tank temperatures or ambient temperatures. We propose a control strategy that uses inventory control to manage the mass flow rate in the thermodynamic cycles, which facilitates load following. We also propose a control strategy for the storage fluid mass flow rates, which are varied to ensure the molten salt is maintained at its design temperature. This maximizes efficiency and minimizes problems with salt freezing or degradation. The cold storage fluid mass flow rate is varied so that the cold tanks have the same state-of-charge as the hot tanks. This leads to variations in cold fluid temperature, but these variations are shown to be acceptably small (e.g. 7.5% increase), and this control method is shown to be simpler and more efficient than an alternative strategy where tanks become unbalanced. The ambient temperature and storage tank temperatures are moved ±50 °C from the design values and the impact on power, duration, and tank temperatures is quantified. Results demonstrate that the proposed control strategy is stable and self-correcting - that is, storage temperatures converge on stable values after two-to-three charge-discharge cycles. When inputs return to design values, the system returns to its design point after two charge-discharge cycles. We also demonstrate that inventory control enables delivery of the target power output even when off-design conditions exist that would normally reduce the power output.

25 ENERGY STORAGE↗

DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Chromium isotopes in an acidic fluvial system: Implications for modern and ancient Cr isotope records

Chromium is an important redox-active element, with significant Cr isotope fractionations driven by redox transformations in the natural environment and quantified through laboratory experiments. Recent work has demonstrated that non-redox-driven transformations can also produce small variations in Cr isotopes, but the role of these processes in controlling environmental Cr isotope variations remains unclear. The highly acidic but biologically diverse acid rock drainage system in Río Tinto, Spain (and affiliated rivers) provides a unique opportunity to examine how non-redox-related processes can substantially affect the isotopic composition of Cr in a natural system. Our results suggest that under conditions where Cr redox transformations are obviated, Cr cycling is largely controlled by formation of Fe-oxyhydroxides and Fe-hydroxysulfates, Cr adsorption, and Cr speciation, with little apparent Cr isotope fractionation between source materials, waters, and sediments. Moderate Cr isotope shifts, however, may be the result of non-redox processes, including adsorption to organic ligands and potentially biological uptake. These results add to a growing body of evidence for the importance of non-redox transformations in contributing to environmental Cr isotope variations, with implications for the use of chromium isotopes in studies of sedimentary Cr isotope records. In particular, this study is a validation of the potential for sulfide weathering in the Proterozoic to deliver relatively unfractionated Cr isotopes to the ocean, despite rising and sustained high oxygen levels.

54 ENVIRONMENTAL SCIENCES↗