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.

Controlling Roll Temperature by Fluid-Solid Coupled Heat Transfer

Currently, when magnesium alloy sheet is rolled, the method of controlling roll temperature is simple and inaccurate. Furthermore, roll temperature has a large influence on the quality of magnesium alloy sheet; therefore, a new model using circular fluid flow control roll temperature has been designed. A fluid heat transfer structure was designed, the heat transfer process model of the fluid heating roll was simplified, and the finite di erence method was used to cal?culate the heat transfer process. Fluent software was used to simulate the fluid?solid coupling heat transfer, and both the trend and regularity of the temperature field in the heat transfer process were identified. The results show that the heating e ciency was much higher than traditional heating methods(when the fluid heat of the roll and tempera?ture distribution of the roll surface was more uniform). Moreover, there was a bigger temperature di erence between the input and the output, and after using reverse flow the temperature di erence decreased. The axial and circum?ferential temperature distributions along the sheet were uniform. Both theoretical calculation results and numerical simulation results of the heat transfer between fluid and roll were compared. The error was 1.8%–12.3%, showing that the theoretical model can both forecast and regulate the temperature of the roll(for magnesium alloy sheets) in the rolling process.

A Mathematical Model of Heat Transfer in Problems of Pipeline Plugging Agent Freezing Induced by Liquid Nitrogen

A mathematical model for one-dimensional heat transfer in pipelines undergoing freezing induced by liquid nitrogen is elaborated.The basic premise of this technology is that the content within a pipeline is frozen to form a plug or two plugs at a position upstream and downstream from a location where work a modification or a repair must be executed.Based on the variable separation method,the present model aims to solve the related coupled heat conduction and moving-boundary phase change problem.An experiment with a 219 mm long pipe,where water was taken as the plugging agent,is presented to demonstrate the relevance and reliability of the proposed model(results show that the error is within 18%).Thereafter,the model is applied to predict the cooling and freezing process of pipelines with different inner diameters at different liquid nitrogen refrigeration temperatures when water is used as the plugging agent.

Analysis of surface cracking in water-cooled rolls and heat transfer between furnace chamber and rolls in a direct flame furnace

In the direct fired furnace of a continuous annealing line, seal rolls are susceptible to deformation that leads to surface defects of steel strips. According to failure analysis, the reasons include improper structural design and heat imbalance. An improved design has been proposed to reduce stress concentration and thermal radiation. A heat transfer model has been employed to determine the proper water flow rate for roll cooling. Industrial application proves that seal rolls with the new design has less deformation and longer service life.

Impact of Separator Thickness on Temperature Distribution in Single Polymer Electrolyte Fuel Cell Based on 1D Heat Transfer

It is known from the New Energy and Industry Technology Development Organization (NEDO) roam map Japan, 2017 that the polymer electrolyte fuel cell (PEFC) power generation system is required to operate at 100°C for application of mobility usage from 2020 to 2025. This study aims to clarify the effect of separator thickness on the distribution of the temperature of reaction surface (Treact) at the initial temperature of cell (Tini) with flow rate, relative humidity (RH) of supply gases as well as RH of air surrounding cell of PEFC. The distribution of Treact is estimated by means of the heat transfer model considering the H2O vapor transfer proposed by the authors. The relationship between the standard deviation of Treact-Tini and total voltage obtained in the experiment is also investigated. We can know the effect of the flow rate of supply gas as well as RH of air surrounding cell of PEFC on the distribution of Treact-Tini is not significant. It is observed the wider distribution of Treact-Tini provides the reduction in power generation performance irrespective of separator thickness. In the case of separator thickness of 1.0 mm, the standard deviation of Treact-Tini has smaller distribution range and the total voltage shows a larger variation compared to the other cases.

Analysis of coupled flow-reaction with heat transfer in heap bioleaching processes

A mathematical model for heap bioleaching is developed to analyze heat transfer, oxygen flow, target ion distribution and oxidation leaching rate in the heap. The model equations are solved with Comsol Multiphysics software. Numerical simulation results show the following facts： Concentration of oxygen is relatively high along the boundary of the slope, and low in the center part where leaching rate is slow. Temper- ature is relatively low along the slope and reaches the highest along the bottom region near the slope, with difference being more than 6℃. Concentration of target mental ions is the highest in the bottom region near the slope. Oxidation leaching rate is relatively large in the bottom and slope part with a fast reaction rate, and small in the other part with low oxygen concentration.

Modeling of Heat Transfer and Solidification of Composite Roll

Modeling of heat transfer and solidification of composite roll was established and used to predict the thermal history and solidification process of roll during spray forming. Evolution of temperature field of the preform and cooling rate in the growing deposit during spray deposition and post-deposition were numerically simulated.

CFD-Based Optimization of a Shell-and-Tube Heat Exchanger

The main objective of this study is the technical optimization of a Shell-and-Tube Heat Exchanger(STHE).In order to do so,a simulation model is introduced that takes into account the related gas-phase circulation.Then,simulation verification experiments are designed in order to validate the model.The results show that the tem-peraturefield undergoes strong variations in time when an inlet wind speed of 6 m/s is considered,while the heat transfer error reaches a minimum of 5.1%.For an inlet velocity of 9 m/s,the heat transfer drops to the lowest point,while the heat transfer error reaches a maximum,i.e.,9.87%.The pressure drop increasesfirst and then decreases with an increase in the wind speed and reaches a maximum of 819 Pa under the 9 m/s wind speed con-dition.Moreover,the pressure drops,and the heat transfer coefficient increases with the Reynolds number.

The fast method and convergence analysis of the fractional magnetohydrodynamic coupled flow and heat transfer model for the generalized second-grade fluid

In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.

CFD study of non-premixed swirling burners: Effect of turbulence models

This research investigates a numerical simulation of swirling turbulent non-premixed combustion.The effects on the combustion characteristics are examined with three turbulence models:namely as the Reynolds stress model,spectral turbulence analysis and Re-Normalization Group.In addition,the P-1 and discrete ordinate(DO)models are used to simulate the radiative heat transfer in this model.The governing equations associated with the required boundary conditions are solved using the numerical model.The accuracy of this model is validated with the published experimental data and the comparison elucidates that there is a reasonable agreement between the obtained values from this model and the corresponding experimental quantities.Among different models proposed in this research,the Reynolds stress model with the Probability Density Function(PDF)approach is more accurate(nearly up to 50%)than other turbulent models for a swirling flow field.Regarding the effect of radiative heat transfer model,it is observed that the discrete ordinate model is more precise than the P-1 model in anticipating the experimental behavior.This model is able to simulate the subcritical nature of the isothermal flow as well as the size and shape of the internal recirculation induced by the swirl due to combustion.

Temperature field model in surface grinding: a comparative assessment

Grinding is a crucial process in machining workpieces because it plays a vital role in achieving the desired precision and surface quality.However,a significant technical challenge in grinding is the potential increase in temperature due to high specific energy,which can lead to surface thermal damage.Therefore,ensuring control over the surface integrity of workpieces during grinding becomes a critical concern.This necessitates the development of temperature field models that consider various parameters,such as workpiece materials,grinding wheels,grinding parameters,cooling methods,and media,to guide industrial production.This study thoroughly analyzes and summarizes grinding temperature field models.First,the theory of the grinding temperature field is investigated,classifying it into traditional models based on a continuous belt heat source and those based on a discrete heat source,depending on whether the heat source is uniform and continuous.Through this examination,a more accurate grinding temperature model that closely aligns with practical grinding conditions is derived.Subsequently,various grinding thermal models are summarized,including models for the heat source distribution,energy distribution proportional coefficient,and convective heat transfer coefficient.Through comprehensive research,the most widely recognized,utilized,and accurate model for each category is identified.The application of these grinding thermal models is reviewed,shedding light on the governing laws that dictate the influence of the heat source distribution,heat distribution,and convective heat transfer in the grinding arc zone on the grinding temperature field.Finally,considering the current issues in the field of grinding temperature,potential future research directions are proposed.The aim of this study is to provide theoretical guidance and technical support for predicting workpiece temperature and improving surface integrity.

Conjugate heat transfer investigations of turbine vane based on transition models

The accurate simulation of boundary layer transition process plays a very important role in the prediction of turbine blade temperature field. Based on the Abu-Ghannam and Shaw （AGS） and c-Re h transition models, a 3D conjugate heat transfer solver is developed, where the fluid domain is discretized by multi-block structured grids, and the solid domain is discretized by unstructured grids. At the unmatched fluid/solid interface, the shape function interpolation method is adopted to ensure the conservation of the interfacial heat flux. Then the shear stress transport （SST） model, SST ＆ AGS model and SST ＆ c-Re h model are used to investigate the flow and heat transfer characteristics of Mark II turbine vane. The results indicate that compared with the full turbulence model （SST model）, the transition models could improve the prediction accuracy of temperature and heat transfer coefficient at the laminar zone near the blade leading edge. Compared with the AGS transition model, the c-Re h model could predict the transition onset location induced by shock/boundary layer interaction more accurately, and the prediction accuracy of temperature field could be greatly improved.

Research on properties of hollow glass microspheres/epoxy resin composites applied in deep rock in-situ temperature-preserved coring

Deep petroleum resources are in a high-temperature environment.However,the traditional deep rock coring method has no temperature preserved measures and ignores the effect of temperature on rock porosity and permeability,which will lead to the distortion of the petroleum resources reserves assessment.Therefore,the hollow glass microspheres/epoxy resin(HGM/EP)composites were innovatively proposed as temperature preserved materials for in-situ temperature-preserved coring(ITP-Coring),and the physical,mechanical,and temperature preserved properties were evaluated.The results indicated that:As the HGM content increased,the density and mechanical properties of the composites gradually decreased,while the water absorption was deficient without hydrostatic pressure.For composites with 50 vol%HGM,when the hydrostatic pressure reached 60 MPa,the water absorption was above 30.19%,and the physical and mechanical properties of composites were weakened.When the hydrostatic pressure was lower than 40 MPa,the mechanical properties and thermal conductivity of composites were almost unchanged.Therefore,the composites with 50 vol%HGM can be used for ITPCoring operations in deep environments with the highest hydrostatic pressure of 40 MPa.Finally,to further understand the temperature preserved performance of composites in practical applications,the temperature preserved properties were measured.An unsteady-state heat transfer model was established based on the test results,then the theoretical change of the core temperature during the coring process was obtained.The above tests results can provide a research basis for deep rock in-situ temperature preserved corer and support accurate assessment of deep petroleum reserves.

Shape of slab solidification end under non-uniform cooling and its influence on the central segregation with mechanical soft reduction

In order to study the effect of continuous casting process parameters on the shape of slab solidification end under non-uniform cooling,a solidification model of a continuous-cast slab with non-uniform cooling condition was established with ProCAST software.The model was verified by the results of nail shooting tests and the infrared temperature measurement equipment.Four characteristic parameters were defined to evaluate the uniformity of the shape of slab solidification end.The results showed that the nonuniformity at the beginning and end of solidification,the solidification end length,and the solidification unevenness increased with the rise of casting speed.For each 10°C increase of superheat,the solidification unevenness increased by about 0.022.However,the effect of superheat on the solidification end length can be ignored.The secondary cooling strength showed minimal effect on the nonuniformity at the beginning and end of solidification.With the increase in secondary cooling intensity,the solidification end length decreased,but the solidification unevenness increased.In addition,the central segregation of the slab produced with and without the mechanical soft reduction(MSR)process was investigated.The transverse flow of molten steel with low solid fraction influenced the central segregation morphology under MSR.

Effect of External Forced Flow and Boiling Film on Heat Transfer of AISI 4140 Steel Horizontal Rod During Direct Quenching

The effects of rod falling and moving, external flow field, boiling film and radiation were investigated on fluid flow and heat transfer of AISI 4140 steel horizontal rod during direct quenching by mathematical modeling. The flow field and heat transfer in quenching tank were simulated by computational fluid dynamics （CFD） method considering falling and moving of rods during process. Therefore, modeling of flow field was done by a fixed-mesh method for general moving objects equations, and then, energy equation was solved with a numerical approach so that effeet of boiling film heat flux was considered as a source term in energy equation for solid-liquid boundary. Simulated results were verified by comparing with published and experimental data and there was a good agreement between them. Also, the effects of external forced flow and film boiling were investigated on heat flux output, temperature distribution and heat transfer coefficient of rod. Also simulated results determined optimum quenching time for this process.

Zonal Metamorphic Complex in the Tongulack Mountain Ridge, Altai, Russia, and Explanation ofIts Origin with the Help of Thermal Modelling

A moderate pressure/high temperature zonal metamorphic complex in the Tongulack Mountain Ridge, Altai, Russia, is described, and the applicability of the models of magmatic intrusion and fluid flow to explanation of its origin discussed. The Precambrian complex was formed at 500–700°C and 3.0–5.5 kbars; it is a linear, 25–30 km wide, thermal anticline with a curved axis showing symmetric metamorphic zoning. The metamorphism was isochemical by its nature, as is corroborated by the chemical compositions of the rocks. Four zones can be recognized within the metamorphic complex: chloritic (on the peripheries), cordieritic, sillimanitic and staurolite-out (in the centre). The zones are separated by successive isograds: cordierite, staurolite-in or sillimanite and staurolite-out. It is argued that the origin of the metamorphic zoning can be explained best by a combined fluid-magmatic model; conductive heat flow from the intrusion predominated considerably over the fluid flux in heat transfer: the fluid flow rate was estimated as about 3 ? 10?9 g/cm2, ? s. The modern position of the axial region of the metamorphic belt is predicted to be lying roughly about 1.5 km above the roof of the intrusive body.

Numerical and Experimental Investigation on the In-Flight Melting Behaviour of Granulated Powders in Induction Thermal Plasmas

An innovative in-flight glass melting technology with thermal plasmas was developed for the purpose of energy conservation and environment protection. In this study, modelling and experiments of argon-oxygen induction thermal plasmas were conducted to investigate the melting behaviour of granulated soda-lime glass powders injected into the plasma. A two-dimensional local thermodynamic equilibrium （LTE） model was performed to simulate the heat and momentum transfer between plasma and particle. Results showed that the particle temperature was strongly affected by the flow rate of carrier gas and the particle size of raw material. A higher flow rate of carrier gas led to lower particle temperature and less energy transferred to particles which resulted in lower vitrification. The incomplete melting of large particles was attributed to the lower central temperature of the particle caused by a larger heat capacity. The numerical analysis explained well the experimental results, which can provide valuable practical guidelines for the process control in the melting process for the glass industry.

Finite Element Modeling of a Solar Box Cooker: Temperature and Fluid Velocity Distributions

We studied the temperature distribution and fluid velocity in a box-type solar cooker by using the Finite Element Method (FEM) in Ziguinchor southern of Senegal. Indeed, this is one of the sunniest countries in the world: more than 3000 hours of sunshine per year with an average temperature of around 30˚C. This abundant and exploitable solar energy contributes to the development of more efficient, profitable and clean sources of energy. This will help to satisfy the increasing demand of energy. This numerical model was validated by comparing the numerical results with those of the experiment carried out on a single day. The relative error obtained is below 3%. The model results confirmed the performance of this cooker as its cooking temperature is available for more than seven hours. They have shown that the temperature and internal fluid velocity fields are not homogeneous. The results, although preliminary and encouraging, are a first step towards the complete simulation of a solar cooker integrated into a drying column.

Preliminary design method for absorber pipe length of tunnel lining ground heat exchangers based on energy efficiency of heat pump

Many ongoing tunnel projects provide a favorable opportunity for the investigation and application of tunnel lining ground heat exchangers(GHEs).Tunnel lining GHEs can be connected to a heat pump to extract geothermal energy for heating and cooling buildings.Numerous studies have focused on the thermal performance of tunnel lining GHEs;however,the studies on the interaction between heat pumps and tunnel lining GHEs are relatively rare.In this study,a coupled heat transfer model of heat pumps and tunnel lining GHEs was proposed and then calibrated based on in situ test results.The model was used to evaluate the energy efficiency of a heat pump with tunnel lining GHEs under different conditions.The results show that the energy efficiency ratio(EER)increases exponentially with the absorber pipe length and thermal conductivity of the surrounding rock.The EER is governed by the convection heat transfer coefficient,which varies exponentially;meanwhile,the EER decreases approximately linearly with the annual average air temperature in the tunnel.Different types of heat pumps affect the EER significantly,and the EER of a Type-3 heat pump is higher than that of a Type-1 heat pump by 27.1%.Based on the aforementioned results,an empirical formula for the EER and absorber pipe length was established.Moreover,a preliminary design method for the absorber pipe length based on this empirical formula was developed.The method was employed to determine the appropriate absorber pipe length for the tunnel lining GHEs in the Shapu tunnel in Shenzhen,China.Finally,groups of absorber pipe layouts with a pipe spacing of 0.5 m,area of 135 m2,and length of 293.5 m were preliminarily determined.

Transient thermal response of micro-thermal conductivity detector (μTCD) for the identification of gas mixtures:An ultra-fast and low power method

Micro-thermal conductivity detector(μTCD)gas sensors work by detecting changes in the thermal conductivity of the surrounding medium and are used as detectors in many applications such as gas chromatography systems.Conventional TCDs use steady-state resistance(i.e.,temperature)measurements of a micro-heater.In this work,we developed a new measurement method and hardware configuration based on the processing of the transient response of a low thermal mass TCD to an electric current step.The method was implemented for a 100-μm-long and 1-μm-thick micro-fabricated bridge that consisted of doped polysilicon conductive film passivated with a 200-nm silicon nitride layer.Transient resistance variations of theμTCD in response to a square current pulse were studied in multiple mixtures of dilute gases in nitrogen.Simulations and experimental results are presented and compared for the time resolved and steady-state regime of the sensor response.Thermal analysis and simulation show that the sensor response is exponential in the transient state,that the time constant of this exponential variation was a linear function of the thermal conductivity of the gas ambient,and that the sensor was able to quantify the mixture composition.The level of detection in nitrogen was estimated to be from 25 ppm for helium to 178 ppm for carbon dioxide.With this novel approach,the sensor requires approximately 3.6 nJ for a single measurement and needs only 300μs of sampling time.This is less than the energy and time required for steady-state DC measurements.