Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “departure time”

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

Towards the Stable Chelation of Radioantimony(V) for Targeted Auger Theranostics

Antimony-119 ( 119 Sb) is one of the most attractive Auger-electron emitters identified to date, but it remains practically unexplored for targeted radiotherapy because no chelators have been identified to stably bind this metalloid in vivo. In a departure from current studies focused on chelator development for Sb(III), we explore the chelation chemistry of Sb(V) using the tris-catecholate ligand TREN-CAM. Through a combination of radiolabeling, spectroscopic, solid-state, and computational studies, the radiochemistry and structural chemistry of TREN-CAM with 1XX/nat Sb(V) were established. The resulting [ 1XX Sb]Sb–TREN-CAM complex remained intact for several days in human serum, signifying high stability under biological conditions. Finally, the first in vivo single photon emission computed tomography and positron emission tomography imaging studies were carried out using 117 Sb, the diagnostic analogue of 119 Sb. These studies revealed marked differences in the uptake and distribution of activity in mice administered unchelated [ 117 Sb]Sb(OH) 6 – versus [ 117 Sb]Sb–TREN-CAM, suggesting that 117 Sb is largely retained by TREN-CAM over the time course of the study. Collectively, these findings demonstrate the most physiologically stable complex of no-carrier-added 1XX Sb yet reported, offering new promise for the clinical implementation of radioantimony in nuclear medicine. In conclusion, our results also establish the feasibility of 117 Sb as an elementally matched partner to 119 Sb for theranostic applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand while simultaneously maintaining air travel as one of the safest forms of transportation. One of the reasons for this success is the ability of the air traffic control system and the operators to adapt and accommodate to situations that routinely disrupt normal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators’ ability to control. These factors can lead to states where automation is unable to properly handle these issues, and therefore air traffic controllers and pilots have to intervene — ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions, complexity increases. This is because, under these conditions humans are required to make tactical decisions in response to external factors. This results in a departure from the original strategic plan where operations would be more efficiently managed. Human operators manage airspace complexity under rigid regulations but in a constantly changing environment. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. Some prior studies devised airspace complexity metrics in commercial aviation and related these metrics to controller workload. The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic — including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to ours that identifies such contributing factors or precursor patterns.

Precursor↗

Predictions of Line Chilldown Boiling Regime Transitions by a Coupled CFD-Sub-Grid Boiling Model Validated against 1G LN2 Experiments

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to capture the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results shows good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements for several different LN2 flowrates in the vertical orientation. The strengths and weaknesses of this modeling approach are presented and discussed.

Cryogenic Fluid Management↗

Advanced Fuels Campaign Time-at-Temperature Irradiation Campaign and Post-Irradiation Examination Plan

With increasing per capita energy consumption and rapidly growing energy consumption predicted for data centers, the US faces challenges in meeting projected energy demands. The current administration is attempting to accelerate the deployment of new nuclear construction by providing significant investments in advanced nuclear and by modernizing regulatory approaches for the licensing new reactors. However, enabling extended power uprates (EPUs) for the current nuclear reactor fleet would provide a shorter-term solution to better calibrate near-term energy generation capacity to the ever-growing energy demands of the modern era. One pathway to achieve power uprates is to increase the operational window by reassessing fuel safety limits around anticipated transients. Existing light-water reactors (LWRs) utilize a targeted operational window with well-defined operational efficiencies, yet this window is also often bound by conservative safety limits that result in suboptimal operational performance. For example, the operation of the existing fleet is highly constrained by safeguards to operation related to anticipated operational occurrences (AOOs), which are moderate-frequency transients expected with a frequency greater than 0.01 per reactor-year. These AOOs include temperature excursions where the critical heat flux for the system is exceeded, resulting in a departure from nucleate boiling (DNB) at the cladding surface. This DNB event is a thermohydraulic condition that, with current conservative fuel safety limits, results in the fuel rods being classified as “failed,” meaning that the cladding may not be returned to service.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Mission Assessment Post Processor (MAPP): A New Tool for Performance Evaluation of Human Lunar Missions

The National Aeronautics and Space Administration s (NASA) Constellation Program paves the way for a series of lunar missions leading to a sustained human presence on the Moon. The proposed mission design includes an Earth Departure Stage (EDS), a Crew Exploration Vehicle (Orion) and a lunar lander (Altair) which support the transfer to and from the lunar surface. This report addresses the design, development and implementation of a new mission scan tool called the Mission Assessment Post Processor (MAPP) and its use to provide insight into the integrated (i.e., EDS, Orion, and Altair based) mission cost as a function of various mission parameters and constraints. The Constellation architecture calls for semiannual launches to the Moon and will support a number of missions, beginning with 7-day sortie missions, culminating in a lunar outpost at a specified location. The operational lifetime of the Constellation Program can cover a period of decades over which the Earth-Moon geometry (particularly, the lunar inclination) will go through a complete cycle (i.e., the lunar nodal cycle lasting 18.6 years). This geometry variation, along with other parameters such as flight time, landing site location, and mission related constraints, affect the outbound (Earth to Moon) and inbound (Moon to Earth) translational performance cost. The mission designer must determine the ability of the vehicles to perform lunar missions as a function of this complex set of interdependent parameters. Trade-offs among these parameters provide essential insights for properly assessing the ability of a mission architecture to meet desired goals and objectives. These trades also aid in determining the overall usable propellant required for supporting nominal and off-nominal missions over the entire operational lifetime of the program, thus they support vehicle sizing.

Williams, Jacob↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Validation of CFD Model Prediction of Flow Boiling Regime Transitions during LN2 Line Chilldown

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to understand the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model developed by Craft Tech that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The model constants are tested for different operating conditions and designated values are reported. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results showing good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements are presented and discussed.

Boiling regimes↗

Small-scale clumping at recombination and the Hubble tension

Despite the success of the standard ΛCDM model of cosmology, recent data improvements have made tensions emerge between low- and high-redshift observables, most importantly in determinations of the Hubble constant H 0 and the (rescaled) clustering amplitude S 8 . The high-redshift data, from the cosmic microwave background (CMB), crucially relies on recombination physics for its interpretation. In this work, we study how small-scale baryon inhomogeneities (i.e., clumping) can affect recombination and consider whether they can relieve both the H 0 and S 8 tensions. Such small-scale clumping, which may be caused by primordial magnetic fields or baryon isocurvature below kpc scales, enhances the recombination rate even when averaged over larger scales, shifting recombination to earlier times. We introduce a flexible clumping model, parametrized via three spatial zones with free densities and volume fractions, and use it to study the impact of clumping on CMB observables. We find that increasing H 0 decreases both Ω m and S 8 , which alleviates the S 8 tension. On the other hand, the shift in Ω m is disfavored by the low-z baryon-acoustic-oscillations measurements. We find that the clumping parameters that can change the CMB sound horizon enough to explain the H 0 tension also alter the damping tail, so they are disfavored by current Planck 2018 data. We test how the CMB damping-tail information rules out changes to recombination by first removing ℓ > 1000 multipoles in Planck data, where we find that clumping could resolve the H 0 tension. Furthermore, we make predictions for future CMB experiments, as their improved damping-tail precision can better constrain departures from standard recombination. Both the Simons Observatory and CMB-S4 will provide decisive evidence for or against clumping as a resolution to the H 0 tension.

79 ASTRONOMY AND ASTROPHYSICS↗

Polarimetry and Flux Distribution in the Debris Disk Around HD 32297

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at > or = 5(sigma) levels from approx. 0.45" to approx.1.7" (50-192AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of approx. 0.75" (NE side) and approx. 0.65" (SW side). Global forward-modeling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from approx. 0.25-1.6", although the central region is quite noisy and high S/N are only found in the range approx. 0.75-1.2". The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from approx. 10% at 0.55" to approx. 25% at 1.6". The maximum is found at scattering angles of 90, either from the main components of the disk or from dust grains blown out to larger radii.

high-contrast angular differential imaging (ADI)↗

Physics of greenhouse effect and convection in warm oceans

Sea surface temperature (SST) in roughly 50% of the tropical Pacific Ocean is warm enough (SST greater than 300 K) to permit deep convection. This paper examines the effects of deep convection on the climatological mean vertical distributions of water vapor and its greenhouse effect over such warm oceans. The study, which uses a combination of satellite radiation budget observations, atmospheric soundings deployed from ships, and radiation model calculations, also examines the link between SST, vertical distribution of water vapor, and its greenhouse effect in the tropical oceans. Since the focus of the study is on the radiative effects of water vapor, the radiation model calculations do not include the effects of clouds. The data are grouped into nonconvective and convective categories using SST as an index for convective activity. On average, convective regions are more humid, trap significantly more longwave radiation, and emit more radiation to the sea surface. The greenhouse effect in regions of convection operates as per classical ideas, that is, as the SST increases, the atmosphere traps the excess longwave energy emitted by the surface and reradiates it locally back to the ocean surface. The important departure from the classical picture is that the net (up minus down) fluxes at the surface and at the top of the atmosphere decrease with an increase in SST; that is, the surface and the surface-troposphere column lose the ability to radiate the excess energy to space. The cause of this super greenhouse effect at the surface is the rapid increase in the lower-troposphere humidity with SST; that of the column is due to a combination of increase in humidity in the entire column and increase in the lapse rate within the lower troposphere. The increase in the vertical distribution of humidity far exceeds that which can be attributed to the temperature dependence of saturation vapor pressure; that is, the tropospheric relative humidity is larger in convective regions. The positive coupling between SST and the radiative warming of the surface by the water vapor greenhouse effect is also shown to exist on interannual time scales.

Inamdar, A. K.↗

The OSIRIS-Rex Asteroid Sample Return: Mission Operations Design

The OSIRIS-REx mission employs a methodical, phased approach to ensure success in meeting the missions science requirements. OSIRIS-REx launches in September 2016, with a backup launch period occurring one year later. Sampling occurs in 2019. The departure burn from Bennu occurs in March 2021. On September 24, 2023, the SRC lands at the Utah Test and Training Range (UTTR). Stardust heritage procedures are followed to transport the SRC to Johnson Space Center, where the samples are removed and delivered to the OSIRIS-REx curation facility. After a six-month preliminary examination period the mission will produce a catalog of the returned sample, allowing the worldwide community to request samples for detailed analysis.Traveling and returning a sample from an Asteroid that has not been explored before requires unique operations consideration. The Design Reference Mission (DRM) ties together space craft, instrument and operations scenarios. The project implemented lessons learned from other small body missions: APLNEAR, JPLDAWN and ESARosetta. The key lesson learned was expected the unexpected and implement planning tools early in the lifecycle. In preparation to PDR, the project changed the asteroid arrival date, to arrive one year earlier and provided additional time margin. STK is used for Mission Design and STKScheduler for instrument coverage analysis.

The Evolution of Tropical Precipitation Patterns During ENSO Events Using 21+ Years of GPCP Merged Data

The ENSO phenomenon is characterized by fluctuations in the climate system of the tropical Pacific. Quantifying changes in the precipitation component of this system is important in understanding the distribution of heating in the atmosphere which drives the large-scale circulation and affects the weather patterns in the mid-latitudes. Monitoring precipitation anomalies in the Pacific is also an important component for tracking the evolution of ENSO. The most timely and complete observations of the earth come from satellite instruments. In this study, the state of the art satellite-gauge merged monthly precipitation data set from the Global Precipitation Climatology Project (GPCP) is used to depict tropical rainfall patterns during ENSO events over the past two decades and quantify these patterns using indices. This analysis will be complemented by daily precipitation data which can resolve the Madden-Julian Oscillation and westerly wind burst events. The 1997-98 El Nino and 1998-2000 La Nina were the best observed ENSO cycle in the historic record. Prior to the El Nino (in terms of anomalous warming of the east Pacific) dry anomalies over the Maritime Continent were observed in February 1997 as a westerly wind burst advected convection to the east. The largest SST anomalies occurred around November-December 1997, which were followed by the largest precipitation anomalies in the beginning of 1998. The largest precipitation departures from normal were not colocated with the SST anomalies, but were further west, In the spring of 1998 negative precipitation anomalies to the north of the equator intensified, signaling the mature phase of the El Nino. A rapid increase in the precipitation-based La Nina index from December-January 1998 to March-April 1998 signaled the coming La Nina. The 1982-1983 El Nino was comparable in strength (according to several indices) and the precipitation patterns evolved in a similar fashion. For the 1998-2000 La Nina, the coldest anomalies, were confined to the central equatorial Pacific, while the driest anomalies were found in the west Pacific,

Curtis, Scott↗

Droplet Growth or Evaporation Does Not Buffer the Variability in Supersaturation in Clean Clouds

Abstract Water vapor supersaturation in clouds is a random variable that drives activation and growth of cloud droplets. The Pi Convection–Cloud Chamber generates a turbulent cloud with a microphysical steady state that can be varied from clean to polluted by adjusting the aerosol injection rate. The supersaturation distribution and its moments, e.g., mean and variance, are investigated for varying cloud microphysical conditions. High-speed and collocated Eulerian measurements of temperature and water vapor concentration are combined to obtain the temporally resolved supersaturation distribution. This allows quantification of the contributions of variances and covariances between water vapor and temperature. Results are consistent with expectations for a convection chamber, with strong correlation between water vapor and temperature; departures from ideal behavior can be explained as resulting from dry regions on the warm boundary, analogous to entrainment. The saturation ratio distribution is measured under conditions that show monotonic increase of liquid water content and decrease of mean droplet diameter with increasing aerosol injection rate. The change in liquid water content is proportional to the change in water vapor concentration between no-cloud and cloudy conditions. Variability in the supersaturation remains even after cloud droplets are formed, and no significant buffering is observed. Results are interpreted in terms of a cloud microphysical Damköhler number (Da), under conditions corresponding to , i.e., the slow-microphysics regime. This implies that clouds with very clean regions, such that is satisfied, will experience supersaturation fluctuations without them being buffered by cloud droplet growth. Significance Statement The saturation ratio (humidity) in clouds controls the growth rate and formation of cloud droplets. When air in a turbulent cloud mixes, the humidity varies in space and time throughout the cloud. This is important because it means cloud droplets experience different growth histories, thereby resulting in broader size distributions. It is often assumed that growth and evaporation of cloud droplets buffers out some of the humidity variations. Measuring these variations has been difficult, especially in the field. The purpose of this study is to measure the saturation ratio distribution in clouds with a range of conditions. We measure the in-cloud saturation ratio using a convection cloud chamber with clean to polluted cloud properties. We found in clouds with low concentrations of droplets that the variations in the saturation ratio are not suppressed.

54 ENVIRONMENTAL SCIENCES↗

Observations of Infrared Radiative Cooling in the Thermosphere on Daily to Multiyear Timescales from the TIMED/SABER Instrument

We present observations of the infrared radiative cooling by carbon dioxide (CO2) and nitric oxide (NO) in Earth s thermosphere. These data have been taken over a period of 7 years by the SABER instrument on the NASA TIMED satellite and are the dominant radiative cooling mechanisms for the thermosphere. From the SABER observations we derive vertical profiles of radiative cooling rates (W/cu m), radiative fluxes (W/sq m), and radiated power (W). In the period from January 2002 through January 2009 we observe a large decrease in the cooling rates, fluxes, and power consistent with the declining phase of solar cycle. The power radiated by NO during 2008 when the Sun exhibited few sunspots was nearly one order of magnitude smaller than the peak power observed shortly after the mission began. Substantial short-term variability in the infrared emissions is also observed throughout the entire mission duration. Radiative cooling rates and radiative fluxes from NO exhibit fundamentally different latitude dependence than do those from CO2, with the NO fluxes and cooling rates being largest at high latitudes and polar regions. The cooling rates are shown to be derived relatively independent of the collisional and radiative processes that drive the departure from local thermodynamic equilibrium (LTE) in the CO2 15 m and the NO 5.3 m vibration-rotation bands. The observed NO and CO2 cooling rates have been compiled into a separate dataset and represent a climate data record that is available for use in assessments of radiative cooling in upper atmosphere general circulation models.

Mlynczak, Martin G.↗

Influences of the El Nino Southern Oscillation and the Pacific Decadal Oscillation on the Timing of the North American Spring

Detrended, modelled first leaf dates for 856 sites across North America for the period 1900-2008 are used to examine how the El Nino Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) separately and together might influence the timing of spring. Although spring (mean March through April) ENSO and PDO signals are apparent in first leaf dates, the signals are not statistically significant (at a 95% confidence level (p <0.05)) for most sites. The most significant ENSO/PDO signal in first leaf dates occurs for El Nino and positive PDO conditions. An analysis of the spatial distributions of first leaf dates for separate and combined ENSO/PDO conditions features a northwest-southeast dipole that is significantly (at p <0.05) different than the distributions for neutral conditions. The nature of the teleconnection between Pacific SST's and first leaf dates is evident in comparable composites for detrended sea level pressure (SLP) in the spring months. During positive ENSO/PDO, there is an anomalous flow of warm air from the southwestern US into the northwestern US and an anomalous northeasterly flow of cold air from polar regions into the eastern and southeastern US. These flow patterns are reversed during negative ENSO/PDO. Although the magnitudes of first leaf date departures are not necessarily significantly related to ENSO and PDO, the spatial patterns of departures are significantly related to ENSO and PDO. These significant relations and the long-lived persistence of SSTs provide a potential tool for forecasting the tendencies for first leaf dates to be early or late.

McCabe, Gregory J.↗

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths can be known, drift in certain instruments and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are in observation-space data formats and have not been made easily available to reanalysis users.A test data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). We refer to this proof-of-concept data as the MERRA-2 Gridded Innovations and Observations (GIO). In this paper, we present the data format and its strengths and limitations with some initial testing and validation of the methodology.

SSMI↗