Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

The Sphene-Centered Ocellar Texture: An Effect of Grain-Supported Flow and Melt Migration in a Hyperdense Magma Mush

The sphene-centered ocellar texture consists of leucocratic ocelli with sphene (titanite) crystals at the center, enclosed in a biotite-rich matrix. This texture has been recognized worldwide in hybrid intermediate rocks. On the basis of structural, petrological, and geochronological data from selected outcrops of the Variscan Ribadelago pluton (NW Iberian Massif), we propose that the ocelli were formed by migration and accumulation of a residual melt through a plagioclase- and biotite-dominated crystalline framework. At the late stage of crystallization, the magma acted as a hyperdense suspension and reacted to the pressure gradient caused by the regional stress field, entering the domain of grain-supported flow. Microstructures reveal that aligned crystal domains arose in the crystal framework from the shearing and compaction of the crystal mush and behaved as magmatic microshears. Relative displacement of adjacent crystal clusters along these microshears corresponded to the onset of Reynolds dilatancy that generated an expansion of the crystal mush, involving melt migration and pore aperture. The mineralogy of the ocelli, dominated by andesine and sphene, represents the composition of the migrating melt. The chemistry of this late, Ti-rich melt stems from the incongruent melting of biotite. Magmatic sphene from the ocelli yields a U-Pb age of 317 +/- 1 Ma, which represents the final crystallization of the hybridized magmatic system. Moreover, this texture offers an opportunity to better understand the rheological behavior of highly crystallized magmas.

Estudo numérico de escoamento turbulento em padrão anular gás-líquido em dutos verticais

Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm

Pressure drop correlation for oil-refrigerant R134a mixture flashing flow in a small diameter tube

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Water-assisted Flow of Heavy Oil in a Vertical Pipe: Pilot-scale Experiments

The use of the core-annular flow pattern, where a thin fluid surrounds a very viscous one, has been suggested as an attractive artificial-lift method for heavy oils in the current Brazilian ultra-deepwater production scenario. This paper reports the pressure drop measurements and the core-annular flow observed in a 2 7/8-inch and 300 meter deep pilot-scale well conveying a mixture of heavy crude oil (2000 mPa.s and 950 kg/m3 at 35 C) and water at several combinations of the individual flow rates. The two-phase pressure drop data are compared with those of single-phase oil flow to assess the gains due to water injection. Another issue is the handling of the core-annular flow once it has been established. High-frequency pressure-gradient signals were collected and a treatment based on the Gabor transform together with neural networks is proposed as a promising solution for monitoring and control. The preliminary results are encouraging. The pilot-scale tests, including long-term experiments, were conducted in order to investigate the applicability of using water to transport heavy oils in actual wells. It represents an important step towards the full scale application of the proposed artificial-lift technology. The registered improvements in terms of oil production rate and pressure drop reductions are remarkable.

Quasi-cylindrical approximation to the swirling flow in an atomizer chamber

A quasi-cylindrical approximation is used to analyse the axisymmetric swirling flow of a liquid with a hollow air core in the chamber of a pressure swirl atomizer. The liquid is injected into the chamber with an azimuthal velocity component through a number of slots at the periphery of one end of the chamber, and flows out as an anular sheet through a central orifice at the other end, following a conical convergence of the chamber wall. An effective inlet condition is used to model the effects of the slots and the boundary layer that develops at the nearby endwall of the chamber. An analysis is presented of the structure of the liquid sheet at the end of the exit orifice, where the flow becomes critical in the sense that upstream propagation of long-wave perturbations ceases to be possible. This nalysis leads to a boundary condition at the end of the orifice that is an extension of the condition of maximum flux used with irrotational models of the flow. As is well known, the radial pressure gradient induced by the swirling flow in the bulk of the chamber causes the overpressure that drives the liquid towards the exit orifice, and also leads to Ekman pumping in the boundary layers of reduced azimuthal velocity at the convergent wall of the chamber and at the wall opposite to the exit orifice. The numerical results confirm the important role played by the boundary layers. They make the thickness of the liquid sheet at the end of the orifice larger than predicted by rrotational models, and at the same time tend to decrease the overpressure required to pass a given flow rate through the chamber, because the large axial velocity in the boundary layers takes care of part of the flow rate. The thickness of the boundary layers increases when the atomizer constant (the inverse of a swirl number, proportional to the flow rate scaled with the radius of the exit orifice and the circulation around the air core) decreases. A minimum value of this parameter is found below which the layer of reduced azimuthal velocity around the air core prevents the pressure from increasing and steadily driving the flow through the exit orifice. The effects of other parameters not accounted for by irrotational models are also analysed in terms of their influence on the boundary layers.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Three dimensional models of water flow in the Persian Gulf

Persian Gulf region is globally of great importance due to its economical and political reasons. The importance lies in oil sources and sea exports. Geophysical phenomena dominated in the water circulation affected this region is called Monsoon it stretches from African coasts to the half way of Red Seal affected all coasts of Persian Gulf and goes toward east to the Indian ocean. Other essential factors in the water circulation in this region are net evaporation (several meters in per year), high density and high salinity. In this article the effects of wind stress and evaporation in the water circulation in the region will be considered and model equations for wind forces, density, pressure, gradient, and bottom friction for Persian Gulf will be discussed.

CFD convective flow simulation of the varying properties of CO2-H20 mixtures in geothermal systems

Numerical simulation of a geothermal reservoir, modelled as a bottom-heated square box, filled with water-CO2 mixture is presented in this work. Furthermore, results for two limiting cases of a reservoir filled with either pure water or CO2 are presented. Effects of different parameters including CO2 concentration as well as reservoir pressure and temperature on the overall performance of the system are investigated. It has been noted that, with a fixed reservoir pressure and temperature, any increase in CO2concentration leads to better performance, i.e. stronger convection and higher heat transfer rates. With a fixed CO2 concentration, however, the reservoir pressure and temperature can significantly affect the overall heat transfer and flow rate from the reservoir. Details of such variations are documented and discussed in the present paper.

Turbulent spot growth in favorable pressure gradients

Experimental results are presented on the lateral growth of turbulent spots in a series of flows with favorable pressure gradients. It is shown that the wedge angle increases slowly with the Reynolds number and that a favorable pressure gradient inhibits the growth of turbulent spots and, in general, results in a nonlinear turbulent wedge. As soon as the pressure gradient decreases to the point where the flow becomes supercritical, however, spot growth increases rapidly and the associated turbulent wedge becomes linear.

Laminar compressible boundary layers with vectored mass transfer

The paper studies the influence of vectored suction or injection on the flow and heat transfer at the stagnation point of a two-dimensional body (a cylinder) and an axisymmetric body (a sphere) with allowance for the effects of variable gas properties. The analysis is based on the boundary-layer equations in dimensionless form for the steady compressible fluid with variable properties in the stagnation region of a two-dimensional or an axisymmetric body with tangential and normal surface mass transfer under similarity requirements. It is shown that the variation of the density-viscosity product across the boundary layer has a strong effect on the skin friction and heat transfer. This gives rise to a point of inflection which can be removed by suction and by increasing the wall temperature. The skin friction and heat transfer are significantly affected by the pressure gradient parameter.

Nonsimilar incompressible laminar boundary layers with magnetic field.

The solution of the steady laminar incompressible nonsimilar magneto-hydrodynamic boundary layer flow and heat transfer problem with viscous dissipation for electrically conducting fluids over two-dimensional and axisymmetric bodies with pressure gradient and magnetic field has been presented. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for flow over a cylinder and a sphere. The results indicate that the magnetic field tends to delay or prevent separation. The heat transfer strongly depends on the viscous dissipation parameter. When the dissipation parameter is positive (i.e. when the temperature of the wall is greater than the freestream temperature) and exceeds a certain value, the hot wall ceases to be cooled by the stream of cooler air because the ‘heat cushion’ provided by the frictional heat prevents cooling whereas the effect of the magnetic field is to remove the ‘heat cushion’ so that the wall continues to be cooled. The results are found to be in good agreement with those of the local similarity and local nonsimilarity methods except near the point of separation, but they are in excellent agreement with those of the difference-differential technique even near the point of separation.

On the ignition and extinction problems in forced convection systems

In this paper the problem of ignition and extinction has been formulated for the flow of a compressible fluid with Prandtl and Schmidt numbers taken as unity. In particular, the problems of (i) a jet impinging on a wall of combustible material and (ii) the opposed jet diffusion flame have been studied. In the wall jet case, three approximations in the momentum equation namely, (i) potential flow, (ii) viscous flow, (ii) viscous incompressible with k = 1 and (iii) Lees' approximation (taking pressure gradient terms zero) are studied. It is shown that the predictions of the mass flow rates at extinction are not very sensitive to the approximations made in the momentum equation. The effects of varying the wall temperature in the case (i) and the jet temperature in the case (ii) on the extinction speeds have been studied. The effects of varying the activation energy and the free stream oxidant concentration in case (ii), have also been investigated.

Studies on helium flow in experimental loops with different impedances operated with fountain effect pumps

An experimental flow loop with He II flow driven by fountain effect pumps (FEPs) is studied with respect to operation at different flow impedances and with thermal loads applied at different positions. The measured values of temperature, flow rate and pressure drop are compared with calculations resulting from a simplified model which assumes ideal performance of the porous plug and of the heat exchangers and which does not take into account Gorter-Mellink (GM) conduction. The main features of the loop are shown to be well described by this model. Refined calculations with a more complex model, including GM conduction of the He II, are only required for predicting the temperature distribution in some discrete regions of the loop.

Measurement of temperature and pressure on the surface of a blunt cone using FBG sensor in hypersonic wind tunnel

Measurement of temperature and pressure exerted on the leeward surface of a blunt cone specimen has been demonstrated in the present work in a hypersonic wind tunnel using fiber Bragg grating (FBG) sensors. The experiments were conducted on a 30 degrees apex-angle blunt cone with 51 mm base diameter at wind flow speeds of Mach 6.5 and 8.35 in a 300 mm hypersonic wind tunnel of Indian Institute of Science, Bangalore. A special pressure insensitive temperature sensor probe along with the conventional bare FBG sensors was used for explicit temperature and aerodynamic pressure measurement respectively on the leeward surface of the specimen. computational fluid dynamics (CFD) simulation of the flow field around the blunt cone specimen has also been carried out to obtain the temperature and pressure at conditions analogous to experiments. The results obtained from FBG sensors and the CFD simulations are found to be in good agreement with each other.