NUMERICAL INVESTIGATION OF MASS AND HEAT TRANSFER IN SODIUM POOL COMBUSTION

Liquid sodium pool combustion is a coupled phenomenon of heat and mass transfer, chemical reaction, and aerosol dynamics. This article presents numerical investigation of the coupled phenomena using a new computational tool. Boundary layer equations are solved to obtain the thermalhydraulic field above the pool surface. The chemical reaction and heat and mass transfer are solved interactively considering radiation heat transfer and behavior of reaction product aerosols. Numerical results are compared with experiments, and agreement is excellent concerning burning rate, flame temperature, and height. It is concluded that the sodium pool combustion is self-limited and negative feedback is at work.     

质子交换膜燃料电池内部传热传质三维模拟

针对常规流场质子交换膜燃料电池提出了三维非等温数学模型。模型考虑了电化学反应动力学以及反应气体在流道和多孔介质内的流动和传递过程,详细研究了水在质子膜内的电渗和扩散作用。计算结果表明,反应气体传质的限制和质子膜内的水含量直接决定了电极局部电流密度的分布和电池输出性能;在电流密度大于0.3~0.4A/cm2时开始出现水从阳极到阴极侧的净迁移;高电流密度时膜厚度方向存在很大的温度梯度,这对膜内传递过程有较大影响。     

NUMERICAL ANALYSIS OF FLUID FLOW AND HEAT TRANSFER CHARACTERISTICS IN THREE-DIMENSIONAL PLATE FIN-AND-TUBE HEAT EXCHANGERS

This numerical study provides three-dimensional (3-D) time-dependent modeling of unsteady laminar flow and heat transfer over single- and multirow plate fin-and-tube heat exchangers. The complex nature of the flow field featuring a horseshoe vortex is investigated for both configurations. The time-dependent evolution of the horseshoe vortex mechanism on the forward part of the tube and its journey to the rear of the tube are studied to provide fundamental information on the local flow structure and the corresponding heat transfer characteristics. The effects of various governing parameters, such as fin spacing, Reynolds number, tube row number, and tube arrangement, on the heat transfer and flow characteristics are also studied for the Reynolds number range investigated. It is found that the local flow structure including formation and evolution of vortex systems and singular-point interactions correlates strongly with the heat transfer characteristics. The numerical results for the integral heat transfer parameters agree well with available experimental measurements.     

THREE-DIMENSIONAL MODELING OF TRANSIENT HEAT TRANSFER ON A SOAKING PIT

A numerical study has been conducted for heat transfer in a soaking pit. A transient three-dimensional radiation-convection-conduction model is presented. Parametric simulations are performed to assess the effect of various operating parameters on the heat transfer to the ingots. Results have shown that, under typical operating conditions, forced convection and radiation account for approximately 75% and 25% of the total heat transfer to the ingots, respectively. Heat transfer is maximized for a specific distance between the ingots. The average heat transfer coefficient to the ingots is strongly influenced by the gas mass flow rate but is nearly independent of the inlet gas temperature and of the ingot length.     

CHARACTERISTICS OF THREE-DIMENSIONAL FLOW,HEAT, AND MASS TRANSFER IN A CHEMICAL VAPOR DEPOSITION REACTOR

The three-dimensional flow, heat, and mass transfer characteristics in a horizontal chemical vapor deposition (CVD) reactor with a tilted susceptor are analyzed numerically. As the physical domain for the CVD reactor has an irregular region, the Cartesian coordinates are transformed into a general curvilinear coordinate system. For calculating the governing equations, the finite volume method (FVM) is adopted and the SIMPLE algorithm is extended and modified to employ the curvilinear system. The effects of flow rate and tilted angle of the susceptor on the transport phenomena (i.e., heat transfer rate, uniformity and growth rate of reactant gas, etc.) at the susceptor are investigated. The results show that the existence of return flows leads locally to improvement of the heat transfer, but it is not good for the uniformity of the deposition. As the tilted angle of the susceptor is increased (from 0^o to 9^o), the amount of heat transfer and growth rate in the main flow direction are improved irrespective of the Reynolds number. The growth rate and uniformity are less sensitive to the inclined angle of the susceptor than that of the Reynolds number.     

除湿条件下波纹通道内对流传热传质的数值研究

通过数值计算对紧凑换热器一种波纹翅片通道内除湿条件下周期性充分发展的对流传热传质情况进行数值研究。计算采用曲线坐标系下压力与速度耦合的SIMPLER算法,湿空气流动Re数的范围为100～1100,Pr数为0.71,Sc数为0.61。讨论了不同波纹高度、波纹间距对阻力与换热的影响,给出了不同Re数下的浓度场,并对动量、传热及传质进行了定量比较分析。计算结果表明,整体Nu数及fRe数随着波纹高度的增加或波纹间距的减小而增加;浓度随着Re数的增加沿着流动方向迅速降低;计算能较好的满足Chilton-Colburn相似,表明传热特性均可类推到传质特性中去。     

NUMERICAL SIMULATION OF THREE-DIMENSIONAL SEPARATED FLOW AND HEAT TRANSFER AROUND STAGGERED SURFACE-MOUNTED RECTANGULAR BLOCKS IN A CHANNEL

ABSTRACT

Numerical results simulating a three-dimensional laminar separated flow and heat transfer around staggered surface-mounted rectangular blocks in a plane channel are presented. Treated in the present study is a case of staggered three-row blocks. The finite-difference method is employed to solve the Navier-Stokes and energy Equations directly, and the resulting finite-difference Equations are solved with the SMAC method for Re = 100–500 and Pr = 0.7. The present numerical results are found to simulate well the visualization results such as horseshoe vortices and recirculating flow. The heat transfer coefficient greatly varies on the different side surfaces of blocks and also with Reynolds number.     

FINITE ELEMENT MODELING OF THREE-DIMENSIONAL TRANSIENT HEAT TRANSFER IN STAINLESS STEEL (304) PULSED GTA WELDMENTS

Pulsed current gas tungsten arc welding involves a large number of process parameters, and hence makes it difficult, if not impractical, to develop the procedure for a given application. A numerical model such as the finite element model will be of great use to study the effect of process parameters on heat transfer and attendant characteristics. A procedure for development of a three-dimensional transient nonlinear finite element model is presented in this article. Sensitivity analysis has been carried out to check the response of the model to constant and temperature-dependent material properties, heat loss due to vaporization of alloying elements, and total number of nodes to model the solution domain. Experimental work on thermal profiles and weld bead dimensions has been carried out to validate the results predicted by the model. Finally, parametric studies have been carried out to study the combined effects of parameters that influence the heat input.     

MODELING SIMULTANEOUS HEAT AND MASS TRANSFER USING THE CONTROL-VOLUME METHOD

A control-volume formulation for the solution of a set of two-way coupled heat and diffusive moisture transfer equations is presented in three dimensions. The three-dimensional model in Cartesian coordinates was extended to other geometries using an equivalent surface area-to-volume ratio index. The solution procedure developed uses a fully implicit time-stepping scheme for the solution of the coupled set of equations to study the drying behavior of barley and starch food. Simulated results of the one-way and the two-way coupled equations were compared with the finite-element analysis (two-way coupled and three-way coupled equation set) and experimental results from the literature. Among the products tested, the maximum deviations between experimental and simulation results were 5.0°C in temperature for barley drying, and 0.2 % moisture for starch drying. The overall predictions agreed well with the available experimental data and show good potential for application in grain and food drying.     

吸附式制冷系统传热传质过程的数值模拟

建立了沸石－水吸附式制冷系统吸附床在非第一类边界条件时的二维传热传质模型。考虑了工质对的吸附特性对脱附过程的影响,讨论了肋片数量、肋片物性参数,吸附床的有效导热系数,接触热阻等参数对脱附时间的影响。结果表明,增加肋片可以改善吸附床的性能,其中肋片的数量和接触热阻的影响较大,肋片的物性参数影响不大。     

太阳能蒸馏器中一组传热传质关系式

建立了盘式太阳能蒸馏器中一组传热传质关系式,并用多级迭盘式蒸馏室进行稳态模拟实验。发现该关系式在较大的雷利数范围及较宽的温度段内能很好地预示蒸馏器的产水量。     

热管式吸热器单元热管传热的数值模拟分析

热管式吸热器的热性能分析对吸热器设计有着重要意义，但由于其相变过程与热管传热的耦合作用十分复杂，至今仍是很少有人深入研究的领域。本文基于焓法建立单元热管耦合传热的物理和数学模型，模拟计算了热管壁温、蓄热容器壁温、循环工质出口温度及相变材料熔化率等参数，并与基本型吸热器进行比较，验证了热管吸热器明显改善了温度分布的均匀性和相变材料的熔化率。     

CONJUGATE HEAT AND MASS TRANSFER IN CONVECTIVE DRYING OF POROUS MEDIA

Abstract

A methodology for the analysis of conjugate problems in the convective drying of porous media is presented. In this study, the interface between porous medium and external convective flow is treated as an internal boundary within a two-phase system rather than a geometric limit. The problems of solid drying and convection boundary layer are connected by expressing the continuity of the state variables and their respective fluxes through the interface. The performance of the proposed methodology is evaluated by applying it to wood-drying problems. The analysis of the drying of porous media as a conjugate problem allows the assessment of the effect of the heat and mass transfer within the solid on the transfer in the adjacent fluid, providing good insight on the complexity of the transfer mechanisms.     

MESHLESS EFG METHOD IN THREE-DIMENSIONAL HEAT TRANSFER PROBLEMS: A NUMERICAL COMPARISON,COST AND ERROR ANALYSIS

The present analysis deals with the numerical solution of three-dimensional heat transfer problems using a meshless element-free Galerkin (EFG) method. The EFG method utilizes the moving least-square (MLS) approximation to approximate the unknown function of temperature T(x) with T ^h (x). The approximants are constructed by using a weight function, a monomial basis, and a set of coefficients that depends on position. A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced using the Lagrange multiplier method. MATLAB codes have been developed to obtain the numerical results for a model problem of three-dimensional heat transfer in orthotropic materials using different EFG weight functions. Three new weight functions, exponential, elliptical, and cosine, are proposed. A numerical comparison is made among the results obtained using proposed (exponential, elliptical, and cosine) and existing (quadratic) weight functions for a model problem. L₂ error is calculated for the proposed and existing EFG weight functions using 125 nodes. FORTRAN software has also been developed and executed on a PARAM 10000 supercomputing machine to obtain the computational cost of the EFG method. The computational cost of the EFG method is obtained for different orders of Gaussian quadrature and for different values of scaling parameter. The effect of scaling parameter on EFG results (temperature values) is also discussed in detail. The effectiveness of EFG method is shown by comparing the EFG results with those obtained by the finite-element method.     

U型管埋地换热器三维传热模型及实验对比分析

针对U型管埋地换热器的传热特点,考虑了两管脚之间的相互影响,建立起非稳态传热模型。采用数值解法对该理论模型进行了求解,并将结果与实验值进行了对比分析,发现两者变化趋势一致,但理论计算结果比实验测试值低,其中的一个重要原因是因为换热器与土壤间的热交换不是纯导热过程,还应该进一步考虑土壤中热湿迁移等因素的影响。     

NUMERICAL INVESTIGATION OF FLOW AND HEAT TRANSFER IN TWISTED CIRCULAR-SECTOR DUCTS

Flow and heat transfer in a twisted circular-sector duct are analyzed numerically for steady, fully developed, and incompressible laminar flow with a uniform-wall-temperature boundary condition. A rotating coordinate system is employed to account for the duct twist. The friction factors and Nusselt numbers are predicted for duct sector angles ranging from 15° to 90°, Reynolds numbers ranging from 1 to 1,000, and Prandtl numbers ranging from 1 to 100. Results show significant influence of duct twist on both friction factors and Nusselt numbers, particularly at large values of Reynolds and Prandtl numbers. Accurate correlations are developed to predict the friction factors and Nusselt numbers for the entire range of geometric and operating conditions studied.     

KRYLOV SUBSPACE SOLVERS IN PARALLEL NUMERICAL COMPUTATIONS OF PARTIAL DIFFERENTIAL EQUATIONS MODELING HEAT TRANSFER APPLICATIONS

In this study, parallel numerical algorithms for Krylov methods such as GMRES(k), Bi-CGM, Bi-CGSTAB, etc., for handling large-scale linear systems resulting from finite-difference analysis (FDA) and finite-element analysis (FEA) of coupled nonlinear partial differential equations (PDEs) describing problems in heat transfer applications are discussed. Parallel code has been successfully implemented on an eight-noded cluster under ANULIB message-passing library environment. Bi-CGM and ILU-GMRES(k) are found to give good performance for linear systems resulting from FEA, whereas Bi-CGSTAB is seen to give good performance with linear systems resulting from FDA.     

NUMERICAL INVESTIGATION OF CONVECTION HEAT TRANSFER IN A HEATED ROOM

We describe numerical investigation of airflow and temperature field in a room with a convective heat source. The simulation involves using computational fluid dynamics (CFD) to validate different turbulence models, i.e., the standard k- k model and the low Reynolds number k- k model. The comparisons between computations and experiments show good and acceptable agreement. It can be concluded that the CFD simulations can capture the main flow features and provide satisfactory results. It can be seen that the thermal wall jet created by the heat source greatly influences airflow pattern and temperature field in the room. It can also be seen that the advanced turbulence model may produce better results than the standard one under a suitable mesh scheme.     

NUMERICAL SIMULATION OF HEAT TRANSFER IN MIST FLOW

A numerical simulation of heat transfer over a row of tubes, in the presence of mist flow, is described. Computations include the solution of the flow field around the tubes, the prediction of the motion of water droplets, and the evaluation of the cooling effect of the water film on the tube surface. The entire analysis is carried out using FENSAP-ICE (Finite Element Navier-Stokes Analysis Package for In-flight icing), a simulation system developed by Newmerical Technologies for icing applications. The numerical model is described, including the Navier-Stokes solution, the water thin film computation, the droplet impingement prediction, and the conjugate heat transfer procedure. The predictions are verified against experimental data for different droplet mass flow rates, showing satisfactory agreement and allowing a useful insight in the physical characteristics of the problem.