Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Outer approximation”

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 1,045 records · Page 58

Locations of termination shock and heliopause based on Voyager plasma and magnetic field data

The locations of the termination shock and the heliopause are studied taking into account the effects of pickup protons. The study uses available plasma and magnetic field data from Voyagers over a 14-year period (1978-1991) and Voyager observation of the 1992-93 radio emission event. Outside 30 AU, pickup protons have a significant influence on dynamical structures of the outer heliosphere. The solar wind is treated as a mixture of electrons, solar wind protons, and interstellar pickup protons. If the magnitude of the interstellar magnetic field B(sub int) is given, one can quantitatively study the motion and location of the termination shock. The location is anti-correlated with the sun spot number and the shock has an average speed of approx. 24 km/s. Because B(sub int) is poorly known, additional information is needed in studying the termination shock. Cummings, et al. have used observations of anomalous cosmic rays to estimate the location of the shock. The observations of the 1991 GMIR and GMIR shock and the 1992-93 radio emission event provide another handle for the study of the termination shock and the heliopause. After its penetration through the termination shock, the GMIR shock continued to propagate in the subsonic region of the solar wind and eventually interacted with the heliopause. This interaction produces a transmitted shock propagating outward in the interstellar medium and a reflected shock propagating inward toward the sun in the subsonic solar wind. The plasma frequencies behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. Taking into account the effects of pickup protons we found that the average locations of the termination shock and the heliopause in 1991-92 are at approximately 66 AU and 150 AU, respectively.

Whang, Y. C.↗

Characterization of Meteorites by Focused Ion Beam Sectioning: Recent Applications to CAIs and Primitive Meteorite Matrices

Focused ion beam (FIB) sectioning has revolutionized preparation of meteorite samples for characterization by analytical transmission electron microscopy (TEM) and other techniques. Although FIB is not "non-destructive" in the purest sense, each extracted section amounts to no more than nanograms (approximately 500 cubic microns) removed intact from locations precisely controlled by SEM imaging and analysis. Physical alteration of surrounding material by ion damage, fracture or sputter contamination effects is localized to within a few micrometers around the lift-out point. This leaves adjacent material intact for coordinate geochemical analysis by SIMS, microdrill extraction/TIMS and other techniques. After lift out, FIB sections can be quantitatively analyzed by electron microprobe prior to final thinning, synchrotron x-ray techniques, and by the full range of state-of-the-art analytical field-emission scanning transmission electron microscope (FE-STEM) techniques once thinning is complete. Multiple meteorite studies supported by FIB/FE-STEM are currently underway at NASA-JSC, including coordinated analysis of refractory phase assemblages in CAIs and fine-grained matrices in carbonaceous chondrites. FIB sectioning of CAIs has uncovered epitaxial and other overgrowth relations between corundum-hibonite-spinel consistent with hibonite preceding corundum and/or spinel in non-equilibrium condensation sequences at combinations of higher gas pressures, dust-gas enrichments or significant nebular transport. For all of these cases, the ability of FIB to allow for coordination with spatially-associated isotopic data by SIMS provides immense value for constraining the formation scenarios of the particular CAI assemblage. For carbonaceous chondrites matrix material, FIB has allowed us to obtain intact continuous sections of the immediate outer surface of Murchison (CM2) after it has been experimentally ion processed to simulate solar wind space weathering. The surface amorphization and loss of OH produced by the irradiation provides important clues regarding space weathering on primitive asteroids such as the OSIRIS-Rex target 101955 Bennu.

Christoffersen, Roy↗

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗

Experimental Determination of the Recovery Factor and Analytical Solution of the Conical Flow Field for a 20 deg Included Angle Cone at Mach Numbers of 4.6 and 6.0 and Stagnation Temperatures to 2600 degree R

The local recovery factor was determined experimentally along the surface of a thin-walled 20 deg included angle cone for Mach numbers near 6.0 at stagnation temperatures between 1200 deg R and 2600 deg R. In addition, a similar cone configuration was tested at Mach numbers near 4.5 at stagnation temperatures of approximately 612 deg R. The local Reynolds number based on flow properties at the edge of the boundary layer ranged between 0.1 x 10(exp 4) and 3.5 x 10(exp 4) for tests at temperatures above 1200 deg R and between 6 x 10(exp 4) and 25 x 10(exp 4) for tests at temperatures near 612 deg R. The results indicated, generally, that the recovery factor can be predicted satisfactorily using the square root of the Prandtl number. No conclusion could be made as to the necessity of evaluating the Prandtl number at a reference temperature given by an empirical equation, as opposed to evaluating the Prandtl number at the wall temperature or static temperature of the gas at the cone surface. For the tests at temperatures above 1200 deg R (indicated herein as the tests conducted in the slip-flow region), two definite trends in the recovery data were observed - one of increasing recovery factor with decreasing stagnation pressure, which was associated with slip-flow effects and one of decreasing recovery factor with increasing temperature. The true cause of the latter trend could not be ascertained, but it was shown that this trend was not appreciably altered by the sources of error of the magnitude considered herein. The real-gas equations of state were used to determine accurately the local stream properties at the outer edge of the boundary layer of the cone. Included in the report, therefore, is a general solution for the conical flow of a real gas using the Beattie-Bridgeman equation of state. The largest effect of temperature was seen to be in the terms which were dependent upon the internal energy of the gas. The pressure and hence the pressure drag terms were unaffected.

Pfyl, Frank A.↗

Fabrication and Characterization of Core–Shell Density-Graded 316L Stainless Steel Porous Structure

In the present work, cylindrical shape 316 L cellular structure with different densities in its core and outer layer (shell) is fabricated by using carbamide as a space holder, via powder metallurgy route and layer-by-layer technique. The arrangement of the created pore is the same as face-centered cubic atomic structure, and different densities are created in two regions by using carbamide particles having two different sizes in the range of 1.7-2 and 2-2.4 mm. The effect of creating a structure with higher porosity (64.5%) in the core and lower porosity (53.8%) in the shell and vice versa and also change in the ratio of the core to the cylinder cross-sectional areas, on the mechanical properties and compaction load bearing of the fabricated foam samples were investigated. In leaching process of the carbamide particles, as an important step of porous structure’s fabrication, it is shown that discontinuous leaching process is more favorable than continuous, by which it would be possible to remove more carbamide particles (around 20%) at the same time. Furthermore, the deformation of the density-graded foam shows the parallel mechanism in the core and shell sections, and the contribution of each part depends on its density and thickness. The energy absorption behavior of the fabricated specimens is evaluated optimally in terms of the energy absorption value associated with the ideal adsorption behavior. The maximum ideal energy absorption efficiency for the samples with more porosity in the shell was approximately equal to 0.91, while for the sample includes the lower amount of porosity in the shell, this value was in the range of 0.85 and 0.89.

36 MATERIALS SCIENCE↗

Large cosmic ray transient decreases observed in the heliosphere at the Earth, Voyagers 1 and 2 and Pioneer 10 from 1987 to 1991

Using the greater than 70-MeV cosmic ray rate data from the Interplanetary Monitoring Platform (IMP), Voyager, and Pioneer spacecraft, we have examined the role of large transient decreases in the overall cosmic ray decrease from 1987 to 1990. At least 12 separate transient decreases with magnitudes greater than 5% were observed at Voyager 2, which moved from a heliocentric radius of approximately 23 to 35 AU during the time period. Many, but not all, of these decreases were observed at IMP 8 near the Earth and at the Pioneer 10 between 42 and 54 AU near the ecliptic plane. However, only five of these decreases were observed at Voyager 1 at north heliographic latitude approximately 30 deg and radial distance between 32 and 46 AU. The period of decreasing cosmic ray intensity from 1987 to 1990 can be divided into two phases. During phase 1, between mid-1987 and early 1989, the average heliospheric current sheet tilt rapidly increased from approximately 8 deg to 65 deg as the solar activity was also increasing rapidly and the intensity at the Earth decreased by about 35%. Four large transient decreases were observed at the Earth during this time; however, transient decreases beyond the Earth were observed only near the ecliptic plane. In phase 2, from early 1989 to mid-1990 when the intensity at the Earth decreased by an additional approximately 60%, five large transient decreases were observed at all latitudes. During this time the current sheet tilt was always greater than 60 deg. Using a procedure to separate these large transient decreases from any possible slower long-term variations, we argue that these decreases appear to be responsible for much of the overall intensity decrease of 20 - 35% observed at all spacecraft near the equatorial plane in phase 1 as well as most of the larger intensity decrease of about 45 - 60% observed at all spacecraft during phase 2. The recovery of the cosmic ray intensity that started at the Earth in early 1990 moved outward in the heliosphere reaching the Voyager and Pioneer spacecraft about 6 months later. This is in contrast to what happened as the intensity began to recover was observed to begin almost simultaneously at all radii less than 25 AU. This difference suggests that during the current cycle, the mechanism responsible for the recovery must have propagated outward from the Sun. This could have been related to the reversal of the solar polar magnetic field which occurred in early 1990. The effects of large transients propagating outward in the heliosphere and becoming ineffective when they reach the heliospheric boundary could not have played a major role in the onset of the recovery in 1990, since in this case the recovery should have been observed first in the outer heliosphere.

Webber, W. R.↗

Advanced Burner Reactor with Breed-and-Burn Thorium Blankets for Improved Economics and Resource Utilization

This study assesses the feasibility of designing Seed and Blanket (S&B) Sodium-cooled Fast Reactor (SFR) to generate a significant fraction of the core power from radial thorium fueled blankets that operate on the Breed-and-Burn (B&B) mode without exceeding the radiation damage constraint of presently verified cladding materials. The S&B core is designed to maximize the fraction of neutrons that radially leak from the seed (or “driver”) into the subcritical blanket and reduce neutron loss via axial leakage. The blanket in the S&B core makes beneficial use of the leaking neutrons for improved economics and resource utilization. A specific objective of this study is to maximize the fraction of core power that can be generated by the blanket without violating the thermal hydraulic and material constraints. Since the blanket fuel requires no reprocessing along with remote fuel fabrication, a larger fraction of power from the blanket will result in a smaller fuel recycling capacity and lower fuel cycle cost per unit of electricity generated. A unique synergism is found between a low conversion ratio (CR) seed and a B&B blanket fueled by thorium. Among several benefits, this synergism enables the very low leakage S&B cores to have small positive coolant voiding reactivity coefficient and large enough negative Doppler coefficient even when using inert matrix fuel for the seed. The benefits of this synergism are maximized when using an annular seed surrounded by an inner and outer thorium blankets. Among the high-performance S&B cores designed to benefit from this unique synergism are: (1) the ultra-long cycle core that features a cycle length of ~7 years; (2) the high-transmutation rate core where the seed fuel features a TRU CR of 0.0. Its TRU transmutation rate is comparable to that of the reference Advanced Burner Reactor (ABR) with CR of 0.5 and the thorium blanket can generate close to 60% of the core power; but requires only one sixth of the reprocessing and fabrication capacity per unit of core power. Nevertheless, these high-performance cores were designed to set upper bounds on the S&B core performance by using larger height and pressure drop than those of typical SFR design. A study was subsequently undertaken to quantify the tradeoff between S&B core design variables and the core performance. This study concludes that a viable S&B core can be designed without significant deviation from SFR core design practices. For example, the S&B core with 120cm active height will be comparable in volume, HM mass and specific power with the S-PRISM core and could fit within the S-PRISM reactor vessel. 43% of this core power will be generated by the once-through thorium blanket; the required capacity for reprocessing and remote fuel fabrication per unit of electricity generated will be approximately one fifth of that for a comparable ABR. The sodium void worth of this 120cm tall S&B core is significantly less positive than that of the reference ABR and the Doppler coefficient is only slightly smaller even though the seed uses a fertile-free fuel. The seed in the high transmutation core requires inert matrix fuel (TRU-40Zr) that has been successfully irradiated by the Fuel Cycle Research & Development program. An additional sensitivity analysis was later conducted to remove the bias introduced by the discrepancy between radiation damage constraints -- 200 DPA applied for S&B cores and fast fluence of 4x1023 n(>0.1MeV)/cm2 applied for ABR core design. Although the performance characteristics of the S&B cores are sensitive to the radiation damage constraint applied, the S&B cores offer very significant performance improvements relative to the conventional ABR core design when using identical constraint.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-E-T-16: A Comparison of Expected Dwell Times and Dose Variations for NAMD Patients Treated with An Episcleral Brachytherapy Device

Purpose: To evaluate the variations in dwell times and doses expected when using an episcleral brachytherapy device for treatment of neovascular agerelated macular degeneration (n-AMD) based on accurate imaging modalities Methods: Data from 40 eyes from 40 subjects with known n- AMD acquired through the Distance of Choroid Study (DOCS) conducted at Moorfields Eye Hospital was used to determine the target depth; the distance from the outer scleral surface of the eye, through the choroid, to the apex of the choroidal neovascularization (CNV). Each subject underwent, in triplicate, enhanced-depth Spectral Domain Optical Coherence Tomography (SD-OCT), Swept Source Optical Coherence Tomography, (SS-OCT) and Ocular Ultrasound (O-US). These data are the most comprehensive and accurate measurements of the dimensions of the CNV and adjacent layers of the eye for this cohort of patients. During treatment of n-AMD, patients receive a dose of 24Gy to the apex at the target depth. Using the percentage depth dose for a Sr-90 episcleral brachytherapy device, dwell times and doses to the apex were computed to determine the expected variations. Results: The mean target depth and the 95% confidence interval (CI) determined by combining O-US with SD-OCT were 1326 (956,1696)µm and with SS-OCT were 1332 (970,1693)µm. The calculated corresponding mean dwell times and 95% (CI) were 334 (223,445)s and 335 (226,445)s for SD-OCT and SS-OCT determined depths, respectively. The corresponding mean apex dose and 95% (CI) were 24 (35.9,18.0)Gy (SD-OCT) and 24 (35.6,18.1)Gy (SS-OCT). Conclusion: For episcleral brachytherapy treatment of n-AMD, using a patient population average target depth for treatment planning is inadequate, resulting in dose variations of a factor of approximately two over the 95% CI and larger variations for a nontrivial segment of the population. Each patient should have individualized imaging studies to determine the target depth for use in the dwell time calculation. Study was sponsored by Salutaris Medical Devices, Ltd., a subsidiary of Salutaris Medical Devices, Inc. Hamilton and Marsteller are founders of Salutaris Medical Devices, Inc. Drew, McGovern and Vitali are minor equity holders in Salutaris Medical Devices, Inc.

60 APPLIED LIFE SCIENCES↗

Where's the Beaverhead beef?

Only rare quartz grains with single-set planar (1013) deformation features (PDF's) are present in breccia dikes found in association with uniformly oriented shatter cones that occur over an area 8 x 25 km. This suggests that the Beaverhead shocked rocks come from only the outer part of the central uplift of what must have been a large (greater than 100 km diameter) complex impact structure. An impact event of this magnitude on continental crust (thought to have occurred in late Precambrian or ealy Paleozoic time) could be expected to punctuate local geologic history. Furthermore, although it may now be covered, its scar should remain despite all the considerable subsequent erosion/deposition and tectonism since the impact. The following are three large-scale singularities or anomalies that may reflect the event and mark its source. (1) The Lemhi Arch is a major structural uplift that occurred in late Proterozoic-early Paleozoic time in East Central Idaho and caused the erosion of at least 4 km of sedimentary cover. This may be directly related to the impact. (2) Of the many thrust sheets comprising the Cordilleran belt, the Cabin plate that carries the shocked rocks is unique in that it alone intersected the crystalline basement. It also now marks the apex of the Southwest Montana Recess in the Sevier front. The basement uplift remaining from the impact may have constituted a mechanical obstacle to the advancing thrust sheets in Cretaceous time, causing the recess. (3) What could be interpreted as a roughly circular aeromagnetic anomaly approx. 70 km in diameter can be discerned in the state aeromagnetic map centered about 20 km southeast of Challis, Idaho, in the Lost River range. It is in approximately the right place, and ignoring the possibility that the anomalies have diverse causes and the circular pattern is coincidental, it may mark what remains of the buried central uplift structure.

Hargraves, R. B.↗

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Cassini atmospheric chemistry mapper. Volume 1. Investigation and technical plan

The Cassini Atmospheric Chemistry Mapper (ACM) enables a broad range of atmospheric science investigations for Saturn and Titan by providing high spectral and spatial resolution mapping and occultation capabilities at 3 and 5 microns. ACM can directly address the major atmospheric science objectives for Saturn and for Titan, as defined by the Announcement of Opportunity, with pivotal diagnostic measurements not accessible to any other proposed Cassini instrument. ACM determines mixing ratios for atmospheric molecules from spectral line profiles for an important and extensive volume of the atmosphere of Saturn (and Jupiter). Spatial and vertical profiles of disequilibrium species abundances define Saturn's deep atmosphere, its chemistry, and its vertical transport phenomena. ACM spectral maps provide a unique means to interpret atmospheric conditions in the deep (approximately 1000 bar) atmosphere of Saturn. Deep chemistry and vertical transport is inferred from the vertical and horizontal distribution of a series of disequilibrium species. Solar occultations provide a method to bridge the altitude range in Saturn's (and Titan's) atmosphere that is not accessible to radio science, thermal infrared, and UV spectroscopy with temperature measurements to plus or minus 2K from the analysis of molecular line ratios and to attain an high sensitivity for low-abundance chemical species in the very large column densities that may be achieved during occultations for Saturn. For Titan, ACM solar occultations yield very well resolved (1/6 scale height) vertical mixing ratios column abundances for atmospheric molecular constituents. Occultations also provide for detecting abundant species very high in the upper atmosphere, while at greater depths, detecting the isotopes of C and O, constraining the production mechanisms, and/or sources for the above species. ACM measures the vertical and horizontal distribution of aerosols via their opacity at 3 microns and, particularly, at 5 microns. ACM recovers spatially-resolved atmospheric temperatures in Titan's troposphere via 3- and 5-microns spectral transitions. Together, the mixing ratio profiles and the aerosol distributions are utilized to investigate the photochemistry of the stratosphere and consequent formation processes for aerosols. Finally, ring opacities, observed during solar occultations and in reflected sunlight, provide a measurement of the particle size and distribution of ring material. ACM will be the first high spectral resolution mapping spectrometer on an outer planet mission for atmospheric studies while retaining a high resolution spatial mapping capability. ACM, thus, opens an entirely new range of orbital scientific studies of the origin, physio-chemical evolution and structure of the Saturn and Titan atmospheres. ACM provides high angular resolution spectral maps, viewing nadir and near-limb thermal radiation and reflected sunlight; sounds planetary limbs, spatially resolving vertical profiles to several atmospheric scale heights; and measures solar occultations, mapping both atmospheres and rings. ACM's high spectral and spatial resolution mapping capability is achieved with a simplified Fourier Transform spectrometer with a no-moving parts, physically compact design. ACM's simplicity guarantees an inherent stability essential for reliable performance throughout the lengthy Cassini Orbiter mission.

Smith, William Hayden↗

Multi-Shock Shield Performance at 15 MJ for Catalogued Debris

While orbital debris of ten centimeters or more are tracked and catalogued, the difficulty of finding and accurately accounting for forces acting on the objects near the ten centimeter threshold results in both uncertainty of their presence and location. These challenges result in difficult decisions for operators balancing potential costly operational approaches with system loss risk. In this paper, the assessment of the feasibility of protecting a spacecraft from this catalogued debris is described using numerical simulations and a test of a multi-shock shield system against a cylindrical projectile impacting normal to the surface with approximately 15 MJ of kinetic energy. The hypervelocity impact test has been conducted at the Arnold Engineering Development Complex (AEDC) with a 598 g projectile at 6.905 km/s on a NASA supplied multi-shock shield. The projectile used is a hollow aluminum and nylon cylinder with an outside diameter of 8.6 cm and length of 10.3 cm. Figure 1 illustrates the multi-shock shield test article, which consisted of five separate bumpers, four of which are fiberglass fabric and one of steel mesh, and two rear walls, each consisting of Kevlar fabric. The overall length of the test article was 2.65 m. The test article was a 5X scaled-up version of a smaller multi-shock shield previously tested using a 1.4 cm diameter aluminum projectile for an inflatable module project. The distances represented by S1 and S1/2 in the figure are 61 cm and 30.5 cm, respectively. Prior to the impact test, hydrodynamic simulations indicated that some enhancement to the standard multi-shock system is needed to address the effects of the cylindrical shape of the projectile. Based on the simulations, a steel mesh bumper has been added to the shield configuration to enhance the fragmentation of the projectile. The AEDC test occurred as planned, and the modified NASA multi-shock shield successfully stopped 598 g projectile using 85.6 kg/m(exp 2). The fifth bumper layer remained in tact, although it was torn free from its support structure and thrown into the first rear wall. The outer Kevlar layer of the first rear wall tore likely from the impact of the fifth bumper's support structure, but the back of the rear wall was intact. No damage occurred to the second rear wall, or to the witness plate behind the target.

Miller, J. E.↗

Microstructural Observations on Advanced Multilayer Protective Coatings for Nuclear Thermal Propulsion

The present research discusses the microstructures of a new multilayered coating concept proposed to solve the “midrange corrosion” problem observed in NbC and ZrC-coated graphite (Gr) nuclear fuels during the NERVA/Rover programs. The concept envisions designing a compliant multilayered coatings architecture to accommodate the thermal expansion mismatch strains between the ZrC outer coating and the Gr/(U,Zr)C fuel matrix along with using a diffusion barrier to minimize carbon and hydrogen diffusion, respectively. Multilayered coatings were deposited on Gr substrates by chemical vapor deposition (CVD). Graphite disks, as well as the inner channels of a 19-hole hexagonal rod 152.4 mm long, were coated by CVD in proof-of concept process demonstration studies. Coated disk specimens were thermally-cycled between ambient and 1900 K for as long as 20 cycles with a 1h hold time at the temperature using both stepped and non-interrupted thermal cycling methods. Detailed cross-sectional microstructural analyses were conducted on the thermally-cycled coated disks. These results showed that the coatings were intact after 20 thermal cycles at 1900 K although the bonding of the overlay Mo coating appeared to be influenced by the oxygen content in the layer. Seven transverse sections were cut along the length from the coated hexagonal rod for microstructural studies. All the nineteen channels were observed to have been coated in all the sections. This paper reports the results of detailed microstructural analyses, compositional line scans and x-ray dot maps conducted on three coated channels from a mid-section of the rod. The measured values of carbon (C) were observed to be exceptionally high in all the specimens, which were attributed to hydrocarbon contamination of the specimen surfaces. It was observed that the oxygen (O) levels were high in all the coating layers even after approximately correcting for the observed high values of C. The source of O contamination of the CVD reactor system was attributed to either the presence of oxychloride impurities in the precursor salts or a pin hole leak in the coating system. It is concluded that high values of O result in poor bonding of the Mo overlay to the ZrC layer. These trials have established that it is possible to deposit multilayered coatings by CVD in 19-hole hexagonal rods 152.4 mm long similar to the NERVA fuel rod design although the quality of the coatings and the uniformity of their thicknesses require further improvement.

Nuclear thermal propulsion, multilayer coatings, m↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercial CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and or these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were run to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire left wing without the wheel. well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was nom1 to the spar. It included all o f the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas H.↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercia2 CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and for these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were m to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire lee wing without the wheel well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was normal to the spar. It included all of the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was Found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas↗

Studies using IMP, Voyager and Pioneer cosmic ray data to determine the size of the heliosphere

The purpose of this project was to use the cosmic ray data from the IMP, Voyager and Pioneer spacecraft in the heliosphere out to approximately 65 AU to estimate the size of the heliosphere. We used several techniques to develop a consistent picture of the size of the heliosphere. The first method used a response function approach which determined the intensity as a function of time by scaling the modulation effect as they move outward and eventually reach the boundary of the heliosphere. In this model the effects of transient cosmic ray disturbances is included. A second approach using the perturbation method in which drifts are considered as a perturbation to the standard diffusion-convection modulation models was not fully developed. In a third approach the location of the modulation boundary beyond the termination shock was estimated using observations of the intensity and radial gradients between Voyager 2 and Pioneer 10 along with new estimates of the interstellar intensity of more than 70 MeV galactic cosmic rays. Using this method we found that for 7 years, from 1983 to 1990, the modulation boundary remained constant at 83 +/- 5 AU. We infer from these studies that a modulation boundary can be estimated only by extrapolating the observed radial gradients when the solar magnetic field polarity is such that cosmic-ray particles are drifting in the heliosphere inward toward the Earth along the neutral sheet. The boundary distance is larger than the estimates of the location of the termination shock at 67 +/- 5 AU using the same method. Two other studies partially supported by this grant are attached. The first deals with the recovery period of the greater than 70 MeV cosmic rays in the outer heliosphere from 1992-1995. In the second paper we compare the rigidity dependence of the 11-year cosmic ray variation at the Earth in two cycles of opposite solar magnetic field polarity.

Lockwood, John A.↗

Estimating the Libration of Mercury by Remote Sensing of Gravity and Altimetry

One of the most intriguing scientific questions about Mercury is whether the planet possesses a liquid core. One way that this can be ascertained is through the measurement of the longitudinal libration of Mercury about its rotation axis. The MESSENGER (Mercury Surface, Space, Environment, Geochemistry, and Ranging) mission to Mercury will attempt this measurement by deriving the rotation of the planet from remotely sensed observations of the Mercury gravity field and the planet's shape. For a fully fluid core the libration is expected to be an approximately 350-m-amplitude oscillation (at the equator) about the mean rotation of the planet. Thus, the shape of the planet undergoes a librational oscillation in the same way that features on the surface participate in the motion. In this sense altimetric topography can be considered analogous to imaging, with the important exception that for altimetry, the long wavelengths corresponding to low degrees and orders in a spherical harmonic expansion, can be precisely determined from orbital measurements. The rotation of the planet thus introduces a time dependence of the topographic image in an inertial (i.e. fixed) reference frame and the libration is a variation in the planetary rotation rate, manifest as an oscillation of the topographic image. Similarly, gravity measures the distribution of mass within a planet and as such can be viewed as providing a vertically integrated 'image' of the internal density structure. For a librating planet the gravitational field will also reveal a variation in the rotation rate, manifest as an oscillation of the field with respect to an inertial reference frame. While the long-wavelength gravity field contains contributions from the radial distribution of mass with depth, the long wavelength terms are most sensitive to the mass distribution of the deep interior. Thus, both the topography and the gravity independently can be used to determine the (irregular) rotation of a planet, although the topographical method measures the rotation of the crust while the gravity measures the rotation of the distributed mass, and these are not necessarily identical. For example, for a planet that contains a fluid outer core that is effectively decoupled from the mantle the gravity field will not librate in the manner of a solid body. The libration obtained from the variable rotation of the gravity field will represent a combination of contributions due to a librating mantle and a differentially rotating fluid core. Differential rotation of the Earth's inner core had been theoretically predicted from three dimensional geodynamical models and subsequently reported from seismic observations. A difference between the libration of the lithosphere as determined from altimetry and the deep interior as determined from gravity, if it could be detected, might provide further insight into the nature of core-mantle coupling. Additional information is contained in the original extended abstract.

Smith, D. E.↗

Hydrogen and related materials at high density: Physics, chemistry and planetary implications

Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.

Hemley, R. J.↗