A simulation of solar convection at supergranulation scale

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$\omega$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations.Comment: Contribution to the proceedings of the workshop entitled "THEMIS and the new frontiers of solar atmosphere dynamics" (March 2001), 6 pages, to appear in Nuovo Cimento

Waves attractors in rotating fluids: a paradigm for ill-posed Cauchy problems

In the limit of low viscosity, we show that the amplitude of the modes of oscillation of a rotating fluid, namely inertial modes, concentrate along an attractor formed by a periodic orbit of characteristics of the underlying hyperbolic Poincar\'e equation. The dynamics of characteristics is used to elaborate a scenario for the asymptotic behaviour of the eigenmodes and eigenspectrum in the physically relevant r\'egime of very low viscosities which are out of reach numerically. This problem offers a canonical ill-posed Cauchy problem which has applications in other fields.Comment: 4 pages, 5 fi

Area evolution, bulk viscosity and entropy principles for dynamical horizons

We derive from Einstein equation an evolution law for the area of a trapping or dynamical horizon. The solutions to this differential equation show a causal behavior. Moreover, in a viscous fluid analogy, the equation can be interpreted as an energy balance law, yielding to a positive bulk viscosity. These two features contrast with the event horizon case, where the non-causal evolution of the area and the negative bulk viscosity require teleological boundary conditions. This reflects the local character of trapping horizons as opposed to event horizons. Interpreting the area as the entropy, we propose to use an area/entropy evolution principle to select a unique dynamical horizon and time slicing in the Cauchy evolution of an initial marginally trapped surface.Comment: Some references added, 5 pages, Phys. Rev. D, in pres

Percolation properties of 3-D multiscale pore networks: how connectivity controls soil filtration processes

Quantifying the connectivity of pore networks is a key issue not only for modelling fluid flow and solute transport in porous media but also for assessing the ability of soil ecosystems to filter bacteria, viruses and any type of living microorganisms as well inert particles which pose a contamination risk. Straining is the main mechanical component of filtration processes: it is due to size effects, when a given soil retains a conveyed entity larger than the pores through which it is attempting to pass. We postulate that the range of sizes of entities which can be trapped inside soils has to be associated with the large range of scales involved in natural soil structures and that information on the pore size distribution has to be complemented by information on a critical filtration size (CFS) delimiting the transition between percolating and non percolating regimes in multiscale pore networks. We show that the mass fractal dimensions which are classically used in soil science to quantify scaling laws in observed pore size distributions can also be used to build 3-D multiscale models of pore networks exhibiting such a critical transition. We extend to the 3-D case a new theoretical approach recently developed to address the connectivity of 2-D fractal networks (Bird and Perrier, 2009). Theoretical arguments based on renormalisation functions provide insight into multi-scale connectivity and a first estimation of CFS. Numerical experiments on 3-D prefractal media confirm the qualitative theory. These results open the way towards a new methodology to estimate soil filtration efficiency from the construction of soil structural models to be calibrated on available multiscale data

Mode identification in rapidly rotating stars

Context: Recent calculations of pulsation modes in rapidly rotating polytropic models and models based on the Self-Consistent Field method have shown that the frequency spectrum of low degree pulsation modes can be described by an empirical formula similar to Tassoul's asymptotic formula, provided that the underlying rotation profile is not too differential. Aims: Given the simplicity of this asymptotic formula, we investigate whether it can provide a means by which to identify pulsation modes in rapidly rotating stars. Methods: We develop a new mode identification scheme which consists in scanning a multidimensional parameter space for the formula coefficients which yield the best-fitting asymptotic spectra. This mode identification scheme is then tested on artificial spectra based on the asymptotic formula, on random frequencies and on spectra based on full numerical eigenmode calculations for which the mode identification is known beforehand. We also investigate the effects of adding random frequencies to mimic the effects of chaotic modes which are also expected to show up in such stars. Results: In the absence of chaotic modes, it is possible to accurately find a correct mode identification for most of the observed frequencies provided these frequencies are sufficiently close to their asymptotic values. The addition of random frequencies can very quickly become problematic and hinder correct mode identification. Modifying the mode identification scheme to reject the worst fitting modes can bring some improvement but the results still remain poorer than in the case without chaotic modes

Balltracking: an highly efficient method for tracking flow fields

We present a method for tracking solar photospheric flows that is highly efficient, and demonstrate it using high resolution MDI continuum images. The method involves making a surface from the photospheric granulation data, and allowing many small floating tracers or balls to be moved around by the evolving granulation pattern. The results are tested against synthesised granulation with known flow fields and compared to the results produced by Local Correlation tracking (LCT). The results from this new method have similar accuracy to those produced by LCT. We also investigate the maximum spatial and temporal resolution of the velocity field that it is possible to extract, based on the statistical properties of the granulation data. We conclude that both methods produce results that are close to the maximum resolution possible from granulation data. The code runs very significantly faster than our similarly optimised LCT code, making real time applications on large data sets possible. The tracking method is not limited to photospheric flows, and will also work on any velocity field where there are visible moving features of known scale length

A simulation of solar convection of supergranulation scale

We present here numerical simulations of surface solar convection which cover a box of 30×30×3.2 Mm3 with a resolutionof 315×315×82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges inthe k-ω diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations

Identification of gravity waves in hydrodynamical simulations

The excitation of internal gravity waves by an entropy bubble oscillating in an isothermal atmosphere is investigated using direct two-dimensional numerical simulations. The oscillation field is measured by a projection of the simulated velocity field onto the anelastic solutions of the linear eigenvalue problem for the perturbations. This facilitates a quantitative study of both the spectrum and the amplitudes of excited g-modes.Comment: 12 pages, 11 figures, Appendices only available onlin

Differential rotation in fully convective stars

Under the assumption of thermal wind balance and effective entropy mixing in constant rotation surfaces, the isorotational contours of the solar convective zone may be reproduced with great fidelity. Even at this early stage of development, this helioseismology fit may be used to put a lower bound on the midlatitude {\em radial} solar entropy gradient, which in good accord with standard mixing length theory. In this paper, we generalize this solar calculation to fully convective stars (and potentially planets), retaining the assumptions of thermal wind balance and effective entropy mixing in isorotational surfaces. It is found that each isorotation contour is of the form $R^2 = A+B\Phi(r)$, where $R$ is the radius from the rotation axis, $\Phi(r)$ is the (assumed spherical) gravitational potential, and $A$ and $B$ are constant along the contour. This result is applied to simple models of fully convective stars. Both solar-like surface rotation profiles (angular velocity decreasing toward the poles) as well as "antisolar" profiles (angular velocity increasing toward the poles) are modeled; the latter bear some suggestive resemblance to numerical simulations. We also perform exploratory studies of zonal surface flows similar to those seen in Jupiter and Saturn. In addition to providing a practical framework for understanding the results of large scale numerical simulations, our findings may also prove useful in dynamical calculations for which a simple but viable model for the background rotation profile in a convecting fluid is needed. Finally, our work bears directly on an important goal of the CoRoT program: to elucidate the internal structure of rotating, convecting stars.Comment: 21 pages, 20 figures. Accepted for publication in MNRA