Hostname: page-component-848d4c4894-x24gv Total loading time: 0 Render date: 2024-05-30T23:11:25.932Z Has data issue: false hasContentIssue false

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

Published online by Cambridge University Press:  20 August 2015

Arne Morten Kvarving*
Department of Mathematical Sciences, Norwegian University of Science and Technology (NTNU), Norway
Tormod Bjøntegaard*
Department of Mathematical Sciences, Norwegian University of Science and Technology (NTNU), Norway
Einar M. Rønquist*
Department of Mathematical Sciences, Norwegian University of Science and Technology (NTNU), Norway
Get access


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 explore the possible patterns the system may end up in for a given set of parameters. 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. For a given set of parameters 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.

Research Article
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)


[1]Arrow, K., Hurwicz, L., and Uzawa, H.Studies in Nonlinear programming. Standford University Press, 1958.Google Scholar
[2]Bénard, H.Les tourbillons cellulaires dans une nappe liquide transportant de la chaleur par convection en régime permanent. Annales de Chimie et de Physique, 23:62144, 1901.Google Scholar
[3]Bestehorn, M.Phase and amplitude instabilities for Bénard-Marangoni convection in fluid layers with large aspect ratio. Physical Review E, 48:36223634, 1993.Google Scholar
[4]Bestehorn, M.Square patterns in Bénard-Marangoni convection. Physical Review Letters, 76:4649, 1996.Google Scholar
[5]Bjøntegaard, T., Maday, Y., and Rønquist, E. M.Fast tensor-product solvers: Partially deformed three-dimensional domains. J. Sci. Comput., 39:2848, 2009.Google Scholar
[6]Bjøntegaard, T. and Rønquist, E. M.A high order splitting method for time-dependent domains. Comput. Meth. Appl. Mech. Engr., 197:47634773, 2008.Google Scholar
[7]Bjøntegaard, T. and Rønquist, E.M.Simulation of three-dimensional Bénard-Marangoni flows including deformed surfaces. Commun. Comput. Phys., 5(24):273295, 2009.Google Scholar
[8]Block, M. J.Surface tension as the cause of Bénard cells and surface deformation in a liquid film. Nature, 178:650651, 1956.Google Scholar
[9]Bragard, J. and Lebon, G.Non-linear Marangoni convection in a layer of finite depth. Europhys. Lett., pages 831838, 1993.Google Scholar
[10]Cerisier, P., Rahal, S., and Billia, B.Extrinsic effects on the disorder dynamics of Bénard-Marangoni patterns. Physical Review E, 54:35083517, 1996.Google Scholar
[11]Cerisier, P., Rahal, S., and Rivier, N.Topological correlations in Bénard-Marangoni convective structures. Physical Review E, 54:50865094, 1996.Google Scholar
[12]Chandrasekhar, S.Hydrodynamic and Hydromagnetic Stability. Dover Publications, Inc., 1961.Google Scholar
[13]Cloot, A. and Lebon, G.A nonlinear stability analysis of the Bénard-Marangoni problem. J. Fluid Mech., 145:447469, 1984.Google Scholar
[14]Dauby, P. C. and Lebon, G.Bénard-Marangoni instability in rigid rectangular containers. J. Fluid Mech., 329:2564, 1996.Google Scholar
[15]Dauby, P. C., Lebon, G., and Bouhy, E.Linear Bénard-Marangoni instability in rigid circular containers. Physical Review E, 56(1), 1997.Google Scholar
[16]Davis, S. H. and Homsy, G.Energy stability for free-surface problems: Buoyancy-Thermocapillary layers. J. Fluid Mech., 98:527, 1980.Google Scholar
[17]Van Dyke, M.An Album of Fluid Motion. Parabolic Press, Stanford, Ca., 1982.Google Scholar
[18]Gear, C. W.Numerical Initial value problems in ordinary differential equations. Prentice-Hall, 1971.Google Scholar
[19]Gordon, W. J. and Hall, C. A.Construction of curvilinear co-ordinate systems and applications to mesh generation. Intl. J. Numer. Meth. Eng, 7:461477, 1973.Google Scholar
[20]Guermond, J. L. and Shen, J.On the error estimates for the rotational pressure-correction projection methods. Math. Comp., 73:17191737, 2003.Google Scholar
[21]Koschmieder, E. L.Bénard convection. Adv. Chem. Phys., 26:177212, 1974.Google Scholar
[22]Koschmieder, E. L.Bénard Cells and Taylor Vortices. Cambridge University Press, 1993.Google Scholar
[23]Koschmieder, E. L. and Switzer, D. W.The wavenumbers of supercritical surface-tension-driven Bénard convection. J. Fluid Mech., 240:533548, 1992.Google Scholar
[24]Koschmieder, E. L. and Prahl, S.A.Surface-tension-driven Bénard convection in small containers. J. Fluid Mech., 215:571583, 1990.Google Scholar
[25]Kraska, J. and Sani, R.Finite amplitude Bénard-Rayleigh convection. Intl. J. Heat Mass Transfer, 22:535546, 1979.Google Scholar
[26]Kvarving, A. M.A fast tensor-product solver for incompressible fluid flow in partially deformed three-dimensional domains: Parallel implementation. Submitted to Computers & Fluids.Google Scholar
[27]Kvarving, A. M., Bjøntegaard, T., and Rønquist, E. M.A fast tensor-product solver for incom-pressible fluid flow in partially deformed three-dimensional domains. Technical Report 9, Department for Mathematical Sciences, Numerical Group, NTNU, 2010.Google Scholar
[28]Maday, Y. and Patera, A. T.Spectral element methods for the incompressible Navier-Stokes equations. In State-of-the-art surveys on computational mechanics (A90-47176 21-64). New York, American Society of Mechanical Engineers, pages 71143, 1989.Google Scholar
[29]Maday, Y., Patera, A.T., and Rønquist, E. M.An Operator-Integration-Factor Splitting Method for Time-Dependent Problems: Application to Incompressible Fluid Flow. J. Sci. Comput., 5(4), 1990.Google Scholar
[30]Medale, M. and Cerisier, P.Numerical simulation of Bénard-Marangoni convection in small aspect ratio containers. Numerical Heat Transfer, Part A, 42:5572, 2002.Google Scholar
[31]Nitschke, K. and Thess, A.Secondary instability in surface-tension-driven Bénard convection. Physical Review E, 52:57725775, 1995.Google Scholar
[32]Pearson, J. R. A.On convection cells induced by surface tension. J. Fluid Mech., 4:489500, 1958.Google Scholar
[33]Ramón, M. L., Maza, D., and Mancini, H. L.Patterns in small aspect ratio Bénard-Marangoni convection. Physical Review E, 60, 1999.Google Scholar
[34]Rayleigh, Lord. On convection currents in a horizontal layer of fluid when the higher temperature is on the under side. Phil. Mag., 32:529546, 1916.Google Scholar
[35]Rosenblat, S., Davis, S. H., and Homsy, G. M.Nonlinear Marangoni convection in bounded layers – 1. circular cylindrical containers. J. Fluid Mech., 120, 1982.Google Scholar
[36]Scanlon, J. W. and Segel, L. A.Finite amplitude cellular convection induced by surface tension. Journal of Fluid Mechanics, 30:149162, 1967.Google Scholar
[37]Schriven, L. and Sternling, C.On cellular convection driven by surface-tension gradients; effects of mean surface tension and surface viscosity. J. Fluid Mech., 21:321340, 1964.Google Scholar
[38]Smith, K. A.On convective instability induced by surface-tension gradients. J. Fluid Mech., 24:401414, 1966.Google Scholar
[39]Thess, A. and Orszag, S. A.Surface-tension-driven Bénard convection at infinite Prandtl number. J. Fluid Mech., 283:201230, 1995.Google Scholar
[40]Timmermans, L. J. P., Minev, P. D., and Van De Vosse, F. N.An approximate projection scheme for incompressible flow using spectral elements. Intl. J. Numer. Met. Fluids, 22, 1996.Google Scholar
[41]Tritton, D. J.Physical Fluid Dynamics. Oxford Science Publications, 1988.Google Scholar
[42]van Kan, J.A second-order accurate pressure-correction scheme for viscous incompressible flow. J. Sci. Stat. Comput., 3, 1986.Google Scholar
[43]Yu, S. T., Jiang, B. N., Wu, J., and Duh, J. C.Three-dimensional simulations of Marangoni-Bénard convection in small containers by the least-squares finite element method. Comput. Meth. Appl. Mech. Engr., 160:7188, 1998.Google Scholar