Phase-field simulation of dendritic growth in a binary alloy with thermodynamics data

This paper simulates the dendrite growth process during non-isothermal solidification in the Al-Cu binary alloy by using the phase-field model. The heat transfer equation is solved simultaneously. The thermodynamic and kinetic parameters are directly obtained from existing database by using the Calculation of Phase Diagram （CALPHAD） method. The effects of the latent heat and undercooling on the dendrite growth, solute and temperature profile during the solidification of binary alloy are investigated. The results indicate that the dendrite growing morphologies could be simulated realistically by linking the phase-field method to CALPHAD. The secondary arms of solidification dendritic are better developed with the increase of undercooling. Correspondingly, the tip speed and the solute segregation in solid-liquid interface increase, but the tip radius decreases.

Phase-field modeling of dendritic growth of magnesium alloys with a parallel-adaptive mesh refinement algorithm

A two-dimensional phase field(PF)model was developed to simulate the dendritic solidification in magnesium alloy with hcp crystal structure.By applying a parallel-adaptive mesh refinement(Para-AMR)algorithm,the computational efficiency of the numerical model was greatly improved.Based on the PF model,a series of simulation cases were conducted and the results showed that the anisotropy coefficient and coupling coefficient had a great influence on the dendritic morphology of magnesium alloy.The dendritic growth kinetics was determined by the undercooling and equilibrium solute partition coefficient.A significant finding is acquired that with a large undercooling,the maximum solute concentration is located on both sides of the dendrite tip in the liquid,whereas the maximum solute concentration gradient is located right ahead of the dendrite tip in the liquid.The dendrite tip growth velocity decreases with the increase of the equilibrium solute partition coefficient,while the variation trend of the dendrite tip radius is the opposite.Quantitative analysis was carried out relating to the dendritic morphology and growth kinetics,and the simulated results are consistent with the theoretical models proposed in the previously published works.

Numerical simulation on dendritic growth of Al-Cu alloy under convection based on the cellular automaton lattice Boltzmann method

A numerical model is developed by coupling the cellular automaton(CA)method and the lattice Boltzmann method(LBM)to simulate the dendritic growth of Al-Cu alloy in both two and three dimensions.An improved decentered square algorithm is proposed to overcome the artificial anisotropy induced by the CA cells and to realize simulation of dendritic growth with arbitrary orientations.Based on the established CA-LBM model,effects of forced convection and gravity-driven natural convection on dendritic growth are studied.The simulation results show that the blocking effect of dendrites on melt flow is advanced with a larger number of seeds.The competitive growth of the converging columnar dendrites is determined by the interaction between heat flow and forced convection.Gravity-driven natural convection leads to highly asymmetric growth of equiaxed dendrites.With sinking downwards of the heavy solute,chimney-like or mushroom-like solute plumes are formed in the melt in front of the columnar dendrites when they grow along the gravitational direction.More details on dendritic growth of Al-Cu alloy under convection are revealed by 3D simulations.

Simulation of Dendritic Growth with Melt Convection in Solidification of Ternary Alloys

A cellular automaton-lattice Boltzmann coupled model is extended to study the dendritic growth with melt convection in the solidification of ternary alloys. With a CALPHAD-based phase equilibrium engine, the effects of melt convection, solutal diffusion, interface curvature and preferred growth orientation are incorporated into the coupled model. After model validation, the multi dendritic growth of the Al-4.0 wt%Cu-1.0 wt%Mg alloy is simulated under the conditions of pure diffusion and melt convection. The result shows that the dendritic growth behavior, the final microstructure and microsegregation are significantly influenced by melt convection in the solidification.

Three-dimensional interfacial wave theory of dendritic growth: (Ⅰ). multiple variables expansion solutions

Dendritic pattern formation at the interface between liquid and solid is a commonly observed phenomenon in crystal growth and solidification process. The theoretical investigation of dendritic growth is one of the most profound and highly challenging subjects in the broad areas of interfacial pattern formation, condensed matter physics and materials science, preoccupying many researchers from various areas. Some longstanding key issues on this subject finally gained a breakthrough in the late of last century, via the ＇Interracial Wave （IFW） Theory＇ on the ground of systematical global stability analysis of the basic state of dendritic growth. The original form of the IFW theory mainly focus on the investigation of various axi-symmetric unsteady perturbed modes solutions around the axi-symmetrie basic state of system of dendritic growth. In reality, the system may allow various non-axi-symmetric, unsteady perturbed states. Whether or not the system of dendritic growth allows some growing non-axi-symmetric modes？ Will the stationary dendritic pattern be destroyed by some of such non-axi- symmetric modes？ Or, in one word, what is the stability property of the system, once the non-axi-symmetric modes can be evoked？ The answers for these questions are important for the solid foundation of IFW theory. The present work attempts to settle down these issues and develop a three-dimensional （3D） interfacial wave theory of dendritic growth. Our investigations verify that dendritic growth indeed allows a discrete set of non-axi-symmetric unstable global wave modes, which gives rise to a set of multiple arms spiral waves propagating along the Ivantsov＇s paraboloid.

Three-dimensional interfacial wave theory of dendritic growth: (Ⅱ). non-axi-symmetric global wave modes and selection criterion of pattern formation

This paper is ttie continuation of part （Ⅰ）, which completes the derivations of the 3D global wave modes solutions, yields the stability criterion and, on the basis of the results obtained, demonstrates the selection criterion of pattern formation.

Simulation of dendritic growth of magnesium al oys with fluid flow

Fluid flow has a significant impact on the microstructure evolution of alloys during solidification. Based on the previous work relating simulation of the dendritic growth of magnesium alloys with hcp(hexagonal closepacked) structure, an extension was made to the formerly established CA(cellular automaton) model with the purpose of studying the effect of fluid flow on the dendritic growth of magnesium alloys. The modified projection method was used to solve the transport equations of flow field. By coupling the flow field with the solute field, simulation results of equiaxed and columnar dendritic growth of magnesium alloys with fluid flow were achieved. The simulated results were quantitatively compared with those without fluid flow. Moreover, a comparison was also made between the present work and previous works conducted by others. It can be concluded that a deep understanding of the dendritic growth of magnesium alloys with fluid flow can be obtained by applying the present numerical model.

Effect of buoyancy-driven convection on steady state dendritic growth in a binary alloy

The effect of buoyancy-driven convection on the steady state dendritic growth in an undercooled binary alloy is studied. For the case of the moderate modified Grashof number, the uniformly valid asymptotic solution in the entire region of space is obtained by means of the matched asymptotic expansion method. The analytical results show that the buoyancy- driven convection has a significant effect on the needle-like interface of dendritic growth. Due to the buoyancy-driven convection, the needle-like interface shape of the crystal is changed. When the Peclet number that is not affected by the buoyant flow is less than a certain critical value, the interface shape of the dendrite becomes thinner as the Grashof number increases; when it is larger than the critical value, the interface shape becomes fatter as the Grashof number increases. In the undercooled binary alloy the morphology number plays an active role in the interface shape and leads to the buoyancy effect that is different from the situation for the pure melt. The smaller the morphology number is, the more significant change the interface shape has. As the Peclet number further increases, the effect of buoyancy on the interface diminishes eventually.

A Diffuse-Interface Lattice Boltzmann Method for the Dendritic Growth with Thermosolutal Convection

In this work,we proposed a diffuse-interface model for the dendritic growth with thermosolutal convection.In this model,the sharp boundary between the fluid and solid dendrite is firstly replaced by a thin but nonzero thickness diffuse interface,which is described by the order parameter,and the diffuse-interface based governing equations for the dendritic growth are presented.To solve the model for the dendritic growth with thermosolutal convection,we also developed a diffuse-interface multirelaxation-time lattice Boltzmann(LB)method.In this method,the order parameter in the phase-field equation is combined into the force caused by the fluid-solid interaction,and the treatment on the complex fluid-solid interface can be avoided.In addition,four LB models are designed for the phase-field equation,concentration equation,temperature equation and the Navier-Stokes equations in a unified framework.Finally,we performed some simulations of the dendritic growth to test the present diffuse-interface LB method,and found that the numerical results are in good agreements with some previous works.

Growth and inhibition of zinc anode dendrites in Zn-air batteries:Model and experiment

Zinc(Zn)-air batteries are widely used in secondary battery research owing to their high theoretical energy density,good electrochemical reversibility,stable discharge performance,and low cost of the anode active material Zn.However,the Zn anode also leads to many challenges,including dendrite growth,deformation,and hydrogen precipitation self-corrosion.In this context,Zn dendrite growth has a greater impact on the cycle lives.In this dissertation,a dendrite growth model for a Zn-air battery was established based on electrochemical phase field theory,and the effects of the charging time,anisotropy strength,and electrolyte temperature on the morphology and growth height of Zn dendrites were studied.A series of experiments was designed with different gradient influencing factors in subsequent experiments to verify the theoretical simulations,including elevated electrolyte temperatures,flowing electrolytes,and pulsed charging.The simulation results show that the growth of Zn dendrites is controlled mainly by diffusion and mass transfer processes,whereas the electrolyte temperature,flow rate,and interfacial energy anisotropy intensity are the main factors.The experimental results show that an optimal electrolyte temperature of 343.15 K,an optimal electrolyte flow rate of 40 ml·min^(-1),and an effective pulse charging mode.

Phase-field lattice-Boltzmann study on fully coupled thermal-solute-convection dendrite growth of Al-Cu alloy

Dendrite growth is a complex liquid-solid phase transition process involving multiple physical factors.A phase-field lattice-Boltzmann method was developed to simulate the two-and three-dimension dendrite growth of Al-Cu alloy.The effect of fully coupled thermal-solute-convection interaction on the dendrite growth was investigated by incorporating a parallel-adaptive mesh refinement algorithm into the numerical model.By accurately reproducing the latent heat release,solute diffusion and convective transport behaviors at the liquidsolid interface,the interaction mechanism among thermal-solute-convection transport as well as their coupling effects on the dendrite growth dynamics were discussed.The simulation results show that the release of latent heat slows down the dendrite growth rate,and both natural and forced convection disrupt the symmetrical growth of dendrites.Their combination makes the growth of dendrites more complex,capturing important physical aspects such as recalescence,dendrite tip splitting,dendrite tilting,dendrite remelting,and solute plume in the simulation case.Based on the robustness and powerful ability of the numerical model,the formation mechanisms of these physical aspects were revealed.

A Modified Cellular Automaton Model for the Quantitative Prediction of Equiaxed and Columnar Dendritic Growth

Since the characteristic of dendrite is an important factor determining the performance of castings, a twodimensional cellular automaton model with decentered square algorithm is developed for quantitatively predicting the dendritic growth during solidification process. The growth kinetics of solid/liquid interface are determined by the local equilibrium composition and local actual liquid composition, and the calculation of the solid fraction increment is based on these two compositions to avoid the solution of growth velocity. In order to validate the developed model, quantitative simulations of steady-state dendritic features over a range of undercooling was performed and the results exhibited good agreement with the predictions of LGK（Liptone Glicksman-Kurz） model. Meanwhile, it is demonstrated that the proposed model can be applied to simulate multiple equiaxed dendritic growth, as well as columnar dendritic growth with or without equiaxed grain formation in directional solidification of AleC u alloys. It has been shown that the model is able to simulate the growth process of multi-dendrites with various preferential orientations and can reproduce a wide range of complex dendritic growth phenomena such as nucleation, coarsening of dendrite arms, side branching in dendritic morphologies, competitive growth as well as the interaction among surrounding dendrites.

Study on Dendritic Growth in Pressurized Solidification of Mg–Al Alloy Using Phase Field Simulation

A multi-grain phase field model coupled with thermodynamic calculation was adopted to describe the dendritic growth in pressurized solidification of Mg-A1 alloy during squeeze casting, in which the effects of the pressure on the Gibbs free energy and chemical potential of solid and liquid phases, the solute diffusion coefficient, and the solute partition coefficient were considered. The individual effect of solute diffusion coefficient, and the Gibbs free energy on the dendritic growth was studied. With the compar- ison of the dendritic growth under atmospheric and elevated pressures, the effect of pressure on the microstructure evolution was discussed. The results showed that the grains are refined, the dendritic growth rate tends to increase and the secondary dendrite arms are more developed as the pressure is increased from 0.1 to 100 MPa, which showed a good agreement with the experimental results of direct squeeze casting of Mg-AI alloy. As the pressure increases, the largest dendritic growth rate can be obtained under the pressure between 200 and 250 MPa, while the growth rate decreases with a further increase of pressure.

Dendritic Growth Pattern and Dendritic Network Distortion in the Platform of a Ni-based Superalloy

The pattern of dendritic growth and distortion of dendritic network in the platform have been investigated by one mold casting with different platform length during directional solidification. As the platform length elongates, the symmetry of dendritic growth along left and right edges gradually worsens in platform base. While the platform length reaches 14 mm, the distortion of dendritic network is first observed in outward platform. It is found that the distortion of dendritic network along platform inside is more serious than that along platform edges. Both [001] deviation and accumulated misorientation along platform inside, up to 9~ and 16.3~, respectively, are far greeter than those along left-outward-right edges. The deformation of dendritic network in a platform may be caused by the asymmetry of the solidification front at the mush zone.

A three-dimensional cellular automaton model for simulation of dendritic growth of magnesium alloy

A numerical model based on the cellular automaton method for the three-dimensional simulation of dendritic growth of magnesium alloy was developed. The growth ki- netics was calculated from the complete solution of the transport equations. By con- structing a three-dimensional anisotropy model with the cubic CA cells, simulation of dendritic growth of magnesium alloy with six-fold symmetry in the basal plane was achieved. The model was applied to simulate the equiaxed dendritic growth and columnar dendritic growth under directional solidification, and its capability was addressed by comparing the simulated results to experimental results and those in the previously published works. Meanwhile, the three-dimensional simulated results were also compared with that of in two dimensions, offering a deep insight into the microstructure formation of magnesium alloy during solidification.

Phase-field simulation of tip splitting in dendritic growth of Fe-C alloy

Two types of dendrite tip splitting including dendrite orientation transition and twinned-like dendrites in Fe-C alloys were investigated by phase-field method. In equiaxed growth, the possible dendrite growth directions and the effect of supersaturation on tip splitting were discussed; the dendrite orientation transition was observed, and it was found that the orientation regions of anisotropy parameters were reduced from three to two with increasing the supersaturation, which was due to the effect of interracial anisotropy controlled by the solute in front of S/L interface changing with the increase of supersaturation. In directional solidification, it was found that the twinned like dendrites were formed with the fixed anisotropy couples and no seaweed dendrites were observed; these were concluded from the results of competition between process anisotropy and inherent anisotropy. The formation process of twinned-like dendrite was investigated by tip splitting phenomenon, which was related to the chan ges of dendrite tips growth velocity. Then, the critical speed of tips splitting and solute concentration of twinned-like dendrites were investigated, and a new type of microsegregation in Fe-C alloys was proposed to supplement the dendrite growth theories.

Phase-Field Modeling of Free Dendritic Growth in Binary Alloy under Forced Flow

A phase-field model(PFM)coupling with phase field,flow field and diffuse equation is presented for simulating isothermal dendrite growth of a nickel-copper alloy under a forced flow.Based on the finite difference method with uniform grid,the C programming code is implemented to complete the phase-field simulations.The simulation results indicate that the interfacial morphology,the symmetry of dendrite formation,the tip growth velocity and the concentration distribution are strongly influenced by the fluid flow.

Possibility of Solid Hydrometeor Growth Zone Identification Using Radar Spectrum Width

In this study,the correlation between simulated and measured radar velocity spectrum width(σ_(v))is investigated.The results show that the dendrites growth zones(DGZs)and needles growth zones(NGZs)mostly contain dendrites(DN)and needles(NE),respectively.Clearσ_(v) zones(1.1<σ_(v)(m s^(-1))<1.3 and 0.3<σ_(v)(m s^(-1))<0.7 for the DGZ and NGZ,respectively)could be identified in the case studies(27 and 28 February 2016)near altitudes corresponding to temperatures of–15°C and–5°C,according to the Japan Meteorological Agency and mesoscale model reanalysis data.Oblate particles with diverse particle shapes were observed in the DGZ withσ_(v)>1.2 m s^(-1),a differential reflectivity(ZDR)higher than 0 dB,and a cross-correlation coefficient(ρhv)less than 0.96.In contrast,prolate particles with relatively uniform shapes were observed in the NGZ withσ_(v)<0.6 m s^(-1),a ZDR less than 0 dB,andρhv higher than 0.97.The simulation results show that the DN exhibited a largerσ_(v) compared to the NE,and this observedσ_(v) was strongly dependent on the wind fluctuations(v’)due to turbulence or wind shear.In contrast,the NE exhibited a significantly smallσ_(v)~0.55 m s^(-1),which converges irrespective of v’.In addition,a strong correlation between the measuredσ_(v) values at five radar elevation angles(θ=6.2°,9.1°,13.1°,19°,and 80°)and those simulated in this study confirmed the significance of the analysis results.

Effect of traveling-wave magnetic field on dendrite growth of high-strength steel slab: Industrial trials and numerical simulation

The dendrite growth behavior of high-strength steel during slab continuous casting with a traveling-wave magnetic field was studied in this paper. The morphology of the solidification structure and composition distribution were analyzed. Results showed that the columnar crystals could deflect and break when the traveling-wave magnetic field had low current intensity. With the increase in current intensity, the secondary dendrite arm spacing and solute permeability decreased, and the columnar crystal transformed into an equiaxed crystal. The electromagnetic force caused by the traveling-wave magnetic field changed the temperature gradient and velocity magnitude and promoted the breaking and fusing of dendrites. Dendrite compactness and composition uniformity were arranged in descending order as follows:columnar-toequiaxed transition (high current intensity), columnar crystal zone (low current intensity), columnar-to-equiaxed transition (low current intensity), and equiaxed crystal zone (high current intensity). Verified numerical simulation results combined with the boundary layer theory of solidification front and dendrite breaking–fusing model revealed the dendrite deflection mechanism and growth process. When thermal stress is not considered, and no narrow segment can be found in the dendrite, the velocity magnitude on the solidification front of liquid steel can reach up to 0.041 m/s before the dendrites break.

Dendritic grain growth simulation in weld pool of nickel base alloy

Dendritic grain growth at the edge of the weld pool is simulated using a stochastic numerical model of cellular automaton algorithm. The grain growth miodel is established based upon the balance of solute in the solid/liquid interface of the dendrite tip. Considering the complicated nucleation condition and competitive growth, the dendrite moiphologies of different nucleation condition are simulated. The simulated results reproduced the dendrite grain evolution process at the edge of the weld pool. It is indicated that the nucleation condition is an important factor influencing the grain morphologies especially the morphologies of secondary and tertiary arms.