The smelting reduction process of the ilmenite in an electric arc furnace(EAF)is a commonly used technology for producing titanium slag in the world.It has particular significance to analyze the velocity-temperature-e...The smelting reduction process of the ilmenite in an electric arc furnace(EAF)is a commonly used technology for producing titanium slag in the world.It has particular significance to analyze the velocity-temperature-electromagnetics multi-physical field in an EAF for improving its productivity and reducing energy consumption.A transient three-dimensional mathematical model was developed to characterize the flow,heat transfer,and electromagnetic behavior in a titanium slag EAF.For describing the electromagnetic field and its effects on velocity and temperature distribution in the furnace,magnetohydrodynamic equations and conservation equations for mass,momentum,and energy were solved simultaneously by compiling the user-defined function program.The numerical model was verified by comparing with the literature data.The results indicate that the Lorentz force is the main driving force of the velocity and temperature distribution.Moreover,the influence of input current and location of electrodes on the multi-physical field distribution was also investigated.It is found that the appropriate range of input current and diameter of pitch circle are about 30,000 A and 3000-3500 mm,respectively.The mathematical model established can characterize the multi-physical field more accu-rately than before,which can provide valuable guidance for the operation improvement and design optimization of the EAF for producing titanium slag.展开更多
Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
Gassy soils are distributed in relatively shallow layers the Quaternary deposit in Hangzhou Bay area. The shallow gassy soils significantly affect the construction of underground projects. Proper characterization of s...Gassy soils are distributed in relatively shallow layers the Quaternary deposit in Hangzhou Bay area. The shallow gassy soils significantly affect the construction of underground projects. Proper characterization of spatial distribution of shallow gassy soils is indispensable prior to construction of underground projects in the area. Due to the costly conditions required in the site investigation for gassy soils, only a limited number of gas pressure data can be obtained in engineering practice, which leads to the uncertainty in characterizing spatial distribution of gassy soils. Determining the number of boreholes for investigating gassy soils and their corresponding locations is pivotal to reducing construction risk induced by gassy soils. However, this primarily relies on the engineering experience in the current site investigation practice. This study develops a probabilistic site investigation optimization method for planning investigation schemes (including the number and locations of boreholes) of gassy soils based on the conditional random field and Monte Carlo simulation. The proposed method aims to provide an optimal investigation scheme before the site investigation based on prior knowledge. Finally, the proposed approach is illustrated using a case study.展开更多
With the evolution of nickel-based single crystal superalloys,there is an increase in heavy elements such as Re and Ru.This has made solutal convection more pronounced during the directional solidification process,lea...With the evolution of nickel-based single crystal superalloys,there is an increase in heavy elements such as Re and Ru.This has made solutal convection more pronounced during the directional solidification process,leading to solute redistribution and increasing the risk of casting defects such as low-angle grain boundaries.To avoid casting defects,downward directional solidification(DWS)method is adopted to eliminate solutal convection and change solute redistribution.However,there is currently no in-situ characterization or quantitative simulation studying the solute redistribution during DWS and upward directional solidification(UWS)processes.A multicomponent phase field simulation coupled with lattice Boltzmann method was employed to quantitatively investigate changes in dendrite morphology,solutal convection and deviation of dendrite tips from the perspective of solute redistribution during UWS and DWS processes.The simulation of microstructure agrees well with the experimental results.The mechanism that explains how solutal convection affects side branching behavior is depicted.A novel approach is introduced to characterize dendrite deviation,elucidating the reasons why defects are prone to occur under the influence of natural convection and solute redistribution.展开更多
A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radia...A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]展开更多
This article deals with generation and application of three-dimensional (3D) atmospheric turbulence field in large aircraft real-time flight simulation. The modeling requirements for the turbulence field of large airc...This article deals with generation and application of three-dimensional (3D) atmospheric turbulence field in large aircraft real-time flight simulation. The modeling requirements for the turbulence field of large aircraft flight simulation are analyzed here. The spatial turbulence field is generated in the frequency domain by using the Monte Carlo method,and then transformed back to the time domain with the 3D inverse Fourier transform. The von Karman model is adopted for an accurate description of the turb...展开更多
This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle mo...This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle modulation of turbulence. The particles mainly excite the initial instability of the jet and bring about the earlier breakup of vortex rings in the near-field. The flow fluc- tuating intensity either in the axial or in the radial directions is hence increased by particles. The article also describes the mean velocity modulated by particles. The changing statistical velocity induced by particle modulation implies the effects of modulation of the local flow structures. This study is expected to be useful to the control of two-phase turbulent jets.展开更多
Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the...Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.展开更多
Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short ti...Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given.展开更多
To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS proce...To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.展开更多
In this paper, integrating the Yangtze Estuary with the Hangzhou Bay, a 2-D velocity field model is established. In the model, fine self-adaptive grids are employed to adapt to the complicated coastal shape. The hydro...In this paper, integrating the Yangtze Estuary with the Hangzhou Bay, a 2-D velocity field model is established. In the model, fine self-adaptive grids are employed to adapt to the complicated coastal shape. The hydrodynamic equations satisfied by two contravariant components of velocity vector and surface elevation in non-orthogonal curvilinear coordinates are used. In each momentum equation the coefficients before the two partial derivatives of surface elevation with respect to variables of alternative direction coordinates have different orders of magnitude, i. e., the derivative with the larger coefficient may play a more important role than that with the smaller one. With this advantage, the ADI scheme can be easily employed. The hydrodynamic factors include tidal current, river runoff and wind-induced current. In terms of tidal current, seven main constituents in the area are considered in the open boundaries. The verifications of surface elevation process and current velocity process in the spring tide and in the neap tide show that the model can preferably reflect current fields in the area. Through the simulation of Lagrangian residual current fields in summer and in winter, the paths of the exchange of water and sediment between the Yangtze Estuary and the Hangzhou Bay are elementarily discussed.展开更多
The temperature field of stainless steel during twin-roll strip casting was simulated by experiment and a finite element （FE） model. By comparing the measured result with the simulated values, it is found that they ...The temperature field of stainless steel during twin-roll strip casting was simulated by experiment and a finite element （FE） model. By comparing the measured result with the simulated values, it is found that they fit close to each other, which indicates this FE model is effective. Based on this model, the effects of roll gap （t） and roll radius （R） on solidification were simulated. The simulated results give the relationship between t or R and the position of the freezing point. The larger the t is and the smaller the R is, the closer the position of the freezing point is to the exit.展开更多
Study on tectonic fractures based on the inversion of tectonic stress fields is an effective method. In this study, a geological model was set up based on geological data from the Hudi Coal Mine, Qinshui Basin, a mech...Study on tectonic fractures based on the inversion of tectonic stress fields is an effective method. In this study, a geological model was set up based on geological data from the Hudi Coal Mine, Qinshui Basin, a mechanical model was established under the condition of rock mechanics and geostress, and the finite element method was used to simulate the paleotectonic stress field. Based on the Griffith and Mohr-Coulomb criterion, the distribution of tectonic fractures in the Shanxi Formation during the Indosinian, Yanshanian, and Himalayan period can be predicted with the index of comprehensive rupture rate. The results show that the acting force of the Pacific Plate and the India Plate to the North China Plate formed the direction of principal stress is N-S, NW - SE, and NE - SW, respectively, in different periods in the study area. Changes in the direction and strength of the acting force led to the regional gradients of tectonic stress magnitude, which resulted in an asymmetrical distribution state of the stress conditions in different periods. It is suggested that the low-stress areas are mainly located in the fault zones and extend along the direction of the fault zones. Furthermore, the high-stress areas are located in the junction of fold belts and the binding site of multiple folds. The development of tectonic fractures was affected by the distribution of stress intensity and the tectonic position of folds and faults, which resulted in some developed areas with level I and II. There are obvious differences in the development of tectonic fractures in the fold and fault zones and the anticline and syncline structure at the same fold zones. The tectonic fractures of the Shanxi Formation during the Himalayan period are more developed than those during the Indosinian and Yanshanian period due to the superposition of the late tectonic movement to the early tectonic movement and the differences in the magnitude and direction of stress intensity.展开更多
Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperatu...Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperature variation along the pipe, was proposed for simulating the temperature field of early-age concrete structures containing cooling pipes. The improved model was verified with an engineering example. Then, the p-version self-adaption algorithm for the improved embedded model was deduced, and the initial values and boundary conditions were examined. Comparison of some numerical samples shows that the proposed model can provide satisfying precision and a higher efficiency. The analysis efficiency can be doubled at the same precision, even for a large-scale element. The p-version algorithm can fit grids of different sizes for the temperature field simulation. The convenience of the proposed algorithm lies in the possibility of locating more pipe segments in one element without the need of so regular a shape as in the explicit model.展开更多
Mold-filling process of thin-walled castings under the condition of traveling magnetic field has been studied by physical simulation method using gallium melt and fast speed photography. Flow morphology and its format...Mold-filling process of thin-walled castings under the condition of traveling magnetic field has been studied by physical simulation method using gallium melt and fast speed photography. Flow morphology and its formation mechanism were obtained and discussed for thin-walled casting. The influences of magnetic field density on the filling ability, filling velocity and mold filling time have been studied. The differences in filling capability between gravity casting and casting under the traveling magnetic field have been compared. The results indicate that the mold filling ability of the gallium melt increases greatly under the condition of traveling magnetic field; the filling time is shortened from 18 s under gravity field to 3 s under the traveling magnetic field and average flow rate of the melt increases from 1.6 to 8.68 cm3/s; the change law of the cross-section morphology of the gallium melt during the mold filling is that at first, the cross-section area does not change, then it decreases gradually. When the front of the melt reaches the end of the mold cavity, the front melt will backfill the mold; the wider the width of mold cavity, the better the mold filling ability. The mold filling ability of gallium melt in mold with upper magnetic conductor is better than that without upper magnetic conductor.展开更多
The thermal model of FSW based on the thermal elastic-plastic finite element method, and the transient temperature distribution of FS welded joints of 2024-T4 Al was simulated by using this model, which provides usefu...The thermal model of FSW based on the thermal elastic-plastic finite element method, and the transient temperature distribution of FS welded joints of 2024-T4 Al was simulated by using this model, which provides useful information for the investigation of FSW process. Simulation resuits show that the temperature distribution of the weld gradually decreases toward periphery in a radiate format, whose center is the probe, and the highest temperature in the weld can reach about 400℃. The initial terminal of the weld is a zone, in which the temperature gradient is great, and defects of the welding are easily produced in this zone. Temperature change at the end of the welded joint is as layer variation, the local serious defects are not easy to produce in this zone.展开更多
Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distrib...Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distribution of a mold wall, a simulated apparatus of mold was designed and experiments were performed by it. The measured results indicated that the mold wall temperature approaches the temperature of cooling-water. An equivalent thermal-conductivity coefficient was proposed and deduced on the basis of the conclusion of the experiments. This coefficient was applied to solve the heat transfer between the melt and cooling water, and to characterize the heat transfer capacity of the mold. By this equivalent thermal-conductivity coefficient, it is very easy and convenient to numerically simulate the solidification process of continuous casting. And the calculation results are in agreement with the experiments. The effects of custing speed and water flow rate on the mold temperature field were also discussed.展开更多
Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovati...Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovative form of the POD based on cross power spectral density matrices. By introducing a discretizing scheme, the CPOD-based spectral representation method is obtained for use in stochastic simulation. Moreover, some criteria are proposed that allow the truncation order of CPOD to be conveniently determined. A numerical example to illustrate the application of the proposed method for the simulation of a wind velocity field is provided.展开更多
The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickn...The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickness on the dendrite growth were studied. Results indicate that the grain grows into an equiaxial dendrite during free flow and into a typical branched structure under forced flow. The radius of an initial crystal can affect the growth of side-branches but not the stability of the dendrite＇s tip when an appropriate value is assigned to it. With an increase in space steps, side-branches appear at the upstream of the longitudinal principal branch and they grow rapidly. With an increase in the interface thickness, the trunk of the longitudinal upstream and lateral principal branches grow longer and become more slender while the number of secondary branches increases.展开更多
Based on the Cahn-Hilliard phase field model, a three-dimensional multiple-field coupling model for simulating the motion characteristics of a rising bubble in a liquid is established in a gas-liquid two-phase flow. T...Based on the Cahn-Hilliard phase field model, a three-dimensional multiple-field coupling model for simulating the motion characteristics of a rising bubble in a liquid is established in a gas-liquid two-phase flow. The gas-liquid interface motion is simulated by using a phase-field method, and the effect of the electric field intensity on bubble dynamics is studied without electric field, or with vertical electric field or horizontal electric field. Through the coupling effect of electric field and flow field, the deformation of a single rising bubble and the formation of wake vortices under the action of gravity and electric field force are studied in detail. The correctness of the results is verified by mass conservation, and the influences of different electric field directions and different voltages on the movement of bubbles in liquid are considered. The results show that the ratio of the length to axis is proportional to the strength of the electric field when the air bubble is stretched into an ellipsoid along the electric field line under the action of electrostatic gravity and surface tension. In addition, the bubble rising speed is affected by the electric field, the vertical electric field accelerates the bubble rise, and the horizontal direction slows it down.展开更多
基金supported by the National Natural Science Foundation of China(No.U2003215).
文摘The smelting reduction process of the ilmenite in an electric arc furnace(EAF)is a commonly used technology for producing titanium slag in the world.It has particular significance to analyze the velocity-temperature-electromagnetics multi-physical field in an EAF for improving its productivity and reducing energy consumption.A transient three-dimensional mathematical model was developed to characterize the flow,heat transfer,and electromagnetic behavior in a titanium slag EAF.For describing the electromagnetic field and its effects on velocity and temperature distribution in the furnace,magnetohydrodynamic equations and conservation equations for mass,momentum,and energy were solved simultaneously by compiling the user-defined function program.The numerical model was verified by comparing with the literature data.The results indicate that the Lorentz force is the main driving force of the velocity and temperature distribution.Moreover,the influence of input current and location of electrodes on the multi-physical field distribution was also investigated.It is found that the appropriate range of input current and diameter of pitch circle are about 30,000 A and 3000-3500 mm,respectively.The mathematical model established can characterize the multi-physical field more accu-rately than before,which can provide valuable guidance for the operation improvement and design optimization of the EAF for producing titanium slag.
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
文摘Gassy soils are distributed in relatively shallow layers the Quaternary deposit in Hangzhou Bay area. The shallow gassy soils significantly affect the construction of underground projects. Proper characterization of spatial distribution of shallow gassy soils is indispensable prior to construction of underground projects in the area. Due to the costly conditions required in the site investigation for gassy soils, only a limited number of gas pressure data can be obtained in engineering practice, which leads to the uncertainty in characterizing spatial distribution of gassy soils. Determining the number of boreholes for investigating gassy soils and their corresponding locations is pivotal to reducing construction risk induced by gassy soils. However, this primarily relies on the engineering experience in the current site investigation practice. This study develops a probabilistic site investigation optimization method for planning investigation schemes (including the number and locations of boreholes) of gassy soils based on the conditional random field and Monte Carlo simulation. The proposed method aims to provide an optimal investigation scheme before the site investigation based on prior knowledge. Finally, the proposed approach is illustrated using a case study.
基金supported by the stable support project and the Major National Science and Technology Project(2017-VII-0008-0101).
文摘With the evolution of nickel-based single crystal superalloys,there is an increase in heavy elements such as Re and Ru.This has made solutal convection more pronounced during the directional solidification process,leading to solute redistribution and increasing the risk of casting defects such as low-angle grain boundaries.To avoid casting defects,downward directional solidification(DWS)method is adopted to eliminate solutal convection and change solute redistribution.However,there is currently no in-situ characterization or quantitative simulation studying the solute redistribution during DWS and upward directional solidification(UWS)processes.A multicomponent phase field simulation coupled with lattice Boltzmann method was employed to quantitatively investigate changes in dendrite morphology,solutal convection and deviation of dendrite tips from the perspective of solute redistribution during UWS and DWS processes.The simulation of microstructure agrees well with the experimental results.The mechanism that explains how solutal convection affects side branching behavior is depicted.A novel approach is introduced to characterize dendrite deviation,elucidating the reasons why defects are prone to occur under the influence of natural convection and solute redistribution.
文摘A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]
文摘This article deals with generation and application of three-dimensional (3D) atmospheric turbulence field in large aircraft real-time flight simulation. The modeling requirements for the turbulence field of large aircraft flight simulation are analyzed here. The spatial turbulence field is generated in the frequency domain by using the Monte Carlo method,and then transformed back to the time domain with the 3D inverse Fourier transform. The von Karman model is adopted for an accurate description of the turb...
基金National Natural Science Foundation of China (50706021)Ph.D.Programs Foundation of Ministry of Education of China (20070003018)TNList Cross-discipline Foundation
文摘This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle modulation of turbulence. The particles mainly excite the initial instability of the jet and bring about the earlier breakup of vortex rings in the near-field. The flow fluc- tuating intensity either in the axial or in the radial directions is hence increased by particles. The article also describes the mean velocity modulated by particles. The changing statistical velocity induced by particle modulation implies the effects of modulation of the local flow structures. This study is expected to be useful to the control of two-phase turbulent jets.
基金the National High Technical Reasearch and Development Programme of China （No. 2003AA327140） the National Natural Science Foundation of China （No. 50374081）.
文摘Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.
文摘Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given.
文摘To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.
基金National Natural Science Foundation of China under contract No.49776279National Excellent Youth Foundation of China under contract No.49825161
文摘In this paper, integrating the Yangtze Estuary with the Hangzhou Bay, a 2-D velocity field model is established. In the model, fine self-adaptive grids are employed to adapt to the complicated coastal shape. The hydrodynamic equations satisfied by two contravariant components of velocity vector and surface elevation in non-orthogonal curvilinear coordinates are used. In each momentum equation the coefficients before the two partial derivatives of surface elevation with respect to variables of alternative direction coordinates have different orders of magnitude, i. e., the derivative with the larger coefficient may play a more important role than that with the smaller one. With this advantage, the ADI scheme can be easily employed. The hydrodynamic factors include tidal current, river runoff and wind-induced current. In terms of tidal current, seven main constituents in the area are considered in the open boundaries. The verifications of surface elevation process and current velocity process in the spring tide and in the neap tide show that the model can preferably reflect current fields in the area. Through the simulation of Lagrangian residual current fields in summer and in winter, the paths of the exchange of water and sediment between the Yangtze Estuary and the Hangzhou Bay are elementarily discussed.
文摘The temperature field of stainless steel during twin-roll strip casting was simulated by experiment and a finite element （FE） model. By comparing the measured result with the simulated values, it is found that they fit close to each other, which indicates this FE model is effective. Based on this model, the effects of roll gap （t） and roll radius （R） on solidification were simulated. The simulated results give the relationship between t or R and the position of the freezing point. The larger the t is and the smaller the R is, the closer the position of the freezing point is to the exit.
基金financially supported by the National Natural Science Foundation of China (Nos.41402138, 41330635, and 41272154)Key Laboratory of Coal-based CO2 Capture and Geological Storage, Jiangsu Province (China University of Mining and Technology) (No.2016B04)
文摘Study on tectonic fractures based on the inversion of tectonic stress fields is an effective method. In this study, a geological model was set up based on geological data from the Hudi Coal Mine, Qinshui Basin, a mechanical model was established under the condition of rock mechanics and geostress, and the finite element method was used to simulate the paleotectonic stress field. Based on the Griffith and Mohr-Coulomb criterion, the distribution of tectonic fractures in the Shanxi Formation during the Indosinian, Yanshanian, and Himalayan period can be predicted with the index of comprehensive rupture rate. The results show that the acting force of the Pacific Plate and the India Plate to the North China Plate formed the direction of principal stress is N-S, NW - SE, and NE - SW, respectively, in different periods in the study area. Changes in the direction and strength of the acting force led to the regional gradients of tectonic stress magnitude, which resulted in an asymmetrical distribution state of the stress conditions in different periods. It is suggested that the low-stress areas are mainly located in the fault zones and extend along the direction of the fault zones. Furthermore, the high-stress areas are located in the junction of fold belts and the binding site of multiple folds. The development of tectonic fractures was affected by the distribution of stress intensity and the tectonic position of folds and faults, which resulted in some developed areas with level I and II. There are obvious differences in the development of tectonic fractures in the fold and fault zones and the anticline and syncline structure at the same fold zones. The tectonic fractures of the Shanxi Formation during the Himalayan period are more developed than those during the Indosinian and Yanshanian period due to the superposition of the late tectonic movement to the early tectonic movement and the differences in the magnitude and direction of stress intensity.
基金supported by the National Natural Science Foundation of China(Grant No.51109071)
文摘Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperature variation along the pipe, was proposed for simulating the temperature field of early-age concrete structures containing cooling pipes. The improved model was verified with an engineering example. Then, the p-version self-adaption algorithm for the improved embedded model was deduced, and the initial values and boundary conditions were examined. Comparison of some numerical samples shows that the proposed model can provide satisfying precision and a higher efficiency. The analysis efficiency can be doubled at the same precision, even for a large-scale element. The p-version algorithm can fit grids of different sizes for the temperature field simulation. The convenience of the proposed algorithm lies in the possibility of locating more pipe segments in one element without the need of so regular a shape as in the explicit model.
文摘Mold-filling process of thin-walled castings under the condition of traveling magnetic field has been studied by physical simulation method using gallium melt and fast speed photography. Flow morphology and its formation mechanism were obtained and discussed for thin-walled casting. The influences of magnetic field density on the filling ability, filling velocity and mold filling time have been studied. The differences in filling capability between gravity casting and casting under the traveling magnetic field have been compared. The results indicate that the mold filling ability of the gallium melt increases greatly under the condition of traveling magnetic field; the filling time is shortened from 18 s under gravity field to 3 s under the traveling magnetic field and average flow rate of the melt increases from 1.6 to 8.68 cm3/s; the change law of the cross-section morphology of the gallium melt during the mold filling is that at first, the cross-section area does not change, then it decreases gradually. When the front of the melt reaches the end of the mold cavity, the front melt will backfill the mold; the wider the width of mold cavity, the better the mold filling ability. The mold filling ability of gallium melt in mold with upper magnetic conductor is better than that without upper magnetic conductor.
文摘The thermal model of FSW based on the thermal elastic-plastic finite element method, and the transient temperature distribution of FS welded joints of 2024-T4 Al was simulated by using this model, which provides useful information for the investigation of FSW process. Simulation resuits show that the temperature distribution of the weld gradually decreases toward periphery in a radiate format, whose center is the probe, and the highest temperature in the weld can reach about 400℃. The initial terminal of the weld is a zone, in which the temperature gradient is great, and defects of the welding are easily produced in this zone. Temperature change at the end of the welded joint is as layer variation, the local serious defects are not easy to produce in this zone.
基金the National Natural Science Foundation of China （No. 599995442）.
文摘Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distribution of a mold wall, a simulated apparatus of mold was designed and experiments were performed by it. The measured results indicated that the mold wall temperature approaches the temperature of cooling-water. An equivalent thermal-conductivity coefficient was proposed and deduced on the basis of the conclusion of the experiments. This coefficient was applied to solve the heat transfer between the melt and cooling water, and to characterize the heat transfer capacity of the mold. By this equivalent thermal-conductivity coefficient, it is very easy and convenient to numerically simulate the solidification process of continuous casting. And the calculation results are in agreement with the experiments. The effects of custing speed and water flow rate on the mold temperature field were also discussed.
基金Research Fund for Communications in Western China Under Grant No. 200431800028
文摘Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovative form of the POD based on cross power spectral density matrices. By introducing a discretizing scheme, the CPOD-based spectral representation method is obtained for use in stochastic simulation. Moreover, some criteria are proposed that allow the truncation order of CPOD to be conveniently determined. A numerical example to illustrate the application of the proposed method for the simulation of a wind velocity field is provided.
基金supported by a Program of the Lanzhou Science and Technology Bureau (No.2009-1-9)the Doctoral Fund of Lanzhou University of Technology(No.SB01200606)
文摘The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickness on the dendrite growth were studied. Results indicate that the grain grows into an equiaxial dendrite during free flow and into a typical branched structure under forced flow. The radius of an initial crystal can affect the growth of side-branches but not the stability of the dendrite＇s tip when an appropriate value is assigned to it. With an increase in space steps, side-branches appear at the upstream of the longitudinal principal branch and they grow rapidly. With an increase in the interface thickness, the trunk of the longitudinal upstream and lateral principal branches grow longer and become more slender while the number of secondary branches increases.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51661020,11504149,and 11364024)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China(Grant No.J201304)
文摘Based on the Cahn-Hilliard phase field model, a three-dimensional multiple-field coupling model for simulating the motion characteristics of a rising bubble in a liquid is established in a gas-liquid two-phase flow. The gas-liquid interface motion is simulated by using a phase-field method, and the effect of the electric field intensity on bubble dynamics is studied without electric field, or with vertical electric field or horizontal electric field. Through the coupling effect of electric field and flow field, the deformation of a single rising bubble and the formation of wake vortices under the action of gravity and electric field force are studied in detail. The correctness of the results is verified by mass conservation, and the influences of different electric field directions and different voltages on the movement of bubbles in liquid are considered. The results show that the ratio of the length to axis is proportional to the strength of the electric field when the air bubble is stretched into an ellipsoid along the electric field line under the action of electrostatic gravity and surface tension. In addition, the bubble rising speed is affected by the electric field, the vertical electric field accelerates the bubble rise, and the horizontal direction slows it down.
