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 523 records · Page 29

Potential Nighttime Contamination of CERES Clear-sky Field of View by Optically Thin Cirrus during the CRYSTAL-FACE Campaign

We investigate the outgoing broadband longwave (LW, 5 to approx. 200 microns) and window (WIN, 8 to approx. 12 microns) channel radiances at the top of atmosphere (TOA) under clear-sky conditions, using data acquired by the Cloud and the Earth s Radiant Energy System (CERES) and Moderate-Resolution Imaging Spectroradiometer (MODIS) instruments onboard the NASA Terra satellite platform. In this study, detailed analyses are performed on the CERES Single Scanner Footprint TOA/Surface Fluxes and Clouds product to understand the radiative effect of thin cirrus. The data are acquired over the Florida area during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers Florida Area Cirrus Experiment (CRYSTAL-FACE) field program. Of particular interest is the anisotropy associated with the radiation field. Measured CERES broadband radiances are compared to those obtained from rigorous radiative transfer simulations. Analysis of results from this comparison indicates that the simulated radiances tend to be larger than their measured counterparts, with differences ranging from 2.1% to 8.3% for the LW band and from 1.7% to 10.6% for the WIN band. The averaged difference in radiance is approximately 4% for both the LW and WIN channels. A potential cause for the differences could be the presence of thin cirrus (i.e., optically thin ice clouds with visible optical thicknesses smaller than approximately 0.3). The detection and quantitative analysis of these thin cirrus clouds are challenging even with sophisticated multispectral instruments. While large differences in radiance between the CERES observations and the theoretical calculations are found, the corresponding difference in the anisotropic factors is very small (0.2%). Furthermore, sensitivity studies show that the influence due to a 1 K bias of the surface temperature on the errors of the LW and WIN channel radiances is of the same order as that associated with a 2% bias of the surface emissivity. The LW and WIN errors associated with a 5% bias of water vapor amount in the lower atmosphere in conjunction with a 50% bias of water vapor amount in the upper atmosphere is similar to that of a 1 K bias of the vertical temperature profile. Even with the uncertainties considered for these various factors, the simulated LW and WIN radiances are still larger than the observed radiances if thin cirrus clouds are excluded.

Lee, Yong-Keun↗

Seat Capacity Selection for an Advanced Short-Haul Aircraft Design

A study was performed to determine the target seat capacity for a proposed advanced short-haul aircraft concept projected to enter the fleet by 2030. This analysis projected the potential demand in the U.S. for a short-haul aircraft using a transportation theory approach, rather than selecting a target seat capacity based on recent industry trends or current market demand. A transportation systems model was used to create a point-to-point network of short-haul trips and then predict the number of annual origin-destination trips on this network. Aircraft of varying seat capacities were used to meet the demand on this network, assuming a single aircraft type for the entire short-haul fleet. For each aircraft size, the ticket revenue and operational costs were used to calculate a total market profitability metric for all feasible flights. The different aircraft sizes were compared, based on this market profitability metric and also the total number of annual round trips and markets served. Sensitivity studies were also performed to determine the effect of changing the aircraft cruise speed and maximum trip length. Using this analysis, the advanced short-haul aircraft design team was able to select a target seat capacity for their design.

Marien, Ty V.↗

SAR terrain classifier and mapper of biophysical attributes

In preparation for the launch of SIR-C/X-SAR and design studies for future orbital SAR, a program has made considerable progress in the development of an SAR terrain classifier and algorithms for quantification of biophysical attributes. The goal of this program is to produce a generalized software package for terrain classification and estimation of biophysical attributes and to make this package available to the larger scientific community. The basic elements of the SAR (Synthetic Aperture Radar) terrain classifier are outlined. An SAR image is calibrated with respect to known system and processor gains and external targets (if available). A Level 1 classifier operates on the data to differentiate: urban features, surfaces and tall and short vegetation. Level 2 classifiers further subdivide these classes on the basis of structure. Finally, biophysical and geophysical inversions are applied to each class to estimate attributes of interest. The process used to develop the classifiers and inversions is shown. Radar scattering models developed from theory and from empirical data obtained by truck-mounted polarimeters and the JPL AirSAR are validated. The validated models are used in sensitivity studies to understand the roles of various scattering sources (i.e., surface trunk, branches, etc.) in determining net backscatter. Model simulations of sigma (sup o) as functions of the wave parameters (lambda, polarization and angle of incidence) and the geophysical and biophysical attributes are used to develop robust classifiers. The classifiers are validated using available AirSAR data sets. Specific estimators are developed for each class on the basis of the scattering models and empirical data sets. The candidate algorithms are tested with the AirSAR data sets. The attributes of interest include: total above ground biomass, woody biomass, soil moisture and soil roughness.

Ulaby, Fawwaz T.↗

Comparison of Different Global Information Sources Used in Surface Radiative Flux Calculation: Radiative Properties of the Surface

Direct estimates of surface radiative fluxes that resolve regional and weather-scale variabilty over the whole globe with reasonable accuracy have only become possible with the advent of extensive global, mostly satellite, datasets within the past couple of decades. The accuracy of these fluxes, estimated to be about 10-15 W per square meter is largely limited by the accuracy of the input datasets. The leading uncertainties in the surface fluxes are no longer predominantly induced by clouds but are now as much associated with uncertainties in the surface and near-surface atmospheric properties. This study presents a fuller, more quantitative evaluation of the uncertainties for the surface albedo and emissivity and surface skin temperatures by comparing the main available global datasets from the Moderate-Resolution Imaging Spectroradiometer product, the NASA Global Energy and Water Cycle Experiment Surface Radiation Budget project, the European Centre for Medium-Range Weather Forecasts, the National Aeronautics and Space Administration, the National Centers for Environmental Prediction, the International Satellite Cloud Climatology Project (ISCCP), the Laboratoire de Meteorologie Dynamique, NOAA/NASA Pathfinder Advanced Very High Resolution Radiometer project, NOAA Optimum Interpolation Sea Surface Temperature Analysis and the Tropical Rainfall Measuring Mission (TRMM) Microwave Image project. The datasets are, in practice, treated as an ensemble of realizations of the actual climate such that their differences represent an estimate of the uncertainty in their measurements because we do not possess global truth datasets for these quantities. The results are globally representative and may be taken as a generalization of our previous ISCCP-based uncertainty estimates for the input datasets. Surface properties have the primary role in determining the surface upward shortwave (SW) and longwave (LW) flux. From this study, the following conclusions are obtained. Although land surface albedos in the near near-infrared remain poorly constrained (highly uncertain), they do not cause too much error in total surface SW fluxes; the more subtle regional and seasonal variations associated with vegetation and snow are still on doubt. The uncertainty of the broadband black-sky SW albedo for land surface from this study is about 7%, which can easily induce 5-10 W per square meter uncertainty in (upwelling) surface SW flux estimates. Even though available surface (broadband) LW emissivity datasets differ significantly (3%-5% uncertainty), this disagreement is confined to wavelengths greater than 20 micrometers so that there is little practical effect (1-3 W per square meters) on the surface upwelling LW fluxes. The surface skin temperature is one of two leading factors that cause problems with surface LW fluxes. Even though the differences among the various datasets are generally only 2-4 K, this can easily cause 10-15 W per square meter uncertainty in calculated surface (upwelling) LW fluxes. Significant improvements could be obtained for surface LW flux calculations by improving the retrievals of (in order of decreasing importance): (1) surface skin temperature, (2) surface air and near-surface-layer temperature, (3) column precipitable water amount and (4) broadband emissivity. And for surface SW fluxes, improvements could be obtained (excluding improved cloud treatment) by improving the retrievals of (1) aerosols (from our sensitivity studies but not discussed in this work), and (2) surface (black-sky) albedo, of which, NIR part of the spectrum has much larger uncertainty.

Zhang, Yuanchong↗

Ocean Carbon Flux, Transport, and Burial Within the Western and Eastern US Coastal Zones

This project has been to develop and apply a regional. eddy-resolving circulation and biogeochemistry model of both the western and eastern U.S. coastal regions, capable of simulating the processes that control the carbon cycle. Validation has been by statistical comparison with analyses from various satellite measurements, including those from EOS sensors, as well as from in situ measurements. Sensitivity studies were carried out to investigate how the coastal ecosystem and biogeochemical cycles respond to changes in climate, large-scale eutrophication from indus- trial pollution, and other anthropogenic induced changes. The research has been conducted in collaboration with research groups at UCLA. NASA/GSFC (Wallops), Rutgers, and SIO. Overall. the project was focused on several key modeling issues, each of which tie back into completing the primary task of developing a coastal carbon model for both the eastern and western US. coasts. Individual groups within the entire program are still collaborating to address these specific tasks. These include: implementation of the coupled circulation/biogeochemical model within the U.S. West Coast. including high-resolution, embedded subdomains for the Southern California Bight and Monterey Bay region; development of a biogeochemical model with resolved carbon, nitrogen and oxygen cycles; development of data assimilation techniques for use of satellite data sets; reconfiguration of the model domain to U.S. East Coast; development of coastal forcing fields: development of methods to compare the model against remotely sensed data; and, the test of model sensitivity to environmental conditions. Below, we present a summary of the progress made toward achieving these soak. Because this has been a multi-institutional, collaborative effort, we note the groups involved with particular activities.

McWilliams, James C.↗

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗

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↗

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↗

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↗

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.↗

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↗

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↗

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↗

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.↗

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.↗

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↗

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↗

Investigation of powertrain system decarbonization using electrically assisted turbocharging and hybridization in off-road vehicles

In recent years, the causes, effects, and existential threat of a global anthropogenic climate shift have drawn significant attention and stimulated mitigation efforts from all genres of scientific, political, and industrial bodies. The greenhouse effect and the detrimental environmental impact of excess greenhouse gas (GHG) emissions, like carbon dioxide, is well studied and targeted as the primary culprit for aversive action. However, some fields exist where reducing GHG emissions is met by considerable challenges. One such field is the production and operations of off-road vehicles. Applications of these vehicles are highly diverse and are often characterized by rugged, especially transient, and power intensive duty cycles that make one-fits-all vehicle configuration solutions impractical. This research surveys existing literature to identify the modern technologies packages and challenges that face development and configuration of powertrain systems which currently navigate the regulated emissions and unique duty cycle requirements of this market space. A novel powertrain concept is proposed and evaluated with respect to the pinnacle objectives of load performance improvement, criteria pollutant reduction, and reducing the life cycle GHG intensity of its operation. A sure-fire path to vehicular GHG reduction is through improving the fuel efficiency of internal combustion engines (ICEs), an entrenched component of the off-road vehicle sector. Unfortunately, this is more easily said than done. An approach that has proven successful in this endeavor is engine downsizing and turbocharging, where a larger engine is replaced with a smaller one with added air system boosting via turbocharger to reduce frictional and pumping losses while also enabling access to additional fuel energy. However, practical realization of these potential benefits is often impeded by the transient response capability of the smaller engine across the operating space. For this reason, downsized ICE powertrains have turned to electrified forced induction systems (EFISs) for a decoupling of exhaust energy and engine speed from boost capability. Platformed on 48V hybrid technology, these systems introduce the need for more sophisticated controls around the engine gas-exchange process for the management of boost performance and exhaust gas emissions. On this account, simulation studies utilizing a GT-POWER model of a turbocharged 4.5 L engine outfit with an EFIS using an electrically driven compressor (eBooster®) are conducted to provide insight into the performance of this air handling architecture on an off-road engine. The results show that improvement in transient torque response time in sync with reduction in engine-out soot and NOx emissions are possible with an engine recalibration that leverages the transient air-fuel ratio authority of the eBooster®. Benefits are further demonstrated when duty cycle simulations of the 48V mild-hybrid engine concept are exercised, showing an acute decrease in cumulative fuel usage and soot production. Powertrain hybridization is another technological pathway achieving pronounced GHG reduction successes in modern on-road vehicle applications through integration of Li-ion battery technology. In the off-road vehicle segment, a review of available literature concludes that hybridized architectures are present but generally lack the depth of technologies that have both high specific energy and power capabilities, and thus are limited in their inclusion of Li-ion batteries for ICE assistance and enhanced energy storage capability. Therefore, building on the mild-hybrid engine results, the downsized and eBoosted engine concept was integrated into a larger high-voltage battery-hybrid series-electric powertrain system. Hybrid powertrain parameter sensitivity studies were carried out in a numerical charge-sustaining framework, providing novel insights into power flows between the battery and the engine and how their respective capabilities and operation contribute to GHG and criteria pollutant emissions of diverse duty cycles. Application of supervisory power management introduced robustness into the power sourcing and battery SOC control process and showed that optimum specification of battery properties can yield synergies between GHG emissions and battery energy capacity. Furthermore, examination of recent literature on Li-ion batteries has shown that pack manufacturing is a highly energy intensive process, thereby producing considerable quantities of GHGs that scale with energy storage capacity. Also scaling with a battery’s energy storage capacity is its investment cost. In consideration of these factors, an inclusive technoeconomic and GHG life cycle analysis is conducted. This analysis systematically compares the carbon footprint and total cost of ownership associated with the proposed hybrid powertrain concept to reference and alternative powertrain configurations, facilitating a thorough evaluation of the decarbonization effectiveness and economic viability.

99 GENERAL AND MISCELLANEOUS↗