Numerical simulation of flow and heat transfer of n-decane in sub-millimeter spiral tube at supercritical pressure

The flow and heat transfer characteristics of n-decane in the sub-millimeter spiral tube(SMST) at supercritical pressure(p = 3 MPa) are studied by the RNG k-ε numerical model in this paper. The effects of various Reynolds numbers(Re) and structural parameters pitch(s) and spiral diameter(D) are analyzed.Results indicate that the average Nusselt numberNu and friction factorNu increase with an increase in Re, and decrease with an increase in D/d(tube diameter). In terms of the structural parameter s/d, it is found that as s/d increases, the Nu first increase, and then decrease. and the critical structural parameter is s/d = 4. Compared with the straight tube, the SMST can improve Nu by 34.8% at best, while it can improve Nu by 102.1% at most. In addition, a comprehensive heat transfer coefficient is applied to analyze the thermodynamic properties of SMST. With the optimal structural parameters of D/d = 6 and s/d = 4, the comprehensive heat transfer factor of supercritical pressure hydrocarbon fuel in the SMST can reach 1.074. At last, correlations of the average Nusselt number and friction factor are developed to predict the flow and heat transfer of n-decane at supercritical pressure.

Heat transfer and temperature evolution in underground mininginduced overburden fracture and ground fissures: Optimal time window of UAV infrared monitoring

Heat transfer and temperature evolution in overburden fracture and ground fissures are one of the essential topics for the identification of ground fissures via unmanned aerial vehicle(UAV) infrared imager. In this study, discrete element software UDEC was employed to investigate the overburden fracture field under different mining conditions. Multiphysics software COMSOL were employed to investigate heat transfer and temperature evolution of overburden fracture and ground fissures under the influence of mining condition, fissure depth, fissure width, and month alternation. The UAV infrared field measurements also provided a calibration for numerical simulation. The results showed that for ground fissures connected to underground goaf(Fissure Ⅰ), the temperature difference increased with larger mining height and shallow buried depth. In addition, Fissure Ⅰ located in the boundary of the goaf have a greater temperature difference and is easier to be identified than fissures located above the mining goaf. For ground fissures having no connection to underground goaf(Fissure Ⅱ), the heat transfer is affected by the internal resistance of the overlying strata fracture when the depth of Fissure Ⅱ is greater than10 m, the temperature of Fissure Ⅱ gradually equals to the ground temperature as the fissures’ depth increases, and the fissures are difficult to be identified. The identification effect is most obvious for fissures larger than 16 cm under the same depth. In spring and summer, UAV infrared identification of mining fissures should be carried out during nighttime. This study provides the basis for the optimal time and season for the UAV infrared identification of different types of mining ground fissures.

Numerical Analysis on Temperature Distribution in a Single Cell of PEFC Operated at Higher Temperature by1D Heat Transfer Model and 3D Multi-Physics Simulation Model

This study is to understand the impact of operating conditions, especially initial operation temperature (Tini) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (Treact). In addition, to validate the proposed heat transfer model, Treact obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O2 and H2O from the inlet to the outlet is more even with the increase in Tini due to the lower performance of O2 reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in Tini and the value of current density is smaller with the increase in Tini due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in Treact from the inlet to the outlet is more even with the increase in Tini irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in Tini due to the lower performance of O2 reduction reaction.

Large-eddy Simulation of Fluid Flow and Heat Transfer in a Mixing Tee Junction

The temperature fluctuation caused by thermal striping phenomena of hot and cold fluids mixing results in cyclical thermal stress fatigue failure of the pipe wall. Mean temperature difference between hot and cold fluids was often used as thermal load in previous analysis of thermal fatigue failure, thereby the influences of the amplitude and frequency of temperature fluctuation on thermal fatigue failure were neglected. Based on the mechanism of flow and heat transfer which induces thermal fatigue, the turbulent mixing of hot and cold water in a tee junction is simulated with FLUENT platform by using the Large-eddy simulation（LES） turbulent flow model with the sub-grid scale（SGS） model of Smagorinsky-Lilly（SL） to capture the amplitude and frequency of temperature fluctuation. In a simulation case, hot water with temperature of 343.48 K and velocity of 0.15 m/s enters the horizontal main duct with the side length of 100 mm, while cold water with temperature of 296.78 K and velocity of 0.3 m/s enters the vertical branch duct with the side length of 50 mm. The numerical results show that the mean and fluctuating temperatures are in good agreement with the previous experimental data, which describes numerical simulation with high reliability and accuracy; the power spectrum density（PSD） on top wall is higher than that on bottom wall（as the frequency less than 1 Hz）, while the PSD on bottom wall is relatively higher than that on top wall （as the frequency of 1-10Hz）. The temperature fluctuations in full mixing region of the tee junction can be accurately captured by LES and can provide the theoretical basis for the thermal stress and thermal fatigue analyses.

Numerical Simulation of Fluid Flow and Heat Transfer in Funnel Shaped Mold of Thin Slab Continuous Caster

Thin slab casting is used widely in the world. The control of molten steel flow and solidification in the mold is difficult due to the high casting speed and complicated configuration of the mold. Numerical simulation was carried out to study the fluid flow and heat transfer in the funnel shaped mold. The influence of nozzle design, casting speed and nozzle submersion depth on the flow and temperature fields in the mold was investigated, and guidance for selecting configurations of submerged nozzle was obtained.

Direct numerical simulation of viscoelastic-fluid-based nanofluid turbulent channel flow with heat transfer

Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.

Paraffin–CaCl_(2)·6H_(2)O dosage effects on the strength and heat transfer characteristics of cemented tailings backfill

The challenge of high temperatures in deep mining remains harmful to the health of workers and their production efficiency The addition of phase change materials (PCMs) to filling slurry and the use of the cold storage function of these materials to reduce downhole temperatures is an effective approach to alleviate the aforementioned problem.Paraffin–CaCl_(2)·6H_(2)O composite PCM was prepared in the laboratory.The composition,phase change latent heat,thermal conductivity,and cemented tailing backfill (CTB) compressive strength of the new material were studied.The heat transfer characteristics and endothermic effect of the PCM were simulated using Fluent software.The results showed the following:(1) The new paraffin–CaCl_(2)·6H_(2)O composite PCM improved the thermal conductivity of native paraffin while avoiding the water solubility of CaCl_(2)·6H_(2)O.(2) The calculation formula of the thermal conductivity of CaCl_(2)·6H_(2)O combined with paraffin was deduced,and the reasons were explained in principle.(3) The“enthalpy–mass scale model”was applied to calculate the phase change latent heat of nonreactive composite PCMs.(4)The addition of the paraffin–CaCl_(2)·6H_(2)O composite PCM reduced the CTB strength but increased its heat absorption capacity.This research can give a theoretical foundation for the use of heat storage backfill in green mines.

Reaction characteristics of magnesium production under argon flow by silicothermic reduction and numerical simulation of argon entrainment process

In this study, the reaction characteristics of reduction of calcined dolomite with ferrosilicon under argon flow to produce magnesium were studied by conducting experiments Pidgeon pellets were used to study the effect of reduced temperature, argon flow, and reduced time on the conversion of calcined dolomite reduction by ferrosilicon. The results show that the conversion significantly increases with the increase in the reduction temperature and reduction time. The conversion first increases and then decreases with the increase in argon flow. The highest conversion was obtained when the argon flow rate was 3 L·min^(-1), and a nearly spherical shape, nanoscale magnesium powder was obtained. Then the characters of the circulating argon entrainment process were numerically studied by ANSYS Fluent 17. A physical model of multilayer pellet arrangement was established, and a numerical calculation model of chemical reaction, radiation, heat conduction, and convection heat transfer was constructed. This confirms that high-temperature argon can effectively strengthen the heat exchange between pellets, improve the heat transfer efficiency, and facilitate the pellets to react quickly. When the conversion is 80%, the production efficiency increased by about 28.6%. In addition, the magnesium production efficiency showed an increase tendency with the increase of the argon inlet flow rate.

Fluid Flow and Mixed Heat Transfer in a Horizontal Channel with an Open Cavity and Wavy Wall

Heat exchangers are utilized extensively in different industries and technologies.Consequently,optimizing heat exchangers has been a major concern among researchers.Although various studies have been conducted to improve the heat transfer rate,the use of a wavy wall in the presence of different types of heat transfer mechanisms has not been investigated.This study thus investigates the mixed heat transmission behavior of fluid in a horizontal channel with a cavity and a hot,wavy wall.The fluid flow in the channel is considered laminar,and the governing equations including continuity,momentum,and energy are all solved numerically.The numerical solution is stabilized by using a first-order multi-dimensional characteristic-based scheme in combination with a fifth-order Runge-Kutta method.The flow and heat transfer effects of varying Richardson numbers,Reynolds numbers,wave amplitude,wavelength,channel height,and cavity width are examined.The results indicate that the mean Nusselt number increases with an increase in Reynolds number,wave amplitude,and cavity width,while it decreases with an increase in Richardson number,wavelength,and channel height.The minimum Nusselt number is calculated to be 0.7,whereas the maximum Nusselt number is 27.09.The Nusselt number has only increased by 40%in the higher depths of the cavity,despite the Richardson number being 10,000 times larger.But this figure increases to 130%at lower depths.The mean Nusselt number is thus significantly influenced by channel height and cavity width.The influence of wave amplitude on the mean Nusselt number is twice that of wavelength.

Simulation of heat transfer in steel billets during continuous casting

This work is focused on the development of computational algorithms to create a simulator for solving the heat transfer during the continuous casting process of steel. The temperatures and the solid shell thickness profiles were calculated and displayed on the screen for a billet through a defined continuous casting plant （CCP）. The algorithms developed to calculate billet temperatures, involve the solutions of the corresponding equations for the heat removal conditions such as radiation, forced convection, and conduction according to the billet posi- tion through the CCP. This is done by a simultaneous comparison with the kinematics model previously developed. A finite difference method known as Cmnk-Nicholson is applied to solve the two-dimensional computational array （2D model）. Enthalpy （HIJ） and temperature （TIL） in every node are updated at each step time. The routines to display the results have been developed using a graphical user interface （GUI） in the programming language C＋＋. Finally, the results obtained are compared with those of industrial trials for the surface temperature of three steel casters with different plant configurations in different casting conditions.

A study of fire propagation in coal seam with numerical simulation of heat transfer and chemical reaction rate in mining field

Spontaneous combustion of coal is a problem that affects the mining operation and generates environ-mental,economic,social and geotechnical impacts.This phenomenon has been divided into two pro-cesses:ignition and propagation.Fire propagation develops in coal seams because of a set of factors such as direction and wind speed,fracturing and temperature.In this work,heat transfer and chemical kinetics are studied from conservation equations of energy and species,respectively,using the software COMSOL Multiphysics to simulate the propagation of fires in coal seams.Two possible scenarios were analyzed that usually occur in the walls of the coal seams,such as fire focus and fire complete screens.It was found that the propagation kinetics of the fire changes depending on the temperature,the fractur-ing of rock mass and the area of fire influence.For temperature values lower than 300℃,there is con-sumption around 250 cm^3/h,values around 700℃,the consumption is 1500 cm^3/h,and for fires of 1200℃ have values of 3000 cm^3/h.Depending on the speed of propagation can vary from 4 to 17cm/day,considering on the level and fracturing of the final wall of the open pit.

Numerical Simulation of Heat Transfer and Deformation of Initial Shell in Soft Contact Continuous Casting Mold Under High Frequency Electromagnetic Field

Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromagnetic parameters were obtained from the previous studies. Owing to the induction heating of a high frequency electromagnetic field （20 kHz）, the thickness of initial solidification shell decreases, and the temperature of strand surface and slit copper mold increases when compared with the case without the electromagnetic filed. The viscosity of flux de- creases because of the induction heating of the high frequency electromagnetic field, and the dimension of the flux channel increases with electromagnetic pressure; thus, the deformation behavior of initial solidification shell was different before and after the action of high frequency electromagnetic field. Furthermore, the abatement mechanism of oscillation marks under high frequency electromagnetic field was explained.

Study on Numerical Simulation of Mold Filling and Heat Transfer in Die Casting Process

A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.

Numerical Simulation of Three-dimensional Heat and Mass Transfer in Spray Cooling of Converter Gas in a Venturi Scrubber

In order to predict the pressure drop, collection efficiency, velocity, temperature and mole fraction of vapor in an industrial venturi scrubber with water spraying for converter gas cooling, a three-dimensional model of heat and mass transfer with phase change is established. The gas flow and liquid droplets are treated as a continuous phase with a Eulerian approach and as a discrete phase with a Lagrangian approach, respectively. The coupled problem of heat, force, and mass transfers between gas flow and liquid droplets is solved by a commercial computational fluid dynamics（CFD） package, FLUENT. The numerical results show that the water injections have an important influence on the distributions of pressure, velocity, temperature, and mole fraction of vapor, especially for the spraying region in the throat. In the spraying region, the pressure drop is higher and the velocity is lower than in other regions due to the gas-droplet drag, while the temperature is lower because the droplet absorbs large amounts of heat from the high temperature gas and the mole fraction of vapor is higher due to the phase change of the liquid droplet. A number of cases with different Water-to-gas volume flow ratios and baffle openings were simulated. The dependence of pressure drop, velocity, temperature, mole fraction of vapor, and collection efficiency on both the water-to-gas volume flow ratio and baffle opening are analyzed. The good agreements between simulation results and experiment data of pressure drop, temperature, and collection efficiency validate the model. The model should facilitate optimization of the venturi scrubber design in order to give better performance with lower pressure drops and higher collection efficiency.

SIMULATION OF STEEL COIL HEAT TRANSFER IN HPH FURNACE

The mathematical model has been estublished for the simulation of steel coil＇s heat transfer during annealing thermal process in HPH （high performance hydrogen） furnace. The equivalent radial thermal conductivity is adopted by statistical analysis regression approach through the combination of a large quantity of production data collected in practice and theoretical analyses. The effect of the number of coils on circulating flow gas is considered for calculating the convection heat transfer coefficient, The temperature within the coil is predicted with the developed model during the annealing cycle including heating process and cooling process. The good consistently between the predicted results and the experimental data has demonstrated that the mathematical model established and the parameters identified by this paper are scientifically feasible and the effective method of calculation for coil equivalent radial heat transfer coefficient and circulating gas flow has been identified successfully, which largely enhances the operability and feasibility of the mathematic- model. This model provides a theoretical basis and an effective means to conduct studies on the impact that foresaid factors may imposed on the steel coil＇s temperature field, to analyze the stress within coils, to realize online control and optimal production and to increase facilily output by increasing heating and cooling rates of coils without producing higher thermal stress.

A review on metal additive manufacturing:modeling and application of numerical simulation for heat and mass transfer and microstructure evolution

Metal additive manufacturing technology has been widely used in prototyping,parts manufacturing and repairing.Metal additive manufacturing is a multi-scale and multi-physical coupling process with complex physical phenomena of heat and mass transfer and microstructure evolution.It is hard to directly observe the dynamic behavior and microstructure evolution of molten pool during additive manufacturing.Therefore,numerical simulation of additive manufacturing process is significant since it can efficiently and pertinently predict and analyze the physical phenomena in the process of metal additive manufacturing,and provide a reference for technological parameters selection.In this review,the research progress of numerical simulation of metal additive manufacturing is discussed.Various aspects of numerical simulation models are reviewed,including:(1)Introduction of basic control method and physical description of numerical simulation models;(2)Comparison of various heat and mass transfer models based on different physical assumptions(heat conduction model;heat flux coupling model;discrete powder particle heat flux coupling model);(3)Applications of various microstructure evolution models[phase field(PF),cellular automata(CA),and Monte Carlo(MC)].Finally,the development trend of numerical simulation of metal additive manufacturing,including the thermal-flow-solid coupling model and deep learning for numerical model,is analyzed.

Numerical simulation of aluminum holding furnace with fluid-solid coupled heat transfer

To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.

NUMERICAL SIMULATION OF HEAT,MASS AND MOMENTUM TRANSFERS IN METALLIC INGOTS AND PREDICTION OF MACROSEGREGATION

A continuum model and numerical methods were established for description of heat,mass and momentum transfers as well as macrosegregation formations in metallic ingots.Numerical simulation of temperature,composition and liquid flow fields during the solidification of an Al-4.5% Cu ingot were performed on an IBM personal computer.The model and numerical methods were verified through two experiments.

A COMPUTER CODE FASTOR-3D FOR TRANSIENT THREE-DIMENSIONAL SIMULATION ON FLOW AND HEAT TRANS-FER WITH POROSITY APPROXIMATION

ＡＣＯＭＰＵＴＥＲＣＯＤＥＦＡＳＴＯＲ┐３ＤＦＯＲＴＲＡＮＳＩＥＮＴＴＨＲＥＥ┐ＤＩＭＥＮＳＩＯＮＡＬＳＩＭＵＬＡＴＩＯＮＯＮＦＬＯＷＡＮＤＨＥＡＴＴＲＡＮＳ┐ＦＥＲＷＩＴＨＰＯＲＯＳＩＴＹＡＰＰＲＯＸＩＭＡＴＩＯＮ＊ＬｕＷａｎｃｈｅｎｇ（卢万成）Ｄ...

Investigations on Structural Safety of Office Room Based on Fire Simulation and Transient Heat Transfer Analysis

This study aims at investigating heat propagations inside the structural members due to fire using fire simulation and transient heat transfer analysis. Toward that goal, fire simulation and transient heat transfer analysis for a 5-story building are carried out sequentially using Fire Dynamics Simulator (FDS) and ABAQUS 6.10-3, respectively. As results from fire simulation, temporal temperature information is obtained depending on various locations of the building, which is used as boundary condition for the structural elements generated in transient heat transfer analysis. Predictions from the transient heat transfer analysis show that the structural members are exposed spatially non-uniform temperatures which can cause significant eccentric deformation and aceleration of structural damages.