Phase field modeling of multiple dendrite growth of Al-Si binary alloy under isothermal solidification

Phase field method offers the prospect of being able to perform realistic numerical experiments on dendrite growth in metallic systems. In this study, the growth process of multiple dendrites in AI-2-mole-%-Si binary alloy under isothermal solidification was simulated using phase field model. The simulation results showed the impingement of arbitrarily oriented crystals and the competitive growth among the grains during solidification. With the increase of growing time, the grains begin to coalesce and impinge the adjacent grains. When the dendrites start to impinge, the dendrite growth is obviously inhibited.

Phase field modeling of dendrite growth

Single dendrite and multi-dendrite growth for Al-2 mol pct Si alloy during isothermal solidification are simulated by phase field method. In the case of single equiaxed dendrite growth, the secondary and the necking phenomenon can be observed. For multi-dendrite growth, there exists the competitive growth among the dendrites during solidification. As solidification proceeds, growing and coarsening of the primary arms occurs, together with the branching and coarsening of the secondary arms. When the diffusion fields of dendrite tips come into contact with those of the branches growing from the neighboring dendrites, the dendrites stop growing and being to ripen and thicken.

Reversed-phase fused-core HPLC modeling of peptides

Different fused-core stationary phase chemistries（C18,Amide,Phenyl-hexyl and Peptide ES-C18） were used for the analysis of 21 structurally representative model peptides.In addition,the effects of the mobile phase composition（ACN or MeOH as organic modifier;formic acid or acetic acid,as acidifying component） on the column selectivity,peak shape and overall chromatographic performance were evaluated.The RP-amide column,combined with a formic acid-acetonitrile based gradient system,performed as best.A peptide reversed-phase retention model is proposed,consisting of 5 variables：log SumAA,log Sv,clog P,log nHDon and log nHAcc.Quantitative structure-retention relationship（QSRR） models were constructed for 16 different chromatographic systems.The accuracy of this peptide retention model was demonstrated by the comparison between predicted and experimentally obtained retention times,explaining on average 86% of the variability.Moreover,using an external set of 5 validation peptides,the predictive power of the model was also demonstrated.This peptide retention model includes the novel in-silico calculated amino acid descriptor,AA,which was calculated from log P,3D-MoRSE,RDF and WHIM descriptors.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Phase field modeling of the ring-banded spherulites of crystalline polymers: The role of thermal diffusion

The ring-banded spherulite is a special morphology of polymer crystals and has attracted considerable attention over recent decades. In this study, a new phase field model with polymer characteristics is established to investigate the emergence and formation mechanism of the ring-banded spherulites of crystalline polymers. The model consists of a nonconserved phase field representing the phase transition and a temperature field describing the diffusion of the released latent heat. The corresponding model parameters can be obtained from experimentally accessible material parameters.Two-dimensional calculations are carried out for the ring-banded spherulitic growth of polyethylene film under a series of crystallization temperatures. The results of these calculations demonstrate that the formation of ring-banded spherulites can be triggered by the self-generated thermal field. Moreover, some temperature-dependent characteristics of the ring-banded spherulites observed in experiments are reproduced by simulations, which may help to study the effects of crystallization temperature on the ring-banded structures.

Phase field model for electric-thermal coupled discharge breakdown of polyimide nanocomposites under high frequency electrical stress

In contrast to conventional transformers, power electronic transformers, as an integral component of new energy power system, are often subjected to high-frequency and transient electrical stresses, leading to heightened concerns regarding insulation failures. Meanwhile, the underlying mechanism behind discharge breakdown failure and nanofiller enhancement under high-frequency electrical stress remains unclear. An electric-thermal coupled discharge breakdown phase field model was constructed to study the evolution of the breakdown path in polyimide nanocomposite insulation subjected to high-frequency stress. The investigation focused on analyzing the effect of various factors, including frequency, temperature, and nanofiller shape, on the breakdown path of Polyimide(PI) composites. Additionally, it elucidated the enhancement mechanism of nano-modified composite insulation at the mesoscopic scale. The results indicated that with increasing frequency and temperature, the discharge breakdown path demonstrates accelerated development, accompanied by a gradual dominance of Joule heat energy. This enhancement is attributed to the dispersed electric field distribution and the hindering effect of the nanosheets. The research findings offer a theoretical foundation and methodological framework to inform the optimal design and performance management of new insulating materials utilized in high-frequency power equipment.

Bifurcation Analysis Reveals Solution Structures of Phase Field Models

The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems.Here,we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models.Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions.To elucidate the idea,we apply this analytical approach to three representative phase field equations:the Allen-Cahn equation,the Cahn-Hilliard equation,and the Allen-Cahn-Ohta-Kawasaki system.The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

Thermodynamic consistent phase field model for sintering process with multiphase powders

A thermodynamic consistent phase field model is developed to describe the sintering process with multiphase powders. In this model, the interface region is assumed to be a mixture of different phases with the same chemical potential, but with different compositions. The interface diffusion and boundary diffusion are also considered in the model. As an example, the model is applied to the sintering process with Fe-Cu powders. The free energy of each phase is described by the well-developed thermodynamic models, together with the published optimized parameters. The microstructure and solute distribution during the sintering process can both be obtained quantitively.

Phase field calculation of interface mobility in a ternary alloy

A novel approach based on the quantitative phase field model was proposed to calculate the interface mobility and applied to the α/β interface of a ternary Ti-6Al-4V alloy.Phase field simulations indicate that the higher interface mobility leads to the faster transformation rate,but only a unique value of interface mobility matches the diffusion equation under the diffusion-controlled condition.By comparing the transformation kinetics from phase field simulations with that from classical diffusion equation,the interface mobility at different temperatures can be obtained.The results show that the calculated interface mobility increases with increasing temperature and accords with Arrhenius equation very well.

Phase-field study for effect of chemical composition on pre-precipitation process of γ′ phase in Ni-Al-Ti alloy

The pre-precipitation processes of γ＇ phase in Ni-Al-Ti alloy with different chemical compositions are investigated by microscopic phase-field dynamics model. The simulation results show that the early precipitated phase with L10 structure is formed before the formation of γ＇ phase. With the increase of Al content, the composition of precipitated phase is transformed from the co-existence of γ＇ phase and material phase into single y phase. Furthermore, excessive Al content can stimulate the production of γ＇ phase, and the conversion time from L10 phase to L12 phase will lag behind other samples with lower Al con-tent. By analyzing the occupation probability of Ti, it can be found that the solid solution strengthening effect of Ti is reduced with the increase of Al component.

Phase-field-lattice Boltzmann study for lamellar eutectic growth in a natural convection melt

In the present study, the influence of natural convection on the lamellar eutectic growth is determined by a phase-field-lattice Boltzmann study for Al-Cu eutectic alloy. The mass difference resulting from concentration difference led to the fluid flow, and a robust parallel and adaptive mesh refinement algorithm was employed to improve the computational efficiency without any compromising accuracy. Results show that the existence of natural convection would affect the growth undercooling and thus control the interface shape by adjusting the lamellar width. In particular, by alternating the magnitude of the solute expansion coefficient, the strength of the natural convection is changed. Corresponding microstructure patterns are discussed and compared with those under no-convection conditions.

Phase Equilibrium of Isobutanol in Supercritical CO2

Vapor-liquid phase equilibrium data including composition,densities,molar volume and equilibrium constant of isobutanol in supercritical carbon dioxide from 313.2K to 353.2K were measured in a variable-volume visual cell.The properties of critical point were obtained by extrapolation.The results showed that critical temperature,critical pressure and critical compressibility factor of CO2-isobutanol system decreased with the increase of critical CO2 content.The phase equilibrium model was established by Peng-Robinson equation of state and van der Waals-2 mixing regulation,and model parameters were determined by optimization calculation of nonlinear least square method.The correlation between calculated values and the experimental data showed good agreement.

On the Liquid-Vapor Phase-Change Interface Conditions for Numerical Simulation of Violent Separated Flows

Numerous models have been proposed in the literature to include phase change into numerical simulations of two-phase flows.This review paper presents the modeling options that have been taken in order to obtain a model for violent separated flows with application to sloshing wave impacts.A relaxation model based on linear non-equilibrium thermodynamics has been chosen to compute the rate of phase change.The integration in the system of partial differential equations is done through a non-conservative advection term.For each of these modelling choices,some alternative models from the literature are presented and discussed.The theoretical framework for all phase change model(conservation equations and entropy growth)is also summarized.

DETERMINATION OF EFFECTIVE MODULI FOR FOAM PLASTICS BASED ON THREE PHASE SPHEROIDAL MODEL

The relations of bulk modulus, shear modulus, Young's modulus and the Poisson's ratio with porosity of foam plastics are determined by a three phase spheroidal model commonly used in Composite Mechanics. The results are compared with those using differential scheme. It is shown that the material properties derived from the present model normally are larger than those obtained by differential scheme for foam plastics with identical porosity. The differences in shear moduli and Young's moduli obtained by the two methods are small but they are larger for bulk moduli of incompressible matrix and Poisson's ratios. The Young's moduli of high density foam plastics derived by the present model agree better with experimental ones.

A Phase-field Model to Simulate Recrystallization in an AZ31 Mg Alloy in Comparison of Experimental Data

A model has been established to simulate the realistic spatio-temporal microstructure evolution in recrystallization of a magnesium alloy using the phase field approach. A set of rules have been proposed to decide the real physical value of all parameters in the model. The thermodynamic software THERMOCALC is applied to determine the local chemical free energy and strain energy, which is added to the free energy density of grains before recrystallization. The Arrhenius formula is used to describe boundary mobility and the activity energy is suggested with a value of zinc segregation energy at the boundary. However, the mobility constant in the formula was found out by fitting to a group of grain size measurements during recrystallization of the alloy. The boundary range is suggested to decide the gradient parameters in addition of fitting to the experimental boundary energy value. These parameter values can be regarded as a database for other similar simulations and the fitting rules can also be applied to build up databases for any other alloy systems. The simulated results show a good agreement with reported experimental measurement of the alloy at the temperatures from 300 to 400℃ for up to 100 min but not at 250℃. This implies a mechanism variation in activity energy of the boundary mobility in the alloy at low temperature.

Modeling the influence of forced ventilation on the dispersion of droplets ejected from roadheader-mounted external sprayer

In order to reveal the influence of forced ventilation on the dispersion of droplets ejected from roadheader-mounted external sprayer,the paper studies the air-flowing field and the droplet distribution under the condition of gentle breeze and normal forced ventilation in heading face using the particle tracking technology of computational fluid dynamics(CFD).The results show that air-flowing tendency in the same section presents great comparability in the period of gentle breeze and forced ventilation,and the difference mainly embodies in the different wind velocity.The influence of ventilation on the dispersion of droplets is faint under the gentle breeze condition.The droplet can be evenly distributed around the cutting head.However,under the normal forced ventilation,a large number of droplets will drift to the return air side.At the same time,droplet clusters are predominantly presented in the lower part of windward side and the middle of the leeward side around the cutting head.In contrast,the droplet concentration in other parts around cutting head decreases a lot and the droplets are unable to form close-grained mist curtain.So the dust escape channel is formed.In addition,the simulation results also reveal that the disturbance of air flow on the droplet distribution can be effectively relieved when using ventilation duct with Coanda effect(VDCE).Field experiment results show that the dust suppression efficiency of total dust and respirable dust increases respectively by 10.5%and 9.3%when using VDCE,which proves that it can weaken the influence of airflow on droplet dispersion.

Phase field simulation of the columnar dendritic growth and microsegregation in a binary alloy

This paper applies a phase field model for polycrystalline solidification in binary alloys to simulate the formation and growth of the columnar dendritic array under the isothermal and constant cooling conditions. The solidification process and microsegregation in the mushy zone are analysed in detail. It is shown that under the isothermal condition solidification will stop after the formation of the mushy zone, but dendritic coarsening will progress continuously, which results in the decrease of the total interface area. Under the constant cooling condition the mushy zone will solidify and coarsen simultaneously. For the constant cooling solidification, microsegregation predicted by a modified Brody- Flemings model is compared with the simulation results. It is found that the Fourier number which characterizes microsegregation is different for regions with different microstructures. Dendritic coarsening and the larger area of interface should account for the enhanced Fourier number in the region with well developed second dendritic arms.

Fracture propagation laws of staged hydraulic fracture in fractured geothermal reservoir based on phase field model

Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.

MICRO-DESCRIPTION OF THE SOLUTE-FIELD AND THE PHASE-FIELD MODEL FOR ISOTHERMAL PHASE TRANSITION IN BINARY ALLOYS

A new phase field method for two-dimensional simulations of binary alloy solidification was studied. A model basing on solute conservative in every unit was developed for solving the solute diffusion equation during solidification. Two-dimensional computations were performed for ideal solutions and Ni-Cu dendritic growth into an isothermal and highly supersaturated liquid phase.

Prediction and programming of microemulsion phase behavior simulation

In the process of microemulsion flooding,microemulsion phase state may be affected by the adsorption of core and the distribution of oil and water,and the upper phase,the middle phase or the lower phase microemulsion will appear.Accurate description of microemulsion phase state and quantitative discrimination of equilibrium phase composition are of great theoretical value and research significance for the design of microemulsion flooding system and improving the recovery efficiency of microemulsion.Therefore,in this paper,microemulsion phase model is deeply studied.Aiming at the difficulty of solving the existing Hand model and the unknown parameters of the improved HLD-NAC model,we introduce the binodal curve range parameter D and the asymmetric migration degree parameter B,and establish a new phase behavior description method of the binodal curve,the two-phase plait point lines and the III-phase node line.Then the phase discrimination programming of microemulsion was carried out,and the change laws of phase number,phase type and mass fraction of equilibrium phases of any microemulsion system under different salinity were revealed.The results show that the simulation results of microemulsions accord with the salinity scanning law,and can accurately identify the complex microemulsion phase states.