Search NASA⌕ Search

SEARCH · Search NASA

Results for “lateral mixing”

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 271 records · Page 15

Arctic chlorine monoxide observations during spring 1993 over Thule, Greenland, and implications for ozone depletion

We have determined the vertical distribution of chlorine monoxide (ClO), from measurements of pressure-broadened molecular-emission spectra made over Thule, Greenland, during the 1993 Arctic spring. The measurements show a weak lower stratospheric layer of chlorine monoxide inside the vortex in late February, which was, however, significantly greater in mixing ratio than that seen in observations we made in the spring of 1992. ClO was also observed in much smaller quantities in early to mid-March 1993 when Thule was outside the vortex. The amount of ClO within the vortex was severely reduced by the time it returned over Thule in late March. This reduction occurred several weeks earlier relative to the winter solstice than the decline of ClO inside the Antarctic vortex in 1993. The enhanced Arctic lower stratospheric layer seen in late February 1993 at a nearly equivalent photochemical period, and beyond. We have calculated daily ozone loss rates, due primarily to the dimer chlorine catalytic cycle, from both sets of measurements. The vertical integral of the Arctic daily percentage ozone loss when the largest ClO levels were present, at the end of February, is found to be approximately one quarter of that in the Antarctic at a photochemical period only 1 week later. The relative weakness of daily ozone depletion, combined with the early disappearance of ClO in the Arctic, suggests that hemispheric dilution by ozone-poor air from within the Arctic vortex is unlikely to be sufficient to explain the historically extreme loss of midlatitude northern hemisphere ozone which began in 1992 and persisted throughout 1993.

Shindell, D. T.↗

Stratigraphy of the Apollo 17 Landslide Core 73002 From FMR Maturity and VNIR and Mössbauer Spectroscopy

A suite of samples from the Apollo missions to the Moon (1969-1973) were set aside and stored under controlled conditions to have unexamined lunar samples available decades later for analyses that take advantage of evolved sample handling techniques, maturation of then existing instrumentation, and development of new analytical techniques and instrumentation [e.g., 1]. One preserved sample is the Apollo 17 double drive tube core (73001/2) that was driven into the lunar surface on the landslide deposit at Station 3 on the South Massif in the Taurus-Littrow valley [2]. The deeper section (73001) was stored frozen in a Core Sample Vacuum Container (CSVC) to maximize preservation of lunar volatiles. We report here, as a part of the Apollo Next Generation Sample Analysis (ANGSA) Initiative [1,3], stratigraphy for the upper core section (73002) with respect to maturity (ferromagnetic resonance (FMR) maturity index Is/FeO [4-6]), visible-near-IR (VNIR) spectroscopy, and, for representative samples, grain-size analysis and Mössbauer spectroscopy. Stratigraphy provides data to model the dynamics of lunar landslide deposits and their post-emplacement evolution in comparison to impact-driven mixing on an airless body and stratigraphic context for volatile and (if any) organic behavior. Reported by [7] are preliminary results for multispectral imaging and hyperspectral scanning of the first dissection pass of 73002.

R. V. Morris↗

Magnetotail Structure and its Internal Particle Dynamics During Northward IMF

This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawnside are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Ashour-Abdalia, M.↗

Magnetotail Structure and its Internal Particle Dynamics During Northward IMF

This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawn side are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Ashour-Abdalla, M.↗

Drop Size Distribution - Based Separation of Stratiform and Convective Rain

For applications in hydrology and meteorology, it is often desirable to separate regions of stratiform and convective rain from meteorological radar observations, both from ground-based polarimetric radars and from space-based dual frequency radars. In a previous study by Bringi et al. (2009), dual frequency profiler and dual polarization radar (C-POL) observations in Darwin, Australia, had shown that stratiform and convective rain could be separated in the log10(Nw) versus Do domain, where Do is the mean volume diameter and Nw is the scaling parameter which is proportional to the ratio of water content to the mass weighted mean diameter. Note, Nw and Do are two of the main drop size distribution (DSD) parameters. In a later study, Thurai et al (2010) confirmed that both the dual-frequency profiler based stratiform-convective rain separation and the C-POL radar based separation were consistent with each other. In this paper, we test this separation method using DSD measurements from a ground based 2D video disdrometer (2DVD), along with simultaneous observations from a collocated, vertically-pointing, X-band profiling radar (XPR). The measurements were made in Huntsville, Alabama. One-minute DSDs from 2DVD are used as input to an appropriate gamma fitting procedure to determine Nw and Do. The fitted parameters - after averaging over 3-minutes - are plotted against each other and compared with a predefined separation line. An index is used to determine how far the points lie from the separation line (as described in Thurai et al. 2010). Negative index values indicate stratiform rain and positive index indicate convective rain, and, moreover, points which lie somewhat close to the separation line are considered 'mixed' or 'transition' type precipitation. The XPR observations are used to evaluate/test the 2DVD data-based classification. A 'bright-band' detection algorithm was used to classify each vertical reflectivity profile as either stratiform or convective, depending on whether or not a clearly-defined melting layer is present at an expected height, and if present, maximum reflectivity within the melting layer as well as the corresponding height are determined. We will present results of quantitative comparisons between the XPR observations-based classifications and the simultaneous 2DVD data-based classifications. Time series comparisons will be presented for thirteen events in Huntsville.

Thurai, Merhala↗

Sampling of Trace Atmospheric Constituents Above the Surface of Mars from a Montgolfiere Balloon

We suggest that use of an aerially deployed, solar-heated balloon, or Montgolfiere, during a Martian polar summer could provide an optimal platform for atmospheric measurements near the surface. Solar Montgolfieres are empty balloons with a hole in the bottom that fill with ambient atmosphere while failing and are rapidly heated in by the sun, thus providing buoyancy. Such a balloon would be driven by the Martian winds in the constant sunlight of a high latitude entry and could provide a traverse of thousands of kilometers. It would encircle the Martian north polar region over a period,likely lasting over one earth month, and its altitude could be controlled as low as 100 meters to as high as 4000 meters above the surface. The significant payload of tens of kilograms allowed by this technology could readily accommodate the mass of instrumentation and the power source needed for a sensitive search for disequilibrium species that might provide the atmospheric signature of localized subsurface biological or geological activity. A mass spectrometer used in combination with gas enrichment techniques can provide sensitivities to sub parts per billion for several species of interest. Additional highly specific and sensitive detectors can search for other species or measurements of their mixing ratio variation along the path of the balloon. It would be useful to combine these atmospheric measurements with electromagnetic sounding to search for variations in subsurface conductivity that might indicate the presence of subsurface water. After a month of operation, the balloon could also soft land the DBMS to provide additional data on soil samples from a 1-2 meter subsurface drill. In order to float a 30-kg payload at altitudes between 100 m and 4000 m above the Martian north polar surface, an aluminized polyethylene Montgolfiere would be approximately 29 meters in diameter and have a mass of 27 kg. Altitude variations could be accomplished by means of an upper vent to release hot gas, as has been tested at JPL, or be means of a simple, internal winch mechanism to decrease the volume of the balloon. JPL has successfully deployed three polyethylene Montgolfieres at 4-6 mbar pressure in the Earth's stratosphere and will be deploying a 20-m diameter polyethylene Montgolfiere later in 2001.

Jones, Jack A.↗

An atlas of mixed-morphology pairs

Current theories of galaxy formation imply that environment is the most important factor in deciding whether a collapsing region becomes a spiral galaxy or an elliptical galaxy. If this is the case, then isolated pairs comprising a galaxy of each morphology provide an especially interesting test, since, strictly speaking, they should not exist. The first phase of studying E-S pairs is to define a sample which is genuinely isolated and genuinely E-S. The catalog compiled by Karachentsev is the most useful Northern Hemisphere starting point, containing 603 pairs which obeyed certain well-defined angular isolation criteria. All of these pairs now have velocities, enabling further removal of at least some of the non-physical binaries. The conventional cutoff at 500 km s(exp -1) was used for this purpose. However, the morphological types were estimated from Palomar Sky Survey prints, and are consequently not very accurate. It was therefore decided to obtain good charge coupled device (CCD) images in the B-band of all pairs from the Karachentsev sample which were classified as either E-S (73 pairs) or E-SO (91 pairs), for a total of 164 pairs. This includes the 24 pairs with known redshift differences greater than 500 km s(exp -1) (9 E-S and 15 E-SO), which were also imaged in this survey, not only for completeness but also in case the criteria were later changed, or the redshifts were found to be in error. A total of 198 CCD frames were obtained between February 1986 and the present, with most of the observations being completed by the summer of 1987. A number of missing stragglers and replacements for poor quality images were taken more recently. Most observations used the number 1-0.9m telescope at Kitt Peak, with a few frames from the 2.1m. Some pairs required more than one image due to their separation. The morphological study did not require photometric conditions, but calibrations were made for almost all images anyway, and the weather was probably photometric at least half of the time, at least to an accuracy of around 0.1 magnitudes. The results of this first phase will be published in the form of an Atlas, from which any interested readers will be able to select their own sub-sample using their preferred morphological criteria.

Sharp, N. A.↗

Assessment of Microbiologically Influenced Corrosion Potential in the International Space Station Internal Active Thermal Control System Heat Exchanger Materials: A 6-Momths Study

The fluid in the Internal Active Thermal Control System (IATCS) of the International Space Station (ISS) is water based. The fluid in the ISS Laboratory Module and Node 1 initially contained a mix of water, phosphate (corrosion control), borate (pH buffer), and silver sulfate (Ag2SO4) (microbial control) at a pH of 9.5+/-0.5. Over time, the chemistry of the fluid changed. Fluid changes included a pH drop from 9.5 to 8.3 due to diffusion of carbon dioxide (CO2) through Teflon(reistered Trademark) (DuPont) hoses, increases in dissolved nickel (Ni) levels, deposition of silver (Ag) to metal surfaces, and precipitation of the phosphate (PO4) as nickel phosphate (NiPO4). The drop in pH and unavailability of a antimicrobial has provided an environment conducive to microbial growth. Microbial levels in the fluid have increased from >10 colony-forming units (CFUs)/100 ml to 10(exp 6) CFUs/100 ml. The heat exchangers in the IATCS loops are considered the weakest point in the loop because of the material thickness (=7 mil). It is made of a Ni-based braze filler/CRES 347. Results of a preliminary test performed at Hamilton Sundstrand indicated the possibility of pitting on this material at locations where Ag deposits were found. Later, tests have confirmed that chemical corrosion of the materials is a concern for this system. Accumulation of micro-organisms on surfaces (biofilm) can also result in material degradation and can amplify the damage caused by the chemical corrosion, known as microbiologically influenced corrosion (MIC). This paper will discuss the results of a 6-mo test performed to characterize and quantify the damage from microbial accumulation on the surface of the ISS/ATCS heat exchanger materials. The test was designed to quantify the damage to the materials under worst-case conditions with and without micro-organisms present at pH 8.3 and 9.5.

Roman, Monsi C.↗

Definition of experiments to investigate fire suppressants in microgravity

Defined and justified here are the conceptual design and operation of a critical set of experiments expected to yield information on suppressants and on suppressant delivery systems under realistic spacecraft-fire conditions (smoldering). Specific experiment parameters are provided on the solid fuel (carbon), oxidants (habitable spacecraft atmospheres), fuel/oxidant supply, mixing mode, and rate (quiescent and finite; ventilated and replenishable), ignition mode, event, and reignition tendency, fire-zone size, fire conditions, lifetime, and consequences (toxicity), suppressants (CO2, H2O, N2) and suppressant delivery systems, and diagnostics. Candidate suppressants were identified after an analysis of how reduced gravity alters combustion, and how these alterations may influence the modes, mechanisms, and capacities of terrestrial agents to suppress unwanted combustion, or fire. Preferred spacecraft suppression concepts included the local, near-quiescent application of a gas, vapor, or mist that has thermophysical fire-suppression activity and is chemically inert under terrestrial (normal gravity) combustion conditions. The scale, number, and duration (about 1 hour) of the proposed low-gravity experiments were estimated using data not only on the limitations imposed by spacecraft-carrier (Shuttle or Space Station Freedom) accommodations, but also data on the details and experience of standardized smolder-suppression experiments at normal gravity. Deliberately incorporated into the conceptual design was sufficient interchangeability for the prototype experimental package to fly either on Shuttle now or Freedom later. This flexibility is provided by the design concept of up to 25 modular fuel canisters within a containment vessel, which permits both integration into existing low-gravity in-space combustion experiments and simultaneous testing of separate experiments to conserve utilities and time.

Reuther, James J.↗

A Mixed Integer Linear Program for Airport Departure Scheduling

Aircraft departing from an airport are subject to numerous constraints while scheduling departure times. These constraints include wake-separation constraints for successive departures, miles-in-trail separation for aircraft bound for the same departure fixes, and time-window or prioritization constraints for individual flights. Besides these, emissions as well as increased fuel consumption due to inefficient scheduling need to be included. Addressing all the above constraints in a single framework while allowing for resequencing of the aircraft using runway queues is critical to the implementation of the Next Generation Air Transport System (NextGen) concepts. Prior work on airport departure scheduling has addressed some of the above. However, existing methods use pre-determined runway queues, and schedule aircraft from these departure queues. The source of such pre-determined queues is not explicit, and could potentially be a subjective controller input. Determining runway queues and scheduling within the same framework would potentially result in better scheduling. This paper presents a mixed integer linear program (MILP) for the departure-scheduling problem. The program takes as input the incoming sequence of aircraft for departure from a runway, along with their earliest departure times and an optional prioritization scheme based on time-window of departure for each aircraft. The program then assigns these aircraft to the available departure queues and schedules departure times, explicitly considering wake separation and departure fix restrictions to minimize total delay for all aircraft. The approach is generalized and can be used in a variety of situations, and allows for aircraft prioritization based on operational as well as environmental considerations. We present the MILP in the paper, along with benefits over the first-come-first-serve (FCFS) scheme for numerous randomized problems based on real-world settings. The MILP results in substantially reduced delays as compared to FCFS, and the magnitude of the savings depends on the queue and departure fix structure. The MILP assumes deterministic aircraft arrival times at the runway queues. However, due to taxi time uncertainty, aircraft might arrive either earlier or later than these deterministic times. Thus, to incorporate this uncertainty, we present a method for using the MILP with "overlap discounted rolling planning horizon". The approach is based on valuing near-term decision results more than future ones. We develop a model of taxitime uncertainty based on real-world data, and then compare the baseline FCFS delays with delays using the above MILP in a simple rolling-horizon method and in the overlap discounted scheme.

Gupta, Gautam↗

Aerogel/Particle Composites for Thermoelectric Devices

Optimizing solution chemistry and the addition of titania and fumed silica powder reduces shrinkage. These materials would serve to increase thermal efficiency by providing thermal insulation to suppress lateral heat leaks. They would also serve to prolong operational lifetime by suppressing sublimation of certain constituents of thermoelectric materials (e.g., sublimation of Sb from CoSb3) at typical high operating temperatures. [The use of pure silica aerogels as cast-in-place thermal-insulation and sublimation-suppression materials was described in "Aerogels for Thermal Insulation of Thermoelectric Devices" (NPO-40630), NASA Tech Briefs, Vol. 30, No. 7 (July 2006), page 50.] A silica aerogel is synthesized in a solgel process that includes preparation of a silica sol, gelation of the sol, and drying of the gel in a solvent at a supercritical temperature and pressure. The utility of pure silica aerogel is diminished by a tendency to shrink (and, therefore, also to crack) during the gelation and supercritical-drying stages. Moreover, to increase suppression of sublimation, it is advantageous to make an aerogel having greater density, but shrinkage and cracking tend to increase with density. A composite material of the type under investigation consists mostly of titania oxide powder particles and a small addition of fumed silica powder, which are mixed into the sol along with other ingredients prior to the gelation stage of processing. The silica aerogel and fumed silica act as a binder, gluing the titania particles together. It is believed that the addition of fumed silica stiffens the aerogel network and reduces shrinkage during the supercritical-drying stage. Minimization of shrinkage enables establishment of intimate contact between thermoelectric legs and the composite material, thereby maximizing the effectiveness of the material for thermal insulation and suppression of sublimation. To some extent, the properties of the composite can be tailored via the proportions of titania and other ingredients. In particular (see figure), the addition of a suitably large proportion of titania (e.g., 0.6 g/cu cm) along with a 10-percent increase in the amount of tetraethylorthosilicate [TEOS (an ingredient of the sol)] to an aerogel component having a density 40 mg/cm3makes it possible to cast a high-average-density (>0.1 g/cm3) aerogel/particle composite having low shrinkage (2.3 percent).

Paik, Jong-Ah↗

Interactive Ozone and Methane Chemistry in GISS-E2 Historical and Future Climate Simulations

The new generation GISS climate model includes fully interactive chemistry related to ozone in historical and future simulations, and interactive methane in future simulations. Evaluation of ozone, its tropospheric precursors, and methane shows that the model captures much of the largescale spatial structure seen in recent observations. While the model is much improved compared with the previous chemistry-climate model, especially for ozone seasonality in the stratosphere, there is still slightly too rapid stratospheric circulation, too little stratosphere-to-troposphere ozone flux in the Southern Hemisphere and an Antarctic ozone hole that is too large and persists too long. Quantitative metrics of spatial and temporal correlations with satellite datasets as well as spatial autocorrelation to examine transport and mixing are presented to document improvements in model skill and provide a benchmark for future evaluations. The difference in radiative forcing (RF) calculated using modeled tropospheric ozone versus tropospheric ozone observed by TES is only 0.016W/sq. m. Historical 20th Century simulations show a steady increase in whole atmosphere ozone RF through 1970 after which there is a decrease through 2000 due to stratospheric ozone depletion. Ozone forcing increases throughout the 21st century under RCP8.5 owing to a projected recovery of stratospheric ozone depletion and increases in methane, but decreases under RCP4.5 and 2.6 due to reductions in emissions of other ozone precursors. RF from methane is 0.05 to 0.18W/ sq. m higher in our model calculations than in the RCP RF estimates. The surface temperature response to ozone through 1970 follows the increase in forcing due to tropospheric ozone. After that time, surface temperatures decrease as ozone RF declines due to stratospheric depletion. The stratospheric ozone depletion also induces substantial changes in surface winds and the Southern Ocean circulation, which may play a role in a slightly stronger response per unit forcing during later decades. Tropical precipitation shifts south during boreal summer from 1850 to 1970, but then shifts northward from 1970 to 2000, following upper tropospheric temperature gradients more strongly than those at the surface.

Ozone↗

Expanding the REE Partitioning Database for Lunar Materials

Positive europium anomalies are ubiquitous in the plagioclase-rich rocks of the lunar highlands, and complementary negative Eu anomalies are found in most lunar basalts. This is taken as evidence of a large-scale differentation event, with crystallization of a global-scale lunar magma ocean (LMO) resulting in a plagioclase flotation crust and a mafic lunar interior from which mare basalts were later derived. However, the extent of the Eu anomaly in lunar rocks is variable. Some plagioclase grains in a lunar impact rock (60635) have been reported to display a negative Eu anomaly, or in some cases single grains display both positive and neagtive anomalies. Cathodoluminescence images reveal that some crystals have a negative anomaly in the core and positive at the rim, or vice versa, and the negative anomalies are not associated with crystal overgrowths. Oxygen fugacity is known to affect Eu partitioning into plagioclase, as under low fO2 conditions Eu can be divalent, and has an ionic radius similar to Ca2+ - significant in lunar samples where plagioclase compositions are predominantly anorthitic. However, there are very few experimental studies of rare earth element (REE) partitioning in plagioclase relevant to lunar magmatism, with only two plagioclase DEu measurements from experiments using lunar materials, and little data in low fO2 conditions relevant to the Moon. We report on REE partitioning experiments on lunar compositions. We investigate two lunar basaltic compositions, high-alumina basalt 14072 and impact melt breccia 60635. These samples span a large range of lunar surface bulk compositions. The experiments are carried out at variable fO2 in 1 bar gas mixing furnaces, and REE are analysed by and LA-ICP-MS. Our results not only greatly expand the existing plagioclase DREE database for lunar compositions, but also investigate the significance of fO2 in Eu partitioning, and in the interpretation of Eu anomalies in lunar materials.

Rapp, Jennifer F.↗

Downscaling GISS ModelE Boreal Summer Climate over Africa

The study examines the perceived added value of downscaling atmosphere-ocean global climate model simulations over Africa and adjacent oceans by a nested regional climate model. NASA/Goddard Institute for Space Studies (GISS) coupled ModelE simulations for June- September 1998-2002 are used to form lateral boundary conditions for synchronous simulations by the GISS RM3 regional climate model. The ModelE computational grid spacing is 2deg latitude by 2.5deg longitude and the RM3 grid spacing is 0.44deg. ModelE precipitation climatology for June-September 1998-2002 is shown to be a good proxy for 30-year means so results based on the 5-year sample are presumed to be generally representative. Comparison with observational evidence shows several discrepancies in ModelE configuration of the boreal summer inter-tropical convergence zone (ITCZ). One glaring shortcoming is that ModelE simulations do not advance the West African rain band northward during the summer to represent monsoon precipitation onset over the Sahel. Results for 1998-2002 show that onset simulation is an important added value produced by downscaling with RM3. ModelE Eastern South Atlantic Ocean computed sea-surface temperatures (SST) are some 4 K warmer than reanalysis, contributing to large positive biases in overlying surface air temperatures (Tsfc). ModelE Tsfc are also too warm over most of Africa. RM3 downscaling somewhat mitigates the magnitude of Tsfc biases over the African continent, it eliminates the ModelE double ITCZ over the Atlantic and it produces more realistic orographic precipitation maxima. Parallel ModelE and RM3 simulations with observed SST forcing (in place of the predicted ocean) lower Tsfc errors but have mixed impacts on circulation and precipitation biases. Downscaling improvements of the meridional movement of the rain band over West Africa and the configuration of orographic precipitation maxima are realized irrespective of the SST biases.

Regional climate models↗

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗

Ar-Ar and I-Xe Ages and the Thermal History of IAB Meteorites

Studies of several samples of the large Caddo County IAB iron meteorite reveal andesitic material, enriched in Si, Na, Al and Ca, which is essentially unique among meteorites. This material is believed to have formed from a chondritic source by partial melting and to have further segregated by grain coarsening. Such an origin implies extended metamorphism of the IAB parent body. New Ar-39- Ar-40 ages for silicate from three different Caddo samples are consistent with a common age of 4.50-4.51 Gyr ago. Less well defined Ar-Ar degassing ages for inclusions from two other IABs, EET8333 and Udei Station, are approx.4.32 Gyr, whereas the age for Campo del Cielo varies considerably over approx.3.23-4.56 Gyr. New I-129-Xe-129 ages for Caddo County and EET8333 are 4557.9+/-0.1 Myr and 4557-4560 Myr, respectively, relative to an age of 4562.3 Myr for Shallowater. Considering all reported Ar-Ar degassing ages for IABs and related winonaites, the range is approx.4.32-4.53 Gyr, but several IABs give similar Ar ages of 4.50-4.52 Gyr. We interpret these older Ar ages to represent cooling after the time of last significant metamorphism on the parent body, and the younger ages to represent later 40Ar diffusion loss. The older Ar-Ar ages for IABs are similar to Sm-Nd and Rb-Sr isochron ages reported in the literature for Caddo County. Considering the possibility that IAB parent body formation was followed by impact disruption, reassembly, and metamorphism (e.g., Benedix et al. 2000), the Ar-Ar ages and IAB cooling rates deduced from Ni concentration profiles in IAB metal (Herpfer et al., 1994) are consistent if the time of the post-assembly metamorphism was as late as approx.4.53 Gyr ago. However, I-Xe ages reported for some IABs define much older ages of approx.4558-4566 Myr, which cannot easily be reconciled with the much younger Ar-Ar and Sm-Nd ages. An explanation for the difference in radiometric ages of IABs may reside in combinations of the following: a) I-Xe ages have very high closure temperatures and were not reset during metamorphism approx.4.53 Gyr ago; b) a bias exists in the 40K decay constants which makes these Ar-Ar ages approx.30 Myr too young; c) the reported Sm-Nd and Rb-Sr ages for Caddo are in error by amounts equal to or exceeding their reported 2-sigma uncertainties; and d) about 30 Myr after the initial heating that produced differentiation of Caddo silicate and mixing of silicate and metal, a mild metamorphism of the IAB parent body reset the Ar-Ar ages.

Bogard, Donald D.↗

Subharmonic resonance in a mixing layer

The subharmonic resonance phenomenon in a spatially evolving mixing layer is studied using direct simulations of the 2-D Navier-Stokes equations. The computational domain extends to + or - infinity in the cross-stream direction with U(1) = 1.25 and U(2) = 0.25 imposed at + infinity and - infinity respectively. The domain is finite in the streamwise direction with inflow and outflow boundary conditions imposed at x/delta sub omega = 0 and 100, respectively. A hyperbolic-tangent mean velocity profile is assumed at the inlet and the Reynolds number based on the inlet vorticity thickness and velocity difference is Re = 600. It is observed that the phase angle between the fundamental and its subharmonic plays a key role in the spatial development of these modes. Contour plots of vorticity show that varying the phase will have a dramatic effect on the dynamics of the vortices. Pairing or shredding is observed depending on the phase. Fourier decomposition of the time traces show that the fundamental grows, saturates and decays with the downstream distance. The subharmonic has a similar behavior. However, the level at which the modes will saturate is affected by the phase. At 0 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is enhanced. At 90 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is inhibited. In the later case, the growth rate of the subharmonic recovers after saturation of the fundamental. These results are in qualitative agreement with experimental data.

Mansour, N. N.↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗