Universality of anisotropic fluctuations from numerical simulations of turbulent flows

We present new results from a direct numerical simulation of a three dimensional homogeneous Rayleigh-Benard system (HRB), i.e. a convective cell with an imposed linear mean temperature profile along the vertical direction. We measure the SO(3)-decomposition of both velocity structure functions and buoyancy terms. We give a dimensional prediction for the values of the anisotropic scaling exponents in this Rayleigh-Benard systems. Measured scaling does not follow dimensional estimate, while a better agreement can be found with the anisotropic scaling of a different system, the random-Kolmogorov-flow (RKF). Our findings support the conclusion that scaling properties of anisotropic fluctuations are universal, i.e. independent of the forcing mechanism sustaining the turbulent flow.Comment: 4 pages, 3 figure

Velocity gradients statistics along particle trajectories in turbulent flows: the refined similarity hypothesis in the Lagrangian frame

We present an investigation of the statistics of velocity gradient related quantities, in particluar energy dissipation rate and enstrophy, along the trajectories of fluid tracers and of heavy/light particles advected by a homogeneous and isotropic turbulent flow. The Refined Similarity Hypothesis (RSH) proposed by Kolmogorov and Oboukhov in 1962 is rephrased in the Lagrangian context and then tested along the particle trajectories. The study is performed on state-of-the-art numerical data resulting from numerical simulations up to Re~400 with 2048^3 collocation points. When particles have small inertia, we show that the Lagrangian formulation of the RSH is well verified for time lags larger than the typical response time of the particle. In contrast, in the large inertia limit when the particle response time approaches the integral-time-scale of the flow, particles behave nearly ballistic, and the Eulerian formulation of RSH holds in the inertial-range.Comment: 7 pages, 7 figures; Physical Review E (accepted Dec 7, 2009

Non-Oberbeck-Boussinesq effects in two-dimensional Rayleigh-Benard convection in glycerol

We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of "non-Oberbeck-Boussinesqness'' or "NOBness'' ($\Delta$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$.Comment: 6 pages, 7 figure

Copepods encounter rates from a model of escape jump behaviour in turbulence

A key ecological parameter for planktonic copepods studies is their interspecies encounter rate which is driven by their behaviour and is strongly influenced by turbulence of the surrounding environment. A distinctive feature of copepods motility is their ability to perform quick displacements, often dubbed jumps, by means of powerful swimming strokes. Such a reaction has been associated to an escape behaviour from flow disturbances due to predators or other external dangers. In the present study, the encounter rate of copepods in a developed turbulent flow with intensity comparable to the one found in copepods' habitat is numerically investigated. This is done by means of a Lagrangian copepod (LC) model that mimics the jump escape reaction behaviour from localised high-shear rate fluctuations in the turbulent flows. Our analysis shows that the encounter rate for copepods of typical perception radius of ~ {\eta}, where {\eta} is the dissipative scale of turbulence, can be increased by a factor up to ~ 100 compared to the one experienced by passively transported fluid tracers. Furthermore, we address the effect of introducing in the LC model a minimal waiting time between consecutive jumps. It is shown that any encounter-rate enhancement is lost if such time goes beyond the dissipative time-scale of turbulence, {\tau}_{\eta}. Because typically in the ocean {\eta} ~ 0.001m and {\tau}_{\eta} ~ 1s, this provides stringent constraints on the turbulent-driven enhancement of encounter-rate due to a purely mechanical induced escape reaction.Comment: 11 pages, 10 figure

Axially-homogeneous Rayleigh-Benard convection in a cylindrical cell

Previous numerical studies have shown that the "ultimate regime of thermal convection" can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers are eliminated by replacing the walls with periodic boundary conditions (homogeneous Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^{1/2} and turbulence intensity as Re ~ Ra^{1/2}, where the Rayleigh number Ra indicates the strength of the driving force. However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally confined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic boundary layers, we still find Nu ~ Re ~ Ra^{1/2}. The system supports exact solutions composed of modes of exponentially growing vertical velocity and temperature fields, with Ra as the critical parameter determining the properties of these modes. Counterintuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions which can dominate the dynamics and lead to very high and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased

Exponentially growing solutions in homogeneous Rayleigh-Benard convection

It is shown that homogeneous Rayleigh-Benard flow, i.e., Rayleigh-Benard turbulence with periodic boundary conditions in all directions and a volume forcing of the temperature field by a mean gradient, has a family of exact, exponentially growing, separable solutions of the full non-linear system of equations. These solutions are clearly manifest in numerical simulations above a computable critical value of the Rayleigh number. In our numerical simulations they are subject to secondary numerical noise and resolution dependent instabilities that limit their growth to produce statistically steady turbulent transport.Comment: 4 pages, 3 figures, to be published in Phys. Rev. E - rapid communication

Matched filters for coalescing binaries detection on massively parallel computers

We discuss some computational problems associated to matched filtering of experimental signals from gravitational wave interferometric detectors in a parallel-processing environment. We then specialize our discussion to the use of the APEmille and apeNEXT processors for this task. Finally, we accurately estimate the performance of an APEmille system on a computational load appropriate for the LIGO and VIRGO experiments, and extrapolate our results to apeNEXT.Comment: 19 pages, 6 figure

Lagrangian filtered density function for LES-based stochastic modelling of turbulent dispersed flows

The Eulerian-Lagrangian approach based on Large-Eddy Simulation (LES) is one of the most promising and viable numerical tools to study turbulent dispersed flows when the computational cost of Direct Numerical Simulation (DNS) becomes too expensive. The applicability of this approach is however limited if the effects of the Sub-Grid Scales (SGS) of the flow on particle dynamics are neglected. In this paper, we propose to take these effects into account by means of a Lagrangian stochastic SGS model for the equations of particle motion. The model extends to particle-laden flows the velocity-filtered density function method originally developed for reactive flows. The underlying filtered density function is simulated through a Lagrangian Monte Carlo procedure that solves for a set of Stochastic Differential Equations (SDEs) along individual particle trajectories. The resulting model is tested for the reference case of turbulent channel flow, using a hybrid algorithm in which the fluid velocity field is provided by LES and then used to advance the SDEs in time. The model consistency is assessed in the limit of particles with zero inertia, when "duplicate fields" are available from both the Eulerian LES and the Lagrangian tracking. Tests with inertial particles were performed to examine the capability of the model to capture particle preferential concentration and near-wall segregation. Upon comparison with DNS-based statistics, our results show improved accuracy and considerably reduced errors with respect to the case in which no SGS model is used in the equations of particle motion

Statistical properties of an ideal subgrid-scale correction for Lagrangian particle tracking in turbulent channel flow

One issue associated with the use of Large-Eddy Simulation (LES) to investigate the dispersion of small inertial particles in turbulent flows is the accuracy with which particle statistics and concentration can be reproduced. The motion of particles in LES fields may differ significantly from that observed in experiments or direct numerical simulation (DNS) because the force acting on the particles is not accurately estimated, due to the availability of the only filtered fluid velocity, and because errors accumulate in time leading to a progressive divergence of the trajectories. This may lead to different degrees of inaccuracy in the prediction of statistics and concentration. We identify herein an ideal subgrid correction of the a-priori LES fluid velocity seen by the particles in turbulent channel flow. This correction is computed by imposing that the trajectories of individual particles moving in filtered DNS fields exactly coincide with the particle trajectories in a DNS. In this way the errors introduced by filtering into the particle motion equations can be singled out and analyzed separately from those due to the progressive divergence of the trajectories. The subgrid correction term, and therefore the filtering error, is characterized in the present paper in terms of statistical moments. The effects of the particle inertia and of the filter type and width on the properties of the correction term are investigated.Comment: 15 pages,24 figures. Submitted to Journal of Physics: Conference Serie