Numerical studies on quasilinear and linear elliptic equations

We explore here the acceleration of convergence of iterative methods for the solution of a class of quasilinear and linear algebraic equations. The specific systems are the finite difference form of the Navier-Stokes equations and the energy equation for recirculating flows. The acceleration procedures considered are: the successive over relaxation scheme; several implicit methods; and a second-order procedure. A new implicit method—the alternating direction line iterative method—is proposed in this paper. The method combines the advantages of the line successive over relaxation and alternating direction implicit methods. The various methods are tested for their computational economy and accuracy on a typical recirculating flow situation. The numerical experiments show that the alternating direction line iterative method is the most economical method of solving the Navier-Stokes equations for all Reynolds numbers in the laminar regime. The usual ADI method is shown to be not so attractive for large Reynolds numbers because of the loss of diagonal dominance. This loss can however be restored by a suitable choice of the relaxation parameter, but at the cost of accuracy. The accuracy of the new procedure is comparable to that of the well-tested successive overrelaxation method and to the available results in the literature. The second-order procedure turns out to be the most efficient method for the solution of the linear energy equation.

An HI study of three long-tailed irregular galaxies in the cluster Abell 1367

We present the results on the distribution and kinematics of HI gas with higher sensitivity and in one case of higher spectral resolution as well than reported earlier, of three irregular galaxies CGCG 097073, 097079 and 097087 (UGC 06697) in the cluster Abell 1367. These galaxies are known to exhibit long (50 - 75 kpc) tails of radio continuum and optical emission lines (H alpha) pointing away from the cluster centre and arcs of starformation on the opposite sides of the tails, These features as well as the HI properties, with two of the galaxies (CGCG 097073 and 097079) exhibiting sharper gradients in HI intensity on the side of the tails, are consistent with the HI gas being affected by the ram pressure of the intracluster medium. However the HI emission in all the three galaxies extends to much smaller distances than the radio-continuum and H alpha tails, and are possibly still bound to the parent galaxies. Approximately 20 - 30 per cent of the HI mass is seen to accumulate on the downstream side due to the effects of ram pressure.

The Formation, Numbers And Radio Output Of Giant Radio Galaxies

The numbers and mean radio luminosities of giant radio galaxies (GRGs) have been calculated for redshifts up to z = 0.6, assuming a sensitivity limit of 1 Jy at 1 GHz for the observations. The estimates are obtained with a model for the beam propagation, first through the hot gaseaous halo around the parent galaxy, and thereafter, through the even hotter but less dense intergalactic medium. The model is able to accurately reproduce the observed numbers and mean radio luminosities of GRGs at redshifts of less than 0.1, and it predicts that a somewhat larger number of GRGs should be found at redshifts of greater than 0.1.

PIV studies of large scale structures in the near field of small aspect ratio elliptic jets

The near flow field of small aspect ratio elliptic turbulent free jets (issuing from nozzle and orifice) was experimentally studied using a 2D PIV. Two point velocity correlations in these jets revealed the extent and orientation of the large scale structures in the major and minor planes. The spatial filtering of the instantaneous velocity field using Gaussian convolution kernel shows that while a single large vortex ring circumscribing the jet seems to be present at the exit of nozzle, the orifice jet exhibited a number of smaller vortex ring pairs close to jet exit. The smaller length scale observed in the case of the orifice jet is representative of the smaller azimuthal vortex rings that generate axial vortex field as they are convected. This results in the axis-switching in the case of orifice jet and may have a mechanism different from the self induction process as observed in the case of contoured nozzle jet flow.

Hot gaseous coronae of early-type galaxies and their radio luminosity function

Recent X-ray observations have revealed that early-type galaxies (which usually produce extended double radio sources) generally have hot gaseous haloes extending up to approx102kpc1,2. Moreover, much of the cosmic X-ray background radiation is probably due to a hotter, but extremely tenuous, intergalactic medium (IGM)3. We have presented4–7 an analytical model for the propagation of relativistic beams from galactic nuclei, in which the beams' crossing of the pressure-matched interface between the IGM and the gaseous halo, plays an important role. The hotspots at the ends of the beams fade quickly when their advance becomes subsonic with respect to the IGM. This model has successfully predicted (for typical double radio sources) the observed8 current mean linear-size (approx2Dsime350 kpc)4,5, the observed8–11 decrease in linear-size with cosmological redshift4–6 and the slope of the linear-size versus radio luminosity10,12–14 relation6. We have also been able to predict the redshift-dependence of observed numbers and radio luminosities of giant radio galaxies7,15. Here, we extend this model to include the propagation of somewhat weaker beams. We show that the observed flattening of the local radio luminosity function (LRLF)16–20 for radio luminosity Papproximately 1024 W Hz-1 at 1 GHz can be explained without invoking ad hoc a corresponding break in the beam power function Phi(Lb), because the heads of the beams with Lb < 1025 W Hz-1 are decelerated to sonic velocity within the halo itself, which leads to a rapid decay of radio luminosity and a reduced contribution of these intrinsically weaker sources to the observed LRLF.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

H I deficiency in cluster spiral galaxies - Dependence on galaxy size

From the available H I data on spiral galaxies in three rich Abell clusters and the Virgo Cluster, it is shown that galaxies with medium to large optical sizes tend to be more severely deficient in atomic hydrogen than the small galaxies. This is so both in terms of the fractional number of galaxies that are deficient and the amount of gas lost by a galaxy. The fraction of H I-deficient galaxies increases with size over most of the size range, saturating or dropping only for the largest galaxies. A comparative study is made of various currently accepted gas removal mechanisms, namely those which are a result of galaxy-intracluster medium interactions, e.g., ram pressure stripping, as well as those due to galaxy-galaxy interactions, i.e., collisions and tidal interactions. It is shown that, with the exception of tidal interactions, all of these mechanisms would produce a size dependence in H I deficiency that is the opposite of that observed. That is, the gas in the largest galaxies would be the least affected by these mechanisms. However, if there is significant mass segregation, these processes may give the trends observed in the size dependence of H I deficiency.

Some nonstandard error analysis of discontinuous Galerkin methods for elliptic problems

An a priori error analysis of discontinuous Galerkin methods for a general elliptic problem is derived under a mild elliptic regularity assumption on the solution. This is accomplished by using some techniques from a posteriori error analysis. The model problem is assumed to satisfy a GAyenrding type inequality. Optimal order L (2) norm a priori error estimates are derived for an adjoint consistent interior penalty method.

Elliptic flow of thermal photons in heavy-ion collisions at Relativistic Heavy Ion Collider and Large Hadron Collider

We calculate the thermal photon transverse momentum spectra and elliptic flow in $\sqrt{s_{NN}} = 200$ GeV Au+Au collisions at RHIC and in $\sqrt{s_{NN}} = 2.76$ TeV Pb+Pb collisions at the LHC, using an ideal-hydrodynamical framework which is constrained by the measured hadron spectra at RHIC and LHC. The sensitivity of the results to the QCD-matter equation of state and to the photon emission rates is studied, and the photon $v_2$ is discussed in the light of the photonic $p_T$ spectrum measured by the PHENIX Collaboration. In particular, we make a prediction for the thermal photon $p_T$ spectra and elliptic flow for the current LHC Pb+Pb collisions.

Free convection about cylinders of elliptic cross section embedded in a porous medium

Numerical solutions are presented for the free convection boundary layers over cylinders of elliptic cross section embedded in a fluid-saturated porous medium. The transformed conservation equations of the nonsimilar boundary layers are solved numerically by an efficient finite-difference method. The theory was applied to a number of cylinders and the results compared very well with published analytical solutions. The results are of use in the design of underground electrical cables, power plant steam, and water distribution lines, among others.

Seyfert Galaxies: Nuclear Radio Structure And Unification

A radio study of a carefully selected sample of 20 Seyfert galaxies that are matched in orientation-independent parameters, which are measures of intrinsic active galactic nucleus power and host galaxy properties, is presented to test the predictions of the unified scheme hypothesis. Our sample sources have core flux densities greater than 8 mJy at 5 GHz on arcsec scales due to the feasibility requirements. These simultaneous parsec-scale and kiloparsec-scale radio observations reveal (1) that Seyfert 1 and Seyfert 2 galaxies have an equal tendency to show compact radio structures on milliarcsecond scales, (2) the distributions of parsec-scale and kiloparsec-scale radio luminosities are similar for both Seyfert 1 and Seyfert 2 galaxies, (3) there is no evidence for relativistic beaming in Seyfert galaxies, (4) similar distributions of source spectral indices in spite of the fact that Seyferts show nuclear radio flux density variations, and (5) the distributions of the projected linear size for Seyfert 1 and Seyfert 2 galaxies are not significantly different as would be expected in the unified scheme. The latter could be mainly due to a relatively large spread in the intrinsic sizes. We also find that a starburst alone cannot power these radio sources. Finally, an analysis of the kiloparsec-scale radio properties of the CfA Seyfert galaxy sample shows results consistent with the predictions of the unified scheme.