SEARCH · Engineering Papers

Results for “variable selection”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Automated characterization of perceptual quality of clinical chest radiographs: Validation and calibration to observer preference

Purpose: The authors previously proposed an image-based technique [Y. Lin et al. Med. Phys. 39, 7019–7031 (2012)] to assess the perceptual quality of clinical chest radiographs. In this study, an observer study was designed and conducted to validate the output of the program against rankings by expert radiologists and to establish the ranges of the output values that reflect the acceptable image appearance so the program output can be used for image quality optimization and tracking. Methods: Using an IRB-approved protocol, 2500 clinical chest radiographs (PA/AP) were collected from our clinical operation. The images were processed through our perceptual quality assessment program to measure their appearance in terms of ten metrics of perceptual image quality: lung gray level, lung detail, lung noise, rib–lung contrast, rib sharpness, mediastinum detail, mediastinum noise, mediastinum alignment, subdiaphragm–lung contrast, and subdiaphragm area. From the results, for each targeted appearance attribute/metric, 18 images were selected such that the images presented a relatively constant appearance with respect to all metrics except the targeted one. The images were then incorporated into a graphical user interface, which displayed them into three panels of six in a random order. Using a DICOM calibrated diagnostic display workstation and under low ambient lighting conditions, each of five participating attending chest radiologists was tasked to spatially order the images based only on the targeted appearance attribute regardless of the other qualities. Once ordered, the observer also indicated the range of image appearances that he/she considered clinically acceptable. The observer data were analyzed in terms of the correlations between the observer and algorithmic rankings and interobserver variability. An observer-averaged acceptable image appearance was also statistically derived for each quality attribute based on the collected individual acceptable ranges. Results: The observer study indicated that, for each image quality attribute, the averaged observer ranking strongly correlated with the algorithmic ranking (linear correlation coefficient R > 0.92), with highest correlation (R = 1) for lung gray level and the lowest (R = 0.92) for mediastinum noise. There was a strong concordance between the observers in terms of their rankings (i.e., Kendall’s tau agreement > 0.84). The observers also generally indicated similar tolerance and preference levels in terms of acceptable ranges, as 85% of the values were close to the overall tolerance or preference levels and the differences were smaller than 0.15. Conclusions: The observer study indicates that the previously proposed technique provides a robust reflection of the perceptual image quality in clinical images. The results established the range of algorithmic outputs for each metric that can be used to quantitatively assess and qualify the appearance quality of clinical chest radiographs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

SU-E-T-369: Evaluating Intensity Modulated Proton Therapy Relative to Passive Scattering Proton Therapy for Increased Vertebral Column Sparing in CSI of Pediatric Patients

Purpose: To investigate the trade-off between vertebral column sparing and thecal-sac target coverage in craniospinal irradiation (CSI) of pediatric patients treated with passive-scattering (PS) and intensity modulated (IMPT) proton therapy. Methods: We selected 2 pediatric patients treated with PS CSI for medulloblastoma. Spinal irradiation was re-planned with IMPT. For all cases, we assumed prescription dose of 23.4 Gy(RBE), with the spinal canal receiving at least 95% of 23.4 Gy(RBE). PS planning was performed using the commercial system XiO. IMPT planning was done using the Astroid planning system. Beam arrangements consisted of (a) PS posterior-anterior (PA) field, PS-PA, (b) IMPT PA field, IMPT-PA, and (c) two posterior oblique IMPT fields, IMPT2 (-35°, 35°). Dose distributions were re-calculated using TOPAS Monte Carlo, along with LET distributions, to investigate LET variations within the target and vertebra anatomy. Variable RBE-weighed dose distributions were also calculated based on a dose and LET-dependent biophysical model. Dosimetric data were compared among the plans for the target volume, spinal cord and adjacent critical organs (thecal-sac and cauda equina). Results: IMPT2 resulted in better sparing of the posterior vertebral column (entrance region posterior to thecal-sac), where planned dose was approximately 6–8Gy(RBE). For IMPT-PA and PS-PA the MC-calculated dose to the posterior vertebral column was, on average, 20Gy and 18Gy respectively. For IMPT2 higher mean-LET (5keV/µm/(g/cm3)) values were observed in anterior vertebral column (beyond the thecal-sac) relative to IMPT-PA and PS-PA, where mean-LET was 3.5keV/µm/(g/cm3) and 2.5keV/µm/(g/cm3) respectively. The higher LET region observed for both IMPT plans was in the distal end of treatment fields, where dose delivered was less 5Gy(RBE). Conclusion: The two-oblique proton beams IMPT2 best spared the spinal column, while reducing the dose to the posterior spinal column from 18–20 to 6–8 Gy(RBE). The best LET distribution was obtained with the PS-PA fields.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Watershed zonation through hillslope clustering for tractably quantifying above- and below-ground watershed heterogeneity and functions

Abstract. In this study, we develop a watershed zonation approach for characterizing watershed organization and functions in a tractable manner by integrating multiple spatial data layers. We hypothesize that (1) a hillslope is an appropriate unit for capturing the watershed-scale heterogeneity of key bedrock-through-canopy properties and for quantifying the co-variability of these properties representing coupled ecohydrological and biogeochemical interactions, (2) remote sensing data layers and clustering methods can be used to identify watershed hillslope zones having the unique distributions of these properties relative to neighboring parcels, and (3) property suites associated with the identified zones can be used to understand zone-based functions, such as response to early snowmelt or drought and solute exports to the river. We demonstrate this concept using unsupervised clustering methods that synthesize airborne remote sensing data (lidar, hyperspectral, and electromagnetic surveys) along with satellite and streamflow data collected in the East River Watershed, Crested Butte, Colorado, USA. Results show that (1) we can define the scale of hillslopes at which the hillslope-averaged metrics can capture the majority of the overall variability in key properties (such as elevation, net potential annual radiation, and peak snow-water equivalent – SWE), (2) elevation and aspect are independent controls on plant and snow signatures, (3) near-surface bedrock electrical resistivity (top 20 m) and geological structures are significantly correlated with surface topography and plan species distribution, and (4) K-means, hierarchical clustering, and Gaussian mixture clustering methods generate similar zonation patterns across the watershed. Using independently collected data, we show that the identified zones provide information about zone-based watershed functions, including foresummer drought sensitivity and river nitrogen exports. The approach is expected to be applicable to other sites and generally useful for guiding the selection of hillslope-experiment locations and informing model parameterization.

58 GEOSCIENCES↗

DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗

DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

DOE OSTI2024

National Climate Database (NCDB)

The National Climate Database (NCDB) is a high resolution, bias-corrected climate dataset consisting of the three most widely used variables of solar radiation- global horizontal (GHI), direct normal (DNI), and diffuse horizontal irradiance (DHI)- as well as other meteorological data. The goal of the NCDB is to provide unbiased high temporal and spatial resolution climate data needed for renewable energy modeling. The NCDB is modeled using a statistical downscaling approach with Regional Climate Model (RCM)-based climate projections obtained from the North American Coordinated Regional Climate Downscaling Experiment (NA-CORDEX; linked below). Daily climate projections simulated by the Canadian Regional Climate Model 4 (CanRCM4) forced by the second-generation Canadian Earth System Model (CanESM2) for two Representative Concentration Pathways (RCP4.5 or moderate emissions scenario and RCP8.5 or highest baseline emission scenario) are selected as inputs to the statistical downscaling models. The National Solar Radiation Database (NSRDB) is used to build and calibrate statistical models.

Array↗

DOE OSTI2015

Modular Soft X-Ray Emission Spectrograph

We have been developing the Hettrick-Underwood type soft x-ray spectrograph for over a decade. This optical design uses one spherical mirror for focusing x-rays onto the CCD detector and a set of variable line-spacing gratings for energy analysis. Unlike the previous effort carried out at the ALS, what we have achieved in the current design is to incorporate the concept of modularity. The model of spectrograph is shown below. It contains three major components: the optics chamber and downstream section (detector plus driving mechanism using a swing arm pivoted at the grating location, see figure) for inside and outside order (translation stage without swing arm, not shown here). Inside order design is often used when source size is larger than the detector pixel resolution, while the outside order is used when source size is smaller than the detector pixel resolution. Irrespective to which order is selected, there is no change in the mechanism of optics chamber. The location of optics chamber relative to the source and detector can also be changed to optimize the energy resolution at target photon energy. The incident angles on mirror and gratings can be varied, in conjunction with the grating line densities, to balance the resolution and throughput. The internal components of optics chamber are also modular and changing the dimension of optics can be done by using different holders. The modularity not only greatly reduces the design and assembling effort, but also offers the cost benefit. Although the cost of this type of spectrograph is expected to be higher than XES350, it is still significant lower than the one-of-the-kind design used at every synchrotron facilities. The unique flexibility and modularity not only allows the users to tailor the design of spectrograph for the characteristic of their sources, but also to have cost down in standardizing the fabrication parts.

Hussain, Zahid↗

DOE OSTI2022

Hydrogen Storage for Flexible Fossil Fuel Power Generation: Integration of Underground Hydrogen Storage with Gas Turbine (Final Report)

As the nation continues to encourage, through market structures and financial incentives, the proliferation of intermittent renewable electricity, how to optimize the ever-changing electric grid and identify means to retain and improve resilience, while ensuring continued reductions in GHG emissions, will be critical. According to Bloomberg, wind & solar generated 10.5% of US electricity in 2020 and that percentage continues to grow. In support of expanding renewable energy use, and to address its intermittent nature, this project will develop the Hydrogen Storage for Flexible Fossil Fuel Power Generation platform that is dispatchable, reliable, repeatable and have the ability to produce zero or negative carbon power while interfacing with geology capable of CO2 and hydrogen storage. GTI Energy (GTIE) and team members Illinois State Geological Survey (ISGS), Mitsubishi Heavy Industries America (MHIA), Ameren Illinois, Hexagon Purus, and the Low Carbon Resources Initiative (LCRI) completed a Phase I Conceptual Study under contract DE-FE0032012 for Hydrogen Storage for Flexible Fossil Fuel Power Generation: Integration of Underground Hydrogen Storage with Gas Turbine. The Hydrogen Storage for Flexible Fossil Fuel Power Generation platform addresses the intermittent nature of the expanding use of Variable Renewable Energy (VRE) generation. The low cost of the electricity (COE) generated results in greater dispatch and more operation at higher power levels (higher efficiency), fewer short intervals, and fewer start/stop cycles. The reliable, resilient system can produce zero carbon power and store hydrogen. It will demonstrate hydrogen storage in geologic formations like those used in natural gas underground storage thus enabling large scale storage of hydrogen in sedimentary strata across the United States rather than in geographically restricted salt caverns. The Phase I study confirmed the system is feasible and generates power at lower cost than other low carbon approaches. The demonstration defines the pathway for broad commercial application and will accelerate the development of larger systems suitable for centralized utility scale electricity production. The study advanced the maturity of the H 2 storage-based system with flexible power generation by completing a Pre-FEED study (Phase II). The Pre-FEED focused on the selected Energy Farm on the University of Illinois Urbana-Champaign (UIUC) site that includes above ground and underground hydrogen storage, low-carbon hydrogen production (GTI’s Compact Hydrogen Generator, CHG) with underground CO2 sequestration, and a 40-MW class gas turbine. The Pre-FEED addressed the entire system and its interconnection to the natural gas and electric grid and mitigation of key risks, such as storage behavior, load-following, and system operation. During Phase 1 of the project, the team completed key tasks, which moved the entire demonstration project, specific components and approaches closer to commercialization. These Phase I Accomplishments include: Completing System Requirements Review; Completing System Layout and Modeling - Heat & Mass Balance and Process Flow Diagram; Completing modelling of 9 turbine performance cases; Evaluating rock strata for underground storage of hydrogen and sequestration of carbon dioxide; Completing initial modelling of underground storage of hydrogen and withdrawal with evaluation of loss and water production; Identifying roadable storage for above ground hydrogen storage; Identifying existing electrical infrastructure for receiving/delivering electricity; Identifying existing gas supply infrastructure for receiving natural gas; Document concept design/development plans in required reports. Conclusions: The 12-month Feasibility study in Phase I study was completed and confirmed the system is feasible and generates power at lower cost than other low carbon approaches and even lower cost than the reference NGCC plant without carbon capture when taking advantage of 45Q carbon credits. The study enabled the fidelity of the concept to be improved and allowed identification of the requirements for the system. Defining the individual system and component requirements was performed via the system requirements review with the whole team. These requirements were then incorporated into and iterated with our Heat & Mass Balance process model and process flow diagrams were generated to reflect the overall system. This information was then used to complete the TEA and show economic feasibility. Large scale non-salt geologic storage of hydrogen is an enabling technology for a hydrogen-fired turbine that can be retrofitted into large-scale electric generating units (EGU). Our demonstration will include 428 MWh or ~4 hours full load of hydrogen storage (above and underground). Carbon capture inherent to the CHG process can capture 90% CO 2 (with upgrades to >98%). This system provides a COE of 23% savings relative to an NGCC with a post combustion amine system. Our proposed storage system decouples carbon capture and hydrogen production from power production; therefore, we expect our proposed system’s efficiency and variable COE to be superior resulting in overall higher dispatch and reduced deep cycling. Our demonstration will be full to multi-day hydrogen storage and has the potential for longer (seasonal) duration commercially. The demonstration defines the pathway for broad commercial application and will accelerate the development of larger systems suitable for centralized utility scale electricity production.

03 NATURAL GAS↗

DOE OSTI2021

Explaining and predicting human behavior and social dynamics in simulated virtual worlds: reproducibility, generalizability, and robustness of causal discovery methods

Ground Truth program was designed to evaluate social science modeling approaches using simulation test beds with ground truth intentionally and systematically embedded to understand and model complex Human Domain systems and their dynamics Lazer et al. (Science 369:1060–1062, 2020). Our multidisciplinary team of data scientists, statisticians, experts in Artificial Intelligence (AI) and visual analytics had a unique role on the program to investigate accuracy, reproducibility, generalizability, and robustness of the state-of-the-art (SOTA) causal structure learning approaches applied to fully observed and sampled simulated data across virtual worlds. In addition, we analyzed the feasibility of using machine learning models to predict future social behavior with and without causal knowledge explicitly embedded. In this paper, we first present our causal modeling approach to discover the causal structure of four virtual worlds produced by the simulation teams—Urban Life, Financial Governance, Disaster and Geopolitical Conflict. Our approach adapts the state-of-the-art causal discovery (including ensemble models), machine learning, data analytics, and visualization techniques to allow a human-machine team to reverse-engineer the true causal relations from sampled and fully observed data. We next present our reproducibility analysis of two research methods team’s performance using a range of causal discovery models applied to both sampled and fully observed data, and analyze their effectiveness and limitations. We further investigate the generalizability and robustness to sampling of the SOTA causal discovery approaches on additional simulated datasets with known ground truth. Our results reveal the limitations of existing causal modeling approaches when applied to large-scale, noisy, high-dimensional data with unobserved variables and unknown relationships between them. We show that the SOTA causal models explored in our experiments are not designed to take advantage from vasts amounts of data and have difficulty recovering ground truth when latent confounders are present; they do not generalize well across simulation scenarios and are not robust to sampling; they are vulnerable to data and modeling assumptions, and therefore, the results are hard to reproduce. Finally, when we outline lessons learned and provide recommendations to improve models for causal discovery and prediction of human social behavior from observational data, we highlight the importance of learning data to knowledge representations or transformations to improve causal discovery and describe the benefit of causal feature selection for predictive and prescriptive modeling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Encapsulation in a Bacterial Microcompartment Shell Improves Thermal Stability of a Glycolytic Enzyme

Selective encapsulation of target enzymes is an increasingly well-studied field, with a host of potential applications for biotechnology. Natively, many bacteria utilize bacterial microcompartments (BMCs) for enzyme encapsulation to enhance catalysis. BMCs are protein shells that enable selective localization of targeted metabolic enzymes and may improve catalytic rates by colocalizing pathway enzymes and/or serve to sequester toxic or volatile intermediates. The microcompartment shell of Haliangium ochraceum (HO) is a notable BMC chassis because of its modularity and versatility; it is easily expressed and assembled outside its native host and can accept a wide array of cargo. Recently, it was demonstrated that assembly of HO BMC shells can be easily achieved in vitro. Following up on our previous work on in vivo assembly of HO-BMCs with triose phosphate isomerase (TPI) as a model enzyme cargo, here we have demonstrated the advantages of in vitro assembly (IVA) for targeted enzyme encapsulation. We achieved variable loading of BMC shells with targeted amounts of TPI and demonstrated enhanced thermal stability of encapsulated TPI versus free TPI up to 62 °C.

assays↗

DOE OSTI2025

Exploring single-cell biosynthetic noise and dynamics for enhanced betaxanthin production in Escherichia coli

Cell-to-cell variability often limits the efficiency of microbial bioproduction, yet how individual cells fluctuate over time and how these fluctuations shape population-level output remain unclear. To address this issue, we tracked a heterologous betaxanthin pathway in Escherichia coli using microfluidics-assisted time-lapse microscopy, allowing simultaneous measurement of fluctuations in betaxanthin, its biosynthetic enzyme DOD and growth across generations. Here we show that over 50% of high betaxanthin producers become medium or low producers after two divisions. Betaxanthin variation primarily originates from DOD noise, with a smaller contribution from growth rate fluctuations. We further develop a stochastic model to explore various control circuits and find that pathway enzyme or metabolite-based growth selection strategies are most effective in enhancing production. We experimentally validate the model by coupling enzyme expression to nutrient availability, which enriches high producers and boosts titer by 4.4-fold. Our results highlight key sources of metabolic heterogeneity and provide a framework for designing robust microbial processes.

Bacterial systems biology↗

DOE OSTI2019

Evaluation of new approaches to the access of official monitoring results for live bivalves molluscs

Abstract Background Live bivalve molluscs, echinoderms, tunicates and marine gastropod are referred in EU food laws, and require member states to implement official controls in classified production areas, with the monitoring and classification of those areas. If, due to contaminant tests results, a production area is closed, any product from there is prohibited to be commercialized. Mobile applications optical character recognition (OCR) functionalities could ease the access to contaminant levels and production area classifications. This study verifies what information is available in live bivalves’ labels, describes an OCR algorithm for those labels and evaluates it. Methods 86 labels were selected from four sale points in Lisbon, and photographed using smartphones. Each label was evaluated by a human to determine what data was available (either required or not). An OCR algorithm was developed and applied on the collected labels and validated against the data extracted by the human analysis. Results The analysis shows that all the labels included the required information, and 63% of the labels included the identifier for the production zone. The label-reading algorithm performs with an accuracy of 79.85% for the individual values. Conclusion High accuracy of the developed label-reading algorithm shows potential for providing instant automatic access to the date and production area, but is affected by the variability on the label structure. Although not required by food laws, the majority of the sampled labels included complementary information (classified production area) that will allow access to more precise information about the existing biotoxin tests and analysis results.

da Costa, Manuel Beja↗

DOE OSTI2018

Sensitivity of Turbine-Height Wind Speeds to Parameters in the Planetary Boundary-Layer Parametrization Used in the Weather Research and Forecasting Model: Extension to Wintertime Conditions

This work extends the model sensitivity analysis of Yang et al. (Boundary-Layer Meteorology, 2017, Vol. 162, 117-142) to include results for February 2011, in addition to May of the same year. We investigated the sensitivity of simulated hub-height wind speeds to the selection of 26 parameters applied in the Mellor-Yamada-Nakanishi-Niino (MYNN) planetary boundary-layer parameterization in the Weather Research and Forecasting (WRF) model, including parameters used to represent the dissipation of turbulence kinetic energy (TKE), Prandtl number (Pr), and turbulence length scales. Differences in the sensitivity of the ensemble of simulated wind speed to the various parameters can largely be explained by changes in the static stability. The largest monthly differences were found during the day, while the sensitivity to many of the parameters during the night were similar regardless of the month. This finding is consistent with an increased frequency of day-time stable conditions in February compared to May. The spatial variability of the sensitivity to TKE dissipation, Pr, and k can also be attributed to variability in the static stability across the domain at any point in time.

parametrization schemes, Parametric sensitivity, P↗

DOE OSTI2020

Elevation distributed micro-climatology data in a coastal glaciated watershed

Air temperature, ground temperature and relative humidity data were collected in a longitudinal transect of the Nooksack watershed at varying elevations from 500 to 1800 m above sea level. Data were collected by anchoring sensors from trees above winter snow levels and shaded from direct solar radiation. Paired sensors were also buried 3 cm under ground near each air temperature sensor to determine snow absence or presence. Select sites included relative humidity sensors to indicate whether precipitation was occurring. Data were collected every 3-4 h from December 2015 to Sept 2018 (with ongoing collection). Code for analysis of daily mean, minimum, maximum, and temperature change with elevation (lapse rates) are available on Github (https://doi.org/10.5281/zenodo.3239539). The sensor download and intermediate data products are available on HydroShare at (http://www.hydroshare.org/resource/222e832d3df24dea9bae9bbeb6f4219d) with publicly accessible visualization available from the Nooksack Observatory at data.cuahsi.org. Hydrologic models are generally structured with a single annual average lapse rate parameter which assumes a linear temperature gradient with elevation. The daily data (2016-2018) is used as part of ongoing studies on the non-linear dynamics and temporal variability of temperature with elevation to improve assessments of watershed function and salmon habitat.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Gas Hydrate Systems Along the Eastern Continental Margin of India

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Collett, Tim↗

DOE OSTI2016

Genotype × environment interaction analysis of North American shrub willow yield trials confirms superior performance of triploid hybrids

Abstract Development of dedicated bioenergy crop production systems will require accurate yield estimates, which will be important for determining many of the associated environmental and economic impacts of their production. Shrub willow ( Salix spp) is being promoted in areas of the USA and Canada due to its adaption to cool climates and wide genetic diversity available for breeding improvement. Willow breeding in North America is in an early stage, and selection of elite genotypes for commercialization will require testing across broad geographic regions to gain an understanding of how shrub willow interacts with the environment. We analyzed a dataset of first‐rotation shrub willow yields of 16 genotypes across 10 trial environments in the USA and Canada for genotype‐by‐environment interactions using the additive main effects and multiplicative interactions ( AMMI ) model. Mean genotype yields ranged from 5.22 to 8.58 oven‐dry Mg ha −1 yr −1 . Analysis of the main effect of genotype showed that one round of breeding improved yields by as much as 20% over check cultivars and that triploid hybrids, most notably Salix viminalis × S. miyabeana , exhibited superior yields. We also found important variability in genotypic response to environments, which suggests specific adaptability could be exploited among 16 genotypes for yield gains. Strong positive correlations were found between environment main effects and AMMI parameters and growing environment temperatures. These findings demonstrate yield improvements are possible in one generation and will be important for developing cultivar recommendations and for future breeding efforts.

09 BIOMASS FUELS↗

DOE OSTI2016

WE-FG-202-12: Investigation of Longitudinal Salivary Gland DCE-MRI Changes

Purpose: To determine the correlation between dose and changes through treatment in dynamic contrast enhanced (DCE) MRI voxel parameters (Ktrans, kep, Ve, and Vp) within salivary glands of head and neck oropharyngeal squamous cell carcinoma (HNSCC) patients. Methods: 17 HNSCC patients treated with definitive radiation therapy completed DCE-MRI scans on a 3T scanner at pre-treatment, mid-treatment, and post-treatment time points. Mid-treatment and post-treatment DCE images were deformably registered to pre-treatment DCE images (Velocity software package). Pharmacokinetic analysis of the DCE images used a modified Tofts model to produce parameter maps with an arterial input function selected from each patient’s perivertebral space on the image (NordicICE software package). In-house software was developed for voxel-by-voxel longitudinal analysis of the salivary glands within the registered images. The planning CT was rigidly registered to the pre-treatment DCE image to obtain dose values in each voxel. Voxels within the lower and upper dose quartiles for each gland were averaged for each patient, then an average of the patients’ means for the two quartiles were compared. Dose-relationships were also assessed by Spearman correlations between dose and voxel parameter changes for each patient’s gland. Results: Changes in parameters’ means between time points were observed, but inter-patient variability was high. Ve of the parotid was the only parameter that had a consistently significant longitudinal difference between dose quartiles. The highest Spearman correlation was Vp of the sublingual gland for the change in the pre-treatment to mid-treatment values with only a ρ=0.29. Conclusion: In this preliminary study, there was large inter-patient variability in the changes of DCE voxel parameters with no clear relationship with dose. Additional patients may reduce the uncertainties and allow for the determination of the existence of parameter and dose relationships.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

WE-H-BRA-06: Experimental Investigation of RBE for Lung Cancer Cell Lines as a Function of Dose and LET in Proton, Helium and Carbon Beams

Purpose: Investigate and quantify the effect of dose and LET on the RBE of protons, helium and carbon ions. Methods: High throughput, high accuracy experimental setups were custom designed to investigate the Relative Biological Effectiveness (RBE) dependence on the dose and Linear Energy Transfer (LET) values for proton, helium and carbon ion beams. The experiment was conducted at the HIT facility in collaboration with the DKFZ in Heidelberg/Germany. Clonogenic assays of two human lung cancer cell lines, H460 and H1437, were investigated in this study. γH2AX foci staining on the H460 cell line was also undertaken to facilitate the study of differential DNA double-strand break induction and repair between low-design available at the HIT facility. Specific points along the Bragg curve corresponding to well-defined doses and LET values were chosen by appropriate selection of the pre-absorber thicknesses. With a setup design for horizontal beam lines we were able to minimize ion scattering in the cell plate, resulting in narrower energy spectra and hence LET distributions in the Bragg peak and in the distal falloff regions, compared to the earlier experiments. Results: Approximately 16,000 samples of cancer cells were irradiated during 23 hours of beam time. The preliminary results of the survival curves for both cell lines show a distinct dependence on LET for a given dose with decreased survival fractions at increasing LET values, encountered at the Bragg peak and in the distal falloff. Conclusion: Our preliminary findings are indicative of the importance of novel variable-RBE models for proton therapy and provide insight into the RBE of heavy ions for possible future heavy ion therapy facilities in the US. Funding support: SINF 2015/16.

60 APPLIED LIFE SCIENCES↗