Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stratified”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Planetary Geophysics and Tectonics

Research supported by grant NAGW-1928 has addressed a variety of problems related to planetary evolution. One important focus has been on questions related to the role of chemical buoyancy in planetary evolution with application to both Venus and the Moon. We have developed a model for the evolution of the Moon (Hess and Parmentier, 1995) in which dense, highly radioactive, late stage magma ocean cumulates sink forming a core. This core heats the overlying, chemically layered mantle giving rise to a heated, chemically well-mixed layer that thickens with time. This Mixed layer eventually becomes hot enough and thick enough that its top begins to melt at a pressure low enough that melt is buoyant, thus creating mare basalts from a high pressure source of the correct composition and at an appropriate time in lunar evolution. In work completed during the last year, numerical experiments on convection in a chemically stably stratified fluid layer heated from below have been completed. These results show us how to calculate the evolution of a mixed layer in the Moon, depending on the heat production in the ilmenite- cumulate core and the chemical stratification of the overlying mantle. Chemical stratification of the mantle after its initial differentiation is would trap heat in the deep interior and prevent the rapid rise of plumes with accompanying volcanism. This trapping of heat in the interior can explain the thickness of the lunar lithosphere as a function of time as well as the magmatic evolution. We show that heat transported to the base of the lithosphere at a rate determined by current estimates of radioactivity in the Moon would not satisfy constraints on elastic lithosphere thickness from tectonic feature associated with basin loading. Trapping heat at depth by a chemically stratified mantle may also explain the absence of global compressional features on the surface that previous models predict for an initially hot lunar interior. For Venus, we developed a model in which the chemical buoyancy of crust and a depleted mantle layer stabilizes the lithosphere for long periods of time and provides a mechanism of episodic planetary evolution (Parmentier and Hess, 1992). Continued thickening of a residual depleted mantle layer eventually suppresses pressure release melting and the creation of depleted mantle. Continued cooling then allows the lithosphere to become heavier than the underlying hotter, undepleted mantle. This repeated instability can occur on time scales appropriate for episodic global resurfacing on Venus. We have also examined the role of the gabbro-eclogite phase transformation on crust and lithosphere stability and as a mechanism of crustal recycling in the absence of plate tectonics. Our work thus far concentrates on the scale of instability that would occur due to cooling or crustal thickening associated with horizontal shortening. Whether repeated overturn can explain the evolution of Venus depends in part on whether sufficient heat transfer can occur between overturns and on constraints provided by understanding observed surface features and evolution.

Parmentier, E. M.↗
NASA NTRS2002

Partial Wave Analysis of Coupled Photonic Structures

The very high quality factors sustained by microcavity optical resonators are relevant to applications in wavelength filtering, routing, switching, modulation, and multiplexing/demultiplexing. Increases in the density of photonic elements require that attention be paid to how electromagnetic (EM) coupling modifies their optical properties. This is especially true when cavity resonances are involved, in which case, their characteristics may be fundamentally altered. Understanding the optical properties of microcavities that are near or in contact with photonic elements---such as other microcavities, nanostructures, couplers, and substrates---can be expected to advance our understanding of the roles that these structures may play in VLSI photonics, biosensors and similar device technologies. Wc present results from recent theoretical studies of the effects of inter- and intracavity coupling on optical resonances in compound spherical particles. Concentrically stratified spheres and bispheres constituted from homogeneous and stratified spheres are subjects of this investigation. A new formulation is introduced for the absorption of light in an arbitrary layer of a multilayered sphere, which is based on multiple reflections of the spherical partial waves of the Lorenz-Mie solution for scattering by a sphere. Absorption efficiencies, which can be used to profile cavity resonances and to infer fluorescence yields or the onset of nonlinear optical processes in the microcavities, are presented. Splitting of resonances in these multisphere systems is paid particular attention, and consequences for photonic device development and possible performance enhancements through carefully designed architectures that exploit EM coupling are considered.

Fuller, Kirk A.↗
NASA NTRS2003

Cox Proportional Hazards Models for Modeling the Time to Onset of Decompression Sickness in Hypobaric Environments

In this paper we fit Cox proportional hazards models to a subset of data from the Hypobaric Decompression Sickness Databank. The data bank contains records on the time to decompression sickness (DCS) and venous gas emboli (VGE) for over 130,000 person-exposures to high altitude in chamber tests. The subset we use contains 1,321 records, with 87% censoring, and has the most recent experimental tests on DCS made available from Johnson Space Center. We build on previous analyses of this data set by considering more expanded models and more detailed model assessments specific to the Cox model. Our model - which is stratified on the quartiles of the final ambient pressure at altitude - includes the final ambient pressure at altitude as a nonlinear continuous predictor, the computed tissue partial pressure of nitrogen at altitude, and whether exercise was done at altitude. We conduct various assessments of our model, many of which are recently developed in the statistical literature, and conclude where the model needs improvement. We consider the addition of frailties to the stratified Cox model, but found that no significant gain was attained above a model that does not include frailties. Finally, we validate some of the models that we fit.

Thompson, Laura A.↗
NASA NTRS2003

Coherent Effects in Tiny Optics: Tunneling Through the Looking Glass

I will discuss two types of one-dimensional photonic bandgap (PBG) effects that can arise in systems of coupled spherical resonators: (1) nearly-free-photon Fabry-Perot photonic bands that arise in quarter-wave concentrically stratified spheres and, (2) tight- binding photonic bands that arise in weakly-coupled mutually-resonant spheres as a result of whispering-gallery mode splitting. These effects can be derived directly from Mie theory, in a more straightforward manner, by exploiting an analogy with stratified planar systems. For odd numbers of mutually-resonant lossless coupled ring resonators, the circulating intensity can increase exponentially with the number of resonators, which can potentially be exploited for the development of advanced sensors. For even numbers of resonators, mode splitting and classical destructive interference lead to a cancellation of absorption and slow light on-resonance, reminiscent of electromagnetic induced transparency. The analogy between these coherent photon trapping effects and population trapping in an atomic system will be explored.

Smith, David D.↗
NASA NTRS2003

Validation of an In-Water, Tower-Shading Correction Scheme

Large offshore structures used for the deployment of optical instruments can significantly perturb the intensity of the light field surrounding the optical measurement point, where different portions of the visible spectrum are subject to different shadowing effects. These effects degrade the quality of the acquired optical data and can reduce the accuracy of several derived quantities, such as those obtained by applying bio-optical algorithms directly to the shadow-perturbed data. As a result, optical remote sensing calibration and validation studies can be impaired if shadowing artifacts are not fully accounted for. In this work, the general in-water shadowing problem is examined for a particular case study. Backward Monte Carlo (MC) radiative transfer computations- performed in a vertically stratified, horizontally inhomogeneous, and realistic ocean-atmosphere system are shown to accurately simulate the shadow-induced relative percent errors affecting the radiance and irradiance data profiles acquired close to an oceanographic tower. Multiparameter optical data processing has provided adequate representation of experimental uncertainties allowing consistent comparison with simulations. The more detailed simulations at the subsurface depth appear to be essentially equivalent to those obtained assuming a simplified ocean-atmosphere system, except in highly stratified waters. MC computations performed in the simplified system can be assumed, therefore, to accurately simulate the optical measurements conducted under more complex sampling conditions (i.e., within waters presenting moderate stratification at most). A previously reported correction scheme, based on the simplified MC simulations, and developed for subsurface shadow-removal processing of in-water optical data taken close to the investigated oceanographic tower, is then validated adequately under most experimental conditions. It appears feasible to generalize the present tower-specific approach to solve other optical sensor shadowing problems pertaining to differently shaped deployment platforms, and also including surrounding structures and instrument casings.

Hooker, Stanford B.↗
NASA NTRS1998

T wave alternans as a predictor of recurrent ventricular tachyarrhythmias in ICD recipients: prospective comparison with conventional risk markers

INTRODUCTION: The current standard for arrhythmic risk stratification is electrophysiologic (EP) testing, which, due to its invasive nature, is limited to patients already known to be at high risk. A number of noninvasive tests, such as determination of left ventricular ejection fraction (LVEF) or heart rate variability, have been evaluated as additional risk stratifiers. Microvolt T wave alternans (TWA) is a promising new risk marker. Prospective evaluation of noninvasive risk markers in low- or moderate-risk populations requires studies involving very large numbers of patients, and in such studies, documentation of the occurrence of ventricular tachyarrhythmias is difficult. In the present study, we identified a high-risk population, recipients of an implantable cardioverter defibrillator (ICD), and prospectively compared microvolt TWA with invasive EP testing and other risk markers with respect to their ability to predict recurrence of ventricular tachyarrhythmias as documented by ICD electrograms. METHODS AND RESULTS: Ninety-five patients with a history of ventricular tachyarrhythmias undergoing implantation of an ICD underwent EP testing, assessment of TWA, as well as determination of LVEF, baroreflex sensitivity, signal-averaged ECG, analysis of 24-hour Holter monitoring, and QT dispersion from the 12-lead surface ECG. The endpoint of the study was first appropriate ICD therapy for electrogram-documented ventricular fibrillation or tachycardia during follow-up. Kaplan-Meier survival analysis revealed that TWA (P < 0.006) and LVEF (P < 0.04) were the only significant univariate risk stratifiers. EP testing was not statistically significant (P < 0.2). Multivariate Cox regression analysis revealed that TWA was the only statistically significant independent risk factor. CONCLUSIONS: Measurement of microvolt TWA compared favorably with both invasive EP testing and other currently used noninvasive risk assessment methods in predicting recurrence of ventricular tachyarrhythmias in ICD recipients. This study suggests that TWA might also be a powerful tool for risk stratification in low- or moderate-risk patients, and needs to be prospectively evaluated in such populations.

Non-NASA Center↗
NASA NTRS2008

Model Effects on GLAS-Based Regional Estimates of Forest Biomass and Carbon

ICESat/GLAS waveform data are used to estimate biomass and carbon on a 1.27 million sq km study area. the Province of Quebec, Canada, below treeline. The same input data sets and sampling design are used in conjunction with four different predictive models to estimate total aboveground dry forest biomass and forest carbon. The four models include nonstratified and stratified versions of a multiple linear model where either biomass or (square root of) biomass serves as the dependent variable. The use of different models in Quebec introduces differences in Provincial biomass estimates of up to 0.35 Gt (range 4.942+/-0.28 Gt to 5.29+/-0.36 Gt). The results suggest that if different predictive models are used to estimate regional carbon stocks in different epochs, e.g., y2005, y2015, one might mistakenly infer an apparent aboveground carbon "change" of, in this case, 0.18 Gt, or approximately 7% of the aboveground carbon in Quebec, due solely to the use of different predictive models. These findings argue for model consistency in future, LiDAR-based carbon monitoring programs. Regional biomass estimates from the four GLAS models are compared to ground estimates derived from an extensive network of 16,814 ground plots located in southern Quebec. Stratified models proved to be more accurate and precise than either of the two nonstratified models tested.

Nelson, Ross↗
NASA NTRS2010

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗
NASA NTRS2010

Model Effects on GLAS-Based Regional Estimates of Forest Biomass and Carbon

Ice, Cloud, and land Elevation Satellite (ICESat) / Geosciences Laser Altimeter System (GLAS) waveform data are used to estimate biomass and carbon on a 1.27 X 10(exp 6) square km study area in the Province of Quebec, Canada, below the tree line. The same input datasets and sampling design are used in conjunction with four different predictive models to estimate total aboveground dry forest biomass and forest carbon. The four models include non-stratified and stratified versions of a multiple linear model where either biomass or (biomass)(exp 0.5) serves as the dependent variable. The use of different models in Quebec introduces differences in Provincial dry biomass estimates of up to 0.35 G, with a range of 4.94 +/- 0.28 Gt to 5.29 +/-0.36 Gt. The differences among model estimates are statistically non-significant, however, and the results demonstrate the degree to which carbon estimates vary strictly as a function of the model used to estimate regional biomass. Results also indicate that GLAS measurements become problematic with respect to height and biomass retrievals in the boreal forest when biomass values fall below 20 t/ha and when GLAS 75th percentile heights fall below 7 m.

Nelson, Ross F.↗
NASA NTRS2010

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗
NASA NTRS2005

Using Cloud-to-Ground Lightning Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourn, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in central Florida, especially during the warm season months of May-September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC (0700 AM EST) each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to increase consistency between forecasters while enabling them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) for which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities from National Lightning Detection Network (NLDN) data. The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at the two organizations provided this data and supporting software for the work performed by the AMU. The densities were first stratified by flow regime, then by time in 1-, 3-, 6-, 12-, and 24-hour increments while maintaining the 2.5 km x 2.5 km grid resolution. A CG frequency of occurrence was calculated for each stratification and grid box by counting the number of days with lightning and dividing by the total number of days in the data set. New CG strike densities were calculated for each stratification and grid box by summing the strike number values over all warm seasons, then normalized by dividing the summed values by the number of lightning days. This makes the densities conditional on whether lightning occurred. The frequency climatology values will be used by forecasters as proxy inputs for lightning prObability, while the density climatology values will be used for CG amount. In addition to the benefits outlined above, these climatologies will provide improved temporal and spatial resolution, expansion of the lightning threat area to include adjacent coastal waters, and potential to extend the forecast to include the day-2 period. This presentation will describe the lightning threat index map, discuss the work done to create the maps initialized with climatological guidance, and show examples of the climatological CG lightning densities and frequencies of occurren based on flow regime.

Lambert, Winnie↗
NASA NTRS2015

Laser Radar Through the Window (LRTW) Coordinate Correction Method

A method for corrections of measurements of points of interests measured by beams of radiation propagating through stratified media including performance of ray-tracing of at least one ray lunched from a metrology instrument in a direction of an apparent point of interest, calculation a path length of the ray through stratified medium, and determination of coordinates of true position of the point interest using the at least one path length and the direction of propagation of the ray.

Hayden, Joseph Ethan↗
NASA NTRS2014

Global Characteristics of Porosity and Density Stratification Within the Lunar Crust from GRAIL Gravity and Lunar Orbiter Laser Altimeter Topography Data

The Gravity Recovery and Interior Laboratory (GRAIL) mission is providing unprecedentedly high-resolution gravity data. The gravity signal in relation to topography decreases from 100 km to 30 km wavelength, equivalent to a uniform crustal density of 2450 kg/cu m that is 100 kg/cu m smaller than the density required at 100 km. To explain such frequency-dependent behavior, we introduce rock compaction models under lithostatic pressure that yield radially stratified porosity (and thus density) and examine the depth extent of porosity. Our modeling and analysis support the assertion that the crustal density must vary from surface to deep crust by up to 500 kg/cu m. We found that the surface density of mega regolith is around 2400 kg/cu m with an initial porosity of 10-20%, and this porosity is eliminated at 10-20 km depth due to lithostatic overburden pressure. Our stratified density models provide improved fits to both GRAIL primary and extended mission data.

Density↗
NASA NTRS2021

Mangrove Carbon Stocks in Pongara National Park, Gabon

Mangroves are recognized for their valued ecosystem services to coastal areas, and the functional linkages between those services and ecosystem carbon stocks have been established. However, spatially explicit inventories are necessary to facilitate management and protection of mangroves, as well as providing a foundation for payment for ecosystem service programs such as REDD+. We conducted an inventory of carbon stocks in mangroves within Pongara National Park (PNP), Gabon using a stratified random sampling design based on forest canopy height derived from TanDEM-X remote sensing data. Ecosystem carbon pools, including aboveground and belowground biomass and necromass, and soil carbon to a depth of 2 m were assessed using measurements and samples from plots distributed among three canopy height classes within the park. There were two mangrove species within the inventory area in PNP, Rhizophora racemosa and R. harrisonii. R. harrisonii was predominant in the sparse, low-stature stands that dominated the west side of the park. In the east side of the park, both species occurred in tall-stature stands, with tree height often exceeding 30 m. Canopy height was an effective means to stratify the inventory area, as biomass was significantly different among the height classes. Despite those differences in aboveground biomass, the soil carbon density was not significantly different among height classes. Soils were the main component of the ecosystem carbon stock, accounting for over 84% of the total. The ecosystem carbon density ranged from 644 to 943 Mg C ha−1 among the three height classes. The ecosystem carbon stock within PNP is estimated to be 40,588 Gg C. The combination of pre-inventory information about stand conditions and their spatial distribution within the assessment area obtained from remote sensing data and a spatial decision support system were fundamental to implementing this relatively large-scale field inventory. This work exemplifies how mangrove carbon stocks can be quantified to augment national C reporting statistics, provide a baseline for projects involving monitoring, reporting and verification (i.e., MRV), and provide data on the forest composition and structure for sustainable management and conservation practices.

Carl C Trettin↗
NASA NTRS2022

A Dynamo Simulation Generating Saturn-Like Small Magnetic Dipole Tilts

Abstract Among planetary dynamos, the magnetic field of Saturn stands out in its exceptional level of axisymmetry. One of its peculiar features is that the magnetic dipole mode is tilted with respect to the planetary rotation axis by only ≈0.007° or less. Numerical dynamo simulations performed in this context have had great difficulty in producing such small dipole tilt angles without introducing ad hoc ingredients such as a latitudinally varying heat flux pattern in the outer layers or stably stratified layers (SSLs). Here we present a numerical dynamo simulation that generates a highly axisymmetric dynamo with a dipole tilt of about ≈0.0008° on average. The model consists of a deep dynamo layer and an overlying low-conductivity layer but without any SSLs. We highlight a novel mechanism where strong differential rotation generated in the atmospheric layer penetrates into the dynamo region, helping to maintain a very small magnetic dipole tilt. Plain Language Summary Saturn's dipole-dominant magnetic field exhibits a very peculiar feature: the dipole component of the planetary magnetic field is tilted by less than ≈0.007° with respect to the planetary spin axis. Numerical simulations performed in this context suggest that if a spatial heat-flux variation is imposed, along with a stably stratified region, on top of an active dynamo layer, then small dipole tilt values can be realized. Here we present a model where extremely small dipole tilt values can be achieved without these ad hoc ingredients. Our simulations demonstrate that dynamo theory allows extremely small dipole tilt values in a relatively simple model configuration.

Rakesh Kumar Yadav↗
DOE OSTI2023

High-fidelity Experiments of Turbulent Buoyant Jets from Simultaneous Particle Image Velocimetry and Laser Induced Fluorescence

Accurate models of turbulent buoyant flows are essential for the design of nuclear reactors thermal hydraulics and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We performed experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous Particle Image Velocimetry and Laser Induced Fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The Refractive Index Matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy driven flows (e.g., natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as basis for the development of a new one that captures turbulent mixing.

42 ENGINEERING↗
DOE OSTI2021

Parametric Studies of Heat Pump Water Heater Using Low GWP Refrigerants

A storage heat pump water heater (HPWH) comprises a vapor compression system and a hot water tank. It can lead to more than 70% energy saving, compared with electric resistance heating. HPWH usually uses a wrapped-tank coil as the water heater. The wrapped coil has a strong interaction with the water tank and impose a new challenge for the engineering design and simulation study. Based on the US DOE/ORNL Heat Pump Design Model (HPDM), i.e., a public-domain, component-based thermal system modelling platform, a segment-to-segment wrapped-tank coil model coupled with a one-dimensional stratified water tank model were developed, which facilitates modeling and design HPWHs having various configurations and components. Additionally, the hardware-based design model was used to evaluate refrigerants having lower global warming potentials to replace R-134a, including R-1234yf, R-1234ze, R-290, R-513A and R-450A, when effects of compressor size, refrigerant charge on 24-hour uniform energy factor (UEF) and first-hour rating (FHR) were investigated.In the drop-in replacements, R-1234ze achieved the best UEF, and R-290 resulted in the lowest UEF. When downsizing the compressor to maximize the UEFs, the UEFs of R-134a, R-290, R1234yf, R-450A and R-513A ascended to the range of 4.4 to 4.53, with R-513A having the highest UEF. The R-1234ze ended up having the lowest UEF of 4.27. Refrigerants having larger volumetric vaporization heat led to higher FHRs if using the fixed-size compressor, that R-290 resulted in the highest FHR.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Isothermal Versus Sensible Cold Storage: A Model-Based Performance Comparison for Pumped Thermal Electricity Storage

Pumped-thermal energy storage (PTES) systems consume and produce electrical energy using thermal storage media as an intermediate stage. PTES lends itself to long-duration energy storage to facilitate high penetration of intermittent electricity generation. This study presents a model-based comparison of two thermal storage types within a PTES system: a conventional, single-phase, stratified water-glycol sensible storage system (SGS), and an ideal isothermal, two-phase heat exchanger that freezes a water reservoir (isothermal heat exchanger (IHEX)). The SGS thermal storage capacity is based on the liquid’s sensible heat change with temperature, whereas the capacity of the IHEX is based on the latent heat of isothermally freezing and melting water. The idealized IHEX modeled here undergoes steady-state melting and freezing (in contrast to transient rates, as observed with ice-on-coil storage). A computational model of a complete PTES system is presented and used to evaluate the PTES system-level performance with each type of cold storage. Compared to SGS-based PTES, under nominal operating conditions, the IHEX-based PTES increased electrical round-trip efficiency from 61% to 82% and increased energy density from 1.13 to 8.09 kWh/m 3 In conclusion, the performance of the PTES configured with IHEX storage was also analyzed under varying operating parameters.

25 ENERGY STORAGE↗