Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern formation”

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.

1,017 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Composition of Hydrothermal Vent Microbial Communities as Revealed by Analyses of Signature Lipids, Stable Carbon Isotopes and Aquificales Cultures

Extremely thermophilic microbial communities associated with the siliceous vent walls and outflow channel of Octopus Spring, Yellowstone National Park, have been examined for lipid biomarkers and carbon isotopic signatures. These data were compared with that obtained from representatives of three Aquificales genera. Thermocrinis ruber. "Thermocrinis sp. HI", Hydrogenobacter thermophilus TK-6, Aquifex pyrophilus and Aquifex aeolicus all contained phospholipids composed not only of the usual ester-linked fatty acids, but also ether-linked alkyls. The fatty acids of all cultured organisms were dominated by a very distinct pattern of n-C-20:1 and cy-C-21 compounds. The alkyl glycerol ethers were present primarily as CIS() monoethers with the expection of the Aquifex spp. in which dialkyl glycerol ethers with a boarder carbon-number distribution were also present. These Aquificales biomarker lipids were the major constituents in the lipid extracts of the Octopus Spring microbial samples. Two natural samples, a microbial biofilm growing in association with deposition of amorphous silica on the vent walls at 92 C, and the well-known 'pink-streamers community' (PSC), siliceous filaments of a microbial consortia growing in the upper outflow channel at 87 C were analyzed. Both the biofilm and PSC samples contained mono and dialkyl glycerol ethers with a prevalence of C-18 and C-20 alkyls. Phospholipid fatty acids were comprised of both the characteristic Aquificales n-C-20:1 and cy-C-21, and in addition, a series of iso-branched fatty acids from i-C-15:0 to i-C-21:0, With i-C-17:0 dominant in the PSC and i-C-19:0 in the biofilm, suggesting the presence of two major bacterial groups. Bacteriohopanepolyols were absent and the minute quantities of archaeol detected showed that Archaea were only minor constituents. Carbon isotopic compositions of the PSC yielded information about community structure and likely physiology. Biomass was C-13-depleted (10.9%) relative to available CO2 from the source water inorganic carbon pool with lipids further depleted by 6.3% relative to biomass The C-20-21 Aquificales fatty acids of the PSC were somewhat heavier than the iso-branched fatty acids. The carbon isotopic signatures of lipid biomarkers were also explored using a pure culture, T ruber, previously isolated from the PSC. Cells grown on C02 with O2 and both H2 and thiosulfate as electron donors were only slightly depleted (3.3%) relative to the C-source while cells grown on formate with O2 showed a major discrimination (19.7%), possibly the result of a metabolic branch point involving the assimilation of C-formate to biomass and the dissimilation to CO2 associated with energy production. T. ruber lipids were slightly heavier than biomass (+1.3%) whether cells were grown using CO2 or formate. Fatty acids from CO2 grown T. ruber cells were a so slightly heavier (average +2.1%) than biomass. The relatively depleted PSC C-20-21 fatty acids suggest that any associated Thermocrinis biomass would also be similarly depleted and much too light to be explained by growth on CO2. The C-fractionations determined with the pure culture suggest that growth of Thermocrinis in the PSC is more likely to occur on formate, presumably generated by geothermal activity. This study points to the value of the analysis of the structural and isotopic composition of lipid blomarkers both in pure culture studies, and in establishing community structure and physiology, as a complement to genomic profiles of microbial diversity. This is especially so when the members of the microbial community are novel and difficult to cultivate in the laboratory.

Jahnke, Linda L.↗
NASA NTRS2020

Durability of YSZ Coated Ti2AlC in 1300 °C Mach 0.3 Burner Rig Tests

A thermal barrier coating system survived burner rig testing at 1300 °C for 500 h. A 160 µm thick yttria stabilized zirconia (YSZ) coating was applied to a Ti2AlC MAX phase bar sample by plasma spray physical vapor deposition (PS-PVD) and tested face-on in an atmospheric Mach 0.3 jet fuel burner, using 5-h thermal cycles. No thermal barrier coating (TBC) spallation or recession was observed, only a 2.4 mg/cm2 mass gain. The modest weight gain precluded severe volatility losses under high velocity burner conditions. The coating surface exhibited colonies of (111)flourite fiber-textured columns separated by craze patterns, with no visible moisture attack. The metastable tetragonal t' YSZ phase was obtained initially, transitioning to equilibrium teq and cubic YSZ, but with little detrimental monoclinic. The thickness of the alumina TGO was ~21 to 23 μm under the heated YSZ face and ~13 to 15 μm on the uncoated, cooler backside. The backside exhibited removal of initial transient TiO2 nodules and partial etching of the underlying Al2O3 scale by volatile hydroxides formed in high temperature, high velocity water vapor. Aerodynamic forces produced some bending of the cantilevered sample via creep. The test indicated exceptional stability of YSZ coatings on Ti2AlC under turbine conditions, with thermal expansion matching playing a key role. The purpose of this study was to demonstrate long term durability of YSZ/MAX phase system in aggressive high temperature burner rig testing. MAX phases have been keenly studied because of their unique crystal structure and intriguing properties (Refs. 1 and 2). Having Mn+1(Al,Si)(C,N)n general composition, they are defined as ceramics, but possess unusual desirable attributes such as high conductivity, thermal shock resistance, easy machinability, and deformation tolerance. The mechanical properties derive from weak M-(Al,Si) bonding in the basal plane that leads to sliding and kinking in preference to catastrophic crack growth. Like most ceramics they are phase stable at high temperatures, generally up to 1500 °C. High temperature oxidation resistance is excellent for alumina-forming Ti3AlC2, Ti2AlC, and Cr2AlC, as reviewed by Tallman, et al. (Ref. 3). Compatibility with α-Al2O3 scales is further enhanced in cyclic exposures by a close matching of thermal expansion coefficients, (Ref. 4) i.e., (~9.3, 10.2, 11.3×10–6/K for Al2O3, Ti2AlC, and YSZ, to be discussed). Turbine environments generally contain 10 percent water vapor in the combustion gases, therefore moisture effects can be a concern for some materials (Ref. 5). Furnace tests of MAX phases in high temperature steam generally showed little effect on Al2O3 scale growth (Ref. 6). However, high velocity and high pressure gas can influence scale losses by the formation of volatile reaction products, such as TiO(OH)2 and Al(OH)3 (Refs. 7 to 10). This phenomenon had been discussed for 1100 to 1300 °C high pressure burner rig tests of Ti2AlC (Ref. 11). A single cubic growth rate parameter kcubic was measurably lower than comparable furnace TGA data, but it could be matched reasonably well if corrected for a slight volatility term. In general, a two-parameter cubic-linear growth-volatility law was believed to apply. Corresponding scale volatility loss rates, directly measured at 1300 °C on a pre-oxidized sample, were moderate (0.012 mg/cm2/h) and largely attributed to removal of the initial TiO2 transient scale. A related CH4 burner study of high purity Cr2AlC MAX phase demonstrated 1200 °C durability after 500 rapid (5 min. heat and 2 min. cool) thermal shock cycling (29 h hot time) (Ref. 12). Heating and cooling rates were ~1000 and 500 °C per minute, with a gas velocity of 5 m/s, producing a 75 °C/mm gradient. A 7 μm Al2O3 surface scale and a 13 μm Cr7C3 depletion zone formed with no signs of failure. No evidence of scale volatility was evident, although weight change was not provided, the velocity was moderate, and the total hot time was not extensive. The same high gradient BRT was used to produce 1400 °C surface temperatures for a YSZ/Cr2AlC/IN738 system in the first study of MAX phases used as bond coats for thermal barrier coatings (TBC) (Ref. 25). Here TBC failure was reported after 745 cycles, with only a 1.5 μm Al2O3 scale entrained within a porous, Cr7C3 bondcoat depletion phase. YSZ thermal barrier coatings have been considered to be a compatible complement to Al-MAX phases because of thermal expansion matching and extremely low volatility in water vapor. Initial studies showed superior oxidative stability up to 1300 °C, for long times (at least 500 h) for Ti2AlC substrates and less (268 h) for Cr2AlC, while withstanding large alumina TGO scale thickness (~35 to 40 μm) (Refs. 13 and 14). By comparison, typical superalloy systems can only survive 1150 °C maximum interface temperatures for extended periods, with a maximum sustained TGO below 10 μm (Ref. 15). High temperature SiC based systems are known to form slow-growing SiO2 scales. But these are subject to rate enhancement and volatile Si(OH)4 products in the presence of water vapor, as described comprehensively by Opila, et al. (Refs. 5, 16 to 19). Net weight losses are generally observed in high velocity, high pressure burner rig studies (e.g., 0.084 mg/cm2/h at 1300 °C) (Ref. 20). Furthermore, the loss rates have been shown from chemical physics to scale with v1/2 and pH2O 2 (Ref. 16). Low activity, moisture-resistant environmental barrier coatings (EBC), such as rare earth silicates, are needed to prevent substrate recession under turbine conditions (Refs. 21 to 23).

Smialek, James L.↗
NASA NTRS2012

Engineered Multifunctional Surfaces for Fluid Handling

Designs incorporating variations in capillary geometry and hydrophilic and/or antibacterial surface properties have been developed that are capable of passive gas/liquid separation and passive water flow. These designs can incorporate capillary grooves and/or surfaces arranged to create linear and circumferential capillary geometry at the micro and macro scale, radial fin configurations, micro holes and patterns, and combinations of the above. The antibacterial property of this design inhibits the growth of bacteria or the development of biofilm. The hydrophilic property reduces the water contact angle with a treated substrate such that water spreads into a thin layer atop the treated surface. These antibacterial and hydrophilic properties applied to a thermally conductive surface, combined with capillary geometry, create a novel heat exchanger capable of condensing water from a humid, two-phase water and gas flow onto the treated heat exchanger surfaces, and passively separating the condensed water from the gas flow in a reduced gravity application. The overall process to generate the antibacterial and hydrophilic properties includes multiple steps to generate the two different surface properties, and can be divided into two major steps. Step 1 uses a magnetron-based sputtering technique to implant the silver atoms into the base material. A layer of silver is built up on top of the base material. Completion of this step provides the antibacterial property. Step 2 uses a cold-plasma technique to generate the hydrophilic surface property on top of the silver layer generated in Step 1. Completion of this step provides the hydrophilic property in addition to the antibacterial property. Thermally conductive materials are fabricated and then treated to create the antibacterial and hydrophilic surface properties. The individual parts are assembled to create a condensing heat exchanger with antibacterial and hydrophilic surface properties and capillary geometry, which is capable of passive phase separation in a reduced gravity application. The plasma processes for creating antibacterial and hydrophilic surface properties are suitable for applications where water is present on an exposed surface for an extended time, such that bacteria or biofilms could form, and where there is a need to manage the water on the surface. The processes are also suitable for applications where only the hydrophilic property is needed. In particular, the processes are applicable to condensing heat exchangers (CHXs), which benefit from the antibacterial properties as well as the hydrophilic properties. Water condensing onto the control surfaces of the CHX will provide the moist conditions necessary for the growth of bacteria and the formation of biofilms. The antibacterial properties of the base layer (silver) will mitigate and prevent the growth of bacteria and formation of biofilms that would otherwise reduce the CHX performance. In addition, the hydrophilic properties reduce the water contact angle and prevent water droplets from bridging between control surfaces. Overall, the hydrophilic properties reduce the pressure drop across the CHX.

Thomas, Chris↗
NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗
NASA NTRS1997

Biogeochemistry of DMS in Surface Waters

Dimethylsulfide (DMS) is important in influencing the formation of aerosols in the troposphere over large areas of the world's oceans. Understanding the dynamics of aerosols is important to understanding the earth's radiation balance. In evaluating the factors controlling DMS in the troposphere it is vital to understand the dynamics of DMS in the surface ocean. The biogeochemical processes controlling DMS concentration in seawater are myriad; modeling and theoretical estimation are problematic. At the beginning of this project we believed that we were on the verge of simplifying the ship-track measurement of DMS, and we proposed to deploy such a system to develop a database relating high frequency DMS measurements to biological and physicochemical and optical properties of surface water that can be quantified by remote sensing techniques. We designed a system to measure DMS concomitantly with other basic chemical and biological data in a flow-through system. The project was collaborative between Woods Hole Oceanographic Institution (WHOI) and Bermuda Biological Station for Research (BBSR). The project on which we are reporting was budgeted for only one year with a one year no-cost extension. At WHOI our effort was directed towards designing traps which would be used to concentrate DMS from seawater and allow storage for subsequent analysis. At that time, GC systems were too large for easy long-term deployment on a research vessel like R/V Weatherbird, so we focused on simplifying the shipboard sampling procedure. Initial studies of sample recovery with high levels of DMS suggested that Carboxen 1000, a relatively new carbon molecular sieve, could be used as a stable storage medium. The affinity of Carboxen for DMS is several orders of magnitude higher than gold wool (another adsorbent used for DMS collection) on a weight or volume basis. Furthermore, Carboxen's affinity for DMS is also far less susceptible to humidity than gold wool. Unfortunately, further experiments with low level DMS indicated that recovery of DMS after storage was not quantitative. The material has proven to be completely acceptable for short term storage and has been incorporated into a micro-GC system. Since working on this project, we have collaborated with RVM Scientific in Santa Barbara in the design and construction of small portable micro-GC's that will make feasible at-sea measurement in moving ships, making rapid gas analysis and quantification feasible in a ship-track mode. Throughout this period at both WHOI and BBSR, we continued to analyze field data to understand that patterns of time and space variability in DMS and the processes that govern it. These insights will be crucial to determining the specifications for our automated sampling program. The data from this, the longest continuous sampling program for ocean DMS, provided insights into year to year and short-term variability.

Dacey, J. W. H.↗
NASA NTRSDate not supplied

Trace Elemental Abundances in Calcium-Aluminum-Rich Inclusions in CV Chondrites

Introduction: Calcium-aluminum-rich inclusions (CAIs), are the first formed solids that define the age of the Solar System [1,2]. CAIs are thought to have condensed from nebular gas [3,4] within the first <1 Ma of Solar System formation [5,6]. CAIs have experienced numerous early Solar System processes including condensation, evaporation, melting, recrystallization, and aqueous alteration [e.g., 7]. The chemical, mineralogical, and textural diversity among CAIs results from a range of chemical and physical processes recorded during nebular and parent body epoch. This study aims to explore the mineralogical, textural, and chemical compositions of CAIs including the trace elemental abundances in CAI phases to determine the early Solar System processes recorded in them. Samples and Analytical Methods: We analyzed one CAI each from CV3 chondrites Northwest Africa (NWA) 5508 designated as ‘Saguaro’, and Northwest Africa (NWA) 12772 designated as ‘Hoopoe’. Back-scatter electron (BSE) images were collected using a Phenom XL scanning electron microscope (SEM) at the Lunar and Planetary Institute (LPI) and the JEOL JXA-8530F electron probe microanalyzer (EPMA) at Johnson Space Center (JSC)-NASA. Additionally, energy dispersive X-ray spectrometry (EDS) elemental maps of select areas for these samples were collected using a 15.0kV beam energy and a 40µA emission current. Using the EPMA, wavelength-dispersive X-ray spectroscopy (WDS) quantitative data were collected. In-situ trace element measurements for both CAIs were determined at JSC-NASA using a Photon Machines 193nm laser ablation system and a Thermo-Scientific Element-XR inductively coupled plasma mass spectrometer (ICP-MS). Analyses consisted of 30s ablations at 10Hz, spot sizes of 20-25µm, and a fluence of 6.0 J/cm2 for anorthite and melilite, and a 3.5 J/cm2 fluence for all other phases. NIST612 was used to correct for instrument drift, while BHVO-2g was used as a primary calibration standard. BCR-2g and in-house mineral standards were regularly measured as unknowns to ensure accuracy. Results: Saguaro is a coarse-grained CAI, ~11 x 6 mm in dimensions. Saguaro contains spinel, Al-rich pyroxene, anorthite, Mg-rich melilite, and minor perovskite in its interior and is therefore classified as a Type B CAI. Individual melilite grains shows normal compositional zoning with an Ak content ranging from ~24 to 54 with no apparent trend from the core to the edge of the CAI. The spinel appears euhedral and occurs both as clusters and as spinel palisades [8]. Two rim sequences surround most of the sample: the inner rim being a Wark-Lovering (WL) rim (~10-35 µm) containing pyroxene, spinel, and melilite (or anorthite), and the outer rim is a finer-grained, thicker (~100 µm), accretionary rim (Fig. 1). The mineral phases in Saguaro record an overall flat REE pattern with an average negative Eu anomaly in pyroxene, and an average positive Eu anomaly in anorthite and melilite respectively. Anorthite, melilite, and pyroxene have a minor depletion in Tm (Fig. 2). The Hoopoe CAI is a compact, coarse-grained ~6 × 4 mm in size. The major mineralogy includes hibonite, spinel, melilite, anorthite, and perovskite. Therefore, it is classified as a compact transitional type A and B. (?)zoning was observed in some hibonites. Individual melilite grains show both reverse and normal zoning, where the Ak content ranges from ~6- to 28. Melilite shows two distinct textures. One texture consisted of smooth melilite that appeared homogenous, while the second appeared to consist of many fine fractures. The spinel also often appears clustered. The WL-rim sequence surrounding Hoopoe is ~25 µm thick and composed of spinel, perovskite, hibonite, and melilite/anorthite. It is then partially surrounded by an outer accretionary rim (~75µm). Like before, refractory metal nuggets appeared concentrated near the WL rims. Other metal assemblages rich in Fe and Ni were also observed. All major mineral phases in Hoopoe display relatively flat REE patterns, except for varying Eu and Tm between phases (Fig. 2). There is a prominent negative Eu anomaly in perovskite and an average positive Eu anomaly in spinel, anorthite, and hibonite respectively (Fig. 2, 3). The mixed phases along the rim of the CAI also display a negative Eu anomaly, and all phases the CAI were depleted in Pb. Discussion: The CV3 CAIs analyzed in this study were classified based on their mineralogy and textures into Type A versus Type B CAIs [10]. Hibonite appears to be pseudomorphically replacing the spinel, (i.e., is hibonite in composition, but appears in the shape of spinel). Spinel palisades. The presence of spinel palisades present in Saguaro are consistent with the melting and recrystallization experienced by this CAI. Trace elemental analyses. Saguaro and Hoopoe display similar trace element patterns to each other, with both appearing generally flat, with anomalies in Eu, and Tm. Melilite and anorthite display positive Eu anomalies in both CAIs, in addition to the hibonite in Hoopoe (Fig. 2). The phases that are depleted in Eu are pyroxene and perovskite in both Saguaro and Hoopoe, respectively (Fig. 3). Given that Eu is volatile in reducing environments [11], this could possibly indicate reducing conditions at the time anorthite and melilite crystallized, with the gas they formed from containing Eu. As these CAIs continued to form, this gas as a result would become depleted in Eu. This also could be supported by the propensity of anorthite and melilite to take up Eu from its surroundings and incorporate it into their structure [12]. In addition, analyzing the assemblage of the phases in the Saguaro, melilite and anorthite (Eu enriched) often surround the pyroxene (Eu depleted) as they are crystallized. This intergrowth of phases and the proximity of the phases would support that the Eu is being incorporated into some phases, preventing it from incorporating into other. Trace elemental analyses of the CAI rims will be evaluated in more detail, as they are complicated by the transient signal being composed of a mixture of mineral phases. Broadly, however, the patterns in the rims of both CAIs are comparable to each other, and for Hoopoe, to the mixed phase patterns in the core (Fig. 3). Other studies have found that CAI rims can be depleted in Ce and Yb [13], however we did not observe these anomalies in the two CAIs discussed here. Given their similarity, the trace elemental analyses of the mixed interior (i.e. core) and rim phases could be interpreted as forming from similar, if not the same, reservoirs. The REE abundance between the rim and core of Hoopoe are also similar, indicating they may have formed from a gas of the same or similar composition. Acknowledgments: We thank the ASU Center for Meteorite Studies for loaning the samples used in this work and Tabb Prissel for his assistance with the analysis. Mouti Al-Hashimi thanks Sam Crossley and Cyrena Goodrich for their help with the LPI SEM training. This work was supported by the LPI Summer Intern Program in Planetary Science and the LPI Cooperative Agreement. References: [1] Connelly J.N. (2012) Science, 338, 651-655. [2] MacPherson G. J. (2014) Treatise on Geochem., 2, 139-179. [3] Grossman L. (1972) GCA, 36, 597-619. [4] Ebel, D.S. (2006) Meteorites and the Early Solar System II (D. S. Lauretta & H. Y. McSween, Eds.) 253-277. [5] MacPherson G. J. (2012) Earth Planet. Sci. Lett., 331-332, 43-54. [6] MacPherson G.J. (2017) GCA, 201, 65-82. [7] Krot A.N. (1995) Meteoritics & Planet. Sci., 30, 748-775. [8] Wark and Lovering (1982) GCA, 46, 2595-2607. [9] Palme H. and Jones A. (2003) Treatise on Geochemistry (H. D. Holland and K. K. Turekian Eds.), 1, 41-61. [10] Grossman L. (1980) Ann. Rev. Earth Planet. Sci., 8, 559-608. [11] Floss C. et al. (1996) GCA, 60, 1975-1997. [12] Mason B. and Martin P. M. (1974) Earth Planet. Sci. Lett., 22, 141-144. [13] Wark B. and Boynton W. V. (2001) Meteoritics & Planet. Sci., 36, 1135-1166.

X Mouti↗
NASA NTRSDate not supplied

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗
NASA NTRSDate not supplied

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗
NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗