Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

Impact of electrode porosity architecture on electrochemical performances of 1 mm-thick LiFePO4 binder-free Li-ion electrodes fabricated by Spark Plasma Sintering

Thick electrodes with high active material loadings have been intensively studied over the last couple of decades in pursuit of achieving high energy density systems. To optimize and enhance the electrochemical performance of such electrodes, one has to control the pore morphology by, for example, varying the pore size and shape, and the level of porosity. In the present work, the fabrication of thick binder-free LiFePO4 (LFP) electrodes with two different pore sizes (12 and 20 mu m) and porosities (21 vol% and 44 vol%) using Spark Plasma Sintering (SPS) and templating approach is reported. The well-controlled porous architecture inside the thick electrodes is realized by fine-tuning experimental parameters. The impact of porosity architecture on electrochemical performance is quantified and correlated with the 3D tortuosity values determined from both micro-computed tomography and electrochemical impedance-based experimental methods. Based on the micro-computed tomography data analysis, estimated tortuosity values along X, Y, and Z axes reveal an anisotropic effect perpendicularly to the SPS compression axis (Z-direction). This is particularly profoundly observed in the samples with larger pores (20 mu m). The correlation between morphological properties and the rate capability performance is established indicating that the capacity loss happens mainly due to the Li-ion diffusion limitations.

Cathode material↗

Photoexcited energy relaxation and vibronic couplings in π-conjugated carbon nanorings

Conjugated carbon nanorings exhibit unique photophysical properties that, combined with their tunable sizes and conformations, make them suitable for a variety of practical applications. These properties are intimately associated to their strained, bent and sterically hindered cyclic structures. Herein we perform a comparative analysis of the photoinduced dynamics in carbon nanorings composed of nine phenyl units([9]CPP) and nine naphthyl units ([9]CN) respectively. The sterically demanding naphthyl units lead to large dihedral angles between neighboring units. Nevertheless, the ultrafast electronic and vibrational energy relaxation and redistribution is found to be similar for both systems. We observe that vibronic couplings, introduced by nonadiabatic energy transfer between electronic excited states, ensure the intramolecular vibrational energy redistribution through specific vibrational modes. The comparative impact of the internal conversion process on the exciton spatial localization and intra-ring migration indicates that naphthyl units in [9]CN achieve more efficient but less dynamical self-trapping compared to that of phenyl units in [9]CPP. That is, during the photoinduced process, the exciton in [9]CN is more static and localized than the exciton in [9]CPP. Furthermore, the internal conversion processes take place through a specific set of middle- to high-frequency normal modes, which directly influence the spatial exciton redistribution during the internal conversion, self-trapping and intra-ring migration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of Hybrid FPOG Applications in Regulated and Deregulated Markets Using HERON

Recent changes in the U.S. energy market, such as low natural gas prices and increased electricity production for variable renewable energy (VRE) sources, have led to an economic crisis for existing light-water reactor (LWR) nuclear power plants (NPP). Many owners and operators of LWRs have elected to decommission these plants rather than continue using them as consistent sources of clean baseload power. This has led to exploration of various possibilities to increase the economic viability of these units, including market restructuring to monetize benefits LWRs already provide to the grid through ancillary markets, load following and economic dispatch, and possible integration of secondary systems directly to the NPP for production of additional products through technologies such as hydrogen electrolysis or water desalination. Previous studies have considered the technologies associated with these Integrated Energy Systems (IES) activities, and the analysis of markets for these secondary products. To analyze the economic viability of various system configurations including IES, especially given the uncertainty surrounding load demand, electricity prices, and the availability of VRE resources, the stochastic technoeconomic analysis package HERON (Heuristic Energy Resource Optimization Network) was released earlier this year as an extension of the risk analysis framework RAVEN (Risk Analysis Virtual Environment). HERON focuses foremost on making the complex uncertainty quantification analysis tools approachable for energy systems analysts, also providing general dispatch optimization algorithms for those workflows. HERON continues to be improved and tested as a significant part of the IES viability analyses performed in this work. HERON is not a capacity expansion model. To consider market and grid energy system development in a variety of scenarios, HERON is best used in coupling with modelling tools such as US-REGEN, which sacrifice some of the uncertainty analysis and resolution of HERON's modelling for the ability to efficiently predict the change in the grid energy system's profile due to economic drivers over decades. HERON can then use this information to explore the economic viability of introducing changes to the predicted outcomes, such as the introduction of an IES. In this work, experts at EPRI using US-REGEN provide six projection scenarios for use in HERON stochastic technoeconomic analysis (STEA) in considering the options available for increasing LWR economic viability through introduction of a hydrogen-centric IES using a high-temperature steam electrolysis plant (HTSE), hydrogen storage, and a constant-rate contracted hydrogen consumer. The results obtained are differential in nature; they do not report expected profits for any configuration, but rather report on the possible increase in the NPV of a configuration with respect to a baseline no-IES configuration. Due to the uncertainty captured in the variable net load of the systems, there is likewise uncertainty in the mean values reported. We consider this viability both in terms of a regulated market, where the energy producers and IES are owned and operated by single entity, as well as a deregulated market, where the IES chooses its bid for electricity generation and is then dispatched by the grid system operator. Results indicate that for deregulated markets, the inclusion of the IES is often statistically beneficial. This is especially true in policies that are not favorable towards nuclear, as nuclear is less often dispatched and is forced to deal with frequent idle capacity. In the nominal case as well as the case of carbon tax policies, inclusion of the IES clearly benefited the economic performance of the NPP. In the regulated case, however, there was a trend towards minimizing the IES, likely due to the optimal sizing performed by US-REGEN of the NPP within the system as well as the lack of penalty for idle capacity at the NPP in the regulated market analyses.

99 GENERAL AND MISCELLANEOUS↗

Influence of Mid-Latitude Cyclones on European Background Surface Ozone Investigated in Observations, MACC and MERRA-2 Reanalyses

The relationship between springtime mid-latitude cyclones and background ozone at two rural monitoring sites on the west coast of Europe -- Mace Head, Ireland and Monte Velho, Portugal -- is explored using a combination of observations and three reanalyses: 1) the European Centre for Medium-Range Weather Forecasts' (ECMWF) ERA-Interim reanalysis, 2) the Monitoring Atmospheric Composition and Climate (MACC) reanalysis and 3) NASA's Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) reanalysis. The ERA-Interim cyclone tracks are used here to establish the long-term relationship between cyclones and ozone observations (since 1988). The MACC reanalysis data set, which covers the period 2003-2012, is produced with the ECMWF integrated forecast system (IFS) model two-way coupled to a chemistry transport model (CTM). Since the MACC reanalysis uses a similar atmospheric model to ERA-Interim, MACC is used to explore the mechanisms within the case study cyclones that can influence surface ozone concentrations at Mace Head and Monte Velho. The MERRA-2 reanalysis also provides 3D distributions of ozone, although less ideal for analysis of surface ozone concentrations since MERRA-2 ozone under represents ozone variability outside the stratosphere as it does not have a detailed chemistry scheme or emission sources for the troposphere. The MERRA-2 reanalysis, which has the potential to identify more features within the cyclones as the resolution is higher than the MACC reanalysis, is used in conjunction with the MACC reanalysis to provide a measure of uncertainty to the case study analysis. We found the main source of high ozone to these two sites is from the stratosphere, which is well represented in both the MERRA-2 and the MACC reanalyses, either from direct injection into the cyclone or associated with aged airstreams from decaying downstream cyclones that can become entrained and descend toward the surface within new cyclones over the NA region.

MERRA-2↗

Enabling residential heating decarbonization through hydronic low-temperature thermal distribution using forced-air assistive devices

Space heating represents approximately one-tenth of the United States’ energy use and has a breadth of potential for emission reduction. An element of space heating, hydronic heat distribution methods, use water supply temperatures up to 82.2 °C (180 °F). However, this operating temperature can be incompatible with high-efficiency heat generation systems, which typically provide heating temperatures of up to 60 °C (140 °F). In this work, a low-cost retrofit solution is developed using experimentally validated computational tools by incorporating an airflow distributor that preferentially directs the airflow from a fan to enhance heat transfer over the finned-tube heat exchanger found in conventional baseboards. The same heat output of traditional baseboard heat distribution systems operating at higher water supply temperatures (71.1 – 82.2 °C) can also be achieved at lower temperatures (≤60 °C). Specifically, results indicate that the technology can produce more than a 46.7% improvement in the heat transfer output at temperatures as low as 60 °C (140 °F), effectively matching the same output range achieved by 71.1 – 82.2 °C (160–180 °F) water supply temperatures. An important benefit of this solution is that it utilizes the existing infrastructure as opposed to requiring the replacement of the entire distribution system. Additionally, this technology enables existing building infrastructure to be coupled with newer, high-efficiency heating systems such as condensing boilers, solar-thermal systems, and geothermal/air-to-water heat pumps that may produce lower water supply temperatures (60 °C). A geospatially resolved techno-economic and environmental analysis is completed and presented to further understand the equivalent carbon footprint of enabling higher efficiency heat generation systems with the improved efficiency heat distribution system disclosed. Using 2021 grid-average emission factors, an emission reduction of up to a 67.5% decrease (2658 kgCO 2 /yr) for a single-family home, depending on state and climate region, could be realized by replacing a traditional natural gas-fired boiler with an air-to-water heat pump coupled to the same high-temperature heat distribution system along with the low-cost retrofit solution. A complete CO 2 emission reduction of 6688 kgCO 2 /yr, depending on state and climate region, could be realized if all electricity is further renewably sourced. Thus, this study provides a possible pathway towards enabling the reduction of operational emissions in space heating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mineralogical, Elemental and Tomographic Reconnaissance Investigation for CLPS (“METRIC”)

METRIC comprises a suite of two instruments in a Lunar lander that perform X-ray diffraction (XRD) for mineral structure, X-ray fluorescence (XRF) for elemental composition, and X-ray Micro Computed Tomography (XCT) for 3D internal micromorphology. The instruments are accompanied by optical-near IR cameras to provide local geologic context. The Honeybee Robotics PlanetVac pneumatic sampling and transfer system [1],positioned on a lander footpad, will deliver sieved regolith to the X-ray instruments for analysis. The instrument suite is intended for delivery to the lunar surface on a Commercial Lunar Payload Services (CLPS) lander. The METRICXRD/F instrument draws on heritage from the Mars Science Laboratory CheMin instrument [2] and improves upon the design in multiple ways [3]. Like CheMin, Rietveld refinement and full-pattern fitting of METRIC XRD data can identify minerals at a detection limit of ~1 wt.%, quantify their abundances when present at >3 wt.%, and determine mineral composition (e.g., Fo#in olivine) from lattice parameters for minerals present at >5 wt.%[4,5]. An optimized XRF geometry provides for improved detection and quantification of major, minor and trace elements. The METRIC XCT instrument is a miniaturized X-ray Computed Tomography scanner[6]. Lunar regolith is delivered to a 3 mm diameter, 10 mm long graphite tube inside the instrument. The tube is rotated through 360°in 0.9deg.incrementsand a divergent point source X-ray beam is directed through the material. A Charge Coupled Device (CCD)records attenuation images whose brightness and contrast area function of average atomic number and density. Quantitative data, including particle and void sizes, 3D particle shape parameters, modal volumes and pore geometry can be derived from the resulting 3D reconstructions(voxel resolution: 30 μm). Crystal morphologies derived from METRIC XCT data complement the bulk mineralogy determined by the METRIC XRD/F and provide a measure of grain size distribution for the different phases. Taken together, the METRIC instrument suite determines crystal structure, elemental composition and morphology, three principal characteristics of geological materials that are highly useful in determining the origin and subsequent processing of lunar regolith.

Moon↗

Seismoacoustic analysis of underground explosions using the Rayleigh integral

SUMMARY A coupled seismoacoustic model is developed for the analysis of acoustic signals produced by underground explosive events with an aim to develop a means of improving estimated depth and yield for explosion monitoring. A ground spall model is used to predict surface motion characteristics produced by an underground explosion and the Rayleigh integral is applied to relate the surface motion to the acoustic signal some distance from surface ground zero. The low-frequency component of the ground motion associated with the prolonged free fall of lofted material during spall is found to dominate the acoustic signal propagating away from surface ground zero at shallow angles. The model is applied to observed ground motion and acoustic signals recorded during the Source Physics Experiment (SPE) with promising results. In addition to accurately predicting characteristics of the observed acoustic signals from several of the SPE events, the model provides a means of explaining the lack of signals observed during several events in the SPE due to the directionality of the higher frequency acoustic signals associated with the uplift and closure components as well as the lack of a spatially localized, longer duration dwell.

58 GEOSCIENCES↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗

Structural insights on ligand recognition at the human leukotriene B4 receptor 1

The leukotriene B4 receptor 1 (BLT1) regulates the recruitment and chemotaxis of different cell types and plays a role in the pathophysiology of infectious, allergic, metabolic, and tumorigenic human diseases. Here we present a crystal structure of human BLT1 (hBLT1) in complex with a selective antagonist MK-D-046, developed for the treatment of type 2 diabetes and other inflammatory conditions. Comprehensive analysis of the structure and structure-activity relationship data, reinforced by site-directed mutagenesis and docking studies, reveals molecular determinants of ligand binding and selectivity toward different BLT receptor subtypes and across species. The structure helps to identify a putative membrane-buried ligand access channel as well as potential receptor binding modes of endogenous agonists. These structural insights of hBLT1 enrich our understanding of its ligand recognition and open up future avenues in structure-based drug design.

59 BASIC BIOLOGICAL SCIENCES↗

A Technology Pathway for Airbreathing, Combined-Cycle, Horizontal Space Launch Through SR-71 Based Trajectory Modeling

Access to space is in the early stages of commercialization. Private enterprises, mainly under direct or indirect subsidy by the government, have been making headway into the LEO launch systems infrastructure, of small-weight-class payloads of approximately 1000 lbs. These moderate gains have emboldened the launch industry and they are poised to move into the middle-weight class (roughly 5000 lbs). These commercially successful systems are based on relatively straightforward LOX-RP, two-stage, bi-propellant rocket technology developed by the government 40 years ago, accompanied by many technology improvements. In this paper we examine a known generic LOX-RP system with the focus on the booster stage (1st stage). The booster stage is then compared to modeled Rocket-Based and Turbine-Based Combined Cycle booster stages. The air-breathing propulsion stages are based on/or extrapolated from known performance parameters of ground tested RBCC (the Marquardt Ejector Ramjet) and TBCC (the SR-71/J-58 engine) data. Validated engine models using GECAT and SCCREAM are coupled with trajectory optimization and analysis in POST-II to explore viable launch scenarios using hypothetical aerospaceplane platform obeying the aerodynamic model of the SR-71. Finally, and assessment is made of the requisite research technology advances necessary for successful commercial and government adoption of combined-cycle engine systems for space access.

Kloesel, Kurt J.↗

Product-specific kinetics reveal effect of epoxy crystallization on thermoset thermal degradation

Crystallization is a common problem for epoxy resins, which are ubiquitous in industrial and commercial use. Integration of crystallized epoxy monomers into cured thermosets has been shown to alter the thermosets’ final mechanical properties. However, no studies have investigated the impact of these crystals on the thermal stability of the thermoset. Here we investigate the degradation kinetics of a bisphenol F–based epoxy thermoset with and without crystallized monomers using Product-Specific Kinetic (PSK) analysis coupled with Evolved Gas Analysis-Mass Spectrometry (EGA-MS). PSK analysis revealed significant differences in evolved product ion kinetics, suggesting changes in the degradation kinetics between thermoset configurations. Further, it was concluded that early stages of degradation are influenced most by crystal presence due to the high concentration of unreacted epoxy monomers and lower cross-linking density of the cured network. After post-cure annealing, significant changes are observed in the degradation kinetics of thermosets without crystal inclusions. Conversely, post-cure annealing procedures of crystal integrated thermosets showed little change in the thermoset degradation kinetics across all conversion extents. These findings suggest that post-cure annealing of thermosets with crystals present at the onset does not alter the cross-linking density of the polymer network enough to significantly change the degradation kinetics. We hypothesize this is because the excess monomers from the melted crystals are unable to find suitable reaction sites for complete binding into the polymer network, which has direct implications for the material properties and final thermal stability of the thermoset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic and Metallurgical Aspects of the Space Shuttle Columbia Accident Investigation and Reconstruction

The Space Shuttle Columbia was descending for a landing at the Kennedy Space Center (KSC) on February 1, 2003. Approximately 20 minutes prior to touchdown, the Columbia began disintegrating over the western United States; the majority of debris eventually impacted in eastern Texas and western Louisiana. A monumental effort eventually recovered approximately 84,000 pieces of debris, approximately 38% of the Orbiter's original dry weight. The debris was transported to KSC, where the items were catalogued and evaluated. Critical areas of interest, such as the left and right leading edge surfaces and the underside of the ship, were placed upon a grid to aid in the reconstruction. Items of interest included metallic structures, reinforced carbon-carbon composites, and ceramic heat insulation tiles. Many of the leading edge elements had re-solidified metallic deposits spattered on them. These deposits became known as slag and were one of the main focuses of the investigation. In order to help determine the sequence of events inside the left wing during the accident, the slag's composition, layering order, and directionality of deposition were studied. A myriad of analytical tests were performed in an attempt to ascertain the compositional and depositional characteristics of selected slag deposits, including the ordering of deposited layers within each individual slag deposit harvested. Initially, Scanning Electron Microscopy and Energy Dispersive X-Ray Spectroscopy (SEM/EDX) were performed to quickly characterize the overall composition of individual slag deposits: SEM utilizes a narrowlyfocused high-energy electron beam impinging upon a specimen. The incident beam excites and liberates lower energy secondary electrons, which are detected and analyzed, providing a visual representation of the sample's surface topography. EDX also relies on an incident electron beam, except an EDX unit measures X-ray energies generated by the impinging beam. Each element generates a unique X-ray signature; the EDX detector measures these discreet energies. EDX actually penetrates approximately 2 microns into the bulk of the sample. However, random examination of various portions of slag, coupled with the semiquantitative nature of the SEM/EDX analysis, did not yield convincingly pertinent data. Therefore, X-ray dot mapping was conducted, which provided more understandable data, both in terms of slag layering and composition. An X-ray dot map is generated by performing numerous EDX scans for individual elements, then compiling the scans in a visual representation. Eventually, specimens consisting of not only the slag, but of the adjacent RCC substrate as well were cross-sectioned. X-Ray dot mapping of the materialographicallymounted and -polished cross- sections provided a visual representation of both the layering sequence and compositional characteristics of the slag. Contemporaneously, Electron Spectroscopy for Chemical Analysis/X-Ray Photoelectron Spectroscopy (ESCA/XPS) and powdered X-Ray Diffraction (XRD) were performed to further characterize the deposits and to attempt to identify what, if any, compounds were present. The ESCA/XPS analysis allowed the analyst to "sputter" into the sample with an electron gun, aiding in the identification of the layering sequence. XPS uses photons, rather than electrons, which impinge upon the surface of the sample. XPS measures the electrons emitted from within the first 5 nm of the sample's surface. The XRD measures the scatter angles of incident X-rays; the angle and intensity of scatter depend upon the crystalline structure of the pulverized sample. XRD is considered a qualitative rather than quantitative technique. ESCA/XPS revealed that the final layer to deposit was predominantly carbonaceous. XRD was successful in identifying specific compounds, such as Al 2O3, Al and/or Al3 21SiO47, mullite (3(Al2)O3 -SiO2), and nickel-aluminides. Eventually, Electron MicroProbe Analysis (EMPA) was conducted on the marialographically-prepared cross- sections of selected slag deposits. Microprobe combines SEM and Wavelength Dispersive X-Ray Spectroscopy (WDS), and, like EDX, uses a narrowly-focused high-energy electron beam impinging upon a specimen to elicit, in the case of EPMA, characteristic X-rays with specific wavelengths. This quantitative, analytical tool proved the most useful in determining depositional layering and composition of the slag deposits. This information was utilized in verifying the location of the breach in the left leading edge of the wing of the Columbia.

McDaniels, Steven J.↗

Infrared Spectroscopy as a Chemical Fingerprinting Tool

Infrared (IR) spectroscopy is a powerful analytical tool in the chemical fingerprinting of materials. Any sample material that will interact with infrared light produces a spectrum and, although normally associated with organic materials, inorganic compounds may also be infrared active. The technique is rapid, reproducible and usually non-invasive to the sample. That it is non-invasive allows for additional characterization of the original material using other analytical techniques including thermal analysis and RAMAN spectroscopic techniques. With the appropriate accessories, the technique can be used to examine samples in liquid, solid or gas phase. Both aqueous and non-aqueous free-flowing solutions can be analyzed, as can viscous liquids such as heavy oils and greases. Solid samples of varying sizes and shapes may also be examined and with the addition of microscopic IR (microspectroscopy) capabilities, minute materials such as single fibers and threads may be analyzed. With the addition of appropriate software, microspectroscopy can be used for automated discrete point or compositional surface area mapping, with the latter providing a means to record changes in the chemical composition of a material surface over a defined area. Due to the ability to characterize gaseous samples, IR spectroscopy can also be coupled with thermal processes such as thermogravimetric (TG) analyses to provide both thermal and chemical data in a single run. In this configuration, solids (or liquids) heated in a TG analyzer undergo decomposition, with the evolving gases directed into the IR spectrometer. Thus, information is provided on the thermal properties of a material and the order in which its chemical constituents are broken down during incremental heating. Specific examples of these varied applications will be cited, with data interpretation and method limitations further discussed.

Huff, Timothy L.↗

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗

What the Long-Term Sunspot Record Tells Us About Space Climate

Direct observations of sunspots span the nearly 400 years since the time of Galileo. Dedicated observing programs at several observatories over the last 150 years have provided detailed information not only on the number of sunspots but on their sizes and positions as well. The data acquired by those original observers, and by those who have more recently brought those observations to light, provide important clues about the nature of the solar cycle and its contribution to space climate. The period of the cycle, the equator-ward drift of the active latitudes, the asymmetry between the rise to maximum and the fill to minimum, shifting asymmetries between northern and southern hemisphere activity, the tilt of active regions, and the increasing amplitude of the cycles since the Maunder Minimum are all well established. Other, less well established characteristics such as multi-cycle and short-term periodicities, often depend upon the method of data analysis. The strong correlation between sunspot statistics and other measures of solar activity, coupled with the length of the sunspot record, make these observations extremely valuable for characterizing and understanding space climate.

Hathaway, D. H.↗

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗