Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sensitivity study”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Modeling surface water dynamics in the Amazon Basin using MOSART-Inundation v1.0: impacts of geomorphological parameters and river flow representation

Abstract. In the Amazon Basin, floodplain inundation is a key component of surface water dynamics and plays an important role in water, energy and carbon cycles. The Model for Scale Adaptive River Transport (MOSART) was extended with a macroscale inundation scheme for representing floodplain inundation. The extended model, named MOSART-Inundation, was used to simulate surface hydrology of the entire Amazon Basin. Previous hydrologic modeling studies in the Amazon Basin identified and addressed a few challenges in simulating surface hydrology of this basin, including uncertainties of floodplain topography and channel geometry, and the representation of river flow in reaches with mild slopes. This study further addressed four aspects of these challenges. First, the spatial variability of vegetation-caused biases embedded in the HydroSHEDS digital elevation model (DEM) data was explicitly addressed. A vegetation height map of about 1 km resolution and a land cover dataset of about 90 m resolution were used in a DEM correction procedure that resulted in an average elevation reduction of 13.2 m for the entire basin and led to evident changes in the floodplain topography. Second, basin-wide empirical formulae for channel cross-sectional dimensions were refined for various subregions to improve the representation of spatial variability in channel geometry. Third, the channel Manning roughness coefficient was allowed to vary with the channel depth, as the effect of riverbed resistance on river flow generally declines with increasing river size. Lastly, backwater effects were accounted for to better represent river flow in mild-slope reaches. The model was evaluated against in situ streamflow records and remotely sensed Envisat altimetry data and Global Inundation Extent from Multi-Satellites (GIEMS) inundation data. In a sensitivity study, seven simulations were compared to evaluate the impacts of the five modeling aspects addressed in this study. The comparisons showed that representing floodplain inundation could significantly improve the simulated streamflow and river stages. Refining floodplain topography, channel geometry and Manning roughness coefficients, as well as accounting for backwater effects had notable impacts on the simulated surface water dynamics in the Amazon Basin. The understanding obtained in this study could be helpful in improving modeling of surface hydrology in basins with evident inundation, especially at regional to continental scales.

58 GEOSCIENCES↗
DOE OSTI2016

A LASER INTERFERENCE-BASED SURFACE TREATMENT OF AL AND CARBON FIBER POLYMER COMPOSITES FOR ENHANCED BONDING

Due to its increased use in the automotive and aerospace industries, joining of Carbon Fiber-reinforced Polymer matrix Composites (CFPC) to metals demands enhanced surface preparation and control of surface morphology prior to joining. In this study, surfaces of both composite and aluminum were prepared for joining using a new laser based technique, in which the laser interference power profile was created by splitting the beam and guiding those beams to the sample surface by overlapping each other with defined angles to each other. Results were presented for the overlap shear testing of single-lap joints made with Al 5182 and CFPC specimens whose surfaces prepared by (a) surface abrasion and solvent cleaning; and (b) laser-interference structured surfaces by rastering with a 4 mm laser beam at approximately 3.5 W power. CFPC specimens of T700S carbon fiber, Prepreg T70 epoxy, 4 or 5 ply thick, 0/90o plaques were used. Adhesive DP810 was used to bond Al and CFPC. The bondline was 0.25mm and the bond length was consistent among all joints produced. First, the effect of the laser speed on the joint performance was evaluated by laser-interference structure Al and CFPC surfaces with a beam angle of 3o and laser beam speeds of 3, 5, and 10 mm/s. For this sensitivity study, 3 joint specimens were used per each joint type. Based on the results for minimum, maximum, and mean values for the shear lap strength and maximum load for all the 9 joint types, two joint types were selected for further evaluations. Six additional joint specimens were prepared for these two joint types in order to obtain better statistics and the shear test data was presented for the range, mean, and standard deviation. The results for the single-lap shear tests obtained for six joint specimens, indicate that the shear lap strength, maximum load, and displacement at maximum load for those joints made with laser-interference structured surfaces were increased by approximately 14.8%, 16%, and 100%, respectively over those measured for the baseline joints.

Sabau, Adrian S↗
NASA NTRS1993

The three-dimensional structure of cumulus clouds over the ocean. 1: Structural analysis

Thermal channel (channel 6, 10.4-12.5 micrometers) images of five Landsat thematic mapper cumulus scenes over the ocean are examined. These images are thresholded using the standard International Satellite Cloud Climatology Project (ISCCP) thermal threshold algorithm. The individual clouds in the cloud fields are segmented to obtain their structural statistics which include size distribution, orientation angle, horizontal aspect ratio, and perimeter-to-area (PtA) relationship. The cloud size distributions exhibit a double power law with the smaller clouds having a smaller absolute exponent. The cloud orientation angles, horizontal aspect ratios, and PtA exponents are found in good agreement with earlier studies. A technique also is developed to recognize individual cells within a cloud so that statistics of cloud cellular structure can be obtained. Cell structural statistics are computed for each cloud. Unicellular clouds are generally smaller (less than or equal to 1 km) and have smaller PtA exponents, while multicellular clouds are larger (greater than or equal to 1 km) and have larger PtA exponents. Cell structural statistics are similar to those of the smaller clouds. When each cell is approximated as a quadric surface using a linear least squares fit, most cells have the shape of a hyperboloid of one sheet, but about 15% of the cells are best modeled by a hyperboloid of two sheets. Less than 1% of the clouds are ellipsoidal. The number of cells in a cloud increases slightly faster than linearly with increasing cloud size. The mean nearest neighbor distance between cells in a cloud, however, appears to increase linearly with increasing cloud size and to reach a maximum when the cloud effective diameter is about 10 km; then it decreases with increasing cloud size. Sensitivity studies of threshold and lapse rate show that neither has a significant impact upon the results. A goodness-of-fit ratio is used to provide a quantitative measure of the individual cloud results. Significantly improved results are obtained after applying a smoothing operator, suggesting the eliminating subresolution scale variations with higher spatial resolution may yield even better shape analyses.

Kuo, Kwo-Sen↗
DOE OSTI2024

Cost Analysis and Projections for U.S.-Manufactured Automotive Lithium-ion Batteries

This document reports on a study conducted to estimate the cost of U.S-produced automotive battery packs for model years (MY) 2023 to 2035, using Argonne National Laboratory’s BatPaC tool. The costs were estimated by designing batteries for several classes of vehicles for four discrete model years (2023, 2026, 2030, and 2035), where a representative battery technology and material prices are selected based on information available today. Correlations were developed from the four discrete years to enable annual pack cost estimates as a function of pack size (kWh) and model year. A consolidated cost curve was then developed that includes battery size, technology by model year, and the anticipated sales volumes of each class of vehicles over the years. This cost curve estimates the volume-averaged, U.S.-manufactured battery pack cost of PHEVs and BEVs in the United States to be $\$140$/kWh for the model year 2023, which will reduce to $\$86$/kWh in MY2035. Applying tax credits from section 45X of the Inflation Reduction Act can further reduce the average pack cost to as low as $\$56$/kWh in MY2029. The report also includes several sensitivity studies that investigate the effect of pack production volume, material prices, fast charge requirements, and labor rates.

25 ENERGY STORAGE↗
DOE OSTI2018

Comparative environmental and human health evaluations of thermolysis and solvolysis recycling technologies of carbon fiber reinforced polymer waste

Highlights: • Performed life cycle impact assessment (LCIA) of pyrolysis for CFRP waste recycling. • Performed LCIA of solvolysis using supercritical water (SCW) for CFRP recycling. • Performed comparative LCIA of pyrolysis vs SCW-based solvolysis for CFRP recycling. • Closed current knowledge gap about sustainability of SCW-based solvolysis. • Compared LCIA of pyrolyzer’s heating using natural gas in lieu of electricity. - Abstract: This quantitative research aims to compare environmental and human health impacts associated with two recycling technologies of CFRP waste. The ‘baseline’ recycling technology is the conventional thermolysis process via pyrolysis and the ‘alternative’ recycling technology is an emerging chemical treatment via solvolysis using supercritical water (SCW) to digest the thermoset matrix. Two Gate-to-Gate recycling models are developed using GaBi LCA platform. The selected functional unit (FU) is 1 kg CFRP waste and the geographical boundary of this comparative LCIA is defined to be within the U.S. The results of this comparative assessment brought to light new insights about the environmental and human health impacts of CFRP waste recycling via solvolysis using SCW and, therefore, helped close a gap in the current state of knowledge about sustainability of SCW-based solvolysis as compared to pyrolysis. Two research questions are posed to identify whether solvolysis recycling offers more environmental and human health gains relative to the conventional pyrolysis recycling. These research questions lay the basis for formulating two null hypotheses (H{sub 0,1} and H{sub 0,2}) and their associated research hypotheses (H{sub 1,1} and H{sub 1,2}). LCIA results interpretation included ‘base case’ scenarios, ‘sensitivity studies,’ and ‘scenarios analysis.’ The results revealed that: (a) recycling via solvolysis using SCW exhibits no gains in environmental and human health impacts relative to those impacts associated with recycling via pyrolysis and (b) use of natural gas in lieu of electricity for pyrolyzer’s heating reduces the environmental and human health impacts by 37% (lowest) and up to 95.7% (highest). It is recommended that on-going experimental efforts that focus only on identifying the best solvent for solvolysis-based recycling should also consider quantification of the energy intensity as well as environmental and human health impacts of the proposed solvents.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

On the Radiation Effects of Strontium Ions on Satellite Solar Cells in Low Earth Orbits

This study focuses on the radiation effects 19 of Sr + ions—generated from high-altitude nuclear explosions (HANE)—on satellite solar cells in low-Earth orbits (LEO). Along four selected satellite orbits, ion fluences are sampled inside the evolving Sr + distributions for days, determined from our newly developed HANE environment model. These fluences, along with the help of radiation transport codes including the MULASSIS and SRIM models, enable us to quantify the radiation damages by determining the values of total ionizing doses and the equivalent 1 MeV electron fluences for displacement damages. Comparing the dose values to existing experimental data, we conclude that HANE-generated Sr + ions have limited darkening effects to quartz solar cell coverglasses in LEO with apogees of 100s to 1000 km. In addition, with the extremely high equivalent fluences, we also conclude that these Sr+ ions may cause severe or even fatal displacement damage to exposed solar photovoltaic (PV) cells on satellites in LEO. The radiation effects of Sr + ions are much less significant for the orbits with high apogees beyond ten thousand km. We also conducted model parameter sensitivity studies on the charge exchange cross-sections, neutral atmosphere density profiles and explosion local time positions, and the above conclusions stay unchanged. The methodology developed in this study can be extended to other HANE-generated heavy ion species in the future.

14 SOLAR ENERGY↗
NASA NTRS2002

Effects of Volcanic Eruptions on Stratospheric Ozone Recovery

The effects of the stratospheric sulfate aerosol layer associated with the Mt. Pinatubo volcano and future volcanic eruptions on the recovery of the ozone layer is studied with an interactive two-dimensional photochemical model. The time varying chlorine loading and the stratospheric cooling due to increasing carbon dioxide have been taken into account. The computed ozone and temperature changes associated with the Mt. Pinatubo eruption in 1991 agree well with observations. Long model runs out to the year 2050 have been carried out, in which volcanoes having the characteristics of the Mount Pinatubo volcano were erupted in the model at 10-year intervals starting in the year 2010. Compared to a non-volcanic run using background aerosol loading, transient reductions of globally averaged column ozone of 2-3 percent were computed as a result of each of these eruptions, with the ozone recovering to that computed for the non-volcanic case in about 5 years after the eruption. Computed springtime Arctic column ozone losses of from 10 to 18 percent also recovered to the non-volcanic case within 5 years. These results suggest that the long-term recovery of ozone would not be strongly affected by infrequent volcanic eruptions with a sulfur loading approximating Mt. Pinatubo. Sensitivity studies in which the Arctic lower stratosphere was forced to be 4 K and 10 K colder resulted in transient ozone losses of which also recovered to the non-volcanic case in 5 years. A case in which a volcano five times Mt. Pinatubo was erupted in the year 2010 led to maximum springtime column ozone losses of 45 percent which took 10 years to recover to the background case. Finally, in order to simulate a situation in which frequent smaller volcanic eruptions result in increasing the background sulfate loading, a simulation was made in which the background aerosol was increased by 10 percent per year. This resulted in a delay of the recovery of column ozone to 1980 values of more than 10 years.

Rosenfield, Joan E.↗
DOE OSTI2016

Standalone BISON Fuel Performance Results for Watts Bar Unit 1, Cycles 1-3

The Consortium for Advanced Simulation of Light Water Reactors (CASL) is moving forward with more complex multiphysics simulations and increased focus on incorporating fuel performance analysis methods. The coupled neutronics/thermal-hydraulics capabilities within the Virtual Environment for Reactor Applications Core Simulator (VERA-CS) have become relatively stable, and major advances have been made in analysis efforts, including the simulation of twelve cycles of Watts Bar Nuclear Unit 1 (WBN1) operation. While this is a major achievement, the VERA-CS approaches for treating fuel pin heat transfer have well-known limitations that could be eliminated through better integration with the BISON fuel performance code. Several approaches are being implemented to consider fuel performance, including a more direct multiway coupling with Tiamat, as well as a more loosely coupled one-way approach with standalone BISON cases. Fuel performance typically undergoes an independent analysis using a standalone fuel performance code with manually specified input defined from an independent core simulator solution or set of assumptions. This report summarizes the improvements made since the initial milestone to execute BISON from VERA-CS output. Many of these improvements were prompted through tighter collaboration with the BISON development team at Idaho National Laboratory (INL). A brief description of WBN1 and some of the VERA-CS data used to simulate it are presented. Data from a small mesh sensitivity study are shown, which helps justify the mesh parameters used in this work. The multi-cycle results are presented, followed by the results for the first three cycles of WBN1 operation, particularly the parameters of interest to pellet-clad interaction (PCI) screening (fuel-clad gap closure, maximum centerline fuel temperature, maximum/minimum clad hoop stress, and cumulative damage index). Once the mechanics of this capability are functioning, future work will target cycles with known or suspected PCI failures to determine how well they can be estimated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Long-term single-column model intercomparison of diurnal cycle of precipitation over midlatitude and tropical land

General Circulation Models (GCMs) have for decades exhibited difficulties in modelling the diurnal cycle of precipitation (DCP). This issue can be related to inappropriate representation of the processes controlling sub-diurnal phenomena like convection. In this study, 11 single-column versions of GCMs are used to investigate the interactions between convection and environmental conditions, processes that control nocturnal convections, and the transition from shallow to deep convection on a diurnal time-scale. Long-term simulations are performed over two continental land sites: the Southern Great Plains (SGP) in the USA for 12 summer months from 2004 to 2015 and the Manacapuru site at the central Amazon (MAO) in Brazil for two full years from 2014 to 2015. The analysis is done on two regimes: afternoon convective regime and nocturnal precipitation regime. Most models produce afternoon precipitation too early, likely due to the missing transition of shallow-to-deep convection in these models. At SGP, the unified convection schemes better simulate the onset time of precipitation. At MAO, models produce the heating peak in a much lower level compared with observation, indicating too shallow convection in the models. For nocturnal precipitation, models that produce most of nocturnal precipitation all allow convection to be triggered above the boundary layer. This indicates the importance of model capability to detect elevated convection for simulating nocturnal precipitation. Sensitivity studies indicate that (a) nudging environmental variables towards observations has a minor impact on DCP, (b) unified treatment of shallow and deep convection and the capability to capture mid-level convection can help models better capture DCP, and (c) the interactions of the atmosphere with other components in the climate system (e.g. land) are also important for DCP simulations in coupled models. Finally, these results provide long-term statistical insights on which physical processes are essential in climate models to simulate DCP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
NASA NTRS2012

Compression After Impact on Honeycomb Core Sandwich Panels with Thin Facesheets, Part 2: Analysis

A two part research study has been completed on the topic of compression after impact (CAI) of thin facesheet honeycomb core sandwich panels. The research has focused on both experiments and analysis in an effort to establish and validate a new understanding of the damage tolerance of these materials. Part 2, the subject of the current paper, is focused on the analysis, which corresponds to the CAI testings described in Part 1. Of interest, are sandwich panels, with aerospace applications, which consist of very thin, woven S2-fiberglass (with MTM45-1 epoxy) facesheets adhered to a Nomex honeycomb core. Two sets of materials, which were identical with the exception of the density of the honeycomb core, were tested in Part 1. The results highlighted the need for analysis methods which taken into account multiple failure modes. A finite element model (FEM) is developed here, in Part 2. A commercial implementation of the Multicontinuum Failure Theory (MCT) for progressive failure analysis (PFA) in composite laminates, Helius:MCT, is included in this model. The inclusion of PFA in the present model provided a new, unique ability to account for multiple failure modes. In addition, significant impact damage detail is included in the model. A sensitivity study, used to assess the effect of each damage parameter on overall analysis results, is included in an appendix. Analysis results are compared to the experimental results for each of the 32 CAI sandwich panel specimens tested to failure. The failure of each specimen is predicted using the high-fidelity, physicsbased analysis model developed here, and the results highlight key improvements in the understanding of honeycomb core sandwich panel CAI failure. Finally, a parametric study highlights the strength benefits compared to mass penalty for various core densities.

Mcquigg, Thomas D.↗
DOE OSTI2016

Influence of plasma pedestal profiles on access to ELM-free regimes in ITER

The influence of current density and pressure gradient profiles in the pedestal on the access to the regimes free from edge localized modes (ELMs) like quiescent H-mode in ITER is investigated. Using the simulator of MHD modes localized near plasma boundary based on the KINX code, calculations of the ELM stability were performed for the ITER plasma in scenarios 2 and 4 under variations of density and temperature profiles with the self-consistent bootstrap current in the pedestal. Low pressure gradient values at the separatrix, the same position of the density and temperature pedestals and high poloidal beta values facilitate reaching high current density in the pedestal and a potential transition into the regime with saturated large scale kink modes. New version of the localized MHD mode simulator allows one to compute the growth rates of ideal peeling-ballooning modes with different toroidal mode numbers and to determine the stability region taking into account diamagnetic stabilization. The edge stability diagrams computations and sensitivity studies of the stability limits to the value of diamagnetic frequency show that diamagnetic stabilization of the modes with high toroidal mode numbers can help to access the quiescent H-mode even with high plasma density but only with low pressure gradient values at the separatrix. The limiting pressure at the top of the pedestal increases for higher plasma density. With flat density profile the access to the quiescent H-mode is closed even with diamagnetic stabilization taken into account, while toroidal mode numbers of the most unstable peeling-ballooning mode decrease from n = 10−40 to n = 3−20.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Major improvements on the longwave radiative interactions between surface and clouds in the polar regions in atmospheric global circulation model

Motivations: One of the most important tasks in the study of physical and socioeconomic aspects of climate change is the simulation of polar climate and how it changes in response to global warming. One import component of the polar energy budget (and thus of polar climate) is the radiative interaction between surface and clouds. As of today, the vast majority of state-of-the-art global climate models (GCMs), including all three US flagship GCMs (NCAR CESM, GFDL CM3, and GISS E2-H/E2-R), still assume constant surface emissivity overall spectral bands in their longwave radiation treatment. Similarly, a majority of state-of-the-art GCMs assume non-scattering clouds in the longwave spectrum, including both the NCAR CESM and GFDL CM3. The issues of ignoring spectral variation in surface emissivity and scattering by clouds manifest themselves in the far-infrared (IR) over the polar continents because (1) there is so little water vapor in such areas that surface far-IR emission can reach clouds; and (2) the imaginary part of the refractive index of ice has a local minimum over 350-550 cm-1, which implies possible strong scattering effects over this spectral region both within and between the snow surface and ice clouds. Our sensitivity studies using CloudSat-retrieved hydrometeor profiles over the high-elevation Antarctic continent show that, in winter, an appropriate inclusion of snow surface spectral emissivity and ice cloud scattering in radiative transfer calculations can noticeably reduce the monthly-mean surface net downward far-IR flux and net atmospheric far-IR emission over the entire region. For the far-IR alone, the magnitude of these effects is ~1Wm-2 or even larger. Projection Objectives: We propose to carry out the following studies: (1) Develop or update schemes that can be used in the CESM to compute the optical properties of the snow surface and ice clouds consistently across the longwave and the shortwave spectrum; (2) Modify the longwave radiation scheme in the CESM such that it can incorporate surface spectral emissivity and scattering of clouds; (3) Validate schemes developed for (1) and (2) against measurements or benchmark models; (4) Evaluate the effects of these new treatments on the simulated polar energy budget and climate by the CESM. Technical Approaches:1.Use a detailed snowpack radiative transfer model to simulate spectrally resolved emissivity of snow surfaces as a function of snow effective radius. The model will be validated against available high spectral-resolution measurements of snow spectral emissivity (ASTER spectra library v2.0). Currently the CESM prognoses vertically-resolved snow effective radius, and snow shortwave albedo is parameterized based on snow effective radius. Our treatment of snow spectral emissivity here will be consistent with the treatment of snow albedo in the shortwave and seamlessly connected with the snow physical module in the CESM. 2. Update the parameterization of ice cloud optical properties in the CESM, for both longwave and shortwave, based on the latest ice-cloud scattering database developed by co-I Prof. Yang. 3. Modify RRTMG_LW, the current longwave radiation scheme in the CESM, to include the surface spectral emissivity and the scattering of clouds in the longwave calculation. Validate modified scheme against a line-by-line radiative transfer model with multiple scattering capability (LBLRTM/CHARTS). 4. Implement products from Steps 1-3 into the CESM and conduct simulations to evaluate their effects on the simulated polar climate and energy budget. These runs will be compared to control runs and evaluated using routine model-observation diagnostics, including data from the ARM NSA (North Slope of Alaska) site. Moreover, a novel aspect in this step is that we will evaluate the simulated longwave band-by-band fluxes and cloud radiative effect against observational counterparts over both polar continents, using a product recently derived by the lead PI Prof. Huang from collocated AIRS and CERES satellite observations from 2002 to 2013. Benefits and outcomes: This project will significantly improve the longwave radiation scheme in the CESM and its successors and provide more realistic treatment of radiative interactions between polar surfaces and clouds. In the future, it will also be straightforward to extend this framework to include spectral emissivity of all surface types. This is responsive to identified DoE Earth System Modeling Science Drivers to improve the performance of the fully coupled climate system modeling and to help address questions on how rapid changes in cryospheric systems interact with the climate system.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

The Role of Clear Sky Identification in the Study of Cloud Radiative Effects: Combine Analysis from ISCCP and the Scanner of Radiation Budget (ScaRaB)

Since the effect of clouds on the earth's radiation balance is often estimated as the difference of net radiative fluxes at the top of the atmosphere between all situations and monthly averaged clear sky situations of the same regions, a reliable identification of clear sky is important for the study of cloud radiative effects. The Scanner for Radiation Balance (ScaRaB) radiometer on board the Russian Meteor-3/7 satellite provided earth radiation budget observations from March 1994 to February 1995 with two ERBE-Re broad-band longwave and shortwave channels. Two narrow-band channels, in the infrared atmospheric window and in the visible band, have been added to the ScaRaB instrument to improve the cloud scene identification. The International Satellite Cloud Climatology Project (ISCCP) method for cloud detection and determination of cloud and surface properties uses the same narrow-band channels as ScaRaB, but is employed to a collection of measurements at a better spatial resolution of about 5 km. By applying the original ISCCP algorithms to the ScaRaB data, the clear sky frequency is about 5% lower than the one over quasi-simultaneous original ISCCP data, an indication that the ISCCP cloud detection is quite stable. However, one would expect an about 10 to 20% smaller clear sky occurrence over the larger ScaRaB pixels. Adapting the ISCCP algorithms to the reduced spatial resolution of 60 km and to the different time sampling of the ScaRaB data leads therefore to a reduction of a residual cloud contamination. A sensitivity study with time-space collocated ScaRaB and original ISCCP data at a spatial resolution of 1deg longitude x 1deg latitude shows that the effect of clear sky identification method plays a higher role on the clear sky frequency and therefore on the statistics than on the zonal mean values of the clear sky fluxes. Nevertheless, the zonal outgoing longwave fluxes corresponding to ERBE clear sky are in general about 2 to 10 W/sq m higher than those obtained from the ScaRaB adapted ISCCP clear sky identifications. The latter are close to (about 1 W/sq m higher) fluxes corresponding to clear sky regions from original ISCCP data, whereas ScaRaB clear sky LW fluxes obtained with the original ISCCP identification lie about 1 to 2 W/sq m below. Especially in the tropics where water vapor abundance is high, the ERBE clear sky LW fluxes seem to be systematically overestimated by about 4 W/sq m, and SW fluxes are lower by about 5 to 10 W/sq m. However, the uncertainty in the analysis of monthly mean zonal cloud radiative effects is also produced by the low frequency of clear sky occurrence, illustrated when averaging over pixels or even over regions of 4deg longitude x 5deg latitude, corresponding to the spatial resolution of General Circulation Models. The systematic bias in the clear sky fluxes is not reflected in the zonal cloud radiative effects, because the clear sky regions selected by the different algorithms can occur in different geographic regions with different cloud properties.

Rossow, W. B.↗
DOE OSTI2026

Geothermal Direct-Use Applications for the District Energy System in Bucharest, Romania

The city of Bucharest, Romania, hosts the second-largest district energy system (DES) in the world. Geothermal resources can be considered as a supplementary heat source to support the demand for domestic hot water and space heating in the winter and shoulder seasons. The National Laboratory of the Rockies (NLR) has conducted a study that considers geothermal energy to serve a fraction of the existing district heating network operated by Electrocentrale Bucure?ti (ELCEN), the utility operating the DES. Lower Cretaceous and Jurassic limestones make up the main geothermal aquifer underlying Bucharest, which hosts temperatures suitable for district heating (up to 90 degrees C to the north of the city). Anomalous geothermal gradients have been observed to the north of the city, where a pumped well has produced 82 degrees C brine at the wellhead to feed the Therme Bucharest Spa. An anomalous gradient has also been reported to the southeast of the city (35 degrees C/km). NLR modeled the building heating loads of a small portion of the DES (a block of nine prototypical buildings) in its Urban Renewable Building and Neighborhood Optimization (URBANopt ) platform. To simulate meeting a baseload benchmark of 20 MWth deliverable to a small portion of the DES, the NLR team used GEOPHIRES to model production scenarios for (1) hydrothermal systems coupled with heat pumps targeting the main geothermal aquifer in the north, (2) enhanced geothermal systems targeting hot dry rock in the southeast, and (3) huff-and-puff systems targeting a gradient of 25 degrees C/km. Finally, NLR conducted a high-level sensitivity study around the techno-economics of these systems. The outcomes of this work indicate that hot dry rock geothermal resources that can deliver at least 90 degrees C hot water to the Geothermal District Energy System (GeoDES) offer a possible solution for supplemental geothermal heat delivered to the existing DES.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

On the Realism of the Rain Microphysics Representation of a Squall Line in the WRF Model. Part I: Evaluation with Multifrequency Cloud Radar Doppler Spectra Observations

This study investigates how multifrequency cloud radar observations can be used to evaluate the representation of rain microphysics in the WRF Model using two bulk microphysics schemes. A squall line observed over Oklahoma on 12 June 2011 is used as a case study. A recently developed retrieval technique combining observations of two vertically pointing cloud radars provides quantitative description of the drop size distribution (DSD) properties of the transition and stratiform regions of the squall-line system. For the first time, the results of this multifrequency cloud radar retrieval are compared to more conventional retrievals from a nearby polarimetric radar, and a supplementary result of this work is that this new methodology provides a much more detailed description of the DSD vertical and temporal variations. While the extent and evolution of the squall line is well reproduced by the model, the 1-h low-reflectivity transition region is not. In the stratiform region, simulations with both schemes are able to reproduce the observed downdraft and the associated significative subsaturation below the melting level, but with a slight overestimation of the relative humidity. Under this subsaturated air, the simulated rain mixing ratio continuously decreases toward the ground, in agreement with the observations. Conversely, the profiles of the mean volume diameter and the concentration parameter of the DSDs are not well reproduced. These discrepancies pinpoint at an issue in the representation of rain microphysics. The companion paper, investigates the sources of the biases in the microphysics processes in the rain layer by performing numerical sensitivity studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗