Rarefied hypersonic flow in a plane channel with a surface glow discharge

The gasdynamic structure of a hypersonic molecular nitrogen flow in a plane channel whose opposite surfaces are segmented electrodes for generating a continuous surface glow discharge is investigated using a two-dimensional computational model. The electrodynamic structure of the surface glow discharge in the hypersonic rarefied gas flow (distributions of the charged particle concentrations, current density, and electric potential) is studied. A two-dimensional conjugate electrical-gasdynamic model consisting of the continuity, Navier-Stokes, and energy conservation equations and the chargedparticle continuity equations in the ambipolar approximation is proposed. The real thermophysical and transport properties of molecular nitrogen are taken into account. It is shown that using a surface glow discharge in a hypersonic rarefied gas flow makes it possible effectively to modify the shock-wave flow structure and hence to consider this type of discharge as additional tool for controlling rarefied gas flows.     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

Parametric investigation of particle acceleration in high enthalpy conical nozzle flows for coating applications

A modified cold gas-dynamic spray technique is under development by using shock tunnel technology, which can enhance the coating quality by increasing the solid particle velocity up to 1,500 m/s. The particle diameter typically amounts to 10 μm. A theoretical model based on gas-particle flows is employed to describe the behaviour of the flow and the solid particles. This quasi-1D model is capable to consider non-equilibrium effects of the gas phase due to high reservoir temperatures, and the influence of wall friction and heat transfer averaged over the nozzle cross section. This model is used for the design and optimization of the nozzle geometry by a parametric study, which results in a conical nozzle with a half opening angle of 2.8° and a length of 325 mm. Particles for coating are injected at about 55 mm downstream of the throat. A shock tunnel facility has been set up at the Shock Wave Laboratory for performing an experimental study of this new technique. The theoretical performance of this setup is evaluated by the KASIMIR simulation software and the quasi-1D method described in this paper. The high reservoir conditions required to achieve particle velocities of 1,500 m/s can be realized by using either a very high driver pressure of about 600 bar for air as driver gas or a relatively low driver pressure of about 200 bar for helium as driver gas.      

Characteristics of turbulence transport for momentum and heat in particle-laden turbulent vertical channel flows

The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, with a focus on the contributions of particle feedback effect to momentum and heat transfer of turbulence. We take into account the effects of particles on flow drag and Nusselt number and explore the possibility of drag reduction in con-junction with heat transfer enhancement in particle-laden turbulent flows.The effects of particles on momentum and heat transfer are analyzed,and the possibility of drag reduc-tion in conjunction with heat transfer enhancement for the prototypical case of particle-laden turbulent channel flows is addressed.We present results of turbulence modification and heat transfer in turbulent particle-laden channel flow,which shows the heat transfer reduction when large inertial parti-cles with low specific heat capacity are added to the flow. However,we also found an enhancement of the heat transfer and a small reduction of the flow drag when particles with high specific heat capacity are involved.The present results show that particles,which are active agents,interact not only with the velocity field,but also the temperature field and can cause a dissimilarity in momentum and heat transport.This demonstrates that the possibility to increase heat transfer and suppress friction drag can be achieved with addition of par-ticles with different thermal properties.     

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.     

Heat Transfer and Gas Flow through Feed Stream within Horizontal Pipe

In the feeding process, the feed stream forms a moving packed bed of particle from the feedstock in the feed channel. When the feeding is at emergency interruption especially in the case of flooding and uncontrollable discharge, the hot gases from reactor would infiltrate into the feed stream. The high heat penetration into feed stream would affect the feeder performance. In this paper, transient thermal response of feed stream within horizontal pipe is described mathematically with a gas flow and heat transfer model. Influences of varied factors on the thermal penetration into feed stream are examined for different conditions. The temperature of the packed-bed particles and the gas velocity distribution curves are obtained for the feeding service at interruption and at normal operating conditions. The numerical results show that the thermal penetration to the packed-bed particles by the seepage flow fluid is high only in the position near the gas entrance. The thermal penetration depth tends to increase with the seepage flow velocity and decrease with feeding rate. There is no appreciable thermal penetration in the feed stream when the feeding service is at normal running. The operating conditions and the porosity of solid bed have importance effects on the gas velocity and temperature field in the thermal penetration zone. A test system is set up to determine the transient thermal response experimentally for the packed bed of particles within a horizontal pipe. The model results are found to compare favorably with the experimental data.     

Asymptotic solutions of the two-dimensional equations for a thin viscous shock layer in a rarefied gas

Two-dimensional hypersonic rarefied gas flow around blunt bodies is investigated for the continuum to free-molecular transition regime. In [1], as a result of an asymptotic analysis, three rarefied gas flow regimes, depending on the relationship between the problem parameters, were detected and one of these regimes was investigated. In the present study, asymptotic solutions of the thin viscous shock layer equations at small Reynolds numbers are obtained for the other two flow regimes. Analytical expressions for the heat transfer, friction and pressure coefficients are obtained as functions of the incident flow parameters and the body geometry and temperature. As the Reynolds number tends to zero, the values of these coefficients approach their values in free-molecular flow. The scaling parameters of hypersonic rarefied gas flow around bodies are determined for different regimes. The asymptotic solutions are compared with the results of direct Monte Carlo simulation.     

Steady Natural Convection Flow of a Particulate Suspension Through a Parallel-Plate Channel

A continuum model for two-phase (fluid/particle) flow induced by natural convection is developed and applied to the problem of steady natural convention flow of a particulate suspension through an infinitely long channel. The walls of the channel are maintained at constant but different temperatures. The two-phase model accounts for particle-phase viscous effects. Boundary conditions borrowed from rarefied gas dynamics are employed for the particle-phase wall conditions. Various closed-form solutions for different special cases are obtained. A parametric study of the physical parameters involved in the problem are performed to illustrate the influence of these parameters on the flow and heat transfer aspects of the problem.     

A novel multiphase DNS approach for handling solid particles in a rarefied gas

A comprehensive model is proposed for multiphase DNS simulations of gas–solid systems involving particles of size comparable to the mean free path of the gas and to that of the bounding geometry. The model can be implemented into any multiphase Direct Numerical Simulation (DNS) method. In the current work, the Volume of Fluid (VOF) method is used, and it is extended to allow for the incorporation of rarefaction effects. For unbounded flow, the model is in excellent agreement with experimental data from the literature. For flows in closed conduits, the model outperforms the alternate approach of using a slip boundary condition at the particle surface for the most relevant degrees of rarefaction and confinement. The proposed model is also able to correctly handle particle–particle interception. The model is intended for low particle Reynolds number flows, and can be applied to resolve in great detail phenomena in a large number of industrial applications (such as filtration of fine particles in porous media).     

Shock-induced Near-wall Two-phase Flow Structure Over a Micron-sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micron-sized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed-phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail. The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.     

Experimental validation of a two equation RANS transitional turbulence model for compressible microflows

Laminar-to-turbulent flow transition in microchannels can be useful to enhance mixing and heat transfer in microsystems. Typically, the small characteristic dimensions of these devices hinder in attaining higher Reynolds numbers to limit the total pressure drop. This is true especially in the presence of a liquid as a working medium. On the contrary, due to lower density, Reynolds number larger than 2000 can be easily reached for gas microflows with an acceptable pressure drop. Since microchannels are used as elementary building blocks of micro heat exchangers and micro heat-sinks, it is essential to predict under which conditions, the laminar-to-turbulent flow transition inside such geometries can be expected. In this paper, experimental validation of a two equations transitional turbulence model, capable of predicting the laminar-to-turbulent flow transition for internal flows as proposed by Abraham etal. (2008), is presented for the first time for microchannels. This is done by employing microchannels in which Nitrogen gas is used as a working fluid. Two different cross-sections namely circular and rectangular are utilized for numerical and experimental investigations. The inlet mass flow rate of the gas is varied to cover all the flow regimes from laminar to fully turbulent flow. Pressure loss experiments are performed for both cross-sectional geometries and friction factor results from experiments and numerical simulations are compared. From the analysis of the friction factor as a function of the Reynolds number, the critical value of the Reynolds number linked to the laminar-to-turbulent transition has been determined. The experimental and numerical critical Reynolds number for all the tested microchannels showed a maximum deviation of less than 12%. These results demonstrate that the transitional turbulence model proposed by Abraham etal. (2008) for internal flows can be extended to microchannels and proficiently employed for the design of micro heat exchangers in presence of gas flows.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

Prediction of turbulent gas–solids flow in curved ducts using the Eulerian–Lagrangian method

The flow of particulate two‐phase flow mixtures occur in several components of solid fuel combustion systems, such as the pressurised fluidised bed combustors (PFBC) and suspension‐fired coal boilers. A detailed understanding of the mixture characteristics in the conveying component can aid in refining and optimising its design. In this study, the flow of an isothermal, dilute two‐phase particulate mixture has been examined in a high curvature duct, which can be representative of that transporting the gas–solid mixture from the hot clean‐up section to the gas turbine combustor in a PFBC plant. The numerical study has been approached by utilising the Eulerian–Lagrangian methodology for describing the characteristics of the fluid and particulate phases. By assuming that the mixture is dilute and the particles are spherical, the governing particle momentum equations have been solved with appropriately prescribed boundary conditions. Turbulence effects on the particle dispersion were represented by a statistical model that accounts for both the turbulent eddy lifetime and the particle transit time scales. For the turbulent flow condition examined it was observed that mixtures with small particle diameters had low interphase slip velocities and low impaction probability with the pipe walls. Increasing the particle diameters (>50 μm) resulted in higher interphase slip velocities and, as expected, their impaction probability with the pipe walls was significantly increased. The particle dispersion is significant for the smaller sizes, whereas the larger particles are relatively insensitive to the gas turbulence. The main particle impaction region, and locations most prone to erosion damage, is estimated to be within an outer duct length of two to six times the duct diameter, when the duct radius of curvature to the duct diameter ratio is equal to unity. Copyright © 1999 John Wiley & Sons, Ltd.     

Effect of inertial particles with different specific heat capacities on heat transfer in particle-laden turbulent flow

The effect of inertial particles with different specific heat on heat transfer in particle-laden turbulent channel flows is studied using the direct numerical simulation(DNS) and the Lagrangian particle tracking method. The simulation uses a two-way coupling model to consider the momentum and thermal interactions between the particles and turbulence. The study shows that the temperature fields display differences between the particle-laden flow with different specific heat particles and the particle-free flow,indicating that the particle specific heat is an important factor that affects the heat transfer process in a particle-laden flow. It is found that the heat transfer capacity of the particle-laden flow gradually increases with the increase of the particle specific heat. This is due to the positive contribution of the particle increase to the heat transfer. In addition,the Nusselt number of a particle-laden flow is compared with that of a particle-free flow.It is found that particles with a large specific heat strengthen heat transfer of turbulent flow, while those with small specific heat weaken heat transfer of turbulent flow.     

Hydrodynamic and thermal fields analysis in gas-solid two-phase flow

The present work aims to investigate numerically the flowfield and heat transfer process in gas-solid suspension in a vertical pneumatic conveying pipe. The Eulerian-Lagrangian model is used to simulate the flow of the two-phases. The gas phase is simulated based on Reynolds Average Navier-Stokes equations (RANS) with low Reynolds number k-ε model, while particle tracking procedure is used for the solid phase. An anisotropic model is used to calculate the Reynolds stresses and the turbulent Prandtl number is calculated as a function of the turbulent viscosity. The model takes into account the lift and drag forces and the effect of particle rotation as well as the particles dispersion by turbulence effect. The effects of inter-particles collisions and turbulence modulation by the solid particles, i.e. four-way coupling, are also included in the model. Comparisons between different models for turbulence modulation with experimental data are carried out to select the best model. The model is validated against published experimental data for velocities of the two phases, turbulence intensity, solids concentration, pressure drop, heat transfer rates and Nusselt number distribution. The comparisons indicate that the present model is able to predict the complex interaction between the two phases in non-isothermal gas-solid flow in the tested range. The results indicate that the particle-particle collision, turbulence dispersion and lift force play a key role in the concentration distribution. In addition, the heat transfer rate increases as the mass loading ratio increases and Nusselt number increases as the pipe diameter increases.     

Modeling gas dynamic particle interaction in rarefied gas flows

In the case of slow, so-called creeping viscous flow a considerable amount of information on the interaction between individual particles or groups of particles has been obtained theoretically by means of the Stokes equations [1]. Regimes in which it is necessary to take inertia, compressibility and the rarefaction of the medium into account (see, for example, [2]) have received much less attention, especially from the experimental standpoint. As a rule, previous experiments have involved freely falling particles in a viscous fluid, but the experimental possibilities of determining the forces and moments exerted on a particle by a liquid or gas at small Reynolds numbers can be considerably expanded by investigating subsonic, flows of rarefied gas obtained with the aid of porous media. The results of such studies of the flow past a spherical particle and its interaction with another particle are presented in this paper.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 152–157, January–February, 1987.     

On the deposition of particles in liquid metals onto vertical ceramic walls

The deposition of non-metallic particles in liquid-metal flows is a serious industrial problem because the build-up of particles on ceramic walls clogs the flow path and interrupts the production, and this leads to large economic losses. This paper is an effort to extend the current state-of-the-art knowledge of particle deposition in air in order to predict particle deposition rates in liquid-metal flows using an improved Eulerian deposition model and considering Brownian and turbulent diffusion, turbophoresis and thermophoresis as transportation mechanisms. The model was used to predict the rate of deposition of particles in an air flow, and the predictions were compared to published measurements to demonstrate its performance. The model was then modified to take into account the differences in properties between air and liquid metals and thereafter applied to liquid-metal flows. Effects on the deposition rate of parameters such as steel flow rate, particle diameter, particle density, wall roughness and temperature gradient near the wall were investigated. It is shown that the steel flow rate has a very important influence on the rate of deposition of large particles, for which turbophoresis is the main deposition mechanism. For small particles, both wall roughness and thermophoresis have a significant influence on the particle deposition rate. Particle deposition rates under various conditions were successfully predicted.     

Heat transfer through rarefied gases between coaxial cylindrical surfaces with arbitrary temperature difference

The problem of nonlinear heat transfer through a rarefied gas confined between concentric cylinders maintained at different temperatures is investigated. The formulation is based on the nonlinear Shakhov kinetic model subject to Cercignani–Lampis boundary conditions, while molecular interaction is modelled by the inverse power law. The detailed behaviour of the radial heat flow, density, temperature and pressure distributions in terms of the normalized temperature difference between the cylindrical walls, the ratio of the two cylindrical radii and the gas rarefaction is investigated and certain interesting characteristics are revealed. The study includes small, moderate and large temperature differences and various radius ratios and is extended in the whole range of the Knudsen number. It is verified that the type of molecular interaction plays an important role when the heat transfer configuration becomes strongly nonlinear, while the influence of the gas–surface scattering law has similar effects both in linear and nonlinear conditions. By comparing linear and nonlinear results corresponding to the same conditions, it is concluded that linearized analysis can capture the correct behaviour of the heat flow configuration not only for infinitesimally small but also for finite temperature differences and that its range of applicability is wider than expected.     

Simulations of dispersed turbulent multiphase flow

Direct numerical simulations of homogeneous isotropic turbulence are used to investigate the effects of turbulence on the transport of particles in gas flows or bubbles in liquid flows. The inertia associated with the bubbles or the particles leads to locally strong concentrations of these in regions of instantaneously strong vorticity for bubbles or strain-rate for particles. This alters the average settling rates and other processes. If the mass-loading of the dispersed phase is significant a random “turbulent” flow is generated by the particle settling. A simple demonstration of this is given, showing the statistically axisymmetric character of this flow and how it can modify an ambient turbulent flow.