Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed predictive control”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Kinetic Modeling of Secondary Organic Aerosol in a Weather-Chemistry Model: Parameterizations, Processes, and Predictions for GOAmazon

Secondary organic aerosol (SOA) forms and evolves in the atmosphere through many pathways and processes, over diverse spatial and time scales. Hence, there is a need to represent these widely-varying kinetic processes in large-scale atmospheric models to allow for accurate predictions of the abundance, properties, and impacts of SOA. In this work, we integrated a kinetic, process-level model (simpleSOM-MOSAIC) into a weather-chemistry model (WRF-Chem) to simulate the oxidation chemistry and microphysics of atmospheric SOA. simpleSOM-MOSAIC simulates multigenerational gas-phase chemistry, autoxidation reactions, heterogeneous oxidation, oligomerization, and phase-state-influenced gas/particle partitioning of SOA. As a case study, the integrated WRF-Chem-simpleSOM-MOSAIC (WC-SSM) model was used to simulate the photochemical evolution downwind of a large city (Manaus, Brazil) in the Amazon and, in turn, study the anthropogenic and biogenic interactions in an otherwise pristine environment. Consistent with previous work, we found that OA was enhanced by up to a factor of four in the urban plume due to elevated hydroxyl radical (OH) concentrations, relative to the background, and that this OA was dominated by SOA from biogenic precursors (80%). Further, in addition to accurately simulating the OA enhancement in the urban plume, the model reproduced the magnitude of the OA oxygen-to-carbon (O:C) ratio and broadly tracked the evolution of the aerosol size distribution. Our work highlights the importance of including an integrated, kinetic representation of SOA processes in an atmospheric model

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)

Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2022

Distributed Real-time Plume Monitoring for Deep Sea Mineral Extraction​

In the emerging industry of deep-sea mining for minerals and deposits (e.g. polymetallic nodules for nickel, cobalt, copper, and manganese), more data is required to understand the effects of sediment plume generation and predict the distribution of disturbed sediment. There are two main sources of plume generation, the first being at the active mining site where the “collector” directly removes the top layer of the sea floor. The other is the “midwater plume” consisting of unwanted sediment that was collected during extraction that is pumped back into the aphotic zone. The vast majority of plume generation is caused by the collector, causing detrimental and long-lasting impacts on seafloor ecosystems due to the lack of wave activity or strong currents at the sea floor. Therefore, it is crucial to invest in the infrastructure to support the study and constant monitoring over a large area of the sea floor where plume generation is present. Due to the limited number of usable channels and power requirements, current subsea wireless communications technologies are not well suited to instrumenting the large areas of the sea floor needed to monitor plume migration. The scope of this effort is to transition experimental demonstrations of high-bandwidth, full-duplex scalable underwater laser communications to the seafloor in an open ocean environment. Specifically tackling challenges associated with the dynamic nature of the subsea world, including but not limited to, deployment logistics, sustainability, and range. The goal is to enable the internet of underwater things for deep sea industries by broadening the capabilities of subsea communications. By using high-precision laser transmitters, many of the challenges current subsea optical systems face can be circumvented, such as power consumption, interference, and bandwidth limitations. This approach lends itself to wireless interlinking multi-node networks, in series or parallel, facilitating the implementation of a wide array of sensor types. This interlinking allows all the data gathered from the network to be processed through a single hardline uplink to the surface, lowering the complexity required for near real-time data processing. Additionally, the laser control systems produce metadata that can be used to help characterize the water column between the nodes. Combining data from various sensors such as turbidity, temperature, current velocity with metadata such as beam attenuation and deflection can produce a high-resolution model of sea floor conditions around an active mining zone. The resulting near real-time model can be used to optimize location and flow rate of the mining operation to minimize and quantify the environmental impact.

Mons, Ishan↗
DOE OSTI2023

Peak Reduction Using Mode Adjustment of Heat Pump Water Heaters in a Residential Neighborhood

Building electrification is putting pressure on distribution grid worldwide. Peak reduction is an important concern that can help reduce the growing stress and allow to defer investments in new capacity. Water heaters represent a convenient way of reducing peak because they are less dependent on weather, and their storage volume allows for asynchronous water heating and hot water use. Previous empirical studies investigated the ability of water heaters to reduce peak through the adjustment of the temperature setpoint. However, not all equipment vendors offer this option. This study aims at understanding the feasibility of peak reduction with an alternative configuration available in the market - by adjusting the device mode rather than temperature setpoint. The peak reduction methodology is tested in an occupied 46-townhome neighborhood located in Atlanta, GA. We find that peak shifting is possible with the adjustable mode approach, with the change in the peak load by 30-60%.

demand response↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electron interference in atomic ionization by two crossing polarized ultrashort pulses

Formation of geometrically regular interference patterns in the photoelectron momentum distributions (PMDs) corresponding to the photoionization of atoms by two single-color, crossing ultrashort pulses is investigated both analytically and numerically. It is shown that, in contrast to the photoionization by monochromatic pulses, PMDs for the ionization by crossing and co-propagating broadband pulses are essentially different (unless both pulses are linearly polarized), namely, when one pulse is linearly polarized along the propagation direction, ˆk, of the circularly polarized (CP) pulse, then interference maxima (minima) of the ionization probability have the form of three-dimensional single-arm regular spirals which are wound along ˆk. Next, the interference maxima (minima) of the ionization probability by a pair of crossing elliptically polarized pulses have the form of either Newton's rings or two-arm Fermat's spirals, depending on the position of a detection plane. Remarkably, these regular patterns occur only for certain values of the pulse ellipticities, and they become distorted for CP pulses. For both above-mentioned pulse configurations, the features of interference patterns depend on the time delay between pulses, their relative electric field amplitude, and relative carrier-envelope phase. Our predictions, illustrated by the numerical results for the ionization of H and He atoms by two orthogonal pulses, are quite general and we expect them to be valid for the ionization of any randomly oriented atomic or molecular target.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Theory and Simulation of Ultrafast Multidimensional Nonlinear X-ray Spectroscopy of Molecules (Final Report)

Emerging X-ray free electron laser (XFEL) beam sources offer new types of probes of matter with unprecedented spatial and temporal resolutions. These experimental advances must be met by robust theoretical and computational tools that provide predictive modeling capacity of the underlining electronic and structural dynamics. The latter will be essential for the design of sophisticated multi-pulse experiments and for their interpretation. The proposed research effort will focus on developing cutting- edge simulation tools for nonlinear multidimensional X-ray/optical spectroscopies and aims to address key questions in Priority Research Opportunities 1 (Probing and controlling electron motion within a molecule) and 3 (Capturing rare events and intermediate states in the transformation of matter) as mentioned in the BES Roundtable Report “Opportunities for Basic Research at the Frontiers of XFEL Ultrafast Science”. XFEL multidimensional nonlinear techniques, which combine sequences of X-ray and possibly optical pulses, provide a unique experimental toolbox for probing the dynamics of core and valence electronic excitations, as well as material structure. Predictive modeling of these dynamical processes requires the combination of analytical theory for nonlinear interactions of light and matter, robust quantum-chemical methodologies for the accurate description of electronic structure of various materials, and multiscale ab initio electron and nuclear dynamics techniques operating beyond Born-Oppenheimer approximation. These challenges will be addressed with three research thrusts (i) Develop and implement theoretical apparatus for modeling a broad range of multidimensional spectroscopic techniques enabled by present and upcoming XFEL facilities. This thrust also includes the incorporation of a computational module in the DOE supported open-source NWChem computational chemistry package as well as the development of other open-source codes ready for dissemination across a broad user base; (ii) Propose and design new multi-pulse experiments that make use of the capabilities of the incoming LCLS-II facility; (iii) Perform selected applications to specific molecular systems that can be carried out at LCLS-II and demonstrate how these X-ray sources may be used to study nonadiabatic dynamics through conical interactions, electronic correlations in multi-core excitons, and charge transfer/energy transfer processes. The proposed research will be carried out by a multi-disciplinary four-institution research team which combines academia and national laboratories and spans the broad and necessary expertise in theoretical spectroscopy, nonlinear optics, quantum chemistry, molecular non-adiabatic dynamics and code development. The work will be performed in a highly interactive team environment with junior researchers shared between institutions thus cementing cross- disciplinary interactions. The developed simulation tools will be immediately deployed for XFEL facility users, both experimentalists and theorists, via freely distributed codes and databases. Altogether, this project will facilitate establishing XFEL-based multidimensional spectroscopies as a novel diagnostic tool for monitoring electronic and structural dynamics in molecular materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗
DOE OSTI2024

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗
DOE OSTI2023

Predictive performance of multi-model ensemble forecasts of COVID-19 across European nations

Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022. Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1–4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models’ predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models’ forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models’ past predictive performance. Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models’ forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models’ forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models’ forecasts of deaths (N=763 predictions from 20 models). Across a 1–4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models. Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks. Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 ( https://www.nfdi4health.de/task-force-covid-19-2 ) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Microsimulation of Energy and Flow Effects from Optimal Automated Driving in Mixed Traffic

In this paper we study the energy and traffic impact of a proposed Anticipative Cruise Controller in a PTV VISSIM microsimulation environment. We dissect our controller into two parts: 1. the unconnected mode, active when following a human-driven vehicle, and 2. the connected mode, active when following another automated vehicle equipped with connectivity. Probabilistic constraints balance safety considerations with inter-vehicle compactness, and vehicle constraints for acceleration capabilities are expressed through the use of powertrain maps. Emergent highway traffic scenarios are then modeled using time headway distributions from empirical traffic data. To study the impact of automation over a range of demands of free-flow to stop-and-go, we vary vehicle flux from low to high and vary automated vehicle penetration from low to high. When examining all-human driving scenarios, network capacity failed to meet demand in high-volume scenarios, such as rush-hour traffic. We further find that with connected automated vehicles introduced, network capacity was improved to support the high-volume scenarios. Finally, we examine energy efficiencies of the fleet for conventional, electric, and hybrid vehicles. We find that automated vehicles perform at a 10%-20% higher energy efficiency over human drivers when considering conventional powertrains, and find that automated vehicles perform at a 3%-9% higher energy efficiency over human drivers when considering electric and hybrid powertrains. Due to secondary effects of smoothing traffic flow and reducing unnecessary braking, energy benefits also apply to human-driven vehicles that interact with automated ones. Such simulated humans were found to drive up to 10% more energy-efficiently than they did in the baseline all-human scenario.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Synchronous Machine Governor Upgrade

Conventional generation sources play a critical role in the stability and reliability of the electrical grid, particularly as we transition towards more renewable energy sources. To understand and accurately emulate their behavior for optimizing grid operations and ensuring seamless integration with renewable technologies, it is essential to better emulate the grid- and plant-level impacts of conventional generation sources, such as natural gas (NG) driven heat recovery steam generators (HRSGs) and combustion turbines (CTs). Therefore, a governor model is developed in a programmable logic controller (PLC) to investigate the performance of the conventional generator under various dynamic operating conditions and to identify the impact on grid stability in a controlled environment. The governor model aims to enable the hardware-in-the-loop (HIL) based emulation of these conventional generation sources using the existing 2 MVA synchronous machine/generator that is driven by a flexible 2.5 MW variable speed drive. This setup will allow us to replicate the dynamic characteristics and response behaviors of NG-driven HRSGs and CTs. The controls for the emulated conventional plants follow the industry standard and are adjustable, ensuring they accurately reflect the operational capabilities and limitations of real-world systems. These controls include load-following capabilities, ramp rates, startup and shutdown sequences, and emissions characteristics. By incorporating these adjustable controls, we aim to capture the nuanced impacts of conventional generation, such as their ability to provide ancillary services like frequency regulation, voltage support, and spinning reserve. In this report, we simulate two types of dynamic operations: grid-connected and islanding. For each dynamic operation, representative starting sequences are tested, including turbine purge, ignition, speed ramping up, generator excitation and synchronizing, and breaker close. The HIL based tests provides insights for field deployment, specifically the high-fidelity governor model provides results to predict the potential stability and reliability risk and suggest possible integration measures (e.g., generation and load balancing, tuning of governor control parameters). Ultimately, this enhanced emulation capability will be integrated into our Advanced Research on Integrated Energy Systems (ARIES), enabling us to conduct comprehensive studies on the interactions between conventional and renewable energy sources. By better understanding these interactions, we can develop strategies to optimize the overall performance and reliability of the grid. This will support the deployment of advanced grid management techniques, such as demand response, grid-forming inverters, and energy storage systems. The main contributions are summarized as follows: (1) This report introduces a PLC-based governor model for gas turbines. This model accurately simulates the dynamic behavior of conventional generation sources under various operational scenarios; (2) The model is integrated with an HIL testbed that includes a 2.5 MW variable speed drive and a 2 MVA synchronous machine. This setup enables realistic, real-time emulation of conventional power plants, particularly NG driven HRSGs and CTs; (3) The developed model is adaptable to various gas turbine configurations and allows for precise control over parameters such as MW ramp rates. This flexibility makes it a valuable tool for future research and industry collaboration; and (4) By incorporating the model into the National Renewable Energy Laboratory's Advanced Research on Integrated Energy Systems, the report lays the groundwork for future studies on interactions between conventional and renewable energy sources, enhancing the ability to develop advanced grid management strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗
DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗