On Pattern Selection in Three-Dimensional Bénard-Marangoni Flows

preprint numerics no. 4/2010 norwegian university of science and technology trondheim, norway In this paper we study Bénard-Marangoni convection in confined containers where a thin fluid layer is heated from below. We consider containers with circular, square and hexagonal cross-sections. For Marangoni numbers close to the critical Marangoni number, the flow patterns are dominated by the appearance of the well-known hexagonal convection cells. The main purpose of this computational study is to demonstrate the importance of the initial conditions on the pattern selection. In a series of numerical experiments, the coupled fluid-thermal system is started with a zero initial condition for the velocity and a random initial condition for the temperature. We demonstrate that the system can end up in more than one state. For example, the final state of the system may be dominated by a steady convection pattern with a fixed number of cells, however, the same system may occasionally end up in a steady pattern involving a slightly different number of cells, or it may end up in a state where most of the cells are stationary, while one or more cells end up in an oscillatory state. For larger aspect ratio containers, we are also able to reproduce dislocations in the convection pattern, which have also been observed experimentally. It has been conjectured that such imperfections (e.g., a localized star-like pattern) are due to small irregularities in the experimental setup (e.g., the geometry of the container). However, we show, through controlled numerical experiments, that such phenomena may appear under otherwise ideal conditions. By repeating the numerical experiments for the same non-dimensional numbers, using a different random initial condition for the temperature in each case, we are able to get an indication of how rare such events are. Next, we study the effect of symmetrizing the initial conditions. Finally, we study the effect of selected geometry deformations on the resulting convection patterns.

[1]  Einar M. Rønquist,et al.  A high order splitting method for time-dependent domains , 2008 .

[2]  Nitschke,et al.  Secondary instability in surface-tension-driven Bénard convection. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[3]  E. Koschmieder,et al.  Bénard cells and Taylor vortices , 1993 .

[4]  J. M. Watt Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[5]  E. Koschmieder,et al.  The wavenumbers of supercritical surface-tension-driven Bénard convection , 1992, Journal of Fluid Mechanics.

[6]  Jie Shen,et al.  On the error estimates for the rotational pressure-correction projection methods , 2003, Math. Comput..

[7]  Pierre Dauby,et al.  Linear Benard-Marangoni instability in rigid circular containers , 1997 .

[8]  Stephen H. Davis,et al.  NONLINEAR MARANGONI CONVECTION IN BOUNDED LAYERS. , 1982 .

[9]  R. Sani,et al.  Finite amplitude bénard-rayleigh convection , 1979 .

[10]  L. Scriven,et al.  On cellular convection driven by surface-tension gradients: effects of mean surface tension and surface viscosity , 1964, Journal of Fluid Mechanics.

[11]  Pierre Dauby,et al.  Benard-Marangoni instability in rigid rectangular containers , 1996 .

[12]  D Maza,et al.  Patterns in small aspect ratio Bénard-Marangoni convection. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[13]  K. A. Smith On convective instability induced by surface-tension gradients , 1966, Journal of Fluid Mechanics.

[14]  Bestehorn Phase and amplitude instabilities for Bénard-Marangoni convection in fluid layers with large aspect ratio. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[15]  Einar M. Rønquist,et al.  An Operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow , 1990 .

[16]  Einar M. Rønquist,et al.  A fast tensor-product solver for incompressible fluid flow in partially deformed three-dimensional domains: Parallel implementation , 2011 .

[17]  M. Block,et al.  Surface Tension as the Cause of Bénard Cells and Surface Deformation in a Liquid Film , 1956, Nature.

[18]  S. H. Davis,et al.  Energy stability theory for free-surface problems: buoyancy-thermocapillary layers , 1980, Journal of Fluid Mechanics.

[19]  M. V. Dyke,et al.  An Album of Fluid Motion , 1982 .

[20]  Rivier,et al.  Topological correlations in Bénard-Marangoni convective structures. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[21]  W. J. Gordon,et al.  Construction of curvilinear co-ordinate systems and applications to mesh generation , 1973 .

[22]  R. A. Wentzell,et al.  Hydrodynamic and Hydromagnetic Stability. By S. CHANDRASEKHAR. Clarendon Press: Oxford University Press, 1961. 652 pp. £5. 5s. , 1962, Journal of Fluid Mechanics.

[23]  Bestehorn Square patterns in Bénard-Marangoni convection. , 1996, Physical review letters.

[24]  Frans N. van de Vosse,et al.  An approximate projec-tion scheme for incompressible ow using spectral elements , 1996 .

[25]  Steven A. Orszag,et al.  Surface-tension-driven Bénard convention at infinite Prandtl number , 1995, Journal of Fluid Mechanics.

[26]  A. Patera,et al.  Spectral element methods for the incompressible Navier-Stokes equations , 1989 .

[27]  S. H. Davis,et al.  On the principle of exchange of stabilities , 1969, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[28]  G. Lebon,et al.  A nonlinear stability analysis of the Bénard–Marangoni problem , 1984, Journal of Fluid Mechanics.

[29]  J. Kan A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[30]  R. Sani,et al.  On finite amplitude cellular convection induced by surface tension , 1974 .

[31]  Lord Rayleigh,et al.  LIX. On convection currents in a horizontal layer of fluid, when the higher temperature is on the under side , 1916 .

[32]  Marc Medale,et al.  NUMERICAL SIMULATION OF BE´NARD-MARANGONI CONVECTION IN SMALL ASPECT RATIO CONTAINERS , 2002 .

[33]  Yvon Maday,et al.  Fast Tensor-Product Solvers: Partially Deformed Three-dimensional Domains , 2009, J. Sci. Comput..

[34]  S. A. Prahl,et al.  Surface-tension-driven Bénard convection in small containers , 1990, Journal of Fluid Mechanics.

[35]  J. Pearson,et al.  On convection cells induced by surface tension , 1958, Journal of Fluid Mechanics.

[36]  Einar M. Rønquist,et al.  Simulation of three-dimensional Benard-Marangoni flows including deformed surfaces , 2008 .