On the evolution of flow topology in turbulent Rayleigh-Bénard convection

Copyright 2016 AIP Publishing. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing.Small-scale dynamics is the spirit of turbulence physics. It implicates many attributes of flow topology evolution, coherent structures, hairpin vorticity dynamics, and mechanism of the kinetic energy cascade. In this work, several dynamical aspects of the small-scale motions have been numerically studied in a framework of Rayleigh-Benard convection (RBC). To do so, direct numerical simulations have been carried out at two Rayleigh numbers Ra = 10(8) and 10(10), inside an air-filled rectangular cell of aspect ratio unity and pi span-wise open-ended distance. As a main feature, the average rate of the invariants of the velocity gradient tensor (Q(G), R-G) has displayed the so-calledPeer ReviewedPostprint (author's final draft

LISA observations of supermassive black holes: parameter estimation using full post-Newtonian inspiral waveforms

We study parameter estimation of supermassive black hole binary systems in the final stage of inspiral using the full post-Newtonian gravitational waveforms. We restrict our analysis to systems in circular orbit with negligible spins, in the mass range 10^8\Ms-10^5\Ms, and compare the results with those arising from the commonly used restricted post-Newtonian approximation. The conclusions of this work are particularly important with regard to the astrophysical reach of future LISA measurements. Our analysis clearly shows that modeling the inspiral with the full post-Newtonian waveform, not only extends the reach to higher mass systems, but also improves in general the parameter estimation. In particular, there are remarkable improvements in angular resolution and distance measurement for systems with a total mass higher than 5\times10^6\Ms, as well as a large improvement in the mass determination.Comment: Final version. Accepted for publication in Phys. Rev.

Studying stellar binary systems with the Laser Interferometer Space Antenna using Delayed Rejection Markov chain Monte Carlo methods

Bayesian analysis of LISA data sets based on Markov chain Monte Carlo methods has been shown to be a challenging problem, in part due to the complicated structure of the likelihood function consisting of several isolated local maxima that dramatically reduces the efficiency of the sampling techniques. Here we introduce a new fully Markovian algorithm, a Delayed Rejection Metropolis-Hastings Markov chain Monte Carlo method, to efficiently explore these kind of structures and we demonstrate its performance on selected LISA data sets containing a known number of stellar-mass binary signals embedded in Gaussian stationary noise.Comment: 12 pages, 4 figures, accepted in CQG (GWDAW-13 proceedings

New subgrid-scale models for large-eddy simulation of Rayleigh-Bénard convection

Published under licence in Journal of Physics: Conference Series by IOP Publishing Ltd. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.At the crossroad between flow topology analysis and the theory of turbulence, a new eddy-viscosity model for Large-eddy simulation has been recently proposed by Trias et al.[PoF, 27, 065103 (2015)]. The S3PQR-model has the proper cubic near-wall behaviour and no intrinsic limitations for statistically inhomogeneous flows. In this work, the new model has been tested for an air turbulent Rayleigh-Benard convection in a rectangular cell of aspect ratio unity and n span-wise open-ended distance. To do so, direct numerical simulation has been carried out at two Rayleigh numbers Ra = 108 and 1010, to assess the model performance and investigate a priori the effect of the turbulent Prandtl number. Using an approximate formula based on the Taylor series expansion, the turbulent Prandtl number has been calculated and revealed a constant and Ra-independent value across the bulk region equals to 0.55. It is found that the turbulent components of eddy-viscosity and eddy-diffusivity are positively prevalent to maintain a turbulent wind essentially driven by the mean buoyant force at the sidewalls. On the other hand, the new eddy-viscosity model is preliminary tested for the case of Ra = 108 and showed overestimation of heat flux within the boundary layer but fairly good prediction of turbulent kinetics at this moderate turbulent flow.Peer ReviewedPostprint (published version