Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

Dependence of the large-scale vortex instability on latitude, stratification and domain size

In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions for the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain sizes, increasing stratification, and place the box at different latitudes. The effect of the increasing box size is to increase the sizes of the generated structures, so that the principal vortex always fills roughly half of the computational domain. The instability becomes stronger in the sense that the temperature anomaly and change in the radial velocity are observed to be enhanced. The model with the smallest box size is found to be stable against the instability, suggesting that a sufficient scale separation between the convective eddies and the scale of the domain is required for the instability to work. The instability can be seen upto the colatitude of 30 degrees, above which value the flow becomes dominated by other types of mean flows. The instability can also be seen in a model with larger stratification. Unlike the weakly stratified cases, the temperature anomaly caused by the vortex structures is seen to depend on depth.

The shear dynamo problem for small magnetic Reynolds numbers

We study large-scale kinematic dynamo action due to turbulence in the presence of a linear shear flow in the low-conductivity limit. Our treatment is non-perturbative in the shear strength and makes systematic use of both the shearing coordinate transformation and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed to have low magnetic Reynolds number (Re-m), but could have arbitrary fluid Reynolds number. The equation for the magnetic fluctuations is expanded perturbatively in the small quantity, Re-m. Our principal results are as follows: (i) the magnetic fluctuations are determined to the lowest order in Rem by explicit calculation of the resistive Green's function for the linear shear flow; (ii) the mean electromotive force is then calculated and an integro-differential equation is derived for the time evolution of the mean magnetic field. In this equation, velocity fluctuations contribute to two different kinds of terms, the 'C' and 'D' terms, respectively, in which first and second spatial derivatives of the mean magnetic field, respectively, appear inside the space-time integrals; (iii) the contribution of the D term is such that its contribution to the time evolution of the cross-shear components of the mean field does not depend on any other components except itself. Therefore, to the lowest order in Re-m, but to all orders in the shear strength, the D term cannot give rise to a shear-current-assisted dynamo effect; (iv) casting the integro-differential equation in Fourier space, we show that the normal modes of the theory are a set of shearing waves, labelled by their sheared wavevectors; (v) the integral kernels are expressed in terms of the velocity-spectrum tensor, which is the fundamental dynamical quantity that needs to be specified to complete the integro-differential equation description of the time evolution of the mean magnetic field; (vi) the C term couples different components of the mean magnetic field, so they can, in principle, give rise to a shear-current-type effect. We discuss the application to a slowly varying magnetic field, where it can be shown that forced non-helical velocity dynamics at low fluid Reynolds number does not result in a shear-current-assisted dynamo effect.

The modulation of radiation in an electron-positron plasma

The modulational instability of a large-amplitude, linearly polarized electromagnetic wave propagating in an electron-positron plasma is considered, including the combined effect of relativistic mass variation of the plasma particles, harmonic generation, and the non-resonant, finite-frequency electrostatic density perturbations, all caused by the large-amplitude radiation field. The radiation from many strong sources, such as AGN and pulsars, has been observed to vary over a host of time-scales. It is possible that the extremely rapid variations in the non-thermal continuum of AGN, as well as in the non-thermal radio radiation from pulsars, can be accounted for by the modulational instabilities to which radiation may be subjected during its propagation out of the emission region.

A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.

Proposed visible emission line space solar coronagraph

The outer atmosphere of the sun called the corona has been observed during total solar eclipse for short periods (typically <6 min), from as early as the eighteenth century. In the recent past, space-based instruments have permitted us to study the corona uninterruptedly. In spite of these developments, the dynamic corona and its high temperature (1-2 million K) are yet to be Ally understood. It is conjectured that their dynamic nature and associated energetic events are possible reasons behind the high temperature. In order to study these in detail, a visible emission line space solar coronagraph is being proposed as a payload under the small-satellite programme of the Indian Space Research Organisation. The satellite is named as Aditya-1 and the scientific objectives of this payload are to study: (i) the existence of intensity oscillations for the study of wave-driven coronal heating; (ii) the dynamics and formation of coronal loops and temperature structure of the coronal features; (iii) the origin, cause and acceleration of coronal mass ejections (CMEs) and other solar active features, and (iv) coronal magnetic field topology and three-dimensional structures of CMEs using polarization information. The uniqueness of this payload compared to previously flown space instruments is as follows: (a) observations in the visible wavelength closer to the disk (down to 1.05 solar radii); (b) high time cadence capability (better than two-images per second), and (c) simultaneous observations of at least two spectral windows all the time and three spectral windows for short durations.

Gain calibration of CCD systems at VBO

The system gain of two CCD systems in regular use at the Vainu Bappu Observatory, Kavalur, is determined at a few gain settings. The procedure used for the determination of system gain and base-level noise is described in detail. The Photometrics CCD system at the 1-m reflector uses a Thomson-CSF TH 7882 CDA chip coated for increased ultraviolet sensitivity. The gain is programme-selected through the parameter 'cgain' varying between 0 and 4095 in steps of 1. The inverse system gain for this system varies almost linearly from 27.7 electrons DN-1 at cgain = 0 to 1.5 electrons DN-1 at cgain = 500. The readout noise is less than or similar 11 electrons at cgain = 66. The Astromed CCD system at 2.3-m Vainu Bappu Telescope uses a GEC P8603 chip which is also coated for enhanced ultraviolet sensitivity. The amplifier gain is selected in discrete steps using switches in the controller. The inverse system gain is 4.15 electrons DN-1 at the gain setting of 9.2, and the readout noise approximately 8 electrons.

Flat rotation curves: A result of the helical inverse cascade in turbulent media

We have modeled the rotation curves of 21 galaxies observed by Amram et al. (1992), by combining the effects of rigid rotation, gravity, and turbulence. The main motivation behind such modeling is to study the formation of coherent structures in turbulent media and explore its role in the large-scale structures of the universe. The values of the parameters such as mass, turbulent velocity, and angular velocity derived from the rotation curve fits are in good agreement with those derived from the prevalent models.

Perturbative growth of cosmological clustering. I: Formalism

Here we rederive the hierarchy of equations for the evolution of distribution functions of various orders using a convenient parameterization. We use this to obtain equations for two- and three-point correlation functions in powers of a small parameter, viz., the initial density contrast. The correspondence of the lowest order solutions of these equations to the results from the linear theory of density perturbations is shown for an OMEGA = 1 universe. These equations are then used to calculate, to the lowest order, the induced three-point correlation function that arises from Gaussian initial conditions in an OMEGA = 1 universe. We obtain an expression which explicitly exhibits the spatial structure of the induced three-point correlation function. It is seen that the spatial structure of this quantity is independent of the value of OMEGA. We also calculate the triplet momentum. We find that the induced three-point correlation function does not have the ''hierarchical'' form often assumed. We discuss possibilities of using the induced three-point correlation to interpret observational data. The formalism developed here can also be used to test a validity of different schemes to close the

Propagating intensity disturbances in polar corona as seen from AIA/SDO

Context. Polar corona is often explored to find the energy source for the acceleration of the fast solar wind. Earlier observations show omni-presence of quasi-periodic disturbances, traveling outward, which is believed to be caused by the ubiquitous presence of outward propagating waves. These waves, mostly of compressional type, might provide the additional momentum and heat required for the fast solar wind acceleration. It has been conjectured that these disturbances are not due to waves but high speed plasma outflows, which are difficult to distinguish using the current available techniques. Aims. With the unprecedented high spatial and temporal resolution of AIA/SDO, we search for these quasi-periodic disturbances in both plume and interplume regions of the polar corona. We investigate their nature of propagation and search for a plausible interpretation. We also aim to study their multi-thermal nature by using three different coronal passbands of AIA. Methods. We chose several clean plume and interplume structures and studied the time evolution of specific channels by making artificial slits along them. Taking the average across the slits, space-time maps are constructed and then filtration techniques are applied to amplify the low-amplitude oscillations. To suppress the effect of fainter jets, we chose wider slits than usual. Results. In almost all the locations chosen, in both plume and interplume regions we find the presence of propagating quasi-periodic disturbances, of periodicities ranging from 10-30 min. These are clearly seen in two channels and in a few cases out to very large distances (approximate to 250 `') off-limb, almost to the edge of the AIA field of view. The propagation speeds are in the range of 100-170 km s(-1). The average speeds are different for different passbands and higher in interplume regions. Conclusions. Propagating disturbances are observed, even after removing the effects of jets and are insensitive to changes in slit width. This indicates that a coherent mechanism is involved. In addition, the observed propagation speed varies between the different passpands, implying that these quasi-periodic intensity disturbances are possibly due to magneto-acoustic waves. The propagation speeds in interplume region are higher than in the plume region.

The gasification of wood-char spheres in CO2---N2 mixtures: analysis and experiments

The gasification of charcoal spheres in an atmosphere of carbon-dioxide-nitrogen mixture involving diffusion and reactions in the pores is modelled and the results are compared with experiments of Standish and Tanjung and those performed in the laboratory on wood-char spheres to determine the effects of diameter, density, gas composition and flow. The results indicate that the conversion time, t(c) approximately d1.03 for large particles (> 5 mm), departing substantially from the t(c) approximately d2 law valid for diffusion limited conditions. The computational studies indicate that the kinetic limit for the particle is below 100 mum. The conversion time varies inversely as the initial char density as expected in the model. Predictions from the model show that there is no significant change in conversion time up to 60% N2 consistent with the CO2-N2 experiments. The variation of diameter and density with time are predicted. The peculiar dependence of conversion time on flow velocity in the experiments is sought to be explained by opposing free and forced convection heat transfer and the attempt is only partly successful. The studies also indicate that the dependence on the CO concentration with low CO2 is significant, indicating the need for multistep reaction mechanism against the generally accepted single-step reaction.