多孔介质在压力梯度作用下的热质耦合数值模拟

多孔介质干燥导致热质耦合传输过程。本文基于连续介质力学的宏观尺度,对多孔介质的热、湿和气三者耦合迁移进行数值模拟,研究压力梯度对热质传输的影响。多孔介质传质机理主要为水汽和空气的对流和扩散传输、吸附水在含湿量梯度作用下的自由扩散和其在温度梯度即Soret效应驱动下的流动。采用Galerkin加权余量的有限元方法,提出了...     

Numerical simulation and experimental study of heat and mass transfer in fuel droplet evaporation

Heat and mass transfer in fuel droplet evaporation are investigated through numerical simulation and experimental study. The effect of liquid-phase heat transfer is studied using the temperature difference between the gas- and liquid-phase droplets, different turbulent intensity and oscillatory flow frequency. For the two-droplet array, some differences in heat and mass transfer mechanisms are found. For different spacing of the two-droplet array, the downstream droplet evaporation is affected by the lead droplet.     

Conjugate model for heat and mass transfer of porous wall in the high temperature gas flow

Heat and mass transfer of a porous permeable wall in a high temperature gas dynamical flow is considered. Numerical simulation is conducted on the ground of the conjugate mathematical model which includes filtration and heat transfer equations in a porous body and boundary layer equations on its surface. Such an approach enables one to take into account complex interaction between heat and mass transfer in the gasdynamical flow and in the structure subjected to this flow. The main attention is given to the impact of the intraporous heat transfer intensity on the transpiration cooling efficiency. The project supported by the National Natural Science Foundation of China (19889209) and Russian Foundation for Basic Research (97-02-16943)     

Numerical simulation of turbulent heat transfer from perforated plate-fin heat sinks

Excessive heat from microelectronic components is essential to remove to increase the reliability of the system. In this paper, various types of perforations in the form of small channels such as square, circular, triangular and hexagonal cross sections are introduced and thermal performances are compared to improve the cooling performance of heat sink. The governing equations are solved by adopting a control volume based finite element method with an unstructured non-uniform grid system. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 10⁴ to 4 × 10⁴ based on the fin length and Prandtl number is taken as Pr = 0.71. RANS based k-ε turbulence model is used to predict the turbulent flow parameters. The predicted results are validated by the previously published experimental data and in reasonable agreement with the experiment. Results show that fins having circular perforations have better thermal and fluid dynamic performances than the other types of fins considered here.     

Numerical simulation of effects of flow maldistribution on heat and mass transfer in a PEM fuel cell stack

A 3D numerical study was carried out to analyze flow, heat and mass transfer first in a single half-cell cathode channel of proton exchange membrane (PEM) fuel cell. From practical point of view, it is necessary to put the appropriate number of cells in a stack. Hence, the above study on a single half-cell is extended to a stack of channels. Due to stacking, the assumption of uniform flow distribution would no longer hold true. Therefore, the channel flow-maldistribution is considered. The water formed at the active surface due to the electrochemical reaction diffuses through the porous layer and eventually enters the gas flow duct. The higher gas velocities in the duct result in faster water vapour removal which leads to a lower value of water vapour into the duct and hence a lower Nusselt number.     

Numerical simulation of heat transfer in a micro channel heat sinks using nanofluids

In this study, a numerical simulation of copper microchannel heatsink (MCHS) using nanofluids as coolants is presented. The nanofluid is a mixture of pure water and nanoscale metallic or nonmetallic particles with various volume fractions. Also, the effects of various volume fractions, volumetric flow rate and various materials of nanoparticles on the performance of MCHS have been developed. A three-dimensional computational fluid dynamics model was developed using the commercial software package FLUENT, to investigate the conjugate fluid flow and heat transfer phenomena in micro channel heatsinks. The results show that the cooling performance of a microchannel heat sink with water based nanofluid containing Al₂O₃ (vol 8%) is enhanced by about 4.5% compared with micro channel heatsink with pure water. Nanofluids reduce both the thermal resistance and the temperature difference between the top (heated) surface of the MCHS and inlet nanofluid compared with that pure water. The cooling performance of a micro channel heat sink with metal nanofluids improves compared with that of a micro channel heat sink with oxide metal nanofluids because the thermal conductivity of metal nanofluid is higher than oxide metal nanofluids. Micro channel heat sinks with nanofluids are expected to be good candidates as the next generation cooling devices for removing ultra high heat flux.     

高温下带金属壳PBX炸药低速撞击敏感性数值模拟

炸药撞击感度和热安全性是评价炸药安全性能的重要指标。为了对高温下炸药撞击敏感性变化规律进行可靠预测,本文中通过数值模拟,研究不同预加热温度下带壳PBX炸药装药在小弹丸低速撞击下的热力学响应,得到炸药点火前至点火阶段局部高温区的位置、形态、温度和应变随时间在炸药中分布的变化。结果显示,炸药发生点火的撞击阈值速度与烤燃温度的关系并非单一随温度升高而降低,而是在加热至348.15 K时达到最高;根据温度和应力应变云图分析可得,随着烤燃温度的提高,炸药强度下降,PBX炸药装药局部高温区快速升温的主导因素由局部剪切变为压缩。热软化对炸药的撞击敏感性起重要作用。     

高温高压气体在低含水率介质中迁移的数值模拟

针对气液两相非等温渗流模型高度非线性的特点,发展了适宜的数值离散方法。根据相态转换准则和控制方程的性质,采用最低饱和度法简化算法。空间离散方面,使用有限体积法;时间离散方面,设计了一套包含合理求解顺序的Picard迭代法,解决了方程组强耦合的问题。利用上述数值方法对高温高压气体的迁移行为进行数值模拟,证明了气体在低含水率介质和等效孔隙度的干燥介质内的运动基本一致,并分析了空腔内的气液相态转变过程。在此基础上,研究了多孔介质孔隙度和渗透率对气体压强演化和示踪气体迁移的影响。研究表明,孔隙度越小（相同渗透率）、渗透率越高（相同孔隙度）,示踪气体的迁移距离越远,并给出了估算不同孔隙度和渗透率下迁移距离的半经验公式。     

Experimental heat and mass transfer of the separated and coupled rotating desiccant wheel and heat wheel

The experimental evaluation of the separated and coupled rotating desiccant wheel and heat wheel is reported. The study aims to investigate the performance of the desiccant wheel and of the heat wheel both when operated separately and jointly. The performance evaluation of the desiccant wheel is based on its moisture removal capacity (MRC), moisture removal regeneration (MRR), and moisture mass balance (MMB). In addition, the study used the total energy balance (TEB), sensible coefficient of performance (COP_Sensible), latent coefficient of performance (COP_Latent) and, total coefficient of performance (COP_Total). The performance of the heat wheel is based on its effectiveness. The COP_Sensible, COP_Latent and, COP_Total are used in the performance evaluation of the coupled desiccant wheel and heat wheel. The general results of the study show that the MRC, MRR and MMB coupled with the TEB, COP_Latent, COP_Sensible and COP_Total predict adequately the performance of the desiccant wheel. In addition, the coupled operation of the desiccant wheel and heat wheel, contributed to the reduction of the external thermal energy requirement for the regeneration of the desiccant wheel. This study can be applied in other researches seeking evaluation of the desiccant wheel, heat wheel, and their combined operation. Moreover, the data presented here are significant for the desiccant wheel benchmarking and for evaluation of the desiccant wheel models.     

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering.CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed.3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique.Investigation was focused on low Reynolds number flow(Re=4.6-56.2),which typically occurs in packed bed reactors in bio-chemical fields.Detailed temperature field information was obtained.Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.     

CFD modeling and experimental validation of heat and mass transfer in wood poles subjected to high temperatures: a conjugate approach

In this article, a coupling method is presented in the case of high thermal treatment of a wood pole and a three-dimensional numerical simulation is proposed. The conservation equations for the wood sample are obtained using diffusion equation with variables diffusion coefficients and the incompressible Reynolds averaged Navier–Stokes equations have been solved for the flow field. The connection between the two problems is achieved by expressing the continuity of the state variables and their respective fluxes through the interface. Turbulence closure is obtained by the use of the standard k–ɛ model with the usual wall function treatment. The model equations are solved numerically by the commercial package ANSYS-CFX10. The wood pole was subjected to high temperature treatment under different operating conditions. The model validation is carried out via a comparison between the predicted values with those obtained experimentally. The comparison of the numerical and experimental results shows good agreement, implying that the proposed numerical algorithm can be used as a useful tool in designing high-temperature wood treatment processes. A parametric study was also carried out to determine the effects of several parameters such as initial moisture content, wood aspect ratio and final gas temperature on temperature and moisture content distributions within the samples during heat treatment.     

燃气射流冲击传热特性的数值模拟

针对射流传热问题,利用基于RNGk-ε湍流模型的数值方法模拟了射流垂直冲击平板的流动过程,并与实验数据比较,验证了模型的可行性。在此基础上,以火箭喷管入口参数为入口条件,建立了超音速燃气射流垂直冲击平板和冲击浸没平板的计算模型,分析了不同冲击条件下努塞尔数分布规律和温度分布规律, 论述了超音速射流传热的特性及影响传热特性的因素。得到了冲击距离为(14~18)D的努塞尔数取值范围,并表明冲击距离和射流温度是影响传热效率的关键因素;冲击距离增加,传热效率降低,冲击平板表面的射流温度越高,传热效率越高。     

Forced convective boiling heat transfer in microtubes at low mass and heat fluxes

Convective boiling of HCFC123 and FC72 in 0.19, 0.3 and 0.51 mm ID tubes is investigated. The experimental setup as well as the data reduction procedure has carefully been designed, so that the relative uncertainty interval of the measured heat transfer coefficient in microtubes is kept within ±10%. Up to 70 K liquid superheat over the saturation temperature is observed at low heat and mass fluxes. The onset of the superheat is found to be dependent on the mass flux and the boiling number of the refrigerant examined. In the saturated boiling regime, the heat transfer characteristics are much different from those in conventional-size tubes. The heat transfer coefficient is monotonically decreased with increasing the vapor quality, and becomes independent of the mass flux. Most empirical formulas are not in accordance with the present experimental data. Since the prediction using the nucleate boiling term of Kandlikar’s empirical correlations coincides with the present results, the convection effect should be minor in microtubes. On the other hand, the pressure loss characteristics are qualitatively in accordance with the conventional correlation formula while quantitatively much lower. These phenomena can be explained by the fact that the annular flow prevails in microtubes.     

Numerical study of heat and mass transfer of ammonia-water in falling film evaporator

To investigate the performance of the heat and mass transfer of ammonia water during the process of falling film evaporation in vertical tube evaporator, a mathematical model of evaporation process was developed and solved based on stream function. Then an experimental study of falling film evaporation was carried out in order to validate the mathematical model. A series of parameters, such as velocity, film thickness and concentration, etc., were obtained from the mathematical model. The calculated results show that the average velocity and the film thickness change sharp at the entrance region when x?x?>?100 mm. The film thickness depends largely on the flow rate of solution. It is observed that the heating power and mass flow of solution significantly affect the concentration difference between the inlet and outlet of evaporation tube. The calculated results reveal that the tube length has a significant impact on the amounts of ammonia vapor evaporated. It is suggested that the roll-worked enhanced tube should be used in order to decrease the concentration gradient in the film thickness direction and enhance the heat and mass transfer rate. Furthermore, the experimental and calculated results indicate that the inlet solution concentration has a great influence on the heat exchange capacity, the amounts of ammonia vapor evaporated and the evaporation pressure.     

Numerical investigation of transient heat and mass transfer in a parallel-flow liquid-desiccant absorber

Liquid desiccant systems have received significant attention as a way to reduce latent loads. Tests of liquid desiccant systems in humid climates have shown significant reductions in energy consumption. As moisture in the air is absorbed at the dehumidifier, the temperature of the liquid desiccant increases due to the addition of heat from the enthalpy of condensation of the water vapor. Thus, the coupled effects of heat and mass transfer are relevant for these applications. A two-dimensional mathematical model of the transient heat and mass transfer for an absorber where a thin film of liquid desiccant flows down its walls and dehumidifies the air in parallel-flow configuration is developed and the dynamics of the modeled system are analyzed.     

Numerical simulation of fully developed turbulent flow and heat transfer in annular-sector ducts

 The mixing length theory is employed to simulate the fully developed turbulent heat transfer in annular-sector ducts with five apex angles (θ₀=18,20,24,30,40^∘) and four radius ratios (R _o/R _i=2,3,4,5). The Reynolds number range is 10⁴10⁵. The numerical results agree well with an available correlation which was obtained in following parameter range: θ₀=18,20,24,30,40^∘, R _o/R _i=4 and Re=10⁴5×10⁴. The present work demonstrates that the application range of the correlation can be much extended. Apart from the mixing length theory, the kɛ model with wall function and the Reynolds stress model are also employed. None of the friction factor results predicted by the three models agrees well with the test data. For the heat transfer prediction the mixing length theory seems the best for the cases studied. Received on 17 July 2000 / Published online: 29 November 2001