Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Microphysical evidence of the transition between predominant convective/stratiform rainfall associated with the intraseasonal oscillation in the Southwest Amazon

The distinction between convective and stratiform precipitation profiles around various precipitating systems existent in tropical regions is very important to the global atmospheric circulation, which is extremely sensitive to vertical latent heat distribution. In South America, the convective activity responds to the Intraseasonal Oscillation (IOS). This paper analyzes a disdrometer and a radar profiler data, installed in the Ji-Paraná airport, RO, Brazil, for the field experiment WETAMC/LBA & TRMM/LBA, during January and February of 1999. The microphysical analysis of wind regimes associated with IOS showed a large difference in type, size and microphysical processes of hydrometeor growth in each wind regime: easterly regimes had more turbulence and consequently convective precipitation formation, and westerly regimes had a more stratiform precipitation formation.

O uso de cálculos teóricos como ferramenta para a detecção de moléculas no meio interestelar: o caso dos carbenos CnNH (n=1, 3 e 5)

This work reviews some applications of ab initio molecular orbital calculations to the elucidation of structures of interstellar molecules. The case of the CnNH (n=1, 3 and 5 ) carbenes is extensively analyzed and discussed. Theoretical conformational analysis and predicted values for the rotational constants and dipole moments of the singlet state of C5NH are reported for the first time and a comparison is performed with results previously obtained for C3NH and CNH.

Uma revisão sobre o efeito Toms: o fenômeno onde macromoléculas atenuam a turbulência em um líquido

It is very well known that the addition of polymers to a liquid increases the shear viscosity of the solution. In other words, the polymer increases the dissipation of the flow energy. Contrarily, in turbulent flow, some particular macromolecules in very low concentration are able to produce large attenuation in the turbulence and thus, decreasing the dissipation of the energy. This article present a brief revision about macroscopic and molecular models used to explain this dynamic effect. Some of the experimental techniques used to quantify the attenuation of the turbulence and the main active substances are also discussed.

Chumbo e zinco em águas e sedimentos de área de mineração e metalurgia de metais

Lead metallurgy at Adrianópolis is the largest environmental problem of Paraná, Brazil. The objective of this study was to determine Pb and Zn levels in water and sediment in two catchments by different extraction methods. The high levels of lead in water in most of samples do not allow the human use. Total Pb concentration as high as 795.3 µg L-1 was observed in Ribeira River bank, in a pluvial water stream flowering from a abandoned factory. Due to the high Pb levels in sediments from some sites (maximum of 24,300 mg kg-1) is recommended to avoid the water turbulence.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

The pulsed wire anemometer: review and further developments

The purpose of the present paper is to review work that has been done on the pulsed wire anemometer technique and also suggest further developments that could be made in its range of application. The aper discusses the three types of probes that have been used in pulsed wire anemometry: the crossed wire velocity probe, the parallel wire wall shear stress probe and the parallel wire velocity probe. The work shows that the crossed wire and the parallel wire techniques can be used to make velocity, turbulence and wall shear stress measurements in highly turbulent flows without any upper restriction on turbulence level. Comments are also made on the potential of a parallel wire probe for use in highly turbulent flows that would enable higher order velocity cross-product terms to be measured.

Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.

A law of the wall formulation for recirculating flows

This work presents a new law of the wall formulation for recirculating turbulent flows. An alternative expression for the internal length which can be applied in the separated region is also presented. The formulation is implemented in a numerical code which solves the k-epsilon model through a finite volume method. The theoretical results are compared with the experimental data of Vogel and Eaton (J. of Heat Transfer, Transactions of ASME, vol.107, pp. 922-929, 1985). The paper shows that the present formulation furnishes better results than the standard k-epsilon formulation.

Dynamics of coherent vortices in mixing layers using direct numerical and large-eddy simulations

Coherent vortices in turbulent mixing layers are investigated by means of Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES). Subgrid-scale models defined in spectral and physical spaces are reviewed. The new "spectral-dynamic viscosity model", that allows to account for non-developed turbulence in the subgrid-scales, is discussed. Pseudo-spectral methods, combined with sixth-order compact finite differences schemes (when periodic boundary conditions cannot be established), are used to solve the Navier- Stokes equations. Simulations in temporal and spatial mixing layers show two types of pairing of primary Kelvin-Helmholtz (KH) vortices depending on initial conditions (or upstream conditions): quasi-2D and helical pairings. In both cases, secondary streamwise vortices are stretched in between the KH vortices at an angle of 45° with the horizontal plane. These streamwise vortices are not only identified in the early transitional stage of the mixing layer but also in self-similar turbulence conditions. The Re dependence of the "diameter" of these vortices is analyzed. Results obtained in spatial growing mixing layers show some evidences of pairing of secondary vortices; after a pairing of the primary Kelvin-Helmholtz (KH) vortices, the streamwise vortices are less numerous and their diameter has increased than before the pairing of KH vortices.

An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

Instabilities in electrochemical systems with a rotating disc electrode

Polarization curves experimentally obtained in the electro-dissolution of iron in a 1 M H2SO4 solution using a rotating disc as the working electrode present a current instability region within the range of applied voltage in which the current is controlled by mass transport in the electrolyte. According to the literature (Barcia et. al., 1992) the electro-dissolution process leads to the existence of a viscosity gradient in the interface metal-solution, which leads to a velocity field quantitatively different form the one developed in uniform viscosity conditions and may affect the stability of the hydrodynamic field. The purpose of this work is to investigate whether a steady viscosity profile, depending on the distance to the electrode surface, affects the stability properties of the classic velocity field near a rotating disc. Two classes of perturbations are considered: perturbations monotonically varying along the radial direction, and perturbations periodically modulated along the radial direction. The results show that the hydrodynamic field is always stable with respect to the first class of perturbations and that the neutral stability curves are modified by the presence of a viscosity gradient in the second case, in the sense of reducing the critical Reynolds number beyond which perturbations are amplified. This result supports the hypothesis that the current oscillations observed in the polarization curve may originate from a hydrodynamic instability.

Wind tunnel simulation of atmospheric boundary layer flows

The present work shows how thick boundary layers can be produced in a short wind tunnel with a view to simulate atmospheric flows. Several types of thickening devices are analysed. The experimental assessment of the devices was conducted by considering integral properties of the flow and the spectra: skin-friction, mean velocity profiles in inner and outer co-ordinates and longitudinal turbulence. Designs based on screens, elliptic wedge generators, and cylindrical rod generators are analysed. The paper describes in detail the experimental arrangement, including the features of the wind tunnel and of the instrumentation. The results are compared with experimental data published by other authors and with naturally developed flows.