Convective flow during dendritic growth

A review is presented of the major experimental findings obtained from recent ground-based research conducted under the SPAR program. Measurements of dendritic growth at small supercoolings indicate that below approximately 1.5 K a transition occurs from diffusive control to convective control in succinonitrile, a model system chosen for this study. The key theoretical ideas concerning diffusive and convective heat transport during dendritic growth are discussed, and it is shown that a transition in the transport control should occur when the characteristic length for diffusion becomes larger than the characteristic length for convection. The experimental findings and the theoretical ideas discussed suggest that the Fluid Experiment System could provide appropriate experimental diagnostics for flow field visualization and quantification of the fluid dynamical effects presented here

Spontaneous deterministic side-branching behavior in phase-field simulations of equiaxed dendritic growth

The accepted view on dendritic side-branching is that side-branches grow as the result of selective amplification of thermal noise and that in the absence of such noise dendrites would grow without the development of side-arms. However, recently there has been renewed speculation about dendrites displaying deterministic side-branching [see, e.g., M. E. Glicksman, Metall. Mater. Trans A 43, 391 (2012)]. Generally, numerical models of dendritic growth, such as phase-field simulation, have tended to display behaviour which is commensurate with the former view, in that simulated dendrites do not develop side-branches unless noise is introduced into the simulation. However, here, we present simulations that show that under certain conditions deterministic side-branching may occur. We use a model formulated in the thin interface limit and a range of advanced numerical techniques to minimise the numerical noise introduced into the solution, including a multigrid solver. Spontaneous side-branching seems to be favoured by high undercoolings and by intermediate values of the capillary anisotropy, with the most branched structures being obtained for an anisotropy strength of 0.03. From an analysis of the tangential thermal gradients on the solid-liquid interface, the mechanism for side-branching appears to have some similarities with the deterministic model proposed by Glicksman

Anisotropic diffusion limited aggregation in three dimensions : universality and nonuniversality

We explore the macroscopic consequences of lattice anisotropy for diffusion limited aggregation (DLA) in three dimensions. Simple cubic and bcc lattice growths are shown to approach universal asymptotic states in a coherent fashion, and the approach is accelerated by the use of noise reduction. These states are strikingly anisotropic dendrites with a rich hierarchy of structure. For growth on an fcc lattice, our data suggest at least two stable fixed points of anisotropy, one matching the bcc case. Hexagonal growths, favoring six planar and two polar directions, appear to approach a line of asymptotic states with continuously tunable polar anisotropy. The more planar of these growths visually resembles real snowflake morphologies. Our simulations use a new and dimension-independent implementation of the DLA model. The algorithm maintains a hierarchy of sphere coverings of the growth, supporting efficient random walks onto the growth by spherical moves. Anisotropy was introduced by restricting growth to certain preferred directions

Evidance for an Oxygen Diffusion Model for the Electric Pulse Induced Resistance Change Effect in Oxides

Electric pulse induced resistance (EPIR) switching hysteresis loops for Pr0.7Ca0.7MnO3 (PCMO) perovskite oxide films were found to exhibit an additional sharp "shuttle peak" around the negative pulse maximum for films deposited in an oxygen deficient ambient. The device resistance hysteresis loop consists of stable high resistance and low resistance states, and transition regions between them. The resistance relaxation of the "shuttle peak" and its temperature behavior as well as the resistance relaxation in the transition regions were studied, and indicate that the resistance switching relates to oxygen diffusion with activation energy about 0.4eV. An oxygen diffusion model with the oxygen ions (vacancies) as the active agent is proposed for the non-volatile resistance switching effect in PCMO.Comment: 7 pages, 5 figure

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results

Finite-Difference Lattice Boltzmann Methods for binary fluids

We investigate two-fluid BGK kinetic methods for binary fluids. The developed theory works for asymmetric as well as symmetric systems. For symmetric systems it recovers Sirovich's theory and is summarized in models A and B. For asymmetric systems it contributes models C, D and E which are especially useful when the total masses and/or local temperatures of the two components are greatly different. The kinetic models are discretized based on an octagonal discrete velocity model. The discrete-velocity kinetic models and the continuous ones are required to describe the same hydrodynamic equations. The combination of a discrete-velocity kinetic model and an appropriate finite-difference scheme composes a finite-difference lattice Boltzmann method. The validity of the formulated methods is verified by investigating (i) uniform relaxation processes, (ii) isothermal Couette flow, and (iii) diffusion behavior.Comment: RevTex, 3 figures. Phys. Rev. E (2005, in press

Undercooling of acoustically levitated molten drops

It was observed that the acoustically levitated molten SCN (succinonitrile) drops can generally be undercooled to a degree where the impurities in the drop are responsible for the nucleation of the solid phase. However, it was also observed that ultrasound occasionally terminates undercooling of the levitated drops by initiating the nucleation of the solid at an undercooling level which is lower than that found for the nucleation catalyzed by the impurities in the drop. This premature nucleation can be explained by thermodynamic considerations which predict an increase in effective undercooling of the liquid upon the collapse of cavities. Pre-existing gas microbubbles which grow under the influence of ultrasound are suggested as the source of cavitation. The highly undercooled SCN drops can be utilized to measure the growth velocity of the solid in the deeply undercooled region including the hypercooled region

Local and Global Bifurcations of Flow Fields During Physical Vapor Transport: Application to a Microgravity Experiment

The local bifurcation of the flow field, during physical vapor transport for a parametric range of experimental interest, shows that its dynamical state ranges from steady to aperiodic. Comparison of computationally predicted velocity profiles with laser doppler velocimetry measurements shows reasonable agreement in both magnitude and planform. Correlation of experimentally measured crystal quality with the predicted dynamical state of the flow field shows a degradation of quality with an increase in Rayleigh number. The global bifurcation of the flow field corresponding to low crystal quality indicates the presence of a traveling wave for Ra = 1.09 x 10(exp 5). For this Rayleigh number threshold a chaotic transport state occurs. However, a microgravity environment for this case effectively stabilizes the flow to diffusive-advective and provides the setting to grow crystals with optimal quality

Advanced wet-dry cooling tower concept

Thesis. 1977. M.S.--Massachusetts Institute of Technology. Dept. of Mechanical Engineering.The purpose of this years' work has been to test and analyze the new dry cooling tower surface previously developed. The model heat transfer test apparatus built last year has been instrumented for temperature, humidity and flow measurement and performance has been measured under a variety of operating conditions. Tower Tests showed approximately 40-50% of the total energy transfer as taking place due to evaporation. This can be compared to approximately 80 to 85% for a conventional wet cooling tower. Comparison of the model tower test results with those of a computer simulation has demonstrated the validity of that simulation and its use as a design tool. Computer predictions have been made for a full-size tower system operating at several locations. Experience with this counterflow model tower has suggested that several design problems may be avoided by blowing the cooling air horizontally through the packing section. This crossflow concept was built from the previous counterflow apparatus and included the design and fabrication of new packing plates. Instrumentation and testing of the counterflow model produced data with an average experimental error of 10%. These results were compared to the predictions of a computer model written for the crossflow configuration. In 14 test runs the predicted total heat transfer differed from the measured total heat transfer by no more than 8% with most runs coming well within 5%. With the computer analogy's validity established, it may now be used to help predict the performance of fullscale wet-dry towers