Processes controlling atmospheric dispersion through city centres

We develop a process-based model for the dispersion of a passive scalar in the turbulent flow around the buildings of a city centre. The street network model is based on dividing the airspace of the streets and intersections into boxes, within which the turbulence renders the air well mixed. Mean flow advection through the network of street and intersection boxes then mediates further lateral dispersion. At the same time turbulent mixing in the vertical detrains scalar from the streets and intersections into the turbulent boundary layer above the buildings. When the geometry is regular, the street network model has an analytical solution that describes the variation in concentration in a near-field downwind of a single source, where the majority of scalar lies below roof level. The power of the analytical solution is that it demonstrates how the concentration is determined by only three parameters. The plume direction parameter describes the branching of scalar at the street intersections and hence determines the direction of the plume centreline, which may be very different from the above-roof wind direction. The transmission parameter determines the distance travelled before the majority of scalar is detrained into the atmospheric boundary layer above roof level and conventional atmospheric turbulence takes over as the dominant mixing process. Finally, a normalised source strength multiplies this pattern of concentration. This analytical solution converges to a Gaussian plume after a large number of intersections have been traversed, providing theoretical justification for previous studies that have developed empirical fits to Gaussian plume models. The analytical solution is shown to compare well with very high-resolution simulations and with wind tunnel experiments, although re-entrainment of scalar previously detrained into the boundary layer above roofs, which is not accounted for in the analytical solution, is shown to become an important process further downwind from the source.

Secondary cyclogenesis along an occluded front leading to damaging wind gusts: windstorm Kyrill, January 2007

Windstorm Kyrill affected large parts of Europe in January 2007 and caused widespread havoc and loss of life. In this study the formation of a secondary cyclone, Kyill II, along the occluded front of the mature cyclone Kyrill and the occurrence of severe wind gusts as Kyrill II passed over Germany are investigated with the help of high-resolution regional climate model simulations. Kyrill underwent an explosive cyclogenesis south of Greenland as the storm crossed polewards of an intense upper-level jet stream. Later in its life cycle secondary cyclogenesis occurred just west of the British Isles. The formation of Kyrill II along the occluded front was associated (a) with frontolytic strain and (b) with strong diabatic heating in combination with a developing upper-level shortwave trough. Sensitivity studies with reduced latent heat release feature a similar development but a weaker secondary cyclone, revealing the importance of diabatic processes during the formation of Kyrill II. Kyrill II moved further towards Europe and its development was favored by a split jet structure aloft, which maintained the cyclone’s exceptionally deep core pressure (below 965 hPa) for at least 36 hours. The occurrence of hurricane force winds related to the strong cold front over North and Central Germany is analyzed using convection-permitting simulations. The lower troposphere exhibits conditional instability, a turbulent flow and evaporative cooling. Simulation at high spatio-temporal resolution suggests that the downward mixing of high momentum (the wind speed at 875 hPa widely exceeded 45 m s-1) accounts for widespread severe surface wind gusts, which is in agreement with observed widespread losses.

Changes on emitter discharge under different water temperature and pressure

The study was conducted at the Research Laboratory of Hydraulic and Irrigation Group in the Rural Engineering Department, Technical University of Madrid (Universidad Politecnica de Madrid), Madrid, Spain. Water temperatures of 20, 30, 40 degrees C and system pressures often encountered in irrigation practices of 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190 and 200 k Pa were applied to determine the effects of different water temperatures and pressures on emitter discharge. Non-pressure compensating in-line emitter which has turbulent flow regime with a long-path (labyrinth), emitter discharge was 4 L h(-1) at system pressure of 100 kPa according to the manufacturer recommended, was used. Emitters were spaced 20 cm along the drip laterals with 16 mm diameter. Discharge equations and coefficients of variation related to temperatures of 20, 30 and 40 degrees C were obtained as q = 0.375H(0.51), q = 0.358H(0.52), q = 0.346H(0.53) and 2.68, 2.09, 3.65, respectively. Discharge of the emitter was affected by different system pressures and increased as potentially (R = 0993-0996). In general. the emitter discharge increased with increasing temperature. However, especially in the common system pressures of 90-120 k Pa, differences of obtained emitter discharges between the different water temperatures were not significant (1%).

Desgaste corrosivo-cavitativo-erosivo de um aço-carbono em meio aquoso com frações de sal (NaCl), CO2 e particulados sólidos (SiO2)

A batch of eighty-four coupons of low carbon steel were investigated at laboratory conditions under a corrosive, cavitative-corrosive (CO2) and corrosive-erosive (SiO2 + CO2) in an aqueous salt solution and two levels of temperature. The following measurements were made on Vickers (HV0,05, HV0,10, HV0,20) Microhardness tests at three levels of subsurface layer. A turbulent flow collided on the cylindrical sample, with and without mechanical stirring and gas bubbling, with and without fluid contamination by solid particles of SiO2, at two temperatures. Surface Roughness and Waviness, under two conditions "as received, after machining" and "after worn out", as well as gravimetric and electrochemical parameter were measured on the two opposite generatrices of each cylindrical sample, on the flow upstream (0°) and downstream (180°) by Profilometry, Mass Variation and Linear Polarization Resistance (LPR). The results of the Microhardness and Surface Texture of all coupons were subjected to statistical comparison, using the software Statgraphics® Centurion XVI, 95% statistical certainty, and significant differences were observed in some arrays of measurements. The corrosive wear rate measured by LPR and mass variation shown to be sensitive to the presence of bubbles and hydrodynamic fluctuations inside the cell, considering the temperature and contamination of corrosive fluid by solid particles. The main results of visual inspection relative to some topologies of the surface damages involving different mechanisms that were seen to give explanation for some fluctuations in wear rates of the steel experimentally investigated

Modelagem e simulação da deposição de parafinas em escoamento turbulento

The flow assurance has become one of the topics of greatest interest in the oil industry, mainly due to production and transportation of oil in regions with extreme temperature and pressure. In these operations the wax deposition is a commonly problem in flow of paraffinic oils, causing the rising costs of the process, due to increased energy cost of pumping, decreased production, increased pressure on the line and risk of blockage of the pipeline. In order to describe the behavior of the wax deposition phenomena in turbulent flow of paraffinic oils, under different operations conditions, in this work we developed a simulator with easy interface. For that we divided de work in four steps: (i) properties estimation (physical, thermals, of transport and thermodynamics) of n-alkanes and paraffinic mixtures by using correlations; (ii) obtainment of the solubility curve and determination the wax appearance temperature, by calculating the solid-liquid equilibrium of parafinnic systems; (iii) modelling wax deposition process, comprising momentum, mass and heat transfer; (iv) development of graphic interface in MATLAB® environment for to allow the understanding of simulation in different flow conditions as well as understand the matter of the variables (inlet temperature, external temperature, wax appearance temperature, oil composition, and time) on the behavior of the deposition process. The results showed that the simulator developed, called DepoSim, is able to calculate the profile of temperature, thickness of the deposit, and the amount of wax deposited in a simple and fast way, and also with consistent results and applicable to the operation

Processamento e interpretação de dados 2D e 3D de GPR :aplicações no imageamento de feições kársticas e estruturas de dissolução no campo de petróleo de Fazenda Belém-CE

In Fazenda Belém oil field (Potiguar Basin, Ceará State, Brazil) occur frequently sinkholes and sudden terrain collapses associated to an unconsolidated sedimentary cap covering the Jandaíra karst. This research was carried out in order to understand the mechanisms of generation of these collapses. The main tool used was Ground Penetrating Radar (GPR). This work is developed twofold: one aspect concerns methodology improvements in GPR data processing whilst another aspect concerns the geological study of the Jandaíra karst. This second aspect was strongly supported both by the analysis of outcropping karst structures (in another regions of Potiguar Basin) and by the interpretation of radargrams from the subsurface karst in Fazenda Belém. It was designed and tested an adequate flux to process GPR data which was adapted from an usual flux to process seismic data. The changes were introduced to take into account important differences between GPR and Reflection Seismic methods, in particular: poor coupling between source and ground, mixed phase of the wavelet, low signal-to-noise ratio, monochannel acquisition, and high influence of wave propagation effects, notably dispersion. High frequency components of the GPR pulse suffer more pronounced effects of attenuation than low frequency components resulting in resolution losses in radargrams. In Fazenda Belém, there is a stronger need of an suitable flux to process GPR data because both the presence of a very high level of aerial events and the complexity of the imaged subsurface karst structures. The key point of the processing flux was an improvement in the correction of the attenuation effects on the GPR pulse based on their influence on the amplitude and phase spectra of GPR signals. In low and moderate losses dielectric media the propagated signal suffers significant changes only in its amplitude spectrum; that is, the phase spectrum of the propagated signal remains practically unaltered for the usual travel time ranges. Based on this fact, it is shown using real data that the judicious application of the well known tools of time gain and spectral balancing can efficiently correct the attenuation effects. The proposed approach can be applied in heterogeneous media and it does not require the precise knowledge of the attenuation parameters of the media. As an additional benefit, the judicious application of spectral balancing promotes a partial deconvolution of the data without changing its phase. In other words, the spectral balancing acts in a similar way to a zero phase deconvolution. In GPR data the resolution increase obtained with spectral balancing is greater than those obtained with spike and predictive deconvolutions. The evolution of the Jandaíra karst in Potiguar Basin is associated to at least three events of subaerial exposition of the carbonatic plataform during the Turonian, Santonian, and Campanian. In Fazenda Belém region, during the mid Miocene, the Jandaíra karst was covered by continental siliciclastic sediments. These sediments partially filled the void space associated to the dissolution structures and fractures. Therefore, the development of the karst in this region was attenuated in comparison to other places in Potiguar Basin where this karst is exposed. In Fazenda Belém, the generation of sinkholes and terrain collapses are controlled mainly by: (i) the presence of an unconsolidated sedimentary cap which is thick enough to cover completely the karst but with sediment volume lower than the available space associated to the dissolution structures in the karst; (ii) the existence of important structural of SW-NE and NW-SE alignments which promote a localized increase in the hydraulic connectivity allowing the channeling of underground water, thus facilitating the carbonatic dissolution; and (iii) the existence of a hydraulic barrier to the groundwater flow, associated to the Açu-4 Unity. The terrain collapse mechanisms in Fazenda Belém occur according to the following temporal evolution. The meteoric water infiltrates through the unconsolidated sedimentary cap and promotes its remobilization to the void space associated with the dissolution structures in Jandaíra Formation. This remobilization is initiated at the base of the sedimentary cap where the flow increases its abrasion due to a change from laminar to turbulent flow regime when the underground water flow reaches the open karst structures. The remobilized sediments progressively fill from bottom to top the void karst space. So, the void space is continuously migrated upwards ultimately reaching the surface and causing the sudden observed terrain collapses. This phenomenon is particularly active during the raining season, when the water table that normally is located in the karst may be temporarily located in the unconsolidated sedimentary cap

Friction factors, convective heat transfer coefficients and the Colburn analogy for industrial sugarcane juices

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Friction losses in valves and fittings for liquid food products

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

Transition to turbulence in the Reynolds' experiment

In this study we simulate numerically the Reynolds' experiment for the transition from laminar to turbulent flow in a pipe. We present a discussion of the results from a dynamical systems perspective when a control parameter, the Reynolds number, is increased. The Landau scenario, where the transition is described by the excitation of infinite oscillatory modes within the fluid, is not observed. Instead what happens is best explained by the Ruelle-Takens scenario in terms of strange attractors. The Lyapunov exponent and fractal dimension for the attractor are calculated together with a measure of complex behaviour called the Lempel-Ziv complexity. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Friction losses in valves and fittings for power-law fluids

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

The gradually truncated Levy flight for systems with power-law distributions

Power-law distributions have been observed in various economical and physical systems. Levy flights have infinite variance which discourage a physical approach. We introduce a class of stochastic processes, the gradually truncated Levy flight in which large steps of a Levy flight are gradually eliminated. It has finite variance and the system can be analyzed in a closed form. We applied the present method to explain the distribution of a particular economical index. The present method can be applied to describe time series in a variety of fields, i.e. turbulent flow, anomalous diffusion, polymers, etc. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Friction losses in valves and fittings for viscoplastic fluids

Data on pressure drop were obtained in stainless steel, sanitary fittings and valves during laminar and turbulent flow of aqueous suspensions of sucrose and bentonite. The rheological properties of these suspensions were determined and the Bingham model provided the best fitting with the experimental data. Friction losses were measured in fully- and partially-open butterfly and plug valves, bends and union. Values of loss coefficients (k(f)) were calculated and correlated as functions of the classical Reynolds number and the Reynolds number proposed by Govier and Aziz (1972) for viscoplastic fluids. The two-k method and a new proposed model presented the best adjustments for the Govier and Aziz Reynolds number, and Hedstrom and classical Reynolds numbers, respectively.

Continuous iteration of dynamical maps

A precise meaning is given to the notion of continuous iteration of a mapping. Usual discrete iterations are extended into a dynamical flow which is a homotopy of them all. The continuous iterate reveals that a dynamic map is formed by independent component modes evolving without interference with each other. An application to turbulent flow suggests that the velocity field assumes nonseparable values. © 1998 American Institute of Physics.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Multiloop controller and reference generator for a dynamic voltage restorer implementation

This paper deals with the design and analysis of a Dynamic Voltage Restorer output voltage control. Such control is based on a multiloop strategy, with an inner current PID regulator and an outer P+Resonant voltage controller. The inner regulator is applied on the output inductor current. It will be also demonstrated how the load current behavior may influence in the DVR output voltage, which justifies the need for the resonant controller. Additionally, it will be discussed the application of a modified algorithm for the identification of the DVR voltage references, which is based on a previously presented positive sequence detector. Since the studied three-phase DVR is assumed to be based on three identical H-bridge converters, all the analysis and design procedures were realized by means of single-phase equivalent circuits. The discussions and conclusions are supported by theoretical calculations, nonlinear simulations and some experimental results. ©2008 IEEE.