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.

Extended Generalized Riccati Equation Mapping for Thermal Traveling-Wave Distribution in Biological Tissues through a Bio-Heat Transfer Model with Linear/Quadratic Temperature-Dependent Blood Perfusion

Analytical thermal traveling-wave distribution in biological tissues through a bio-heat transfer (BHT) model with linear/quadratic temperature-dependent blood perfusion is discussed in this paper. Using the extended generalized Riccati equation mapping method, we find analytical traveling wave solutions of the considered BHT equation. All the travelling wave solutions obtained have been used to explicitly investigate the effect of linear and quadratic coefficients of temperature dependence on temperature distribution in tissues. We found that the parameter of the nonlinear superposition formula for Riccati can be used to control the temperature of living tissues. Our results prove that the extended generalized Riccati equation mapping method is a powerful tool for investigating thermal traveling-wave distribution in biological tissues.

An Improved Heat Transfer Model for Flame Resistant Fabric Used in Thermal Protective Clothing

An improved numerical heat transfer model considering pyrolysis effect is proposed to predict thermal performance of heat-resistant fabric subjected to radiant heat flux. The model incorporates the heat-induced changes in fabric thermophysical properties. The new model has been validated with data from modified Radiant Protective Performance (RPP) tests of flame-resistant cotton fabrics. Comparison with experimental data shows that the predictions of mass loss rates and temperature profiles within the charring material and skin simulant are in reasonably good agreement with the experiments. Results from the numerical model contribute to a better understanding of the heat transfer process within flame-resistant fabrics under high heat flux conditions, and also to establish a systematic method for analyzing heat transfer in other fibrous materials applications.

Surface-particle-emulsion heat transfer model between fluidized bed and horizontal immersed tube

A mathematical model, surface-particle-emulsion heat transfer model, is presented by considering voidage variance in emulsion in the vicinity of an immersed surface. Heat transfer near the surface is treated by dispersed particles touching the surface and through the emulsion when the distance from the surface is greater than the diameter of a particle. A film with an adjustable thickness which separates particles from the surface is not introduced in this model. The coverage ratio of particles on the surface is calculated by a stochastic model of particle packing density on a surface. By comparison of theoretical solutions with experimental data from some References, the mathematical model shows better qualitative and quantitative prediction for local heat transfer coefficients around a horizontal immersed tube in a fluidized bed.

A Revised CFB Wall-to-suspension Heat Transfer Model

Based on the Cluster Renewal Model of the particle motion in a CFB riser, a revised heat transfer model is developed, which introduces the latest research results of the hydrodynamics of the suspension flow in CFB. This model divides the heat transfer into two parts, which are due to the transient heat conduction by the covered clusters and the convection between the uncovered wall and the dispersed phase. Radiation at high temperature is regarded as being additive. The fraction of the covered wall by clusters is revised by a new formula, which is a function of the operating condition and the particle properties. The radiation between the dispersed phase and the uncovered wall includes not only the direct radiation to the uncovered wall, but also the radiation to the clusters and then reflected to the uncovered wall. Calculation was carried out for the CFB heat transfer model. The results were compared with the published typical experimental data of other researchers and showed a good agreement between them.

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.

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.

A soil water and heat transfer model including changes in soil frost and thaw fronts

Freeze-thaw processes in soils,including changes in frost and thaw fronts(FTFs),are important physical processes.The movement of FTFs affects soil hydrothermal characteristics,as well as energy and water exchanges between the land surface and the atmosphere and hydrothermal processes in the land surface.This paper reduces the issue of soil freezing and thawing to a multiple moving-boundary problem and develops a soil water and heat transfer model which considers the effects of FTF on soil hydrothermal processes.A local adaptive variable-grid method is used to discretize the model.Sensitivity tests based on the hierarchical structure of the Community Land Model(CLM)show that multiple FTFs can be continuously tracked,which overcomes the difficulties of isotherms that cannot simultaneously simulate multiple FTFs in the same soil layer.The local adaptive variable-grid method is stable and offers computational efficiency several times greater than the high-resolution case.The simulated FTF depths,soil temperatures,and soil moisture values fit well with the observed data,which further demonstrates the potential application of this simulation to the land-surface process model.

Modelling Study to Compare the Flow and Heat Transfer Characteristics of Low-Power Hydrogen,Nitrogen and Argon Arc-Heated Thrusters

A modelling study is performed to compare the plasma flow and heat transfer characteristicsof low-power arc-heated thrusters (arcjets) for three different propellants: hydrogen,nitrogen and argon.The all-speed SIMPLE algorithm is employed to solve the governing equations,which take into account the effects of compressibility,Lorentz force and Joule heating,aswell as the temperature- and pressure-dependence of the gas properties.The temperature,velocityand Mach number distributions calculated within the thruster nozzle obtained with differentpropellant gases are compared for the same thruster structure,dimensions,inlet-gas stagnantpressure and arc currents.The temperature distributions in the solid region of the anode-nozzlewall are also given.It is found that the flow and energy conversion processes in the thrusternozzle show many similar features for all three propellants.For example,the propellant is heatedmainly in the near-cathode and constrictor region,with the highest plasma temperature appearingnear the cathode tip; the flow transition from the subsonic to supersonic regime occurs withinthe constrictor region; the highest axial velocity appears inside the nozzle; and most of the inputpropellant flows towards the thruster exit through the cooler gas region near the anode-nozzlewall.However,since the properties of hydrogen,nitrogen and argon,especially their molecularweights,specific enthalpies and thermal conductivities,are different,there are appreciable differencesin arcjet performance.For example,compared to the other two propellants,the hydrogenarcjet thruster shows a higher plasma temperature in the arc region,and higher axial velocitybut lower temperature at the thruster exit.Correspondingly,the hydrogen arcjet thruster has thehighest specific impulse and arc voltage for the same inlet stagnant pressure and arc current.Thepredictions of the modelling are compared favourably with available experimental results.

Optimising Gas Quenching Technology through Modelling of Heat Transfer

Gas Quenching represents an environmentally friendly alternative to more commonly-used oil quenching. Yet, the performances of this technology remain limited in terms of cooling rates reached compared to oil quenching. Distortion and process homogeneity also have to be controlled carefully. The efficiency of the gas quenching process fully depends on the heat transfer between the gas and the quenched parts. The goal of this study is the optimisation of the gas quenching process efficiency through a better understanding of the heat transfer phenomena involved. The study has been performed with modelling means and validated by an experimental approach. The configuration of the gas flow has a major influence on the heat transfer phenomena between the gas and the parts. The fluid dynamics modelling approach performed in this study allows to optimise the heat transfer phenomena. New gas quenching processes allowing enhanced gas quenching performance through higher cooling rates can be thereby identified. The new solutions have been validated in experimental and industrial conditions. Results obtained allow to expect significant improvement of high pressure gas quenching technology.

Heat Transfer Investigation and Modeling of Heat Integrated Distillation Column

The high degree of reversibility of heat integrated distillation column(HIDiC) has been thermodynamically interpreted by the entropy method. In this paper, a heat transfer model and a more universal method were proposed, through which the overall heat transfer coefficient at different height of column under different operating conditions could be obtained before the experiment. Then the separation of a binary ethanol-water system was carried out experimentally as a case study to verify the heat transfer model and the aforementioned calculation method. The close results between the calculation, the simulation, and the experiments suggested that the proposed model and the calculation method in this paper were accurate and applicable. Meanwhile, it was demonstrated that the HIDiC shows obvious effect of reducing entropy increase and improving thermodynamic efficiency as compared to conventional distillation column.

Modeling and Simulation of Heat Transfer in Loaded Continuous Heat Treatment Furnace

Continuous furnaces are widely used in the heat treatment of mass-produced parts. However, the heating up process of parts in continuous furnace is still decided by experience. In this paper the heat transfer in the continuous furnace is formulated firstly. The heat balance in each zone is discussed and equations are given. Coupled with the model for heat transfer between workpieces and furnace and the heat transfer in the workload as well presented in the former developed CHT-6/ for batch furnaces, a program CHT- for continuous furnaces was developed. The model deals with two typical movements of parts: continuous or step by step. The moving speed of parts and load pattern can be optimized based on the calculated temperature distributions and curves, especially, the fastest heated and slowest-heated temperature-distance profiles. A case study is carried out for the heat treatment of a kind of hook-shaped part. The calculated results are analyzed and in good agreement with the measured ones.

Fluid Flow and Heat Transfer Modeling of AC Arc in Ferrosilicon Submerged Arc Furnace

A two-dimensional mathematical model was developed to describe the heat transfer and fluid flow in an AC arc zone of a ferrosilicon submerged arc furnace.In this model,the time-dependent conservation equations of mass,momentum,and energy in the specified domain of plasma zone were numerically solved by coupling with the Maxwell and Laplace equations for magnetic filed and electric potential,respectively.A control volume-based finite difference method was used to solve the governing equations in cylindrical coordinates.The reliability of the developed model was checked by experimental data from the previous available literature.The results of present model were in good agreement with the given data comparing with other models,because of solving the Maxwell and Laplace equations simultaneously in order to calculate current density.In addition,parametric studies were carried out to evaluate the effects of electrical current and arc length on flow field and temperature distribution within the arc.According to the computed results,a lower power input led to a higher arc efficiency.

A coupled model on fluid flow,heat transfer and solidification in continuous casting mold

Fluid flow, heat transfer and solidification of steel in the mold are so complex but crucial, determining the surface quality of the continuous casting slab. In the current study, a 2D numerical model was established by Fluent software to simulate the fluid flow, heat transfer and solidification of the steel in the mold. The VOF model and k-ε model were applied to simulate the flow field of the three phases(steel, slag and air), and solidification model was used to simulate the solidification process. The phenomena at the meniscus were also explored through interfacial tension between the liquid steel and slag as well as the mold oscillation. The model included a 20 mm thick mold to clarify the heat transfer and the temperature distribution of the mold. The simulation results show that the liquid steel flows as upper backflow and lower backflow in the mold, and that a small circulation forms at the meniscus. The liquid slag flows away from the corner at the meniscus or infiltrates into the gap between the mold and the shell with the mold oscillating at the negative strip stage or at the positive strip stage. The simulated pitch and the depth of oscillation marks approximate to the theoretical pitch and measured depth on the slab.

Modelling Study on the Plasma Flow and Heat Transfer in a Laminar Arc Plasma Torch Operating at Atmospheric and Reduced Pressure

A modelling study is performed to investigate the characteristics of both plasma flowand heat transfer of a laminar non-transferred arc argon plasma torch operated at atmosphericand reduced pressure.It is found that the calculated flow fields and temperature distributions arequite similar for both cases at a chamber pressure of 1.0 atm and 0.1 atm.A fully developed flowregime could be achieved in the arc constrictor-tube between the cathode and the anode of theplasma torch at 1.0 atm for all the flow rates covered in this study.However the flow field couldnot reach the fully developed regime at 0.1 atm with a higher flow rate.The arc-root is alwaysattached to the torch anode surface near the upstream end of the anode,i.e.the abruptly expandedpart of the torch channel,which is in consistence with experimental observation.The surroundinggas would be entrained from the torch exit into the torch interior due to a comparatively largeinner diameter of the anode channel compared to that of the arc constrictor-tube.

Establishment and Application of UFC-ACC Heat Transfer Coefficient Model

Based on medium plate runout table ultra-fast cooling( UFC)-accelerated cooling equipment( ACC) system,a heat transfer coefficient model was constructed. Firstly,according to the measured data,heat transfer coefficients under different roll speed and water volume were calculated by using an inverse heat conduction method. Secondly,a monofactorial heat transfer coefficient calculation formula was obtained. Finally,the heat transfer coefficient model based on medium plate runout table UFC-ACC system was constructed by intercept function,slope function,interaction influence function and linear or nonlinear influencing factors. The precision of these models was validated by comparing model prediction value with measured data,and the results were in good agreement with practical needs,and the average deviation was less than 5%.

Coupled Transfer of Water and Heat in Red Soil: Experiment and Numerical Modelling

Coupled transfer of soil water and heat in closed columns of homogeneous red soil was studied under laboratory conditions. A coupled model was constructed using soil physical theory, empirical equations and experimental data to predict the coupled transfer. The results show that transport of soil water was affected by temperature gradient, and the largest net water transport was found in the soil column with initial water content of 0.148 m3 m-3. At the same time, temperature changes with the transport of soil water was in a nonlinear shape as heat parameters were function of water content, and the changes of temperature were positively correlated with the net amount of water transported. Numerical modelling results show that the predicted values of temperature distribution were close to the observed values, while the predicted values of water content exhibited limited deviation at both ends of the soil column due to the slight temperature changes at both ends. It was indicated that the model proposed here was applicable.

A Mathematical Model to Solve Bio-Heat Transfer Problems through a Bio-Heat Transfer Equation with Quadratic Temperature-Dependent Blood Perfusion under a Constant Spatial Heating on Skin Surface

We consider the one-dimensional bio-heat transfer equation with quadratic temperature-dependent blood perfusion, which governs the temperature distribution inside biological tissues. Using an extended mapping method with symbolic computation, we obtain the exact analytical thermal traveling wave solution, which describes the non-uniform temperature distribution inside the bodies. The found exact solution is used to investigate the temperature distribution in the tissues. It is found that the surrounding medium with higher temperature does not necessarily imply that the tissue will quickly (after a short duration of heating process) reach the desired temperature. It is also found that increased perfusion causes a decline in local temperature.

Modelling of heat transfer for progressive freeze concentration process by spiral finned crystallizer

This study presents a novel design for a spiral finned crystallizer which is the primary element of progressive freeze concentration(PFC) system, which simplifies the setup of the conventional system. After the crystallizer has been designed, the research experiments have been conducted and evaluated through a thorough analysis of its performance by developing a mathematical model that can be used to predict the productivity of ice crystal at a range of coolant temperature. The model is developed based on the basic heat transfer equation, and by considering the solution's and the coolant's convective heat transfer coefficient(h) under the forced flow condition.The model's accuracy is verified by making comparison between the ice crystal mass' experimental value and the values predicted by the model. Consequently, the study found that the model helps in enhancing the PFC system.