Diffusion wave phenomena
Thermal and bleaching waves as diffusion phenomena in opposite limits, discussing laser heating
SEARCH · Engineering Papers
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.
Thermal and bleaching waves as diffusion phenomena in opposite limits, discussing laser heating
This work is based on the study of the photochemical balance of molecular oxygen in the martian atmosphere by using a one-dimensional model of photochemical reactions involving species derived from CO2 and H2O. The model is basically similar to one used previously for the study of the regulation of CO on a global scale, but the chemical rates are taken from another source. In the present scheme, the regulation of molecular oxygen is studied over timescales of the order of its photochemical lifetime (approximately equals 30 yr), which is much shorter than typical escape timescales. Thus, the escape fluxes are fixed to the values given by 3 and 4. We examine the calculated equilibrium abundances of O2 for given thermal, eddy diffusion coefficients and H2O profiles. The thermal profile is taken from in the lower atmosphere. At higher levels, in order to include the diurnal and seasonal thermal profile variability, we have also used the IRTM data. In order to study the influence of both temperature and pressure profiles on the O2 mixing ratio, we have made several tests corresponding to different martian seasons. The results show that the influence of pressure and temperature is quantitatively weak compared to the one of K and of the water vapor density (H2O). Thus, in the following we have fixed the pressure at the surface to a value of 7 mbar and we have used unique standard thermal profile corresponding to a profile roughly averaged over the year, the season, and the day: T equal 205 K at 0 km altitude, 175 K at 25 km, and 145 K at 50 km.
An optical model of solar sail material originally derived at JPL in 1978 has since served as the de facto standard for NASA and other solar sail researchers. The optical model includes terms for specular and diffuse reflection, thermal emission, and non-Lambertian diffuse reflection. The standard coefficients for these terms are based on tests of 2.5 micrometer Kapton sail material coated with 100 nm of aluminum on the front side and chromium on the back side. The original derivation of these coefficients was documented in an internal JPL technical memorandum that is no longer available. Additionally more recent optical testing has taken place and different materials have been used or are under consideration by various researchers for solar sails. Here, where possible, we re-derive the optical coefficients from the 1978 model and update them to accommodate newer test results and sail material. The source of the commonly used value for the front side non-Lambertian coefficient is not clear, so we investigate that coefficient in detail. Although this research is primarily designed to support the upcoming NASA NEA Scout and Lunar Flashlight solar sail missions, the results are also of interest to the wider solar sail community.
Heat transfer, and thermal control characteristics of diffuse bladed louver system for spacecraft
We report on the diffusive-ballistic thermal conductance of multi-moded single-crystal silicon beams measured below 1 K. It is shown that the phonon mean-free-path is a strong function of the surface roughness characteristics of the beams. This effect is enhanced in diffuse beams with lengths much larger than, even when the surface is fairly smooth, 510 nm rms, and the peak thermal wavelength is 0.6 microns. Resonant phonon scattering has been observed in beams with a pitted surface morphology and characteristic pit depth of 30 nm. Hence, if the surface roughness is not adequately controlled, the thermal conductance can vary significantly for diffuse beams fabricated across a wafer. In contrast, when the beam length is of order, the conductance is dominated by ballistic transport and is effectively set by the beam cross-sectional area. We have demonstrated a uniformity of +/-8% in fractional deviation for ballistic beams, and this deviation is largely set by the thermal conductance of diffuse beams that support the micro-electro-mechanical device and electrical leads. In addition, we have found no evidence for excess specific heat in single-crystal silicon membranes. This allows for the precise control of the device heat capacity with normal metal films. We discuss the results in the context of the design and fabrication of large-format arrays of far-infrared and millimeter wavelength cryogenic detectors.
A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.
The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.
A totally automated facility for semiconductor oxidation and diffusion was developed using a state-of-the-art diffusion furnace and high temperature grown oxides. Major innovations include: (1) a process controller specifically for semiconductor processing; (2) an automatic loading system to accept wafers from an air track, insert them into a quartz carrier and then place the carrier on a paddle for insertion into the furnace; (3) automatic unloading of the wafers back onto the air track, and (4) boron diffusion using diborane with plus or minus 5 percent uniformity. Processes demonstrated include Wet and dry oxidation for general use and for gate oxide, boron diffusion, phosphorous diffusion, and sintering.
The energy balance in the lower transition region is analyzed by constructing theoretical models which satisfy the energy balance constraint. The energy balance is achieved by balancing the radiative losses and the energy flowing downward from the corona. This energy flow is mainly in two forms: conductive heat flow and hydrogen ionization energy flow due to ambipolar diffusion. Hydrostatic equilibrium is assumed, and, in a first calculation, local mechanical heating and Joule heating are ignored. In a second model, some mechanical heating compatible with chromospheric energy-balance calculations is introduced. The models are computed for a partial non-LTE approach in which radiation departs strongly from LTE but particles depart from Maxwellian distributions only to first order. The results, which apply to cases where the magnetic field is either absent, or uniform and vertical, are compared with the observed Lyman lines and continuum from the average quiet sun. The approximate agreement suggests that this type of model can roughly explain the observed intensities in a physically meaningful way, assuming only a few free parameters specified as chromospheric boundary conditions.
A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.
This paper investigates the role played by vapor-phase chemical reactions on CVD rates by comparing the results of two extreme theories developed to predict CVD mass transport rates in the absence of interfacial kinetic barrier: one based on chemically frozen boundary layer and the other based on local thermochemical equilibrium. Both theories consider laminar convective-diffusion boundary layers at high Reynolds numbers and include thermal (Soret) diffusion and variable property effects. As an example, Na2SO4 deposition was studied. It was found that gas phase reactions have no important role on Na2SO4 deposition rates and on the predictions of the theories. The implications of the predictions of the two theories to other CVD systems are discussed.
The seasonal behavior of the global ionosphere was studied using a time-dependent three-dimensional physical model (developed by Shunk and his coworkers) of the ionosphere at altitudes between 120 and 800 km. This model accounts for field-aligned diffusion, cross-field electrodynamic drifts both the equatorial region and at high latitudes, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, and local heating and cooling processes. The model studies were carried out for both June and December solstice conditions at solar maximum and for low geomagnetic activity. The ionospheric features predicted by the model agreed qualitatively with the available measurements.
A physical model for the effective thermal conductivity of water frost is proposed for application to the full range of frost density. The proposed model builds on the Zehner-Schlunder one-dimensional formulation for porous media appropriate for solid-to-fluid thermal conductivity ratios less than about 1000. By superposing the effects of mass diffusion and eddy convection on stagnant conduction in the fluid, the total effective thermal conductivity of frost is shown to be satisfactorily described. It is shown that the effects of vapor diffusion and eddy convection on the frost conductivity are of the same order. The results also point out that idealization of the frost structure by cylindrical inclusions offers a better representation of the effective conductivity of frost as compared to spherical inclusions. Satisfactory agreement between the theory and the measurements for the effective thermal conductivity of frost is demonstrated for a wide range of frost density and frost temperature.
A time-dependent, three-dimensional, multi-ion model of the ionospheric F region at 120-800 km altitude is presented. Account is taken of field-aligned diffusion, cross-field electrodynamic drifts in equatorial and high latitude regions, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions and neutral composition changes. Attention is also given to the effects of ion production by solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, local heating and cooling processes, offsets between the geomagnetic and geographic poles, and bending of field lines near the magnetic equator. The model incorporates all phenomena described by previous models and can be applied to tracing magnetic storm and substorm disturbances from high to low latitudes on a global scale. Sample results are provided for ionospheric features during a June solstice, the solar maximum and in a period of low geomagnetic activity. The model will eventually be used to study coupled ionosphere-thermosphere activity.
Finite laser pulse time effects on temperature wave retardation in flash diffusivity techniques
It is demonstrated that diffusion induced stresses in low resistivity silicon solar cells can significantly reduce both the open-circuit voltage and collection efficiency. The degradation mechanism involves stress induced changes in both the minority carrier mobility and the diffusion length. Thermal recovery characteristics indicate that the stresses are relieved at higher temperatures by divacancy flow (silicon self diffusion). The level of residual stress in as-fabricated cells was found to be negligible in the cells tested.
It is demonstrated that diffusion induced stresses in low resistivity silicon solar cells can significantly reduce both the open-circuit voltage and collection efficiency. The degradation mechanism involves stress induced changes in both the minority carrier mobility and the diffusion length. Thermal recovery characteristics indicate that the stresses are relieved at higher temperatures by divacancy flow (silicon self diffusion). The level of residual stress in as-fabricated cells was found to be negligible in the cells tested.
Presently favored mechanisms of coronal heating (current sheet dissipation and Alfven wave resonant heating) deposit heat in thin layers. Classical thermal conduction cannot explain how heat is transported across the magnetic field. If heating occurs in thin layers, large pressure gradients can be created which can give rise to ballooning modes. These instabilities are caused by the pressure gradient and the curvature of the magnetic field, and are stabilized by magnetic tension. The modes are broad band in wavelength and should produce turbulence. A mixing length expression for the turbulent heat transport shows that it is more than adequate to rapidly convect heat into much broader layers. Furthermore, the turbulent resistivity implies that heating occurs over most of the width of these broadened layers. The broadening also implies that much shorter time scales are required for heating. The beta values in the corona suggest that 1-10 turbulent layers are formed in typical loop or arch structures.