Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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 883 records · Page 49

A Linearized Prognostic Cloud Scheme in NASAs Goddard Earth Observing System Data Assimilation Tools

A linearized prognostic cloud scheme has been developed to accompany the linearized convection scheme recently implemented in NASA's Goddard Earth Observing System data assimilation tools. The linearization, developed from the nonlinear cloud scheme, treats cloud variables prognostically so they are subject to linearized advection, diffusion, generation, and evaporation. Four linearized cloud variables are modeled, the ice and water phases of clouds generated by large-scale condensation and, separately, by detraining convection. For each species the scheme models their sources, sublimation, evaporation, and autoconversion. Large-scale, anvil and convective species of precipitation are modeled and evaporated. The cloud scheme exhibits linearity and realistic perturbation growth, except around the generation of clouds through large-scale condensation. Discontinuities and steep gradients are widely used here and severe problems occur in the calculation of cloud fraction. For data assimilation applications this poor behavior is controlled by replacing this part of the scheme with a perturbation model. For observation impacts, where efficiency is less of a concern, a filtering is developed that examines the Jacobian. The replacement scheme is only invoked if Jacobian elements or eigenvalues violate a series of tuned constants. The linearized prognostic cloud scheme is tested by comparing the linear and nonlinear perturbation trajectories for 6-, 12-, and 24-h forecast times. The tangent linear model performs well and perturbations of clouds are well captured for the lead times of interest.

Precipitation↗

Burn rates of TiH2/KClO4/Viton and output testing of NASA SKD26100098-301 pressure cartridges

The burn rates of the pyrotechnic TiH2/KClO4/Viton with a mass ratio of 30/65/5 have been measured as a function of pressure in nitrogen up to 312 MPa(45 Kpsi). The burn rates were fit to R = a pn, with a = 2.055 cm/sec/MPan and n = 0.472 between 0.15 MPa (22 psi) and 21.6 MPa (3.13 Kpsi) and a = 4.38 cm/sec/MPan and n = 0.266 between 70 MPa (10.15 Kpsi) and 312 MPa (45.25 Kpsi). The decrease in slope at the higher pressures is attributed to a diffusion limited reaction. No acoustically driven flame instabilities or large conductive-to-convective burn transitions were observed. Solid reaction products were analyzed by x-ray diffraction and scanning electron microscopy (SEM). X-ray diffraction detected only TiO2 and KC1. SEM showed that the particle size of the reaction products increased as the nitrogen pressure increased. There were no anomalous characteristics of the burn of this pyrotechnic that could be interpreted as a cause of the o-ring blow-by problem in the forward shear bolt assembly. Three NASA SKD26100098-301 pressure cartridges were fired into a fixed volume vessel that was sealed with an O-ring. A maximum pressure of 181.7 MPa(26,350 psi) was reached in around 100 ,mu sec for two shots fired into a volume of 16.3 cm3(0.996 in3). A maximum pressure of 33,460 psi was reached for one shot fired into a volume of 9.55 cm3(0.583 in3). The O-ring burned through on one shot in the larger volume and leaked on the other two thereby simulating the effects of an O-ring leak. The results imply that the piston in the shear bolt assembly would receive a large impulse even if there was a leak in an O-ring seal.

Holy, John A.↗

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

Experimental/Analytical Characterization of the RBCC Rocket-Ejector Mode

The experimental/analytical research work described here addresses the rocket-ejector mode (Mach 0-2 operational range) of the RBCC engine. The experimental phase of the program includes studying the mixing and combustion characteristics of the rocket-ejector system utilizing state-of-the-art diagnostic techniques. A two-dimensional variable geometry rocket-ejector system with enhanced optical access was utilized as the experimental platform. The goals of the experimental phase of the research being conducted at Penn State are to: (a) systematically increase the range of rocket-ejector understanding over a wide range of flow/geometry parameters and (b) provide a comprehensive data base for evaluating and anchoring CFD codes. Concurrent with the experimental activities, a CFD code benchmarking effort at Marshall Space Flight Center is also being used to further investigate the RBCC rocket-ejector mode. Experiments involving the single rocket based optically-accessible rocket-ejector system have been conducted for Diffusion and Afterburning (DAB) as well as Simultaneous Mixing and Combustion configurations. For the DAB configuration, air is introduced (direct-connect) or ejected (sea-level static) into a constant area mixer section with a centrally located gaseous oxygen (GO2)/gaseous hydrogen (GH2) rocket combustor. The downstream flowpath for this configuration includes a diffuser, an afterburner and a final converging nozzle. For the SMC configuration, the rocket is centrally located in a slightly divergent duct. For all tested configurations, global measurements of the axial pressure and heat transfer profiles as well as the overall engine thrust were made. Detailed measurements include major species concentration (H2 O2 N2 and H2O) profiles at various mixer locations made using Raman spectroscopy. Complementary CFD calculations of the flowfield at the experimental conditions also provide additional information on the physics of the problem. These calculations are being conducted at Marshall Space Flight Center to benchmark the FDNS code for RBCC engine operations for such configurations. The primary fluid physics of interests are the mixing and interaction of the rocket plume and secondary flow, subsequent combustion of the fuel rich rocket exhaust with the secondary flow and combustion of the injected afterburner flow. The CFD results are compared to static pressure along the RBCC duct walls, Raman Spectroscopy specie distribution data at several axial locations, net engine thrust and entrained air for the SLS cases. The CFD results compare reasonably well with the experimental results.

Ruf, J. H.↗

COSIM: A Finite-Difference Computer Model to Predict Ternary Concentration Profiles Associated With Oxidation and Interdiffusion of Overlay-Coated Substrates

A finite-difference computer program (COSIM) has been written which models the one-dimensional, diffusional transport associated with high-temperature oxidation and interdiffusion of overlay-coated substrates. The program predicts concentration profiles for up to three elements in the coating and substrate after various oxidation exposures. Surface recession due to solute loss is also predicted. Ternary cross terms and concentration-dependent diffusion coefficients are taken into account. The program also incorporates a previously-developed oxide growth and spalling model to simulate either isothermal or cyclic oxidation exposures. In addition to predicting concentration profiles after various oxidation exposures, the program can also be used to predict coating life based on a concentration dependent failure criterion (e.g., surface solute content drops to 2%). The computer code is written in FORTRAN and employs numerous subroutines to make the program flexible and easily modifiable to other coating oxidation problems.

Nesbitt, James A.↗

Superconducting Mixers for Far-Infrared Spectroscopy

The goal of this project was to fabricate and test planar arrays of superconducting mixers for the 2-6 THz band. The technology is intended for multi-beam receivers aboard Explorer-class missions and the SOFIA Airborne Observatory. The mixer technology is the superconducting transition-edge microbolometer, which is more commonly known as the Hot-Electron micro-Bolometer (HEB). As originally proposed, two superconducting technologies were to be developed: (1) low-Tc niobium HEBs which could approach quantum-noise-limited sensitivities but require cooling to 2- 4 K, and (2) high-Tc YBCO HEBs with sensitivities 10 times worse but with a relaxed cooling requirement of 30-60 K. The low-Tc devices would be best for astronomy applications on SOFIA, whereas the high-Tc devices would be more suitable for planetary missions using systems without stored cryogens. The work plan called for planar micro-fabrication and initial testing of HEB devices at the NIST Boulder clean-room facility. Subsequent assembly and RF testing of selected devices would be done at the CASA laboratory at U. Colorado. Approximately 1-year after work began on this project, Dr. Eyal Gerecht joined the NIST group, and assumed day-to-day responsibility for Nb-HEB development at NIST outside of micro-fabrication. The YBCO-HEB work was to be guided by Dr. Ron Ono, who was the NIST expert in YBCO technology. Unfortunately, recurrent health problems limited the time Ron could devote to the project in its first year. These problems became aggravated in early 2001, and sadly led to Ron's death in October, 2001. His loss was not only a blow to his friends and associates at NIST, but was mounted by the US superconductivity community at large. With his passing, work on high-Tc HEBs ceased at NIST. There was no one to replace him or his expertise. Our work subsequently shifted solely to Nb-HEB devices. In the sections which follow, our progress in the development of diffusion-cooled Nb-HEB mixers is detailed. To simplify the terminology, these devices will subsequently be called DHEB mixers to distinguish them from phonon-cooled devices (PHEBs).

Betz, A. L.↗

Atmospheric Correction for Satellite Ocean Color Radiometry

This tutorial is an introduction to atmospheric correction in general and also documentation of the atmospheric correction algorithms currently implemented by the NASA Ocean Biology Processing Group (OBPG) for processing ocean color data from satellite-borne sensors such as MODIS and VIIRS. The intended audience is graduate students or others who are encountering this topic for the first time. The tutorial is in two parts. Part I discusses the generic atmospheric correction problem. The magnitude and nature of the problem are first illustrated with numerical results generated by a coupled ocean-atmosphere radiative transfer model. That code allow the various contributions (Rayleigh and aerosol path radiance, surface reflectance, water-leaving radiance, etc.) to the topof- the-atmosphere (TOA) radiance to be separated out. Particular attention is then paid to the definition, calculation, and interpretation of the so-called "exact normalized water-leaving radiance" and its equivalent reflectance. Part I ends with chapters on the calculation of direct and diffuse atmospheric transmittances, and on how vicarious calibration is performed. Part II then describes one by one the particular algorithms currently used by the OBPG to effect the various steps of the atmospheric correction process, viz. the corrections for absorption and scattering by gases and aerosols, Sun and sky reflectance by the sea surface and whitecaps, and finally corrections for sensor out-of-band response and polarization effects. One goal of the tutorial-guided by teaching needs- is to distill the results of dozens of papers published over several decades of research in atmospheric correction for ocean color remote sensing.

MODIS↗

Observational Signatures of Coronal Heating in Magnetohydrodynamic Simulations without Radiation or a Lower Atmosphere

It is extremely difficult to simulate the details of coronal heating and also make meaningful predictions of the emitted radiation. Thus, testing realistic models with observations is a major challenge. Observational signatures of coronal heating depend crucially on radiation, thermal conduction, and the exchange of mass and energy with the transition region and chromosphere below. Many magnetohydrodynamic simulation studies do not include these effects, opting instead to devote computational resources to the magnetic aspects of the problem. We have developed a simple method of accounting approximately for the missing effects. It is applied to the simulation output ex post facto and therefore may be a valuable tool for many studies. We have used it to predict the emission from a model corona that is driven by vortical boundary motions meant to represent photospheric convection. We find that individual magnetic strands experience short-term brightenings, both scattered throughout the computational volume and in localized clusters. The former may explain the diffuse component of the observed corona, while the latter may explain bright coronal loops. Several observed properties of loops are reproduced reasonably well: width, lifetime, and quasi-circular cross section (aspect ratio not high). Our results lend support to the idea that loops are multistranded structures heated by "storms" of nanoflares.

James A. Klimchuk↗

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

Aerodynamic and gasdynamic effects in cosmogony

An improved numerical code was constructed to model coagulation and settling of particles in disk nebula containing generic turbulence with arbitrary velocities in the gas. The turbulence is assumed to have a Kolmogorov eddy spectrum. Relative velocities of particles, which lead to collisions and possible coagulation, are computed as due all significant causes in their appropriate regimes: thermal motion, shear and inertial effects in turbulent eddies, and systematic motions due to settling and non-keplerian rotation of the gas. Significant improvements to this program were produced. One significant problem was the disparity of timescales for turbulent mixing and coagulation. To accurately compute the former, the timestep must be shorter than the smallest spatial scale (layer thickness) divided by the turbulent velocity. However, the size distribution often varies due to coagulation on much longer timescales. To minimize the computational overhead associated with collisions between particles of all sizes, a dual timestep was introduced. Collisional changes in the size distribution are computed once in every N substeps, where the substep is controlled by the turbulent diffusion velocity, and N is determined by the rate of collisions. This algorithm allowed simulations to be extended to longer times and later stages.

Weidenschilling, S. J.↗

NASA Solar Array Demonstrates Commercial Potential

A state-of-the-art solar-panel array demonstration site at NASA's Dryden Flight Research Center provides a unique opportunity for studying the latest in high-efficiency solar photovoltaic cells. This five-kilowatt solar-array site (see Figure 1) is a technology-transfer and commercialization success for NASA. Among the solar cells at this site are cells of a type that was developed in Dryden Flight Research Center s Environmental Research Aircraft and Sensor Technology (ERAST) program for use in NASA s Helios solar-powered airplane. This cell type, now denoted as A-300, has since been transferred to SunPower Corporation of Sunnyvale, California, enabling mass production of the cells for the commercial market. High efficiency separates these advanced cells from typical previously commercially available solar cells: Whereas typical previously commercially available cells are 12 to 15 percent efficient at converting sunlight to electricity, these advanced cells exhibit efficiencies approaching 23 percent. The increase in efficiency is due largely to the routing of electrical connections behind the cells (see Figure 2). This approach to increasing efficiency originated as a solution to the problem of maximizing the degree of utilization of the limited space available atop the wing of the Helios airplane. In retrospect, the solar cells in use at this site could be used on Helios, but the best cells otherwise commercially available could not be so used, because of their lower efficiencies. Historically, solar cells have been fabricated by use of methods that are common in the semiconductor industry. One of these methods includes the use of photolithography to define the rear electrical-contact features - diffusions, contact openings, and fingers. SunPower uses these methods to produce the advanced cells. To reduce fabrication costs, SunPower continues to explore new methods to define the rear electrical-contact features. The equipment at the demonstration site includes two fixed-angle solar arrays and one single-axis Sun-tracking array. One of the fixed arrays contains typical less-efficient commercial solar cells and is being used as a baseline for comparison of the other fixed array, which contains the advanced cells. The Sun-tracking array tilts to follow the Sun, using an advanced, real-time tracking device rather than customary pre-programmed mechanisms. Part of the purpose served by the demonstration is to enable determination of any potential advantage of a tracking array over a fixed array. The arrays are monitored remotely on a computer that displays pertinent information regarding the functioning of the arrays.

Creech, Gray↗

Astronomical Problems and Laboratory Solutions

Studying the chemical and isotopic composition of interstellar ice and dust provides insight into the composition and chemical history of the solid bodies in the solar nebula and the nature of the material subsequently brought into the inner part of the solar system by comets and meteorites. It is now possible to probe the composition of these microscopic interstellar particles (some hundreds of light years away), thanks to substantial progress in two areas: astronomical spectroscopic techniques in the middle-infrared, the spectral region most diagnostic of composition; and laboratory simulations which realistically reproduce the critical conditions in various interstellar environments. High quality infrared spectra of many different astronomical sources, some associated with dark molecular clouds, and others in the diffuse interstellar medium (DISM) are now available. What comparisons of these spectra with laboratory spectra tell us about the complex organic components of these materials is the subject of this talk. Most interstellar material is concentrated in large molecular clouds where simple molecules are formed by gas phase and dust grain surface reactions. Gaseous species (except H2) striking the cold (10K) dust will stick, forming an icy grain mantle. This accretion, coupled with energetic particle bombardment and UV photolysis, will produce a complex chemical mixture containing volatile, non-volatile, and isotopically fractionated species. One can compare spectra of the diffuse and dense interstellar medium with the spectra of analogs produced in the laboratory under conditions which mimic those in these different environments. In this way one can determine the composition and abundances of the major constituents present and place general constraints on the types and relative abundances of organics coating the grains. Ices in dense clouds contain H2O, CH3OH, CO, perhaps some NH3 and H2CO, as well as nitriles and ketones or esters. there is some evidence that the later, more complex species, are also present on the grains in the DISM. The evidence for these materials, in addition to carbon rich materials such as amorphous carbon, microdiamonds, and polycyclic aromatic hydrocarbons will be reviewed and the possible connection with meteoritic organics will be discussed.

Allamandola, L. J.↗

Optical design of the Diffuse Infrared Background Experiment for NASA's Cosmic Background Explorer

The conceptual design for a ten-band absolute filter photometer (the Diffuse Infrared Background Experiment) to operate at 2 K and measure galactic and extragalactic infrared radiation in the 1 to 300-micron range and polarization in the 1 to 3.5-micron range is presented, as part of the NASA Cosmic Background Explorer. The telescope optical design, a Gregorian design incorporating bafffles and shades to provide high stray-light rejection, is described. Pupil nonuniformity in the detector-assembly optical design has been limited. It is determined that detector sensitiity requirements can be met, and that the problem of radiation-induced responsivity variations can be solved by minimizing detector-assembly size, providing for in situ thermal annealing, and allowing for frequent detector calibration. Limitations on mirror performance are to be met by fabricating mirrors and structure from the same aluminum 6061 ingot.

Miller, M. S.↗

A time-accurate implicit method for chemically reacting flows at all Mach numbers

The objective of this work is to develop a unified solution algorithm capable of treating time-accurate chemically reacting flows at all Mach numbers, ranging from molecular diffusion velocities to supersonic speeds. A rescaled pressure term is used in the momentum equation to circumvent the singular behavior of pressure at low Mach numbers. A dual time-stepping integration procedure is established. The system eigenvalues become well behaved and have the same order of magnitude, even in the very low Mach number regime. The computational efficiency for moderate and high speed flow is competitive with the conventional density-based scheme. The capabilities of the algorithm are demonstrated by applying it to selected model problems including nozzle flows and flame dynamics.

Withington, J. P.↗

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗

COSIM: A Finite-Difference Computer Model to Predict Ternary Concentration Profiles Associated with Oxidation and Interdiffusion of Overlay-Coated Substrates

A finite-difference computer program (COSIM) has been written which models the one-dimensional, diffusional transport associated with high-temperature oxidation and interdiffusion of overlay-coated substrates. The program predicts concentration profiles for up to three elements in the coating and substrate after various oxidation exposures. Surface recession due to solute loss is also predicted. Ternary cross terms and concentration-dependent diffusion coefficients are taken into account. The program also incorporates a previously-developed oxide growth and spalling model to simulate either isothermal or cyclic oxidation exposures. In addition to predicting concentration profiles after various oxidation exposures, the program can also be used to predict coating fife based on a concentration dependent failure criterion (e.g., surface solute content drops to two percent). The computer code, written in an extension of FORTRAN 77, employs numerous subroutines to make the program flexible and easily modifiable to other coating oxidation problems.

Nesbitt, James A.↗

X-Ray Dust Scattering At Small Angles: The Complete Halo Around GX13+1

The exquisite angular resolution available with Chandra should allow precision measurements of faint diffuse emission surrounding bright sources, such as the X-ray scattering halos created by interstellar dust. However, the ACIS CCDs suffer from pileup when observing bright sources, and this creates difficulties when trying to extract the scattered halo near the source. The initial study of the X-ray halo around GX13+1 using only the ACIS-I detector done by Smith, Edgar & Shafer (2002) suffered from a lack of sensitivity within 50" of the source, limiting what conclusions could be drawn. To address this problem, observations of GX13+1 were obtained with the Chandra HRC-I and simultaneously with the RXTE PCA. Combined with the existing ACIS-I data, this allowed measurements of the X-ray halo between 2-1000". After considering a range of dust models, each assumed to be smoothly distributed with or without a dense cloud along the line of sight, the results show that there is no evidence in this data for a dense cloud near the source, as suggested by Xiang et al. (2005). In addition, although no model leads to formally acceptable results, the Weingartner & Draine (2001) and all but one of the composite grain models from Zubko, Dwek & Arendt (2004) give particularly poor fits.

Smith, Randall K.↗

Integrated design of structures, controls, and materials

In this talk we shall discuss algorithms and CAD tools for the design and analysis of structures for high performance applications using advanced composite materials. An extensive mathematical theory for optimal structural (e.g., shape) design was developed over the past thirty years. Aspects of this theory have been used in the design of components for hypersonic vehicles and thermal diffusion systems based on homogeneous materials. Enhancement of the design methods to include optimization of the microstructure of the component is a significant innovation which can lead to major enhancements in component performance. Our work is focused on the adaptation of existing theories of optimal structural design (e.g., optimal shape design) to treat the design of structures using advanced composite materials (e.g., fiber reinforced, resin matrix materials). In this talk we shall discuss models and algorithms for the design of simple structures from composite materials, focussing on a problem in thermal management. We shall also discuss methods for the integration of active structural controls into the design process.

Blankenship, G. L.↗