Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model Project”

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 487 records · Page 27

Chromatographic Modeling of 1-Foot Detritiation Test Column Results for the Optimization of Regeneration Time

This report documents the methods and results of a modeling study on a series of data generated by 1-foot test columns for the Water Detritiation project. In order to compliment modeling work already in progress on this project, a purely empirical curve-fitting approach is used in this study to determine chromatographic performance indicators such as retention time, curve width, and the number of chromatographic stages. The data used in this study comes from four sets of chromatographic test separations performed in March through May of 2017 using 1-foot columns. All are using the same flow rates and temperature (room temperature, nominally 25 C°). The first set of data (referred to in this report as the Column 1 data) is a series of twelve breakthrough curve measurements. Prior to each of these twelve runs, the column was regenerated by flowing hydrogen gas through the column, for increasing amounts of time, ranging from 2 hours to 25 hours. The second set of data (referred to as the Column 2 data) is comprised of five measurements of breakthrough curves measured in similar fashion as the Column 1 data, across a range of regeneration times from 2 to 21 hours. The third set of data (referred to as the Purge data) is of three breakthrough curves, measured after a regeneration time of 8 hours, which takes place after a nitrogen purge of 15, 30, or 45 minutes. The fourth data set is a series of five breakthrough curve measurements using a 2 nd batch of stationary phase, taken after regeneration times ranging from 2 to 20 hours. The overall modeling approach is as follows. First, the data is pre-processed for formatting, artifact removal, and is normalized. The data is then fit by integrated versions of four common chromatographic models. The best model is chosen based on the predictive performance of the model, and chromatographic performance values are calculated based on the best model’s fit parameters. Each of these processes are detailed below.

42 ENGINEERING↗

Data-driven projection pursuit adaptation of polynomial chaos expansions for dependent high-dimensional parameters

Uncertainty quantification (UQ) and inference involving a large number of parameters are valuable tools for problems associated with heterogeneous and non-stationary behaviors. The difficulty with these problems is exacerbated when these parameters are statistically dependent requiring statistical characterization over joint measures. Probabilistic modeling methodologies stand as effective tools in the realms of UQ and inference. Among these, polynomial chaos expansions (PCE), when adapted to low-dimensional quantities of interest (QoI), provide effective yet accurate approximations for these QoI in terms of an adapted orthogonal basis. These adaptation techniques have been cast as projection pursuits in Gaussian Hilbert space in what has been referred to as a projection pursuit adaptation (PPA) by Xiaoshu Zeng and Roger Ghanem (2023). The PPA method efficiently identifies an optimal low-dimensional space for representing the QoI and simultaneously evaluates an optimal PCE within that space. The quality of this approximation clearly depends on the size of the training dataset, which is typically a function of the adapted reduced dimension. Here, the complexity of the problem is thus mediated by the complexity of the low-dimensional quantity of interest and not the complexity of the high-dimensional parameter space.

Data-driven↗

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗

Clusters have edges: the projected phase-space structure of SDSS redMaPPer clusters

ABSTRACT We study the distribution of line-of-sight velocities of galaxies in the vicinity of Sloan Digital Sky Survey (SDSS) red-sequence Matched-filter Probabilistic Percolation (redMaPPer) galaxy clusters. Based on their velocities, galaxies can be split into two categories: galaxies that are dynamically associated with the cluster, and random line-of-sight projections. Both the fraction of galaxies associated with the galaxy clusters, and the velocity dispersion of the same, exhibit a sharp feature as a function of radius. The feature occurs at a radial scale Redge ≈ 2.2Rλ, where Rλ is the cluster radius assigned by redMaPPer. We refer to Redge as the ‘edge radius’. These results are naturally explained by a model that further splits the galaxies dynamically associated with a galaxy cluster into a component of galaxies orbiting the halo and an infalling galaxy component. The edge radius Redge constitutes a true ‘cluster edge’, in the sense that no orbiting structures exist past this radius. A companion paper tests whether the ‘halo edge’ hypothesis holds when investigating the full three-dimensional phase-space distribution of dark matter substructures in numerical simulations, and demonstrates that this radius coincides with a suitably defined splashback radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing transmit field inhomogeneity of parallel RF transmit design in 7T MRI using deep learning

Ultrahigh field (UHF) Magnetic Resonance Imaging (MRI) provides a higher signal-to-noise ratio and, thereby, higher spatial resolution. However, UHF MRI introduces challenges such as transmit radiofrequency (RF) field (B+1) inhomogeneities, leading to uneven flip angles and image intensity anomalies. These issues can significantly degrade imaging quality and its medical applications. This study addresses B+1 field homogeneity through a novel deep learning-based strategy. Traditional methods like Magnitude Least Squares (MLS) optimization have been effective but are time-consuming and dependent on the patient’s presence. Recent machine learning approaches, such as RF Shim Prediction by Iteratively Projected Ridge Regression and deep learning frameworks, have shown promise but face limitations like extensive training times and oversimplified architectures. We propose a two-step deep learning strategy. First, we obtain the desired reference RF shimming weights from multi-channel B+1 fields using random-initialized Adaptive Moment Estimation. Then, we employ Residual Networks (ResNets) to train a model that maps B+1 fields to target RF shimming outputs. Our approach does not rely on pre-calculated reference optimizations for the testing process and efficiently learns residual functions. Comparative studies with traditional MLS optimization demonstrate our method’s advantages in terms of speed and accuracy. The proposed strategy achieves a faster and more efficient RF shimming design, significantly improving imaging quality at UHF. This advancement holds potential for broader applications in medical imaging and diagnostics.

Lu, Zhengyi [Vanderbilt University]↗

NGEE Arctic Plant Traits: Vegetation Plot Locations, Ecotypes, and Photos, Kougarok Road Mile Marker 64, Seward Peninsula, Alaska, 2016

This dataset provides field site documentation for related datasets including GPS waypoints for intensive plots, reference points, vegetation plots, and soil temperature/moisture monitoring stations that were established in July 2016 at the Kougarok hill slope located at Kougarok Road, Mile Marker 64. At all intensive plots and reference points, photographs of the surrounding area were also taken using an Olympus Stylus TG-3 digital camera. Photographs were July 2nd through July 7th looking in each of the four cardinal directions. GPS waypoints were taken using a Garmin eTrek 10 handheld GPS and World Geodetic System 1984 (WGS-84). The error associated with the Garmin GPS coordinates was 9 foot or approximately 3 m. When marking waypoints, users waited until the instrument reached this error level before recording. Coordinate units are decimal degrees. This dataset is comprised of one comma-separated (*.csv) file containing GPS waypoints locating intensive plots, reference points, soil temperature/moisture monitoring stations, and vegetation plots. The dataset also contains one folder of photographs (*.jpg) taken at intensive stations, one folder of photographs (*.jpg) taken at reference points, and one User file (*.pdf). The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

Solar Resource Calibration, Measurement, and Dissemination: Final Report FY 2016-FY 2018

The Solar Resource Calibration, Measurement, and Dissemination project supports the U.S. Department of Energy's Solar Energy Technologies Office (SETO) initiative to make solar energy cost-competitive with other forms of electricity by improving the tools and methods to measure and model solar radiation and therefore reduce uncertainty in predicting solar output and improve the bankability of solar projects. This project has three tasks that conduct research on advancing solar resource measurements and modeling: Task 1: Applied solar Radiation Measurements. Task 2: National Solar Radiation Data Base (NSRDB). Task 3: Knowledge Sharing. The National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) activities ensure traceable solar measurements throughout the United States that are essential for project feasibility, due diligence, financing, and plant operations. The SRRL houses NREL's Baseline Measurement System (BMS), which provides a platform for improving the accuracy of solar measurements through research to improve the calibration and characterization of instruments and enable the development of new measurement technology and standards. The platform contains a variety of instruments measuring solar radiation and other meteorological parameters. Following best practices, BMS instruments are calibrated to ensure traceability to the World Radiometric Reference in compliance with ISO/IEC 17025 accreditation requirements. The BMS data sets are quality-controlled and made publicly available through the Measurement and Instrumentation Data Center website. The NSRDB is the most accessed public database of solar irradiance and associated weather parameters for use in energy modeling. It supports SETO's goals of reducing barriers to high penetrations of solar technologies by providing easy access to high-quality, foundational data essential for innovative product development and downstream modeling through the NSRDB website. Data are currently available for the contiguous United States and beyond for 20 years covering 1998-2017. The Physical Solar Model (PSM) that underpins the current NSRDB data sets is a unique physics-based model that has been developed at NREL and has opened the door to the use of next-generation satellite data sets for solar resource assessment and forecasting. Through the American Society for Testing Materials (ASTM) International and other international standards organizations, NREL continues to provide leadership in the development of standards that are relevant to the measurement and modeling of solar radiation for solar energy applications. Standards developed cover the calibration method of pyrheliometers (direct normal irradiance measurements), pyranometers (global horizontal irradiance measurements), spectroradiometers, ultraviolet radiometers, correct field deployment of these instruments, spectral distribution reference tables for photovoltaic modeling, and transfer of indoor calibrations. Further, NREL leads the International Energy Agency (IEA) Photovoltaic Power Systems Programme Task 16: Subtask 1, which focuses on the evaluation of current and emerging resource assessment methodologies. NREL and the IEA updated the Best Practices Handbook for the Collection and Use of Solar Resource Data for Solar Energy Applications, which is widely used by the solar industry.

14 SOLAR ENERGY↗

Empowering Geothermal Research: The Geothermal Data Repository's New AI Research Assistant

The Department of Energy's (DOE) Geothermal Data Repository (GDR) team has integrated a Large Language Model (LLM) with the metadata and supporting documents associated with GDR datasets to create an Artificially Intelligent (AI) research assistant. By leveraging work done to make GDR metadata machine-readable and an open-source LLM integration model called the Energy Language Model, developed by the National Renewable Energy Laboratory, AskGDR serves as a virtual research assistant to GDR users. It provides answers to a variety of user-provided questions using natural language processing and generative machine learning. Users can get answers to questions about specific datasets, including inquiries about the equipment, assumptions and methodologies used in the origination of the data; or more abstract questions, such as the applicability of data to specific research fields. AskGDR improves the discoverability of geothermal data by helping guide users to datasets beyond simple keyword searches. It enables users to find data based on properties of the data, discover information contained within supporting documents, and explore data from projects related to their research objectives. This paper will outline the development, integration, output, and efficacy of the AskGDR LLM, including adherence to scientific rigor through improvements designed to increase the accuracy of generated answers, avoid speculation, and provide proper references for all resources used.

access↗

Modelling of Wastewater Heat Recovery Heat Pump Systems

Wastewater heat recovery is currently an underutilized technology that could be part of solving the climate crisis. A large portion of the heat that leaves a building in the form of wastewater is potentially recoverable for pre-heating domestic hot water or other service water systems. While there are several different approaches to wastewater heat recovery, this project focused on creating detailed, integrated building models for wastewater heat recovery heat pump systems. EnergyPlus models were developed featuring inputs and assumptions corresponding to manufacturers’ specifications, performance lab test data and feedback from engineering consultants. EnergyPlus’s supervisory control Energy Management System objects were heavily relied upon to overcome modelling challenges. The developed EnergyPlus model was integrated into U.S. Department of Energy New Construction Reference Building models for various climate zones and building types to assess potential energy use, energy cost and greenhouse gas emission reductions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Crystal structures of Salmonella enterica FraB deglycase reveal a conformational heterodimer with remarkable structural plasticity at the active site

Abstract Thefralocus ofSalmonella entericaencodes five genes for metabolism of fructose‐asparagine, an Amadori product formed by condensation of asparagine with glucose. In the last step of this pathway, the FraB deglycase cleaves 6‐phospho‐fructose‐aspartate into glucose‐6‐phosphate and aspartate. In homology models, FraB forms a homodimer with two equivalent active sites located at the dimer interface. E214 and H230, two invariant residues essential for catalysis, project into each active site cleft from opposing subunits of the dimer. Here, we have determined six crystal structures of FraB, three of a variant containing an N‐terminal His 6 tag and two mutations needed for crystallization (hereafter referred to as WT′), two with additional mutations to active site residues (E214A and P232A), and one of a variant with C‐terminal residues 313–325 deleted. Surprisingly, in the WT′ FraB structure, the two catalytic residues, E214 (general base) and H230 (general acid), are positioned ~22 Å apart. In the E214A and C‐terminus‐truncated FraB variants, however, a conformational change in the E214‐residing helix brings E214 and H230* to ~7 Å (* indicates residue from the second protomer that creates the inter‐subunit catalytic center). The loop bearing H230 also exhibits significant variation, ranging from being completely disordered to adopting open or closed states, with the nearby P232* residue being eithercisortrans. The C‐terminal residues 313–325 form a flexible “C‐tail” that can be fully disordered, bind in the active site to block access of substrate, or angle across the active site to wrap across the other subunit of the dimer and potentially close over substrate. Collectively, these structures reveal that FraB is a conformational heterodimer with two chemically identical subunits that are constrained to adopt different structures as they come together for catalysis. This plasticity likely involves correlated opening and closure of the two active sites for their respective binding and release of substrates and ligands.

Biochemistry & Molecular Biology↗

The Implications of Global Change for the Co-Evolution of Argentina's Integrated Energy-Water-Land Systems

This study seeks to understand how Argentina’s energy, water, and land (EWL) systems will co-evolve under a representative array of human and earth system influences, including socioeconomic change, climate change, and climate policy. To capture Argentina’s sub-national EWL dynamics in the context of global change, we couple the Global Change Analysis Model with a suite of consistent, gridded sectoral downscaling models to explore multiple stakeholder-engaged scenarios. Across scenarios, Argentina has the economic opportunity to use its vast land resources to satisfy growing domestic and international demand for crops, such as oil (e.g., soy) and biomass. The human (rather than earth) system produces the most dominant changes in mid-century EWL resource use. A Reference scenario characterized by modest socioeconomic growth projects a 40% increase in Argentina’s agricultural production by 2050 (relative to 2020) by using 50,000 km 2 of additional cropland and 40% more water. A Climate Policy scenario designed to achieve net-zero carbon emissions globally shortly after mid-century projects that Argentina could use 100,000 km 2 of additional land (and 65% more water) to grow biomass and other crops. The burden of navigating these national opportunities and challenges could fall disproportionately on a subset of Argentina’s river basins. The Colorado and Negro basins could experience moderate-to-severe water scarcity as they simultaneously navigate substantial irrigated crop demand growth and climate-induced declines in natural water availability. Argentina serves as a generalizable testbed to demonstrate that multi-scale EWL planning challenges can be identified and managed more effectively via integrated analysis of coupled human-earth systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Simultaneous 3D–2D image registration and C-arm calibration: Application to endovascular image-guided interventions

Purpose: Three-dimensional to two-dimensional (3D–2D) image registration is a key to fusion and simultaneous visualization of valuable information contained in 3D pre-interventional and 2D intra-interventional images with the final goal of image guidance of a procedure. In this paper, the authors focus on 3D–2D image registration within the context of intracranial endovascular image-guided interventions (EIGIs), where the 3D and 2D images are generally acquired with the same C-arm system. The accuracy and robustness of any 3D–2D registration method, to be used in a clinical setting, is influenced by (1) the method itself, (2) uncertainty of initial pose of the 3D image from which registration starts, (3) uncertainty of C-arm’s geometry and pose, and (4) the number of 2D intra-interventional images used for registration, which is generally one and at most two. The study of these influences requires rigorous and objective validation of any 3D–2D registration method against a highly accurate reference or “gold standard” registration, performed on clinical image datasets acquired in the context of the intervention. Methods: The registration process is split into two sequential, i.e., initial and final, registration stages. The initial stage is either machine-based or template matching. The latter aims to reduce possibly large in-plane translation errors by matching a projection of the 3D vessel model and 2D image. In the final registration stage, four state-of-the-art intrinsic image-based 3D–2D registration methods, which involve simultaneous refinement of rigid-body and C-arm parameters, are evaluated. For objective validation, the authors acquired an image database of 15 patients undergoing cerebral EIGI, for which accurate gold standard registrations were established by fiducial marker coregistration. Results: Based on target registration error, the obtained success rates of 3D to a single 2D image registration after initial machine-based and template matching and final registration involving C-arm calibration were 36%, 73%, and 93%, respectively, while registration accuracy of 0.59 mm was the best after final registration. By compensating in-plane translation errors by initial template matching, the success rates achieved after the final stage improved consistently for all methods, especially if C-arm calibration was performed simultaneously with the 3D–2D image registration. Conclusions: Because the tested methods perform simultaneous C-arm calibration and 3D–2D registration based solely on anatomical information, they have a high potential for automation and thus for an immediate integration into current interventional workflow. One of the authors’ main contributions is also comprehensive and representative validation performed under realistic conditions as encountered during cerebral EIGI.

60 APPLIED LIFE SCIENCES↗

Analysis of Subject Matter Expert (SME) Feedback on the First Version of the Biosurveillance Analytics Resource directory (BaRD)

As part of a recently funded project, LANL was tasked with the development of the BaRD, a reference tool in support of the Biosurveillance Ecosystem (BSVE) that will serve to enhance situational awareness during an infectious disease outbreak and provide timely and relevant information to support effective decision making to prevent, stop or control an outbreak or epidemic. The BaRD is being designed to enable rapid selection of appropriate models for use in disease prediction, monitoring, or forecasting at various population levels. It is a searchable, relational database that catalogs and classifies model specific information to allow facile reporting on tools of interest.

60 APPLIED LIFE SCIENCES↗

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Commercial Building Technology Demonstration (Final Scientific/Technical Report)

The background for this project was a unique and innovative concept that had been proposed and demonstrated by enVerid Systems, reducing the energy consumption in HVAC systems in commercial and public buildings by up to 50% at peak loads. The concept involves selective scrubbing of gas contaminants from indoor air as a practical substitute for outside air ventilation, the latter being the conventional practice today. This scrubbing enables recycling of most of the indoor air and significantly decreasing the amount of outside air, resulting in a reduced HVAC load and commensurate energy savings for the facility. Key to our approach is low cost scrubbing of all gas contaminants from indoor air using novel, efficient and regenerable sorbent materials. The total solution, labeled HLR (HVAC Load Reduction) is implemented by means of compact modules that can be retrofitted onto virtually any building, and utilizes novel sorbents that have been developed and brought to scale production by enVerid. Although the technology had been deployed and demonstrated, the US commercial building and HVAC markets, famously conservative with respect to innovation and new technologies, had yet to embrace the idea or the product. The overarching goal of the project was to help overcome market hesitancy and enable widespread adoption, by creating a critical mass of success stories across different regions and building types, all under oversight and verification of the DOE and NREL, to ensure credibility and accuracy. The project, overall, was a resounding success. Several installations were deployed and tested comprehensively for the impact on both energy consumption and indoor air quality, and results were not only excellent, but quantitatively in line with expectations. A total number of eight (8) installations were attempted, with 2 aborted due to building issues. Of the six successful installations, however, complete M&V was performed by NREL on 2 (Miami and NYC) and partial analysis on a third site (ArcBest), and these three, along with their results, were therefore presented from a financial and reporting standpoint. Although the number of case studies fell short of the original target of 10, we nevertheless had excellent demonstrations in vastly different climate zones: Southeast, Northeast and Midwest, and in very different building types. This included a modern high-rise office building in New York City and a University Campus health center in Miami, and a mid-rise high-density office in Arkansas. The cooling load reduction was shown to be significant, averaging 37% during the final cooling season in the Miami site. In the NYC site, the cooling energy savings were generally less, at around 12%, as expected. While this is a lower number for average/total savings, two comments are important to note: (i) the outside airflow (OA) reduction was not consistent during the data collection period, at least in part understating the attainable savings through incorrect measurement of the “high energy” reference state; and (ii) the peak impact in NYC is almost certainly as high as Miami, with critical ramifications for peak load reduction and equipment downsizing ability. Importantly, the measured energy savings in all locations were related to the outdoor climate conditions and the actual achievable reduction in outside air, very much in line with models and calculations; this leaves us with increased confidence in our understanding of the interaction between the technology and the building, and our ability to project the impact of this innovation in the future, through broad adoption. Simply put: the energy savings and load reduction impact of the HLR solution has been confirmed quantitatively. The results also confirmed the air quality objective. This may be even more significant, in that air quality is more complex and harder to measure and quantify, and the performance of the HLR as an air scrubber is the linchpin of the entire solution. In all sites, indoor air quality met or exceeded the required parameters, in terms of achieving the target (low) levels of indoor pollutants. Furthermore, the concentrations of outdoor sourced pollutants were markedly lower, as expected, mostly due to the reduction in outside air intake although in part also due to the scrubbing. The success was not immediate nor easy, and early challenges with the systems, especially (but not only) around controls and communication, caused delays in commissioning and undermined some of our M&V schedules. The good news is no fundamental issues were encountered and the number of such glitches was declining quickly by the end of the project. Other challenges had to do with coordinating the control of the building OA dampers, which the HLR needs to have in order to deliver the load reduction, but which the facility manager controls regularly and independently. This is not necessarily a problem for normal HLR operation but interferes with accurate M&V, where switching OA settings back and forth on a precise timetable is essential for clean data. The project provided valuable experience in retrofitting buildings that had been originally designed to operate with less recirculation and more outside air, and some key learnings related to the cyclical regeneration of the sorbents in the operational environment of the building. Some of these learnings gave already been incorporated into ongoing improvements in the design of the HLR module as well as the software that is used to control its operation. The cost and time of retrofit installations is much better understood; so is the test and balance procedure of the building that determines the available range of outside air reduction, which in turn determines the amount of load reduced. We had also implemented at least one annual cycle of sorbent cartridge replacement in the field, after which we tested the used sorbent to assess the rate of degradation of sorbent potency – an important factor affecting the long-term economics of HLR deployments. A notable shortfall of the project, other than the fact that the number of sites was below target, was that we had not (yet) completed testing to measure and verify (M&V) the winter energy savings. This is primarily due to late start on the NYC sites (Miami does not portend any heating savings associated with reduced outside air), combined with the issues of OA damper control. It is by no means suggested that such savings cannot be captured, only that the M&V with respect to winter in NYC has been incomplete. The value of this project and its interest to the general public are hard to overstate. The case studies created from Miami and especially from NYC have, without a doubt, opened up the interest of the broad HVAC industry to HLR in the past year. Major manufacturer reps have signed up to represent and promote the product, and dozens of top tier engineering consultants have introduced HLR as a base of design in upcoming construction and renovation projects – precisely the overarching goal of this project. These would not have happened without the highly visible installation in NYC and its ultimately satisfied customer. Furthermore, the DOE itself issued a report in late 2017 analyzing hundreds of candidate technologies for reducing energy consumption in HVAC; remarkably, the budding HLR solution ranked among the top 3(!); no doubt, its status as a viable and significant contender on this list is informed by the success of these demonstrations. And furthermore, in the wake of these demos, several major energy utilities, led by ConEdison in NYC, have endorsed the HLR solution and now routinely provide substantial cash rebates for customers choosing to incorporate the HLR in their buildings and reap the benefit of the reduced peak load as well as reduced cumulative energy demand. Last but not least, in January of 2019, as this report was being written, the HLR module won the AHR Expo Innovation Award for Green Buildings and, remarkably, the prestigious and only Product of the Year 2019, being selected by an independent professional panel of judges from AHR and ASHRAE. Given the confirmation of the enormous energy saving potential of HLR in widescale deployment, as well as ancillary but significant benefits in terms of equipment downsizing and indoor air quality, the benefit of this effort to the public speaks for itself. Accelerating adoption of HLR, without dependence on subsidies and driven by the private sector, is well under way. There is still much work to be done in terms of educating the HVAC and construction ecosystems and improving the product’s robustness and simplicity, to further facilitate the adoption of this technology and the capture of its potential value. But the American public will greatly benefit both directly and indirectly in the immediate future and beyond.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Verification & Validation of High-Order Short-Characteristics-Based Deterministic Transport Methodology on Unstructured Grids

The research team has developed a practical, high-order, discrete-ordinates, short characteristics neutron transport code for three-dimensional configurations represented on unstructured tetrahedral grids that can be used for realistic reactor physics applications at both the assembly and core levels. This project will perform a comprehensive verification and validation of this new computational tool against both a continuous-energy Monte Carlo simulation (e.g. MCNP) and experimentally measured data, an essential prerequisite for its deployment in reactor core modeling. Verification is divided into three phases. The team will first conduct spatial mesh and expansion order refinement studies to monitor convergence of the numerical solution to reference solutions. This is quantified by convergence rates that are based on integral error norms computed from the cell-by-cell difference between the code’s numerical solution and its reference counterpart. The latter is either analytic or very fine- mesh numerical solutions from independent computational tools. For the second phase, the team will create a suite of code-independent benchmark configurations to enable testing the theoretical order of accuracy of any particular discretization of the discrete ordinates approximation of the transport equation. For each tested case (i.e. mesh and spatial approximation order), researchers will execute the code and compare the resulting numerical solution to the exact solution on a per cell basis to determine the distribution of the numerical error. The final activity comprises a comparison to continuous-energy Monte Carlo solutions for zero-power critical configuration measurements at Idaho National Laboratory’s Advanced Test Reactor (ATR). Results of this comparison will allow the investigators to distinguish between modeling errors and the above-listed discretization errors introduced by the deterministic method, and to separate the sources of uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗