Numerical modeling of heat and mass transfer in laser surface alloying: Non-equilibrium solidification effects

A transient macroscopic model is developed for studying heat and mass transfer in a single-pass laser surface alloying process, with particular emphasis on non-equilibrium solidification considerations. The solution for species concentration distribution requires suitable treatment of non-equilibrium mass transfer conditions. In this context, microscopic features pertaining to non-equilibrium effects on account of solutal undercooling are incorporated through the formulation of a modified partition-coefficient. The effective partition-coefficient is numerically modeled by Means of a number of macroscopically observable parameters related to the solidifying domain. The numerical model is so developed that the modifications on account of non-equilibrium solidification considerations can be conveniently implemented in existing numerical codes based on equilibrium solidification considerations.

Symmetry-breaking transitions in equilibrium shapes of coherent precipitates: Effect of elastic anisotropy and inhomogeneity

We examine the symmetry-breaking transitions in equilibrium shapes of coherent precipitates in two-dimensional (2-D) systems under a plane-strain condition with the principal misfit strain components epsilon(xx)*. and epsilon(yy)*. For systems with cubic elastic moduli, we first show all the shape transitions associated with different values of t = epsilon(yy)*/epsilon(xx)*. We also characterize each of these transitions, by studying its dependence on elastic anisotropy and inhomogeneity. For systems with dilatational misfit (t = 1) and those with pure shear misfit (t = -1), the transition is from an equiaxed shape to an elongated shape, resulting in a break in rotational symmetry. For systems with nondilatational misfit (-1 < t < 1; t not equal 0), the transition involves a break in mirror symmetries normal to the x- and y-axes. The transition is continuous in all cases, except when 0 < t < 1. For systems which allow an invariant line (-1 less than or equal to t < 0), the critical size increases with an increase in the particle stiffness. However, for systems which do not allow an invariant line (0 < t less than or equal to 1), the critical size first decreases, reaches a minimum, and then starts increasing with increasing particle stiffness; moreover, the transition is also forbidden when the particle stiffness is greater than a critical value.

Equilibrium and relaxation in turbulent wakes

In order to study the memory of the larger eddies in turbulent shear flow, experiments have been conducted on plane turbulent wakes undergoing transition from an initial (carefully prepared) equilibrium state to a different final one, as a result of a nearly impulsive pressure gradient. It is shown that under the conditions of the experiments the equations of motion possess self-preserving solutions in the sense of Townsend (1956), but the observed behaviour of the wake is appreciably different when the pressure gradient is not very small, as the flow goes through a slow relaxation process before reaching final equilibrium. Measurements of the Reynolds stresse show that the approach to a new equilibrium state is exponential, with a relaxation length of the order of 103 momentum thicknesses. It is suggested that a flow satisfying the conditions required by a self-preservation analysis will exhibit equilibrium only if the relaxation length is small compared with a characteristic streamwise length scale of the flow.

Turbulence measurements in an ‘equilibrium’ axisymmetric wall jet Journal

This paper reports measurements of turbulent quantities in an axisymmetric wall jet subjected to an adverse pressure gradient in a conical diffuser, in such a way that a suitably defined pressure-gradient parameter is everywhere small. Self-similarity is observed in the mean velocity profile, as well as the profiles of many turbulent quantities at sufficiently large distances from the injection slot. Autocorrelation measurements indicate that, in the region of turbulent production, the time scale of ν fluctuations is very much smaller than the time scale of u fluctuations. Based on the data on these time scales, a possible model is proposed for the Reynolds stress. One-dimensional energy spectra are obtained for the u, v and w components at several points in the wall jet. It is found that self-similarity is exhibited by the one-dimensional wavenumber spectrum of $\overline{q^2}(=\overline{u^2}+\overline{v^2}+\overline{w^2})$, if the half-width of the wall jet and the local mean velocity are used for forming the non-dimensional wavenumber. Both the autocorrelation curves and the spectra indicate the existence of periodicity in the flow. The rate of dissipation of turbulent energy is estimated from the $\overline{q^2}$ spectra, using a slightly modified version of a previously suggested method.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Three-dimensional double-diffusive convection and macrosegregation during non-equilibrium solidification of binary mixtures

A three-dimensional transient mathematical model (following a fixed-grid enthalpy-based continuum formulation) is used to study the interaction of double-diffusive natural convection and non-equilibrium solidification of a binary mixture in a cubic enclosure cooled from a side. Investigations are carried out for two separate test systems, one corresponding to a typical model "metal-alloy analogue" system and other corresponding to a real metal-alloy system. Due to stronger effects of solutal buoyancy in actual metal-alloy systems than in corresponding analogues, the convective transport mechanisms for the two cases are quite different. However, in both cases, similar elements of three-dimensionality are observed in the curvature and spacing of the projected streamlines. As a result of three-dimensional convective flow patterns, a significant solute macrosegregation is observed across the transverse sections of the cavity, which cannot be captured by two-dimensional simulations. (C) 2003 Elsevier Science Ltd. All rights reserved.

Closed form solutions for element equilibrium and flexibility matrices of eight node rectangular plate bending element using integrated force method

Closed form solutions for equilibrium and flexibility matrices of the Mindlin-Reissner theory based eight-node rectangular plate bending element (MRP8) using integrated Force Method (IFM) are presented in this paper. Though these closed form solutions of equilibrium and flexibility matrices are applicable to plate bending problems with square/rectangular boundaries, they reduce the computational time significantly and give more exact solutions. Presented closed form solutions are validated by solving large number of standard square/rectangular plate bending benchmark problems for deflections and moments and the results are compared with those of similar displacement-based eight-node quadrilateral plate bending elements available in the literature. The results are also compared with the exact solutions.

Activities and their relation to cation distribution in NiAl2O4 MgAl2O4 spinel solid solutions

The activity of NiAl2O4 in NiAl2O4MgAl2O4 solid solutions has been measured by using a solid oxide galvanic cell of the type, Pt, Ni + NiAl2O4 + Al2O3(α)/CaOZrO2/Ni + NixMg1−xAl2O4 + Al2O3(α). Pt, in the temperature range 750–1150°C. The activities in the spinel solid solutions show negative deviations from Raoult's law. The cation distribution in the solid solutions has been calculated using site preference energies independent of composition for Ni2+, Mg2+, and Al3+ ions obtained from crystal field theory and measured cation disorder in pure NiAl2O4 and MgAl2O4, and assumi g ideal mixing of cations on the tetrahedral and octahedral positions. The calculated values correctly predict the decrease in the fraction, α, of Ni2+ ions on tetrahedral sites for 1>x>0.25, observed by Porta et al. [J. Solid State Chem.11, 135 (1974)] but do not support their tentative evidence for an increase in α for x < 0.25. The measured excess free energy of mixing can be completely accounted for by using either the calculated or the measured cation distributions. This suggests that the Madelung energy is approximately a linear function of composition in the solid solutions. The composition of NiOMgO solid solutions in equilibrium with NiAl2O4MgAl2O4 solid solutions has been calculated from the results and information available in literature.

Thermodynamic study of Fe2O3 Fe2(SO4)3 equilibrium using an oxyanionic electrolyte (Na2SO4 I)

The emf of the cell, Pt, Ar + O2 + SO2 + SO3/Na2SO4-I/Fe2O2 + Fe2(SO4)3, Pt, has been measured in the temperature range 800 to 1000 K, using a gas mixture of known input composition as the reference electrode. The equilibrium composition of the reference gas at the measuring temperatures was computed using the thermodynamic data on the gaseous species reported in the literature. A mixture of ferric oxide and sulfate was kept in a closed system to ensure establishment of equilibrium partial pressure at the electrode. The cell was designed to avoid physical contact between Fe2(SO4)3 and Na2SO4 electrolyte. Uncertainties arising from the formation of sulfate solid solution were thus eliminated. The Gibbs’ energy of formation of ferric sulfate calculated from the emf is discussed in comparison with data reported in the literature. There is no evidence for the formation of oxysulfates in the Fe-S-0 system. Based on the results obtained in the present study for Fe2(SO4)3 and literature data for other phases, chemical potential diagrams have been constructed for the Fe-S-O system at 900 and 1100 K.

Oxygen content of liquid cobalt in equilibrium with CoO.(1+x)Al2O3 and -Al2O3

The oxygen concentration of liquid cobalt in equilibrium with cobalt aluminate and a-alumina has been measured by suction sampling and crucible quenching techniques at temperatures between 1770 and 1975 K. Experiments were made with cobalt of high and low initial oxygen contents, and with and without the addition of cobalt aluminate. The effect of temperature on the equilibrium oxygen content is represented by the equation, log (at.% 0) = -10,4001T(K) + 4.64 (±0.008). The composition of the spinel phase, CoO.(1+x)AI20 3, saturated with alumina, has been determined by electron probe microanalysis. The values of x are 0.22 at 1770 Kand 0.28 at 1975 K. The oxygen potential corresponding to the three-phase equilibrium between cobalt, aluminate and alumina, and the standard Gibbs' energy of formation of nonstoichiometric cobalt aluminate are evaluated by combining the results of this study with recently published data on the activity of oxygen in liquid cobalt. Implications of the present results to aluminium deoxidation of liquid cobalt are discussed.

Solubility and activity of oxygen in liquid nickel in equilibrium with -Al2O3 and NiO.(1+x)Al2O3

The limiting solubility of oxygen in liquid nickel in equilibrium withα-alumina and nickel aluminate has been measured by inert gas fusion analysis of suction samples in the temperature range 1730 to 1975 K. The corresponding oxygen potential has been monitored by a solid electrolyte cell consisting of calcia stabilized zirconia as the electrolyte and Mo + MoO2 as the reference electrode. The results can be summarized by the following equations: log(at. pct O) = \frac - 10,005T + 4.944 ( ±0.015)log(atpctO)=T−10005+4944(0015) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn DmO2 /4.606RT = log P O2 1/2 = \frac - 13,550T + 4.411 ( ±0.009)O24606RT=logPO212=T−13550+4411(0009) From simultaneous measurements of the potential and concentration of oxygen in melts, not in thermodynamic equilibrium with alumina and aluminate phases, information on the composition dependence of the activity coefficient and the standard free energy of solution of oxygen is obtained. For the reaction, $\frac{1}{2} O_2 \to \underset{\raise0.3em\hbox{$Missing close brace ΔG o = -72,930 - 7.11T (±840) J gr.at.–1 = + 0.216 at. pct OlogfO=T−500+0216atpctO where the standard state for dissolved oxygen is that which makes the value of activity equal to the concentration (in at. pct) in the limit as concentration approaches zero. The oxygen solubility in liquid nickel in equilibrium with solid NiO, evaluated from thermodynamic data, is compared with information reported in the literature. Implications of the results to the deoxidation equilibria of aluminum in nickel are discussed.

Measurement and modeling of Alloy-Spinel-Corundum equilibrium in the Ni-Mn-Al-0 system at 1873 K

The oxygen content of liquid Ni-Mn alloy equilibrated with spinel solid solution, (Ni,Mn)O. (1 +x)A12O3, and α-Al2O3 has been measured by suction sampling and inert gas fusion analysis. The corresponding oxygen potential of the three-phase system has been determined with a solid state cell incorporating (Y2O3)ThO2 as the solid electrolyte and Cr + Cr2O3 as the reference electrode. The equilibrium composition of the spinel phase formed at the interface of the alloy and alumina crucible was obtained using EPMA. The experimental data are compared with a thermodynamic model based on the free energies of formation of end-member spinels, free energy of solution of oxygen in liquid nickel, interaction parameters, and the activities in liquid Ni-Mn alloy and spinel solid solution. Mixing properties of the spinel solid solution are derived from a cation distribution model. The computational results agree with the experimental data on oxygen concentration, potential, and composition of the spinel phase.

Reactive interdiffusion: deviation from local equilibrium in the Fe-NiTiO3 couple

Deviation from local equilibrium between Fe–Ni alloy and (Fe,Ni)TiO3 solid solution in the reaction–diffusion zone of the Fe–NiTiO3 couple at 1273 K is evaluated by comparing the measured compositions in the zone with experimentally determined equilibrium tie-lines. The deviation is quantified by computing the Gibbs energy change for the reaction, Fe + NiTiO3 → FeTiO3 + Ni, from measured compositions in the zone and activity data available in the literature. Except near the extremities of the zone, the computed Gibbs energy change is constant, 8.2 kJ mol−1 higher than the standard Gibbs energy change for the reaction.