Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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.

572 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2022

Local magnetic moments due to loop currents in metals

We present Hartree-Fock calculations on a simple model to obtain the conditions of formation of local magnetic moments due to loop currents $L$ $o$ and spin-loop currents $L$ $s$ and compare them to the conditions of formation of local spin moments $M$ which were given long ago in a similar approximation by Anderson. A model with three degenerate orbitals sitting on an equilateral triangle, with on-site and nearest-neighbor repulsions $U$ and $V$, respectively, and intersite kinetic energy, hybridizing with conduction electrons with a parameter $Δ$ is investigated. $L$ $o$ and $L$ $s$ are promoted by large $V$/$Δ$ and their magnitude is relatively unaffected by $U$/$Δ$. Spin-magnetic moments $M$ promoted by large $U$/$Δ$ on the other hand are adversely affected by $V$/$Δ$. In this model, $L$ $o$ for $V$ multiplied by the number of neighbors is approximately the same as the $M$ promoted by $U$ in Anderson’s local model for $M$. $L$ $o$ and $L$ $s$ are degenerate if exchange interactions and Hund’s rule are neglected but $L$ $o$ is favored when they are included. Many of the qualitative results are visible in an expression for the Hartree-Fock ground state energy derived as a function of small $L$ $o$ , $L$ $s$ , and $M$. Numerical minimization of the Hartree-Fock energy is presented for larger values. Additionally, we also briefly discuss the connection and differences of the interaction-generated orbital currents and spin currents discussed here and generalized to a lattice with the topological states in metals and semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nanoporous weakly coordinating anionic frameworks (Final Technical Report)

Under this award, we synthesized and characterized new weakly coordinating, anionic covalent organic frameworks. The motivation was to expand the class of covalent organic frameworks though the introduction of the first weakly coordinating anionic ones. The significance of the weakly coordinating anionic character of these frameworks is associated with the character of the counter (metal) cations in their pores. Those counter-cations, due to the inability to effectively coordinate to the pore walls, have a “naked”, and thus highly electrophilic, Lewis acidic character. Such properties are expected to translate to a number of interesting physical properties that can be potentially used for a range of energy-relevant applications. One of these applications is gas adsorption, separation, and storage. The electrophilic, Lewis acidic character of the counter-cations are expected to strongly bind to gas molecules entering the pores leading to higher gas sorption capacities. Other applications include heterogeneous catalysis, ion exchange, and solid state ion conduction. Within the project we explored a range of chemical blocks, and tested synthesis conditions to polymerize these building blocks to form the target materials. In particular, we focused on materials with group 15 elements like phosphorus in the oxidation state +5 and the coordination number 6, that have a formal -1 charge. The synthesis of such materials is challenging because its success depends on the absence of residual monomeric and oligomeric side-products in the pores, and the absence of significant framework interpenetration. Crystalline forms of these materials require sufficient reversibility of the polymerization reaction (which is difficult to achieve because group 15 elements like P tend to make strong covalent bonds with neighbored elements), as well as suitable polymerization rates that prevent the formation of kinetically controlled products. We had success with the preparation of anionic porous organic frameworks with anionic building blocks containing P(V)O6 octahedra synthesized through Yamamoto coupling reactions. These materials are non-crystalline because of the irreversible nature of the Yamamoto coupling reaction, but show porosity. The porosity and surface area can be significantly increased by co-reaction with electrostatically neutral co-reactants. Very high porosity and surface areas can be achieved when the materials are washed with hydrochloric acid. This method is commonly used to free products from Yamamoto coupling reactions from by-products. Indeed, this method removed by products from the pores of our materials, however it also led to the hydrolysis of the PO6 octahedra to form pending phosphonic acid groups (-PO3H)-groups in the frameworks (PA-POFs) which was actually not desired. Nonetheless, these materials showed high CO2 sorption capacities. In addition, they show extremely high (actually record) sorption capacities for the removal of a range of common micropollutants such as bisphenol A (BPA) and 4-nitrophenol in water. BPA is a wide-spread organic pollutant and a known endocrine disruptor. The maximum adsorption capacity of BPA at equilibrium is found to be as high as 3,366 mg g-1 by Langmuir adsorption model, which is more than 10 times greater than peer materials. The polymer also rapidly removes various other organic micropollutants with more than 90% removal efficiencies. In addition, the PA-POF material can be regenerated at least five times by mild washing using methanol without significant loss in removal efficiency. This is significant advantage over commonly used activated carbon materials (that are in Brita water faucet filters, for example). Comparison studies showed that only materials with the phosphonate groups show extreme sorption although the number of PA group is actually small according to solid state NMR and EDX data. More research would be needed to investigate why this is (a possibility is that the phosphonate groups have a multiplying effect in which a PA-induced initial adsorption leads to the adsorption of much more adsorbate). The extreme sorption capabilities can also be visually seen when colored adsorbates such as 4-nitrophenol are used, as the adsorption of the nitrobenzene leads to the decoloration of the solution. We have produced three short videos that illustrate the sorption capabilities compared to two peer materials (activated carbon from dismantled Brita water filters, and a electrostatically neutral, chemically related, porous organic polymer, respectively).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Overcoming timestep limitations in boosted-frame particle-in-cell simulations of plasma-based acceleration

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Exploration of relationships between low-Ti and high-Ti pristine lunar glasses using an armalcolite assimilation model

The pristine glasses of Delano are the most primitive lunar basaltic magma compositions discovered to date. They are grouped into two (and possibly three) arrays: a low-alumina array and a high alumina array. These glasses are very olivine normative and are multiply saturated at pressures of approximately 20 kbar, implying a depth of origin of 400 to 500 km in the Moon. Thus, these glasses appear to be the best candidates for primitive partial melts of the upper lunar mantle. One of the most perplexing characteristics of the pristine glasses is a positive correlation between Ni and SiO2 within each array. This is contrary to the terrestrial experience, where Ni is observed to positively correlate with MgO and negatively correlate with SiO2. These systematics are believed to be due to the depletion of Ni by olivine fractionation. The difference between the lunar and terrestrial Ni vs. SiO2 trends may be partially ascribed to the Ti-rich component. In the case of the pristine glasses, SiO2 increases not because of olivine fractionation, but because they contain less of the high-Ti component. An attempt was made to model this variation in Ni and SiO2 with a simple assimilation-fractional crystallization (AFC) model. Silica and Ni both decreased dramatically as the AFC process proceeded. Only 15 to 20 percent AFC was necessary to produce the observed variation, and the SiO2 vs. Ni variation was modeled quite well. The D(Ni) for olivine/liquid in this model was taken to be 10 and the olivine was assumed to be Fe sub 80. However, the results of this model for Ti and Mg were less than satisfactory. It seemed difficult to achieve the high TiO2 contents of some glasses (16 to 17 wt. percent) by this method. Continual addition of ilmenite by AFC could indeed raise the titania concentrations to the necessary levels, but only by enriching the magma in FeO and greatly depleting the magma in MgO. An attempt was made to circumvent this problem by using armalcolite, (Fe, Mg)Ti2O5, in the AFC model, and the results are presented.

Jones, John H.↗
NASA NTRS2009

A First Look at Carbon and Oxygen Stable Isotope Measurements of Martian Atmospheric C02 by the Phoenix Lander

Precise stable isotope measurements of the CO2 in the martian atmosphere have the potential to provide important constraints for our understanding of the history of volatiles, the carbon cycle, current atmospheric processes, and the degree of water/rock interaction on Mars. The isotopic composition of the martian atmosphere has been measured using a number of different methods (Table 1), however a precise value (<1%) has yet to be achieved. Given the elevated 13C values measured in carbonates in martian meteorites it has been supposed that the martian atmosphere was enriched in delta(sup 13)C. This was supported by measurements of trapped CO2 gas in EETA 79001[2] which showed elevated delta(sup 13)C values (Table 1). More recently, Earth-based spectroscopic measurements of the martian atmosphere have measured the martian CO2 to be depleted in delta(sup 13)C relative to CO2 in the terrestrial atmosphere. The spectroscopic measurements performed by Krasnopolsky et al. were reported with approx.2% uncertainties which are much smaller than the Viking measurements, but still remain very large in comparison to the magnitude of carbon and oxygen isotope fractionations under martian surface conditions. The Thermal Evolved Gas Analyzer (TEGA) instrument on the Mars Phoenix Lander included a magnetic sector mass spectrometer (EGA) which had the goal of measuring the isotopic composition of martian atmospheric CO2 to within 0.5%. The mass spectrometer is a miniature magnetic sector instrument intended to measure both the martian atmosphere as well as gases evolved from heating martian soils. Ions produced in the ion source are drawn out by a high voltage and focused by a magnetic field onto a set of collector slits. Four specific trajectories are selected to cover the mass ranges, 0.7 - 4, 7 - 35, 14 - 70, and 28 - 140 Da. Using four channels reduces the magnitude of the mass scan and provides simultaneous coverage of the mass ranges. Channel electron multiplier (CEM) detectors that operate in the pulse counting mode detect the ion beams.

Niles, P.B.↗
DOE OSTI2023

Study of nanodiamond photocathodes for MPGD-based detectors of single photons

The proposed new Electron–Ion Collider poses a technical and intellectual challenge for the detector design to accommodate the long-term diverse physics goals envisaged by the program. This requires a 4π detector system capable of reconstructing the energy and momentum of final state particles with high precision. The Electron-Ion Collider also requires identification of particles of different masses over a wide momentum range. A diverse spectrum of Particle IDentification detectors has been proposed. Of the four types of detectors for hadron identification, three are based on Ring Imaging Cherenkov Counter technologies, and one is realized by the Time of Flight method. The quest for a novel photocathode, sensitive in the far vacuum ultraviolet wavelength range and more robust than cesium iodide, motivated an R&D programme to explore nano-diamond (ND) based photocathodes, started by a collaboration between INFN and CNR Bari and INFN Trieste. Systematic measurements of the photoemission in different Ar:CH 4 and Ar:CO 2 gas mixtures with various types of ND powders and Hydrogenated ND (H-ND) powders are reported. A first study of the response of THGEMs coated with different photocathode materials is presented. Additionally, the progress of this R&D programme and the results obtained so far by these exploratory studies are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2005

Apparatus for Testing Flat Specimens of Thermal Insulation

An apparatus has been developed to implement an improved method of testing flat-plate specimens of thermal-insulation materials for cryogenic application. The method includes testing under realistic use conditions that could include vacuum and mechanical loading at a pressure up to 70 psi (=0.48 MPa). The apparatus can accommodate a rigid or flexible specimen having thickness up to 1.25 in. (=3.2 cm) and diameters between 6 and 10 in. (about 15.2 and 25.4 cm, respectively). Typical test conditions include boundary temperatures between 77 K and 373 K and vacuum/interstitial gas filling at a pressure between 10(exp -6) torr (=1.3 x 10(exp -4) Pa) and 760 torr (atmospheric pressure =0.1 MPa). The interstitial gas could be N2, He, CO2, or any other suitable gas to which the insulation is expected to be exposed in use. Relative to prior apparatuses and testing methods, this apparatus and the testing method that it implements offer advantages of relative simplicity and ease of use. The basic principle of operation of the apparatus is that of boil-off calorimetry, using liquid nitrogen or any other suitable liquid that boils at a desired temperature below ambient temperature. Comparative rates of flow of heat through the thicknesses of the specimens (heat-leak rates) and apparent-thermal-conductivity values are obtained from tests of specimens. Absolute values of heat-leak rates and apparent thermal conductivities are computed from a combination of (1) the aforementioned comparative values and (2) calibration factors obtained by testing reference specimens of materials that have known thermal-insulation properties. The apparatus includes a full complement of temperature sensors, a vacuum pump and chamber, a monitoring and control system, and tools and fixtures that enable rapid and reliable installation and removal of specimens. A specimen is installed at the bottom of the vacuum chamber, and a cold-mass assembly that includes a tank is lowered into position above and around the specimen (see figure). A spring-based compensating fixture helps to ensure adequate thermal contact with possibly irregular specimen surfaces. For a high-compression test, the springs can be replaced with spacers. A flat circular load cell at the bottom of the chamber measures the compressive load on the specimen. Once the desired compressive-load, temperature, and vacuum/gas-filling conditions are established, testing begins. During a test, all measurements are recorded by use of a portable data-acquisition system and a computer. The total heat-leak rate is measured and calculated as the boil-off flow rate multiplied by the latent heat of vaporization. The parasitic heat leak (to the side of the specimen and to the top and side of the cold-mass tank) is reduced to a small fraction of the total heat leak by use of a combination of multilayer-insulation (MLI) shield rings, reflective film, a fiberglass/epoxy centering ring, and a bulk fill of aerogel beads. This combination eliminates the need for a cryogenic guard chamber used in a typical prior apparatus to reduce the parasitic heat leak.

Fesmire, James E.↗
DOE OSTI2022

Investigating the AGN-201M research reactor's unique dominance ratio

The AGN-201M reactor is a small, thermal spectrum reactor located at the University of New Mexico (UNM). It is moderated by polyethylene, reflected by graphite, with fuel comprised of uranium microspheres embedded in polyethylene plates that are separated by an aluminum baffle. The dominance ratio of a multiplying system is the ratio of the first higher-mode eigenvalue to the fundamental, k{sub 1}/k{sub 0}. It is a key parameter in characterizing the rate of convergence, analyzing system stability, and determining how tightly coupled the system's neutron fission is. The purpose of this paper is to investigate the system's unusually low dominance ratio of 0.632. This investigation varied the reactor's physical parameters, such as the fuel geometry, fuel density, and reflector thickness, and examined their impact on the reactor's dominance ratio. Additionally, this study examined neutronically similar systems in hopes of identifying common causes for systems with low dominance ratios. The reason why the dominance ratio of the reactor is small compared to larger thermal reactors was determined to be because of its significantly smaller geometry and the composition of the fuel plates. The reflector's effect on the dominance ratio is small in comparison to the other factors, but was found to have a non-zero effect. Furthermore, the AGN-201M was found to have a significantly lower dominance ratio than systems in which it shares a very high (c{sub k} >95%) degree of neutronic similarity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Dendrite-resolved, full-melt-pool phase-field simulations to reveal non-steady-state effects and to test an approximate model

In this report we study the epitaxial, columnar growth of (multiply oriented) dendrites/cells for a spot melt in a polycrystalline Al–Cu substrate using two-dimensional, phase-field, direct numerical simulations (DNS) at the full-melt-pool scale. Our main objective is to compare the expensive DNS model to a much cheaper but approximate “line” model in which a single-crystal phase-field simulation is confined to a narrow rectangular geometry. To perform this comparison, we develop algorithms that automatically extract quantities of interest (QoIs) from both DNS and line models. These QoIs allow us to quantitatively assess the assumptions in the line model and help us analyze its discrepancy with the DNS model. We consider four sets of heat source parameters, mimicking welding and additive manufacturing conditions, that create a combination shallow and deep melt pools. Our largest DNS simulation used 16K x 14K grid points in space. Our main findings can be summarized as follows. Under AM conditions, the QoIs of line models are in excellent agreement with the full DNS results for both shallow and deep melt pools. Under welding conditions, the primary spacing of the DNS model is smaller than the prediction of line model. We identify a geometric crowding effect that accounts for the discrepancies between the DNS and line models. We propose two potential mechanisms that determine the response of the microstructure to geometric crowding.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High Thermoelectric Performance and Defect Energetics of Multipocketed Full Heusler Compounds

We report a first-principles density-functional study of electron-phonon interactions in and thermoelectric transport properties of the full Heusler compounds Sr 2 BiAu and Sr 2 SbAu. Our results show that ultrahigh intrinsic bulk thermoelectric performance across a wide range of temperatures is physically possible and point to the presence of multiply degenerate and highly dispersive carrier pockets as the key factor for achieving this. Sr 2 BiAu, which features ten energy-aligned low-effective-mass pockets (six along Γ-X and four at L), is predicted to deliver n-type zT=0.4-4.9 at T=100-700 K. Comparison with the previously investigated compound Ba 2 BiAu shows that the additional L pockets in Sr 2 BiAu significantly increase its low-temperature power factor to a maximum value of 12 mW m -1 K -2 near T=300 K. However, at high temperatures the power factor of Sr 2 BiAu drops below that of Ba 2 BiAu because the L states are heavier and subject to strong scattering by phonon deformation, as opposed to the lighter Γ-X states, which are limited by polar-optical scattering. Sr2SbAu is predicted to deliver a lower n-type zT=3.4 at T=750 K due to appreciable misalignment between the L and Γ-X carrier pockets, generally heavier scattering, and a slightly higher lattice thermal conductivity. Soft acoustic modes, which are responsible for the low lattice thermal conductivity, also increase the vibrational entropy and high-temperature stability of these Heusler compounds, suggesting that their experimental synthesis may be feasible. The dominant intrinsic defects are found to be Au vacancies, which drive the Fermi level towards the conduction band and work in favor of n-doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Accuracy and precision of ESI-Orbitrap-IRMS observations of hours to tens of hours via reservoir injection

Here, orbitrap isotope ratio mass spectrometry (Orbitrap-IRMS) has recently been applied to high-precision, natural-abundance isotope ratio measurements of a diverse range of compounds, including amino acids, oxyanions, fatty acids, and metals. These measurements can characterize many isotope ratios simultaneously at high (≈1.0‰) precision. In a successful experiment, observed precision will track the shot-noise limit and be limited by experimental time. Some isotope ratios, for example those involving 17 O in organic compounds or multiply-substituted (‘clumped’) isotopologues, require experimental times of hours to tens of hours to achieve desired precision, while current sample introduction techniques focus on observations on the order of seconds to tens of minutes. In this study, we characterize Orbitrap-IRMS performance for three long duration measurements (individual acquisitions ≥1 h and as long as 24 h) using an automated reservoir injection system coupled to a Q Exactive HF Orbitrap with an electrospray ionization (ESI) source. First, we characterize long-term intra-measurement stability through a 24-h long measurement of acetone. We report the following isotope ratios and precisions (as acquisition errors, errors on the observed ratio within this measurement (σ AE )): 13 C/ 12 C (σ AE = 0.07‰), 17 O/ 16 O (σ AE = 1.1‰), 18 O/ 16 O (σ AE = 0.3‰), and 13 C 13 C/ 12 C (σ AE = 0.65‰). The σ AE of each tracks the shot noise limit throughout and is limited by the challenging conditions (high resolution and low numbers of ions per scan) required for 17 O/ 16 O measurement in the presence of 13 C via Orbitrap. Second, we characterize inter-measurement stability via a sequence of seven 75-min analyses of perchlorate. We observe the following ratios and acquisition errors: 37 Cl/ 35 Cl (σ AE = 0.09‰); 17 O/ 16 O (σ AE = 1.6‰); 18 O/ 16 O (σ AE = 0.7‰), 37 Cl 17 O/ 35 Cl 16 O (σ AE = 2.7‰), and 37 Cl 18 O/ 35 Cl 16 O (σ AE = 1.2‰). However, we find that inter-measurement drift between acquisitions limits our accuracy and precision for standardized measurements (i.e., error on reported δ values) to ≈1‰ for the 37 Cl/ 35 Cl measurement. Hence, the benefits of low σ AE may not be fully realized. Third, we demonstrate accuracy via sample/standard comparisons of a methionine sample with 13 C enrichment of ≈20‰ relative to a known standard. Using a sequence of seven 60-min analyses, we recover the following isotope ratios and standardized precisions (i.e., error on reported δ values, denoted propagated acquisition errors, σ PAE ): 33 S/ 32 S (σ PAE = 1.0‰), 34 S/ 32 S (σ PAE = 0.7‰), 15 N/ 14 N (σ PAE = 2.1‰), 2 H/ 1 H (σ PAE = 3.2‰), 13 C/ 12 C (σ PAE = 0.4‰), 18 O/ 16 O (σ PAE = 1.6‰), & 13 C 13 C/ 12 C (σ PAE = 2.8‰) with confirmation of accurate results for the known 13 C/ 12 C and 13 C 13 C/ 12 C enrichments. Together, our results demonstrate the viability of Orbitrap-IRMS for long duration measurements of diverse sample types via an automated reservoir injection system. Inter-measurement stability remains a challenge; we expect our methods to be most applicable to extended measurements of hard-to-observe properties, such as 17 O in organics and clumped isotopologues.

58 GEOSCIENCES↗
NASA NTRS2005

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗