Lyapunov exponents for small aspect ratio Rayleigh-Bénard convection

Leading order Lyapunov exponents and their corresponding eigenvectors have been computed numerically for small aspect ratio, three-dimensional Rayleigh-Benard convection cells with no-slip boundary conditions. The parameters are the same as those used by Ahlers and Behringer [Phys. Rev. Lett. 40, 712 (1978)] and Gollub and Benson [J. Fluid Mech. 100, 449 (1980)] in their work on a periodic time dependence in Rayleigh-Benard convection cells. Our work confirms that the dynamics in these cells truly are chaotic as defined by a positive Lyapunov exponent. The time evolution of the leading order Lyapunov eigenvector in the chaotic regime will also be discussed. In addition we study the contributions to the leading order Lyapunov exponent for both time periodic and aperiodic states and find that while repeated dynamical events such as dislocation creation/annihilation and roll compression do contribute to the short time Lyapunov exponent dynamics, they do not contribute to the long time Lyapunov exponent. We find instead that nonrepeated events provide the most significant contribution to the long time leading order Lyapunov exponent

Scaling laws for rotating Rayleigh-Bénard convection

Numerical simulations of large aspect ratio, three-dimensional rotating Rayleigh-Bénard convection for no-slip boundary conditions have been performed in both cylinders and periodic boxes. We have focused near the threshold for the supercritical bifurcation from the conducting state to a convecting state exhibiting domain chaos. A detailed analysis of these simulations has been carried out and is compared with experimental results, as well as predictions from multiple scale perturbation theory. We find that the time scaling law agrees with the theoretical prediction, which is in contradiction to experimental results. We also have looked at the scaling of defect lengths and defect glide velocities. We find a separation of scales in defect diameters perpendicular and parallel to the rolls as expected, but the scaling laws for the two different lengths are in contradiction to theory. The defect velocity scaling law agrees with our theoretical prediction from multiple scale perturbation theory

THE "GREAT DECISIONS" PROGRAM IN OREGON

Teaching/Communication/Extension/Profession,

Resolving the fine-scale structure in turbulent Rayleigh-Benard convection

We present high-resolution direct numerical simulation studies of turbulent Rayleigh-Benard convection in a closed cylindrical cell with an aspect ratio of one. The focus of our analysis is on the finest scales of convective turbulence, in particular the statistics of the kinetic energy and thermal dissipation rates in the bulk and the whole cell. The fluctuations of the energy dissipation field can directly be translated into a fluctuating local dissipation scale which is found to develop ever finer fluctuations with increasing Rayleigh number. The range of these scales as well as the probability of high-amplitude dissipation events decreases with increasing Prandtl number. In addition, we examine the joint statistics of the two dissipation fields and the consequences of high-amplitude events. We also have investigated the convergence properties of our spectral element method and have found that both dissipation fields are very sensitive to insufficient resolution. We demonstrate that global transport properties, such as the Nusselt number, and the energy balances are partly insensitive to insufficient resolution and yield correct results even when the dissipation fields are under-resolved. Our present numerical framework is also compared with high-resolution simulations which use a finite difference method. For most of the compared quantities the agreement is found to be satisfactory.Comment: 33 pages, 24 figure

Gauge drivers for the generalized harmonic Einstein equations

The generalized harmonic representation of Einstein's equations is manifestly hyperbolic for a large class of gauge conditions. Unfortunately most of the useful gauges developed over the past several decades by the numerical relativity community are incompatible with the hyperbolicity of the equations in this form. This paper presents a new method of imposing gauge conditions that preserves hyperbolicity for a much wider class of conditions, including as special cases many of the standard ones used in numerical relativity: e.g., K freezing, Gamma freezing, Bona-Massó slicing, conformal Gamma drivers, etc. Analytical and numerical results are presented which test the stability and the effectiveness of this new gauge-driver evolution system

Quantum local-field corrections and spontaneous decay

A recently developed scheme [S. Scheel, L. Knoll, and D.-G. Welsch, Phys. Rev. A 58, 700 (1998)] for quantizing the macroscopic electromagnetic field in linear dispersive and absorbing dielectrics satisfying the Kramers-Kronig relations is used to derive the quantum local-field correction for the standard virtual-sphere-cavity model. The electric and magnetic local-field operators are shown to be consistent with QED only if the polarization noise is fully taken into account. It is shown that the polarization fluctuations in the local field can dramatically change the spontaneous decay rate, compared with the familiar result obtained from the classical local-field correction. In particular, the spontaneous emission rate strongly depends on the radius of the local-field virtual cavity.Comment: 7 pages, using RevTeX, 4 figure