Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 523 records · Page 29

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗

Modeling of the Influence of a Damaged Thermally Grown Oxide (TGO) Layer in an Environmental Barrier Coating System

Environmental barrier coatings (EBC) are an enabling technology for the successful application of ceramic matrix composites (CMCs) in air-breathing gas turbine engines. EBCs are susceptible to a number of failure modes including oxidation/delamination, recessions, chemical attack and dissolution, thermomechanical degradation, erosion, and foreign object impact damage in a combustion environment. Spallation of environmental barrier coating (EBC), induced by a thermally grown oxide (TGO) layer, is a key EBC failure mode. The TGO layer, resulting from steam oxidation, grows either from a silicon bond coat layer (if present) or from a silicon carbide (SiC) based substrate itself. The TGO layer evolves (i.e., thickness increases with time) as water vapor and oxygen gradually diffuse through the EBC, and the EBC spalls off once the TGO layer reaches some critical thickness. The critical thickness of the TGO layer for failure is in the range of 20-30 microns, but it can vary due to exposure temperature, microstructure etc. Current work at NASA Glenn Research Center, under the Revolutionary Tools and Methods (RTM) subproject is aimed at addressing associated failure modes in EBC systems and developing robust analysis tools to aid in the design/analysis of these systems. The objective of the current work is to conduct a sensitivity study to examine the influence of uniformly and non-uniformly grown oxide layers on the associated driving forces leading to spallation of the EBC when subjected to isothermal loading. The effect of damage in the TGO layer on the resulting stress states is also assessed both in uniform and non-uniform TGO layers.

Mital, Subodh K.↗

A model of temporal variations in ozone density in the Martian atmosphere

Temporal variations of the ozone density profile in the Martian atmosphere at high latitudes are calculated for the course of a Martian year, taking into account seasonal and diurnal variations in temperature, water vapor and solar radiation. Calculations are based on a model including 35 neutral photochemical reactions, and vertical eddy diffusion using a time step of 12 min for the region from the surface to 240 km altitude. Results are found to be in better agreement with Mariner 9 observations of the time and magnitude of the seasonal maximum than previous model calculations. The diurnal variation is predicted to be small near the solstices, with the nighttime ozone density greater than the daytime and the magnitude of the difference dependent on season. Opposite temporal variations are predicted for ozone densities above and below about 25 km, and an ozone density maximum at 35-40 km is obtained. It is suggested that the effects of an aerosol layer may not be important in enhancing predicted ozone concentration, and may even decrease it.

Shimazaki, T.↗

Global Mapping of Underwater UV Irradiances and DNA-Weighted Exposures using TOMS and SeaWiFS Data Products

The global stratospheric ozone-layer depletion results In an increase in biologically harmful ultraviolet (UV) radiation reaching the surface and penetrating to ecologically significant depths in natural waters. Such an increase can be estimated on a global scale by combining satellite estimates of UV irradiance at the ocean surface from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWIFS satellite ocean-color measurements in the visible spectral region. In this paper we propose a model of seawater optical properties in the UV spectral region based on the Case I water model in the visible range. The inputs to the model are standard monthly SeaWiFS products: chlorophyll concentration and the diffuse attenuation coefficient at 490nm. Penetration of solar UV radiation to different depths in open ocean waters is calculated using the RT (radiative transfer) quasi-single scattering approximation (QSSA). The accuracy of the QSSA approximation in the water is tested using more accurate codes. The sensitivity study of the underwater UV irradiance to atmospheric and oceanic optical properties have shown that the main environmental parameters controlling the absolute levels of the UVB (280-320nm) and DNA-weighted irradiance underwater are: solar-zenith angle, cloud transmittance, water optical properties, and total ozone. Weekly maps of underwater UV irradiance and DNA-weighted exposure are calculated using monthly-mean SeaWiFS chlorophyll and diffuse attenuation coefficient products, daily SeaWiFS cloud fraction data, and the TOMS-derived surface UV irradiance daily maps. The final products include global maps of weekly-average UVB irradiance and DNA-weighted daily exposures at 3m and 10m, and depths where the UVB irradiance and DNA-weighted dose rate at local noon are equal to 10% of their surface values.

Vasilkov, Alexander↗

Rapid Numerical Simulation of Viscous Axisymmetric Flow Fields

A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmetric flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thin-layer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results are in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.

Tweedt, Daniel L.↗

Structure of Microgravity Transitional and Pulsed Jet Diffusion Flames

This paper describes results obtained in a study of pulsed gas jet diffusion flames to better characterize the recently observed vortex/flame interactions in microgravity transitional and turbulent diffusion flames, and to improve the understanding of large-scale structures in corresponding normal-gravity flames. In preparation for a space experiment, tests were conducted in the 5.18-Second Zero-Gravity Facility of the NASA Lewis Research Center. Both unpulsed and pulsed laminar flames were studied and numerical modeling of these flames was carried out for data comparison and model validation. In addition, complementary tests for a series of unpulsed flames were conducted on-board the NASA KC-135 research aircraft. The microgravity transitional and turbulent gas-jet diffusion flames have been observed to be dominated by large-scale disturbances, or structures. These structures first appear intermittently in the flame at Reynolds numbers (based on the cold jet injection properties) of about 2100. With increase in injection Reynolds number, the rate of intermittent disturbances increases until the generation becomes continuous at Reynolds numbers of 3000 and higher. The behavior of these structures depends upon the velocity and temperature characteristics of the jet/flame shear layer. These characteristics are different in normal gravity and microgravity.

Bahadori, M. Yousef↗

Damped Windows for Aircraft Interior Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission into aircraft. To improve the acoustical performance, damped windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. In this paper, numerical and experimental results are used to evaluate the acoustic benefits of damped windows. Tests were performed in the Structural Acoustic Loads and Transmission Facility at NASA Langley Research Center to measure the transmission loss for diffuse acoustic excitation and radiated sound power for point force excitation. Comparisons between uniform and damped plexiglas windows showed increased transmission loss of 6 dB at the first natural frequency, 6 dB at coincidence, and 4.5 dB over a 50 to 4k Hz range. Radiated sound power was reduced up to 7 dB at the lower natural frequencies and 3.7 dB over a 1000 Hz bandwidth. Numerical models are presented for the prediction of radiated sound power for point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. A parametric study is presented that evaluates the optimum configuration of the damped plexiglas windows for reducing the radiated sound power.

Buehrle, Ralph D.↗

Transport properties for a mixture of the ablation products C, C2, and C3

The ablation of carbon-phenolic heat shields upon entry into the atmosphere of one of the outer planets leads to the injection of large amounts of C, C2, and C3 into the shock layer. These species must be included in the calculation of transport properties in the shock layer. The kinetic theory of gases has been used to obtain accurate results for the transport properties of monatomic carbon. The Hulburt-Hirschelder potential, the most accurate general purpose atom-atom potential for states with an attractive minimum, was used to represent such states and repulsive states were represented by fitting quantum mechanical potential energy curves with the exponential repulsive potential. These results were orientation averaged according to the peripheral force model to obtain transport collision integrals for the C-C2 and C2-C2 interaction. Results for C3 were obtained by ignoring the presence of the central carbon atom. The thermal conductivity, viscosity, and diffusion coefficients for pure C, C2, and C3, and for mixtures of these gases, were then calculated from 1000 K - 25,000 K.

Biolsi, L.↗

Characteristics of laminar gas jet diffusion flames under the influence of elevated gravity

Laminar gas jet diffusion flames of hydrogen, methane, ethane, and propane, stabilized on small cylindrical burners, were studied at elevated gravity in a centrifuge to investigate the role of buoyancy in such flames. Photographic observation of the hydrocarbon flames indicated that length and carbon luminosity decreased with increasing buoyancy. These flames eventually separated from the burner rim and finally extinguished as gravity increased. Hydrogen flames, however, did not separate or extinguish up to about 31 times earth normal gravity. Length data are reported and interpreted with the aid of a simplified model which describes the flame as a heated vertical cylindrical sheet to which oxygen is supplied by a free convection boundary layer.

Altenkirch, R. A.↗

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

Wind-US Results for the AIAA 1st Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 1st Propulsion Aerodynamics Workshop. The The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from representatives from government, industry, academia, and commercial software companies. Participants were were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a series of convergent convergent nozzles for a range of nozzle pressure ratios between 1.4 and 7.0. These configurations included a included a reference axisymmetric nozzle as well as 15deg , 25deg , and 40deg conical nozzles. Participants were also asked also asked to examine the plume shock structure for two cases where the 25deg conical nozzle was bifurcated by a bifurcated by a solid plate. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total and total pressure at downstream rake locations were examined.

Yoder, Dennis↗

Photospheric Magnetic Flux Transport - Supergranules Rule

Observations of the transport of magnetic flux in the Sun's photosphere show that active region magnetic flux is carried far from its origin by a combination of flows. These flows have previously been identified and modeled as separate axisymmetric processes: differential rotation, meridional flow, and supergranule diffusion. Experiments with a surface convective flow model reveal that the true nature of this transport is advection by the non-axisymmetric cellular flows themselves - supergranules. Magnetic elements are transported to the boundaries of the cells and then follow the evolving boundaries. The convective flows in supergranules have peak velocities near 500 m/s. These flows completely overpower the superimposed 20 m/s meridional flow and 100 m/s differential rotation. The magnetic elements remain pinned at the supergranule boundaries. Experiments with and without the superimposed axisymmetric photospheric flows show that the axisymmetric transport of magnetic flux is controlled by the advection of the cellular pattern by underlying flows representative of deeper layers. The magnetic elements follow the differential rotation and meridional flow associated with the convection cells themselves -- supergranules rule!

Hathaway, David H.↗

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbo-machinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved by a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage scheme. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

Estimation of minority carrier diffusion lengths in InP/GaAs solar cells

Minority carrier diffusion length is one of the most important parameters affecting the solar cell performance. An attempt is made to estimate the minority carrier diffusion lengths is the emitter and base of InP/GaAs heteroepitaxial solar cells. The PC-1D computer model was used to simulate the experimental cell results measured at NASA Lewis under AMO (air mass zero) spectrum at 25 C. A 16 nm hole diffusion length in the emitter and a 0.42 micron electron diffusion length in the base gave very good agreement with the I-V curve. The effect of varying minority carrier diffusion lengths on cell short current, open circuit voltage, and efficiency was studied. It is also observed that the front surface recombination velocity has very little influence on the cell performance. The poor output of heteroepitaxial cells is caused primarily by the large number of dislocations generated at the interfaces that propagate through the bulk indium phosphide layers. Cell efficiency as a function of dislocation density was calculated and the effect of improved emitter bulk properties on cell efficiency is presented. It is found that cells with over 16 percent efficiencies should be possible, provided the dislocation density is below 10(exp 6)/sq cm.

Jain, R. K.↗

Dark Halo Craters and the Thickness of Grooved Terrain on Ganymede

Dark halo craters on grooved terrain on Ganymede represent potential probes of the subsurface geology. These craters are surrounded by a broad, diffuse low albedo annulus, or halo and have apparently lost the bright rim and ray deposits associated with very young craters. The recognition of dark halo craters, almost all greater than 12 km in diameter, indicates that material darker and redder than grooved terrain material forms a stratigraphy horizon at about 1 km depth the present surface of grooved terrain in Uruk Sulcus. This material is most likely downdropped cratered terrain material, which only larger craters have excavated into. This is most consistent with extensional tectonic models for grooved terrain formation whereby the lithosphere is stretched and blocks of ancient cratered terrain are downdropped along bounding faults, and subsequently resurfaced by a shallow layer of relatively clean icy material. Evidence for the preservation of cratered terrain material at shallow depth beneath grooved terrain material poses great difficulties for alternative grooved and smoothed terrain formation mechanisms.

Schenk, P. M.↗

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗