Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the inf...Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.展开更多
To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dime...To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.展开更多
This paper analyzes the sources of heat losses in magnetic fluid bearings,proposes various cou-pling relationships of physical fields,divides the coupled heat transfer surfaces while ensuring the continuity of heat fl...This paper analyzes the sources of heat losses in magnetic fluid bearings,proposes various cou-pling relationships of physical fields,divides the coupled heat transfer surfaces while ensuring the continuity of heat flux density,and analyzes the overall heat dissipation pathways of the bearings.By changing parameters such as input current,rotor speed,and inlet oil flow rate,the study applies a multi-physics field coupling method to investigate the influence of different parameters on the temper-ature field and heat dissipation patterns of the bearings,which is then validated through experi-ments.This research provides a theoretical basis for the optimal design of magnetic fluid bearing sys-tems.展开更多
A transient multi-physics model incorporated with an electromagneto-thermomechanical coupling is developed to capture the multi-field behavior of a single-pancake(SP)insert no-insulation(NI)coil in a hybrid magnet dur...A transient multi-physics model incorporated with an electromagneto-thermomechanical coupling is developed to capture the multi-field behavior of a single-pancake(SP)insert no-insulation(NI)coil in a hybrid magnet during the charging and discharging processes.The coupled problem is resolved by means of the finite element method(FEM)for the magneto-thermo-elastic behaviors and the Runge-Kutta method for the transient responses of the electrical circuits of the hybrid superconducting magnet system.The results reveal that the transient multi-physics responses of the insert NI coil primarily depend on the charging/discharging procedure of the hybrid magnet.Moreover,a reverse azimuthal current and a compressive hoop stress are induced in the insert NI coil during the charging process,while a forward azimuthal current and a tensile hoop stress are observed during the discharging process.The induced voltages in the insert NI coil can drive the currents flowing across the radial turns where the contact resistance exists.Therefore,it brings forth significant Joule heat,causing a temperature rise and a uniform distribution of this heat in the coil turns.Accordingly,a thermally/mechanically unstable or quenching event may be encountered when a high operating current is flowing in the insert NI coil.It is numerically predicted that a quick charging will induce a compressive hoop stress which may bring a risk of buckling instability in the coil,while a discharging will not.The simulations provide an insight of hybrid superconducting magnets under transient start-up or shutdown phases which are inevitably encountered in practical applications.展开更多
We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numeri...We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numerical simulation and experimental validation. The distribution characteristics and interaction mechanism of the flow field and electromagnetic field of the ICP wind tunnel are investigated using the multi-field coupling method of flow, electromagnetic, chemical, and thermodynamic field. The accuracy of the numerical simulation is validated by comparing the experimental results with the simulation results. Thereafter, the wind tunnel pressure, air velocity, electron density, Joule heating rate, Lorentz force, and electric field intensity obtained using the simulation are analyzed and discussed. The results indicate that for the 1.2-MW ICP wind tunnel, the maximum values of temperature, pressure, electron number density, and other parameters are observed during coil heating. The influence of the radial Lorentz force on the momentum transfer is stronger than that of the axial Lorentz force. The electron number density at the central axis and the amplitude and position of the Joule heating rate are affected by the radial Lorentz force. Moreover, the plasma in the wind tunnel is constantly in the subsonic flow state, and a strong eddy flow is easily generated at the inlet of the wind tunnel.展开更多
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.展开更多
The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and nu...The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and numerical calculations.It is found that the longitudinal magnetic field can affect the heating by changing the level of the pressure heating along the longitudinal direction and that of the Ohmic heating along the direction which is perpendicular to both driving electric field and the applied transverse magnetic field,and a continuously increased longitudinal magnetic field can induce pressure heating to become dominant.Moreover,the electron temperature as well as proportion of some low energy electrons will increase if a small longitudinal magnetic field is introduced,which is attributed to the increased average electron energy.We believe that the research will provide guidance for optimizing the magnetic field configuration of some discharge systems having both transverse and longitudinal magnetic field.展开更多
The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field...The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.展开更多
In order to study the temperature distribution of deep field,mathematical mod- els of temperature field in field and surrounding rock were built based on heat transfer and seepage theory.Combined test data with mathem...In order to study the temperature distribution of deep field,mathematical mod- els of temperature field in field and surrounding rock were built based on heat transfer and seepage theory.Combined test data with mathematical model,the temperature distribution under heat-transfer and underground-water coupling was studied by using Golden Soft- ware Surfer and Matlab.The results show that distribution law of most isothermal lines is very similar in deep field,and temperature gradient is equal in general.At the same time, temperature distribution is influenced by underground-water and fault.In surrounding rock, seepage changes symmetrical distribution of temperature field and vector,and the tem- perature field may divide into inward-diffusion area and outward-diffusion area.Peripheral temperature of working will approach to the temperature of airflow.In inward diffusion area the distribution of temperature and temperature vector is symmetric,and the direction of temperature vector point to the center of working.The action of airflow is stronger than seepage in inward diffusion area,however,the case opposite is true in outward diffusion area.展开更多
In this paper, the development status of casting numerical simulation technology is introduced. In additional, mathematical model, solution method, initial condition, boundary condition and defect predicting scheme of...In this paper, the development status of casting numerical simulation technology is introduced. In additional, mathematical model, solution method, initial condition, boundary condition and defect predicting scheme of foundry process are also analyzed, which include the mold filling process, solidification process and the process coupling fluid flow with heat transfer. Finally, a practical casting is taken out to show how to predict defects and optimize foundry process with numerical simulation technology.展开更多
Coupled turbulent flow, temperature fields of the twin-roll casting strip process were simulated by three-dimensional finite element method. Based on the heat balance calculation and using inverse methods between the ...Coupled turbulent flow, temperature fields of the twin-roll casting strip process were simulated by three-dimensional finite element method. Based on the heat balance calculation and using inverse methods between the simulations and real experiments, the relational models among casting speed, location, and coefficient of heat transfer between molten metal and rolls in different regions are given. In the simulation, the calculated surface temperatures are in good agreement with the measured values. An on-line model of kiss point is derived by simulations and the geometry of molten pool, corresponding control strategy is also proposed.展开更多
To prepare Carbon/Carbon (C/C) composites with advanced performance, the thermal gradient chemical vapor infiltration (TCVI) process has been optimized by simulation. A 2D axisymmetric unstable model was built, which ...To prepare Carbon/Carbon (C/C) composites with advanced performance, the thermal gradient chemical vapor infiltration (TCVI) process has been optimized by simulation. A 2D axisymmetric unstable model was built, which included convection, conduction, diffusion, densification reactions in the pores and the evolution of the porous medium. The multi-physical field coupling model was solved by finite element method (FEM) and iterative calculation. The time evolution of the fluid, temperature and preform density field were obtained by the calculation. It is indicated that convection strongly affects the temperature field. For the preform of carbon/carbon composites infiltrated for 100 h by TCVI, the radial average densities from simulation agrees well with those from experiment. The model is validated to be reliable and the simulation has capability of forecasting the process.展开更多
The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, ...The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, error function and normal distribution function, and a series of FEM equations of coupled fields of temperature and seepage were deduced and put forward. With the example of shaft ground freezing, the formation of freezing wall in seepage ground was simulated.展开更多
<div style="text-align:justify;"> Rainfall infiltration is a porous medium flow problem with variable saturation. Based on the theoretical analysis of the flow field, electrical conductivity of rocks, ...<div style="text-align:justify;"> Rainfall infiltration is a porous medium flow problem with variable saturation. Based on the theoretical analysis of the flow field, electrical conductivity of rocks, the electrical field, the paper simulates the coupling relationship between the water saturation in soil and the apparent resistivity distribution. It combines the Richards equation, the Archie formula and the Laplace equation. The experiment simulates the potential field data by the Wenner setting in electrical exploration on a two-layer geologic model with continuous rainfall during 5 days, which shows that the effective saturation in soil is increasing with the rainfall time, while the apparent resistivity is decreasing. This can provide a theoretical basis for the analyzing the rainfall infiltration and porosity of the soil by using high-density electrical method in the future. </div>展开更多
Purpose-It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.Design/methodology/approach-T...Purpose-It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.Design/methodology/approach-The paper builds up the model for the hydrothermal coupling calculation of permafrost using finite element software COMSOL to study how permafrost temperature field changes in the short term after subgrade filling,on which basis it proposes the method of calculation for the concave distortion of freezing front in the subgrade-covered area.Findings-The results show that the freezing front below the subgrade center sinks due to the thermal effect of subgrade filling,which will trigger hydrothermal erosion in case of sufficient moisture inflows,leading to the thawing settlement or the cracking of the subgrade,etc.The heat output of soil will be hindered the most in case of July filling,in which case the sinking and the distortion of the freezing front is found to be the most severe,which the recovery of the permafrost temperature field,the slowest,constituting the most unfavorable working condition.The concave distortion of the freezing front in the subgrade area increases with the increase in temperature difference between the filler and ground surface,the subgrade height,the subgrade width and the volumetric thermal capacity of filler,while decreases with the increase of the thermal conductivity of filler.Therefore,the filler chose for engineering project shall be of small volumetric thermal capacity,low initial temperature and high thermal conductivity whenever possible.Originality/value-The concave distortion of the freezing front under different working conditions at different times after filling can be calculated using the method proposed.展开更多
文摘Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.
文摘To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.
基金the National Natural Science Foundation of China(No.52075468)the Natural Science Foundation of Hebei Province(No.E2020203052)+1 种基金the Key Scientific Research Projects of North China University of Technology(No.ZD-YG-202306-23)the Tangshan Science and Technology Project(No.23130201E).
文摘This paper analyzes the sources of heat losses in magnetic fluid bearings,proposes various cou-pling relationships of physical fields,divides the coupled heat transfer surfaces while ensuring the continuity of heat flux density,and analyzes the overall heat dissipation pathways of the bearings.By changing parameters such as input current,rotor speed,and inlet oil flow rate,the study applies a multi-physics field coupling method to investigate the influence of different parameters on the temper-ature field and heat dissipation patterns of the bearings,which is then validated through experi-ments.This research provides a theoretical basis for the optimal design of magnetic fluid bearing sys-tems.
基金the National Natural Science Foundation of China(Nos.11932008 and 11672120)the Fundamental Research Funds for the Central Universities of China(No.lzujbky-2022-kb01)。
文摘A transient multi-physics model incorporated with an electromagneto-thermomechanical coupling is developed to capture the multi-field behavior of a single-pancake(SP)insert no-insulation(NI)coil in a hybrid magnet during the charging and discharging processes.The coupled problem is resolved by means of the finite element method(FEM)for the magneto-thermo-elastic behaviors and the Runge-Kutta method for the transient responses of the electrical circuits of the hybrid superconducting magnet system.The results reveal that the transient multi-physics responses of the insert NI coil primarily depend on the charging/discharging procedure of the hybrid magnet.Moreover,a reverse azimuthal current and a compressive hoop stress are induced in the insert NI coil during the charging process,while a forward azimuthal current and a tensile hoop stress are observed during the discharging process.The induced voltages in the insert NI coil can drive the currents flowing across the radial turns where the contact resistance exists.Therefore,it brings forth significant Joule heat,causing a temperature rise and a uniform distribution of this heat in the coil turns.Accordingly,a thermally/mechanically unstable or quenching event may be encountered when a high operating current is flowing in the insert NI coil.It is numerically predicted that a quick charging will induce a compressive hoop stress which may bring a risk of buckling instability in the coil,while a discharging will not.The simulations provide an insight of hybrid superconducting magnets under transient start-up or shutdown phases which are inevitably encountered in practical applications.
基金supported by the National Natural Science Foundation of China (Grant No. 11705143)the Open Foundation for Key Laboratories of National Defense Science and Technology of China (Grant No. 6142202031901)the Foundation for Research and Development of Applied Technology in Beilin District of Xi’an,China (Grant No. GX2047)。
文摘We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numerical simulation and experimental validation. The distribution characteristics and interaction mechanism of the flow field and electromagnetic field of the ICP wind tunnel are investigated using the multi-field coupling method of flow, electromagnetic, chemical, and thermodynamic field. The accuracy of the numerical simulation is validated by comparing the experimental results with the simulation results. Thereafter, the wind tunnel pressure, air velocity, electron density, Joule heating rate, Lorentz force, and electric field intensity obtained using the simulation are analyzed and discussed. The results indicate that for the 1.2-MW ICP wind tunnel, the maximum values of temperature, pressure, electron number density, and other parameters are observed during coil heating. The influence of the radial Lorentz force on the momentum transfer is stronger than that of the axial Lorentz force. The electron number density at the central axis and the amplitude and position of the Joule heating rate are affected by the radial Lorentz force. Moreover, the plasma in the wind tunnel is constantly in the subsonic flow state, and a strong eddy flow is easily generated at the inlet of the wind tunnel.
基金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.
基金supported by National Natural Science Foundation of China (Nos. 11975163 and 12175160)a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)
文摘The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and numerical calculations.It is found that the longitudinal magnetic field can affect the heating by changing the level of the pressure heating along the longitudinal direction and that of the Ohmic heating along the direction which is perpendicular to both driving electric field and the applied transverse magnetic field,and a continuously increased longitudinal magnetic field can induce pressure heating to become dominant.Moreover,the electron temperature as well as proportion of some low energy electrons will increase if a small longitudinal magnetic field is introduced,which is attributed to the increased average electron energy.We believe that the research will provide guidance for optimizing the magnetic field configuration of some discharge systems having both transverse and longitudinal magnetic field.
文摘The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.
基金Foundation of Liaoning Educational Committee(2007F049,20060390)
文摘In order to study the temperature distribution of deep field,mathematical mod- els of temperature field in field and surrounding rock were built based on heat transfer and seepage theory.Combined test data with mathematical model,the temperature distribution under heat-transfer and underground-water coupling was studied by using Golden Soft- ware Surfer and Matlab.The results show that distribution law of most isothermal lines is very similar in deep field,and temperature gradient is equal in general.At the same time, temperature distribution is influenced by underground-water and fault.In surrounding rock, seepage changes symmetrical distribution of temperature field and vector,and the tem- perature field may divide into inward-diffusion area and outward-diffusion area.Peripheral temperature of working will approach to the temperature of airflow.In inward diffusion area the distribution of temperature and temperature vector is symmetric,and the direction of temperature vector point to the center of working.The action of airflow is stronger than seepage in inward diffusion area,however,the case opposite is true in outward diffusion area.
文摘In this paper, the development status of casting numerical simulation technology is introduced. In additional, mathematical model, solution method, initial condition, boundary condition and defect predicting scheme of foundry process are also analyzed, which include the mold filling process, solidification process and the process coupling fluid flow with heat transfer. Finally, a practical casting is taken out to show how to predict defects and optimize foundry process with numerical simulation technology.
基金supported by National Key Research Development Planning Project of China (2004CB619108).
文摘Coupled turbulent flow, temperature fields of the twin-roll casting strip process were simulated by three-dimensional finite element method. Based on the heat balance calculation and using inverse methods between the simulations and real experiments, the relational models among casting speed, location, and coefficient of heat transfer between molten metal and rolls in different regions are given. In the simulation, the calculated surface temperatures are in good agreement with the measured values. An on-line model of kiss point is derived by simulations and the geometry of molten pool, corresponding control strategy is also proposed.
基金Supported by the National Natural Science Foundation of China (Grant No. 90716024)
文摘To prepare Carbon/Carbon (C/C) composites with advanced performance, the thermal gradient chemical vapor infiltration (TCVI) process has been optimized by simulation. A 2D axisymmetric unstable model was built, which included convection, conduction, diffusion, densification reactions in the pores and the evolution of the porous medium. The multi-physical field coupling model was solved by finite element method (FEM) and iterative calculation. The time evolution of the fluid, temperature and preform density field were obtained by the calculation. It is indicated that convection strongly affects the temperature field. For the preform of carbon/carbon composites infiltrated for 100 h by TCVI, the radial average densities from simulation agrees well with those from experiment. The model is validated to be reliable and the simulation has capability of forecasting the process.
文摘The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, error function and normal distribution function, and a series of FEM equations of coupled fields of temperature and seepage were deduced and put forward. With the example of shaft ground freezing, the formation of freezing wall in seepage ground was simulated.
文摘<div style="text-align:justify;"> Rainfall infiltration is a porous medium flow problem with variable saturation. Based on the theoretical analysis of the flow field, electrical conductivity of rocks, the electrical field, the paper simulates the coupling relationship between the water saturation in soil and the apparent resistivity distribution. It combines the Richards equation, the Archie formula and the Laplace equation. The experiment simulates the potential field data by the Wenner setting in electrical exploration on a two-layer geologic model with continuous rainfall during 5 days, which shows that the effective saturation in soil is increasing with the rainfall time, while the apparent resistivity is decreasing. This can provide a theoretical basis for the analyzing the rainfall infiltration and porosity of the soil by using high-density electrical method in the future. </div>
基金supported by the Fund of China Academy of Railway Sciences Corporation Limited (2019YJ041).
文摘Purpose-It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.Design/methodology/approach-The paper builds up the model for the hydrothermal coupling calculation of permafrost using finite element software COMSOL to study how permafrost temperature field changes in the short term after subgrade filling,on which basis it proposes the method of calculation for the concave distortion of freezing front in the subgrade-covered area.Findings-The results show that the freezing front below the subgrade center sinks due to the thermal effect of subgrade filling,which will trigger hydrothermal erosion in case of sufficient moisture inflows,leading to the thawing settlement or the cracking of the subgrade,etc.The heat output of soil will be hindered the most in case of July filling,in which case the sinking and the distortion of the freezing front is found to be the most severe,which the recovery of the permafrost temperature field,the slowest,constituting the most unfavorable working condition.The concave distortion of the freezing front in the subgrade area increases with the increase in temperature difference between the filler and ground surface,the subgrade height,the subgrade width and the volumetric thermal capacity of filler,while decreases with the increase of the thermal conductivity of filler.Therefore,the filler chose for engineering project shall be of small volumetric thermal capacity,low initial temperature and high thermal conductivity whenever possible.Originality/value-The concave distortion of the freezing front under different working conditions at different times after filling can be calculated using the method proposed.