Numerical simulations of compressible Rayleigh-Taylor turbulence in stratified fluids

We present results from numerical simulations of Rayleigh-Taylor turbulence, performed using a recently proposed lattice Boltzmann method able to describe consistently a thermal compressible flow subject to an external forcing. The method allowed us to study the system both in the nearly-Boussinesq and strongly compressible regimes. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.Comment: 15 pages, 11 figures. Proceedings of II Conference on Turbulent Mixing and Beyond (TMB-2009

Rain, power laws, and advection

Localized rain events have been found to follow power-law size and duration distributions over several decades, suggesting parallels between precipitation and seismic activity [O. Peters et al., PRL 88, 018701 (2002)]. Similar power laws are generated by treating rain as a passive tracer undergoing advection in a velocity field generated by a two-dimensional system of point vortices.Comment: 7 pages, 4 figure

Simple Viscous Flows: from Boundary Layers to the Renormalization Group

The seemingly simple problem of determining the drag on a body moving through a very viscous fluid has, for over 150 years, been a source of theoretical confusion, mathematical paradoxes, and experimental artifacts, primarily arising from the complex boundary layer structure of the flow near the body and at infinity. We review the extensive experimental and theoretical literature on this problem, with special emphasis on the logical relationship between different approaches. The survey begins with the developments of matched asymptotic expansions, and concludes with a discussion of perturbative renormalization group techniques, adapted from quantum field theory to differential equations. The renormalization group calculations lead to a new prediction for the drag coefficient, one which can both reproduce and surpass the results of matched asymptotics

On visualizing continuous turbulence scales

Turbulent flows are multi‐scale with vortices spanning a wide range of scales continuously. Due to such complexities, turbulence scales are particularly difficult to analyse and visualize. In this work, we present a novel and efficient optimization‐based method for continuous‐scale turbulence structure visualization with scale decomposition directly in the Kolmogorov energy spectrum. To achieve this, we first derive a new analytical objective function based on integration approximation. Using this new formulation, we can significantly improve the efficiency of the underlying optimization process and obtain the desired filter in the Kolmogorov energy spectrum for scale decomposition. More importantly, such a decomposition allows a ‘continuous‐scale visualization’ that enables us to efficiently explore the decomposed turbulence scales and further analyse the turbulence structures in a continuous manner. With our approach, we can present scale visualizations of direct numerical simulation data sets continuously over the scale domain for both isotropic and boundary layer turbulent flows. Compared with previous works on multi‐scale turbulence analysis and visualization, our method is highly flexible and efficient in generating scale decomposition and visualization results. The application of the proposed technique to both isotropic and boundary layer turbulence data sets verifies the capability of our technique to produce desirable scale visualization results

Elastic turbulence in curvilinear flows of polymer solutions

Following our first report (A. Groisman and V. Steinberg, \sl Nature $\bf 405$, 53 (2000)) we present an extended account of experimental observations of elasticity induced turbulence in three different systems: a swirling flow between two plates, a Couette-Taylor (CT) flow between two cylinders, and a flow in a curvilinear channel (Dean flow). All three set-ups had high ratio of width of the region available for flow to radius of curvature of the streamlines. The experiments were carried out with dilute solutions of high molecular weight polyacrylamide in concentrated sugar syrups. High polymer relaxation time and solution viscosity ensured prevalence of non-linear elastic effects over inertial non-linearity, and development of purely elastic instabilities at low Reynolds number (Re) in all three flows. Above the elastic instability threshold, flows in all three systems exhibit features of developed turbulence. Those include: (i)randomly fluctuating fluid motion excited in a broad range of spatial and temporal scales; (ii) significant increase in the rates of momentum and mass transfer (compared to those expected for a steady flow with a smooth velocity profile). Phenomenology, driving mechanisms, and parameter dependence of the elastic turbulence are compared with those of the conventional high Re hydrodynamic turbulence in Newtonian fluids.Comment: 23 pages, 26 figure

Starcounts Redivivus. IV. Density Laws Through Photometric Parallaxes

In an effort to more precisely define the spatial distribution of Galactic field stars, we present an analysis of the photometric parallaxes of 70,000 stars covering nearly 15 square degrees in seven Kapteyn Selected Areas. We address the affects of Malmquist Bias, subgiant/giant contamination, metallicity and binary stars upon the derived density laws. The affect of binary stars is the most significant. We find that while the disk-like populations of the Milky Way are easily constrained in a simultaneous analysis of all seven fields, no good simultaneous solution for the halo is found. We have applied halo density laws taken from other studies and find that the Besancon flattened power law halo model (c/a=0.6, r^-2.75) produces the best fit to our data. With this halo, the thick disk has a scale height of 750 pc with an 8.5% normalization to the old disk. The old disk scale height is 280-300 pc. Corrected for a binary fraction of 50%, these scale heights are 940 pc and 350-375 pc, respectively. Even with this model, there are systematic discrepancies between the observed and predicted density distributions. Our model produces density overpredictions in the inner Galaxy and density underpredictions in the outer Galaxy. A possible solution is modeling the stellar halo as a two-component system in which the halo has a flattened inner distribution and a roughly spherical, but substructured outer distribution. Further reconciliation could be provided by a flared thick disk, a structure consistent with a merger origin for that population. (Abridged)Comment: 66 pages, accepted to Astrophysical journal, some figures compresse

Dynamics and Excitation of Radio Galaxy Emission-Line Regions - I. PKS 2356-61

Results are presented from a programme of detailed longslit spectroscopic observations of the extended emission-line region (EELR) associated with the powerful radio galaxy PKS 2356-61. The observations have been used to construct spectroscopic datacubes, which yield detailed information on the spatial variations of emission-line ratios across the EELR, together with its kinematic structure. We present an extensive comparison between the data and results obtained from the MAPPINGS II shock ionization code, and show that the physical properties of the line-emitting gas, including its ionization, excitation, dynamics and overall energy budget, are entirely consistent with a scenario involving auto-ionizing shocks as the dominant ionization mechanism. This has the advantage of accounting for the observed EELR properties by means of a single physical process, thereby requiring less free parameters than the alternative scheme involving photoionization by radiation from the active nucleus. Finally, possible mechanisms of shock formation are considered in the context of the dynamics and origin of the gas, specifically scenarios involving infall or accretion of gas during an interaction between the host radio galaxy and a companion galaxy.Comment: 35 pages, LaTeX, uses aas2pp4.sty file, includes 9 PostScript figures. Two additional colour plates are available from the authors upon request. Accepted for publication in the Astrophysical Journa

Twenty-five years of multifractals in fully developed turbulence: a tribute to Giovanni Paladin

The paper {\it On the multifractal nature of fully developed turbulence and chaotic systems}, by R. Benzi {\it et al.} published in this journal in 1984 (vol {\bf 17}, page 3521) has been a starting point of many investigations on the different faces of selfsimilarity and intermittency in turbulent phenomena. Since then, the multifractal model has become a useful tool for the study of small scale turbulence, in particular for detailed predictions of different Eulerian and Lagrangian statistical properties. In the occasion of the 50-th birthday of our unforgettable friend and colleague Giovanni Paladin (1958-1996), we review here the basic concepts and some applications of the multifractal model for turbulence

The USNO-B Catalog

USNO-B is an all-sky catalog that presents positions, proper motions, magnitudes in various optical passbands, and star/galaxy estimators for 1,042,618,261 objects derived from 3,643,201,733 separate observations. The data were obtained from scans of 7,435 Schmidt plates taken for the various sky surveys during the last 50 years. USNO-B1.0 is believed to provide all-sky coverage, completeness down to V = 21, 0.2 arcsecond astrometric accuracy at J2000, 0.3 magnitude photometric accuracy in up to five colors, and 85% accuracy for distinguishing stars from non-stellar objects. A brief discussion of various issues is given here, but the actual data are available from http://www.nofs.navy.mil and other sites.Comment: Accepted by Astronomical Journa

Blood flow and coherent vortices in the normal and aneurysmatic aortas: a fluid dynamical approach to intra-luminal thrombus formation

Abdominal aortic aneurysms (AAAs) are frequently characterized by the development of an intra-luminal thrombus (ILT), which is known to have multiple biochemical and biomechanical implications. Development of the ILT is not well understood, and shear–stress-triggered activation of platelets could be the first step in its evolution. Vortical structures (VSs) in the flow affect platelet dynamics, which motivated the present study of a possible correlation between VS and ILT formation in AAAs. VSs educed by the λ2-method using computational fluid dynamics simulations of the backward-facing step problem, normal aorta, fusiform AAA and saccular AAA were investigated. Patient-specific luminal geometries were reconstructed from computed tomography scans, and Newtonian and Carreau–Yasuda models were used to capture salient rheological features of blood flow. Particularly in complex flow domains, results depended on the constitutive model. VSs developed all along the normal aorta, showing that a clear correlation between VSs and high wall shear stress (WSS) existed, and that VSs started to break up during late systole. In contrast, in the fusiform AAA, large VSs developed at sites of tortuous geometry and high WSS, occupying the entire lumen, and lasting over the entire cardiac cycle. Downward motion of VSs in the AAA was in the range of a few centimetres per cardiac cycle, and with a VS burst at that location, the release (from VSs) of shear-stress-activated platelets and their deposition to the wall was within the lower part of the diseased artery, i.e. where the thickest ILT layer is typically observed. In the saccular AAA, only one VS was found near the healthy portion of the aorta, while in the aneurysmatic bulge, no VSs occurred. We present a fluid-dynamics-motivated mechanism for platelet activation, convection and deposition in AAAs that has the potential of improving our current understanding of the pathophysiology of fluid-driven ILT growth