An explicit scheme for coupling temperature and concentration fields in solidification models

A numerical scheme for coupling temperature and concentration fields in a general solidification model is presented. A key feature of this scheme is an explicit time stepping used in solving the governing thermal and solute conservation equations. This explicit approach results in a local point-by-point coupling scheme for the temperature and concentration and avoids the multi-level iteration required by implicit time stepping schemes. The proposed scheme is validated by predicting the concentration field in a benchmark solidification problem. Results compare well with an available similarity solution. The simplicity of the proposed explicit scheme allows for the incorporation of complex microscale models into a general solidification model. This is demonstrated by investigating the role of dendrite coarsening on the concentration field in the solidification benchmark problem.

The multiphysics modeling of solidification and melting processes

Solidification and melting processes involve a range of physical phenomena and their interactions (i.e., multiphysics). Computational modeling of such processes presents a significant challenge, both in representing the physics involved and in handling the resulting coupled behavior. Two methods for the computational modeling of multiphysics processes in complex geometries are highlighted in the context of four challenging applications

Assessing the scalability of multiphysics tools for modeling solidification and melting processes on parallel clusters

A comprehensive solution of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, nonlinear solid mechanics and, possibly, electromagnetics together with their interactions, in what is now known as multiphysics simulation. Such simulations are computationally intensive and the implementation of solution strategies for multiphysics calculations must embed their effective parallelization. For some years, together with our collaborators, we have been involved in the development of numerical software tools for multiphysics modeling on parallel cluster systems. This research has involved a combination of algorithmic procedures, parallel strategies and tools, plus the design of a computational modeling software environment and its deployment in a range of real world applications. One output from this research is the three-dimensional parallel multiphysics code, PHYSICA. In this paper we report on an assessment of its parallel scalability on a range of increasingly complex models drawn from actual industrial problems, on three contemporary parallel cluster systems.

The effects of thermo-electrically induced convection in alloy solidification

The effects of a constant uniform magnetic field on thermoelectric currents during dendritic solidification were investigated using a 2-dimensional enthalpy based numerical model. Using an approximation of the dendrite growing in free space it was found that the resulting Lorentz force generates a circulating flow influencing the solidification pattern. As the magnetic field strength increases it was found that secondary growth on the clockwise side of the primary arm of the dendrite was encouraged, while the anticlockwise side is suppressed due to a reduction in local free energy. The preferred direction of growth rotated in the clockwise sense under an anti-clockwise flow for both the binary alloy and pure material. The tip velociy is significantly increased compared to growth in stagnant flow. This is due to a small recirculation that follows the tip of the dendrite; bringing in colder liquid and lower concentrations of solute. The recirculation being not normally incident on the tip is most likely the cause for the rotation. Grain growth consisting of multiple seeds with the same anisotropy growing in the same plane, gives a competition to release latent heat resulting in stunted growth. The initial growth for each dendrite is very similar to the single seed cases indicating that dendrites must become before the thermoelectric interactions are significant.

Three-dimensional simulation of the influence of convection on dendritic solidification

Bulk and interdendritic flow during solidification alters the microstructure development, potentially leading to the formation of defects. In this paper, a 3D numerical model is presented for the simulation of dendritic growth in the presence of fluid flow in both liquid and semi-solid zones during solidification. The dendritic growth was solved by the combination of a stochastic nucleation approach with a finite difference solution of the solute diffusion equation and. a projection method solution of the Navier-Stokes equations. The technique was applied first to simulate the growth of a single dendrite in 2D and 3D in an isothermal environment with forced fluid flow. Significant differences were found in the evolution of dendritic morphology when comparing the 2D and 3D results. In 3D the upstream arm has a faster growth velocity due to easier flow around the perpendicular arms. This also promotes secondary arm formation on the upstream arm. The effect of fluid flow on columnar dendritic growth and micro-segregation in constrained solidification conditions is then simulated. For constrained growth, 2D simulations lead to even greater inaccuracies as compared to 3D.

Thermoelectric effects on alloy solidification microstructure

Thermoelectric currents in the presence of a magnetic field generate Lorentz forces which can drive fluid flow. In the case of dendritic growth a naturally occurring thermoelectric current exists and in the presence of a high magnetic field micro convections are generated. Experimental evidence has attributed changes in microstructure to this effect. A numerical model has been developed to study the flow field around an unconstricted equiaxed dendrite growing under these conditions. The growth is modeled in 2D and 3D by an enthalpy based method and a complex flow structure has been predicted. Using a pseudo-3D approximation for economy, realistic 2D simulations are obtained where a fully coupled transient scheme reveals significant changes to the dendrite morphology reflecting experimental evidence. There is a rotation of the preferred direction of growth and increased secondary branching.

Application of accelerated carbonation with a combination of Na2CO3 and CO2 in cement-based solidification/stabilization of heavy metal-bearing sediment

The efficient remediation of heavy metal-bearing sediment has been one of top priorities of ecosystem protection. Cement-based solidification/stabilization (s/s) is an option for reducing the mobility of heavy metals in the sediment and the subsequent hazard for human beings and animals. This work uses sodium carbonate as an internal carbon source of accelerated carbonation and gaseous CO2 as an external carbon source to overcome deleterious effects of heavy metals on strength development and improve the effectiveness of s/s of heavy metal-bearing sediment. In addition to the compressive strength and porosity measurements, leaching tests followed the Chinese solid waste extraction procedure for leaching toxicity - sulfuric acid and nitric acid method (HJ/T299-2007), German leaching procedure (DIN38414-S4) and US toxicity characteristic leaching procedures (TCLP) have been conducted. The experimental results indicated that the solidified sediment by accelerated carbonation was capable of reaching all performance criteria for the disposal at a Portland cement dosage of 10 wt.% and a solid/water ratio of 1: 1. The concentrations of mercury and other heavy metals in the leachates were below 0.10 mg/L and 5 mg/L, respectively, complying with Chinese regulatory level (GB5085-2007). Compared to the hydration, accelerated carbonation improved the compressive strength of the solidified sediment by more than 100% and reduced leaching concentrations of heavy metals significantly. It is considered that accelerated carbonation technology with a combination of Na2CO3 and CO2 may practically apply to cement-based s/s of heavy metal-bearing sediment. (C) 2008 Elsevier B.V. All rights reserved.

Immobilisation of heavy metal in cement-based solidification/stabilisation: a review

Heavy metal-bearing waste usually needs solidification/stabilization (s/s) prior to landfill to lower the leaching rate. Cement is the most adaptable binder currently available for the immobilisation of heavy metals. The selection of cements and operating parameters depends upon an understanding of chemistry of the system. This paper discusses interactions of heavy metals and cement phases in the solidification/stabilisation process. It provides a clarification of heavy metal effects on cement hydration. According to the decomposition rate of minerals, heavy metals accelerate the hydration of tricalcium silicate (C3S) and Portland cement, although they retard the precipitation of portlandite due to the reduction of pH resulted from hydrolyses of heavy metal ions. The chemical mechanism relevant to the accelerating effect of heavy metals is considered to be H+ attacks on cement phases and the precipitation of calcium heavy metal double hydroxides, which consumes calcium ions and then promotes the decomposition Of C3S. In this work, molecular models of calcium silicate hydrate gel are presented based on the examination of Si-29 solid-state magic angle spinning/nuclear magnetic resonance (MAS/NMR). This paper also reviews immobilisation mechanisms of heavy metals in hydrated cement matrices, focusing on the sorption, precipitation and chemical incorporation of cement hydration products. It is concluded that further research oil the phase development during cement hydration in the presence of heavy metals and thermodynamic modelling is needed to improve effectiveness of cement-based s/s and extend this waste management technique. (C) 2008 Elsevier Ltd. All rights reserved.

Thermoelectric MHD in dendritic solidification

The effects of a constant uniform magnetic field on dendritic solidification were investigated using an enthalpy based numerical model. The interaction between thermoelectric currents on a growing crystal and the magnetic field generates a Lorentz force that creates flow. The need for very high resolution at the liquid-solid boundary where the thermoelectric source originates plus the need to accommodate multiple grains for a realistic simulation, make this a very demanding computational problem. For practical simulations, a quasi 3-dimensional approximation is proposed which nevertheless retains essential elements of transport in the third dimension. A magnetic field normal to the plane of growth leads to general flow circulation around an equiaxed dendrite, with secondary recirculations between the arms. The heat/solute advection by the flow is shown to cause a change in the morphology of the dendrite; secondary growth is promoted preferentially on one side of the dendrite arm and the tip velocity of the primary arm is increased. The degree of approximation introduced is quantified by extending the model into 3-dimensions, where the full Navier-Stokes equation is solved, and compared against the 2-dimensional solution.

Directional performance in motion transparency

Motion transparency provides a challenging test case for our understanding of how visual motion, and other attributes, are computed and represented in the brain. However, previous studies of visual transparency have used subjective criteria which do not confirm the existence of independent representations of the superimposed motions. We have developed measures of performance in motion transparency that require observers to extract information about two motions jointly, and therefore test the information that is simultaneously represented for each motion. Observers judged whether two motions were at 90 to one another; the base direction was randomized so that neither motion taken alone was informative. The precision of performance was determined by the standard deviations (S.D.s) of probit functions fitted to the data. Observers also made judgments of orthogonal directions between a single motion stream and a line, for one of two transparent motions against a line and for two spatially segregated motions. The data show that direction judgments with transparency can be made with comparable accuracy to segregated (non-transparent) conditions, supporting the idea that transparency involves the equivalent representation of two global motions in the same region. The precision of this joint direction judgment is, however, 2–3 times poorer than that for a single motion stream. The precision in directional judgment for a single stream is reduced only by a factor of about 1.5 by superimposing a second stream. The major effect in performance, therefore, appears to be associated with the need to compute and compare two global representations of motion, rather than with interference between the dot streams per se. Experiment 2tested the transparency of motions separated by a range of angles from 5 to 180 by requiring subjects to set a line matching the perceived direction of each motion. The S.D.s of these settings demonstrated that directions of transparent motions were represented independently for separations over 20. Increasing dot speeds from 1 to 10 deg/s improved directional performance but had no effect on transparency perception. Transparency was also unaffected by variations of density between 0.1 and 19 dots/deg2

Interaction of directional, neuromuscular and egocentric constraints on the stability of preferred bimanual coordination patterns

We investigated how the relative direction of limb movements in external space (iso- and non-isodirectionality), muscular constraints (the relative timing of homologous muscle activation) and the egocentric frame of reference (moving simultaneously toward/away the longitudinal axis of the body) contribute to the stability of coordinated movements. In the first experiment, we attempted to determine the respective stability of isodirectional and non-isodirectional movements in between-persons coordination. In a second experiment, we determined the effect of the relative direction in external space, and of muscular constraints, on pattern stability during a within-person bimanual coordination task. In the third experiment we dissociated the effects on pattern stability of the muscular constraints, relative direction and egocentric frame of reference. The results showed that (1) simultaneous activation of homologous muscles resulted in more stable performance than simultaneous activation of non-homologous muscles during within-subject coordination, and that (2) isodirectional movements were more stable than non-isodirectional movements during between-persons coordination, confirming the role of the relative direction of the moving limbs in the stability of bimanual coordination. Moreover, the egocentric constraint was to some extent found distinguishable from the effect of the relative direction of the moving limbs in external space, and from the effect of the relative timing of muscle activation. In summary, the present study showed that relative direction of the moving limbs in external space and muscular constraints may interact either to stabilize or destabilize coordination patterns. (C) 2003 Published by Elsevier B.V.

The hierarchy of directional interactions in visual motion processing

It is well known that context influences our perception of visual motion direction. For example, spatial and temporal context manipulations can be used to induce two well-known motion illusions: direction repulsion and the direction after-effect (DAE). Both result in inaccurate perception of direction when a moving pattern is either superimposed on (direction repulsion), or presented following adaptation to (DAE), another pattern moving in a different direction. Remarkable similarities in tuning characteristics suggest that common processes underlie the two illusions. What is not clear, however, is whether the processes driving the two illusions are expressions of the same or different neural substrates. Here we report two experiments demonstrating that direction repulsion and the DAE are, in fact, expressions of different neural substrates. Our strategy was to use each of the illusions to create a distorted perceptual representation upon which the mechanisms generating the other illusion could potentially operate. We found that the processes mediating direction repulsion did indeed access the distorted perceptual representation induced by the DAE. Conversely, the DAE was unaffected by direction repulsion. Thus parallels in perceptual phenomenology do not necessarily imply common neural substrates. Our results also demonstrate that the neural processes driving the DAE occur at an earlier stage of motion processing than those underlying direction repulsion.