Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generator parameter estimation”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Laser Wavelength and Argon Ambient Gas Pressure Effects on Metallic Plasma Emission

The effects of laser wavelength and argon ambient gas pressure on ionic and atomic emission from aluminum and copper metal plasmas has been studied. The plasmas were generated using metal targets excited with a Q-switched Nd:YAG laser operating at the first fundamental (1064 nm), second (532 nm), third (355 nm), and fourth harmonics (266 nm) with repetition rate up to 10 Hz. The spectroscopic analysis of the plasma plume emission revealed that the intensities and signal-to-noise (S/N) ratio of various neutral and ion spectral lines, small components, and trace impurities increase sharply when the fourth laser harmonic is used at a wavelength of 266 nm. Based on LTE assumption, the electron excitation temperature (T{sub e}) and electron number density (N{sub e}) of the produced plasma have been estimated from the Boltzmann plot using neutral spectral lines of Al I at 256.8, 257.5, 309.3, and 396.15 nm and Cu I at 427.58, 510.59, 570.07, and 578.58 nm and from Stark broadening width profile of Al II at 281.6 nm and Cu I at 324.75 nm, respectively. The variations of T{sub e}, N{sub e}, and emission intensities of metals as a function of laser irradiance of different laser wavelengths at 1064, 532, 355, and 266 nm were studied. The effect of the argon ambient gas pressure on the plasma parameters was also studied. The diagnostics of the produced metal plasma was undertaken using laser induced breakdown spectroscopy (LIBS).

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

Assessment of Lung Tumour Motion and Volume Size Dependencies Using Various Evaluation Measures

The purpose of the study was to assess internal target volume changes through the breathing cycle and associated tumour motion for lung patients and to establish possible correlations between different parameters. Respiration-induced volume changes with breathing cycle and the associated tumour motion were analyzed for 11 patients. Selected phases were the maximum and average intensity projections and the 10 phases of equal duration and separation obtained through the respiratory cycle. Tumour centre of mass (COM) motion planes were generated using least square fitting, and their angles and orientations were then compared between the cases studied. Trajectories that are composed by the points of COM location in different phases were identified, and their interrelation was assessed through different similarity measures. The results were used to determine if there is any correlation between parameters chosen and if the margins conventionally used for the planning target volume creation successfully encompassed lung tumour motion and volume change. The results show that the extent of tumour motion was related to its volume and location. The tumour displacement was predominantly left and inferior. Planar fitting to COM motion data through respiratory phases demonstrated some correlation in best fit motion plane positions between different data sets. In the plane fit comparison, for each patient, the lower root mean square error values showed that a good planar fit can be achieved for the COM motion path. The evaluation of the inhale and exhale trajectories may allow, for certain tumour locations and size, contouring on only inhalation or exhalation phases, knowing that tumour motion will be adequately covered on the other phases. Taking all the data into account and knowing the tumour size and location, a good estimate can be made of the motion plane position in the three-dimensional space and the required dosimetric margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Allometric and Mobile Terrestrial LiDAR Modeling of Aboveground Woody Biomass of Populus in Coppice Production

Poplars ( Populus spp.) and their hybrids are increasingly being grown in coppice production to generate bioenergy feedstocks at frequent intervals. Allometric equations are re-quired to predict aboveground biomass (AGB) of coppiced individuals with minimal field measurements. Likewise, remote sensing tools like LiDAR (light detection and ranging) can be used if models are available to predict AGB from point cloud data. Therefore, this study sought to develop equations to predict dry woody AGB from field measurements and LiDAR data from coppiced poplar field trials containing eastern cottonwood ( P. del-toides ) and hybrid poplar taxa. We found that taxa-specific allometric models containing the summed basal area of the three largest stems in the coppice provided the best predictive model, with stem height and stem count failing to provide additional explanatory power. The best predictive LiDAR-based model was independent of taxa but had slightly lower adjusted R 2 and higher RMSE than the allometric model. It contained four parameters including crown volume, leaf area index, variance of height returns, and the top point density (i.e., density metric 9 or the proportion of points in the highest point interval when the point cloud is evenly divided into ten vertical intervals). In total, these models can be used to quickly and efficiently estimate dry woody AGB of Populus coppice systems for bioenergy feedstock production.

AGB↗
DOE OSTI2017

Optimizing observations of drizzle onset with millimeter-wavelength radars

Cloud Doppler radars are increasingly used to study cloud and precipitation microphysical processes. Typical bulk cloud properties such as liquid or ice content are usually derived using the first three standard moments of the radar Doppler spectrum. Recent studies demonstrated the value of higher moments for the reduction of retrieval uncertainties and for providing additional insights into microphysical processes. Large effort has been undertaken, e.g., within the Atmospheric Radiation Measurement (ARM) program to ensure high quality of radar Doppler spectra. However, a systematic approach concerning the accuracy of higher moment estimates and sensitivity to basic radar system settings, such as spectral resolution, integration time and beam width, are still missing. Here In this study, we present an approach on how to optimize radar settings for radar Doppler spectra moments in the specific context of drizzle detection. The process of drizzle development has shown to be particularly sensitive to higher radar moments such as skewness. We collected radar raw data (I/Q time series) from consecutive zenith-pointing observations for two liquid cloud cases observed at the cloud observatory JOYCE in Germany. The I/Q data allowed us to process Doppler spectra and derive their moments using different spectral resolutions and integration times during identical time intervals. This enabled us to study the sensitivity of the spatiotemporal structure of the derived moments to the different radar settings. The observed signatures were further investigated using a radar Doppler forward model which allowed us to compare observed and simulated sensitivities and also to study the impact of additional hardware-dependent parameters such as antenna beam width. For the observed cloud with drizzle onset we found that longer integration times mainly modify spectral width (S w ) and skewness (S k ), leaving other moments mostly unaffected. An integration time of 2 s seems to be an optimal compromise: both observations and simulations revealed that a 10 s integration time – as it is widely used for European cloud radars – leads to a significant turbulence-induced increase of S w and reduction of S k compared to 2 s integration time. This can lead to significantly different microphysical interpretations with respect to drizzle water content and effective radius. A change from 2 s to even shorter integration times (0. 4 s) has much smaller effects on S w and S k . We also find that spectral resolution has a small impact on the moment estimations, and thus on the microphysical interpretation of the drizzle signal. Even the coarsest spectral resolution studied, 0. 08 ms -1 , seems to be appropriate for calculation moments of drizzling clouds. Moreover, simulations provided additional insight into the microphysical interpretation of the skewness signatures observed: in low (high)-turbulence conditions, only drizzle larger than 20 µm (40 µm) can generate S k values above the S k noise level (in our case 0.4). Higher S k values are also obtained in simulations when smaller beam widths are adopted.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

LCOE Uncertainty Analysis for Hydropower using Monte Carlo Simulations

Levelized Cost of Energy (LCOE) is an important metric to evaluate the cost and performance of electricity production generation alternatives, and combined with other measures, can be used to assess the economics of future hydropower development. Multiple assumptions on input parameters are required to calculate the LCOE, which each contain some level of uncertainty, in turn affecting the accuracy of LCOE results. This paper explores these uncertainties, their sources, and ultimately the level of variability they introduce at the screening level of project evaluation for non-powered dams (NPDs) across the U.S.Owing to site-specific differences in site design, the LCOE for hydropower varies significantly from project to project unlike technologies with more standardized configurations such as wind and gas. Therefore, to assess the impact of LCOE input uncertainty on the economics of U.S. hydropower resources, these uncertainties must be modeled across the population of potential opportunities. To demonstrate the impact of uncertainty, resource data from a recent nationwide non-powered dam (NPD) resource assessment (Hadjerioua et al., 2012) and screening-level predictive cost equations (O’Connor et al., 2015) are used to quantify and evaluate uncertainties in project capital and operations & maintenance costs, and generation potential at broad scale. LCOE dependence on financial assumptions is also evaluated on a sensitivity basis to explore ownership/investment implications on project economics for the U.S. hydropower fleet.The results indicate that the LCOE for U.S. NPDs varies substantially. The LCOE estimates for the potential NPD projects of capacity greater than 1 MW range from 40 to 182 $/MWh, with average of 106 $/MWh. 4,000 MW could be developed through projects with individual LCOE values below 100 $/MWh. The results also indicate that typically 90 % of LCOE uncertainty can be attributed to uncertainties in capital costs and energy production; however, for small projects (below 10 MW) O&M uncertainty plays an important role. The potential NPDs have costs with uncertainty bands extending well above 100 $/MWh in some cases. A truly representative evaluation should complement LCOE with an assessment of potential revenues from the provision of ancillary services. Moreover, hydroelectric production is often part of multipurpose water resource projects where other purposes (e.g., irrigation, flood control, navigation, recreation, etc.) also provide value. Neither of these considerations are accounted for LCOE.

Chalise, Dol Raj↗
DOE OSTI2015

Improving Limit Surface Search Algorithms in RAVEN Using Acceleration Schemes: Level II Milestone

The RAVEN code is becoming a comprehensive tool to perform Probabilistic Risk Assessment (PRA); Uncertainty Quantification (UQ) and Propagation; and Verification and Validation (V&V). The RAVEN code is being developed to support the Risk-Informed Safety Margin Characterization (RISMC) pathway by developing an advanced set of methodologies and algorithms for use in advanced risk analysis. The RISMC approach uses system simulator codes applied to stochastic analysis tools. The fundamental idea behind this coupling approach to perturb (by employing sampling strategies) timing and sequencing of events, internal parameters of the system codes (i.e., uncertain parameters of the physics model) and initial conditions to estimate values ranges and associated probabilities of figures of merit of interest for engineering and safety (e.g. core damage probability, etc.). This approach applied to complex systems such as nuclear power plants requires performing a series of computationally expensive simulation runs. The large computational burden is caused by the large set of (uncertain) parameters characterizing those systems. Consequently, exploring the uncertain/parametric domain, with a good level of confidence, is generally not affordable, considering the limited computational resources that are currently available. In addition, the recent tendency to develop newer tools, characterized by higher accuracy and larger computational resources (if compared with the presently used legacy codes, that have been developed decades ago), has made this issue even more compelling. In order to overcome to these limitations, the strategy for the exploration of the uncertain/parametric space needs to use at best the computational resources focusing the computational effort in those regions of the uncertain/parametric space that are “interesting” (e.g., risk-significant regions of the input space) with respect the targeted Figures Of Merit (FOM): for example, the failure of the system, subject of the analysis. These methodologies are named, in the RAVEN environment, adaptive sampling strategies. These methodologies infer system responses from surrogate models constructed from already existing samples (produced using high fidelity simulations) and suggest the most relevant location (coordinate in the input space) of the next sampling point to be explored in the uncertain/parametric domain. When using those methodologies, it is possible to understand features of the system response with a small number of carefully selected samples. This report focuses on the development and improvement of the limit surface search. The limit surface is an important concept in system reliability analysis. Without going into the details, which will be covered later in the report, the limit surface could be briefly described as an hyper-surface in the system uncertainty/parametric space separating the regions leading to a prescribed system outcome. For example, if the uncertainty/parametric space is the one generated by the reactor power level and the duration of the batteries, the system is a nuclear power plant and the system outcome discriminating variable is the clad failure in a station blackout scenario, then the limit surface separates the combinations of reactor power level and battery duration that lead to clad failure from the ones that do not.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Communication Network Layer State Estimation Measurement Model for a Cyber-Secure Smart Grid

Network communication has been proven to be a very important tool and a key factor in the recent development and progress of the power grid operation. It is also considered as the foundation for the smart grid because information and communication are integrated into electricity distribution to achieve reliable and accurate knowledge of the power grid. In previous years, absorbing energy from substations and delivering it to customers was the only type of interaction we knew between utility companies and customers. Presently, the growing connections of small distributed generation units caused by the cost reduction of most of the technologies used in generation and storage of electrical energy, along with the potential benefits of renewable energy have pushed many researchers to look into the improvement of information and communication technologies (ICT) in order to ensure a bidirectional flow of power and data. Moreover, the evolution of information and communication technologies and its applications to smart grid have converted the smart grid into a cyber-physical system where vulnerabilities and additional security challenges such as cyber-threats and cyber-attacks have emerged. Previously, we have demonstrated that using machine learning-based processing on data gathered from communication networks and the power grid was a promising solution for detecting cyber threats by implementing a co-simulation of cyber-security for cross-layer strategy. Since the majority of the challenges observed can only be solved in the network communication layer, we present in this work a physics-based state estimation model of the communication network system towards enhanced cyber-physical security of the smart grid. Information integration with the previously developed machine learning model is developed, providing a enhanced cyber-physical security application for the smart grid. Easy-to-implement model, without hard-to-derive parameters, highlight potential aspects of the model for real-life applications.

Mathieu, Reynold↗
DOE OSTI2016

Simulation Model of Microsphere Distribution for Selective Internal Radiation Therapy Agrees With Observations

Purpose: To perform a detailed analysis of microsphere distribution in biopsy material from a patient treated with {sup 90}Y-labeled resin spheres and characterize microsphere distribution in the hepatic artery tree, and to construct a novel dichotomous bifurcation model for microsphere deposits and evaluate its accuracy in simulating the observed microsphere deposits. Methods and Materials: Our virtual model consisted of arteries that successively branched into 2 new generations of arteries at 20 nodes. The artery diameter exponentially decreased from the lowest generation to the highest generation. Three variable parameters were optimized to obtain concordance between simulations and measure microsphere distributions: an artery coefficient of variation (ACV) for the diameter of all artery generations and the microsphere flow distribution at the nodes; a hepatic tree distribution volume (HDV) for the artery tree; and an artery diameter reduction (ADR) parameter. The model was tested against previously measured activity concentrations in 84 biopsies from the liver of 1 patient. In 16 of 84 biopsies, the microsphere distribution regarding cluster size and localization in the artery tree was determined via light microscopy of 30-μm sections (mean concentration, 14 microspheres/mg; distributions divided into 3 groups with mean microsphere concentrations of 4.6, 14, and 28 microspheres/mg). Results: Single spheres and small clusters were observed in terminal arterioles, whereas large clusters, up to 450 microspheres, were observed in larger arterioles. For 14 microspheres/mg, the optimized parameter values were ACV=0.35, HDV = 50 cm{sup 3}, and ADR=6 μm. For 4.6 microspheres/mg, ACV and ADR decreased to 0.26 and 0 μm, respectively, whereas HDV increased to 130 cm{sup 3}. The opposite trend was observed for 28 microspheres/mg: ACV = 0.49, HDV = 20 cm{sup 3}, and ADR = 8 μm. Conclusion: Simulations and measurements reveal that microsphere clusters are larger and more common in volumes with high microsphere concentrations and indicate that the spatial distribution of the artery tree must be considered in estimates of microsphere distributions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.

58 GEOSCIENCES↗
DOE OSTI2021

Using radar observations to evaluate 3-D radar echo structure simulated by the Energy Exascale Earth System Model (E3SM) version 1

Abstract. The Energy Exascale Earth System Model (E3SM) developed by the Department of Energy has a goal of addressing challenges in understanding the global water cycle. Success depends on correct simulation of cloud and precipitation elements. However, lack of appropriate evaluation metrics has hindered the accurate representation of these elements in general circulation models. We derive metrics from the three-dimensional data of the ground-based Next-Generation Radar (NEXRAD) network over the US to evaluate both horizontal and vertical structures of precipitation elements. We coarsened the resolution of the radar observations to be consistent with the model resolution and improved the coupling of the Cloud Feedback Model Intercomparison Project Observation Simulator Package (COSP) and E3SM Atmospheric Model Version 1 (EAMv1) to obtain the best possible model output for comparison with the observations. Three warm seasons (2014–2016) of EAMv1 simulations of 3-D radar reflectivity features at an hourly scale are evaluated. A general agreement in domain-mean radar reflectivity intensity is found between EAMv1 and NEXRAD below 4 km altitude; however, the model underestimates reflectivity over the central US, which suggests that the model does not capture the mesoscale convective systems that produce much of the precipitation in that region. The shape of the model-estimated histogram of subgrid-scale reflectivity is improved by correcting the microphysical assumptions in COSP. Different from previous studies that evaluated modeled cloud top height, we find the model severely underestimates radar reflectivity at upper levels – the simulated echo top height is about 5 km lower than in observations – and this result is not changed by tuning any single physics parameter. For more accurate model evaluation, a higher-order consistency between the COSP and the host model is warranted in future studies.

58 GEOSCIENCES↗
NASA NTRS2016

A Method for Improving Hotspot Directional Signatures in BRDF Models Used for MODIS

The semi-empirical, kernel-driven, linear RossThick-LiSparseReciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is used to generate the routine MODIS BRDFAlbedo product due to its global applicability and the underlying physics. A challenge of this model in regard to surface reflectance anisotropy effects comes from its underestimation of the directional reflectance signatures near the Sun illumination direction; also known as the hotspot effect. In this study, a method has been developed for improving the ability of the RTLSR model to simulate the magnitude and width of the hotspot effect. The method corrects the volumetric scattering component of the RTLSR model using an exponential approximation of a physical hotspot kernel, which recreates the hotspot magnitude and width using two free parameters (C(sub 1) and C(sub 2), respectively). The approach allows one to reconstruct, with reasonable accuracy, the hotspot effect by adjusting or using the prior values of these two hotspot variables. Our results demonstrate that: (1) significant improvements in capturing hotspot effect can be made to this method by using the inverted hotspot parameters; (2) the reciprocal nature allow this method to be more adaptive for simulating the hotspot height and width with high accuracy, especially in cases where hotspot signatures are available; and (3) while the new approach is consistent with the heritage RTLSR model inversion used to estimate intrinsic narrowband and broadband albedos, it presents some differences for vegetation clumping index (CI) retrievals. With the hotspot-related model parameters determined a priori, this method offers improved performance for various ecological remote sensing applications; including the estimation of canopy structure parameters.

airborne measurements↗
DOE OSTI2022

Development and Validation of Algorithms That Analyze Communicating Thermostat Data to Identify Enclosure Retrofit Opportunities

Annual energy savings of up to $\$ 4$ to $\$ 5$ billion could be achieved nationwide through basic insulation and heating system retrofits of existing homes. However, current utility energy efficiency programs are costly and challenging to scale. Customer acquisition occurs primarily through energy bill mailers, mass media, and online advertising that lack specificity about home-specific retrofit opportunities, expected energy savings, and cost-effectiveness. Specific retrofit opportunities are identified via on-site home energy assessments (HEAs) that are inconvenient to homeowners, expensive, and of variable accuracy. We developed computational algorithms that automatically analyze communicating thermostat (CT) heating data that could be used to increase the customer uptake of insulation and air sealing energy conservation measures (ECMs) by identifying homes with the most significant retrofit opportunities, estimating post-retrofit energy savings, and formulating home-specific outreach. The algorithms are based on an extended second-order grey-box model that characterizes a building’s thermal response using lumped elements, coupled with an empirical model of infiltration that accounts for both wind and stack effects. The basic parameters of the model correspond to actual physical parameters of the home, i.e., the home’s overall R-value of and the building envelope ACH50. Unlike the conventional approach, which estimates model parameters based on the best fit to the observed time-dependent room temperature, our approach derives correlations between the daily heating system runtime and temperature difference (indoor-outdoor) that are more robust to data quality issues in real-world applications. We also used HEA data for algorithm development and validation. With the help of our utility partners, Eversource and National Grid, we obtained data sets for hundreds of Massachusetts homes. For each home, these data sets included three sets of information anonymized by the utility: (1) CT data (HVAC runtime, room temperature, and, for some vendors, outdoor temperature and wind speed) collected by the CT vendor (one of three) over a heating season, (2) HEA report performed by the HEA vendor (same vendor for all homes), (3) Monthly utility gas bills coincident with the CT data (3 to 24 per home, depending on availability). For some homes, we also obtained blower-door test results. Initially, we applied the algorithms developed to homes with a single CT and then extended them to homes with two CTs by using an equivalent home approach. Finally, we developed algorithms for prediction of energy savings and a methodology of comparing our predictions with those generated by HEAs. The main technical results indicate that we can reliably identify homes with insulation and/or air sealing retrofit opportunities and provide accurate savings predictions. Our hypothesis is that the algorithms could be applied to utility energy efficiency programs to identify homes that could realize significant energy savings from insulation and/or air sealing retrofits. This information could then be used to reach out to those homes with highly customized outreach, thereby delivering increased program energy savings and cost-effectiveness. This would: Significantly increase the uptake rate of on-site HEAs, and Significantly increase the fraction of HEAs resulting in ECM implementation. To test these hypotheses, we designed and conducted a randomized controlled trial (RCT). The RCT results suggest that personal messaging leads to a two- to five-fold increase in the HEA uptake rate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Thermalized axion inflation

We analyze the dynamics of inflationary models with a coupling of the inflaton φ to gauge fields of the form φ F F-tilde / f , as in the case of axions. It is known that this leads to an instability, with exponential amplification of gauge fields, controlled by the parameter ξ= φ-dot /(2 fH ), which can strongly affect the generation of cosmological perturbations and even the background. We show that scattering rates involving gauge fields can become larger than the expansion rate H , due to the very large occupation numbers, and create a thermal bath of particles of temperature T during inflation. In the thermal regime, energy is transferred to smaller scales, radically modifying the predictions of this scenario. We thus argue that previous constraints on ξ are alleviated. If the gauge fields have Standard Model interactions, which naturally provides reheating, they thermalize already at ξ∼>2.9, before perturbativity constraints and also before backreaction takes place. In absence of SM interactions (i.e. for a dark photon), we find that gauge fields and inflaton perturbations thermalize if ξ∼>3.4; however, observations require ξ∼>6, which is above the perturbativity and backreaction bounds and so a dedicated study is required. After thermalization, though, the system should evolve non-trivially due to the competition between the instability and the gauge field thermal mass. If the thermal mass and the instabilities equilibrate, we expect an equilibrium temperature of T {sub eq} ≅ ξ H / g-bar where g-bar is the effective gauge coupling. Finally, we estimate the spectrum of perturbations if φ is thermal and find that the tensor to scalar ratio is suppressed by H /(2 T ), if tensors do not thermalize.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Integrated Heat Exchanger-Phase Change Material Thermal Energy Storage System

The purpose of this study is to experimentally investigate the thermal performance of an innovative thermal energy storage (TES) system that combines the advantages of the phase-change material (PCM)/graphite foam latent heat TES medium developed at Argonne National Laboratory (Argonne) and the internally supported plate-fin (ISPE) cell architecture heat transfer fluid (HTF) flow channels developed at Brayton Energy (Brayton). Several essential tasks were accomplished: (1) Thermal property characterization. Thermal properties of the graphite foam were characterized, providing necessary data for experimental result analysis and numerical simulation. (2) Design and optimization of lab-scale test module. Based on Brayton’s full-scale heat exchanger (HX)-TES system, the experimental test module was designed, optimized, and fabricated. (3) Thermal performance testing and data analysis. Five cycle tests were successfully conducted—including one with approximately 3.5 psig of pressure applied to the diaphragms—to investigate the thermal performance of the experimental test module for charging and discharging. Temperature profiles were generated for each charging test and discharging test as a function of time. The temperature profiles clearly show three TES stages: sensible heat (temperature increase), latent heat (melting), and sensible heat (temperature increase) for the charging process. Similarly, the temperature profiles clearly show three thermal energy release stages: sensible heat (temperature decrease), latent heat (solidification), and sensible heat (temperature decrease). Melting and solidification of the PCM generally occurred in relatively narrow temperature ranges, indicated by the flattened temperature regions in the temperature profiles. These phase changes ranged approximately 3°C for melting and 3.5°C for solidification. The charging and discharging temperature profiles were similar for similar experimental parameter tests whether or not pressure was applied to the diaphragm to eliminate the gap between the HX surface and the TES subsystem. This indicates that the effect of a small gap between the HX surface and the TES subsystem is insignificant for charging and discharging. (4) Comparison of experimental data and simulation results. We compared the experimental data to the numerical simulation results. Numerical simulations were conducted by using the ANSYS FLUENT 2019 R3 commercial computational fluid dynamics software. The predicted phase-change times agreed reasonably well with those from the experimental data. In most cases, the estimated time differences between the relative phase changes were within 16%. The predicted start and end times for the charging process agreed well with those from the experimental data. However, the simulation results showed earlier start and end times than the experimental data for the discharging process. Overall, the experimental data and its comparison with the simulation predictions verified the technical viability of the integrated ISPF HX-PCM/graphite foam latent-heat TES system.

25 ENERGY STORAGE↗
DOE OSTI2018

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development

This DOE GTO funded project was a collaborative effort between Array Information Technology and the University of California at Berkeley. The project goal is to improve technology to assess in-situ stress magnitude and orientation, kinematic fracture parameters, rupture size, rupture orientation as well as temporal changes and volumetric distribution of the injected fluid during Enhanced Geothermal System (EGS) resource development. The project leveraged high-frequency seismic data recorded by the LBNL 34-station permanent geophone network and seismic broadband data recorded by a temporary 33-station seismometer network that operated in The Geysers during the injection phase of the DOE GTO funded EGS demonstration project at the well Prati-32. The operation of the broadband network with high station density and good azimuthal coverage to monitor the injection phase produced an unprecedented dataset that is typically unavailable for EGS operations. The two datasets offer the advantages of studying the parameters under investigation over a broader frequency band and the higher sensitivity of the broadband seismometers resulted in a richer dataset with a lower magnitude of completeness. Moment tensor (MT) inversion have been employed to develop a catalog of earthquake source solutions including moment tensor and finite source rupture area. This resulted in a 168-event waveform-based seismic moment tensor catalog for events 0.7≤Mw≤3.9 based on a semi-automated moment tensor analysis technique. The MT solutions were used for in-situ stress estimation during the EGS injection phases. We found an approximate 15-degree counterclockwise rotation of the least compressive stress sigma-3, and a rotation of the maximum compressive stress sigma-1 toward the vertical as the injected volume of water increased. We developed a rupture-area magnitude scaling relationship of The Geysers earthquakes obtained from finite-source inversion for fault slip that reveals a very high correlation to published results. Based on the rupture area-magnitude relationship we mapped Mw from the LBNL earthquake catalog to estimate the activated fracture area, its location and orientation within the injection volume of Prati-32. The double-difference Wadati (DDW) technique was applied to high-accuracy P- and S-wave differential travel times derived from waveform cross correlation to image fluid saturation in the subsurface based on Vp/Vs-ratio. The recently developed technique has been applied to estimate temporal changes of fluid distribution in fault zones and in volcanic regions but not yet in engineered geo-reservoirs. The improvement of the technology includes high-resolution waveform cross correlation applied to the high-quality datasets with lower magnitude of completeness resulting in a higher number of recorded events. The improvement allowed to estimate both, the spatial and temporal distribution of injected fluids and the state of the injected fluid during reservoir stimulation. Results of the DDW analysis revealed the development of an injection-derived steam plume during the early injection phase, particularly in the northern part of the Prati-32 study area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Influence of regional nighttime atmospheric regimes on canopy turbulence and gradients at a closed and open forest in mountain-valley terrain

Stable stratification of the nocturnal lower boundary layer inhibits convective turbulence, such that turbulent vertical transfer of ecosystem carbon dioxide (CO 2 ), water vapor (H 2 O) and energy is driven by mechanically forced turbulence, either from frictional forces near the ground or top of a plant canopy, or from shear generated aloft. The significance of this last source of turbulence on canopy flow characteristics in a closed and open forest canopy is addressed in this paper. We present micrometeorological observations of the lower boundary layer and canopy air space collected on nearly 200 nights using a combination of atmospheric laser detection and ranging (lidar), eddy covariance (EC), and tower profiling instrumentation. Two AmeriFlux/Fluxnet sites in mountain-valley terrain in the Western U.S. are investigated: Wind River, a tall, dense conifer canopy, and Tonzi Ranch, a short, open oak canopy. On roughly 40% of nights lidar detected down-valley or downslope flows above the canopy at both sites. Nights with intermittent strong bursts of “top-down” forced turbulence were also observed above both canopies. The strongest of these bursts increased sub-canopy turbulence and reduced canopy virtual potential temperature (θv) gradient at Tonzi, but did not appear to change the flow characteristics within the dense Wind River canopy. At Tonzi we observed other times when high turbulence (via friction velocity, u*) was found just above the trees, yet CO2 and θv gradients remained large and suggested flow decoupling. These events were triggered by regional downslope flow. Lastly, a set of turbulence parameters is evaluated for estimating canopy turbulence mixing strength. The relationship between turbulence parameters and canopy θv gradients was found to be complex, although better agreement between the canopy θv gradient and turbulence was found for parameters based on the standard deviation of vertical velocity, or ratios of 3-D turbulence to mean flow, than for u*. These findings add evidence that the relationship between canopy turbulence, static stability, and canopy mixing is far from straightforward even within an open canopy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Accounting for partiality in serial crystallography using ray-tracing principles

Serial crystallography generates partial reflections from still diffraction images. Partialities are estimated with EVAL ray-tracing simulations, thereby improving merged reflection data to a similar quality as conventional rotation data. Serial crystallography generates ‘still’ diffraction data sets that are composed of single diffraction images obtained from a large number of crystals arbitrarily oriented in the X-ray beam. Estimation of the reflection partialities, which accounts for the expected observed fractions of diffraction intensities, has so far been problematic. In this paper, a method is derived for modelling the partialities by making use of the ray-tracing diffraction-integration method EVAL. The method estimates partialities based on crystal mosaicity, beam divergence, wavelength dispersion, crystal size and the interference function, accounting for crystallite size. It is shown that modelling of each reflection by a distribution of interference-function weighted rays yields a ‘still’ Lorentz factor. Still data are compared with a conventional rotation data set collected from a single lysozyme crystal. Overall, the presented still integration method improves the data quality markedly. The R factor of the still data compared with the rotation data decreases from 26% using a Monte Carlo approach to 12% after applying the Lorentz correction, to 5.3% when estimating partialities by EVAL and finally to 4.7% after post-refinement. The merging R{sub int} factor of the still data improves from 105 to 56% but remains high. This suggests that the accuracy of the model parameters could be further improved. However, with a multiplicity of around 40 and an R{sub int} of ∼50% the merged still data approximate the quality of the rotation data. The presented integration method suitably accounts for the partiality of the observed intensities in still diffraction data, which is a critical step to improve data quality in serial crystallography.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗