[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driv...[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driving forces derived from the volumetric-deviatoric strain decomposition strategy,incorporating distinct tension,compression,and shear degradation mechanisms.Inertia effects capture compaction-band formation driven by wave-like disturbances,grain crushing,and frictional rearrangement.A monolithic algorithm ensures numerical stability and rapid convergence.[Results]The framework reproduces tensile,shear,mixed tensile-shear,and compressive-shear failures using the Benzeggagh-Kenane criterion.Validation against benchmark simulations-including uniaxial compression of rock-like and triaxial compression of V-notched sandstone specimens-demonstrates accurate predictions of crack initiation stress,localization orientation,and energy dissipation.[Conclusions]The framework provides a unified and robust numerical tool for analyzing the spatiotemporal evolution of strain localization and fracture in geomaterials.[Significance]By linking microscale fracture dynamics with macroscale failure within a thermodynamically consistent scheme,this study advances predictive modeling of rock stability,slope failure,and subsurface energy systems,contributing to safer and more sustainable geotechnical practice.展开更多
Based on the principles of thermodynamics, we elucidate the fundamental reasons behind the hysteresis of spontaneous polarization in ferroelectric materials during heating and cooling processes. By utilizing the effec...Based on the principles of thermodynamics, we elucidate the fundamental reasons behind the hysteresis of spontaneous polarization in ferroelectric materials during heating and cooling processes. By utilizing the effective Hamiltonian method in conjuction with the phase-field model, we have successfully reproduced the thermal hysteresis observed in ferroelectric materials during phase transitions. The computational results regarding the electrocaloric effect from these two different computational scales closely align with experimental measurements. Furthermore, we analyze how the first-order ferroelectric phase transition gradually diminishes with an increasing applied electric field, exhibiting characteristics of second-order-like phase transition. By employing the characteristic parameters of thermal hysteresis, we have established a pathway for calculations across different computational scales, thereby providing theoretical support for further investigations into the properties of ferroelectric materials through concurrent multiscale simulations.展开更多
Concrete materials are employed extensively in a variety of large-scale structures due to their economic viability and superior mechanical properties.During the service life of concrete structures,they are inevitably ...Concrete materials are employed extensively in a variety of large-scale structures due to their economic viability and superior mechanical properties.During the service life of concrete structures,they are inevitably subjected to damage from impact loading from natural disasters,such as earthquakes and storms.In recent years,the phasefield model has demonstrated exceptional capability in predicting the stochastic initiation,propagation,and bifurcation of cracks in materials.This study employs a phase-field model to focus on the rate dependency and failure response of concrete under impact deformation.A viscosity coefficient is introduced within the phase-field model to characterize the viscous behavior of dynamic crack propagation in concrete.The rate-dependent cohesive strength is defined within the yield function of concrete,where the rate sensitivity of cohesive strength facilitates the accumulation of the plastic driving force in the phase-field model.This process effectively captures the impact failure response of concrete.The applicability of the model was validated through unit cell experiments and numerical simulations of concrete under impact compression.Furthermore,the mechanical response and damage evolution mechanisms of concrete under impact loading were analyzed.It was observed that crack propagation in concrete initiates at material defects and,with increasing load,eventually develops in a direction perpendicular to the loading axis.展开更多
A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution t...A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.展开更多
A phase-field model integrated with the thermodynamic databases was constructed to investigate the impact of Ni content on the precipitation kinetics and phase transformation of the Cu-rich phase in Fe-Cu-Ni alloy at ...A phase-field model integrated with the thermodynamic databases was constructed to investigate the impact of Ni content on the precipitation kinetics and phase transformation of the Cu-rich phase in Fe-Cu-Ni alloy at 773 K.The results demonstrated that the Cu core-Ni shell structures form via the decomposition of Cu-Ni co-clusters,which is consistent with previous experimental results.As the Ni content increases,both the volume fraction and number density of Cu-rich precipitates increase,while their size decreases.With the increase in Ni content,the transformation from a Cu to 9R Cu is accelerated,which is the opposite to the result of increasing Mn content.Magnetic energy can increase the nucleation rate of the Cu-rich phase,but it does not affect the phase transformation driving force required for its crystal structure transformation.展开更多
Magnesium is distinguished by its highly anisotropic inelastic deformation involving a profuse activity of deformation twinning.Instrumented micro/nano-indentation technique has been widely applied to characterize the...Magnesium is distinguished by its highly anisotropic inelastic deformation involving a profuse activity of deformation twinning.Instrumented micro/nano-indentation technique has been widely applied to characterize the mechanical properties of magnesium,typically through the analysis of the indentation load-depth response,surface topography,and less commonly,the post-mortem microstructure within the bulk material.However,experimental limitations prevent the real-time observation of the evolving microstructure.To bridge this gap,we employ a recently-developed finite-strain model that couples the phase-field method and conventional crystal plasticity to simulate the evolution of the indentation-induced twin microstructure and its interaction with plastic slip in a magnesium single-crystal.Particular emphasis is placed on two aspects:orientation-dependent inelastic deformation and indentation size effects.Several outcomes of our 2D computational study are consistent with prior experimental observations.Chief among them is the intricate morphology of twin microstructure obtained at large spatial scales,which,to our knowledge,represents a level of detail that has not been captured in previous modeling studies.To further elucidate on size effects,we extend the model by incorporating gradient-enhanced crystal plasticity,and re-examine the notion of‘smaller is stronger’.The corresponding results underscore the dominant influence of gradient plasticity over the interfacial energy of twin boundaries in governing the size-dependent mechanical response.展开更多
Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilt...Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.展开更多
Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field con...Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field controlled equations. And the alternating direction implicit(ADI) algorithm for solving temperature field controlled equation was also employed to avoid the restriction of time step. Some characteristics of the Ni-Cu alloy were captured in the process of non-isothermal solidification, and the comparative analysis of the isothermal and the non-isothermal solidification was investigated. The simulation results indicate that the non-isothermal model is favorable to simulate the real solidification process of binary alloys, and when the thermal diffusivity decreases, the non-isothermal phase-field model is gradually consistent with the isothermal phase-field model.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
By incorporating two different fracture mechanisms and salient unilateral effects in rock materials,we propose a thermomechanical phase-field model to capture thermally induced fracture and shear heating in the proces...By incorporating two different fracture mechanisms and salient unilateral effects in rock materials,we propose a thermomechanical phase-field model to capture thermally induced fracture and shear heating in the process of rock failure.The heat conduction equation is derived,from which the plastic dissipation is treated as a heat source.We then ascertain the effect of the non-associated plastic flow on frictional dissipation and show how it improves the predictive capability of the proposed model.Taking advantage of the multiscale analysis,we propose a phase-field-dependent thermal conductivity with considering the unilateral effect of fracture.After proposing a robust algorithm for solving involved three-field coupling and damage-plasticity coupling problems,we present three numerical examples to illustrate the abilities of our proposed model in capturing various thermo-mechanically coupled behaviors.展开更多
A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the dr...A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.展开更多
The early precipitation process of Ni（75）Al（14）Mo（11） alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and oc...The early precipitation process of Ni（75）Al（14）Mo（11） alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and occupation probability of the three kinds of atoms were investigated.It is indicated that the non-stoichiometric Ll0（Ⅰ/Ⅱ） phases are found in the precipitation process.With the temperature increasing,the appearance time of Ll0 is brought forward.The Ll0（Ⅱ） structure always precipitates earlier than the Ll0（Ⅰ） structure.Compared with lower temperature,higher temperature brings the formation time of Ll0 phase forward and makes Ll0 phase have a higher order degree.But lower temperature shortens the process time of the Ll0 phase to the Ll2 phase.Al and Mo atoms tend to occupy γ site,Ni atom tends to occupy a and β sites.At the same temperature,Al atom has stronger occupation ability than Mo atom in the same site.Ni,Al and Mo collectively form the composited Ll2 structure.展开更多
Solidification structure is a key aspect for understanding the mechanical performance of metal alloys,wherein composition and casting parameters considerably influence solidification and determine the unique microstru...Solidification structure is a key aspect for understanding the mechanical performance of metal alloys,wherein composition and casting parameters considerably influence solidification and determine the unique microstructure of the alloys.By following the principle of free energy minimization,the phase-field method eliminates the need for tracking the solid/liquid phase interface and has greatly accelerated the research and development efforts geared toward optimizing metal solidification microstructures.The recent progress in the application of phasefield simulation to investigate the effect of alloy composition and casting process parameters on the solidification structure of metals is summarized in this review.The effects of several typical elements and process parameters,including carbon,boron,silicon,cooling rate,pulling speed,scanning speed,anisotropy,and gravity,on the solidification structure are discussed.The present work also addresses the future prospects of phase-field simulation and aims to facilitate the widespread applications of phase-field approaches in the simulation of microstructures during solidification.展开更多
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 comput...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.展开更多
A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is...A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is applied to the grain growth process of polycrystalline Mg(ZK60)alloy to investigate the evolution characteristics in different systems with varying proportions of low-angle grain boundary(LAGB)caused by different texture levels.It is found that the GB energy anisotropy can cause the grain growth kinetics to change,namely,higher texture levels(also means higher LAGB proportion)result in lower kinetics,and vice versa.The simulation results also show that the topological characteristics,such as LAGB proportion and distribution of grain size,undergo different evolution characteristics in different systems,and a more serious grain size fluctuation can be caused by a higher texture level.The mechanism is mainly the slower evolution of textured grains in their accumulation area and the faster coarsening rate of non-textured grains.Therefore,weakening the texture level is an effective way for implementing a desired homogenized microstructure in ZK60 Mg alloy.The rules revealed by the simulation results should be of great significance for revealing how the GB anisotropy affects the evolution of polycrystalline during the grain growth after recrystallization and offer the ideas for processing the alloy and optimizing the microstructure.展开更多
Hydrogen embrittlement of steels is directly linked to hydrogen diffusion and trapping in the microstructure,which can hardly be precisely measured by modern experimental techniques.A phase-field model,in which a chem...Hydrogen embrittlement of steels is directly linked to hydrogen diffusion and trapping in the microstructure,which can hardly be precisely measured by modern experimental techniques.A phase-field model,in which a chemical potential well of hydrogen in the grain boundaries is introduced,is proposed to simulate hydrogen diffusion and trapping in the polycrystalline iron.It was interestingly found that grain boundaries,as connected trap sites,have a complex influence on the effective diffusivity of hydrogen,which are strongly linked to grain boundary diffusivity and binding energy.展开更多
The objective of the work is focused on predictions of microsegregation,solidification speed,dendritic arm spacings and dendrite morphology by phase-field model.The numerical results were compared with experimental da...The objective of the work is focused on predictions of microsegregation,solidification speed,dendritic arm spacings and dendrite morphology by phase-field model.The numerical results were compared with experimental data.The experimental values for cooling rates and effective partition coefficient were adopted during calculations.The results of microsegregation through phase-field model show excellent agreement with the experimental data.Such excellent agreement is because cooling rates,effective partition coefficient and back-diffusion of solute are considered in the model.For solidification speed,the calculation results show good agreement with the experimental data.Tertiary dendritic arm spacing calculated with phase-field model is compared with experimental data.The results show good agreement between them.The dendrite arm spacing varies with position because high cooling rates are responsible for the refinement effect on microstructure.Finally,two-dimensional simulation produced a dendrite that is similar to that found in the experiment.展开更多
The quantitative phase-field simulations were reviewed on the processes of solidification of pure metals and alloys.The quantitative phase-field equations were treated in a diffuse thin-interface limit,which enabled t...The quantitative phase-field simulations were reviewed on the processes of solidification of pure metals and alloys.The quantitative phase-field equations were treated in a diffuse thin-interface limit,which enabled the quantitative links between interface dynamics and model parameters in the quasi-equilibrium simulations.As a result,the quantitative modeling is more effective in dealing with microstructural pattern formation in the large scale simulations without any spurious kinetic effects.The development of the quantitative phase-field models in modeling the formation of microstructures such as dendritic structures,eutectic lamellas,seaweed morphologies,and grain boundaries in different solidified conditions was also reviewed with the purpose of guiding to find the new prospect of applications in the quantitative phase-field simulations.展开更多
The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and...The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and most importantly, no paradoxical parameters appeared in the model. Moreover, it can simulate different microstructure patterns owing to the use of a new different free energy function for the simulation of morphologies of polymer. The new free energy function considers both the cases of T〈Tm and T≥Tm, which is more reasonable than that in published literatures that all ignored the T≥Tm case. In order to show the validity of the modified model, the finite difference method is used to solve the model and different crystallization morphologies during the solidification process of isotactic polystyrene are obtained under different conditions. Numerical results show that the growth rate of the initial secondary arms is obviously increased as the anisotropy strength increases. But the anisotropy strength seems to have no apparent effect on the global growth rate. The whole growth process of the dendrite depends mainly upon the latent heat and the latent heat has a direct effect on the tip radius and tip velocity of side branches.展开更多
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 meth...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.展开更多
文摘[Objective]This study aims to develop a thermodynamically consistent phase-field framework for modeling the initiation and evolution of discontinuous structures in geomaterials.[Methods]Our model introduces crack driving forces derived from the volumetric-deviatoric strain decomposition strategy,incorporating distinct tension,compression,and shear degradation mechanisms.Inertia effects capture compaction-band formation driven by wave-like disturbances,grain crushing,and frictional rearrangement.A monolithic algorithm ensures numerical stability and rapid convergence.[Results]The framework reproduces tensile,shear,mixed tensile-shear,and compressive-shear failures using the Benzeggagh-Kenane criterion.Validation against benchmark simulations-including uniaxial compression of rock-like and triaxial compression of V-notched sandstone specimens-demonstrates accurate predictions of crack initiation stress,localization orientation,and energy dissipation.[Conclusions]The framework provides a unified and robust numerical tool for analyzing the spatiotemporal evolution of strain localization and fracture in geomaterials.[Significance]By linking microscale fracture dynamics with macroscale failure within a thermodynamically consistent scheme,this study advances predictive modeling of rock stability,slope failure,and subsurface energy systems,contributing to safer and more sustainable geotechnical practice.
基金Project supported financially by the National Natural Science Foundation of China (Grant No. 52372100)the National Key Research and Development Program of China (Grant No. 2019YFA0307900)。
文摘Based on the principles of thermodynamics, we elucidate the fundamental reasons behind the hysteresis of spontaneous polarization in ferroelectric materials during heating and cooling processes. By utilizing the effective Hamiltonian method in conjuction with the phase-field model, we have successfully reproduced the thermal hysteresis observed in ferroelectric materials during phase transitions. The computational results regarding the electrocaloric effect from these two different computational scales closely align with experimental measurements. Furthermore, we analyze how the first-order ferroelectric phase transition gradually diminishes with an increasing applied electric field, exhibiting characteristics of second-order-like phase transition. By employing the characteristic parameters of thermal hysteresis, we have established a pathway for calculations across different computational scales, thereby providing theoretical support for further investigations into the properties of ferroelectric materials through concurrent multiscale simulations.
文摘Concrete materials are employed extensively in a variety of large-scale structures due to their economic viability and superior mechanical properties.During the service life of concrete structures,they are inevitably subjected to damage from impact loading from natural disasters,such as earthquakes and storms.In recent years,the phasefield model has demonstrated exceptional capability in predicting the stochastic initiation,propagation,and bifurcation of cracks in materials.This study employs a phase-field model to focus on the rate dependency and failure response of concrete under impact deformation.A viscosity coefficient is introduced within the phase-field model to characterize the viscous behavior of dynamic crack propagation in concrete.The rate-dependent cohesive strength is defined within the yield function of concrete,where the rate sensitivity of cohesive strength facilitates the accumulation of the plastic driving force in the phase-field model.This process effectively captures the impact failure response of concrete.The applicability of the model was validated through unit cell experiments and numerical simulations of concrete under impact compression.Furthermore,the mechanical response and damage evolution mechanisms of concrete under impact loading were analyzed.It was observed that crack propagation in concrete initiates at material defects and,with increasing load,eventually develops in a direction perpendicular to the loading axis.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.
基金supported by the National Natural Science Foundation of China(Grant No.51871086).
文摘A phase-field model integrated with the thermodynamic databases was constructed to investigate the impact of Ni content on the precipitation kinetics and phase transformation of the Cu-rich phase in Fe-Cu-Ni alloy at 773 K.The results demonstrated that the Cu core-Ni shell structures form via the decomposition of Cu-Ni co-clusters,which is consistent with previous experimental results.As the Ni content increases,both the volume fraction and number density of Cu-rich precipitates increase,while their size decreases.With the increase in Ni content,the transformation from a Cu to 9R Cu is accelerated,which is the opposite to the result of increasing Mn content.Magnetic energy can increase the nucleation rate of the Cu-rich phase,but it does not affect the phase transformation driving force required for its crystal structure transformation.
文摘Magnesium is distinguished by its highly anisotropic inelastic deformation involving a profuse activity of deformation twinning.Instrumented micro/nano-indentation technique has been widely applied to characterize the mechanical properties of magnesium,typically through the analysis of the indentation load-depth response,surface topography,and less commonly,the post-mortem microstructure within the bulk material.However,experimental limitations prevent the real-time observation of the evolving microstructure.To bridge this gap,we employ a recently-developed finite-strain model that couples the phase-field method and conventional crystal plasticity to simulate the evolution of the indentation-induced twin microstructure and its interaction with plastic slip in a magnesium single-crystal.Particular emphasis is placed on two aspects:orientation-dependent inelastic deformation and indentation size effects.Several outcomes of our 2D computational study are consistent with prior experimental observations.Chief among them is the intricate morphology of twin microstructure obtained at large spatial scales,which,to our knowledge,represents a level of detail that has not been captured in previous modeling studies.To further elucidate on size effects,we extend the model by incorporating gradient-enhanced crystal plasticity,and re-examine the notion of‘smaller is stronger’.The corresponding results underscore the dominant influence of gradient plasticity over the interfacial energy of twin boundaries in governing the size-dependent mechanical response.
基金Project supported by the National Natural Science Foundation of China(Nos.12072297 and12202370)the Natural Science Foundation of Sichuan Province of China(No.24NSFSC4777)。
文摘Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of China
文摘Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field controlled equations. And the alternating direction implicit(ADI) algorithm for solving temperature field controlled equation was also employed to avoid the restriction of time step. Some characteristics of the Ni-Cu alloy were captured in the process of non-isothermal solidification, and the comparative analysis of the isothermal and the non-isothermal solidification was investigated. The simulation results indicate that the non-isothermal model is favorable to simulate the real solidification process of binary alloys, and when the thermal diffusivity decreases, the non-isothermal phase-field model is gradually consistent with the isothermal phase-field model.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
基金funding provided by the National Natural Science Foundation of China(No.12202137)TY's contribution is funded by the China and Germany Postdoctoral Exchange Program(Grant No.ZD202137).The first author(TY)would like to express his gratitude to Prof.Keita Yoshioka for reviewing this manuscript and for his invaluable feedback.
文摘By incorporating two different fracture mechanisms and salient unilateral effects in rock materials,we propose a thermomechanical phase-field model to capture thermally induced fracture and shear heating in the process of rock failure.The heat conduction equation is derived,from which the plastic dissipation is treated as a heat source.We then ascertain the effect of the non-associated plastic flow on frictional dissipation and show how it improves the predictive capability of the proposed model.Taking advantage of the multiscale analysis,we propose a phase-field-dependent thermal conductivity with considering the unilateral effect of fracture.After proposing a robust algorithm for solving involved three-field coupling and damage-plasticity coupling problems,we present three numerical examples to illustrate the abilities of our proposed model in capturing various thermo-mechanically coupled behaviors.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.
基金Project(51275486)supported by the National Natural Science Foundation of China
文摘The early precipitation process of Ni（75）Al（14）Mo（11） alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and occupation probability of the three kinds of atoms were investigated.It is indicated that the non-stoichiometric Ll0（Ⅰ/Ⅱ） phases are found in the precipitation process.With the temperature increasing,the appearance time of Ll0 is brought forward.The Ll0（Ⅱ） structure always precipitates earlier than the Ll0（Ⅰ） structure.Compared with lower temperature,higher temperature brings the formation time of Ll0 phase forward and makes Ll0 phase have a higher order degree.But lower temperature shortens the process time of the Ll0 phase to the Ll2 phase.Al and Mo atoms tend to occupy γ site,Ni atom tends to occupy a and β sites.At the same temperature,Al atom has stronger occupation ability than Mo atom in the same site.Ni,Al and Mo collectively form the composited Ll2 structure.
基金financially supported by the National Key Research and Development Program of China(No.2021YFB3702401)the National Natural Science Foundation of China(Nos.51901013,52122408,52071023)+3 种基金financial support from the Fundamental Research Funds for the Central Universities,China(University of Science and Technology Beijing(USTB),Nos.FRF-TP-2021-04C1,06500135)financial support from the Qilu Young Talent Program of Shandong University,Zhejiang Lab Open Research Project,China(No.K2022PE0AB05)the Shandong Provincial Natural Science Foundation,China(No.ZR2023MA058)the Guangdong Basic and Applied Basic Research Foundation,China(No.2023A1515011819)。
文摘Solidification structure is a key aspect for understanding the mechanical performance of metal alloys,wherein composition and casting parameters considerably influence solidification and determine the unique microstructure of the alloys.By following the principle of free energy minimization,the phase-field method eliminates the need for tracking the solid/liquid phase interface and has greatly accelerated the research and development efforts geared toward optimizing metal solidification microstructures.The recent progress in the application of phasefield simulation to investigate the effect of alloy composition and casting process parameters on the solidification structure of metals is summarized in this review.The effects of several typical elements and process parameters,including carbon,boron,silicon,cooling rate,pulling speed,scanning speed,anisotropy,and gravity,on the solidification structure are discussed.The present work also addresses the future prospects of phase-field simulation and aims to facilitate the widespread applications of phase-field approaches in the simulation of microstructures during solidification.
基金financially supported by the National Natural Science Foundation of China(No.51805389)the Natural Science Foundation of Hubei Province,China(No.2018CFB210)the 111 Project(No.B17034)。
文摘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.
基金Project supported by the National Key Research and Development Program of China(Grant No.2016YFB0701204)the 111 Project,China(Grant No.B20029)+2 种基金the Fundamental Research Funds for the Central Universities,China(Grant Nos.N2002017 and N2007011)the National Natural Science Foundation of China(Grant No.51571055)the Science Fund from the Science and Technology Bureau of Jiangyin High-Tech Industrial Development Zone,China(Grant No.ZX20200062)。
文摘A three-dimensional(3D)multiple phase field model,which takes into account the grain boundary(GB)energy anisotropy caused by texture,is established based on real grain orientations and Read-Shockley model.The model is applied to the grain growth process of polycrystalline Mg(ZK60)alloy to investigate the evolution characteristics in different systems with varying proportions of low-angle grain boundary(LAGB)caused by different texture levels.It is found that the GB energy anisotropy can cause the grain growth kinetics to change,namely,higher texture levels(also means higher LAGB proportion)result in lower kinetics,and vice versa.The simulation results also show that the topological characteristics,such as LAGB proportion and distribution of grain size,undergo different evolution characteristics in different systems,and a more serious grain size fluctuation can be caused by a higher texture level.The mechanism is mainly the slower evolution of textured grains in their accumulation area and the faster coarsening rate of non-textured grains.Therefore,weakening the texture level is an effective way for implementing a desired homogenized microstructure in ZK60 Mg alloy.The rules revealed by the simulation results should be of great significance for revealing how the GB anisotropy affects the evolution of polycrystalline during the grain growth after recrystallization and offer the ideas for processing the alloy and optimizing the microstructure.
基金the National Key R&D program of China(No.2016YFB0300104)the financial support by the Beijing Natural Science Foundation(No.2182024)+1 种基金the National Natural Science Foundation of China(Grant Nos.51501099 and U1764252)the National Young 1000-Talents Program(No.D1101073)。
文摘Hydrogen embrittlement of steels is directly linked to hydrogen diffusion and trapping in the microstructure,which can hardly be precisely measured by modern experimental techniques.A phase-field model,in which a chemical potential well of hydrogen in the grain boundaries is introduced,is proposed to simulate hydrogen diffusion and trapping in the polycrystalline iron.It was interestingly found that grain boundaries,as connected trap sites,have a complex influence on the effective diffusivity of hydrogen,which are strongly linked to grain boundary diffusivity and binding energy.
基金FAPERJ(Fundacao de Amparo a Pesquisa do Estado do Rio de Janeiro,APQ#1:E-26/010.001942/2019)CNPq(Conselho Nacional de Desenvolvimento Cientifico e Tecnologico)for their financial support。
文摘The objective of the work is focused on predictions of microsegregation,solidification speed,dendritic arm spacings and dendrite morphology by phase-field model.The numerical results were compared with experimental data.The experimental values for cooling rates and effective partition coefficient were adopted during calculations.The results of microsegregation through phase-field model show excellent agreement with the experimental data.Such excellent agreement is because cooling rates,effective partition coefficient and back-diffusion of solute are considered in the model.For solidification speed,the calculation results show good agreement with the experimental data.Tertiary dendritic arm spacing calculated with phase-field model is compared with experimental data.The results show good agreement between them.The dendrite arm spacing varies with position because high cooling rates are responsible for the refinement effect on microstructure.Finally,two-dimensional simulation produced a dendrite that is similar to that found in the experiment.
基金supported by National Natural Science Foundation of China(No.51174177)the Fund of the State Key Solidification Laboratory of Solidification Processing in Northwestern Polytechnical University(No.SKLSP 201714)
文摘The quantitative phase-field simulations were reviewed on the processes of solidification of pure metals and alloys.The quantitative phase-field equations were treated in a diffuse thin-interface limit,which enabled the quantitative links between interface dynamics and model parameters in the quasi-equilibrium simulations.As a result,the quantitative modeling is more effective in dealing with microstructural pattern formation in the large scale simulations without any spurious kinetic effects.The development of the quantitative phase-field models in modeling the formation of microstructures such as dendritic structures,eutectic lamellas,seaweed morphologies,and grain boundaries in different solidified conditions was also reviewed with the purpose of guiding to find the new prospect of applications in the quantitative phase-field simulations.
基金This work is supported by the National Natural Science Foundation of China （No.11402210）, the Natural Science Foundation of Shanxi Province （No.2012011019-2）, and the Doctoral Fund of Taiyuan University of Science and Technology （No.20152024）.
文摘The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and most importantly, no paradoxical parameters appeared in the model. Moreover, it can simulate different microstructure patterns owing to the use of a new different free energy function for the simulation of morphologies of polymer. The new free energy function considers both the cases of T〈Tm and T≥Tm, which is more reasonable than that in published literatures that all ignored the T≥Tm case. In order to show the validity of the modified model, the finite difference method is used to solve the model and different crystallization morphologies during the solidification process of isotactic polystyrene are obtained under different conditions. Numerical results show that the growth rate of the initial secondary arms is obviously increased as the anisotropy strength increases. But the anisotropy strength seems to have no apparent effect on the global growth rate. The whole growth process of the dendrite depends mainly upon the latent heat and the latent heat has a direct effect on the tip radius and tip velocity of side branches.
基金the National Natural Science Foundation of China(Nos.10964004 and 50804019)the Natural Science Foundation of Gansu Province (No.1014RJZA003)
文摘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.
