Onset of Thermal Convection in a Vertical Porous Cylinder with a Partly Conducting and Partly Penetrative Cylinder Wall

The onset of thermal convection in a vertical porous cylinder in three dimensions is investigated analytically. Top and bottom of the cylinder are set to be perfectly heat conducting and impermeable, and is uniformly heated from below. The convection problem is solved for a cylinder wall that is partly conducting and partly penetrative. The expressions for semi-conduction and semi-penetration are based on a porous medium separated from its surroundings by a thin wall. The eigenvalue problem is split into two Helmholtz equations, and the results are expressed by Bessel functions in the radial direction. Comparisons are made with existing solutions for the limit cases of a closed cylinder wall that is either conducting or insulating. Two different models are compared for the kinematic limit condition of an open boundary.     

Two-phase mass injection from the stagnation zone of a blunt body in a hypersonic flow

A mathematical model of the hypersonic steady gas flow over the stagnation zone of an axisymmetric blunt body with given two-phase injection from the surface is proposed. The two-continuum model of a dusty gas [3] is used for describing the flow in the region of the wall. The problem is solved in the boundary layer and thin viscous shock layer approximations. On the basis of the numerical calculations the distribution of the parameters of the carrier and dispersed phases near the axis of symmetry is obtained. The similarity parameters determining the convective heat transfer are found. The stagnation point heat fluxes with and without particles are compared. The range of parameters on which particles can significantly reduce the heat transfer is determined.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 60–66, July–August, 1992.     

激波诱导含灰气体边界层中的粒子分布

本文研究当激波沿着一个固体表面等速地穿越含灰气体运动时所诱导的层流边界层特性。考虑了作用在气体边界层中球形粒子的 Saffman 升力,建议了一种计算近壁区中弥散相密度剖面的方法,并给出了数值计算结果。本文结果表明:在激波后方存在着一个弯曲的薄层区域,其中的粒子密度可以比其波前原始值增加许多倍。这种粒子聚集效应对于工业中粉尘爆炸等实际问题具有重要意义。     

Development and application of wall-function treatments for turbulent forced and mixed convection flows

CFD calculations of turbulent flow near smooth walls generally employ one of two broad strategies to resolve the very influential, complex, but thin near-wall viscosity-affected sub-layer. One approach uses a fine numerical mesh and a turbulence model incorporating viscous influences; the other employs “wall functions”—formulae that attempt to account for the overall resistance of the sublayer to momentum and heat transport. The latter requires only a fraction of the computational effort of the former and is thus strongly favoured for industrial calculations. However, the wall-function performance is often poor, partly because of inappropriate implementations and partly because the schemes themselves have inherent limitations.The present paper reviews the evolution of wall-function strategies. Attention is then given to two new schemes developed by the authors, one based on an analytical treatment and the other on a numerical resolution of the near-wall sub-layer. Several applications are shown of mixed and forced convection.     

The influence of thermal expansion flow on droplet evaporation

It has been demonstrated recently that it follows from conservation of mass that unsteady temperature fields create flow in an incompressible fluid with a temperature-dependent density even in the absence of gravity. The paper studies the influence of thermal expansion flow on spherically symmetric evaporation of an isolated droplet. A model problem of a droplet evaporating at a constant rate is first considered. In this idealized situation one can use the assumption of a thin thermal boundary layer to solve analytically the unsteady moving-boundary heat conduction problem to find the temperature field inside the droplet both with and without the thermal expansion flow. Next evaporation of a fuel droplet in a diesel engine is studied numerically. The heat diffusion equation is solved in the liquid phase while the standard quasi-steady model is used for the gas phase. The results of the calculation show that for high ambient temperatures the influence of the thermal expansion flow on the droplet lifetime can be considerable.     

Effect of a latitudinal temperature gradient on convective motions arising in a rotating spherical layer

The hydrodynamics of planetary atmospheres and Interiors are frequently directly or indirectly connected with convective motions taking place in rotating liquid spherical layers in the field of a central force. Convective stability in a spherical layer at rest, in a central gravity field, was first discussed in [1, 2]. It was shown that the critical Rayleigh number Rao at which convective instability sets in and the wave number of the critical perturbations depend essentially on the thickness of the layer. As in the plane case, the problem of the convective stability of a spherical layer is found to be degenerate, and the form of the critical perturbations cannot be determined from the linear problem. In actuality, minimization of the Rayleigh number permits establishing only the wave numberl for the spherical harmonic Y _l ^m (θ, ?), realized at the limit of stability; the parameter m remains indeterminate and thus 2l+1 independent convective modes correspond to Rao. In [3] a study was made of the convective stability of a liquid in a slowly rotating thin spherical layer. It was shown that the presence of rotation eliminates the degeneracy; at the limit of stability there arise motions corresponding to the Y _l l (θ, ?) -harmonic with a degenerate maximum at the equator, and propagating in a wave manner toward the side opposite to the rotation. In the present work a study is made of the convective stability of a flow of liquid, arising in a rotating spherical layer due to a nonuniform distribution of the temperatures at one of the boundaries of the layer. In such a statement of the problem it is possible to model large-scale motions in the atmospheres of large planets having internal sources of heat and absorbing solar radiation near the cloud cover of the atmosphere. It is established that, depending on the relationships between the parameters imparting the rotation and the inhomogeneous distribution of the temperature, there is either stabilization or destabilization of the layer in comparison with a fixed layer of the same thickness and with the same, but uniformly distributed heat flux supplied to the layer. A study is made of the form of the corresponding critical perturbations.     

On Disturbance Propagation in Internal-Flow Boundary Layers

The unsteady interaction of plane-channel wall boundary layers with a supersonic inviscid flow is investigated. The flow regimes in which disturbances introduced by the boundary layer developing on one wall influence the boundary layer on the other wall are considered. The regime of relatively large pressure disturbance amplitudes generated near the nozzle outlet or by deforming the channel walls is studied. In these conditions, the interaction process is described by a system of Burgers equations with retarded arguments. Numerical solutions of this system are obtained for symmetric and antisymmetric perturbations of the channel walls.     

The distribution of particles in a shock-induced boundary layer of a dusty gas over a solid surface

The laminar boundary layer behind a constant-speed shock wave moving through a dusty gas along a solid surface is studied. The Saffman lift force acting on a spherical particle in a gas boundary layer is taken into account. A method for calculating the density profile of dispersed phase near the wall is proposed and some numerical results are given. It is shown that behind the shock wave, there exists a curved thin layer where the density of particles is many times higher than the original one. This dust collection effect may be of essential importance to the problem of dust explosion in industry.     

High Péclet number heat exchange between cocurrent streams

Summary  This paper concentrates on the analysis of the heat transfer between two cocurrent laminar flows in parallel channels. For high values of the Péclet number Pe, a boundary layer arises near the wall separating the streams. Matched asymptotic expansions (MAE) are used to obtain approximate solutions. We consider arbitrary inlet temperatures and derive higher-order corrections of the boundary problem. The separating wall is supposed to be sufficiently thin to neglect the heat conduction in it. Analyticity and adiabatic conditions at the outer walls impose restrictions on the inlet temperatures. It turns out, however, that only the inlet temperatures at the wall separating the two fluids enter the leading-order problem. The Nusselt numbers thus calculated are in the leading order proportional to (Pe/x)^1/3, where x is the stream-wise coordinate. An estimate of the thickness of the separating wall to validate the MAE approach is obtained. It is also demonstrated that the MAE analysis is unable to describe the heat exchange of counterflowing fluids. Received 9 June 1999; accepted for publication 17 November 1999     

Wall shape effects on multiphase flow in channels

Two-fluid flow is examined analytically and numerically for increased flow rates through a channel with surface roughness or branching or both. The viscosity and density ratios of the fluids are of order unity. There is much concern in terms of applications as well as fluid dynamical phenomena in configurations where one fluid is present only as a thin layer near an outer wall, leaving the other fluid occupying the channel core and part of a viscous wall layer. The interactive dynamics in both regions is studied and numerical and asymptotic analyses are performed. The major situations examined are: the flow to two symmetrically bifurcating daughters and the flow in a single channel over a rough wall, as well as a combination of the two situations. The principal phenomena of interest are conditions for flow reversal, the presence of upstream influence and the trajectories of the injected fluid as the density or viscosity ratio varies. Special relatively thin or thinning wall layers are produced when the core fluid viscosity increases or when the fluid travels downstream into a daughter vessel.     

Evaporation of Falling and Shear-Driven Thin Films on Smooth and Grooved Surfaces

One of the most important tasks in development of modern gas turbine combustors is the reduction of NO_x emissions. An effective way to reduce the NO_x emission is using the lean premixed prevaporization (LPP) concept. An important phenomenon taking place in LPP chambers is the evaporation of thin fuel films. To increase the fuel evaporation rate, the use of microstructured walls has been suggested. The wall microstructures make use of the capillary forces to evenly distribute the liquid fuel over the wall, so that the appearance of uncontrolled dry patches can be avoided. Moreover, the wall structures promote the thin film evaporation characterized by ultra-high evaporation rates. An experimental setup was built for the investigation of thin liquid films falling down on the outer surface of vertical tubes with either a smooth or structured surface. In the first testing phase water is used, fuel like liquids will be used later on. The thin film can be heated from both sides, by hot oil flowing inside the tube, and by hot compressed air flowing in co-current direction to the thin film. The film is partly evaporated along the flow. Results for the wavy film structure at different Reynolds numbers are reported. For theoretical investigations a model describing the hydrodynamics and heat transfer due to evaporation of the gravity- and shear-driven undisturbed liquid film on structured surfaces was developed. For low Reynolds numbers or low liquid mass fluxes the wall surface is only partly covered with liquid and the heat transfer is shown to be governed by the evaporation of the ultra-thin film in the vicinity of the three-phase contact line. A numerical model for the solution of a two-dimensional free-surface flow of a liquid film over a structured wall was also developed. The Navier–Stokes equations are solved using the Volume of Fluid (VOF) technique. The energy equation is included in the model. The model is verified by comparison with data from the literature showing favorable agreement. In particular, the proposed model predicts the formation of capillary waves observed in the experiments. The model is used to investigate the flow of liquid on a structured wall. This calculation is the first step towards the modeling of a three-dimensional wavy flow of a gravity- and shear-driven film along a wall with longitudinal grooves. It is found that due to the Marangoni effect, a circulating flow arises within the cavity, thereby leading to an enhancement in the evaporation rate.     

Effects of vibrations on particle motion near a wall: Existence of attraction force

The effects of small vibrations on a particle oscillating near a solid wall in a fluid cell, relevant to material processing such as crystal growth in space, have been investigated experimentally and theoretically. Assuming the boundary layer around the particle to be thin compared to the particle radius at high vibration frequencies, an inviscid fluid model was developed to predict the motion of a spherical particle placed near a wall of a rectangular liquid-filled cell subjected to a sinusoidal vibration. Under these conditions, a non-uniform pressure distribution around the particle results in an average pressure that gives rise to an attraction force. Theoretical expressions for the attraction force are derived for the particle vibrating normal to and parallel with the nearest cell wall. The magnitude of this attractive force has been verified experimentally by measuring the motion of a steel particle suspended in the fluid cell by a thin wire. Experiments performed at high frequencies showed that the mean particle position, when the particle is brought near a cell wall, shifts towards the same wall, and is dependent on the cell amplitude and frequency, particle and fluid densities.     

The influence of boundary layers on supersonic inlet flow unstart induced by mass injection

A transverse jet is injected into a supersonic model inlet flow to induce unstart. Planar laser Rayleigh scattering from condensed CO₂ particles is used to visualize flow dynamics during the unstart process, while in some cases, wall pressure traces are simultaneously recorded. Studies conducted over a range of inlet configurations reveal that the presence of turbulent wall boundary layers strongly affect the unstart dynamics. It is found that relatively thick turbulent boundary layers in asymmetric wall boundary layer conditions prompt the formation of unstart shocks; in symmetric boundary conditions lead to the propagation of pseudo-shocks; and in both cases facilitate fast inlet unstart, when compared with thin, laminar boundary layers. Incident shockwaves and associated reflections are found to affect the speed of pressure disturbances. These disturbances, which induce boundary layer separation, are found to precede the formation of unstart shocks. The results confirm the importance of and need to better understand shock-boundary layer interactions in inlet unstart dynamics.     

Numerical studies on erosive burning in cylindrical solid propellant grain

This paper addresses erosive burning of a cylindrical composite propellant grain. Equations governing the steady axisymmetric, chemically reacting boundary layer are solved numerically. The turbulence is described by the two equation (k-ɛ) model and Spalding’s eddy break up model is employed for the gas phase reaction rate. The governing equations are transformed and solved in the normalized stream function coordinate system. The results indicate that the dominant reaction zone lies within 20% of the boundary layer thickness close to the wall. The sharp gradient of the temperature profile near the wall is found responsible for bringing the maximum heat release zone near the surface and hence enhancement in the burning rate. The model reproduces the experimental observation that erosive burning commences only above a threshold value of axial velocity.     

On constructing the theory of ductile fracture near friction surfaces

A criterion of ductile fracture is proposed, which takes into account the singular character of theoretical solutions near the maximum friction surfaces and the emergence of a thin layer with intense plastic strains near surfaces with high friction stresses in real processes of metal forming. The equation for the thickness of the layer with intense plastic strains and the fracture criterion include the strain rate intensity factor, apparently, characterizing the intensity of physical processes that occur in a thin material layer near the friction surfaces. Some experimental data are used to determine the thickness of this layer. The ductile fracture criterion is analyzed by solving the problem of strip extrusion under conditions of plane strain deformation.     

哈特曼边界层的初级稳定性分析

导电流体在法向外置磁场的作用下,在贴近壁面处会形成哈特曼边界层.哈特曼边界层的稳定性研究对电磁冶金过程和热核聚变反应冷却系统等相关设备的设计和运行都有着十分重要的意义.本文采用非正则模态稳定性分析方法,对两无限大绝缘平行平板内导电流体流动的稳定性进行了研究.通过在时间上迭代求解扰动变量的控制方程组和伴随控制方程组,获得了在磁场作用下初级扰动的增长情况及其空间分布形式,分析了磁场强度对最优扰动增长倍数G_max、最优展向波数β_opt和最优时刻t_opt的影响,并考察了上下两个哈特曼边界层之间的相互作用.结果表明,最优初始扰动的空间分布形式为沿着流场方向的漩涡,关于法向方向对称或者反对称.当哈特曼数Ha较大时(Ha>10),对称漩涡和反对称漩涡形式的初始扰动增长倍数基本相等;上下两个哈特曼边界层可以认为是彼此独立的,不会相互影响,此时最优扰动增长倍数G_max与局部雷诺数R的平方成正比,相应的最优展向波数β_opt和最优时刻t_opt均正比于哈特曼数Ha.当哈特曼数Ha较小时(Ha<10),反对称漩涡形式的初始扰动更为不稳定,其增长倍数大于对称漩涡的增长倍数,且上下两个边界层之间存在着一定的相互作用,并对整个流场的稳定性产生一定的影响.      

Two-phase flow over a surface with the formation of a liquid film by particle deposition

A consistent asymptotic theory of wall flow with film formation is constructed with reference to subsonic two-phase flow over a blunt body. The external flow problem and the film equations are solved simultaneously. This formulation of the problem supplements the investigation carried out in [4] in which particles deposited on the surface were assumed to disappear from the flow. It is shown that depending on the values of the governing parameters the flow in the film should be described either by the boundary layer equations or by the equations of creeping flow in a layer of unknown thickness. At the outer edge of the film the mass, momentum and energy fluxes found from the numerical solution of the flow problem are given. The case of isothermal film flow on the front of a sphere is investigated. The thickness of the film and the friction and heat transfer coefficients near the axis of symmetry are found for nonisothermal flows. The conditions under which the presence of a film significantly reduces the heat flow to the wall are determined. A similar formulation of the problem (but with another type of mass, momentum and energy sources at the outer edge) is encountered in problems of film condensation on a cold surface [5, 6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 85–92, July–August, 1989.     

高超声速钝锥流场红外辐射实验研究

报道了在JF-10氢氧爆轰驱动高焓激波风洞中开展的再入流场红外辐射实验研究. 风洞的试 验状态为：驻室总压19.6MPa, 总焓15.5MJ/kg, 自由流速度约5 km/s. 实验以锑化铟多元红外成像系统为测量手段，以球头钝锥体为试验模型，测量激波层与近尾流中红外辐射功率 的横向分布剖面. 试验数据呈现明显的规律性. 试验结果表明，激波层内壁面附近的红外辐 射功率较小，中间有一区域辐射较大且相对均匀，激波层外缘辐射单调减小；尾流中红外辐 射功率在轴线附近的核心区最大，随着离轴线距离的增大而单调减小.     

Transient natural convection in a cavity with heat input and a constant temperature wall on opposite sides

The transient natural convection of a fluid with Prandtl number of order 200 in a two-dimensional square cavity has been numerically studied. One of the vertical walls of the cavity is kept at a constant (ambient) temperature and a constant heat flux is applied on the opposite wall. The other walls are adiabatic. Initially, a boundary layer is formed near the heated wall; subsequently, a large vortical structure is generated, together with an upper intrusion layer. As time progresses, the average temperature in the cavity increases, and a descending boundary layer is formed near the constant temperature wall. During the transition to the steady-state regime, a thermal stratification pattern is formed. The results are compared with the scale analysis presented by Patterson and Imberger (1980).