The experiments of laser fusion welding with Al-foil addition was carried out for DP590 dual-phase steel and AZ31B magnesium alloy in an overlap steel-on-magnesium configuration.Temperature field was simulated by COMS...The experiments of laser fusion welding with Al-foil addition was carried out for DP590 dual-phase steel and AZ31B magnesium alloy in an overlap steel-on-magnesium configuration.Temperature field was simulated by COMSOL finite element software for steel/magnesium laser fusion welding.The results show that when Al-foil is added,some defects,such as pores,cracks and softening in heat affected zone(HAZ),can be avoided in welding joint,the bonding strength of steel/magnesium joints is increased,heat transfer between steel and magnesium is regulated.In the case of adding Al-foil,welding pool is divided into two parts,the upper and lower pools contact each other but do not mix,the transition layer at the interface between the upper and lower molten pools mainly contains Al−Fe phases,such as AlFe,Al2Fe and AlFe3,and these new phases are helpful for promoting the metallurgical connection between the upper and lower molten pools.Hence,adding Al-foil laser fusion welding is an effective way in joining steel to magnesium alloy.展开更多
A new severe plastic deformation method for manufacturing tubes made of AZ31 magnesium alloy with a large diameter was developed,which is called the TCESE(tube continuous extrusion−shear−expanding)process.The process ...A new severe plastic deformation method for manufacturing tubes made of AZ31 magnesium alloy with a large diameter was developed,which is called the TCESE(tube continuous extrusion−shear−expanding)process.The process combines direct extrusion with a two-step shear−expanding process.The influences of expanding ratios,extrusion temperatures on the deformation of finite element meshes,strain evolution and flow velocity of tube blanks during the TCESE process were researched based on numerical simulations by using DEFORM-3D software.Simulation results show that the maximum expanding ratio is 3.0 in the TCESE process.The deformation of finite element meshes of tube blanks is inhomogeneous in the shear−expanding zone,and the equivalent strains increase significantly during the TCESE process of the AZ31 magnesium alloy.A extrusion temperature of 380°C and expanding ratio of 2.0 were selected as the optimized process parameters from the numerical simulation results.The average grain size of tubes fabricated by the TCESE process is approximately 10μm.The TCESE process can refine grains of magnesium alloy tubes with the occurrence of dynamic recrystallization.The(0001)basal texture intensities of the magnesium alloy tube blanks decrease due to continuous plastic deformation during the TCESE process.The average hardness of the extruded tubes is approximately HV 75,which is obviously improved.展开更多
With the microscopic phase-field dynamic model, the effects of temperature and concentration on the nucleation incubation time of Ni75AlxV25-x alloy were studied and the relation between the incubation time and precip...With the microscopic phase-field dynamic model, the effects of temperature and concentration on the nucleation incubation time of Ni75AlxV25-x alloy were studied and the relation between the incubation time and precipitation mechanism was investigated by using the atomic occupation probability picture and average order parameter curve. The simulation results demonstrate that there exists the incubation time for different precipitation mechanisms~ such as non-classical nucleation, the mixed style of non-classical nucleation and spinodal decomposition, and spinodal ordering; and the incubation time shortens in turn for the three kinds of mechanisms. With the increase of Al content of Ni75AlxV25-x alloy, the incubation time of Llz phases shortens continuously and that of DOzz phases is prolonged. The effects of temperature on the incubation time of Llz and DOzz phases are accordant, i.e. the incuba- tion time is greatly prolonged with the temperature rising.展开更多
The hot deformation behavior of Ti-6 Al-4 V-0.1 Ru titanium alloy was investigated by isothermal compression tests on a Gleeble-3500 thermal simulator over deformation temperature range of 1023-1423 K and strain rate ...The hot deformation behavior of Ti-6 Al-4 V-0.1 Ru titanium alloy was investigated by isothermal compression tests on a Gleeble-3500 thermal simulator over deformation temperature range of 1023-1423 K and strain rate of 0.01-10 s-1.Arrhenius-type constitutive models were developed for temperature ranges of bothα+βdual phase andβsingle phase at strain of 0.1.Afterwards,a series of material constants(including activation energy Q,material constants n,αand ln A)as polynomial functions of strain were introduced into Arrhenius-type models.Finally,the improved Arrhenius-type models in temperature field ofα+βandβphase were constructed.The results show that the improved Arrhenius-type models contribute to the calculation of Zener-Hollomon(Z)parameter,and the microstructural evolution mechanism is uncovered by combining microstructure observations with Z-parameter.Furthermore,the improved Arrhenius-type models are also helpful to improve the accuracy of finite element method(FEM)simulation in the deformation process of Ti-6 Al-4 V-0.1 Ru titanium alloy.展开更多
The deformation behavior of V-10Cr-5Ti alloy was studied on the Gleeble-1500 thermomechanical simulator at the temperatures of 950-1350℃, and the strain rates of 0.01-10 s^-1. Based on the Arrhenius model, dislocatio...The deformation behavior of V-10Cr-5Ti alloy was studied on the Gleeble-1500 thermomechanical simulator at the temperatures of 950-1350℃, and the strain rates of 0.01-10 s^-1. Based on the Arrhenius model, dislocation density model, nucleation model and grain growth model, a numerical cellular automaton (CA) model coupling simulation of hot deformation is established to simulate and characterize the microstructural evolution during DRX. The results show that the flow stress is fairly sensitive to the strain rate and deformation temperature. The error between the predicted stress by the Arrhenius model and the actual measured value is less than 8%. The initial average grain size calculated by the CA model is 86.25 μm, which is close to the experimental result (85.63 μm). The simulations show that the effect of initial grain size on the dynamic recrystallization microstructure evolution is not significant, while increasing the strain rate or reducing the temperature can refine the recrystallized grains.展开更多
A general CFD-PBE(computational fluid dynamics-population balance equation) solver for gas–liquid poly-dispersed flows of both low and high gas volume fractions is developed in OpenFOAM(open-source field operation an...A general CFD-PBE(computational fluid dynamics-population balance equation) solver for gas–liquid poly-dispersed flows of both low and high gas volume fractions is developed in OpenFOAM(open-source field operation and manipulation) in this work. Implementation of this solver in OpenFOAM is illustrated in detail. The PBE is solved with the cell average technique. The coupling between pressure and velocity is dealt with the transient PIMPLE algorithm, which is a merged PISO-SIMPLE(pressure implicit split operator-semi-implicit method for pressure-linked equations) algorithm. Results show generally good agreement with the published experimental data, whereas the modeling precision could be improved further with more sophisticated closure models for interfacial forces, the models for the bubble-induced turbulence and those for bubble coalescence and breakage.The results also indicate that the PBE could be solved out the PIMPLE loop to save much computation time while still preserving the time information on variables. This is important for CFD-PBE modeling of many actual gas–liquid problems, which are commonly high-turbulent flows with intrinsic transient and 3 D characteristics.展开更多
Influences of tip radius and sampling interval on applying atomic force microscopy (AFM) in quantitative surface evaluations are investigated by numerical simulations and experiments. Several evaluation parameters o...Influences of tip radius and sampling interval on applying atomic force microscopy (AFM) in quantitative surface evaluations are investigated by numerical simulations and experiments. Several evaluation parameters of surfaces ranging from amplitude to functional parameters are studied. Numerical and experimental results are in good agreements. The accuracy of estimating tip radius on random rough surface with Ganssian distribution of heights using a blind reconstruction method is also discussed theoretically. It is found that the accuracy is greatly depending on the ratio of actual tip radius to root-mean-square (rms) radius of curvature. To obtain an accurate estimation of tip radius under Gaussian rough surface, the ratio has to be larger than 3/2.展开更多
The Hadley circulation is one of the most important atmospheric circulations.Widening of the Hadley circulation has drawn extensive studies in the past decade.The key concern is that widening of the Hadley circulation...The Hadley circulation is one of the most important atmospheric circulations.Widening of the Hadley circulation has drawn extensive studies in the past decade.The key concern is that widening of the Hadley circulation would cause poleward shift of the subtropical dry zone.Various metrics have been applied to measure the widening of the tropics.What are responsible for the observed widening trends of the Hadley circulation? How anthropogenic and natural forcings caused the widening? How the widening results in regional climatic effects? These are the major questions in studing the widening of the Hadley circulation.While both observations and simulations all show widening of the Hadley circulation in the past few decades,there are no general agreements of changes in the strength of the Hadley circulation.Although some reanalysis datasets show strengthening of the Hadley circulation,it was shown that the strengthening trend could be artificial,and simulations show weakening of the Hadley circulation for global greenhouse warming.In the present paper,we shall briefly review the major progresses of studies in trends in width and strength of the Hadley circulation.We address answers to these questions,clarify inconsistent results,and propose ideas for future studies.展开更多
Riemann problems for the compressible Euler system in two space dimensions are complicated and difficult,but a viable alternative remains missing.The author lists merits of one-dimensional Riemann problems and compare...Riemann problems for the compressible Euler system in two space dimensions are complicated and difficult,but a viable alternative remains missing.The author lists merits of one-dimensional Riemann problems and compares them with those for the current two-dimensional Riemann problems,to illustrate their worthiness.Two-dimensional Riemann problems are approached via the methodology promoted by Andy Majda in the spirits of modern applied mathematics;that is,simplified model is built via asymptotic analysis,numerical simulation and theoretical analysis.A simplified model called the pressure gradient system is derived from the full Euler system via an asymptotic process.State-of-the-art numerical methods in numerical simulations are used to discern small-scale structures of the solutions,e.g.,semi-hyperbolic patches.Analytical methods are used to establish the validity of the structure revealed in the numerical simulation.The entire process,used in many of Majda's programs,is shown here for the two-dimensional Riemann problems for the compressible Euler systems of conservation laws.展开更多
The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses...The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses.At first,the strong governing formulation of these microstructures is briefly constructed,and the relevant weak formulation is then deduced based on the virtual work theorem.Afterwards,a subdividing approach,which cuts the elements intercepted by the interfaces between distinct phases automatically,is established with the aid of the level set method(LSM),and the discrete counterpart of the governing formula is obtained by combining the weak formulation derived and a discretized model.To be noted,two silent merits are found when the elaborated strategy is applied:(1) the continents constituting the microstructures of different scales can be arbitrarily-shaped and conveniently reproduced;(2) the periodic boundary condition commonly employed can be enforced on the external surfaces of representative unit cells(RUCs) with no difficulty.Besides,a boundary value problem(BVP) involving a simplified HTE nanostructure is designed,analytically solved,and hereafter applied to verify the correctness of the proposed strategy.It is observed that both the displacement and stress predictions by the computational approach are in good agreement with the relevant analytical results irrespective of the material combinations applied.Eventually,discussions are made on the influence of material organizations of both the 2D and 3D HTE microstructures at the ultrastructural and repeated rod levels,and some concluding remarks are drawn.展开更多
In this paper, a time-delayed predator-prey system is considered. The existence of Hopf bifurcations at the positive equilibrium is established by analyzing the distribution of the characteristic values. An explicit a...In this paper, a time-delayed predator-prey system is considered. The existence of Hopf bifurcations at the positive equilibrium is established by analyzing the distribution of the characteristic values. An explicit algorithm for determining the direction of the Hopf bifurcation and the stability of the bifurcating periodic solutions are derived by using the normal form and the center manifold theory. Numerical simulations to support the analytical conclusions are carried out.展开更多
We develop and assess a model of the turbulent burning velocity ST over a wide range of conditions.The aim is to obtain an explicit ST model for turbulent combustion modeling and flame analysis.The model consists of s...We develop and assess a model of the turbulent burning velocity ST over a wide range of conditions.The aim is to obtain an explicit ST model for turbulent combustion modeling and flame analysis.The model consists of sub models of the stretch factor and the turbulent flame area.The stretch factor characterizes the flame response of turbulence stretch and incorporates detailed chemistry and transport effects with a lookup table of laminar counterflow flames.The flame area model captures the area growth based on Lagrangian statistics of propagating surfaces and considers the effects of turbulence length scales and fuel characteristics.The present model predicts sT via an algebraic expression without free parameters.We assess the model using 490 cases of the direct numerical simulation or experiment reported from various research groups on planar and Bunsen flames over a wide range of conditions,covering fuels from hydrogen to n-dodecane,pressures from 1 to 30 atm,lean and rich mixtures,turbulence intensity ratios from 0.1 to 177.6,and turbulence length ratios from 0.5 to 66.7.Despite the scattering sT data in the literature,the comprehensive comparison shows that the proposed ST model has an overall good agreement over the wide range of conditions,with the averaged modeling error of 28.1%.展开更多
基金Projects(51774125,51674112)supported by the National Natural Science Foundation of ChinaProject(2018YFB1107905)supported by the National Key Research and Development Program of China。
文摘The experiments of laser fusion welding with Al-foil addition was carried out for DP590 dual-phase steel and AZ31B magnesium alloy in an overlap steel-on-magnesium configuration.Temperature field was simulated by COMSOL finite element software for steel/magnesium laser fusion welding.The results show that when Al-foil is added,some defects,such as pores,cracks and softening in heat affected zone(HAZ),can be avoided in welding joint,the bonding strength of steel/magnesium joints is increased,heat transfer between steel and magnesium is regulated.In the case of adding Al-foil,welding pool is divided into two parts,the upper and lower pools contact each other but do not mix,the transition layer at the interface between the upper and lower molten pools mainly contains Al−Fe phases,such as AlFe,Al2Fe and AlFe3,and these new phases are helpful for promoting the metallurgical connection between the upper and lower molten pools.Hence,adding Al-foil laser fusion welding is an effective way in joining steel to magnesium alloy.
基金financially supported by the National Natural Science Foundation of China (Nos.52071042,51771038)the Chongqing Talent Plan,China (No.CQYC202003047)Chongqing Natural Science Foundation,China (Nos.cstc2018jcyj AX0249,cstc2018jcyj AX0653)。
文摘A new severe plastic deformation method for manufacturing tubes made of AZ31 magnesium alloy with a large diameter was developed,which is called the TCESE(tube continuous extrusion−shear−expanding)process.The process combines direct extrusion with a two-step shear−expanding process.The influences of expanding ratios,extrusion temperatures on the deformation of finite element meshes,strain evolution and flow velocity of tube blanks during the TCESE process were researched based on numerical simulations by using DEFORM-3D software.Simulation results show that the maximum expanding ratio is 3.0 in the TCESE process.The deformation of finite element meshes of tube blanks is inhomogeneous in the shear−expanding zone,and the equivalent strains increase significantly during the TCESE process of the AZ31 magnesium alloy.A extrusion temperature of 380°C and expanding ratio of 2.0 were selected as the optimized process parameters from the numerical simulation results.The average grain size of tubes fabricated by the TCESE process is approximately 10μm.The TCESE process can refine grains of magnesium alloy tubes with the occurrence of dynamic recrystallization.The(0001)basal texture intensities of the magnesium alloy tube blanks decrease due to continuous plastic deformation during the TCESE process.The average hardness of the extruded tubes is approximately HV 75,which is obviously improved.
文摘With the microscopic phase-field dynamic model, the effects of temperature and concentration on the nucleation incubation time of Ni75AlxV25-x alloy were studied and the relation between the incubation time and precipitation mechanism was investigated by using the atomic occupation probability picture and average order parameter curve. The simulation results demonstrate that there exists the incubation time for different precipitation mechanisms~ such as non-classical nucleation, the mixed style of non-classical nucleation and spinodal decomposition, and spinodal ordering; and the incubation time shortens in turn for the three kinds of mechanisms. With the increase of Al content of Ni75AlxV25-x alloy, the incubation time of Llz phases shortens continuously and that of DOzz phases is prolonged. The effects of temperature on the incubation time of Llz and DOzz phases are accordant, i.e. the incuba- tion time is greatly prolonged with the temperature rising.
基金Projected(51775068)supported by the National Natural Science Foundation of China.
文摘The hot deformation behavior of Ti-6 Al-4 V-0.1 Ru titanium alloy was investigated by isothermal compression tests on a Gleeble-3500 thermal simulator over deformation temperature range of 1023-1423 K and strain rate of 0.01-10 s-1.Arrhenius-type constitutive models were developed for temperature ranges of bothα+βdual phase andβsingle phase at strain of 0.1.Afterwards,a series of material constants(including activation energy Q,material constants n,αand ln A)as polynomial functions of strain were introduced into Arrhenius-type models.Finally,the improved Arrhenius-type models in temperature field ofα+βandβphase were constructed.The results show that the improved Arrhenius-type models contribute to the calculation of Zener-Hollomon(Z)parameter,and the microstructural evolution mechanism is uncovered by combining microstructure observations with Z-parameter.Furthermore,the improved Arrhenius-type models are also helpful to improve the accuracy of finite element method(FEM)simulation in the deformation process of Ti-6 Al-4 V-0.1 Ru titanium alloy.
基金Project(51405110)supported by the National Natural Science Foundation of ChinaProject(20132302120002)supported by the Specialized Research Fund for the Doctoral Program of Higher Education,ChinaProject(LBH-Z14096)supported by Heilongjiang Province Postdoctoral Fund,China
文摘The deformation behavior of V-10Cr-5Ti alloy was studied on the Gleeble-1500 thermomechanical simulator at the temperatures of 950-1350℃, and the strain rates of 0.01-10 s^-1. Based on the Arrhenius model, dislocation density model, nucleation model and grain growth model, a numerical cellular automaton (CA) model coupling simulation of hot deformation is established to simulate and characterize the microstructural evolution during DRX. The results show that the flow stress is fairly sensitive to the strain rate and deformation temperature. The error between the predicted stress by the Arrhenius model and the actual measured value is less than 8%. The initial average grain size calculated by the CA model is 86.25 μm, which is close to the experimental result (85.63 μm). The simulations show that the effect of initial grain size on the dynamic recrystallization microstructure evolution is not significant, while increasing the strain rate or reducing the temperature can refine the recrystallized grains.
基金Supported by the National Key Research and Development Program(2016YFB0301702)National Natural Science Foundation of China(21776284,21476236)+1 种基金Key Research Program of Frontier Sciences,CAS(QYZDJ-SSW-JSC030)Jiangsu National Synergetic Innovation Center for Advanced Materials
文摘A general CFD-PBE(computational fluid dynamics-population balance equation) solver for gas–liquid poly-dispersed flows of both low and high gas volume fractions is developed in OpenFOAM(open-source field operation and manipulation) in this work. Implementation of this solver in OpenFOAM is illustrated in detail. The PBE is solved with the cell average technique. The coupling between pressure and velocity is dealt with the transient PIMPLE algorithm, which is a merged PISO-SIMPLE(pressure implicit split operator-semi-implicit method for pressure-linked equations) algorithm. Results show generally good agreement with the published experimental data, whereas the modeling precision could be improved further with more sophisticated closure models for interfacial forces, the models for the bubble-induced turbulence and those for bubble coalescence and breakage.The results also indicate that the PBE could be solved out the PIMPLE loop to save much computation time while still preserving the time information on variables. This is important for CFD-PBE modeling of many actual gas–liquid problems, which are commonly high-turbulent flows with intrinsic transient and 3 D characteristics.
文摘Influences of tip radius and sampling interval on applying atomic force microscopy (AFM) in quantitative surface evaluations are investigated by numerical simulations and experiments. Several evaluation parameters of surfaces ranging from amplitude to functional parameters are studied. Numerical and experimental results are in good agreements. The accuracy of estimating tip radius on random rough surface with Ganssian distribution of heights using a blind reconstruction method is also discussed theoretically. It is found that the accuracy is greatly depending on the ratio of actual tip radius to root-mean-square (rms) radius of curvature. To obtain an accurate estimation of tip radius under Gaussian rough surface, the ratio has to be larger than 3/2.
基金supported by the National Natural Science Foundation of China(41530423,and 41761144072)
文摘The Hadley circulation is one of the most important atmospheric circulations.Widening of the Hadley circulation has drawn extensive studies in the past decade.The key concern is that widening of the Hadley circulation would cause poleward shift of the subtropical dry zone.Various metrics have been applied to measure the widening of the tropics.What are responsible for the observed widening trends of the Hadley circulation? How anthropogenic and natural forcings caused the widening? How the widening results in regional climatic effects? These are the major questions in studing the widening of the Hadley circulation.While both observations and simulations all show widening of the Hadley circulation in the past few decades,there are no general agreements of changes in the strength of the Hadley circulation.Although some reanalysis datasets show strengthening of the Hadley circulation,it was shown that the strengthening trend could be artificial,and simulations show weakening of the Hadley circulation for global greenhouse warming.In the present paper,we shall briefly review the major progresses of studies in trends in width and strength of the Hadley circulation.We address answers to these questions,clarify inconsistent results,and propose ideas for future studies.
基金supported partially by the National Science Foundation (No.DMS-0603859)
文摘Riemann problems for the compressible Euler system in two space dimensions are complicated and difficult,but a viable alternative remains missing.The author lists merits of one-dimensional Riemann problems and compares them with those for the current two-dimensional Riemann problems,to illustrate their worthiness.Two-dimensional Riemann problems are approached via the methodology promoted by Andy Majda in the spirits of modern applied mathematics;that is,simplified model is built via asymptotic analysis,numerical simulation and theoretical analysis.A simplified model called the pressure gradient system is derived from the full Euler system via an asymptotic process.State-of-the-art numerical methods in numerical simulations are used to discern small-scale structures of the solutions,e.g.,semi-hyperbolic patches.Analytical methods are used to establish the validity of the structure revealed in the numerical simulation.The entire process,used in many of Majda's programs,is shown here for the two-dimensional Riemann problems for the compressible Euler systems of conservation laws.
基金supported by the National Natural Science Foundation of China(Grant Nos.51535010,51305362,11372260&11572266)the Fundamental Research Funds for the Central Universities(Grant Nos.2682014BR016&2682016CX024)the China Scholar Council
文摘The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses.At first,the strong governing formulation of these microstructures is briefly constructed,and the relevant weak formulation is then deduced based on the virtual work theorem.Afterwards,a subdividing approach,which cuts the elements intercepted by the interfaces between distinct phases automatically,is established with the aid of the level set method(LSM),and the discrete counterpart of the governing formula is obtained by combining the weak formulation derived and a discretized model.To be noted,two silent merits are found when the elaborated strategy is applied:(1) the continents constituting the microstructures of different scales can be arbitrarily-shaped and conveniently reproduced;(2) the periodic boundary condition commonly employed can be enforced on the external surfaces of representative unit cells(RUCs) with no difficulty.Besides,a boundary value problem(BVP) involving a simplified HTE nanostructure is designed,analytically solved,and hereafter applied to verify the correctness of the proposed strategy.It is observed that both the displacement and stress predictions by the computational approach are in good agreement with the relevant analytical results irrespective of the material combinations applied.Eventually,discussions are made on the influence of material organizations of both the 2D and 3D HTE microstructures at the ultrastructural and repeated rod levels,and some concluding remarks are drawn.
文摘In this paper, a time-delayed predator-prey system is considered. The existence of Hopf bifurcations at the positive equilibrium is established by analyzing the distribution of the characteristic values. An explicit algorithm for determining the direction of the Hopf bifurcation and the stability of the bifurcating periodic solutions are derived by using the normal form and the center manifold theory. Numerical simulations to support the analytical conclusions are carried out.
基金supported by the National Natural Science Foundation of China(Grant Nos.91841302,11925201,and 11988102)the National Key Research and Development.Program of China(Grant No.2020YFE0204200)。
文摘We develop and assess a model of the turbulent burning velocity ST over a wide range of conditions.The aim is to obtain an explicit ST model for turbulent combustion modeling and flame analysis.The model consists of sub models of the stretch factor and the turbulent flame area.The stretch factor characterizes the flame response of turbulence stretch and incorporates detailed chemistry and transport effects with a lookup table of laminar counterflow flames.The flame area model captures the area growth based on Lagrangian statistics of propagating surfaces and considers the effects of turbulence length scales and fuel characteristics.The present model predicts sT via an algebraic expression without free parameters.We assess the model using 490 cases of the direct numerical simulation or experiment reported from various research groups on planar and Bunsen flames over a wide range of conditions,covering fuels from hydrogen to n-dodecane,pressures from 1 to 30 atm,lean and rich mixtures,turbulence intensity ratios from 0.1 to 177.6,and turbulence length ratios from 0.5 to 66.7.Despite the scattering sT data in the literature,the comprehensive comparison shows that the proposed ST model has an overall good agreement over the wide range of conditions,with the averaged modeling error of 28.1%.