Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic stabilization”

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.

391 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Environmental Stability and Oxidation Behavior of HfO2-Si and YbGd(O) Based Environmental Barrier Coating Systems for SiCSiC Ceramic Matrix Composites

Ceramic environmental barrier coatings (EBC) and SiCSiC ceramic matrix composites (CMCs) will play a crucial role in future aircraft propulsion systems because of their ability to significantly increase engine operating temperatures, improve component durability, reduce engine weight and cooling requirements. Advanced EBC systems for SiCSiC CMC turbine and combustor hot section components are currently being developed to meet future turbine engine emission and performance goals. One of the significant material development challenges for the high temperature CMC components is to develop prime-reliant, environmental durable environmental barrier coating systems. In this paper, the durability and performance of advanced Electron Beam-Physical Vapor Deposition (EB-PVD) NASA HfO2-Si and YbGdSi(O) EBC bond coat top coat systems for SiCSiC CMC have been summarized. The high temperature thermomechanical creep, fatigue and oxidation resistance have been investigated in the laboratory simulated high-heat-flux environmental test conditions. The advanced NASA EBC systems showed promise to achieve 1500C temperature capability, helping enable next generation turbine engines with significantly improved engine component temperature capability and durability.

Ceramic Matrix Composites↗
NASA NTRS1995

Cytoskeletal filament assembly and the control of cell spreading and function by extracellular matrix

This study was undertaken to analyze how cell binding to extracellular matrix produces changes in cell shape. We focused on the initial process of cell spreading that follows cell attachment to matrix and, thus, cell 'shape' changes are defined here in terms of alterations in projected cell areas, as determined by computerized image analysis. Cell spreading kinetics and changes in microtubule and actin microfilament mass were simultaneously quantitated in hepatocytes plated on different extracellular matrix substrata. The initial rate of cell spreading was highly dependent on the matrix coating density and decreased from 740 microns 2/h to 50 microns 2/h as the coating density was lowered from 1000 to 1 ng/cm2. At approximately 4 to 6 hours after plating, this initial rapid spreading rate slowed and became independent of the matrix density regardless of whether laminin, fibronectin, type I collagen or type IV collagen was used for cell attachment. Analysis of F-actin mass revealed that cell adhesion to extracellular matrix resulted in a 20-fold increase in polymerized actin within 30 minutes after plating, before any significant change in cell shape was observed. This was followed by a phase of actin microfilament disassembly which correlated with the most rapid phase of cell extension and ended at about 6 hours; F-actin mass remained relatively constant during the slow matrix-independent spreading phase. Microtubule mass increased more slowly in spreading cells, peaking at 4 hours, the time at which the transition between rapid and slow spreading rates was observed. However, inhibition of this early rise in microtubule mass using either nocodazole or cycloheximide did not prevent this transition. Use of cytochalasin D revealed that microfilament integrity was absolutely required for hepatocyte spreading whereas interference with microtubule assembly (using nocodazole or taxol) or protein synthesis (using cycloheximide) only partially suppressed cell extension. In contrast, cell spreading could be completely inhibited by combining suboptimal doses of cytochalasin D and nocodazole, suggesting that intact microtubules can stabilize cell form when the microfilament lattice is partially compromised. The physiological relevance of the cytoskeleton and cell shape in hepatocyte physiology was highlighted by the finding that a short exposure (6 hour) of cells to nocodazole resulted in production of smaller cells 42 hours later that exhibited enhanced production of a liver-specific product (albumin). These data demonstrate that spreading and flattening of the entire cell body is not driven directly by net polymerization of either microfilaments or microtubules.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Cell Biology↗
NASA NTRS2011

High Power Flex-Propellant Arcjet Performance

A MW-class electrothermal arcjet based on a water-cooled, wall-stabilized, constricted arc discharge configuration was subjected to extensive performance testing using hydrogen and simulated ammonia propellants with the deliberate aim of advancing technology readiness level for potential space propulsion applications. The breadboard design incorporates alternating conductor/insulator wafers to form a discharge barrel enclosure with a 2.5-cm internal bore diameter and an overall length of approximately 1 meter. Swirling propellant flow is introduced into the barrel, and a DC arc discharge mode is established between a backplate tungsten cathode button and a downstream ringanode/ spin-coil assembly. The arc-heated propellant then enters a short mixing plenum and is accelerated through a converging-diverging graphite nozzle. This innovative design configuration differs substantially from conventional arcjet thrusters, in which the throat functions as constrictor and the expansion nozzle serves as the anode, and permits the attainment of an equilibrium sonic throat (EST) condition. During the test program, applied electrical input power was varied between 0.5-1 MW with hydrogen and simulated ammonia flow rates in the range of 4-12 g/s and 15-35 g/s, respectively. The ranges of investigated specific input energy therefore fell between 50-250 MJ/kg for hydrogen and 10-60 MJ/kg for ammonia. In both cases, observed arc efficiencies were between 40-60 percent as determined via a simple heat balance method based on electrical input power and coolant water calorimeter measurements. These experimental results were found to be in excellent agreement with theoretical chemical equilibrium predictions, thereby validating the EST assumption and enabling the utilization of standard TDK nozzle expansion analyses to reliably infer baseline thruster performance characteristics. Inferred specific impulse performance accounting for recombination kinetics during the expansion process implied nearly frozen flow in the nozzle and yielded performance ranges of 800-1100 sec for hydrogen and 400-600 sec for ammonia. Inferred thrust-to-power ratios were in the range of 30-10 lbf/MWe for hydrogen and 60-20 lbf/MWe for ammonia. Successful completion of this test series represents a fundamental milestone in the progression of high power arcjet technology, and it is hoped that the results may serve as a reliable touchstone for the future development of MW-class regeneratively-cooled flex-propellant plasma rockets.

Litchford, Ron J.↗
NASA NTRS1996

A Magnetic Bumper-Tether System Using ZFC Y123

We consider the use of magnetic forces in a bumper system, to soften docking procedures. We investigate a system which exhibits no magnetic field except during the docking process, which, if desired, can automatically tether two craft together, and which provides lateral stability during docking. A system composed of zero field cooled Y(1.7)Ba2Cu3O(7-delta) (Y123) tiles and electromagnets is proposed. The Y123 high temperature superconductor (HTS) is mounted on one craft, and the electromagnet on the other. Results of small prototype laboratory experiments are reported. The electromagnet has, for convenience, been replaced by a permanent SmCo ferromagnet in these measurements. When the two craft approach, a mirror image of the ferromagnet is induced in the Y123, and a repulsive bumper force, F(sub B), results. F(sub B) is velocity dependent, and increases with v. For presently available HTS materials, bumper pressure of approx. 3.7 N/cm(exp 2) is achieved using SmCo. This extrapolates to approx. 18 N/cm(exp 2) for an electromagnet, or a force of up to 20 tons for a 1 m(exp 2) system. After reaching a minimum distance of approach, the two colliding craft begin to separate. However, the consequent change of SmCo magnetic field at the Y123 results in a reversal of current in the Y123 so that the Y123 is attractive to the SmCo. The attractive (tether) force, F(sub T), is a function of R = B(sub Fe)/B(sub t, max), where B(sub Fe) is the field at the surface of the ferromagnet, and B(sub t, max) is the maximum trapped field of the Y123, i.e., the trapped field in the so-called critical state. For R greater than or equal to 2, F(sub T) saturates at a value comparable to F(sub B). For a range of initial approach velocities the two craft are tethered following the bumper sequence. Most of the kinetic energy of the collision is first converted to magnetic field energy in the Y123, and then into heat via the creep mechanism. About 15% of the work done against magnetic forces during collision remains stored as magnetic energy after 1 hour. Experiments have also been conducted on the spatial range of the bumper force for arrays of HTS tiles. For a single HTS tile approx. 2 cm in diameter, the range of F(sub B) is approx. l cm. For a l m(exp 2) array the range of F(sub B) will be circa 50 cm.

Weinstein, Roy↗
NASA NTRS2014

Robonaut 2 - IVA Experiments On-Board ISS and Development Towards EVA Capability

Robonaut 2 (R2) has completed its fixed base activities on-board the ISS and is scheduled to receive its climbing legs in early 2014. In its continuing line of firsts, the R2 torso finished up its on-orbit activities on its stanchion with the manipulation of space blanket materials and performed multiple tasks under teleoperation control by IVA astronauts. The successful completion of these two IVA experiments is a key step in Robonaut's progression towards an EVA capability. Integration with the legs and climbing inside the ISS will provide another important part of the experience that R2 will need prior to performing tasks on the outside of ISS. In support of these on-orbit activities, R2 has been traversing across handrails in simulated zero-g environments and working with EVA tools and equipment on the ground to determine manipulation strategies for an EVA Robonaut. R2 made significant advances in robotic manipulation of deformable materials in space while working with its softgoods task panel. This panel features quarter turn latches that secure a space blanket to the task panel structure. The space blanket covers two cloth cubes that are attached with Velcro to the structure. R2 was able to open and close the latches, pull back the blanket, and remove the cube underneath. R2 simulated cleaning up an EVA worksite as well, by replacing the cube and reattaching the blanket. In order to interact with the softgoods panel, R2 has both autonomously and with a human in the loop identified and localized these deformable objects. Using stereo color cameras, R2 identified characteristic elements on the softgoods panel then extracted the location and orientation of the object in its field of view using stereo disparity and kinematic transforms. R2 used both vision processing and supervisory control to successfully accomplish this important task. Teleoperation is a key capability for Robonaut's effectiveness as an EVA system. To build proficiency, crewmembers have attempted increasingly difficult tasks using R2 inside the Station. After donning motion capture equipment and a virtual reality visor, Expedition 34/35 flight engineer Tom Marshburn began operations with simple hand movements. Having gained confidence, Marshburn guided R2's arms in a leader-follower exercise with crewmate Chris Cassidy. He was also able to use the hand to grab a tumbling roll of tape, a task only demonstrable in microgravity. Later efforts saw Cassidy handle softgoods through shared control with ground operators, mimicking an activity previously achieved using only autonomy. Robotic climbing through the ISS on handrails requires both precision motion and compliant grasps in order to both position grippers on handrails/seat track and prevent large internal forces. R2 climbs using actively controlled compliance and torque limiting to meet both the precision and softness requirements. During a step, the attached leg is controlled to be strong and stiff in order to maintain precision trajectory tracking. The swing leg is controlled to be stiff but weak to minimize unintentional impact forces while maintaining precision. During a simulated dual limb grasp (as shown in Figure 1), the R2 controller maintains one limb rigid and one limb soft to prevent large internal forces from building up. R2's grippers also use a form of force control to limit grip force while not fully closed on either a handrail or seat track thus limiting unintentional forces on cables/objects that may be present in R2's translational path. The on-board torso R2 safety system relies on a single end-effector velocity limit to prevent potential impact forces from exceeding Station maximum load requirements. R2's mobile configuration required modifications to the velocity limiting safety function due to its large, dynamic inertia. R2's legs maneuver the robot's mass creating configuration dependent, joint-relative inertias. A single all-encompassing velocity limit to cover worst case inertia is prohibitively low. The upgraded R2 control and safety systems solve this problem using momentum limiting, momentum control, and kinetic energy minimization. Momentum and kinetic energy take the robot mass into account relieving low velocity restrictions on low inertia end-effectors while ensuring that the overall mass of R2 is limited from hazardous velocities. The momentum of R2's five safety nodes (each of the four end-effectors and the body) is monitored and compared to a single momentum limit. If any of the five nodes exceeds the safety limit, the motor power is removed and the robot comes to a stop. Momentum control/limiting also provides a simple, reliable method to integrate hand held tools into the safety system by providing the tool mass to the control system thus automatically reducing the allowable velocity of the end-effector with the tool. Work on the ground continues to build the skill set for an EVA Robonaut. Recent experiments (Figure 2) demonstrate how a teleoperator can use R2 to manipulate a tether hook, an important safety precaution on spacewalks. Another task displayed Robonaut's ability to pull back a protective jacket over a hose and search for damage, as well as inspect a quick-disconnect fitting for debris. Demonstrations such as these are indicative of EVA work done on ISS, specifically seen during a series of spacewalks over 2012 and 2013 where astronauts searched for an ammonia leak in one of the external cooling loops. Through experiments both on ISS and on the ground, R2 is evolving and providing the information needed to plan out the upgrades that will make an EVA Robonaut an effective tool. With the addition of legs, R2 will start climbing inside the space station and supply invaluable information on how the climbing strategies and task stabilization techniques must be refined. Ground R2 systems will continue to work with additional EVA tools and equipment in preparation for onboard IVA testing and future EVA applications.

Diftler, Myron↗
NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗
NASA NTRS2006

Coating MCPs with AlN and GaN

A development effort underway at the time of reporting the information for this article is devoted to increasing the sensitivity of microchannel plates (MCPs) as detectors of photons and ions by coating the MCPs with nitrides of elements in period III of the periodic table. Conventional MCPs are relatively insensitive to slowly moving, large-mass ions for example, ions of biomolecules under analysis in mass spectrometers. The idea underlying this development is to coat an MCP to reduce its work function (decrease its electron affinity) in order to increase both (1) the emission of electrons in response to impingement of low-energy, large-mass ions and (2) the multiplying effect of secondary electron emission. Of particular interest as coating materials having appropriately low or even negative electron affinities are gallium nitride, aluminum nitride, and ternary alloys of general composition Al(x)Ga(1-x)N (where 0<x<1). These materials exhibit attractively high degrees of chemical, mechanical, and thermal stability plus acceptably high resistance to sputtering. The electron-excitation cross sections of these materials are expected to exceed those of several other materials (including diamond) that are, variously, in use or under development for the same purpose. Moreover, by doping these materials with silicon, one can render them partly electrically conductive, thereby suppressing the undesired accumulation of electric charge that could otherwise occur during bombardment by ions. For experiments, thin films of AlN and GaN both undoped and doped with Si were deposited on commercial MCPs by radio-frequency molecular-beam epitaxy (also known as plasma-assisted molecular-beam epitaxy) at temperatures <200 C. This deposition technique is particularly suitable because (1) MCPs cannot withstand the higher deposition-substrate temperatures used to decompose constituent compounds in some other deposition techniques and (2) in this technique, the constituent Al, Ga, and N are supplied in elemental form, so that there is no need for thermal decomposition at the substrate surface. The nitride films thus formed were, variously, amorphous or polycrystalline. The nitride films were coated with surface layers of gold <100 thick. The MCPs were tested in a standard configuration in which the output stage of a first MCP was coupled to the input stage of a second MCP. Each pair of MCPs was mounted in a standard holder that included front and back contact rings and an anode for collecting the output electrons of the second MCP. The MCP pairs were biased at potentials between 1.7 and 1.9 kV, and count rates measured after preamplification and discrimination. To enable a direct comparison, in one pair, the second MCP was uncoated while the first MCP was coated over half its surface. The coated and uncoated sides of the half-coated MCP were exposed to fluxes of argon ions at kinetic energies of 1.0 and 0.5 keV. At 1.0 keV, the count rate for the coated side was about 2.3 times greater than that for the uncoated side; at 0.5 keV, the count rate for the coated side was about 1.8 times greater than that for the uncoated side.

Bensaoula, Abdelhakim↗
NASA NTRS2016

Ice Melt, Sea Level Rise and Superstorms: Evidence from Paleoclimate Data, Climate Modeling, and Modern Observations that 2C Global Warming Could Be Dangerous

We use numerical climate simulations, paleoclimate data, and modern observations to study the effect of growing ice melt from Antarctica and Greenland. Meltwater tends to stabilize the ocean column, inducing amplifying feedbacks that increase subsurface ocean warming and ice shelf melting. Cold meltwater and induced dynamical effects cause ocean surface cooling in the Southern Ocean and North Atlantic, thus increasing Earth's energy imbalance and heat flux into most of the global ocean's surface. Southern Ocean surface cooling, while lower latitudes are warming, increases precipitation on the Southern Ocean, increasing ocean stratification, slowing deepwater formation, and increasing ice sheet mass loss. These feedbacks make ice sheets in contact with the ocean vulnerable to accelerating disintegration. We hypothesize that ice mass loss from the most vulnerable ice, sufficient to raise sea level several meters, is better approximated as exponential than by a more linear response. Doubling times of 10, 20 or 40 years yield multi-meter sea level rise in about 50, 100 or 200 years. Recent ice melt doubling times are near the lower end of the 10-40-year range, but the record is too short to confirm the nature of the response. The feedbacks, including subsurface ocean warming, help explain paleoclimate data and point to a dominant Southern Ocean role in controlling atmospheric CO2, which in turn exercised tight control on global temperature and sea level. The millennial (500-2000-year) timescale of deep-ocean ventilation affects the timescale for natural CO2 change and thus the timescale for paleo-global climate, ice sheet, and sea level changes, but this paleo-millennial timescale should not be misinterpreted as the timescale for ice sheet response to a rapid, large, human-made climate forcing. These climate feedbacks aid interpretation of events late in the prior interglacial, when sea level rose to C6-9m with evidence of extreme storms while Earth was less than 1 C warmer than today. Ice melt cooling of the North Atlantic and Southern oceans increases atmospheric temperature gradients, eddy kinetic energy and baroclinicity, thus driving more powerful storms. The modeling, paleoclimate evidence, and ongoing observations together imply that 2 C global warming above the preindustrial level could be dangerous. Continued high fossil fuel emissions this century are predicted to yield (1) cooling of the Southern Ocean, especially in the Western Hemisphere; (2) slowing of the Southern Ocean overturning circulation, warming of the ice shelves, and growing ice sheet mass loss; (3) slowdown and eventual shutdown of the Atlantic overturning circulation with cooling of the North Atlantic region; (4) increasingly powerful storms; and (5) nonlinearly growing sea level rise, reaching several meters over a timescale of 50-150 years. These predictions, especially the cooling in the Southern Ocean and North Atlantic with markedly reduced warming or even cooling in Europe, differ fundamentally from existing climate change assessments. We discuss observations and modeling studies needed to refute or clarify these assertions.

Computer simulation↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗
NASA NTRS2017

Correlations Among Microstructure, Morphology, Chemistry, and Isotopic Systematics of Hibonite in CM Chondrites

Introduction: Hibonite is a primary refractory phase occurring in many CAIs, typically with spinel and perovskite. Our microstructural studies of CAIs from carbonaceous chondrites reveal a range of stacking defect densities and correlated non-stoichiometry in hibonite. We also conducted a series of annealing experiments, demonstrating that the Mg-Al substitution stabilized the formation of defect-structured hibonite. Here, we continue a detailed TEM analysis of hibonite-bearing inclusions from CM chondrites that have been well-characterized isotopically. We examine possible correlations of microstructure, morphology, mineralogy, and chemical and isotopic systematics of CM hibonites in order to better understand the formation history of hibonite in the early solar nebula. Methods: Fifteen hibonite-bearing inclusions from the Paris CM chondrite were analyzed using a JEOL 7600F SEM and a JEOL 8530F electron microprobe. In addition to three hibonite-bearing inclusions from the Murchison CM chondrite previously reported, we selected three inclusions from Paris, Pmt1-6, 1-9, and 1-10, representing a range of 26Al/27Al ratios and minor element concentrations for a detailed TEM study. We extracted TEM sections from hibonite grains using a FEI Quanta 3D field emission gun SEM/FIB. The sections were then examined using a JEOL 2500SE field-emission scanning TEM equipped with a Thermo-Noran thin window EDX spectrometer. Results and Discussion: A total of six hibonite-bearing inclusions, including two platy hibonite crystals (PLACs) and four spinel-hibonite inclusions (SHIBs), were studied. There are notable differences in chemical and isotopic compositions between the inclusions (Table 1), indicative of their different formation environment or timing. Our TEM observations show perfectly-ordered, stoichiometric hibonite crystals without stacking defects in two PLACs, 2-7-1 and 2-8-2, and in three SHIBs, Pmt1-6, 1-9, and 1-10. In contrast, SHIB 1-9-5 hibonite grains contain a low density of stacking defects linked to an increase in MgO contents, indicating complex, disordered intergrowths of stoichiometric and MgO-enriched hibonites. From the data collected to date, we find no clear correlation between the microstructures of hibonite and its morphological and mineralogical types that reflect distinct chemical and isotopic systematics [6-8,10]. Interestingly, the presence of no or few stacking defects in hibonite from the PLACs and SHIBs are in contrast to our experimental studies that produced very high densities of stacking defects in hibonite [3-5]. Unlike our experi-ments, electron microprobe data from the PLACs and SHIBs hibonite grains show a strong correlation between (Ti4++Si4+) and Mg2+ cations, suggesting that coupled substitutions of (Ti4++Mg2+) and (Si4++Mg2+) for 2Al3+ inhibit the formation of defect-structured hibonite. However, our experimental studies suggest that kinetics (e.g., cooling rate) or other thermal effects also exert a strong control on the microstructures and chemical compositions of hibonite. In Pmt1-6, elongated perovskite grains present at the hibonite grain boundaries display (121) twinning, indicative of a fast cooling (>50degC/min) after high-temperature events. Therefore, the nebular microstructural characteristics of hibonite, at least in this inclusion, would not have destroyed by subsequent high-temperature annealing. Conclusions: Our TEM observations thus far show no clear correlation in microstructures, morphological and mineralogical characteristics, and chemical and isotopic systematics of hibonites from CM chondrites. The observed variation in stacking defect densities in the hibonites may be controlled by thermal processes in the early solar nebula. A detailed TEM analysis of additional CM hibonite samples is underway to evaluate this hypothesis.

Han, J.↗
NASA NTRS2005

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗