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

Black Hole Boundary Conditions and Coordinate Conditions

This paper treats boundary conditions on black hole horizons for the full 3+1D Einstein equations. Following a number of authors, the apparent horizon is employed as the inner boundary on a space slice. It is emphasized that a further condition is necessary for the system to be well posed; the ``prescribed curvature conditions" are therefore proposed to complete the coordinate conditions at the black hole. These conditions lead to a system of two 2D elliptic differential equations on the inner boundary surface, which coexist nicely to the 3D equation for maximal slicing (or related slicing conditions). The overall 2D/3D system is argued to be well posed and globally well behaved. The importance of ``boundary conditions without boundary values" is emphasized. This paper is the first of a series. This revised version makes minor additions and corrections to the previous version.Comment: 13 pages LaTeX, revtex. No figure

Resonance bifurcations of robust heteroclinic networks

Robust heteroclinic cycles are known to change stability in resonance bifurcations, which occur when an algebraic condition on the eigenvalues of the system is satisfied and which typically result in the creation or destruction of a long-period periodic orbit. Resonance bifurcations for heteroclinic networks are more complicated because different subcycles in the network can undergo resonance at different parameter values, but have, until now, not been systematically studied. In this article we present the first investigation of resonance bifurcations in heteroclinic networks. Specifically, we study two heteroclinic networks in $\R^4$ and consider the dynamics that occurs as various subcycles in each network change stability. The two cases are distinguished by whether or not one of the equilibria in the network has real or complex contracting eigenvalues. We construct two-dimensional Poincare return maps and use these to investigate the dynamics of trajectories near the network. At least one equilibrium solution in each network has a two-dimensional unstable manifold, and we use the technique developed in [18] to keep track of all trajectories within these manifolds. In the case with real eigenvalues, we show that the asymptotically stable network loses stability first when one of two distinguished cycles in the network goes through resonance and two or six periodic orbits appear. In the complex case, we show that an infinite number of stable and unstable periodic orbits are created at resonance, and these may coexist with a chaotic attractor. There is a further resonance, for which the eigenvalue combination is a property of the entire network, after which the periodic orbits which originated from the individual resonances may interact. We illustrate some of our results with a numerical example.Comment: 46 pages, 20 figures. Supplementary material (two animated gifs) can be found on http://www.maths.leeds.ac.uk/~alastair/papers/KPR_res_net_abs.htm

Traveling waves in rotating Rayleigh-Bénard convection: Analysis of modes and mean flow

Numerical simulations of the Boussinesq equations with rotation for realistic no-slip boundary conditions and a finite annular domain are presented. These simulations reproduce traveling waves observed experimentally. Traveling waves are studied near threshhold by using the complex Ginzburg-Landau equation (CGLE): a mode analysis enables the CGLE coefficients to be determined. The CGLE coefficients are compared with previous experimental and theoretical results. Mean flows are also computed and found to be more significant as the Prandtl number decreases (from sigma=6.4 to sigma=1). In addition, the mean flow around the outer radius of the annulus appears to be correlated with the mean flow around the inner radius

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