We simulate three-dimensional, horizontally periodic Rayleigh-B\'enard convection between free-slip horizontal plates, rotating about a distant horizontal axis. When both the temperature difference between the plates and the rotation rate are sufficiently large, a strong horizontal wind is generated that is perpendicular to both the rotation vector and the gravity vector. The wind is turbulent, large-scale, and vertically sheared. Horizontal anisotropy, engendered here by rotation, appears necessary for such wind generation. Most of the kinetic energy of the flow resides in the wind, and the vertical turbulent heat flux is much lower on average than when there is no wind.

The passage from kinetic theory to fluid dynamics as discussed by Hilbert has perplexed Uhlenbeck. Herein, I try to smooth over the discord.

We report on direct numerical simulations of two-dimensional, horizontally periodic Rayleigh-B\'enard convection, focusing on its ability to drive large-scale horizontal flow that is vertically sheared. For the Prandtl numbers ($Pr$) between 1 and 10 simulated here, this large-scale shear can be induced by raising the Rayleigh number ($Ra$) sufficiently, and we explore the resulting convection for $Ra$ up to $10^{10}$. When present in our simulations, the sheared mean flow accounts for a large fraction of the total kinetic energy, and this fraction tends towards unity as $Ra\to\infty$. The shear helps disperse convective structures, and it reduces vertical heat flux; in parameter regimes where one state with large-scale shear and one without are both stable, the Nusselt number of the state with shear is smaller and grows more slowly with $Ra$. When the large-scale shear is present with $Pr\lesssim2$, the convection undergoes strong global oscillations on long timescales, and heat transport occurs in bursts. Nusselt numbers, time-averaged over these bursts, vary non-monotonically with $Ra$ for $Pr=1$. When the shear is present with $Pr\gtrsim3$, the flow does not burst, and convective heat transport is sustained at all times. Nusselt numbers then grow roughly as powers of $Ra$, but the growth rates are slower than any previously reported for Rayleigh-B\'enard convection without large-scale shear. We find the Nusselt numbers grow proportionally to $Ra^{0.077}$ when $Pr=3$ and to $Ra^{0.19}$ when $Pr=10$. Analogies with tokamak plasmas are described.

Recent developments in observational cosmology have led to attempts to make modifications on both sides of the Einstein equation to explain some of the puzzling new findings. What follows is an examination of the source of gravity that we usually find on the right hand side of Einstein's equation. The outcome is a modified version of the stress-energy tensor that is the source of the gravitational field. The derivation is based on the kinetic theory of a gas of identical particles with no internal structure. The presentation here is in two parts. In Part I, I describe the stress tensor that Xinzhong Chen and I have proposed for the matter tensor for a nonrelativistic gas with input from Hongling Rao and Jean-Luc Thiffeault. Our derivation of the equations of fluid dynamics is based on kinetic theory without recourse to the standard Chapman-Enskog approximation. In Part II, I present the analogous derivation of our form for the stress-energy tensor in the relativistic case. Then I exhibit its application to the usual isotropic cosmological model. The result of that, in addition to the Friedmann solution, is a second solution that arises from terms discarded in the usual Chapman-Enskog approximation. The new solution is a temporal analogue of a spatial shock wave. Just as the usual shock waves make transitions in properties within a mean free path, the new solution can change its properties appreciably in a mean flight time. Whereas the Friedmann solution is not dissipative, the new solution produces entropy at a rate that may be of cosmological interest. For the calculation of cosmic entropy production I use a formula derived in the ultrarelativistic limit in which particle masses are negligible. Independently of the cosmological aspects, the fluid dynamical equations that we derive are causal, even for the heat equation (or Fourier equation).

We derive causal relativistic fluid dynamical equations from the relaxation model of kinetic theory as in a procedure previously applied in the case of non-relativistic rarefied gases. By treating space and time on an equal footing and avoiding the iterative steps of the conventional Chapman-Enskog --- CE---method, we are able to derive causal equations in the first order of the expansion in terms of the mean flight time of the particles. This is in contrast to what is found using the CE approach. We illustrate the general results with the example of a gas of identical ultrarelativistic particles such as photons under the assumptions of homogeneity and isotropy. When we couple the fluid dynamical equations to Einstein's equation we find, in addition to the geometry-driven expanding solution of the FRW model, a second, matter-driven nonequilibrium solution to the equations. In only the second solution, entropy is produced at a significant rate.

We consider the effect of radial fluid injection and suction on Taylor-Couette flow. Injection at the outer cylinder and suction at the inner cylinder generally results in a linearly unstable steady spiralling flow, even for cylindrical shears that are linearly stable in the absence of a radial flux. We study nonlinear aspects of the unstable motions with the energy stability method. Our results, though specialized, may have implications for drag reduction by suction, accretion in astrophysical disks, and perhaps even in the flow in the earth's polar vortex.

We introduce on/off intermittency into a mean field dynamo model by imposing stochastic fluctuations in either the alpha effect or through the inclusion of a fluctuating electromotive force. Sufficiently strong small scale fluctuations with time scales of the order of 0.3-3 years can produce long term variations in the system on time scales of the order of hundreds of years. However, global suppression of magnetic activity in both hemispheres at once was not observed. The variation of the magnetic field does not resemble that of the sunspot number, but is more reminiscent of the 10 Be record. The interpretation of our results focuses attention on the connection between the level of magnetic activity and the sunspot number, an issue that must be elucidated if long term solar effects are to be well understood.

The description of a stellar system as a continuous fluid represents a convenient first approximation to stellar dynamics, and its derivation from the kinetic theory is standard. The challenge lies in providing adequate closure approximations for the higher-order moments of the phase-space density function that appear in the fluid dynamical equations. Such closure approximations may be found using representations of the phase-space density as embodied in the kinetic theory. In the classic approach of Chapman and Enskog, one is led to the Navier-Stokes equations, which are known to be inaccurate when the mean free paths of particles are long, as they are in many stellar systems. To improve on the fluid description, we derive here a modified closure relation using a Fokker-Planck collision operator. To illustrate the nature of our approximation, we apply it to the study of gravitational instability. The instability proceeds in a qualitative manner as given by the Navier-Stokes equations but, in our description, the damped modes are considerably closer to marginality, especially at small scales.

We analyze examples of projected two dimensional fractals for the case where the objects studied have a Schechter luminosity function. Our results indicate that two dimensional catalogues of sources could be of value in detecting the effects of luminosity functions as well as of large scale structure.

To investigate the mechanism of mixing in oscillatory doubly diffusive (ODD) convection, we truncate the horizontal modal expansion of the Boussinesq equations to obtain a simplified model of the process. In the astrophysically interesting case with low Prandtl number, large-scale shears are generated as in ordinary thermal convection. The interplay between the shear and the oscillatory convection produces intermittent overturning of the fluid with significant mixing. By contrast, in the parameter regime appropriate to sea water, large-scale flows are not generated by the convection. However, if such flows are imposed externally, intermittent overturning with enhanced mixing is observed.

We use the one-dimensional steady version of the equations derived in paper I to compute the structure of shock waves. The agreement with experiment is good, especially when we retain the experimental value of the Prandtl number adopted in II.

We describe an asymptotic procedure for deriving continuum equations from the kinetic theory of a simple gas. As in the works of Hilbert, of Chapman and of Enskog, we expand in the mean flight time of the constituent particles of the gas, but we do not adopt the Chapman-Enskog device of simplifying the formulae at each order by using results from previous orders. In this way, we are able to derive a new set of fluid dynamical equations from kinetic theory, as we illustrate here for the relaxation model for monatomic gases. We obtain a stress tensor that contains a dynamical pressure term (or bulk viscosity) that is process-dependent and our heat current depends on the gradients of both temperature and density. On account of these features, the equations apply to a greater range of Knudsen number (the ratio of mean free path to macroscopic scale) than do the Navier-Stokes equations, as we see in the accompanying paper. In the limit of vanishing Knudsen number, our equations reduce to the usual Navier-Stokes equations with no bulk viscosity.

The equations of fluid dynamics developed in paper I are applied to the study of the propagation of ultrasound waves. There is good agreement between the predicted propagation speed and experimental results for a wide range of Knudsen numbers.

Viscous resistance to changes in the volume of a gas arises when different degrees of freedom have different relaxation times. Collisions tend to oppose the resulting departures from equilibrium and, in so doing, generate entropy. Even for a classical gas of hard spheres, when the mean free paths or mean flight times of constituent particles are long, we find a nonvanishing bulk viscosity. Here we apply a method recently used to uncover this result for a classical rarefied gas to radiative transfer theory and derive an expression for the radiative stress tensor for a gray medium with absorption and Thomson scattering. We determine the transport coefficients through the calculation of the comoving entropy generation. When scattering dominates absorption, the bulk viscosity becomes much larger than either the shear viscosity or the thermal conductivity.

We use the transfer equation in relativistic form to develop an expansion of the one-photon distribution for a medium with constant photon mean free path, $\epsilon$. The resulting radiative viscosity tensor may be expressed either as a simple integral operator acting on the Thomas stress tensor or as the solution of an inhomogenous, linear partial differential equation. The expression obtained for the radiative viscosity tensor applies for media with long, as well as short, photon mean free paths. We also develop results applicable for relatively smooth flows by using the form of the Thomas stress tensor with generalized transport coefficients derived by the application of a suitable operator to the bare Thomas coefficients.

In the problem of planetary formation one seeks a mechanism to gather small solid particles together into larger accumulations of solid matter. Here we describe a scenario in which turbulence mediates this process by aggregating particles into anticyclonic regions. If, as our simulations suggest, anticyclonic vortices form as long-lived coherent structures, the process becomes more powerful because such vortices trap particles effectively. Even if the turbulence is decaying, following the upheaval that formed the disk, there is enough time to make the dust distribution quite lumpy.

We study the longitudinal instabilities of two interpenetrating fluids interacting only through gravity. When one of the constituents is of relatively low density, it is possible to have a band of unstable wave numbers well separated from those involved in the usual Jeans instability. If the initial streaming is large enough, and there is no linear instability, the indefinite sign of the free energy has the possible consequence of explosive interactions between positive and negative energy modes in the nonlinear regime. The effect of dissipation on the negative energy modes is also examined.

Beginning from a relatively simple set of dynamical equations for a fluid permeated by a radiative field strong enough to produce significant forces, we find the structure of plane-parallel equilibria and study their stability to small acoustic disturbances. In doing this, we neglect viscous effects and complications of nongreyness. We find that acoutic instabilities occur over a wide range of conditions below the Eddington limit. This result is in line with findings reported twenty years ago but it contradicts some more recent reports of the absence of instabilities. We briefly attempt to identify the causes of the discrepancies and then close with a discussion of the possible astrophysical interest of such instabilities.

If we treat the galaxies in published redshift catalogues as point sets, we may determine the generalized dimensions of such sets by standard means, outlined here. For galaxy separations up to about 5 Mpc, we find the dimensions of the galaxy set to be about 1.2, with not a strong indication of multifractality. For larger scales, out to about 30 Mpc, there is also good scaling with a dimension of about 1.8. For even larger scales, the data seem too sparse to be conclusive, but we find that the dimension is climbing as the scales increase. We report simulations that suggest a rationalization of such measurements, namely that in the intermediate range the scaling behavior is dominated by flat structures (pancakes) and that the results on the smallest scales are a reflection of the formation of density singularities.

We study the weakly nonlinear evolution of acoustic instability of a plane- parallel polytrope with thermal dissipation in the form of Newton's law of cooling. The most unstable horizontal wavenumbers form a band around zero and this permits the development of a nonlinear pattern theory leading to a complex Ginzburg-Landau equation (CGLE). Numerical solutions for a subcritical, quintic CGLE produce vertically oscillating, localized structures that resemble the oscillons observed in recent experiments of vibrated granular material.

Rotating, turbulent cosmic fluids are generally pervaded by coherent structures such as vortices and magnetic flux tubes. The formation of such structures is a robust property of rotating turbulence as has been confirmed in computer simulations and laboratory experiments. We defend here the notion that accretion disks share this feature of rotating cosmic bodies. In particular, we show that the intense shears of Keplerian flows do not inhibit the formation of vortices. Given suitable initial disturbances and high enough Reynolds numbers, long-lived vortices form in Keplerian shear flows and analogous magnetic structures form in magnetized disks. The formation of the structures reported here should have significant consequences for the transport properties of disks and for the observed properties of hot disks.

Calculations of the stopping power of a medium lead to divergent integrals in gravitational theory as they do in the analogous electromagnetic problem. In the em case, one introduces (normally in an arbitrary manner) the Debye length as a cutoff at large distances to remove the divergence. We show here that in the Newtonian gravitational analogue the distant cutoff problem is solved by including the self-gravity of the medium. Then the Jeans length appears naturally as a cutoff. Despite the different sign of the coupling constant from the case of the electric plasma, the same sort of theory works in both cases and removes the need of introducing ad hoc cutoffs from both of them.

We show that an array of singular density concentrations accounts for the scaling exponents observed in the luminous matter distribution in the Universe for scales below 10Mpc. This model provides a good description of the matter distribution on those smaller scales.