COMPUTER SIMULATION OF TWO-DIMENSIONAL TEMPERATURE FIELD FOR CASTING SOLIDIFICATION

Based on the mathematical model of Fourier Heat Conductivity Equations, this paper firstchooses an alternating-direction-implicit (ADI) type of finite difference method as the numericalcalculation method which is absolutely convergent and stable for time steps of any size when usedto calculate the temperature field. Then, the generalized programs are designed with FORTRANlanguage which can be used to calculate two-dimensional temperature field of casting solidificationof cast steel or cast iron for any shape. The calculated results about their solidification tendencyobtained from the above programs are found in good agreement with the experimental ones. Inaddition, we have measured the thermal properties of the furane resin sand which is the mouldingmaterial to make steel ingot moulds.

Numerical Study on the Two-Dimensional Temperature Fields of Titanium/Aluminum Double-Layer Target Irradiated by High-Intensity Pulsed Ion Beam

Interaction between high-intensity pulsed ion beam （HIPIB） and a double-layer target with titanium film on top of aluminum substrate was simulated. The two-dimensional nonlinear thermal conduction equations, with the deposited energy in the target taken as source term, were derived and solved by finite differential method. As a result, the two-dimensional spatial and temporal evolution profiles of temperature were obtained for a titanium/aluminum double-layer target irradiated by a pulse of HIPIB. The effects of ion beam current density on the phase state of the target materials near the film and substrate interface were analyzed. Both titanium and aluminum were melted near the interface after a shot when the ion beam current density fell in the range of 100 A/cm2 to 200 A/cm2.

A two-dimensional air streamer discharge modified model based on artificial stability term under non-uniform electric field at low temperature and sub-atmospheric pressure

In this paper, an improved air discharge fluid model under non-uniform electric field is constructed based on the plasma module COMSOL Multiphysics with artificial stability term, and the boundary conditions developed in the previous paper are applied to the calculation of photoionization rate. Based on the modified model, the characteristics of low temperature subatmospheric air discharge under 13 kV direct current voltage are discussed, including needle-plate and needle-needle electrode structures. Firstly, in order to verify the reliability of the model, a numerical example and an experimental verification were carried out for the modified model respectively. Both verification results show that the model can ensure the accuracy and repeatability of the calculation. Secondly, according to the calculation results of the modified model, under the same voltage and spacing, the reduced electric field under low temperature subatmosphere pressure is larger than that under normal temperature and atmospheric pressure. The high electric field leads to the air discharge at low temperature and sub atmospheric pressure entering the streamer initiation stage earlier, and has a faster propagation speed in the streamer development stage, which shortens the overall discharge time. Finally, the discharge characteristics of the two electrode structures are compared, and it is found that the biggest difference between them is that there is a pre-ionization region near the cathode in the needle-needle electrode structure. When the pre-ionization level reaches 1013 cm-3, the propagation speed of the positive streamer remains unchanged throughout the discharge process, and is no longer affected by the negative streamer. The peak value of electric field decreases with the increase of pre-ionization level, and tends to be constant during streamer propagation. Based on the previous paper, this paper constructs the air discharge model under non-uniform electric field, complements with the previous paper, and forms a relatively complete set of air discharge simulation system under low temperature and sub atmospheric pressure, which provides a certain reference for future research.

Numerical simulation of melt flow and temperature field during DC casting 2024 aluminium alloy under different casting conditions

Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during casting,which are crucial for the quality of the ingot and can determine the success or failure of the casting operation.Numerical simulation,with the advantages of low cost,rapid execution,and visualized results,is an important method to study and optimize the DC casting process.In the present work,a simulation model of DC casting 2024 aluminum alloy was established,and the reliability of the model was verified.Then,the influence of casting parameters on flow field and temperature field was studied in detail by numerical simulation method.Results show that with the increase of casting speed,the melt flow becomes faster,the depths of slurry zone and mushy zone increase,and the variation of slurry zone depth is greater than that of mushy zone.With an increase in casting temperature,the melt flow rate increases,the depth of the slurry zone becomes shallower,and the depth of the mushy zone experiences only minor changes.The simulation results further indicate that the increase of the flow rate of the secondary cooling water slightly reduces the depths of both slurry and mushy zone.

Stability analysis of a liquid crystal elastomer self-oscillator under a linear temperature field

Self-oscillating systems abound in the natural world and offer substantial potential for applications in controllers,micro-motors,medical equipments,and so on.Currently,numerical methods have been widely utilized for obtaining the characteristics of self-oscillation including amplitude and frequency.However,numerical methods are burdened by intricate computations and limited precision,hindering comprehensive investigations into self-oscillating systems.In this paper,the stability of a liquid crystal elastomer fiber self-oscillating system under a linear temperature field is studied,and analytical solutions for the amplitude and frequency are determined.Initially,we establish the governing equations of self-oscillation,elucidate two motion regimes,and reveal the underlying mechanism.Subsequently,we conduct a stability analysis and employ a multi-scale method to obtain the analytical solutions for the amplitude and frequency.The results show agreement between the multi-scale and numerical methods.This research contributes to the examination of diverse self-oscillating systems and advances the theoretical analysis of self-oscillating systems rooted in active materials.

Investigation on Temperature Field Calibration and Analysis of Wind Tunnel

For wind tunnels,it is essential to conduct temperature and flow field calibration on their test section,which is an important indicator for evaluating the quality of wind tunnel flow fields.In the paper,a truss composed of temperature sensors was used to calibrate the temperature field of a completed wind tunnel section.By adjusting the distance between the temperature measurement truss and the nozzle,as well as the wind speed,the temperature field distribution data at different positions could be obtained.Analyze these data to identify important factors that affect the distribution of temperature field.Simultaneously,the temperature field of the wind tunnel was simulated accordingly.The purpose is to further analyze the fluid heat transfer between air and wind tunnel walls through numerical simulation.Through the above analysis methods,the quality of the temperature field in the wind tunnel has been further verified,providing reference for future wind tunnel tests of relevant models.

Two-dimensional investigation of characteristic parameters and their gradients for the self-generated electric and magnetic fields of laser-induced zirconium plasma

Two-dimensional diagnosis of laser-induced zirconium(Zr)plasma has been experimentally performed using the time-of-flight method by employing Faraday cups in addition to electric and magnetic probes.The characteristic parameters of laser-induced Zr plasma have been evaluated as a function of different laser irradiances ranging from 4.5 to 11.7 GW cm-2 at different axial positions of 1–4 cm with a fixed radial distance of 2 cm.A well-supporting correlation between the plume parameters and the laser-plasma-produced spontaneous electric and magnetic(E and B)fields was established.The measurements of the characteristic parameters and spontaneously induced fields were observed to have an increasing trend with the increasing laser irradiance.However,when increasing the spatial distance in both the axial and radial directions,the plasma parameters(electron/ion number density,temperature and kinetic energy)did not show either continuously increasing or decreasing trends due to various kinetic and dynamic processes during the spatial evolution of the plume.However,the E and B fields were observed to be always diffusing away from the target.The radial component of electron number densities remained higher than the axial number density component,whereas the axial ion number density at all laser irradiances and axial distances remained higher than the radial ion number density.The higher axial self-generated electric field(SGEF)values than radial SGEF values are correlated with the effective charge-separation mechanism of electrons and ions.The generation of a self-generated magnetic field is observed dominantly in the radial direction at increasing laser irradiance as compared to the axial one due to the deflection of fast-moving electrons and the persistence of two-electron temperature on the radial axis.

An Explanation of the Temperature-Dependent Upper Critical Field Data of H3S on the Basis of the Thermodynamics of a Superconductor in a Magnetic Field

Excellent fits to a couple of the data-sets on the temperature (T)-dependent upper critical field (Hc2) of H3S (critical temperature, Tc ≈ 200 K at pressure ≈ 150 GPa) reported by Mozaffari, et al. (2019) were obtained by Talantsev (2019) in an approach based on an ingenious mix of the Ginzberg-Landau (GL), the Werthamer, Helfand and Hohenberg (WHH), and the Gor’kov, etc., theories which have individually been employed for the same purpose for a long time. Up to the lowest temperature (TL) in each of these data-sets, similarly accurate fits have also been obtained by Malik and Varma (2023) in a radically different approach based on the Bethe-Salpeter equation (BSE) supplemented by the Matsubara and the Landau quantization prescriptions. For T TL, however, while the (GL, WHH, etc.)-based approach leads to Hc2(0) ≈ 100 T, the BSE-based approach leads to about twice this value even at 1 K. In this paper, a fit to one of the said data-sets is obtained for the first time via a thermodynamic approach which, up to TL, is as good as those obtained via the earlier approaches. While this is interesting per se, another significant result of this paper is that for T TL it corroborates the result of the BSE-based approach.

Temperature field calculation of rail flash welding

The forging stage of rail flash welding has a decisive influence on joint strength,and the study of the temperature distribution in the process has an important role in further improving joint strength.In this paper,three calculation methods for the temperature field are given.First,the finite element model of the temperature field before forging rail flash welding is established by using the transient heat module of Ansys software and verified by infrared temperature measurement.Second,the temperature distribution of different parts of the rail before flash welding is obtained by using infrared thermal imaging equipment.Third,Matlab software is used to calculate the temperature of the non-measured part.Finally,the temperature distribution function along the rail axis is fitted through the temperature measurement data.The temperature distribution before the top forging of the rail flash welding can be used to analyze the joint and heat-affected zone organization and properties effectively and to guide the parameter setting and industrial production.

Interfacial built-in electric field and crosslinking pathways enabling WS_(2)/Ti_(3)C_(2)T_(x) heterojunction with robust sodium storage at low temperature

Developing efficient energy storage for sodium-ion batteries(SIBs)by creating high-performance heterojunctions and understanding their interfacial interaction at the atomic/molecular level holds promise but is also challenging.Besides,sluggish reaction kinetics at low temperatures restrict the operation of SIBs in cold climates.Herein,cross-linking nanoarchitectonics of WS_(2)/Ti_(3)C_(2)T_(x) heterojunction,featuring built-in electric field(BIEF),have been developed,employing as a model to reveal the positive effect of heterojunction design and BIEF for modifying the reaction kinetics and electrochemical activity.Particularly,the theoretical analysis manifests the discrepancy in work functions leads to the electronic flow from the electron-rich Ti_(3)C_(2)T_(x) to layered WS_(2),spontaneously forming the BIEF and“ion reservoir”at the heterogeneous interface.Besides,the generation of cross-linking pathways further promotes the transportation of electrons/ions,which guarantees rapid diffusion kinetics and excellent structure coupling.Consequently,superior sodium storage performance is obtained for the WS_(2)/Ti_(3)C_(2)T_(x) heterojunction,with only 0.2%decay per cycle at 5.0 A g^(-1)(25℃)up to 1000 cycles and a high capacity of 293.5 mA h g^(-1)(0.1A g^(-1)after 100 cycles)even at-20℃.Importantly,the spontaneously formed BIEF,accompanied by“ion reservoir”,in heterojunction provides deep understandings of the correlation between structure fabricated and performance obtained.

First-principles investigation of the significant anisotropy and ultrahigh thermoelectric efficiency of a novel two-dimensional Ga_(2)I_(2)S_(2) at room temperature

Two-dimensional(2D)thermoelectric(TE)materials have been widely developed;however,some 2D materials exhibit isotropic phonon,electron transport properties,and poor TE performance,which limit their application scope.Thus,exploring excellent anisotropic and ultrahigh-performance TE materials are very warranted.Herein,we first investigate the phonon thermal and TE properties of a novel 2D-connectivity ternary compound named Ga2I2S2.This paper comprehensively studies the phonon dispersion,phonon anharmonicity,lattice thermal conductivity,electronic structure,carrier mobility,Seebeck coefficient,electrical conductivity,and the dimensionless figure of merit(ZT)versus carrier concentration for 2D Ga_(2)I_(2)S_(2).We conclude that the in-plane lattice thermal conductivities of Ga_(2)I_(2)S_(2) at room temperature(300 K)are found to be 1.55 W mK^(−1) in the X-axis direction(xx-direction)and 3.82 W mK^(−1)in the Y-axis direction(yy-direction),which means its anisotropy ratio reaches 1.46.Simultaneously,the TE performance of p-type and n-type doping 2D Ga2I2S2 also shows significant anisotropy,giving rise to the ZT peak values of p-type doping in xx-and yy-directions being 0.81 and 1.99,respectively,and those of n-type doping reach ultrahigh values of 7.12 and 2.89 at 300 K,which are obviously higher than the reported values for p-type and n-type doping ternary compound Sn2BiX(ZT∼1.70 and∼2.45 at 300 K)(2020 Nano Energy 67104283).This work demonstrates that 2D Ga_(2)I_(2)S_(2) has high anisotropic TE conversion efficiency and can also be used as a new potential room-temperature TE material.

Temperature and stress fields in electron beam welded Ti-15-3 alloy to 304 stainless steel joint with copper interlayer sheet

Electron beam welding of Ti-15-3 alloy to 304 stainless steel （STS） using a copper filler metal was carried out. The temperature fields and stress distributions in the Ti/Fe and Ti/Cu/Fe joint during the welding process were numerically simulated and experimentally measured. The results show that the rotated parabola body heat source is fit for the simulation of the electron beam welding. The temperature distribution is asymmetric along the weld center and the temperature in the titanium alloy plate is higher than that in the 304 STS plate. The thermal stress also appears to be in asymmetric distribution. The residual tensile stress mainly exists in the weld at the 304 STS side. The copper filler metal decreases the peak temperature and temperature grade in the joint as well as the residual stress. The longitudinal and lateral residual tensile strengths reduce by 66 MPa and 31 MPa, respectively. From the temperature and residual stress, it is concluded that copper is a good filler metal candidate for the electron beam welding of Ti-15-3 titanium alloy to 304 stainless steel.

NUMERICAL SIMULATIONS OF TEMPERATURE FIELD IN DIRECT METAL LASER SINTERING PROCESS

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]

Analysis of Static Temperature Field of Vehicle's Solid Rubber Tire

Aim To analyse the static temperature field ofthe solid rubber tire(SRT).Methods The mechanical and thermal FE models were developed and analyzed respectively with the FE software ANSYS.Results The maximum temperature becomes higher with the higher with the higher velocity of tire and scales down slightly with the higher convection coefficients.The mixed models are reasonable.Conclusion The study on static temperature field is important and reasonable.It gives the fundament for life analysis of SRT.

Measurement and theoretical analysis of solar temperature field in steel-concrete composite girder

The solar temperature field of a large three-span continuous bridge with steel-concrete composite box girder and variable cross-section is measured to verify a calculation method for the temperature field of steel-concrete composite beams. The test results show that the temperature of an external steel web- plate is higher than that of an internal web-plate due to the difference in solar radiation. Air temperature inside the box matches the average temperature of the whole steel box. Based on actual measurements, a transient thermal analysis with multiple boundary conditions is also carried out by a software program ANSYS. Convective boundary situation and states of solar radiation on steel web plates in different situations are determined in the analysis. The feature of the temperature field is preliminarily achieved through a comparative study between the actual measurement and the finite element analysis. The computed results are in good consistence with the actual measurement results, with the maximum difference within 2 ℃. This indicates that the theoretical calculation method is reliable and it provides a foundation for further research on temperature field distribution in the steel-concrete composite box girder.

THE ANALYSIS OF TEMPERATURE FIELD OFARMORED VEHICLE

论述了负重轮轮胎温度场有限元分析方法；给出了生热率数学模型；建立了具有三重非线性的负重轮力学有限元分析模型和热学分析模型，利用有限元分析软件进行了温度场的计算分析；研究了负重轮轮胎温度场，分析了温升过程，为负重轮寿命分析提供了理论和数据。

STUDY ON TEMPERATURE FIELD INDUCED IN HIGH FREQUENCY INDUCTION HEATING

A mathematical model was established for the temperature field developed during high frequency induction heating （HFIH） by Maxwell＇s equations. It required solving the coupled equations of the electromagnetic and temperature fields. The numerical simudation was performed using FEMLAB. The comparison of the calculations using the proposed model with experimental results showed a very good correlation. The effects of the heating parameters in high frequency induction such as the distance between the plate and the coil, the applied current, the frequency, and the turns of the coil on the temperature profiles developed in the plate were also discussed using the established model.

Numerical Simulation of Multi-track and Multi-layer Temperature Field on Laser Direct Metal Shaping

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.

Numerical Simulation of Temperature Field and Thermal Stress Field of Work Roll During Hot Strip Rolling

Based on the thermal conduction equations, the three-dimensional （3D） temperature field of a work roll was investigated using finite element method （FEM）. The variations in the surface temperature of the work roll during hot strip rolling were described, and the thermal stress field of the work roll was also analyzed. The results showed that the highest roll surface temperature is 593 ℃, and the difference between the minimum and maximum values of thermal stress of the work roll surface is 145.7 MPa. Furthermore, the results of this analysis indicate that temperature and thermal stress are useful parameters for the investigation of roll thermal fatigue and also for improving the quality of strip during rolling.

COUPLED NUMERICAL SIMULATION ON COLD ROLLER'S TEMPERATURE FIFLD-PHASE TRANSFORMATION - STRESS FIELD DURING ITS QUENCHING PSOCESS

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.