Direct numerical simulation of supercritical annular electroconvection

We use direct numerical simulation to study electrically driven convection in an annular thin film. The simulation models a laboratory experiment that consists of a weakly conducting, submicron thick liquid crystal film suspended between two concentric electrodes. The film is driven to convect by imposing a sufficiently large voltage across it. The flow is driven by a surface charge density inversion which is unstable to the imposed electrical force. This mechanism is closely analogous to the mass density inversion which is unstable to the buoyancy force in conventional, thermally driven Rayleigh-Bénard convection. The simulation uses a pseudospectral method with Chebyshev polynomials in the radial direction and Fourier modes in the azimuthal direction. The numerical results, which are in good agreement with previous experimental data and theoretical predictions, reveal several insights. The mode competition near a codimension-two point exhibits hysteresis. The primary bifurcation is supercritical for a broad range of fluid and geometrical parameters.