Three-dimensional computational modeling of momentum, heat, and mass transfer in a laser surface alloying process

A three- dimensional, transient model is developed for studying heat transfer, fluid flow, and mass transfer for the case of a single- pass laser surface alloying process. The coupled momentum, energy, and species conservation equations are solved using a finite volume procedure. Phase change processes are modeled using a fixed-grid enthalpy-porosity technique, which is capable of predicting the continuously evolving solid- liquid interface. The three- dimensional model is able to predict the species concentration distribution inside the molten pool during alloying, as well as in the entire cross section of the solidified alloy. The model is simulated for different values of various significant processing parameters such as laser power, scanning speed, and powder feedrate in order to assess their influences on geometry and dynamics of the pool, cooling rates, as well as species concentration distribution inside the substrate. Effects of incorporating property variations in the numerical model are also discussed.

Grain Floatation During Equiaxed Solidification of an Al-Cu Alloy in a Side-Cooled Cavity: Part II-Numerical Studies

In this article, a single-phase, one-domain macroscopic model is developed for studying binary alloy solidification with moving equiaxed solid phase, along with the associated transport phenomena. In this model, issues such as thermosolutal convection, motion of solid phase relative to liquid and viscosity variations of the solid-liquid mixture with solid fraction in the mobile zone are taken into account. Using the model, the associated transport phenomena during solidification of Al-Cu alloys in a rectangular cavity are predicted. The results for temperature variation, segregation patterns, and eutectic fraction distribution are compared with data from in-house experiments. The model predictions compare well with the experimental results. To highlight the influence of solid phase movement on convection and final macrosegregation, the results of the current model are also compared with those obtained from the conventional solidification model with stationary solid phase. By including the independent movement of the solid phase into the fluid transport model, better predictions of macrosegregation, microstructure, and even shrinkage locations were obtained. Mechanical property prediction models based on microstructure will benefit from the improved accuracy of this model.

In-situ measurements of concentration and temperature during transient solidification of aqueous solution of ammonium chloride using laser interferometry

The variation in temperature and concentration plays a crucial role in predicting the final microstructure during solidification of a binary alloy. Most of the experimental techniques used to measure concentration and temperature are intrusive in nature and affect the flow field. In this paper, the main focus is laid on in-situ, non-intrusive, transient measurement of concentration and temperature during the solidification of a binary mixture of aqueous ammonium chloride solution (a metal-analog system) in a top cooled cavity using laser based Mach-Zehnder Interferometric technique. It was found from the interferogram, that the angular deviation of fringe pattern and the total number of fringes exhibit significant sensitivity to refractive index and hence are functions of the local temperature and concentration of the NH4Cl solution inside the cavity. Using the fringe characteristics, calibration curves were established for the range of temperature and concentration levels expected during the solidification process. In the actual solidification experiment, two hypoeutectic solutions (5% and 15% NH4Cl) were chosen. The calibration curves were used to determine the temperature and concentration of the solution inside the cavity during solidification of 5% and 15% NH4Cl solution at different instants of time. The measurement was carried out at a fixed point in the cavity, and the concentration variation with time was recorded as the solid-liquid interface approached the measurement point. The measurement exhibited distinct zones of concentration distribution caused by solute rejection and Rayleigh Benard convection. Further studies involving flow visualization with laser scattering confirmed the Rayleigh Benard convection. Computational modeling was also performed, which corroborated the experimental findings. (C) 2011 Elsevier Ltd. All rights reserved.

A coupled VOF-IBM-enthalpy approach for modeling motion and growth of equiaxed dendrites in a solidifying melt

A coupled methodology for simulating the simultaneous growth and motion of equiaxed dendrites in solidifying melts is presented. The model uses the volume-averaging principles and combines the features of the enthalpy method for modeling growth, immersed boundary method for handling the rigid solid-liquid interfaces, and the volume of fluid method for tracking the advection of the dendrite. The algorithm also performs explicit-implicit coupling between the techniques used. A two-dimensional framework with incompressible and Newtonian fluid is considered. Validation with available literature is performed and dendrite growth in the presence of rotational and buoyancy driven flow fields is studied. It is seen that the flow fields significantly alter the position and morphology of the dendrites. (C) 2012 Elsevier Inc. All rights reserved.

A model for growth and engulfment of gas microporosity during aluminum alloy solidification process

A new coupled approach is presented for modeling the hydrogen bubble evolution and engulfment during an aluminum alloy solidification process in a micro-scale domain. An explicit enthalpy scheme is used to model the solidification process which is coupled with a level-set method for tracking the hydrogen bubble evolution. The volume averaging techniques are used to model mass, momentum, energy and species conservation equations in the chosen micro-scale domain. The interaction between the solid, liquid and gas interfaces in the system have been studied. Using an order-of-magnitude study on growth rates of bubble and solid interfaces, a criterion is developed to predict bubble elongation which can occur during the engulfment phase. Using this model, we provide further evidence in support of a conceptual thought experiment reported in literature, with regard to estimation of final pore shape as a function of typical casting cooling rates. The results from the proposed model are qualitatively compared with in situ experimental observations reported in literature. The ability of the model to predict growth and movement of a hydrogen bubble and its subsequent engulfment by a solidifying front has been demonstrated for varying average cooling rates encountered in typical sand, permanent mold, and various casting processes. (C) 2012 Elsevier B.V. All rights reserved.

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Pattern-Formation During Self-organization in Three-Phase Eutectic Solidification

Eutectic growth is an interesting example for exploring the topic of pattern-formation in multi-phase systems, where the growth of the phases is coupled with the diffusive transport of one or more components in the melt. While in the case of binary alloys, the number of possibilities are limited (lamellae, rods, labyrinth etc.), their number rapidly increases with the number of components and phases. In this paper, we will investigate pattern formation during three-phase eutectic solidification using a state-of-the art phase-field method based on the grand-canonical density formulation. The major aim of the study is to highlight the role of two properties, which are the volume fraction of the solid phases and the solid-liquid interfacial energies, in the self-organization of the solid phases during directional growth. Thereafter, we will show representative phase-field simulations of a micro-structure in a real alloy (Ag-Al-Cu) using an asymmetric phase diagram as well as interfacial properties.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

Orientation Selection and Microstructural Evolution in Directionally Solidified Tb0.3Dy0.7Fe1.95

Tb0.3Dy0.7Fe1.95 alloy was directionally solidified by using a modified Bridgman technique at a wide range of growth rates of 5 to 100 cm/h. The directionally grown samples exhibited plane front solidification morphology up to a growth rate of 90 cm/h. Typical island banding feature was observed closer to the chilled end, which eventually gave rise to irregular peritectic coupled growth (PCG). The PCG gained prominence with an increase in the growth rate. The texture study revealed formation of strong aOE (c) 311 > texture in a lower growth rate regime, aOE (c) 110 > and ``rotated aOE (c) 110 > aEuroe in an intermediate growth regime, and aOE (c) 112 > in a higher growth rate regime. In-depth analysis of the atomic configuration of a solid-liquid interface revealed that the growth texture is influenced by the kinetics of atomic attachment to the solid-liquid interface, which is intimately related to a planar packing fraction and an atomic stacking sequence of the interfacial plane. The mechanism proposed in this article is novel and will be useful in addressing the orientation selection mechanism of topologically closed packed intermetallic systems. The samples grown at a higher growth rate exhibit larger magnetostriction (lambda) and d lambda/dH owing to the absence of pro-peritectic (Tb,Dy)Fe-3 and formation of aOE (c) 112 > texture, which lies closer to the easy magnetization direction (EMD).

An electrokinetic sensor for studying immersed surfaces, using focused ultrasound

A novel method for obtaining information on the charge density of an immersed surface is presented. The technique uses focused ultrasound to excite oscillatory fluid motion in the plane of the solid-liquid interface, over a localised area. The displacement current (resulting from the motion of fluid-borne ions in the outer double-layer) is detected by electrodes in the liquid. The method is demonstrated as a means for monitoring protein adsorption, and for monitoring interactions between two different proteins. A second electrokinetic effect at the interface is identified, isolated from the first, and shown to provide additional information on the compressibility and charge density of the double-layer. © 2001 Elsevier Science B.V. All rights reserved.

Influence of Forced Convection on Columnar Microstructure During Directional Solidification of Al-Ni Alloys

This paper presents a summary of cellular and dendritic morphologies resulting from the upward directional solidification of Al - Ni alloys in a cylindrical crucible. We analysed the coupling of solid-liquid interface morphology with natural and forced convection. The influence of natural convection was first analyzed as a function of growth parameters (solute concentration, growth rate and thermal gradient). In a second step, the influence of axial vibrations on solidification microstructure was investigated by varying vibration parameters (amplitude and frequency). Experimental results were compared to preliminary numerical simulations and a good agreement is found for natural convection. In this study, the critical role of the mushy zone in the interaction between fluid flow and solidification microstructure is pointed out.

In situ formation by laser cladding of a TiC composite coating with a gradient distribution

An in situ method was developed to produce an Ni alloy composite coating reinforced by in situ reacted TiC particles with a gradient distribution, using one-step laser cladding with a pre-placed powder mixture on a 5CrMnMo steel substrate. Dispersed and ultra-fine TIC particles were formed in situ in the coating. Most. of the TiC particles, with a marked gradient distribution, were uniformly distributed within interdendritic regions because of the trapping effect of the advancing solid-liquid interface. In addition, the TiC-gamma-Ni interfaces generated in situ were found to be free from any deleterious surface reaction. Finally, the microhardness also showed a gradient variation, with the highest value of 1250 Hv0.2 and the wear properties of the coating were significantly enhanced.

Characteristics Of Energy Dissipation In Hyperconcentrated Flows

An equilibrium equation for the turbulence energy in sediment-laden flows was derived on the basis of solid-liquid two-phase flow theory. The equation was simplified for two-dimensional, uniform, steady and fully developed turbulent hyperconcentrated flows. An energy efficiency coefficient of suspended-load motion was obtained from the turbulence energy equation, which is defined as the ratio of the sediment suspension energy to the turbulence energy of the sediment-laden flows. Laboratory experiments were conducted to investigate the characteristics of energy dissipation in hyperconcentrated flows. A total of 115 experimental runs were carried out, comprising 70 runs with natural sediments and 45 runs with cinder powder. Effects of sediment concentration on sediment suspension energy and flow resistance were analyzed and the relation between the energy efficiency coefficient of suspended-load motion and sediment concentration was established on the basis of experimental data. Furthermore, the characteristics of energy dissipation in hyperconcentrated flows were identified and described. It was found that the high sediment concentration does not increase the energy dissipation; on the contrary, it decreases flow resistance.