Analytical solution of the problem of heat transfer in rarefied gas between two coaxial cylinders

The method of characteristics was used within the framework of the kinetic approach to construct an analytical solution of the problem of heat transfer in a channel whose walls were formed by two coaxial cylinders. The main equation was the Williams kinetic equation, and the boundary condition on the channel walls was the diffusion reflection model. The vector field of the heat flux in the channel was determined, and the specific heat flux through the cross section of the channel was calculated. It was shown that the results obtained for a limiting case, in which the cylinder radii were significantly greater than the mean length of free path of gas molecules, were in good agreement with the results obtained for a plane channel with infinite parallel walls.     

Mass flow rate prediction of pressure–temperature-driven gas flows through micro/nanoscale channels

In this paper, we study mass flow rate of rarefied gas flow through micro/nanoscale channels under simultaneous thermal and pressure gradients using the direct simulation Monte Carlo (DSMC) method. We first compare our DSMC solutions for mass flow rate of pure temperature-driven flow with those of Boltzmann-Krook-Walender equation and Bhatnagar-Gross-Krook solutions. Then, we focus on pressure–temperature-driven flows. The effects of different parameters such as flow rarefaction, channel pressure ratio, wall temperature gradient and flow bulk temperature on the thermal mass flow rate of the pressure–temperature-driven flow are examined. Based on our analysis, we propose a correlated relation that expresses normalized mass flow rate increment due to thermal creep as a function of flow rarefaction, normalized wall temperature gradient and pressure ratio over a wide range of Knudsen number. We examine our predictive relation by simulation of pressure-driven flows under uniform wall heat flux (UWH) boundary condition. Walls under UWH condition have non-uniform temperature distribution, that is, thermal creep effects exist. Our investigation shows that developed analytical relation could predict mass flow rate of rarefied pressure-driven gas flows under UWH condition at early transition regime, that is, up to Knudsen numbers of 0.5.     

A Kinetic Approach to the Plane Poiseuille Flow Over a Porous Matrix

The Poiseuille–Couette gas flow in a channel and the gas flow through an adjacent porous medium are considered when the governing equations are obtained via a molecular kinetic approach based on the Boltzmann equation. The mass continuity, momentum balance and energy conservation are written for the gas in the contiguous regions, whereas the behavior of the solid matrix obeys to the heat diffusion equation. Two different space scalings lead to different forms of the equations for the steady flow through the fully saturated matrix. The boundary conditions at the interface between the two domains are investigated via a matching procedure.     

NUMERICAL STUDY OF TURBULENT FORCED CONVECTION IN A PERIODICALLY RIBBED CHANNEL WITH OSCILLATORY THROUGHFLOW

Numerical computations are performed on the fully developed flow and heat transfer in a periodically ribbed channel with oscillatory throughflow. A uniform heat flux is imposed at the lower plate of the channel. An externally sustained pressure gradient varies sinusoidally in time. A low-turbulent-Reynolds-number version of the k–ϵ two-equation model of turbulence is invoked, together with a preferential dissipation modification, to predict the complex turbulent flow field. Computed results indicate that much heat transfer enhancement is expected by increasing the Womersley number, which measures the relative strength of the oscillatory motion to the viscous effects.     

Heat and mass fluxes upon incomplete accommodation of rarefied gas molecules by the walls of an elliptic channel

The free-molecular rarefied-gas flow in a long channel of elliptic cross-section is considered. The Maxwellian mirror-diffuse reflection model is used as the boundary condition on the channel walls. Under the assumption that a longitudinal temperature gradient is maintained in the channel the mass and heat fluxes are determined as functions of the accommodation coefficient and the ellipse semi-axis ratio. The gas mass velocity profile is obtained. It is shown that for the channels with the same ellipse semi-axis ratios a decrease in the accommodation coefficient leads generally to an increase in the gas mass velocity. The results obtained are compared with those for the case of full gas accommodation on the channel surface.     

Two-phase annular flow and evaporative heat transfer in a microchannel

A physical and mathematical model has been developed to predict the two-phase flow and heat transfer in a microchannel with evaporative heat transfer. Sample solutions to the model were obtained for both constant wall temperature and constant wall heat flux conditions. Results are provided for evaporation rate, liquid film thickness, liquid and vapor phase pressure and temperature distributions. In addition to the sample calculations that were used to illustrate the transport characteristics, computations based on the current model were performed to generate results for comparisons with the experimental results of Qu and Mudawar (2004) where two different mass flow rates of the working fluid were used in the experiment. The comparisons of total pressure drops with the experimental data of Qu and Mudawar (2004) cover the wall heat flux range of 142.71-240 W/cm² with a total channel mass flux of 400.1 kg/m² s and also the wall heat flu range of 99.54-204.39 W/cm² with total channel mass flux of 401.9 kg/m² s. The calculated results from the current model match closely with those of Qu and Mudawar (2004).     

Approximate model for laser ablation of carbon

The paper presents an approximate kinetic theory model of ablation of carbon by a nanosecond laser pulse. The model approximates the process as sublimation and combines conduction heat transfer in the target with the gas dynamics of the ablated plume which are coupled through the boundary conditions at the interface. The ablated mass flux and the temperature of the ablating material are obtained from the assumption that the ablation rate is restricted by the kinetic theory limitation on the maximum mass flux that can be attained in a phase-change process. To account for non-uniform distribution of the laser intensity while keeping the calculation simple the quasi-one-dimensional approximation is used in both gas and solid phases. The results are compared with the predictions of the exact axisymmetric model that uses the conservation relations at the interface derived from the momentum solution of the Boltzmann equation for arbitrary strong evaporation. It is seen that the simpler approximate model provides good accuracy.     

Numerical investigation of rarefied diatomic gas flows through a plane channel into a vacuum

Rarefied flows through a plane microchannel into a vacuum are numerically investigated on the basis of the model kinetic equation for a diatonic gas (nitrogen). The dependence of the gas flow rate through the channel on the Knudsen number, the wall temperature, and the channel length is determined. The energy flux transferred to the cold diatomic gas from the hot channel walls is calculated. The results for diatomic and monatomic gases are compared.     

一种新的玻尔兹曼方程可计算模型构造与分析

临近空间飞行器因各部件尺寸差异较大, 在高空高马赫数条件下飞行会出现多流区共存的多尺度复杂非平衡流动现象, 流场中的气体分子速度分布函数与当地位置、流场分子速度、气体密度、流动速度、温度、热流矢量、应力张量等相关. 通过分析玻尔兹曼方程的一阶查普曼?恩斯科近似解, 构造了一种同时考虑热流矢量和应力张量影响、满足玻尔兹曼方程高阶碰撞矩的跨流域统一可计算模型方程, 并在数学上分析了其守恒性、H定理等基本属性, 证明了新模型方程与玻尔兹曼方程的相容性, 给出了新模型与现有模型如沙克霍夫(Shakhov)模型等的递进关系, 基于碰撞动力学确定了各流域统一气体分子碰撞松弛参数表达式. 在气体动理论统一算法中采用新建模型及现有模型模拟了一维激波结构、二维近空间飞行环境平板和多体圆柱干扰流动, 并与直接模拟蒙特卡洛方法对比分析, 结果表明在流场中粘性效应显著的区域新建模型能更好地捕捉激波位置, 尤其是在激波内部新模型模拟的宏观参数分布与直接模拟蒙特卡洛方法结果符合更好, 验证了新模型的有效性和可靠性, 同时说明在非平衡显著的流动区域碰撞松弛模型受多参数共同作用的影响.      

Numerical simulation of coupled heat and mass transfer in wood dried at high temperature

The mutual effect between heat and mass transfer is investigated for wood dried at high temperature. A numerical model of coupled heat and mass transfer under the effect of the pressure gradient is presented. Based on the macroscopic viewpoint of continuum mechanics, the mathematical model with three independent variables (temperature, moisture content and gas pressure) is constructed. Mass transfer in the pores involves a diffusional flow driven by the gradient of moisture content, convectional flow of gaseous mixture governed by the gradient of gas pressure, the Soret effect and phase change of water. Energy gain or loss due to phase change of water is taken as the heat source. Numerical methods, the finite element method and the finite difference method are used to discretize the spatial and time dimension, respectively. A direct iteration method to solve the nonlinear problem without direct evaluation of the tangential matrix is introduced. The local convergence condition based on the contraction–mapping principle is discussed. The mathematical model is applied to a 3-D wood board dried at high temperature with the Neumann boundary conditions for both temperature and moisture content, and the Dirichlet boundary conditions for gas pressure.     

Fully Developed Forced Convection Heat Transfer in a Porous Channel with Asymmetric Heat Flux Boundary Conditions

An analytical study is performed on steady, laminar, and fully developed forced convection heat transfer in a parallel plate channel with asymmetric uniform heat flux boundary conditions. The channel is filled with a saturated porous medium, and the lower and upper walls are subjected to different uniform heat fluxes. The dimensionless form of the Darcy–Brinkman momentum equation is solved to determine the dimensionless velocity profile, while the dimensionless energy equation is solved to obtain temperature profile for a hydrodynamically and thermally fully developed flow in the channel. Nusselt numbers for the lower and upper walls and an overall Nusselt number are defined. Analytical expressions for determination of the Nusselt numbers and critical heat flux ratio, at which singularities are observed for individual Nusselt numbers, are obtained. Based on the values of critical heat flux ratio and Darcy number, a diagram is provided to determine the direction of heat transfer between the lower or upper walls while the fluid is flowing in the channel.     

Numerical investigation of static flow instability in a low-pressure subcooled boiling channel

A three-dimensional two-fluid model to predict subcooled boiling flow at low pressure is presented. The model is adopted to investigate the two-phase flow and heat transfer characteristics in a heated channel. The presence of bubbles as a consequence of heating flow through a vertical rectangular channel has a significant effect on the overall pressure drop along the channel. Numerical results were compared against a series experimental data performed at various conditions – mass flux, heat flux, inlet temperature and exit pressure. Good agreement on the overall pressure drop was achieved. The onset of flow instability velocity was also accurately determined when compared against measurements. Predicted results of void fraction provided useful information towards a more fundamental understanding of the occurrence of onset of nucleate boiling, onset of significant voiding and onset of flow instability. The phenomenon of boiling onset oscillations was also predicted through the use of the two-fluid model.     

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

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

On the frame dependence of constitutive equations. I. Heat transfer through a rarefied gas between two rotating cylinders

To verify the principle of material frame indifference a numerical calculation of the heat flux field in a rotating gas has been carried out based on the kinetic equation over wide ranges of gas rarefaction and angular velocity. It has been confirmed that a radial gradient of the temperature causes a tangential heat flux. Also, it has been found that the radial heat flux is affected by the rotation.On temporary leave from Department of Physics, Urals State University, 620083 Ekaterinburg, Russia     

Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air

A numerical study of the evaporation in mixed convection of a pure alcohol liquid film: ethanol and methanol was investigated. It is a turbulent liquid film falling on the internal face of a vertical tube. A laminar flow of dry air enters the vertical tube at constant temperature in the downward direction. The wall of the tube is subjected to a constant and uniform heat flux. The model solves the coupled parabolic governing equations in both phases including turbulent liquid film together with the boundary and interfacial conditions. The systems of equations obtained by using an implicit finite difference method are solved by TDMA method. A Van Driest model is adopted to simulate the turbulent liquid film flow. The influence of the inlet liquid flow, Reynolds number in the gas flow and the wall heat flux on the intensity of heat and mass transfers are examined. A comparison between the results obtained for studied alcohols and water in the same conditions is made.     

R410A condensation inside a commercial brazed plate heat exchanger

This paper presents the heat transfer coefficients and the pressure drop measured during HFC-410A condensation inside a commercial brazed plate heat exchanger: the effects of saturation temperature, refrigerant mass flux and vapour super-heating are investigated. The heat transfer coefficients show weak sensitivity to saturation temperature and great sensitivity to refrigerant mass flux and vapour super-heating. At low refrigerant mass flux (<20 kg/m² s) the saturated vapour condensation heat transfer coefficients are not dependent on mass flux and are well predicted by Nusselt [W. Nusselt, Die oberflachenkondensation des wasserdampfes, Energy 60 (1916) 541–546, 569–575] analysis for vertical surface: the condensation process is gravity controlled. For higher refrigerant mass flux (>20 kg/m²s) the saturated vapour condensation heat transfer coefficients depend on mass flux and are well predicted by Akers et al. [W.W. Akers, H.A. Deans, O.K. Crosser, Condensing heat transfer within horizontal tubes, Chem. Eng. Prog. Symp. Series 55 (1959) 171–176] equation: forced convection condensation occurs. In the forced convection condensation region the heat transfer coefficients show a 30% increase for a doubling of the refrigerant mass flux. The condensation heat transfer coefficients of super-heated vapour are 8–10% higher than those of saturated vapour and are well predicted by Webb [R.L. Webb, Convective condensation of superheated vapor, ASME J. Heat Transfer 120 (1998) 418–421] model. A simple linear equation based on the kinetic energy per unit volume of the refrigerant flow is proposed for the computation of the frictional pressure drop.     

Influence of the velocity gradient on the stagnation point heating in hypersonic flow

In a number of experimental and numerical publications a deviation has been found between the measured or computed stagnation point heat flux and that given by the theory of Fay and Riddell. Since the formula of Fay and Riddell is used in many applications to yield a reference heat flux for experiments performed in wind tunnels, for flight testing and numerical simulations, it is important that this reference heat flux is as accurate as possible. There are some shortcomings in experiments and numerical simulations which are responsible in some part for the deviations observed. But, as will be shown in the present paper, there is also a shortcoming on the theoretical side which plays a major role in the deviation between the theoretical and experimental/numerical stagnation point heat fluxes. This is caused by the method used so far to determine the tangential velocity gradient at the stagnation point. This value is important for the stagnation point heat flux, which so far has been determined by a simple Newtonian flow model. In the present paper a new expression for the tangential velocity gradient is derived, which is based on a more realistic flow model. An integral method is used to solve the conservation equations and, for the stagnation point, yields an explicit solution of the tangential velocity gradient. The solution achieved is also valid for high temperature flows with real gas effects. A comparison of numerical and experimental results shows good agreement with the stagnation point heat flux according to the theory of Fay and Riddell, if the tangential velocity gradient is determined by the new theory presented in this paper.This article was processed by the author using theLAT_EX style filepljour2 from Springer-Verlag.     

Numerical simulation of compressible gas flow and heat transfer in a microchannel surrounded by solid media

In this study, a numerical model is developed to investigate the coupled compressible gas flow and heat transfer in a microchannel surrounded by solid media. To accommodate the varying flow cross-section, the compressible gas flow model is established in a non-orthogonal curvilinear coordinate system. An iterative numerical procedure is employed to solve the coupled heat transfer and gas flow equations. The computer code for the compressible gas flow is first validated against two test problems, and then extended by including the heat conduction in the solid media. The effect of the inlet Mach number on the Nusselt number is examined. It is found that the pressure difference from the pyrolysis front to the heated surface is induced essentially by the gas addition from the channel wall, instead from the pyrolysis front. The necessity of accounting for the gas compressibility is clearly demonstrated when severe heating is applied. The pressure distribution obtained along the channel axial direction is useful for further structural analysis of composite materials.     

Certain hydrodynamic characteristics of the evolution of a solitary spherical bubble in the finite compression regime

The characteristics of unsteady incompressible fluid flow around a gas or vapor bubble compressed to a finite size are theoretically investigated. The velocity on the bubble boundary, the integral flow-rate and kinetic energy, and the form of the pressure distribution in the fluid are analyzed. Certain general qualitative regularities independent of viscosity, surface tension, and the intensity of the heat and mass transfer processes are revealed. The general results obtained are illustrated and specified with reference to a simplified polytropic model.     

Mass transfer effects on the non-Newtonian fluids past a vertical plate embedded in a porous medium with non-uniform surface heat flux

 The present study is devoted to investigate the influences of mass transfer on buoyancy induced flow over vertical flat plate embedded in a non-Newtonian fluid saturated porous medium. The Ostwald–de Waele power-law model is used to characterize the non-Newtonian fluid behavior. Similarity solution for the transformed governing equations is obtained with prescribed variable surface heat flux. Numerical results for the details of the velocity, temperature and concentration profiles are shown on graphs. Excess surface temperature as well as concentration gradient at the wall associated with heat flux distributions, which are entered in tables, have been presented for different values of the power-law index n, buoyancy ration B and the exponent λ as well as Lewis number Le. Received on 26 April 2000