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.

Evaluation of ductile-brittle transition temperature (DBTT) of aluminide bond coats by micro-tensile test method

A method has been suggested to accurately determine the DBTT of diffusion aluminide bond coats. Micro-tensile testing of free-standing coating samples has been carried out. The DBTT was determined based on the variation of plastic strain-to-fracture with temperature. The positive features of this method over the previously reported techniques are highlighted. (C) 2010 Elsevier B.V. All rights reserved.

Mutual diffusion coefficients of some binary liquid systems: benzene-n-alkyl alcohol

The mutual diffusion coefficients for binary liquid systems of benzene-n-alkyl alcohol at various compositions have been determined by the diaphragm cell method at 28-degrees-C. The alcohols used were the members of n-paraffinic alcohols ranging from C1 to C8. The maximum possible experimental error is 14%. The data were fitted with a generalized correlation, giving the deviation from the experimental data to within 2.75%, on average.

Protonic conductivity of (NH4) 4Fe(CN)6·1.5H2O by complex admittance method

The protonic conductivity of ammonium ferrocyanide hydrate has been studied by the complex admittance method. The admittance plots show departures from ideal Debye behaviour. The values of ionic conductivity (sigma = 3.7 X 10(-5) (OMEGA-cm)-1) and diffusion coefficient (D = 3.8 X 10(-10) cm2/s) obtained at room temperature are consistent with the corresponding values estimated by an earlier NMR study.

Impurity diffusion in bismuth single-crystals

Measurements of impurity diffusion of 86Rb, 90Sr, 133Ba, and 137Cs in single crystal Bi were carried out. Diffusion samples were prepared from single crystal Bi by ion implantation. About 1012-1013 ions were implanted, resulting in surface activities approx =104 cpm. After implantation, specimens were annealed for specified times at 220-265 deg C, and tracer penetration profiles were determined by an electrolytic method. A typical penetration profile for 137Cs in Bi showed a linear relationship for log C vs x in with Fick's law for volume diffusion. Laws of grain boundary diffusion were not obeyed and the order of magnitude of the penetration distances was much less than on a grain boundary mechanism. Results were interpreted in terms of a modified Fischer analysis using a kinetic trapping term. Effective half lengths for trapping at a twin boundary were determined for each impurity.

Digital simulation of galvanostatic current-potential data for gas-diffusion electrodes and estimation of electrodekinetic parameters

A computerized non-linear-least-squares regression procedure to analyse the galvanostatic current-potential data for kinetically hindered reactions on porous gas-diffusion electrodes is reported. The simulated data fit well with the corresponding measured values. The analytical estimates of electrode-kinetic parameters and uncompensated resistance are found to be in good agreement with their respective values obtained from Tafel plots and the current-interrupter method. The procedure circumvents the need to collect the data in the limiting-current region where the polarization values are usually prone to errors. The polarization data for two typical cases, namely, methanol oxidation on a carbon-supported platinum-tin electrode and oxygen reduction on a Nafion-coated platinized carbon electrode, are successfully analysed.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Unsteady free convection flow over an infinite vertical porous plate due to the combined effects of thermal and mass diffusion, magnetic field and Hall currents

The unsteady free convection flow over an infinite vertical porous plate, which moves with time-dependent velocity in an ambient fluid, has been studied. The effects of the magnetic field and Hall current are included in the analysis. The buoyancy forces arise due to both the thermal and mass diffusion. The partial differential equations governing the flow have been solved numerically using both the implicit finite difference scheme and the difference-differential method. For the steady case, analytical solutions have also been obtained. The effect of time variation on the skin friction, heat transfer and mass transfer is very significant. Suction increases the skin friction coefficient in the primary flow, and also the Nusselt and Sherwood numbers, but the skin friction coefficient in the secondary flow is reduced. The effect of injection is opposite to that of suction. The buoyancy force, injection and the Hall parameter induce an overshoot in the velocity profiles in the primary flow which changes the velocity gradient from a negative to a positive value, but the magnetic field and suction reduce this velocity overshoot.

A Reliable method of obtaining Breakthrough Curves of ions in Soils using Transport Equation

Molecular diffusion plays a dominant role in transport of contaminants through fine-grained soils with low hydraulic conductivity. Attenuation processes occur while contaminants travel through the soils. Effective diffusion coefficient (De) is expected to take into consideration various attenuation processes. Effective diffusion coefficient has been considered to develop a general approach for modelling of contaminant transport in soils.The effective diffusion coefficient of sodium in presence of sulphate has been obtained using the column test.The reliability of De, has been checked by comparing theoretical breakthrough curves of sodium ion in soils obtained using advection diffusion equation with the experimental curve.

Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.

Growth and consumption rates of the phase layers during interdiffusion in a diffusion couple with finite end member

The relations for the growth and consumption rates of a layer with finite thickness as an end member and the product phases in the interdiffusion zone are developed. We have used two different methodologies, the diffusion based and the physico-chemical approach to develop the same relations. We have shown that the diffusion based approach is rather straightforward; however, the physico-chemical approach is much more versatile than the other method. It was found that the position of the marker plane becomes vague in the second stage of the interdiffusion process in pure A thin layer/B couple, where two phases grow simultaneously.

Self-diffusion of colloids at freezing

This is an introduction to the theory of interacting Brownian particles, as applied to charge-stabilised colloidal suspensions near their equilibrium liquid-solid transition. The density functional approach to the statics of the transition is reviewed briefly, and the generalised Langevin equation method for the dynamics presented in detail. Work with A.V. Indrani [1] on a self-consistent approach for calculating the excess single-particle friction is presented, which explains the observed [2] ''universal'' suppression of self-diffusion at freezing as a consequence of the universal structure-factor height at this transition. Criticisms, open questions, and challenges to theory are discussed.

Shear-enhanced diffusion in charge-stabilised colloids

The generalised Langevin equation method for the dynamics of interacting colloids presented in my previous lecture is extended here to the case of a sheared suspension. A calculation of shear-dependent diffusivities using these methods is found to account for puzzling observations in experiments and simulations. The limitations of the method are discussed, and important unresolved questions presented. This lecture summarises work done in collaboration with A.V. Indrani [1].

Reactive Interdiffusion: A Method for Making Composition Graded Metal-ceramic Composites

Presented is a new method for making composition graded metal-ceramic composites using reactive inter-diffusion between a metal and a complex ceramic. Composition variation in both metal and ceramic phases with distance along the direction of diffusion is achieved. The design criteria for developing such composites are discussed. The system should exhibit extensive solid solubility in both metallic and ceramic phases, a defined gradation in the stabilities of the oxides, and mobility of electrons or holes in the oxide solid solution. The complex ceramic used for making the composite should be polycrystalline with sufficient porosity to accommodate the volume expansion caused by alloy precipitation. An inert atmosphere to prevent oxidation and high processing temperature to facilitate diffusive transport are required. The process is illustrated using the reaction couples Fe-NiTiO3, Fe-(Mg,Co)TiO3 and Fe-(Ni,Co)TiO3.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.