Direct Numerical Simulation of a separated channel flow with a smooth profile

A direct numerical simulation (DNS) of a channel flow with one curved surface was performed at moderate Reynolds number (Re_tau = 395 at the inlet). The adverse pressure gradient was obtained by a wall curvature through a mathematical mapping from physical coordinates to Cartesian ones. The code, using spectral spanwise and normal discretization, combines the advantage of a good accuracy with a fast integration procedure compared to standard numerical procedures for complex geometries. The turbulent flow slightly separates on the profile at the lower curved wall and is at the onset of separation at the opposite flat wall. The thin separation bubble is characterized with a reversal flow fraction. Intense vortices are generated near the separation line on the lower wall but also at the upper wall. Turbulent normal stresses and kinetic energy budget are investigated along the channel.Comment: 23 pages, submitted to Journal of Turbulenc

Extented ionized gas emission and kinematics of the compact group galaxies in HCG 16: Signatures of mergers

We report on kinematic observations of Ha emission line from four late-type galaxies of Hickson Compact Group 16 (H16a,b,c and d) obtained with a scanning Fabry-Perot interferometer and samplings of 16 km/s and 1". The velocity fields show kinematic peculiarities for three of the four galaxies: H16b, c and d. Misalignments between the kinematic and photometric axes of gas and stellar components (H16b,c,d), double gas systems (H16c) and severe warping of the kinematic major axis (H16b and c) were some of the peculiarities detected. We conclude that major merger events have taken place in at least two of the galaxies group. H16c and d, based on their significant kinematic peculiarities, their double nuclei and high infrared luminosities. Their Ha gas content is strongly spatially concentred - H16d contains a peculiar bar-like structure confined to the inner $\sim$ 1 h^-1 kpc region. These observations are in agreement with predictions of simulations, namely that the gas flows towards the galaxy nucleus during mergers, forms bars and fuel the central activity. Galaxy H16b, and Sb galaxy, also presents some of the kinematic evidences for past accretion events. Its gas content, however, is very spare, limiting our ability to find other kinematic merging indicators, if they are present. We find that isolated mergers, i.e., they show an anormorphous morphology and no signs of tidal tails. Tidal arms and tails formed during the mergers may have been stripped by the group potential (Barnes & Hernquist 1992) ar alternatively they may have never been formed. Our observations suggest that HCG 16 may be a young compact group in formation throught the merging of close-by objects in a dense environment.Comment: Accepted for publication in ApJ. 35 pages, 13 figures. tar file gzipped and uuencode

Supernova Remnants in the Magellanic Clouds III: An X-ray Atlas of LMC Supernova Remnants

We have used archival ROSAT data to present X-ray images of thirty-one supernova remnants (SNRs) in the Large Magellanic Cloud (LMC). We have classified these remnants according to their X-ray morphologies, into the categories of Shell-Type, Diffuse Face, Centrally Brightened, Point-Source Dominated, and Irregular. We suggest possible causes of the X-ray emission for each category, and for individual features of some of the SNRs.Comment: 27 pages, 6 figures (9 figure files). To appear in the Supplement Series of the Astrophysical Journal, August 1999 Vol. 123 #

Quasilinear theory of the 2D Euler equation

We develop a quasilinear theory of the 2D Euler equation and derive an integro-differential equation for the evolution of the coarse-grained vorticity. This equation respects all the invariance properties of the Euler equation and conserves angular momentum in a circular domain and linear impulse in a channel. We show under which hypothesis we can derive a H-theorem for the Fermi-Dirac entropy and make the connection with statistical theories of 2D turbulence.Comment: 4 page

Kang-Redner Anomaly in Cluster-Cluster Aggregation

The large time, small mass, asymptotic behavior of the average mass distribution \pb is studied in a $d$-dimensional system of diffusing aggregating particles for $1\leq d \leq 2$. By means of both a renormalization group computation as well as a direct re-summation of leading terms in the small reaction-rate expansion of the average mass distribution, it is shown that \pb \sim \frac{1}{t^d} (\frac{m^{1/d}}{\sqrt{t}})^{e_{KR}} for $m \ll t^{d/2}$, where $e_{KR}=\epsilon +O(\epsilon ^2)$ and $\epsilon =2-d$. In two dimensions, it is shown that \pb \sim \frac{\ln(m) \ln(t)}{t^2} for $m \ll t/ \ln(t)$. Numerical simulations in two dimensions supporting the analytical results are also presented.Comment: 11 pages, 6 figures, Revtex

Beyond scaling and locality in turbulence

An analytic perturbation theory is suggested in order to find finite-size corrections to the scaling power laws. In the frame of this theory it is shown that the first order finite-size correction to the scaling power laws has following form $S(r) \cong cr^{\alpha_0}[\ln(r/\eta)]^{\alpha_1}$, where $\eta$ is a finite-size scale (in particular for turbulence, it can be the Kolmogorov dissipation scale). Using data of laboratory experiments and numerical simulations it is shown shown that a degenerate case with $\alpha_0 =0$ can describe turbulence statistics in the near-dissipation range $r > \eta$, where the ordinary (power-law) scaling does not apply. For moderate Reynolds numbers the degenerate scaling range covers almost the entire range of scales of velocity structure functions (the log-corrections apply to finite Reynolds number). Interplay between local and non-local regimes has been considered as a possible hydrodynamic mechanism providing the basis for the degenerate scaling of structure functions and extended self-similarity. These results have been also expanded on passive scalar mixing in turbulence. Overlapping phenomenon between local and non-local regimes and a relation between position of maximum of the generalized energy input rate and the actual crossover scale between these regimes are briefly discussed.Comment: extended versio

Long-time discrete particle effects versus kinetic theory in the self-consistent single-wave model

The influence of the finite number N of particles coupled to a monochromatic wave in a collisionless plasma is investigated. For growth as well as damping of the wave, discrete particle numerical simulations show an N-dependent long time behavior resulting from the dynamics of individual particles. This behavior differs from the one due to the numerical errors incurred by Vlasov approaches. Trapping oscillations are crucial to long time dynamics, as the wave oscillations are controlled by the particle distribution inhomogeneities and the pulsating separatrix crossings drive the relaxation towards thermal equilibrium.Comment: 11 pages incl. 13 figs. Phys. Rev. E, in pres

Scaling laws and vortex profiles in 2D decaying turbulence

We use high resolution numerical simulations over several hundred of turnover times to study the influence of small scale dissipation onto vortex statistics in 2D decaying turbulence. A self-similar scaling regime is detected when the scaling laws are expressed in units of mean vorticity and integral scale, as predicted by Carnevale et al., and it is observed that viscous effects spoil this scaling regime. This scaling regime shows some trends toward that of the Kirchhoff model, for which a recent theory predicts a decay exponent $\xi=1$. In terms of scaled variables, the vortices have a similar profile close to a Fermi-Dirac distribution.Comment: 4 Latex pages and 4 figures. Submitted to Phys. Rev. Let

Relaxation equations for two-dimensional turbulent flows with a prior vorticity distribution

Using a Maximum Entropy Production Principle (MEPP), we derive a new type of relaxation equations for two-dimensional turbulent flows in the case where a prior vorticity distribution is prescribed instead of the Casimir constraints [Ellis, Haven, Turkington, Nonlin., 15, 239 (2002)]. The particular case of a Gaussian prior is specifically treated in connection to minimum enstrophy states and Fofonoff flows. These relaxation equations are compared with other relaxation equations proposed by Robert and Sommeria [Phys. Rev. Lett. 69, 2776 (1992)] and Chavanis [Physica D, 237, 1998 (2008)]. They can provide a small-scale parametrization of 2D turbulence or serve as numerical algorithms to compute maximum entropy states with appropriate constraints. We perform numerical simulations of these relaxation equations in order to illustrate geometry induced phase transitions in geophysical flows.Comment: 21 pages, 9 figure