Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Dynamics of convectively driven banded jets in the laboratory

The banded organization of clouds and zonal winds in the atmospheres of the outer planets has long fascinated observers. Several recent studies in the theory and idealized modeling of geostrophic turbulence have suggested possible explanations for the emergence of such organized patterns, typically involving highly anisotropic exchanges of kinetic energy and vorticity within the dissipationless inertial ranges of turbulent flows dominated (at least at large scales) by ensembles of propagating Rossby waves. The results from an attempt to reproduce such conditions in the laboratory are presented here. Achievement of a distinct inertial range turns out to require an experiment on the largest feasible scale. Deep, rotating convection on small horizontal scales was induced by gently and continuously spraying dense, salty water onto the free surface of the 13-m-diameter cylindrical tank on the Coriolis platform in Grenoble, France. A “planetary vorticity gradient” or “β effect” was obtained by use of a conically sloping bottom and the whole tank rotated at angular speeds up to 0.15 rad s−1. Over a period of several hours, a highly barotropic, zonally banded large-scale flow pattern was seen to emerge with up to 5–6 narrow, alternating, zonally aligned jets across the tank, indicating the development of an anisotropic field of geostrophic turbulence. Using particle image velocimetry (PIV) techniques, zonal jets are shown to have arisen from nonlinear interactions between barotropic eddies on a scale comparable to either a Rhines or “frictional” wavelength, which scales roughly as (β/Urms)−1/2. This resulted in an anisotropic kinetic energy spectrum with a significantly steeper slope with wavenumber k for the zonal flow than for the nonzonal eddies, which largely follows the classical Kolmogorov k−5/3 inertial range. Potential vorticity fields show evidence of Rossby wave breaking and the presence of a “hyperstaircase” with radius, indicating instantaneous flows that are supercritical with respect to the Rayleigh–Kuo instability criterion and in a state of “barotropic adjustment.” The implications of these results are discussed in light of zonal jets observed in planetary atmospheres and, most recently, in the terrestrial oceans.

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

The nature of Arctic polar vortices in chemistry-climate models

The structure of the Arctic stratospheric polar vortex in three chemistry–climate models (CCMs) taken from the CCMVal-2 intercomparison is examined using zonal mean and geometric-based methods. The geometric methods are employed by taking 2D moments of potential vorticity fields that are representative of the polar vortices in each of the models. This allows the vortex area, centroid location and ellipticity to be determined, as well as a measure of vortex filamentation. The first part of the study uses these diagnostics to examine how well the mean state, variability and extreme variability of the polar vortices are represented in CCMs compared to ERA-40 reanalysis data, and in particular for the UMUKCA-METO, NIWA-SOCOL and CCSR/NIES models. The second part of the study assesses how the vortices are predicted to change in terms of the frequency of sudden stratospheric warmings and their general structure over the period 1960–2100. In general, it is found that the vortices are climatologically too far poleward in the CCMs and produce too few large-scale filamentation events. Only a small increase is observed in the frequency of sudden stratospheric warming events from the mean of the CCMVal-2 models, but the distribution of extreme variability throughout the winter period is shown to change towards the end of the twentyfirst century.

A comparison of two identification and tracking methods for polar lows

In this study, we compare two different cyclone-tracking algorithms to detect North Atlantic polar lows, which are very intense mesoscale cyclones. Both approaches include spatial filtering, detection, tracking and constraints specific to polar lows. The first method uses digital bandpass-filtered mean sea level pressure (MSLP) fieldsin the spatial range of 200�600 km and is especially designed for polar lows. The second method also uses a bandpass filter but is based on the discrete cosine transforms (DCT) and can be applied to MSLP and vorticity fields. The latter was originally designed for cyclones in general and has been adapted to polar lows for this study. Both algorithms are applied to the same regional climate model output fields from October 1993 to September 1995 produced from dynamical downscaling of the NCEP/NCAR reanalysis data. Comparisons between these two methods show that different filters lead to different numbers and locations of tracks. The DCT is more precise in scale separation than the digital filter and the results of this study suggest that it is more suited for the bandpass filtering of MSLP fields. The detection and tracking parts also influence the numbers of tracks although less critically. After a selection process that applies criteria to identify tracks of potential polar lows, differences between both methods are still visible though the major systems are identified in both.

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.

Potential vorticity and moisture in extratropical cyclones : climatology and sensitivity experiments

The development of extratropical cyclones can be seen as an interplay of three positive potential vorticity (PV) anomalies: an upper-level stratospheric intrusion, low-tropospheric diabatically produced PV, and a warm anomaly at the surface acting as a surrogate PV anomaly. In the mature stage they become vertically aligned and form a “PV tower” associated with strong cyclonic circulation. This paradigm of extratropical cyclone development provides the basis of this thesis, which will use a climatological dataset and numerical model experiments to investigate the amplitude of the three anomalies and the processes leading in particular to the formation of the diabatically produced low-tropospheric PV anomaly.rnrnThe first part of this study, based on the interim ECMWF Re-Analysis (ERA-Interim) dataset, quantifies the amplitude of the three PV anomalies in mature extratropical cyclones in different regions in the Northern Hemisphere on a climatological basis. A tracking algorithm is applied to sea level pressure (SLP) fields to identify cyclone tracks. Surface potential temperature anomalies ∆θ and vertical profiles of PV anomalies ∆PV are calculated at the time of the cyclones’ minimum SLP and during the intensification phase 24 hours before in a vertical cylinder with a radius of 200 km around the surface cyclone center. To compare the characteristics of the cyclones, they are grouped according to their location (8 regions) and intensity, where the central SLP is used as a measure of intensity. Composites of ∆PV profiles and ∆θ are calculated for each region and intensity class at the time of minimum SLP and during the cyclone intensification phase.rnrnDuring the cyclones’ development stage the amplitudes of all three anomalies increase on average. In the mature stage all three anomalies are typically larger for intense than for weak winter cyclones [e.g., 0.6 versus 0.2 potential vorticity units (PVU) at lower levels, and 1.5 versus 0.5 PVU at upper levels].rnThe regional variability of the cyclones’ vertical structure and the profile evolution is prominent (cyclones in some regions are more sensitive to the amplitude of a particular anomaly than in other regions). Values of ∆θ and low-level ∆PV are on average larger in the western parts of the oceans than in the eastern parts. In addition, a large seasonal variability can be identified, with fewer and weaker cyclones especially in the summer, associated with higher low-tropospheric PV values, but also with a higher tropopause and much weaker surface potential temperature anomalies (compared to winter cyclones).rnrnIn the second part, we were interested in the diabatic low-level part of PV towers. Evaporative sources were identified of moisture that was involved in PV production through condensation. Lagrangian backward trajectories were calculated from the region with high PV values at low-levels in the cyclones. PV production regions were identified along these trajectories and from these regions a new set of backward trajectories was calculated and moisture uptakes were traced along them. The main contribution from surface evaporation to the specific humidity of the trajectories is collected 12-72 hours prior to therntime of PV production. The uptake region for weaker cyclones with less PV in the centre is typically more localized with reduced uptake values compared to intense cyclones. However, in a qualitative sense uptakes and other variables along single trajectories do not vary much between cyclones of different intensity in different regions.rnrnA sensitivity study with the COSMO model comprises the last part of this work. The study aims at investigating the influence of synthetic moisture modification in the cyclone environment in different stages of its development. Moisture was eliminated in three regions, which were identified as important moisture source regions for PV production. Moisture suppression affected the cyclone the most in its early phase. It led to cyclolysis shortly after its genesis. Nevertheles, a new cyclone formed on the other side of a dry box and developed relatively quickly. Also in other experiments, moisture elimination led to strong intensity reduction of the surface cyclone, limited upper-level development, and delayed or missing interaction between the two.rnrnIn summary, this thesis provides novel insight into the structure of different intensity categories of extratropical cyclones from a PV perspective, which corroborates the findings from a series of previous case studies. It reveals that all three PV anomalies are typically enhanced for more intense cyclones, with important regional differences concerning the relative amplitude of the three anomalies. The moisture source analysis is the first of this kind to study the evaporation-condensation cycle related to the intensification of extratropical cyclones. Interestingly, most of the evaporation occurs during the 3 days prior to the time of maximum cyclone intensity and typically extends over fairly large areas along the track of the cyclone. The numerical model case study complements this analysis by analyzing the impact of regionally confined moisture sources for the evolution of the cyclone.

Spatial pattern detection modeling of thrips (Thrips tabaci) on onion fields

Onion (Allium cepa) is one of the most cultivated and consumed vegetables in Brazil and its importance is due to the large laborforce involved. One of the main pests that affect this crop is the Onion Thrips (Thrips tabaci), but the spatial distribution of this insect, although important, has not been considered in crop management recommendations, experimental planning or sampling procedures. Our purpose here is to consider statistical tools to detect and model spatial patterns of the occurrence of the onion thrips. In order to characterize the spatial distribution pattern of the Onion Thrips a survey was carried out to record the number of insects in each development phase on onion plant leaves, on different dates and sample locations, in four rural properties with neighboring farms under different infestation levels and planting methods. The Mantel randomization test proved to be a useful tool to test for spatial correlation which, when detected, was described by a mixed spatial Poisson model with a geostatistical random component and parameters allowing for a characterization of the spatial pattern, as well as the production of prediction maps of susceptibility to levels of infestation throughout the area.

Magnetic fields and acute lymphoblastic leukemia in children: a systematic review of case-control studies

Leukemia incidence in children has increased worldwide in recent decades, particularly due to the rise in acute lymphoblastic leukemia. Studies have associated exposure to non-ionizing radiation generated by low frequency magnetic fields with childhood leukemia. The current article reviews the case-control studies published on this subject. Of 152 articles tracked in different databases, ten studies from North America, Asia, and Europe met the defined selection criteria, with patients diagnosed from 1960 to 2004. Methodological limitations were observed in these articles, including difficulties with the procedures for assessing exposure. An association may exist between exposure to low frequency magnetic fields and acute lymphoblastic leukemia in children, but this association is weak, preventing the observation of consistency in the findings. Future studies from a wider range of geographic regions should focus on the analysis of acute lymphoblastic leukemia, which is the subtype with the greatest impact on the increasing overall incidence of childhood leukemia.

Electronic properties and hyperfine fields of nickel-related complexes in diamond

We carried out a first-principles investigation on the microscopic properties of nickel-related defect centers in diamond. Several configurations, involving substitutional and interstitial nickel impurities, have been considered either in isolated configurations or forming complexes with other defects, such as vacancies and boron and nitrogen dopants. The results, in terms of spin, symmetry, and hyperfine fields, were compared with the available experimental data on electrically active centers in synthetic diamond. Several microscopic models, previously proposed to explain those data, have been confirmed by this investigation, while some models could be discarded. We also provided insights into the microscopic structure of several of those centers.

Origin of primordial magnetic fields

Magnetic fields of intensities similar to those in our galaxy are also observed in high redshift galaxies, where a mean field dynamo would not have had time to produce them. Therefore, a primordial origin is indicated. It has been suggested that magnetic fields were created at various primordial eras: during inflation, the electroweak phase transition, the quark-hadron phase transition (QHPT), during the formation of the first objects, and during reionization. We suggest here that the large-scale fields similar to mu G, observed in galaxies at both high and low redshifts by Faraday rotation measurements (FRMs), have their origin in the electromagnetic fluctuations that naturally occurred in the dense hot plasma that existed just after the QHPT. We evolve the predicted fields to the present time. The size of the region containing a coherent magnetic field increased due to the fusion of smaller regions. Magnetic fields (MFs) similar to 10 mu G over a comoving similar to 1 pc region are predicted at redshift z similar to 10. These fields are orders of magnitude greater than those predicted in previous scenarios for creating primordial magnetic fields. Line-of-sight average MFs similar to 10(-2) mu G, valid for FRMs, are obtained over a 1 Mpc comoving region at the redshift z similar to 10. In the collapse to a galaxy (comoving size similar to 30 kpc) at z similar to 10, the fields are amplified to similar to 10 mu G. This indicates that the MFs created immediately after the QHPT (10(-4) s), predicted by the fluctuation-dissipation theorem, could be the origin of the similar to mu G fields observed by FRMs in galaxies at both high and low redshifts. Our predicted MFs are shown to be consistent with present observations. We discuss the possibility that the predicted MFs could cause non-negligible deflections of ultrahigh energy cosmic rays and help create the observed isotropic distribution of their incoming directions. We also discuss the importance of the volume average magnetic field predicted by our model in producing the first stars and in reionizing the Universe.

Origin of magnetic fields in galaxies

Microgauss magnetic fields are observed in all galaxies at low and high redshifts. The origin of these intense magnetic fields is a challenging question in astrophysics. We show here that the natural plasma fluctuations in the primordial Universe (assumed to be random), predicted by the fluctuation - dissipation theorem, predicts similar to 0.034 mu G fields over similar to 0.3 kpc regions in galaxies. If the dipole magnetic fields predicted by the fluctuation- dissipation theorem are not completely random, microgauss fields over regions greater than or similar to 0.34 kpc are easily obtained. The model is thus a strong candidate for resolving the problem of the origin of magnetic fields in less than or similar to 10(9) years in high redshift galaxies.

Automated supervised classification of variable stars in the CoRoT programme Method and application to the first four exoplanet fields

Aims. In this work, we describe the pipeline for the fast supervised classification of light curves observed by the CoRoT exoplanet CCDs. We present the classification results obtained for the first four measured fields, which represent a one-year in-orbit operation. Methods. The basis of the adopted supervised classification methodology has been described in detail in a previous paper, as is its application to the OGLE database. Here, we present the modifications of the algorithms and of the training set to optimize the performance when applied to the CoRoT data. Results. Classification results are presented for the observed fields IRa01, SRc01, LRc01, and LRa01 of the CoRoT mission. Statistics on the number of variables and the number of objects per class are given and typical light curves of high-probability candidates are shown. We also report on new stellar variability types discovered in the CoRoT data. The full classification results are publicly available.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.