Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

Dynamics and local structure of colossal magnetoresistance manganites

We report Extended X-ray Absorption Fine Structure and anelastic spectroscopy measurements on on hole doped manganese oxides La1-xCaxMnO3 which present the colossal magnetoresistance effect. EXAFS measurements were realized both in the absence and presence of an applied magnetic field of 1.1 Tesla, in a wide temperature range (between 330 and 77 K) and at various dopings (x = 0.25 and x = 0.33). The magnetic field orders the magnetic moments so favouring the electron mobility and the reduction of Mn-O octahedra distortions. We observe the presence of four short and two long Mn-O distances (1.93 and 2.05 Angstrom respectively) above and also below the metal-insulator phase transition. The overall distortion decreases but does not completely disappear in the metallic phase suggesting the possible coexistence of metallic and insulating regions at low temperatures. The magnetic field reduces the lattice distortions showing evidence of a microscopic counterpart of the macroscopic colossal magnetoresistance. We also present preliminary anelastic relaxation spectra in a wide temperature range from 900 K to 1 K on a sample with x = 0.40, in order to study the structural phase transitions and the lattice dynamics. A double peak has been observed at the metal-insulator transition in the imaginary part of Young's modulus. This double peak indicates that the metal-insulator transition could be a more complex phenomenon than a simple second order phase transition. In particular the peak at lower temperatures can be connected with the possible presence of inhomogeneous phase structures. Another intense dissipation peak has been observed corresponding to the structural orthorhombic-trigonal transition around 750 K.

Unsteady MHD forced flow due to a point sink

An analysis is performed to study the unsteady laminar incompressible boundary-layer flow of an electrically conducting fluid in a cone due to a point sink with an applied magnetic field. The unsteadiness in the flow is considered for two types of motion, viz. the motion arising due to the free stream velocity varying continuously with time and the transient motion occurring due to an impulsive change either in the strength of the point sink or in the wall temperature. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The magnetic field increases the skin friction but reduces heat transfer. The heat transfer and temperature field are strongly influenced by the viscous dissipation and Prandtl number. The velocity field is more affected at the early stage of the transient motion, caused by an impulsive change in the strength of the point sink, as compared to the temperature field. When the transient motion is caused by a sudden change in the wall temperature, both skin friction and heat transfer take more time to reach a new steady state. The transient nature of the flow and heat transfer is active for a short time in the case of suction and for a long time in the case of injection. The viscous dissipation prolongs the transient behavior of the flow.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

Unsteady mixed convection with double diffusion over a horizontal cylinder and a sphere within a porous medium

The combined effects of the permeability of the medium, magnetic field, buoyancy forces and dissipation on the unsteady mixed convection flow over a horizontal cylinder and a sphere embedded in a porous medium have been studied. The nonlinear coupled partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The skin friction, heat transfer and mass transfer increase with the permeability of the medium, magnetic field and buoyancy parameter. The heat and mass transfer continuously decrease with the stream-wise distance, whereas the skin friction increases from zero, attains a maximum and then decreases to zero. The skin friction, heat transfer and mass transfer are significantly affected by the free stream velocity distribution. The effect of dissipation parameter is found to be more pronounced on the heat transfer than on the skin friction and mass transfer

Al/SiC carriers for microwave integrated circuits by a new technique of pressureless infiltration

Materials with high thermal conductivity and thermal expansion coefficient matching with that of Si or GaAs are being used for packaging high density microcircuits due to their ability of faster heat dissipation. Al/SiC is gaining wide acceptance as electronic packaging material due to the fact that its thermal expansion coefficient can be tailored to match with that of Si or GaAs by varying the Al:SiC ratio while maintaining the thermal conductivity more or less the same. In the present work, Al/SiC microwave integrated circuit (MIC) carriers have been fabricated by pressureless infiltration of Al-alloy into porous SiC preforms in air. This new technique provides a cheaper alternative to pressure infiltration or pressureless infiltration in nitrogen in producing Al/SiC composites for electronic packaging applications. Al-alloy/65vol% SiC composite exhibited a coefficient of thermal expansion of 7 x 10(-6) K-1 (25 degrees C-100 degrees C) and a thermal conductivity of 147 Wm(-1) K-1 at 30 degrees C. The hysteresis observed in thermal expansion coefficient of the composite in the temperature range 100 degrees C-400 degrees C has been attributed to the presence of thermal residual stresses in the composite. Thermal diffusivity of the composite measured over the temperature range from 30 degrees C to 400 degrees C showed a 55% decrease in thermal diffusivity with temperature. Such a large decrease in thermal diffusivity with temperature could be due to the presence of micropores, microcracks, and decohesion of the Al/SiC interfaces in the microstructure (all formed during cooling from the processing temperature). The carrier showed satisfactory performance after integrating it into a MIC.

Kinetic theory for the density plateau in the granular flow down an inclined plane

The striking lack of observable variation of the volume fraction with height in the center of a granular flow down an inclined plane is analysed using constitutive relations obtained from kinetic theory. It is shown that the rate of conduction in the granular energy balance equation is O(delta(2)) smaller than the rate of production of energy due to mean shear and the rate of dissipation due to inelastic collisions, where the small parameter delta = (d/(1 - e(n))H-1/2), d is the particle diameter, en is the normal coefficient of restitution and H is the thickness of the flowing layer. This implies that the volume fraction is a constant in the leading approximation in an asymptotic analysis in small delta. Numerical estimates of both the parameter delta and its pre-factor are obtained to show that the lack of observable variation of the volume fraction with height can be explained by constitutive relations obtained from kinetic theory.

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

High temperature deformation processing maps for boron modified Ti-6Al-4V alloys

The alloy, Ti-6Al-4V is an alpha + beta Ti alloy that has large prior beta grain size (similar to 2 mm) in the as cast state. Minor addition of B (about 0.1 wt.%) to it refines the grain size significantly as well as produces in-situ TiB needles. The role played by these microstructural modifications on high temperature deformation processing maps of B-modified Ti64 alloys is examined in this paper.Power dissipation efficiency and instability maps have been generated within the temperature range of 750-1000 degrees C and strain rate range of 10(-3)-10(+1) s(-1). Various deformation mechanisms, which operate in different temperature-strain rate regimes, were identified with the aid of the maps and complementary microstructural analysis of the deformed specimens. Results indicate four distinct deformation domains within the range of experimental conditions examined, with the combination of 900-1000 degrees C and 10(-3)-10(-2) s(-1) being the optimum for hot working. In that zone, dynamic globularization of alpha laths is the principle deformation mechanism. The marked reduction in the prior beta grain size, achieved with the addition of B, does not appear to alter this domain markedly. The other domains, with negative values of instability parameter, show undesirable microstructural features such as extensive kinking/bending of alpha laths and breaking of beta laths for Ti64-0.0B as well as generation of voids and cracks in the matrix and TiB needles in the B-modified alloys. (C) 2010 Elsevier B.V. All rights reserved.

Phase-shift randomness in one-dimensional disordered conductors in the quasi-metallic domain

An invariant imbedding method yields exact analytical results for the distribution of the phase theta (L) of the reflection amplitude and for low-order resistance moments (pn) for a disordered conductor of length L in the quasi-metallic regime L<momentum. For realistic values of 2k0L, the resistance moments are dominated by the terms obtained previously by arbitrarily assuming uniformly distributed phases. A direct proof of the validity of the random-phase assumption for studying scaling of resistance in one dimension is thus obtained.

Unsteady three-dimensional boundary layer flow over a flat plate with attached cylinder

The unsteady laminar incompressible three-dimensional boundary layer flow and heat transfer on a flat plate with an attached cylinder have been studied when the free stream velocity components and wall temperature vary inversely as linear and quadratic functions of time, respectively. The governing semisimilar partial differential equations with three independent variables have been solved numerically using a quasilinear finite-difference scheme. The results indicate that the skin friction increases with parameter ? which characterizes the unsteadiness in the free stream velocity and the streamwise distance Image , but the heat transfer decreases. However, the skin friction and heat transfer are found to change little along Image . The effect of the Prandtl number on the heat transfer is found to be more pronounced when ? is small, whereas the effect of the dissipation parameter is more pronounced when ? is comparatively large.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.

Brownian motion of spins; generalized spin Langevin equation

We derive the Langevin equations for a spin interacting with a heat bath, starting from a fully dynamical treatment. The obtained equations are non-Markovian with multiplicative fluctuations and concommitant dissipative terms obeying the fluctuation-dissipation theorem. In the Markovian limit our equations reduce to the phenomenological equations proposed by Kubo and Hashitsume. The perturbative treatment on our equations lead to Landau-Lifshitz equations and to other known results in the literature.

Characterization of hot deformation behavior of brasses using processing maps: Part II. β Brass and α-β brass

The hot deformation behaviors of β brass in the temperature range of 550°C to 800°C and α-β brass in the temperature range of 450°C to 800°C have been characterized in the strain rate range of 0.001 to 100 s−1 using processing maps developed on the basis of the Dynamic Materials Model. The map for β brass revealed a domain of superplasticity in the entire temperature range and at strain rates lower than 1 s−1, with a maximum efficiency of power dissipation of about 68 pct. The temperature variation of the efficiency of power dissipation in the domain is similar to that of the diffusion coefficient for zinc in β brass, confirming that the diffusion-accommodated flow controls the superplasticity. The material undergoes microstructural instability in the form of adiabatic shear bands and strain markings at temperatures lower than 700°C and at strain rates higher than 10 s−1. The map for α-β brass revealed a wide domain for processing in the temperature range of 550°C to 800°C and at strain rates lower than 1 s−1, with a maximum efficiency of 54 pct occurring at about 750°C and 0.001 s−1. In the domain, the α phase undergoes dynamic recrystallization and controls the hot deformation of the alloy, while the β phase deforms superplastically. At strain rates greater than 1 s−1, α-β brass exhibits microstructural instabilities manifested as flow rotations at lower temperatures and localized shear bands at higher temperatures.

Effect of stacking fault energy on the dynamic recrystallization during hot working of FCC metals: A study using processing maps

The influence of stacking fault energy (SFE) on the mechanism of dynamic recrystallization (DRX) during hot deformation of FCC metals is examined in the light of results from the power dissipation maps. The DRX domain for high SFE metals like Al and Ni occurred at homologous temperature below 0·7 and strain rates of 0·001 s−1 while for low SFE metals like Cu and Pb the corresponding values are higher than 0·8 and 100 s−1. The peak efficiencies of power dissipation are 50% and below 40% respectively. A simple model which considers the rate of interface formation (nucleation) involving dislocation generation and simultaneous recovery and the rate of interface migration (growth) occurring with the reduction in interface energy as the driving force, has been proposed to account for the effect of SFE on DRX. The calculations reveal that in high SFE metals, interface migration controls DRX while the interface formation is the controlling factor in low SFE metals. In the latter case, the occurrence of flow softening and oscillations could be accounted for by this model.