Numerical analysis of three-dimensional flow,heat and mass transfer in arbitrary-shaped ducts

In this study, 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 SIMPLE algorithm is extended and modified to employ the curvilinear system. From the above modeling, the effects of flow rate and tilt angle of suscept on the uniformity and growth rate of reactant gas at the susceptor are investigated. The results show that the existence of return flows leads locally to the improvement of heat transfer, but it is not good for the uniformity of the susceptor. As the tilt angle of the susceptor is increased (from 0° to 9°), the amount of heat transfer and growth rate are improved irrespective of Reynolds number.© 1999 Scripta Technica, Heat Trans Asian Res, 28(6): 474–483, 1999     

Numerical analysis of heat flow in contact heat transfer

A finite difference analysis of heat conduction problem in a cylinder terminating in a frustum of a cone is presented. The constriction can be either in vacuum or in a gaseous environment. A fine mesh of 2500 × 800 was used for the construction of the grid such that very small constrictions could be analysed sufficiently accurately. Small constrictions i.e., small contact areas separated by large voids filled with a gas are typical of most practical applications involving contact heat transfer. The result of the finite difference analysis shows that gap conductance is predominant for all the gases considered. Gap-to-solid conductance ratio increases as the cone angle decreases due to the decrease of gap thickness. It also indicates that increase of conductance ratio is less significant at higher constriction angles. Finally, predicted conductance parameters are compared with the experimental results for different interfacial gases and a very good agreement is obtained.     

Numerical analysis for flow,heat and mass transfer in film flow along a vertical fluted tube

We conducted a three-dimensional numerical investigation of the flow, heat and mass transfer characteristics of the fluted evaporating tube with films flowing down on both the inside and outside tube walls. Condensation occurs along the outside wall while evaporation takes place on the free surface of the inside film. The three-dimensional transport equations for momentum and energy were solved by using the finite volume method (FVM). The free-surface shape is tracked by using the moving-grid technique that satisfies the space conservation law (SCL). Because of the secondary motion of the fluid, the film becomes thin at the crest whereas it thickens at the valley. The velocity and temperature fields were successfully predicted for various flute shapes.     

Numerical analysis on flow pattern and heat transfer characteristics of flow boiling in the mini-channels

Abstract

Flow pattern and heat transfer of flow boiling for different flow orientation, mini-channel width and height were presented in this work. The data were obtained by the numerical simulation with the coolant of R141b flow in a vertical mini-channel. Orientation includes upward and downward. A constant heat flux was loaded at the wall of the channel, of which the width ranges from 1?mm to 3?mm, and a length changes from 200?mm to 400?mm. Subsequently the impact of those parameters that referred to heat flux, inlet temperature and inlet temperature of the coolant could be investigated by carrying out the numerical calculation. In addition, a validation for the model was illustrated in comparison with the previously experimental data.     

Flow boiling heat transfer in two-phase micro-channel heat sinks--II. Annular two-phase flow model

This paper is Part II of a two-part study devoted to measurement and prediction of the saturated flow boiling heat transfer coefficient in water-cooled micro-channel heat sinks. Part I discussed the experimental findings from the study, and identified unique aspects of flow boiling in micro-channels such as abrupt transition to the annular flow regime near the point of zero thermodynamic equilibrium quality, and the decrease in heat transfer coefficient with increasing quality. The operating conditions of water-cooled micro-channels fell outside the recommended range for most prior empirical correlations. In this paper, an annular flow model is developed to predict the saturated flow boiling heat transfer coefficient. Features unique to two-phase micro-channel flow, such as laminar liquid and vapor flow, smooth interface, and strong droplet entrainment and deposition effects, are identified and incorporated into the model. The model correctly captures the unique overall trend of decreasing heat transfer coefficient with increasing vapor quality in the low vapor quality region of micro-channels. Good agreement is achieved between the model predictions and heat transfer coefficient data over broad ranges of flow rate and heat flux.     

物性参数对环形空间流动和传热影响的模拟研究

利用有限容积法,建立了环形空间内单相流体竖直向上流动过程中流动和传热的稳态模型。模型将环形空间内管设置为具有固定生热速率的发热体;流体与内管壁之间设置流动和传热边界层,以更精确的描述壁面位置流体与固体之间动量和热量的耦合传递过程。通过与常物性模型的对比,流体密度、导热系数和黏度随温度变化的变物性模型,在传热能力上具有一定的减少,流体与固体传热面之间的界面剪切力稍有下降。通过比较常物性模型和变物性模型的Re和Ri,结果表明,随着流体强制循环速度的加大,流体物性变化对流动和传热过程的影响逐渐减小。     

仿螺旋肋片内冷通道流动与传热数值分析

利用数值分析方法研究了新型仿螺旋肋片内冷通道的传热与流动特性。采用横截面积为矩形、上下表面带有间断性倾斜矩形肋片叉排且对置的仿螺旋肋片内部冷却通道。分析了在通道宽高比AR=2.9、肋化比Ff/F=2.545、肋高与通道当量直径比e/Dh=0.336、肋间距与肋高比p/e=0.6、肋片与轴面的夹角β=15°及R e在1×104~2×105时的非旋转情况下,R e、肋片与主流方向夹角α等参数对内冷通道强化传热与流动阻力特性的影响。计算结果表明,该仿螺旋矩形肋片作为旋流形成装置起到了迫使流体旋转运动、提高流速和减小层流底层厚度的作用,通道内流体流动达到了预期的螺旋流动效果,通道平均换热系数得到了明显的提高,但同时流动阻力也显著增加。     

Numerical analysis of heat transfer and flow field around cross-flow heat exchanger tube with fouling

Fouling is one of the main problems of heat transfer which can be described as the accumulation on the heat exchanger tubes, i.e.; ash deposits on the heat exchanger unit of the boiler. A decrease in heat transfer rate by this deposition causes loss in system efficiency and leads to increasing in operating and maintenance costs. This problem concerns with the coupling among conduction heat transfer mode between solid of different types, conjugate heat transfer at the interface of solid and fluid, and the conduction/convection heat transfer mode in the fluid which can not be solved analytically. In this paper, fouling effect on heat transfer around a cylinder in cross flow has been studied numerically by using conjugate heat transfer approach. Unlike other numerical techniques in existing literatures, an unstructured control volume finite element method (CVFEM) has been developed in this present work. The study deals with laminar flow where the Reynolds number is limited in the range that the flow field over the cylinder is laminar and steady. We concern the fouling shape as an eccentric annulus with constant thermal properties. The local heat transfer coefficient, temperature distribution and mean heat transfer coefficient along the fouling surface are given for concentric and eccentric cases. From the results, we have found that the heat transfer rate of cross-flow heat exchanger depends on the eccentricity and thermal conductivity ratio between the fouling material and fluid. The effect of eccentric is dominant in the region near the front stagnation point due to high temperature and velocity gradients. The mean Nusselt number varies in asymptotic fashion with the thermal conductivity ratio. Fluid Prandtl number has a prominent effect on the distribution of local Nusselt number and the temperature along the fouling surface.     

Numerical study of the influence of the atmospheric pressure on the heat and mass transfer rates of desiccant wheels

A study was conducted to investigate the influence of the atmospheric pressure, from 101,325 to 60,000 Pa (corresponding from 0 to 4217 m of altitude), on the heat and mass transfer rates of desiccant wheels (hygroscopic matrix of silica gel RD), assuming fixed-mass airflow rates and fixed inlet states for the process and regeneration airflows, defined by the respective inlet temperatures and water vapour contents. The results show a generic reduction of the heat and mass transfer rates as the atmospheric pressure decreases, an effect that is more pronounced for short channel wheels and high airflow rates. Correlations based on the numerical results are presented for the correction of the heat and mass transfer rates when the desiccant wheels are operating at non-standard atmospheric pressure.     

气液两相流动及沸腾传热流体模化分析

基于制冷剂类物质热物性表与水和水蒸气热力性质图表,计算并比较了R12、R22、R123、R134a在相同气液密度比条件下模化水热力过程的能力;应用相似原理与量纲理论分析了流体模化方法在研究气液两相流动特性、高温高压传热过程、临界热流密度现象以及近临界和超临界传热方面的优势与方法.结果表明:制冷剂类物质具有良好的热力性质,其中R134a在相同气液密度比条件下的压力值仅约为水的1/5,而且其臭氧消耗潜能值为0,是环境友好的理想模化流体.提出并探讨了从本质上革新准则数以及基于模型探索新的模化条件和发展新的准则关系式的流体模化技术发展思路.     

Numerical investigation on the heat transfer and flow in the mini-fin heat sink for CPU

Fluid flow and heat transfer in the mini-rectangular fin heat sink for CPU of PC using de-ionized water as working fluid are numerically investigated. Based on the real PC operating conditions, the three-dimensional governing equations for fluid flow and heat transfer characteristics are solved using finite volume scheme. The standard k–ε turbulent model is employed to describe the flow structure and behavior. The predicted results obtained from the model are verified by the measured data. There is a reasonable agreement between the predicted results and experiments. The results of this study are expected to lead to guidelines that will allow the design of the cooling system with improved cooling performance of the electronic equipments increasing reliable operation of these devices.     

错列翅片紧凑式换热器湍流流动及换热性能的数值研究

采用高雷诺数κ-ε湍流模型，对中高雷诺数下紧凑式错列翅片换热器的表面换热及流动特性进行数值模拟。结果表明，该种型式的换热器具有良好的流动和换热性能，拓宽了其空调领域的应用。     

燃气轮机回流式燃烧室内两相流动数值模拟及分析

利用流体分析软件STAR—CD对某回流式燃气轮机燃烧室的内部流场进行了三维冷态数值模拟及两相流反应数值模拟。建立了燃烧室的三维计算几何模型及计算网格,计算了燃烧室的单相流场及喷雾两相流场。在计算中气相采用N-S方程求解,采用高雷诺数κ-ε湍流模型及SIMPISO算法;液相采用Lagrange法处理,采用颗粒群轨道模型。根据计算结果进行流动分析,为进一步进行燃烧室内部燃烧过程的数值计算分析及改善燃烧室的结构设计、降低排放奠定了基础。     

Numerical modelling of heat and mass transfer in finned dehumidifier

In several heat exchange devices, phase transition occurs in a small region adjacent to the wall, and the secondary phase is present only in a thin layer running along the wall, allowing for decoupling between the fluid dynamic computation of the core flow and the numerical analysis of the secondary phase. This happens in finned dehumidifier, but also in spray cooling or defogging problems. In a finned dehumidifier, or in air conditioning evaporators, the secondary phase is provided by moist air condensation, and may consist of discrete droplets, continuous film or a collection of rivulets. Several levels of approximation may be adopted, depending on the specific problem: perfect drain assumption requires only the addition of a heat source in the energy equation, otherwise the water layer behaviour has to be taken into account. Furthermore, a heat and mass transfer analogy may or may not be appropriate; in the latter case, the solution of the diffusion equation of humidity is required.Here, different levels of approximation are compared with literature experimental data for condensation over a vertical fin. Results show that thermal resistance and gravity effects, in the considered geometry, are negligible, and the condensate takes the form of a collection of still droplets, rather than a flowing film. This has an effect on the actual heat transfer and water layer build-up, and the variation of temperature along the fin induces some discrepancy with respect to the straightforward application of the heat and mass transfer analogy.     

Numerical analysis on the fluid flow and heat transfer in the channel with V-shaped wavy lower plate

In the present study, the numerical analysis on the heat transfer and flow developments in the channel with one-side corrugated plate under constant heat flux conditions is presented. The corrugated plate with the corrugated tile angles of 40° is simulated with the channel height of 7.5 mm. The flow and heat transfer developments are simulated by the k-ε standard turbulent model. A finite volume method with the structured uniform grid system is employed for solving the model. Effects of relevant parameters on the heat transfer and flow developments are considered. Breaking and destabilizing in the thermal boundary layer are promoted as fluid flowing through the corrugated surface. Therefore, the corrugated surface has a significant effect on the enhancement of heat transfer.     

Numerical analysis of thermal response tests with a groundwater flow and heat transfer model

The Kelvin line-source equation, used to analyze thermal response tests, describes conductive heat transfer in a homogeneous medium with a constant temperature at infinite boundaries. The equation is based on assumptions that are valid for most ground-coupled heat pump environments with the exception of geological settings where there is significant groundwater flow, heterogeneous distribution of subsurface properties, a high geothermal gradient or significant atmospheric temperature variations. To address these specific cases, an alternative method to analyze thermal response tests was developed. The method consists in estimating parameters by reproducing the output temperature signal recorded during a test with a numerical groundwater flow and heat transfer model. The input temperature signal is specified at the entrance of the ground heat exchanger, where flow and heat transfer are computed in 2D planes representing piping and whose contributions are added to the 3D porous medium. Results obtained with this method are compared to those of the line-source model for a test performed under standard conditions. A second test conducted in waste rock at the South Dump of the Doyon Mine, where conditions deviate from the line-source assumptions, is analyzed with the numerical model. The numerical model improves the representation of the physical processes involved during a thermal response test compared to the line-source equation, without a significant increase in computational time.     

Transient two-phase flow and heat transfer with localized heating in porous media

The transient behavior of two-phase flow and heat transfer in a channel filled with porous media was numerically studied in this paper. Based on the two-phase mixture model, numerical solutions were obtained using the Finite-Volume Method (FVM). Two methods to treat the discontinuous diffusion coefficient in the energy equation, i.e. the harmonic mean method and the “modified” Kirchhoff method were compared. It was found that the “modified” Kirchhoff method was better in dealing with the rapid change in the diffusion coefficient. Three different cases, with discrete heat flux applied at (1) the upper wall, (2) lower wall and (3) both the upper and lower walls were studied. The velocity and temperature fields for these cases were discussed. The results show that the liquid and vapor flow fields, as well as the temperature and liquid saturation fields have distinctly different features with the change in heating location. An analysis of the vapor volume fraction indicates that the largest amount of vapor with the highest vapor generation rate was for the case in which the heat flux is applied from the lower wall.