Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution operators”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Sustainable Data Evolution Technology (SDET) for Power Grid Optimization (Final Report)

This project develops an innovative sustainable data evolution technology (SDET) to create open-access power grid datasets and facilitate updates to these datasets by the power grid community. The lack of open-access, realistic large-scale datasets significantly limits the ability of researchers, developers and ultimately end users to develop, benchmark and compare new methods and tools for optimizing the operation and planning of the grid, which leads to slow adoption by end users. Available open datasets include the IEEE sponsored 14-300 bus transmission test cases and distribution test feeders, the Edinburgh Power Systems Test Case Archive (39-2224 buses), and the Polish cases (2383-3375 buses) available with MATPOWER. These datasets serve certain purposes, but for power grid optimization, they are insufficient in one or more of these aspects: 1) they are too small; 2) they are incomplete in both models and scenarios for optimization purposes; and 3) they are static and are not keeping up with grid needs. The SDET technology delivers large-scale realistic datasets and data-creation tools capable of generating new datasets. We work closely with the Category-2 data repository project awarded by the same GRID DATA program. The generated datasets and data creation tools are compatible with and available through the data repository. The objective is to make this a sustained effort within and beyond the ARPA-E GRID DATA program so that the datasets can evolve over time and meet the current and future needs for power grid optimization and potentially other applications in power grid operation and planning. The novel and disruptive features of SDET differentiate itself from prior work in the following aspects: 1) SDET derives features and metrics for both transmission and distribution (T&D) systems by analyzing many public and private datasets provided by our industry partners National Rural Electric Cooperative Association (NRECA) and PJM Interconnection (PJM). 2) SDET develops data-creation tools and use these tools to generate large-scale open-access realistic datasets that comply with the metrics for both T&D systems. 3) SDET validates the created datasets using industry tools provided by our vendor partner GE Energy Solution. The data creation tools, as well as the created datasets, are integrated into the data repository (DR. POWER) funded by ARPA-E GRID DATA program in Category 2. Users can use the datasets and access these data creation tools to create and update datasets, thus datasets can evolve per user requirements and power grid development. The SDET is a novel and disruptive technology compared to current ad hoc and static dataset generation. It is expected the SDET technology will fundamentally enhance the development, assessment, and adoption of new methods and tools for power grid optimization and other applications, which in turn would significantly improve the reliability, resiliency, and efficiency of the power grid – a major ARPA-E mission objective.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

The 13th Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-IV Survey Mapping Nearby Galaxies at Apache Point Observatory

Here, the fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in 2014 July. It pursues three core programs: the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2), Mapping Nearby Galaxies at APO (MaNGA), and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). As well as its core program, eBOSS contains two major subprograms: the Time Domain Spectroscopic Survey (TDSS) and the SPectroscopic IDentification of ERosita Sources (SPIDERS). This paper describes the first data release from SDSS-IV, Data Release 13 (DR13). DR13 makes publicly available the first 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA. It includes new observations from eBOSS, completing the Sloan Extended QUasar, Emission-line galaxy, Luminous red galaxy Survey (SEQUELS), which also targeted variability-selected objects and X-ray-selected objects. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification, and new reductions of the SDSS-III APOGEE-1 data, improving stellar parameters for dwarf stars and cooler stars. DR13 provides more robust and precise photometric calibrations. Value-added target catalogs relevant for eBOSS, TDSS, and SPIDERS and an updated red-clump catalog for APOGEE are also available. This paper describes the location and format of the data and provides references to important technical papers. The SDSS web site provides links to the data, tutorials, examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6 yr operations of SDSS-IV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Soil Moisture Drought Monitoring and Forecasting Using Satellite and Climate Model Data over Southwestern China

In this paper, an experimental soil moisture drought monitoring and seasonal forecasting framework based on the Variable Infiltration Capacity model (VIC) over southwestern China (SW) is presented. Satellite precipitation data are used to force VIC for a near-real-time estimate of land surface hydrologic conditions. Initialized with satellite-aided monitoring (MONIT), the climate model (CFSv2)-based forecast (MONIT+CFSv2) and ensemble streamflow prediction (ESP)-based forecast (MONIT+ESP) are both performed. One dry season drought and one wet season drought are employed to test the ability of this framework in terms of real-time tracking and predicting the evolution of soil moisture (SM) drought, respectively. The results show that the skillful CFSv2 climate forecasts (CFs) are only found at the first month. The satellite-aided monitoring is able to provide a reasonable estimate of forecast initial conditions (ICs) in real-time mode. In the presented cases, MONIT+CFSv2 forecast exhibits comparable performance against the observation-based estimates for the first month, whereas the predictive skill largely drops beyond 1 month. Compared to MONIT+ESP, MONIT+CFSv2 ensembles give more skillful SM drought forecast during the dry season, as indicated by a smaller ensemble range, while the added value of MONIT+CFSv2 is marginal during the wet season. Furthermore, a quantitative attribution analysis of SM forecast uncertainty demonstrates that SM forecast skill is mostly controlled by ICs at the first month and that uncertainties in CFs have the largest contribution to SM forecast errors at longer lead times. This study highlights a value of this framework in generating near-real-time ICs and providing a reliable SM drought prediction with 1 month ahead, which may greatly benefit drought diagnosis, assessment, and early warning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multifold Insights for Power System Dynamics From Data Assimilation: Meeting Current Challenges

The world’s electric power grids have evolved over the last 120 years from a single power line to today’s large networks. The evolution will continue at an accelerated rate with extensive smart grid development worldwide. Here for example, the U.S. government has set a goal of reaching 100 percent carbon pollution-free electricity by 2035, while the Department of Energy (DOE) shared a goal to deploy 30 GW and 110 GW of offshore wind by 2030 and 2050, respectively. To meet such ambitious goals, in the years to come, a significant percentage of electricity will come from intermittent renewable sources, electric vehicles (EVs), and be supplied to a vast number of loads that will actively respond to grid conditions and incentive signals. This development is largely driven by environmental and economic factors, such as reducing carbon emissions and saving electricity cost for consumers. These energy resources lead to new uncertain behaviors and dynamics which the grid has never seen and were not considered in its design. Operating such a dynamic grid with sufficient reliability and efficiency is a monumental challenge.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

2020 Annual Report of Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve

The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of available drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns. The process involves an initial assessment of the pillar-to-diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the structural integrity with more detail using geomechanical models. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent "red flags" for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern's drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report includes an update to the baseline drawdowns for each cavern, and provides an initial assessment of the evolution of drawdown expenditure for several caverns

02 PETROLEUM↗
DOE OSTI2019

Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. Here, the accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics

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

Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve (2021 Annual Report)

The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPRs ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of available drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns (Sobolik et al., 2014; Sobolik 2016). The process involves an initial assessment of the pillar - to - diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the structural integrity with more detail using geomechanical models. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent red flags for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales , purchases, or exchanges authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the caverns drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report is the latest in a series of annual reports, and it includes the baseline available drawdowns for each cavern, and the most recent assessment of the evolution of drawdown expenditure for several caverns .

02 PETROLEUM↗
DOE OSTI2022

Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve (2022 Annual Report)

The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of a vailable drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns (Sobolik et al., 2014; Sobolik 2016). The process involves an initial assessment of the pillar-to-diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the s tructural integrity with more detail using geomechanical models. Finite - element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent "red flags" for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales , purchases, or exchanges authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern's drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report is the latest in a series of annual reports, and it includes the baseline available drawdowns for each cavern, and the most recent assessment of the evolution of drawdown expenditure for several caverns.

02 PETROLEUM↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The evolution of the Nevada National Security site waste acceptance criteria - 15312

The Nevada National Security Site Waste Acceptance Criteria (DOE/NV--325) is an evolution of requirements spanning almost three decades of waste acceptance at the Nevada National Security Site (NNSS). The original waste acceptance criteria (WAC) document was written in 1978 (document NVO-185) and underwent multiple revisions as time passed and disposal experience was gained. This document came into being when the NNSS first started receiving significant quantities of low-level radioactive waste (LLW) from offsite U.S. Department of Energy (DOE) generators, such as the Mound and Rocky Flats Plants. The original WAC was far less prescriptive and rigorous than the current criteria since many of the current DOE order directives, Federal/State laws, and regulations did not exist. The current WAC reflects the results of numerous revisions in order to address the nuances of the varied waste streams and package configurations accepted for disposal at the NNSS. This includes waste streams that exhibit challenges not previously faced or covered by the WAC and therefore, from time to time, require new information and or changes to be made to the WAC. These wastes now represent a growing percentage of the remaining waste at cleanup sites across the DOE Complex. These waste streams present challenges, such as the technical basis for ensuring safe transportation to and disposal at the NNSS. The current WAC (Rev. 10 dated June 2013), is Available online at: http://www.nv.energy.gov/library/publications/Environmental/DOENV{sub 3}25.pdf). Other challenges for the NNSS include incorporating new and emerging technologies and addressing stakeholder concerns. Addressing these challenges sometimes results in changes to the NNSS WAC. Drivers for revising the NNSS WAC over the years have included operating experience improvements, the need to specify waste origin/ownership, and identifying preferred transportation routing. This paper will discuss various aspects of the NNSS WAC, and how it has evolved based upon DOE orders and requirements, acceptance of offsite mixed wastes, acceptance of non-radioactive classified wastes, stakeholder concerns and agreements, and the changing face of the remaining legacy waste at DOE sites across the United States. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Development of Kinetics and Mathematical Models for High Pressure Gasification of Lignite-Switchgrass Blends

The overall objective of the current project was to investigate the high pressure gasification characteristics of a feed containing both coal and biomass. The two feed types differ in their ash contents and ash composition, particularly the alkali content. Gasification of a combined feed of coal and biomass has the potential for considerable synergies that might lead to a dramatic improvement in process economics and flexibility. The proposed study aimed to develop a detailed understanding of the chemistry, kinetics, and transport effects during high pressure gasification of coal-biomass blend feed. Specifically, we studied to develop: (a) an understanding of the catalytic effect of alkali and other inorganic species present in the biomass and coal, (b) an understanding of processing conditions under which synergistic effects of the blending of coal and biomass might be observed. This included the role of particle size, residence time, and proximity of the two feed types, (c) kinetics of high pressure gasification of individual feeds as well as the blends, and (d) development of mathematical models that incorporate kinetics and transport models to enable prediction of gasification rate at a given set of operating conditions, and (e) protocols to extend the results to other feed resources. The goal was to provide a fundamental understanding of the gasification process and guide in optimizing the configurations and design of the next generation of gasifiers. The approach undertaken was centered on two basic premises: (1) the gasification for small particles without internal mass transfer limitations can be treated as the sum of two processes in series (pyrolysis and char gasification) , and (2) the reactivity of the char generated during pyrolysis not only depends on the pressure and temperature but is also affected by the heating rates. Thus low heating rates (10-50 °C/min) typical of PTGA fail to produce char that would typically be formed at high heating rates (~10 4 °C/sec), encountered in entrained flow gasifiers. The char morphology, also a function of the heating rate, would influence the transport rates during the char gasification phase. Thus, heating rate plays a critical role through which both, pyrolysis and char gasification, are interconnected. We utilized two complementary gasification experiments: PEFR (pressurized entrained flow gasifier) and PTGA (pressurized thermo-gravimetric analyzer). The PEFR allowed us to study gasification at pressures, temperatures, and heating rates relevant for coal-biomass gasifiers. The PTGA work was useful in understanding the basic chemistry of the evolution of various gaseous species during pyrolysis. These results helped improved our understanding of the chemistry and chemical changes during pyrolysis. The role alkali metals and other inorganics in char gasification using steam and/or CO 2 was investigated. Finally, the mathematical models for char gasification without the transport effects were developed at commercial operating conditions.

09 BIOMASS FUELS↗
DOE OSTI2021

New insights into fuel blending effects: Intermolecular chemical kinetic interactions affecting autoignition times and intermediate-temperature heat release

Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.

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

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The dynamical state of eROSITA clusters and its impact on the brightest cluster galaxy luminosity

The first Spectrum-Roentgen-Gamma (SRG) eROSITA public release contains 12 247 clusters and groups from its first 6 months of operation. We used the offset between the brightest cluster galaxy (BCG) and the X-ray peak ( D BCG − X ) to classify the cluster dynamical state of 3946 galaxy clusters and groups. We aim to investigate the evolution of the merger and relaxed cluster distributions with redshift and mass, and the distributions’ impact on the BCG. We used the X-ray peak from the eROSITA survey and the BCG position from the LS DR10 optical data, which includes the DECam eROSITA Survey optical data, to measure the D BCG − X offset. We modelled the distribution of D BCG − X , in units of R 500 , as the sum of two Rayleigh distributions representing the cluster’s relaxed and disturbed populations, and explored their evolution with redshift and mass. To explore the impact of the cluster’s dynamical state on the BCG luminosity, we separated the main sample according to the dynamical state. We defined clusters as relaxed if D BCG − X < 0.25R500, disturbed if D BCG − X > 0.5R500, and ‘diverse’ if 0.25R500 < D BCG − X < 0.5R500. We find no evolution of the merging fraction with redshift or mass. The width of the relaxed distribution increases with redshift, while the width of the two Rayleigh distributions decreases with mass. The analysis reveals that BCGs in relaxed clusters are brighter than BCGs in both the disturbed and diverse cluster populations. The most significant differences are found for high-mass clusters at higher redshifts. The results suggest that BCGs in low-mass clusters are less centrally bound than those in high-mass systems, irrespective of the dynamical state. Over time, BCGs in relaxed clusters progressively align with the potential centre. This alignment correlates with their luminosity growth relative to BCGs in dynamically disturbed clusters, underscoring the critical role of the cluster’s dynamical state in regulating BCG evolution.

galaxy clusters↗
DOE OSTI2016

Organ and Tissue-specific Sucrose Transporters. Important Hubs in Gene and Metabolite Networks Regulating Carbon Use in Wood-forming Tissues of Populus

The overall project objective was to probe the relationship between sucrose transporters and plant productivity in the biomass for biofuels woody perennial, Populus. At the time the proposal was written, sucrose transporters had already been investigated in many plant model systems, primarily with respect to the export of photosynthate sucrose from source leaves, and the uptake of sucrose in storage organs and seeds. Preliminary findings by the PI found that in Populus, sucrose transporter genes (SUTs) were well expressed in wood-forming tissues that comprise the feedstock for biofuels production. Because sucrose comprises by far the predominant form in which photosynthate is delivered from source organs to sink organs like roots and wood-forming tissues, SUTs control a gate that nominally at least could impact the allocation or partitioning of sucrose for potentially competing end uses like growth (stem biomass) and storage. In addition, water use might be conditioned by the way in which sucrose is distributed throughout the plant, and/or by the way in which sucrose is partitioned intracellularly. Several dozen transgenic lines were produced in year 1 of the project to perturb the expression ratio of multiple plasma membrane (PM) SUTs (intercellular trafficking), versus the single tonoplast membrane (TM) sucrose transporter that effectively regulates intracellular trafficking of sucrose. It was possible to obtain transgenic lines with dual SUT gene knockdown using the 35S promoter, but not the wood-specific TUA1 promoter. By the end of project year 2, a decision was made to work with the 35S plants while archiving the TUA1 plants. The PhD candidate charged with producing the transgenic lines abandoned the project during its second year, substantially contributing to the decision to operate with just the 35S lines. That student’s interests ranged more toward evolutionary topics, and a report on SUT gene evolution was published (Peng et al 2014).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2016

International Approaches for Nuclear Waste Disposal in Geological Formations: Report on Fifth Worldwide Review

The goal of the Fifth Worldwide Review is to document evolution in the state-of-the-art of approaches for nuclear waste disposal in geological formations since the Fourth Worldwide Review that was released in 2006. The last ten years since the previous Worldwide Review has seen major developments in a number of nations throughout the world pursuing geological disposal programs, both in preparing and reviewing safety cases for the operational and long-term safety of proposed and operating repositories. The countries that are approaching implementation of geological disposal will increasingly focus on the feasibility of safely constructing and operating their repositories in short- and long terms on the basis existing regulations. The WWR-5 will also address a number of specific technical issues in safety case development along with the interplay among stakeholder concerns, technical feasibility, engineering design issues, and operational and post-closure safety. Preparation and publication of the Fifth Worldwide Review on nuclear waste disposal facilitates assessing the lessons learned and developing future cooperation between the countries. The Report provides scientific and technical experiences on preparing for and developing scientific and technical bases for nuclear waste disposal in deep geologic repositories in terms of requirements, societal expectations and the adequacy of cases for long-term repository safety. The Chapters include potential issues that may arise as repository programs mature, and identify techniques that demonstrate the safety cases and aid in promoting and gaining societal confidence. The report will also be used to exchange experience with other fields of industry and technology, in which concepts similar to the design and safety cases are applied, as well to facilitate the public perception and understanding of the safety of the disposal approaches relative to risks that may increase over long times frames in the absence of a successful implementation of final dispositioning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗