Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface mass variations”

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 541 records · Page 30

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗

Supersonic Aerodynamic Characteristics of Blunt Body Trim Tab Configurations

Trim tabs are aerodynamic control surfaces that can allow an entry vehicle to meet aerodynamic performance requirements while reducing or eliminating the use of ballast mass and providing a capability to modulate the lift-to-drag ratio during entry. Force and moment data were obtained on 38 unique, blunt body trim tab configurations in the NASA Langley Research Center Unitary Plan Wind Tunnel. The data were used to parametrically assess the supersonic aerodynamic performance of trim tabs and to understand the influence of tab area, cant angle, and aspect ratio. Across the range of conditions tested (Mach numbers of 2.5, 3.5, and 4.5; angles of attack from -4deg to +20deg; angles of sideslip from 0deg to +8deg), the effects of varying tab area and tab cant angle were found to be much more significant than effects from varying tab aspect ratio. Aerodynamic characteristics exhibited variation with Mach number and forebody geometry over the range of conditions tested. Overall, the results demonstrate that trim tabs are a viable approach to satisfy aerodynamic performance requirements of blunt body entry vehicles with minimal ballast mass. For a 70deg sphere-cone, a tab with 3% area of the forebody and canted approximately 35deg with no ballast mass was found to give the same trim aerodynamics as a baseline model with ballast mass that was 5% of the total entry mass.

Korzun, Ashley M.↗

Measurement of Damage Profiles from Solar Wind Implantation

NASA's Genesis Mission launched from Cape Canaveral in August of 2001 with the goal of collecting solar wind in ultra-pure materials. The samples were returned to Earth more than three years later for subsequent analysis. Although the solar wind is comprised primarily of protons, it also contains ionized species representing the entire periodic table. The Genesis mission took advantage of the natural momentum of these ionized species to implant themselves in specialized collectors including single crystal Si and SiC. The collectors trapped the solar wind species of interest and sustained significant damage to the surface crystal structure as a result of the ion bombardment. In this work, spectroscopic ellipsometry has been used to evaluate the extent of this damage in Si and SiC samples. These results and models are compared for artificially implanted samples and pristine non-flight material. In addition, the flown samples had accumulated a thin film of molecular contamination as a result of outgassing in flight, and we demonstrate that this layer can be differentiated from the material damage. In addition to collecting bulk solar wind samples (continuous exposure), the Genesis mission actually returned silicon exposed to four different solar wind regimes: bulk, high speed, low speed, and coronal mass ejections. Each of these solar wind regimes varies in energy, but may vary in composition as well. While determining the composition is a primary goal of the mission, we are also interested in the variation in depth and extent of the damage layer as a function of solar wind regime. Here, we examine flight Si from the bulk solar wind regime and compare the results to both pristine and artificially implanted Si. Finally, there were four samples which were mounted in an electrostatic "concentrator" designed to reject a large fraction (>85%) of incoming protons while enhancing the concentration of ions mass 4-28 amu by a factor of at least 20. Two of these samples were single crystal 6H silicon carbide. (The others were polycrystalline CVD diamond and amorphous carbon that were not examined in the work.) The ion damaged SiC samples from the concentrator were studied in comparison to the flight Si from the bulk array to understand differences in the extent of the damage.

McNamara, K. M.↗

Apollo 17 Soil Characterization for Reflectance Spectroscopy

It is the fine fractions that dominate the observed spectral signatures of bulk lunar soil, and the next to the smallest size fractions are the most similar to the overall properties of the bulk soil. Thus, our Lunar Soil Characterization Consortium has concentrated on understanding the inter-relations of compositional, mineralogical, and optical properties of the <45-micron size fraction and its component sizes (20-44 micron, 10-20 micron, and <10 micron size fractions). To be able to generalize our results beyond the particular sample set studied, it is necessary to quantitatively identify the observed effects of space weathering and evaluate the processes involved. For this, it is necessary to know the chemistry of each size fraction, modal abundances of each phase, average compositions of the minerals and glasses, I(sub s)/FeO values, reflectance spectra, and the physical makeup of the individual particles and their patinas. This characterization includes the important dissection of the pyroxene minerals into four separate populations, with data on both modes and average chemical compositions. Armed with such data, it should be possible to effectively isolate spectral effects of space weathering from spectral properties related to mineral and glass chemistry. Four mare soils from the Apollo 17 site were selected for characterization based upon similarities in bulk composition and their contrasting maturities, ranging from immature to submature to mature. The methodology of our characterization has been discussed previously. Results of the Apollo 17 mare soils, outlined herein, are being prepared for publication in MAPS. As shown, with decreasing grain size, the agglutinitic (impact) glass content profoundly increases. This is the most impressive change for the mare soils. In several soils we have examined, there is an over two-fold increase in the agglutinitic glass contents between the 90-150- micron and the 10-20-micron size fractions. Accompanying this increase in agglutinitic glass is a definite decrease in pyroxenes and to lesser extents, the oxides (ilmenite), volcanic glass, and olivine. Unexpectedly, however, the absolute plagioclase abundances stay relatively constant throughout the different grain sizes, although the abundance of plagioclase relative to the mafic minerals increases with decreasing particle size. These soils were chosen for study based upon their similarities in FeO and Ti02 content, allowing for direct comparisons between evolutions of chemistry between size fractions and among different maturities of soils. The bulk chemistry of these fractions was determined by EMP analyses of fused glass beads. In contrast to the systematic variations in bulk chemistry discussed below, the relatively uniform composition of agglutinitic glass with grain size and soil maturity is illustrated. The composition of the bulk fraction of each size fraction becomes more feldspathic with increasing maturity, with the effect being most pronounced for the finest fractions. The composition of the agglutinitic glass, however, is relatively invariant and more feldspathic (i.e., rich in Al2O3) than even the <10-micron fraction. This relation not only strengthens the "fusion of the finest fraction" (F(sup 3)) hypothesis, but also highlights the important role of plagioclase in the formation of agglutinitic glass. With decreasing grain size, FeO, MgO, and TiO2 contents decrease, whereas CaO, Na2O, and Al2O3 (plag components) increase for all soils. These chemical variations would appear to be coupled with the significant increase in agglutinitic glass and decrease in oxide (ilmenite),pyroxene, and volcanic glass. These changes in chemistry do not appear to be due to distinct changes in the compositions of individual phases but to their abundances. Values of I(sub s)/FeO increase with decreasing grain size, even though the bulk FeO contents decrease. That is, the percentage of the total Fe that is present as nanophase Fe(sup O) has increased substantially in the smaller size fraction. Note that the increase in nanophase FeO in smaller size fractions is significantly greater than the increase in agglutinitic glass content, with its single-domain FeO component. This would seem to indicate that at least some of the FeO is surface correlated. To illustrate this effect, if it is assumed that the nanophase FeO is entirely surface correlated, then equal masses of 15-micron and 6-micron spheres should have about 3x as much FeO in the finer fraction. The recent findings of Kelleret al. of the major role of vapor-deposited, nanophase FeO-containing patinas on most soil particles is a major breakthrough in our understanding of the distribution of FeO within agglutinitic glass and upon grain surfaces. Bidirectional reflectance spectra for a representative Apollo 17 soil (70181) are shown. The size separates all have similar albedo in the blue and follow a regular sequence in which the continuum slope increases, ferrous bands weaken, and albedo, increases with decreasing particle size. The bulk <45-micron soil is typically close to the 10-20 micron spectrum. It is important to note that although the finest fraction (<10 micron) is close in composition to the abundant agglutinitic glass in each size fraction, this size fraction is relatively featureless and does not dominate the spectrum of the bulk <45-micron soil. It has long been suspected that agglutinitic glass, to a large extent, is the product of melting of the finest fraction of the soils, with a dominance of plagioclase. Given the low abundance of pyroxene in the finest fractions of each soil the source of the FeO in these Apollo 17 agglutinitic glasses is not fully identified. We suspect the abundant volcanic glass in these samples may be a significant contributor and this hypothesis will be tested with the suite under study from other Apollo sites.

Taylor, L. A.↗

A Validated All-Pressure Fluid Drop Model and Lewis Number Effects for a Binary Mixture

The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Harstad, K.↗

Weak-signal Phase Calibration Strategies for Large DSN Arrays

The NASA Deep Space Network (DSN) is studying arrays of large numbers of small, mass-produced radio antennas as a cost-effective way to increase downlink sensitivity and data rates for future missions. An important issue for the operation of large arrays is the accuracy with which signals from hundreds of small antennas can be combined. This is particularly true at Ka band (32 GHz) where atmospheric phase variations can be large and rapidly changing. A number of algorithms exist to correct the phases of signals from individual antennas in the case where a spacecraft signal provides a useful signal-to-noise ratio (SNR) on time scales shorter than the atmospheric coherence time. However, for very weak spacecraft signals it will be necessary to rely on background natural radio sources to maintain array phasing. Very weak signals could result from a spacecraft emergency or by design, such as direct-to-Earth data transmissions from distant planetary atmospheric or surface probes using only low gain antennas. This paper considers the parameter space where external real-time phase calibration will be necessary, and what this requires in terms of array configuration and signal processing. The inherent limitations of this technique are also discussed.

interferometry↗

Preliminary Measurements on the BOLT Geometry in the Supersonic Low Disturbance Tunnel

Experiments were performed in the Mach 3.5 Supersonic Low Disturbance Tunnel on the BOundary-Layer Transition (BOLT) geometry. The goal of this campaign was to assess changes to the transition front by varying the freestream noise and the model surface quality. The model was printed of polycarbonate and is 30% scale of the flight geometry. It was tested under noisy and quiet conditions at different streamwise positions in the tunnel and was also tested before and after improvements were made to the model surface through sanding and gap reduction. Pitch and yaw angles were nominally zero. The model surface temperature was measured using infrared thermography and thermocouples as the unit Reynolds number was swept from 4.75 – 15.8 × 10 6 m −1 . Boundary-layer transition was observed near the midspan of the model for unit Reynolds numbers above 10×10 6 m −1 in quiet flow, while transition was already observed on the shoulders below 5×10 6 m −1 in noisy flow. The model surface enhancements improved the symmetry of the transition front for quiet flow but made only slight differences for noisy flow. In quiet flow, the impact of streamwise positioning on transition was more significant on the model shoulders compared to the central region. The difference in the shoulder heating was likely due to variation in noise radiated from the nozzle sidewalls. The total temperature and initial wall temperatures were near 300 K and 293 K, respectively, which generally resulted in negative convective heat-flux values into the model. Hot-wire anemometry was also utilized to characterize the freestream and to conduct a planar survey near the base of the model for the 7.92 × 10 6 m −1 quiet condition. The survey captured the thickening of the boundary layer at the centerline and the vortical nature of the flow outboard of it. Spectral analysis of the mass-flux fluctuations near the midspan demonstrates that the boundary layer is laminar and suggests that the thick boundary layer caused the relatively warm centerline observed in the infrared images.

boundary layer transition↗

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Variable Emissivity Through MEMS Technology

This paper discusses a new technology for variable emissivity (vari-e) radiator surfaces, which has significant advantages over traditional radiators and promises an alternative design technique for future spacecraft thermal control systems. All spacecraft rely on radiative surfaces to dissipate waste heat. These radiators have special coatings, typically with a low solar absorptivity and a high infrared-red emissivity, that are intended to optimize performance under the expected heat load and thermal sink environment. The dynamics of the heat loads and thermal environment make it a challenge to properly size the radiator and often require some means of regulating the heat rejection rate of the radiators in order to achieve proper thermal balance. Specialized thermal control coatings, which can passively or actively adjust their emissivity offer an attractive solution to these design challenges. Such systems would allow intelligent control of the rate of heat loss from a radiator in response to heat load and thermal environmental variations. Intelligent thermal control through variable emissivity systems is well suited for nano and pico spacecraft applications where large thermal fluctuations are expected due to the small thermal mass and limited electric resources. Presently there are three different types of vari-e technologies under development: Micro ElectroMechanical Systems (MEMS) louvers, Electrochromic devices, and Electrophoretic devices. This paper will describe several prototypes of micromachined (MEMS) louvers and experimental results for the emissivity variations measured on theses prototypes. It will further discuss possible actuation mechanisms and space reliability aspects for different designs. Finally, for comparison parametric evaluations of the thermal performances of the new vari-e technology and standard thermal control systems are presented in this paper.

Darrin, Ann Garrison↗

Preliminary Measurements on the BOLT Geometry in the Supersonic Low Disturbance Tunnel

Experiments were performed in the Mach 3.5 Supersonic Low Disturbance Tunnel on the BOundary-Layer Transition (BOLT) geometry. The goal of this campaign was to assess changes to the transition front by varying the freestream noise and the model surface quality. The model was printed of polycarbonate and is 30% scale of the flight geometry. It was tested under noisy and quiet conditions at different streamwise positions in the tunnel and was also tested before and after improvements were made to the model surface through sanding and gap reduction. Pitch and yaw angles were nominally zero. The model surface temperature was measured using infrared thermography and thermocouples as the unit Reynolds number was swept from 4.75 -- 15.8 x 10^6 m^{-1}. Boundary-layer transition was observed near the midspan of the model for unit Reynolds numbers above 10 x 10^6 m^{-1} in quiet flow, while transition was already observed on the shoulders below 5 x 10^6 m^{-1} in noisy flow. The model surface enhancements improved the symmetry of the transition front for quiet flow but made only slight differences for noisy flow. In quiet flow, the impact of streamwise positioning on transition was more significant on the model shoulders compared to the central region. The difference in the shoulder heating was likely due to variation in noise radiated from the nozzle sidewalls. The total temperature and initial wall temperatures were near 300 K and 293 K, respectively, which generally resulted in negative convective heat-flux values into the model. Hot-wire anemometry was also utilized to characterize the freestream and to conduct a planar survey near the base of the model for the 7.92 x 10^6 m^{-1} quiet condition. The survey captured the thickening of the boundary layer at the centerline and the vortical nature of the flow outboard of it. Spectral analysis of the mass-flux fluctuations near the midspan demonstrates that the boundary layer is laminar and suggests that the thick boundary layer caused the relatively warm centerline observed in the infrared images.

boundary layer transition↗

Updated Human Mars Ascent Vehicle Concept in Support of NASA’s Strategic Analysis Cycle 2021

The NASA Artemis program has brought significant change to the agency’s human exploration strategy over the last several years. The human exploration of Mars remains an ultimate objective of the overall human exploration strategy. However, how the agency intends to execute initial human exploration of Mars has shifted to better align with current policy. These changes find their way into the design of key architecture elements, such as the Mars Ascent Vehicle (MAV). Of the numerous changes in the conceptual human Mars architecture, two in particular have had significant impact on the design of the human MAV: the desire for minimal surface infrastructure for initial human Mars missions, and technology investment timelines in support of initial human missions to Mars. The first leads to a surface mission architecture targeting a 30-day surface duration with two crewmembers. The second drives out in-situ resource utilization (ISRU) from the initial sortie. As a result, the reference MAV design for initial human Mars missions in the Human Exploration and Operations Mission Directorate’s Strategic Analysis Cycle 2021 (SAC21) is based around supporting the return of 2 crew from the surface of Mars, without surface ISRU. This has led to several design changes since 2019, when details of a reference human MAV were last published. The MAV concept for SAC21 consists of a two-stage vehicle capable of supporting two crew for 84 hours. This duration is derived from a baseline ascent trajectory targeting an apoapsis rendezvous with the habitation element in a nominal 5-sol Mars aggregation orbit. With the goal of minimizing technology investments for the initial surface sortie removing ISRU from the trade space, the MAV must now be landed either fully loaded with propellant, or with some portion of the propellant off-loaded. Landing a partially-fueled MAV will require additional surface assets to perform robotic propellant transfer on the surface of Mars from a pre-emplaced propellant depot. Furthermore, cryogenic fluid management and storage for extended duration with no losses represents additional technology investment that would have to be made to support cryogenic-based propellant concepts that have been the baseline in past architectures. As a result, the SAC21 architecture utilizes a storable propellant-based MAV propulsion system to minimize these potential technology investments incurred by cryogenic-based propulsion systems. Recent efforts also focused on evaluating the sensitivity of key design parameters, such as landing site latitude, elevation, and local atmospheric conditions. All play a large role in determining the predicted propellant requirements of the MAV. Understanding their potential impacts is important because of the rippling effects of changes to the MAV design on other architecture elements, such as entry, descent, landing, and transportation systems. Results indicate up to +2% to -4% wet mass variation due to landing site latitude, with higher latitudes greater than 70 degrees North resulting in greater than 6% wet mass increases. The MAV saw reduced impact due to atmospheric design parameters, which only accounting for -0.6% to +0.4% wet mass variation.

Douglas J. Trent↗

Surface conversion techniques for low energy neutral atom imagers

This investigation has focused on development of key technology elements for low energy neutral atom imaging. More specifically, we have investigated the conversion of low energy neutral atoms to negatively charged ions upon reflection from specially prepared surfaces. This 'surface conversion' technique appears to offer a unique capability of detecting, and thus imaging, neutral atoms at energies of 0.01 - 1 keV with high enough efficiencies to make practical its application to low energy neutral atom imaging in space. Such imaging offers the opportunity to obtain the first instantaneous global maps of macroscopic plasma features and their temporal variation. Through previous in situ plasma measurements, we have a statistical picture of large scale morphology and local measurements of dynamic processes. However, with in situ techniques it is impossible to characterize or understand many of the global plasma transport and energization processes. A series of global plasma images would greatly advance our understanding of these processes and would provide the context for interpreting previous and future in situ measurements. Fast neutral atoms, created from ions that are neutralized in collisions with exospheric neutrals, offer the means for remotely imaging plasma populations. Energy and mass analysis of these neutrals provides critical information about the source plasma distribution. The flux of neutral atoms available for imaging depends upon a convolution of the ambient plasma distribution with the charge exchange cross section for the background neutral population. Some of the highest signals are at relatively low energies (well below 1 keV). This energy range also includes some of the most important plasma populations to be imaged, for example the base of the cleft ion fountain.

Quinn, J. M.↗

The Potassium Isotope Composition of Aggregate Material From Asteroid Bennu

NASA’s OSIRIS-REx mission re-turned material from asteroid Bennu on September 24, 2023, marking the first time a U.S. mission has delivered asteroid samples to Earth, and the largest asteroid sample return to date. Due to the pristine nature of these primitive carbonaceous-rich samples, they allow us a rare opportunity to study our solar system’s formation to a degree not previously possible. Isotopes of moderately volatile elements (MVEs) have been recently developed as robust tracers for tracking different volatilization events within our Solar System. Of the MVEs, K has gained significant interest in recent years due to its ideal chemical and physical properties combined with technique improvements. Due to this, several recent studies have investigated the K isotope systematics within bulk chondrites and found an isotopic dichotomy be-tween the carbonaceous and non-carbonaceous chondrites (Figure 1), likely reflecting their different reservoirs which formed in the inner and outer protoplanetary disk respectively. Furthermore, systematic K isotope variations across the different carbonaceous chondrite groups and correlations with mass independent isotope systems such as Cr, Ti, and Ni, have also been observed. As a result, establishing the K isotope composition of bulk Bennu can be used to test the mission hypothesis that “Bennu’s parent asteroid accreted in the outer protoplanetary disk, beyond Jupiter, as recorded by distinct isotopic anomalies in a variety of elements” . In addition to early solar system processes, K isotopes can also be fractionated by space weathering as shown by comparisons between mature and immature lunar regolith. The degree of K isotope fractionation caused by space weathering correlates with the regolith maturity index and is an order of magnitude larger than what is observed in bulk meteorite samples which have not undergone strong space weathering processes. As such, the K isotope systematics of Bennu samples could provide an independent assessment of regolith maturity of the asteroid surface, helping to directly test the mission hypothesis that “Space weathering changed the chemistry and mineralogy of optically active surfaces” . In order to test both this hypothesis, and the hypothesis relating to Bennu’s accretion location, we aim to conduct high-precision K isotope analysis on bulk Bennu aggregates.

K Wang↗

Modeling and Analysis for Spent Nuclear Fuel Seismic Testing

Note: This is the final public release version of PNNL-31671 Draft, which was previously released to the sponsor for review. The intent is the public release version will be PNNL-31671. There are no significant changes from the earlier draft version. This unlimited distribution report is the deliverable for M3SF-21PN010202014: Modeling and Analysis for Spent Nuclear Fuel Seismic Testing. This report summarizes the modeling, analysis, and test plan support completed by Pacific Northwest National Laboratory for the spent nuclear fuel dry storage system seismic test plan through May of 2021. Test plan preparation is planned to continue until the seismic test is completed in July of 2022. This report covers preliminary structural dynamic model development and computer-aided design of test hardware. Based on preliminary modeling, the strongest earthquakes under consideration in this test program provide mechanical loading on the spent nuclear fuel that is comparable to the 30 cm cask drop scenario, although the loads do not appear to be strong enough to cause significant permanent deformation of the fuel assembly spacer grids. The potential for grid deformation during the test will be assessed when the proposed shake table motion becomes available. The weakest earthquakes considered in this test are expected to be comparable to the mechanical loads witnessed in the multimodal transportation test of 2017. The horizontal canister system is predicted to provide nearly-uniform loading condition on the fuel assemblies it contains, while the vertical cask system is predicted to cause non-uniform dynamic loads on the fuel assemblies. In the horizontal case, gravity keeps the fuel assemblies in contact with one basket wall surface unless the loads are strong enough to cause a separation. In the vertical case, the fuel assemblies are relatively long and slender, and seismic motion in the anticipated test range is predicted to cause the assemblies to lean, tilt, and impact the basket walls throughout the seismic event. These gap closures are a nonlinear force transmission condition, so the vertical system is expected to have more variation and variability than the horizontal system. While the horizontal system is expected to have a relatively more linear response than the vertical system, there is still the potential for nonlinear behavior in the horizontal system because the fuel assemblies are free to slide, bounce, and impact the basket walls if the seismic loads are strong enough. One major conclusion of this study is that the use of mixed fuel assemblies in the canister will be acceptable. There are differences in overall system response when the mass, center of gravity, or gaps are changed, but the changes in response are within the bounds of a system that contains completely homogeneous fuel assemblies. One important observation is that the loading for each individual fuel assembly within the vertical canister is expected to be different from that of the others because of nonlinearities in the system. The horizontal canister system is expected to have a more uniform and more predictable response than the vertical canister, making variations in fuel assembly characteristics easier to account for. One important recommendation that comes from this modeling work is to repeat some of the strongest shake tests at different angles, particularly for the vertical cask system. Modeling predicts that strong seismic motion will cause the fuel assemblies to close initial gaps and impact the fuel basket walls when the canister is in the vertical cask configuration. Each shake test will impose a pre-defined three-dimensional time history on the cask system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fine particulate concentrations over East Asia derived from aerosols measured by the Advanced Himawari Imager using machine learning

Fine particulate matter with a diameter below 2.5 μm (PM 2.5 ) is deleterious to the cardiovascular and respiratory systems. It is often difficult to assess the effects of PM 2.5 on human health over regions with limited ground monitoring sites, especially in East Asia. As an alternative, we estimated near-surface PM 2.5 concentrations by analyzing Advanced Himawari Imager (AHI) Yonsei Aerosol Retrieval (YAER) products. This study incorporates daytime data for East Asia covering the Korean Peninsula, China, Japan, Southeast Asia, and southern Mongolia. We collocated AHI YAER product pixels with meteorological, land-cover, and other ancillary data for the period from March 2018 to February 2019. To estimate PM 2.5 concentrations over wide areas spanning many countries displaying various relationships between aerosol optical depth and PM 2.5 , monthly models were developed by considering both the spatial and temporal characteristics of ground-based PM 2.5 measurements. Random forest machine learning model estimated ground-level mass concentrations of PM 2.5 ; subsequent 10-fold cross validation (CV) yielded a CV R 2 value of 0.81 and a CV root mean squared error (RMSE) of 12.3 μg m -3 . We investigated the spatial pattern of PM 2.5 concentrations over multiple countries and seasonal variation in PM 2.5 concentrations. Diurnal variation of a severe PM 2.5 event in the Korean Peninsula was investigated as a case study. The model captured the extremely heterogeneous spatial distribution of PM 2.5 concentrations peaked around local noon. To measure the capability of the developed model to estimate PM 2.5 concentrations in areas with few in-situ data, its predictive performance was evaluated using a dataset independent of the training process with an R 2 of 0.60 and RMSE of 8.18 μg m −3 . This study demonstrates the potential for satellite-based PM 2.5 estimation for areas with insufficient measuring stations.

Pm2.5↗

Impact melt generation and transport

The results from the first two calculations in a series of continuum mechanics computer code calculations, investigating the effects of variations in impactor mass and velocity on the generation and transport of impact melt, are reported. In the present calculations, the impactor is modeled as a spherical iron projectile with a mass of one trillion grams, and the target as a gabbroic anorthosite (GA) half-space, where the cases calculated have impact velocities of 5 and 15.8 km/sec. Early-time ejection velocities are 1-2 km/sec in both cases. The first calculation results in 0.07 projectile masses of GA being partly or completely melted, with all the melted GA being ejected from the crater, and a maximum impact range for the ejected melted material of 30 km. The second calculation yields 10.4 projectile masses of melted GA, 50% of which is ejected from the crater to ranges of up to about 130 km. Peak shock pressure attenuation with depth is reported for both cases, and transient cavity dynamics are described and compared to that for surface and near-surface explosions.

Orphal, D. L.↗

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

What Controls the Vertical Distribution of Aerosol? Relationships Between Process Sensitivity in HadGEM3-UKCA and Inter-Model Variation from AeroCom Phase II

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3-UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3-UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN >3 nm), while the profiles of larger particles (e.g. CN>100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.

dust↗