Statistics of Particle Accumulation in Spatially Developing Turbulent Boundary Layers

We present the results of a Direct Numerical Simulation of a particle-laden spatially developing turbulent boundary layer up to Re _θ ?=?2500. Two main features differentiate the behavior of inertial particles in a zero-pressure-gradient turbulent boundary layer from the more commonly studied case of a parallel channel flow. The first is the variation along the streamwise direction of the local dimensionless parameters defining the fluid-particle interactions. The second is the coexistence of an irrotational free-stream and a near-wall rotational turbulent flow. As concerns the first issue, an inner and an outer Stokes number can be defined using inner and outer flow units. The inner Stokes number governs the near-wall behavior similarly to the case of channel flow. To understand the effect of a laminar-turbulent interface, we examine the behavior of particles initially released in the free stream and show that they present a distinct behavior with respect to those directly injected inside the boundary layer. A region of minimum concentration occurs inside the turbulent boundary layer at about one displacement thickness from the wall. Its formation is due to the competition between two transport mechanisms: a relatively slow turbulent diffusion towards the buffer layer and a fast turbophoretic drift towards the wall.     

Thermal Non-equilibrium Natural Convection in a Square Enclosure with Heat-Generating Porous Layer on Inner Walls

Differentially heated enclosure with heat-generating porous layer on inner walls is studied computationally for non-Darcy flow and thermal non-equilibrium models. In this study, this problem is investigated for different internal and external Rayleigh numbers, Darcy numbers, porosity-scaled thermal conductivity ratio, solid-/fluid-scaled heat transfer coefficient and dimensionless thickness of the porous layer. The results indicate that the dimensionless thickness of the porous layer has an important effect on the heat transfer in the enclosure. It was found that the thermal non-equilibrium model is needed for small values of the porosity-scaled thermal conductivity ratio and the solid-/fluid-scaled heat transfer coefficient. It is shown that the convection of heat due to internal heat generation is increased in the enclosure when the ratio of internal Rayleigh number to external Rayleigh number is larger.     

Regime diagram of the development of long-wave near-wall disturbances in a hypersonic boundary layer

A regime diagram of the development of slow near-wall disturbances induced by an unsteady self-induced pressure perturbation in a hypersonic boundary layer is constructed for a disturbance wavelength greater than the boundary layer thickness. It is shown that the main factors shaping the perturbed flow are the gas enthalpy near the body surface, the intensity of the viscous-inviscid interaction, and the nature (sub- or supersonic) of the main part of the boundary layer. Nonlinear boundary-value problems are formulated for regimes in which the near-wall boundary layer region plays a decisive role. Numerical and analytical solutions are obtained in the linear approximation. It is shown that intensification of the viscous-inviscid interaction or an increase in the role of the supersonic main region of the boundary layer impart generally supersonic properties to the main part of the boundary layer, i.e. the upstream propagation of the disturbances is damped and the disturbance growth downstream becomes more intense. Damping of the viscous-inviscid interaction and an increase in the role of the subsonic main part of the boundary layer have the opposite effect. Surface cooling increases the effect of the main part of the boundary layer on the formation of pressure disturbances and surface heating leads to an increase in the effect of the near-wall boundary layer region. It is also shown that for the regimes considered disturbances propagate in a direction opposite to that of the free stream from the turbulent flow region located downstream of the local disturbance development region.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 59–71. Original Russian Text Copyright © 2004 by Bogolepov and Neiland.     

The flow of a non-Newtonian liquid near a rotating disk

Synopsis The flow of non-Newtonian liquid near a rotating disk has been discussed by using second order stress strain velocity relations of classical hydrodynamics. It is found that the effect of cross-viscosity depends on a non-dimensional number R _c. The boundary layer thickness decreases and the dimensionless moment coefficient increases with the increase of R _c.     

Tollmien-Schlichting Waves in a Hypersonic Boundary Layer Flow in the Neighborhood of a Cooled Surface

A nonlinear time-dependent model of the development of longwave perturbations in a hypersonic boundary layer flow in the neighborhood of a cooled surface is constructed. The pressure in the flow is assumed to be induced the combined variation of the thicknesses of the near-wall and main parts of the boundary layer. Numerical and analytic solutions are obtained in the linear approximation. It is shown that if the main part of the boundary layer is subsonic as a whole, its action reduces the perturbation damping upstream and the perturbation growth downstream, while a supersonic, as a whole, main part of the boundary layer creates the opposite effects. An analysis of the solutions obtained makes it possible to conclude that the asymptotic model proposed can describe the three-dimensional instability of the Tollmien-Schlichting waves.     

Filtration model of longitudinal flow in a finned cylindrical channel

The problem of laminar flow of a viscous incompressible fluid in a finned circular tube is considered. A solution is obtained in the form of series in eigenfunctions of the Laplace operator; the coefficients in the series are found numerically. For the same problem, a simpler filtration approximation is proposed in which the system of fins is modeled by a radially inhomogeneous porous layer, and fluid flow in it is described by the Brinkman equation. A formula for the effective permeability of the porous medium is obtained by varying the number and height of fins. The formula provides an accurate evaluation of the mean flow velocity and viscous drag coefficient in finned channels.     

Steady Flows in an Oscillating Spheroidal Cavity with Elastic Wall

The steady flow arising in a spheroidal cavity with periodically-deformed elastic wall is studied experimentally. It is found that average flows whose intensities and structures depend on the wall oscillation frequency and amplitude can develop in the fluid. The average flow is generated in the Stokes boundary layer whose relative thickness is characterized by the dimensionless frequency of the vibrational action. Flow in the form of a pair of toroidal vortices which occupy the entire cavity volume can be observed over the range of low dimensionless frequencies when the boundary layer thickness is comparable with the characteristic cavity dimension. Increase in the dimensionless frequency (decrease in the relative thickness of the Stokes layers) leads to a displacement of the primary vortices towards the cavity boundary. In this case secondary vortices with opposite swirling are formed in the central part of the cavity above the primary vortices. The further increase in the dimensionless frequency leads to development of the secondary vortices and growth of the flow intensity. The large-scale secondary vortices occupy almost the entire cavity volume over the range of high dimensionless frequencies. The dependences of the regimes of average flows and their intensities on the control dimensionless parameters, the oscillation amplitude and frequency, are found on the basis of the results of the investigation.     

Direct numerical simulation of a particle-laden flow in a flat plate boundary layer

In this paper, a direct numerical simulation of particle-laden flow in a flat plate boundary layer is performed, using the Eulerian–Lagrangian point-particle approach. This is, as far as we know, the first simulation of a particle-laden spatially-developing turbulent boundary layer with two-way coupling. A local minimum of the particle number density is observed in the close vicinity of the wall. The present simulation results indicate that the inertial particles displace the quasi-streamwise vortices towards the wall, which, in turn, enhance the mean streamwise fluid velocity. As a result, the skin-friction coefficient is increased whereas the boundary layer integral thicknesses are reduced. The presence of particles augments the streamwise fluctuating velocity in the near-wall region but attenuates it in the outer layer. Nevertheless, the wall-normal and spanwise velocity fluctuations are significantly damped, and so is the Reynolds stress. In addition, the combined effect of a reduced energy production and an increased viscous dissipation leads to the attenuation of the turbulent kinetic energy.     

Correlation for frost properties on a cold cylinder surface in cross flow

An experimental study was reported earlier on the development of frost formation by humid flow passing over the cylinder. In this study, dimensionless correlations based on previous experimental data and reported empirical correlations of others for frost properties are proposed. The frost conduction coefficient is determined by using an analytical equation. Subsequently, correlations are sought for the heat conduction coefficient, dimensionless temperature, dimensionless thickness and density. The advantages of these correlations to any other proposed correlations are their explicit and dependency on time. Furthermore, an estimation of characteristics of the frost is followed by using the same approach and the established correlations.     

Experimental investigation of helical structures in swirling flows

In this investigation the flow in a generic swirl tube with a tangential double-inlet swirl generator and variable exit orifices was experimentally investigated. Using magnetic resonance velocimetry (MRV) three-dimensional, three-component velocity fields were measured for two different Reynolds numbers: 10,000 and 15,000, and for three different exit orifices. The swirl generator had a fixed geometry producing an initial swirl number of 1.6 for all cases. One major observation is the occurrence of a three-layered flow structure. An annular main flow was surrounded by a recirculation zone, as reported in previous literature. However, this recirculation zone – also of an annular shape – exhibited a third layer inside: a thin, high speed jet in the center of the tube with the same flow direction as the main flow. Therefore, the conventional classification of swirling flows into ring and recirculation zone, has to be extended by a core zone. This three-layered flow structure develops independently of the exit configuration. Helical structures were observed in the near-wall region for all cases investigated. Applying an eccentric exit orifice results in the occurrence of strong stationary helical structures not only in the near-wall region but also in the center of the tube. The results, deviating significantly from previous results in the literature, underline the need for more detailed research on the topic of cyclone type flows.     

High-velocity penetration of a melting solid by a slender body

The high-velocity penetration of a melting solid by a thermally insulated slender body is considered. Under certain constraints on the dimensionless melting parameters the flow in the molten layer can be described within the framework of lubrication theory. The local angle of inclination of the body and the surfaces of the molten layer with respect to the velocity is assumed to be small and is taken into account in the linear approximation. The heat flow into the solid is found by simulating the body and the molten layer by means of a segment with distributed heat sources. Within the framework of this simple formulation a closed solution of the problem of the fusion zone around a moving slender body is constructed. The dependence of the shape of the molten layer and the structure of the temperature and longitudinal velocity fields in the layer on the shape of the body and the other governing parameters of the problem is investigated. The results obtained also give a solution of the problem of the melting of a solid rubbing at high velocity against a thermally insulated rough substrate, when the characteristic height of the roughness is of the order of the thickness of the layer and the characteristic length of the order of the contact length.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 43–48, November–December, 1992.     

Effects of velocity fluctuations on heat transfer enhancement

The three-dimensional velocity fluctuation effects on heat transfer enhancement were experimentally investigated using a wind tunnel system and cylinders placed upstream of the test section in the wind tunnel. The cylinders with different diameters were used as turbulators to generate vortical flow motions with three-dimensional velocity fluctuations. A heated plate, part of the tunnel wall, was placed far downstream of the cylinders such that it was subjected mainly to flows with velocity fluctuations but with negligible steady vortical motions. These studies included three-component velocity measurements to characterize the near-wall and cross-section velocity fields and to obtain the turbulent kinetic energy. The temperatures were measured by thermocouples on the heated plate to obtain the convection heat transfer coefficients and the Nusselt numbers. Results indicate that the heat transfer was enhanced by the velocity fluctuations, which is attributed to the modification of boundary layer velocity profiles without the modification of boundary layer thickness. The resulting normalized Nusselt number was approximately a parabolic function of a dimensionless parameter, the product of Reynolds number and normalized turbulent kinetic energy.     

Formation of particle jetting in a cylindrical shock tube

A dense, solid particle flow is numerically studied at a mesoscale level for a cylindrical shock tube problem. The shock tube consists of a central high pressure gas driver section and an annular solid powder bed with air in void regions as a driven section with its far end adjacent to ambient air. Simulations are conducted to explore the fundamental phenomena, causing clustering of particles and formation of coherent particle jet structures in such a dense solid flow. The influence of a range of parameters is investigated, including driver pressure, particle morphology, particle distribution and powder bed configuration. The results indicate that the physical mechanism responsible for this phenomenon is twofold: the driver gas jet flow induced by the shock wave as it passes through the initial gaps between the particles in the innermost layer of the powder bed, and the chaining of solid particles by inelastic collision. The particle jet forming time is determined as the time when the motion of the outermost particle layer of the powder bed is first detected. The maximum number of particle jets is bounded by the total number of particles in the innermost layer of the powder bed. The number of particle jets is mainly a function of the number of particles in the innermost layer and the mass ratio of the powder bed to the gas in the driver section, or the ratio of powder bed mass (in dimensionless form) to the pressure ratio between the driver and driven sections.     

Gas displacement of viscoplastic liquids in capillary tubes

The displacement of viscoplastic liquids in capillary tubes by gas injection is examined. The viscoplasticity alters the flow kinematics and changes dramatically the amount of mass left attached at the tube wall as compared to the Newtonian case, studied experimentally by G.I. Taylor in 1961 [G.I. Taylor, Deposition of a viscous fluid on the wall of a tube, J. Fluid Mech. 10 (1961) 161–165]. Experiments with Carbopol aqueous solutions were performed for different flow rates. A recently proposed viscosity function for viscoplastic liquids was fitted to the rheological data of the Carbopol solutions. A new dimensionless rheological property – the jump number – arises in the dimensionless version of this viscosity function. The results show the effect of the viscoplastic character of the liquid on the free surface shape and on the thickness of the film of liquid left attached to the wall. This thickness decreases with the jump number and increases with the flow rate. It is also observed that there is a critical dimensionless flow rate below which the displacement is apparently perfect, i.e. there is no observable liquid left attached to the wall. This behavior is shown to be directly related to the fully developed flow far ahead the air–liquid interface.     

Possibility of drag reduction usingd-type roughness

An experiment was carried out in a low-speed wind tunnel to study the turbulence structure of the boundary layer over a two-dimensional square cavity on a flat plate. The main purpose of this investigation is to examine the way a square cavity modifies the near-wall structure of the turbulent boundary layer leading to a possible drag reduction overd-type roughness. The experimental results on pressure coefficient and friction coefficient indicated a small reduction in total drag in this configuration. This seems to be due to the stable vortex flow observed within the cavity which absorbs and reorganizes the incoming turbulence in the cavity, thereby modifying the near-wall turbulence structure of the boundary layer. The resultant turbulence structure was very similar to that over drag-reducing riblets surface.     

Possibility of drag reduction usingd-type roughness

An experiment was carried out in a low-speed wind tunnel to study the turbulence structure of the boundary layer over a two-dimensional square cavity on a flat plate. The main purpose of this investigation is to examine the way a square cavity modifies the near-wall structure of the turbulent boundary layer leading to a possible drag reduction overd-type roughness. The experimental results on pressure coefficient and friction coefficient indicated a small reduction in total drag in this configuration. This seems to be due to the stable vortex flow observed within the cavity which absorbs and reorganizes the incoming turbulence in the cavity, thereby modifying the near-wall turbulence structure of the boundary layer. The resultant turbulence structure was very similar to that over drag-reducing riblets surface.     

Free convection from an arbitrarily inclined plate in a porous medium

A mathematical model for the flow and heat transfer in the free convection from an arbitrary inclined isothermal flat plate embedded in a porous medium is presented, in which the Darcy–Boussinesq approximation is adopted to account for bouyancy force. A novel inclination parameter ξ is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed boundary layer equations. Moreover, the similarity equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting ξ=0 and ξ=1, respectively. Detailed results for the skin friction coefficient and Nusselt number as well as for the dimensionless velocity and temperature profiles are presented for a wide range of the parameter ξ. A comparison with similarity solution shows excellent agreement.     

Numerical Study of the Near-Wall Gas-Droplet Jet in a Tube with a Heat Flux on the Surface

A model for calculating the flow of a turbulent mixture of air and suspended liquid particles injected into the near-wall region is developed within a unified approach of mechanics of heterogeneous media in the two-velocity and two-temperature approximation of the Eulerian approach. The influence of droplet evaporation in the near-wall jet on heat transfer between the two-phase gas-droplet flow and the wall is studied in the case of heat addition to the latter. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 5–17, January–February, 2006.     

Simultaneous convection and radiation in flow between two parallel plates

A numerical solution is described for simultaneous forced convection and radiation in flow between two parallel plates forming ahannel. The front plate is transparent to thermal radiation while the back one is thermally insulated. Analyses for both flow and heat are presented for the case of a non-emitting ‘blackened’ fluid. The governing equations of the stream function and the temperature together with their boundary conditions are presented in non-dimensional expressions. The solution is found to depend on eight dimensionless parameters, namely the ratio of the height of the channel to the distance between the plates, the initial dimensionless temperature, the optical thickness, the absorptivities of both plates, the Reynolds number, the Prandtl number and the heat transfer coefficient from the front plate to the surroundings. The numerical solution is obtained using a finite-difference technique. A study has been made of the effect of the initial temperature of the flow at the channel inlet, the dimensionless loss coefficient from the front plate, the absorptivity of the back plate and the optical thickness, on the temperature distribution in the channel, the heat collection efficiency and the average temperature rise in the channel. Results showed that increasing the optical thickness increases the temperature of the front plate and decreases the temperature of the back plate. Also, increasing the optical thickness increases the efficiency of heat collection, which reaches its maximum asymptotic value at an optical thickness of about 1.5. Moreover, the location of the maximum temperature is found to depend on both the optical thickness and the dimensionless heat loss coefficient from the front plate.     

Free- and forced-convection film condensation from a horizontal elliptic tube with a vertical plate and horizontal tube as special cases

A simple mathematical model is developed for the study of the mixed-convection film condensation with downward flowing vapors onto a horizontal elliptic tube. Analytical analysis for both the local condensate film thickness and heat transfer characteristics under simultaneous effects of interfacial vapor shear and pressure gradient has been performed by adopting a unified geometry parameter, eccentricity e. The present results for two limit cases, e = 0 (circular tube) and e = 1.0 (vertical plate) are in an excellent agreement with the earlier works. For very slow vapor flow, the present result for dimensionless mean heat transfer coefficient reduces to the same form as in the earlier works, , whose value is 0.728 for e = 0 and 0.943 for e = 1.0. As for very fast vapor flow, the dimensionless mean heat transfer coefficient, increase with increasing eccentricity under the effects of pressure gradient caused by potential flow and surface tension.