Wavy stripes and squares in zero P number convection

A simple model to explain numerically observed behaviour of chaotically varying stripes and square patterns in zero Prandtl number convection in Boussinesq fluid is presented. The nonlinear interaction of mutually perpendicular sets of wavy rolls, via higher mode, may lead to a competition between the two sets of wavy rolls. The appearance of square patterns is due to the secondary forward Hopf bifurcation of a set of wavy rolls.Comment: 8 pages and 3 figures, late

Asymmetric Squares as Standing Waves in Rayleigh-Benard Convection

Possibility of asymmetric square convection is investigated numerically using a few mode Lorenz-like model for thermal convection in Boussinesq fluids confined between two stress free and conducting flat boundaries. For relatively large value of Rayleigh number, the stationary rolls become unstable and asymmetric squares appear as standing waves at the onset of secondary instability. Asymmetric squares, two dimensional rolls and again asymmetric squares with their corners shifted by half a wavelength form a stable limit cycle.Comment: 8 pages, 7 figure

Cross-Newell equations for hexagons and triangles

The Cross-Newell equations for hexagons and triangles are derived for general real gradient systems, and are found to be in flux-divergence form. Specific examples of complex governing equations that give rise to hexagons and triangles and which have Lyapunov functionals are also considered, and explicit forms of the Cross-Newell equations are found in these cases. The general nongradient case is also discussed; in contrast with the gradient case, the equations are not flux-divergent. In all cases, the phase stability boundaries and modes of instability for general distorted hexagons and triangles can be recovered from the Cross-Newell equations.Comment: 24 pages, 1 figur

Convection in colloidal suspensions with particle-concentration-dependent viscosity

The onset of thermal convection in a horizontal layer of a colloidal suspension is investigated in terms of a continuum model for binary-fluid mixtures where the viscosity depends on the local concentration of colloidal particles. With an increasing difference between the viscosity at the warmer and the colder boundary the threshold of convection is reduced in the range of positive values of the separation ratio psi with the onset of stationary convection as well as in the range of negative values of psi with an oscillatory Hopf bifurcation. Additionally the convection rolls are shifted downwards with respect to the center of the horizontal layer for stationary convection (psi>0) and upwards for the Hopf bifurcation (psi<0).Comment: 8 pages, 6 figures, submitted to European Physical Journal

Efficient Algorithm on a Non-staggered Mesh for Simulating Rayleigh-Benard Convection in a Box

An efficient semi-implicit second-order-accurate finite-difference method is described for studying incompressible Rayleigh-Benard convection in a box, with sidewalls that are periodic, thermally insulated, or thermally conducting. Operator-splitting and a projection method reduce the algorithm at each time step to the solution of four Helmholtz equations and one Poisson equation, and these are are solved by fast direct methods. The method is numerically stable even though all field values are placed on a single non-staggered mesh commensurate with the boundaries. The efficiency and accuracy of the method are characterized for several representative convection problems.Comment: REVTeX, 30 pages, 5 figure