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.

Criticality of the exponential rate of decay for the largest nearest-neighbor link in random geometric graphs

Let n points be placed independently in d-dimensional space according to the density f(x) = A(d)e(-lambda parallel to x parallel to alpha), lambda, alpha > 0, x is an element of R-d, d >= 2. Let d(n) be the longest edge length of the nearest-neighbor graph on these points. We show that (lambda(-1) log n)(1-1/alpha) d(n) - b(n) converges weakly to the Gumbel distribution, where b(n) similar to ((d - 1)/lambda alpha) log log n. We also prove the following strong law for the normalized nearest-neighbor distance (d) over tilde (n) = (lambda(-1) log n)(1-1/alpha) d(n)/log log n: (d - 1)/alpha lambda <= lim inf(n ->infinity) (d) over tilde (n) <= lim sup(n ->infinity) (d) over tilde (n) <= d/alpha lambda almost surely. Thus, the exponential rate of decay alpha = 1 is critical, in the sense that, for alpha > 1, d(n) -> 0, whereas, for alpha <= 1, d(n) -> infinity almost surely as n -> infinity.

Effect of Strain Rate on Evolution of the Deformation Microstructure and Texture in Polycrystalline Copper and Nickel

The evolution of crystallographic texture in polycrystalline copper and nickel has been studied. The deformation texture evolution in these two materials over seven orders of magnitude of strain rate from 3 x 10(-4) to similar to 2.0 x 10(+3) s(-1) show little dependence on the stacking fault energy (SFE) and the amount of deformation. Higher strain rate deformation in nickel leads to weakerh < 101 > texture because of extensive microband formation and grain fragmentation. This behavior, in turn, causes less plastic spin and hence retards texture evolution. Copper maintains the stable end < 101 > component over large strain rates (from 3 x 10(-4) to 10(+2) s(-1)) because of its higher strain-hardening rate that resists formation of deformation heterogeneities. At higher strain rates of the order of 2 x 10(+3) s(-1), the adiabatic temperature rise assists in continuous dynamic recrystallization that leads to an increase in the volume fraction of the < 101 > component. Thus, strain-hardening behavior plays a significant role in the texture evolution of face-centered cubic materials. In addition, factors governing the onset of restoration mechanisms like purity and melting point govern texture evolution at high strain rates. SFE may play a secondary role by governing the propensity of cross slip that in turn helps in the activation of restoration processes.

High-Temperature (1023 K to 1273 K 750 degrees C to 1000 degrees C]) Plastic Deformation Behavior of B-Modified Ti-6Al-4V Alloys: Temperature and Strain Rate Effects

A minor addition of B to the Ti-6Al-4V alloy, by similar to 0.1 wt pct, reduces its as-cast prior beta grain size by an order of magnitude, whereas higher B content leads to the presence of in situ formed TiB needles in significant amounts. An experimental investigation into the role played by these microstructural modifications on the high-temperature deformation behavior of Ti-6Al-4V-xB alloys, with x varying between 0 wt pct and 0.55 wt pct, was conducted. Uniaxial compression tests were performed in the temperature range of 1023 K to 1273 K (750 degrees C to 1000 degrees C) and in the strain rate range of 10(-3) to 10(+1) s(-1). True stress-true strain responses of all alloys exhibit flow softening at lower strain rates and oscillations at higher strain rates. The flow softening is aided by the occurrence of dynamic recrystallization through lath globularization in high temperature (1173 K to 1273 K 900 degrees C to 1000 degrees C]) and a lower strain rate (10(-2) to 10(-3) s(-1)) regime. The grain size refinement with the B addition to Ti64, despite being marked, had no significant effect on this. Oscillations in the flow curve at a higher strain rate (10(0) to 10(+1) s(-1)), however, are associated with microstructural instabilities such as bending of laths, breaking of lath boundaries, generation of cavities, and breakage of TiB needles. The presence of TiB needles affected the instability regime. Microstructural evidence suggests that the matrix cavitation is aided by the easy fracture of TiB needles.

Effect of favorable pressure gradient on transitional spot formation rate

The calculation of the transitional boundary layer requires estimates of the extent of the transition zone, which in turn depends on the rate at which turbulent spots are formed. This rate has been found to scale with local boundary layer thickness and viscosity, and the resulting nondimensional group (called crumble) is a function of the pressure gradient, among other parameters. Available experimental data are analyzed to show that the crumble increases slowly with increasing favorable pressure gradients, being about four times as large as in constant-pressure flow when the Thwaites pressure gradient parameter at the effective origin of the resulting turbulent boundary layer is 0.1 and when transition is driven by free-stream turbulence.

Measurement of the corrosion rate of aluminium in sodium hydroxide using holographic interferometry

The application of holographic interferometry to the measurement of the corrosion rate of aluminium in sodium hydroxide is investigated. Details of the fabrication of the corrosion cell and the experimental procedure are given. Thickness loss of aluminium was found for different dissolution times and compared with the conventional weight-loss method using a microbalance.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Position and Feed-Rate Control for Contouring Numerical-Control Systems

Numerical control (NC) for contouring operations requires precise control of position and feed rate for approximating the contour by linear moves of the cutter. A control scheme, for generating linear moves with desired slopes for the cutter, is described. This scheme provides for nine successive linear moves, and may be either expanded or implemented in succession, for approximating a contour.

The velocity dispersion of giant molecular clouds. II - Mathematical and numerical refinements

The energy input to giant molecular clouds is recalculated, using the proper linearized equations of motion, including the Coriolis force and allowing for changes in the guiding center. Perturbation theory yields a result in the limit of distant encounters and small initial epicyclic amplitudes. Direct integration of the motion equations allows the strong encounter regime to be studied. The present perturbation theory result differs by a factor of order unity from that of Jog and Ostriker (1988). The result of present numerical integrations for the 2D (planar) velocity dispersion is presented. The accretion rate for a molecular cloud in the Galactic disk is calculated.

Hypersonic stagnation‐point boundary layers with massive blowing in the presence of a magnetic field

The effect of massive blowing rates on the steady laminar hypersonic boundary‐layer flow of an electrically conducting fluid in the stagnation region of an axisymmetric body with an applied magnetic field has been studied. The governing equations have been solved numerically by combining the implicit finite‐difference scheme with the quasi‐linearization technique. It is observed that the effect of massive blowing rates is to remove the viscous layer away from the boundary, whereas the effect of the magnetic field is just the opposite. It is also found that the velocity overshoot increases with blowing rates and also with magnetic field. The effect of the variation of the density‐viscosity product across the boundary layer is strong only when the blowing rate is small, but for the massive blowing rate the effect is negligible.

Nonnegative integral solution of linear equations

A method to obtain a nonnegative integral solution of a system of linear equations, if such a solution exists is given. The method writes linear equations as an integer programming problem and then solves the problem using a combination of artificial basis technique and a method of integer forms.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation