The formation and growth of Kirkendall voids in a binary alloy system during deformation process were investigated byphase field crystal model.The simulation results show that Kirkendall voids nucleate preferentially ...The formation and growth of Kirkendall voids in a binary alloy system during deformation process were investigated byphase field crystal model.The simulation results show that Kirkendall voids nucleate preferentially at the interface,and the averagesize of the voids increases with both the time and strain rate.There is an obvious coalescence of the voids at a large strain rate whenthe deformation is applied along the interface under both constant and cyclic strain rate conditions.For the cyclic strain rate appliedalong the interface,the growth exponent of Kirkendall voids increases with increasing the strain rate when the strain rate is largerthan1.0×10-6,while it increases initially and then decreases when the strain rate is smaller than9.0×10?7.The growth exponent ofKirkendall voids increases initially and then decreases gradually with increasing the length of cyclic period under a square-waveform constant strain rate.展开更多
The anisotropy problem of 3D phase-field model was studied,and various degrees of anisotropy were simulated by numerical calculation method.The results show that with the change of interface anisotropy coefficients,fr...The anisotropy problem of 3D phase-field model was studied,and various degrees of anisotropy were simulated by numerical calculation method.The results show that with the change of interface anisotropy coefficients,from smooth transition to the appearance of angle,equilibrium crystals shape morphology has a critical value,and 3D critical value is 0.3.The growth of dendrites is stable and the interface is smooth when it is less than critical value;the interface is unstable,rolling edge appears and the growth is discontinuous when it is more than critical value.With the increase of anisotropy coefficients,the dendrites grow faster under the same condition.展开更多
In this article, we study the phase-field model of solidification for numerical simulation of dendritic crystal growth that occurs during the casting of metals and alloys. Phase-field model of solidification describes...In this article, we study the phase-field model of solidification for numerical simulation of dendritic crystal growth that occurs during the casting of metals and alloys. Phase-field model of solidification describes the physics of dendritic growth in any material during the process of under cooling. The numerical procedure in this work is based on finite difference scheme for space and the 4th-order Runge-Kutta method for time discretization. The effect of each physical parameter on the shape and growth of dendritic crystal is studied and visualized in detail.展开更多
In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)...In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)equation.In which,the first order linear scheme is based on the invariant energy quadratization approach.The MPFC equation is a damped wave equation,and to preserve an energy stability,it is necessary to introduce a pseudo energy,which all increase the difficulty of constructing numerical methods comparing with the phase field crystal(PFC)equation.Due to the severe time step restriction of explicit timemarchingmethods,we introduce the first order and second order semi-implicit schemes,which are proved to be unconditionally energy stable.In order to improve the temporal accuracy,the semi-implicit spectral deferred correction(SDC)method combining with the first order convex splitting scheme is employed.Numerical simulations of the MPFC equation always need long time to reach steady state,and then adaptive time-stepping method is necessary and of paramount importance.The schemes at the implicit time level are linear or nonlinear and we solve them by multigrid solver.Numerical experiments of the accuracy and long time simulations are presented demonstrating the capability and efficiency of the proposed methods,and the effectiveness of the adaptive time-stepping strategy.展开更多
We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification.Here we focus on Stefan problems,and their quasi-static variants,with ...We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification.Here we focus on Stefan problems,and their quasi-static variants,with applications to crystal growth.New approaches with a high mesh quality for the parametric approximations of the resulting free boundary problems and new stable discretizations of the anisotropic phase field system are taken into account in a comparison involving benchmark problems based on exact solutions of the free boundary problem.展开更多
基金Projects(51275178,51405162,51205135) supported by the National Natural Science Foundation of ChinaProjects(20110172110003,20130172120055) supported by the Doctoral Program of Higher Education of China
文摘The formation and growth of Kirkendall voids in a binary alloy system during deformation process were investigated byphase field crystal model.The simulation results show that Kirkendall voids nucleate preferentially at the interface,and the averagesize of the voids increases with both the time and strain rate.There is an obvious coalescence of the voids at a large strain rate whenthe deformation is applied along the interface under both constant and cyclic strain rate conditions.For the cyclic strain rate appliedalong the interface,the growth exponent of Kirkendall voids increases with increasing the strain rate when the strain rate is largerthan1.0×10-6,while it increases initially and then decreases when the strain rate is smaller than9.0×10?7.The growth exponent ofKirkendall voids increases initially and then decreases gradually with increasing the length of cyclic period under a square-waveform constant strain rate.
基金Project(110218)supported by Program for the Top Young Academic Leaders of Higher Learning Institutions of Shanxi Province,China
文摘The anisotropy problem of 3D phase-field model was studied,and various degrees of anisotropy were simulated by numerical calculation method.The results show that with the change of interface anisotropy coefficients,from smooth transition to the appearance of angle,equilibrium crystals shape morphology has a critical value,and 3D critical value is 0.3.The growth of dendrites is stable and the interface is smooth when it is less than critical value;the interface is unstable,rolling edge appears and the growth is discontinuous when it is more than critical value.With the increase of anisotropy coefficients,the dendrites grow faster under the same condition.
文摘In this article, we study the phase-field model of solidification for numerical simulation of dendritic crystal growth that occurs during the casting of metals and alloys. Phase-field model of solidification describes the physics of dendritic growth in any material during the process of under cooling. The numerical procedure in this work is based on finite difference scheme for space and the 4th-order Runge-Kutta method for time discretization. The effect of each physical parameter on the shape and growth of dendritic crystal is studied and visualized in detail.
基金Research of R.Guo is supported by NSFC grant No.11601490Research of Y.Xu is supported by NSFC grant No.11371342,11626253,91630207.
文摘In this paper,we will develop a first order and a second order convex splitting,and a first order linear energy stable fully discrete local discontinuous Galerkin(LDG)methods for the modified phase field crystal(MPFC)equation.In which,the first order linear scheme is based on the invariant energy quadratization approach.The MPFC equation is a damped wave equation,and to preserve an energy stability,it is necessary to introduce a pseudo energy,which all increase the difficulty of constructing numerical methods comparing with the phase field crystal(PFC)equation.Due to the severe time step restriction of explicit timemarchingmethods,we introduce the first order and second order semi-implicit schemes,which are proved to be unconditionally energy stable.In order to improve the temporal accuracy,the semi-implicit spectral deferred correction(SDC)method combining with the first order convex splitting scheme is employed.Numerical simulations of the MPFC equation always need long time to reach steady state,and then adaptive time-stepping method is necessary and of paramount importance.The schemes at the implicit time level are linear or nonlinear and we solve them by multigrid solver.Numerical experiments of the accuracy and long time simulations are presented demonstrating the capability and efficiency of the proposed methods,and the effectiveness of the adaptive time-stepping strategy.
文摘We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification.Here we focus on Stefan problems,and their quasi-static variants,with applications to crystal growth.New approaches with a high mesh quality for the parametric approximations of the resulting free boundary problems and new stable discretizations of the anisotropic phase field system are taken into account in a comparison involving benchmark problems based on exact solutions of the free boundary problem.