Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 775 records · Page 43

Long-Term Cyclic Oxidation Behavior of Uncoated and Coated Re-108 and In-939 at 980 and 870 C

Very long-term cyclic oxidation behavior of Re108 and In939 with and without a protective coating was evaluated at 980 and 870 C, respectively. Re-108 and In-939 without a protective coating began to show a rapid weight loss at 3000 h due to scale spallation, indicating the need for an oxidation protective coating for longer than thousands of hours of oxidative life. NiAl-base coatings of a vapor phase aluminide (VPA), a pack aluminide (CODEP), and a slurry paint aluminide (SERMALOY J) were applied on Re-108 and In-939. The VPA and CODEP on Re-108 and all three coatings on In-939 showed excellent cyclic oxidation resistance out to 10,000 h. Coated alloys were annealed in an inert atmosphere to determine the loss of Al from the coating into the alloy substrate through diffusion. The Al loss from the coating through diffusion was twice as great as the Al loss through oxidation after 10,000 h of cyclic exposure. The oxidation life of VPA-coated Re-108 was estimated by calculating the amount of Al initially available for protective oxidation and the amount of Al lost through oxidation and diffusion.

Lee, K. N.↗

Small Scale Polygons and the History of Ground Ice on Mars

Recent progress on polygon modeling has focused on the diameter and surface relief that we expect of thermal-contraction polygons in martian permafrost. With this in mind, we developed a finite-element model of thermal-contraction-crack behavior in permafrost in a martian climate. This model was generated from a finite element code by Jay Melosh (called TECTON) originally developed for terrestrial and planetary crustal-deformation studies. We adapted this model to martian permafrost by including time (and temperature) dependent rheologies, boundary conditions, and isotropic thermal-contraction, as well as several small adaptations to a martian environment. We tested our model extensively, including comparison to an analytic solution of pre-fracture stress. We recently published an analysis of two potential sources of water for forming the recent gullies. In this work we first evaluated the potential for near-surface ground ice (in the top meter or so of soil) to melt under conditions of solar heating on sloped surfaces at high obliquity, utilizing both thermal and diffusion-based ground-ice-stability models; our results suggested that the ground ice will sublimate, and the ice table will recede to greater depths before the melting temperature can be reached. An exception can occur only for extremely salt-rich ice, depressing the freezing point.

Mellon, Michael T.↗

Time-dependent modulation of galactic cosmic rays by merged interaction regions

Models that solve the one-dimensional, solar modulation equation have reproduced the 11-year galactic cosmic ray using functional representations of global merged interaction regions (MIRs). This study extends those results to the solution of the modulation equation with explicit time dependence. The magnetometers on Voyagers 1 and 2 provide local magnetic field intensities at regular intervals, from which one calculates the ratio of the field intensity to the average local field. These ratios in turn are inverted to form diffusion coefficients. Strung together in radius and time, these coefficents then fall and rise with the strength of the interplanetary magnetic field, becoming representations of MIRs. These diffusion coefficients, calculated locally, propagate unchanged from approx. 10 AU to the outer boundary (120 AU). Inside 10 AU, all parameters, including the diffusion coefficient are assumed constant in time and space. The model reproduces the time-intensity profiles of Voyager 2 and Pioneer 10. Radial gradient data from 1982-1990 between Pioneer 10 and Voyager 2 are about the same magnitude as those calculated in the model. It is also shows agreement in rough magnitude with the radial gradient between Pioneer 10 and 1 AU. When coupled with enhanced, time-dependent solar wind speed at the probe's high latitude, as measured by independent observers, the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source for the observed negative latitudinal gradients. The model exhibits the 11-year cyclical cosmic ray intensity behavior at all radii, including 1 AU, not just at the location of the spacecraft where the magnetic fields are measured. In addition, the model's point of cosmic ray maximum correctly travels at the solar wind speed, illustrating the well-known propagation of modulation. Finally, at least in the inner heliosphere this model accounts for the delay experienced by lower-rigidity protons in reaching their time-intensity peak. The actual delays in this model, however, are somewhat smaller than the data. In the outer heliosphere the models sees no delays, and the data are ambiguous as to their existence. It appears that strong magnetic field compression regions (merged interaction regions) that are 3-4 times the average field strength can, at least in a helioequatorial band, disrupt effects, such as drifts, that could dominate in quieter magnetic fields. The question remains: Is the heliosphere ever quiet enough to allow such effects to be unambiguously measured, at least in the midlatitudes?

Perko, J. S.↗

Structures Self-Assembled Through Directional Solidification

Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.

Dynys, Frederick W.↗

Improved Detection System Description and New Method for Accurate Calibration of Micro-Channel Plate Based Instruments and Its Use in the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The Fast Plasma Investigation (FPI) on NASAs Magnetospheric MultiScale (MMS) mission employs 16 Dual Electron Spectrometers (DESs) and 16 Dual Ion Spectrometers (DISs) with 4 of each type on each of 4 spacecraft to enable fast (30 ms for electrons; 150 ms for ions) and spatially differentiated measurements of the full 3D particle velocity distributions. This approach presents a new and challenging aspect to the calibration and operation of these instruments on ground and in flight. The response uniformity, the reliability of their calibration and the approach to handling any temporal evolution of these calibrated characteristics all assume enhanced importance in this application, where we attempt to understand the meaning of particle distributions within the ion and electron diffusion regions of magnetically reconnecting plasmas. Traditionally, the micro-channel plate (MCP) based detection systems for electrostatic particle spectrometers have been calibrated using the plateau curve technique. In this, a fixed detection threshold is set. The detection system count rate is then measured as a function of MCP voltage to determine the MCP voltage that ensures the count rate has reached a constant value independent of further variation in the MCP voltage. This is achieved when most of the MCP pulse height distribution (PHD) is located at higher values (larger pulses) than the detection system discrimination threshold. This method is adequate in single-channel detection systems and in multi-channel detection systems with very low crosstalk between channels. However, in dense multi-channel systems, it can be inadequate. Furthermore, it fails to fully describe the behavior of the detection system and individually characterize each of its fundamental parameters. To improve this situation, we have developed a detailed phenomenological description of the detection system, its behavior and its signal, crosstalk and noise sources. Based on this, we have devised a new detection system calibration method that enables accurate and repeatable measurement and calibration of MCP gain, MCP efficiency, signal loss due to variation in gain and efficiency, crosstalk from effects both above and below the MCP, noise margin, and stability margin in one single measurement. More precise calibration is highly desirable as the instruments will produce higher quality raw data that will require less post-acquisition data correction using results from in-flight pitch angle distribution measurements and ground calibration measurements. The detection system description and the fundamental concepts of this new calibration method, named threshold scan, will be presented. It will be shown how to derive all the individual detection system parameters and how to choose the optimum detection system operating point. This new method has been successfully applied to achieve a highly accurate calibration of the DESs and DISs of the MMS mission. The practical application of the method will be presented together with the achieved calibration results and their significance. Finally, it will be shown that, with further detailed modeling, this method can be extended for use in flight to achieve and maintain a highly accurate detection system calibration across a large number of instruments during the mission.

detection system↗

Phase Transformation and Creep Behavior in Ti50Pd30Ni20 High Temperature Shape Memory Alloy in Compression

The creep behavior and the phase transformation of Ti50Pd30Ni20 High Temperature Shape Memory Alloy (HTSMA) is investigated by standard creep tests and thermomechanical tests. Ingots of the alloy are induction melted, extruded at high temperature, from which cylindrical specimens are cut and surface polished. A custom high temperature test setup is assembled to conduct the thermomechanical tests. Following preliminary monotonic tests, standard creep tests and thermally induced phase transformation tests are conducted on the specimen. The creep test results suggest that over the operating temperatures and stresses of this alloy, the microstructural mechanisms responsible for creep change. At lower stresses and temperatures, the primary creep mechanism is a mixture of dislocation glide and dislocation creep. As the stress and temperature increase, the mechanism shifts to predominantly dislocation creep. If the operational stress or temperature is raised even further, the mechanism shifts to diffusion creep. The thermally induced phase transformation tests show that actuator performance can be affected by rate independent irrecoverable strain (transformation induced plasticity + retained martensite) as well as creep. The rate of heating and cooling can adversely impact the actuators performance. While the rate independent irrecoverable strain is readily apparent early in the actuators life, viscoplastic strain continues to accumulate over the lifespan of the HTSMA. Thus, in order to get full actuation out of the HTSMA, the heating and cooling rates must be sufficiently high enough to avoid creep.

Kumar, Parikshith K.↗

Performance, defect behavior and carrier enhancement in low energy, proton irradiated p+nn+ InP solar cells

The highest AMO efficiency (19.1 percent) InP solar cell consisted of an n+pp+ structure epitaxially grown on a p+ InP substrate. However, the high cost and relative fragility of InP served as motivation for research efforts directed at heteroepitaxial growth of InP on more viable substrates. The highest AMO efficiency (13.7 percent) for this type of cell was achieved using a GaAs substrate. Considering only cost and fracture toughness, Si would be the preferred substrate. The fact that Si is a donor in InP introduces complexities which are necessary in order to avoid the formation of an efficiency limiting counterdiode. One method used to overcome this problem lies in employing an n+p+ tunnel junction in contact with the cell's p region. A simpler method consists of using an n+ substrate and processing the cell in the p+ nn+ configuration. This eliminates the need for a tunnel junction. Unfortunately, the p/n configuration has received relatively little attention the best cell with this geometry having achieved an efficiency of 17 percent. Irradiation of these homoepitaxial cells, with 1 Mev electrons, showed that they were slightly more radiation resistant than diffused junction n/p cells. Additional p/n InP cells have been processed by some activity aimed at diffusion. Currently, there has been some activity aimed at producing heteroepitaxial p+nn+ InP cells using n+ Ge substrates. Since, like Si, Ge is an n-dopant in InP, use of this configuration obviates the need for a tunnel junction. Obviously, before attempting to process heteroepitaxial cells, one must produce a reasonably good homoepitaxial cell. In the present case we focus our attention on homoepitaxially on an n+ Ge substrate.

Weinberg, I.↗

Optical studies in the holographic ground station

The Holographic Group System (HGS) Facility in rooms 22 & 123, Building 4708 has been developed to provide for ground based research in determining pre-flight parameters and analyzing the results from space experiments. The University of Alabama, Huntsville (UAH) has researched the analysis aspects of the HGS and reports their findings here. Some of the results presented here also occur in the Facility Operating Procedure (FOP), which contains instructions for power up, operation, and powerdown of the Fluid Experiment System (FES) Holographic Ground System (HGS) Test Facility for the purpose of optically recording fluid and/or crystal behavior in a test article during ground based testing through the construction of holograms and recording of videotape. The alignment of the optical bench components, holographic reconstruction and and microscopy alignment sections were also included in the document for continuity even though they are not used until after optical recording of the test article) setup of support subsystems and the Automated Holography System (AHS) computer. The HGS provides optical recording and monitoring during GCEL runs or development testing of potential FES flight hardware or software. This recording/monitoring can be via 70mm holographic film, standard videotape, or digitized images on computer disk. All optical bench functions necessary to construct holograms will be under the control of the AHS personal computer (PC). These include type of exposure, time intervals between exposures, exposure length, film frame identification, film advancement, film platen evacuation and repressurization, light source diffuser introduction, and control of realtime video monitoring. The completed sequence of hologram types (single exposure, diffuse double exposure, etc.) and their time of occurrence can be displayed, printed, or stored on floppy disk posttest for the user.

Workman, Gary L.↗

Behavior of Particle Depots in Molten Silicon During Float-Zone Growth in Strong Magnetic Fields

Solar cells made from directionally solidified silicon cover 57% of the photovoltaic industry's market [1]. One major issue during directional solidification of silicon is the precipitation of foreign phase particles. These particles, mainly SiC and Si3N4, are precipitated from the dissolved crucible coating, which is made of silicon nitride, and the dissolution of carbon monoxide from the furnace atmosphere. Due to their hardness and size of several hundred micrometers, those particles can lead to severe problems during the wire sawing process for wafering the ingots. Additionally, SiC particles can act as a shunt, short circuiting the solar cell. Even if the particles are too small to disturb the wafering process, they can lead to a grit structure of silicon micro grains and serve as sources for dislocations. All of this lowers the yield of solar cells and reduces the performance of cells and modules. We studied the behaviour of SiC particle depots during float-zone growth under an oxide skin, and strong static magnetic fields. For high field strengths of 3T and above and an oxide layer on the sample surface, convection is sufficiently suppressed to create a diffusive like regime, with strongly dampened convection [2, 3]. To investigate the difference between atomically rough phase boundaries and facetted growth, samples with [100] and [111] orientation were processed.

Jauss, T.↗

Spatially Complete Global Surface Albedos Derived from Terra/MODIS Data

Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. Recent production of land surface anisotropy, diffuse bihemispherical (white-sky) albedo and direct beam directional hemispherical (black-sky) albedo from observations acquired by the MODIS instruments aboard NASA s Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface's radiative characteristics. Cloud cover, which cutails retrievals, and the presence of ephemeral and seasonal snow limit the snow-free data to approximately half the global land surfaces on an annual equal-angle basis. This precludes the MOD43B3 albedo products from being used in some remote sensing and ground-based applications, climate models, and global change research projects. An ecosystem-dependent temporal interpolation technique is described that has been developed to fill missing or seasonally snow-covered data in the official MOD43B3 albedo product. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data. The resulting snow-free value-added products provide the scientific community with spatially and temporally complete global white- and black-sky surface albedo maps and statistics. These products are stored on 1'(approximately 10 km) and coarser resolution equal-angle grids, and are computed for the first seven MODIS wavelengths, ranging from 0.47 through 2.1 microns, and for three broadband wavelengths, 0.3-0.7,0.3-5.0 and 0.7-5.0 microns.

King, Michael D.↗

A comparative study of Full Navier-Stokes and Reduced Navier-Stokes analyses for separating flows within a diffusing inlet S-duct

A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced Navier Stokes (RNS) initial value space marching solution technique has been applied to a class of separated flow problems within a diffusing S-duct configuration characterized by vortex-liftoff. Both the FNS and the RNS solution technique were able to capture the overall flow physics of vortex lift-off, and gave remarkably similar results which agreed reasonably well with the experimental measured averaged performance parameters of engine face total pressure recovery and distortion. However, the Full Navier-Stokes and Reduced Navier-Stokes also consistently predicted separation further downstream in the M2129 inlet S-duct than was indicated by experimental data, thus compensating errors were present in the two Navier-Stokes analyses. The difficulties encountered in the Navier-Stokes separations analyses of the M2129 inlet S-duct center primarily on turbulence model issues, and these focused on two distinct but different phenomena, namely, (1) characterization of low skin friction adverse pressure gradient flows, and (2) description of the near wall behavior of flows characterized by vortex lift-off.

Anderson, B. H.↗

Annealing of irradiated n+p InP buried homojunctions

At the last SPRAT conference, the Naval Research Laboratory (NRL) presented results from two experiments. One studied n+p diffused junction (DJ) InP solar cells, and the other studied n+p shallow homojunction (SHJ) InP mesa diodes grown by metalorganic chemical vapor deposition (MOCVD). The former work showed that a DJ solar cell in which the maximum power P(sub max) had been degraded by nearly 80 percent under irradiation recovered completely under short circuit illumination at 450K. The recovery was accompanied by the removal of all but one of the radiation-induced defect levels. The latter work, on the other hand, showed that the radiation-induced defects in the SHJ diodes did not anneal until the temperature reached 650K. These results suggest that an irradiated DJ solar cell, under illumination, will anneal at a temperature 200K lower than an irradiated SHJ cell. This is an unexpected result considering the similarity of the devices. The goal of the present research is to explain this different behavior. This paper investigates two points which arose from the previous studies. The first point is that the DJ cells were annealed under illumination while the SHJ diodes were annealed without bias. The second point investigated here is that the emitters of the DJ and SHJ devices were significantly different.

Walters, Robert J.↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Elevated temperature deformation of thoria dispersed nickel-chromium

The deformation behavior of thoria nickel-chromium (TD-NiCr) was examined over the temperature range 593 C (1100 F) to 1260 C (2300 F) in tension and compression and at 1093 C (2000 F) in creep. Major emphasis was placed on: (1) the effects of the material and test related variables (grain size, temperature, stress and strain rate) on the deformation process; and (2) the evaluation of single crystal TD-NiCr material produced by a directional recrystallization process. Elevated temperature yield strength levels and creep activation enthalpies were found to increase with increasing grain size reaching maximum values for the single crystal TD-NiCr. Stress exponent of the steady state creep rate was also significantly higher for the single crystal TD-NiCr as compared to that determined for the polycrystalline materials. The elevated temperature deformation of TD-NiCr was analyzed in terms of two concurrent, parallel processes: diffusion controlled grain boundary sliding, and dislocation motion.

Kane, R. D.↗

Evidence for Widespread Cooling in an Active Region Observed with the SDO Atmospheric Imaging Assembly

A well known behavior of EUV light curves of discrete coronal loops is that the peak intensities of cooler channels or spectral lines are reached at progressively later times. This time lag is understood to be the result of hot coronal loop plasma cooling through these lower respective temperatures. However, loops typically comprise only a minority of the total emission in active regions. Is this cooling pattern a common property of active region coronal plasma, or does it only occur in unique circumstances, locations, and times? The new SDO/AIA data provide a wonderful opportunity to answer this question systematically for an entire active region. We measure the time lag between pairs of SDO/AIA EUV channels using 24 hours of images of AR 11082 observed on 19 June 2010. We find that there is a time-lag signal consistent with cooling plasma, just as is usually found for loops, throughout the active region including the diffuse emission between loops for the entire 24 hour duration. The pattern persists consistently for all channel pairs and choice of window length within the 24 hour time period, giving us confidence that the plasma is cooling from temperatures of greater than 3 MK, and sometimes exceeding 7 MK, down to temperatures lower than approx. 0.8 MK. This suggests that the bulk of the emitting coronal plasma in this active region is not steady; rather, it is dynamic and constantly evolving. These measurements provide crucial constraints on any model which seeks to describe coronal heating.

Viall, Nicholeen M.↗

On the Importance of the Dynamics of Discretizations

It has been realized recently that the discrete maps resulting from numerical discretizations of differential equations can possess asymptotic dynamical behavior quite different from that of the original systems. This is the case not only for systems of Ordinary Differential Equations (ODEs) but in a more complicated manner for Partial Differential Equations (PDEs) used to model complex physics. The impact of the modified dynamics may be mild and even not observed for some numerical methods. For other classes of discretizations the impact may be pronounced, but not always obvious depending on the nonlinear model equations, the time steps, the grid spacings and the initial conditions. Non-convergence or convergence to periodic solutions might be easily recognizable but convergence to incorrect but plausible solutions may not be so obvious - even for discretized parameters within the linearized stability constraint. Based on our past four years of research, we will illustrate some of the pathology of the dynamics of discretizations, its possible impact and the usage of these schemes for model nonlinear ODEs, convection-diffusion equations and grid adaptations.

Sweby, Peter K.↗

Solidification (crystal growth) in the presence of gravitational forces

The surface tension behavior of doped and undoped InSb melts was investigated as well as their temperature and composition dependence. Surface tension in InSb melts was determined using the sessile-drop technique covering the temperature range from 530 C to 880 C. A linear regression of the data obtained shows that the temperature dependence of sigma is 392- (T-530) x (7000) plus or minus 10 dyne/cm. The d sigma/d Tau for intrinsics InSb is less than that previously reported. On the basis of the surface tension data obtained, it is concluded that surface tension induced convective flow velocities in InSb under reduced gravity conditions range from zero to at most 1 cm/sec. Accordingly, no interference with dopant segregation can be expected during growth in space because the momentum boundary layer (at the crystal melt interface) associated with any Marangoni-type convective flows would, at the given growth rate, be significantly larger than the predicted diffusion boundary layer thickness.

Gatos, H. C.↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗