Fulltext

Application of the multi distribution function lattice Boltzmann approach to thermal flows

Parmigiani, A. ; Huber, C. ; Chopard, B. ; Latt, J. ; Bachmann, O.

In: The European Physical Journal Special Topics, 2009, vol. 171, no. 1, p. 37-43

Ajouter à la liste personnelle
    Titre
    • Application of the multi distribution function lattice Boltzmann approach to thermal flows
    Auteur
    • Parmigiani, A.. Computer Science Department, University of Geneva, 24 rue du Général Dufour, 1211, Geneva 4, Switzerland
    • Huber, C.. Department of Earth and Planetary Science, University of California - Berkeley, 307 McCone Hall 4767, 94720-4767, Berkeley, CA, USA
    • Chopard, B.. Computer Science Department, University of Geneva, 24 rue du Général Dufour, 1211, Geneva 4, Switzerland
    • Latt, J.. Department of Mathematics, Tufts University, 503 Boston Avenue, 02155, Medford, MA, USA
    • Bachmann, O.. Department of Earth and Space Science, University of Washington, Johnson Hall, 98195-1310, Seattle, WA, USA
    Type de document
    • Postprint
    Langue
    • Anglais
    Publié dans
    • The European Physical Journal Special Topics, 2009, vol. 171, no. 1, p. 37-43. Springer-Verlag
    Autre version électronique
    Classification
    • Physique
    Identifiant OAI-PMH
    • oai:doc.rero.ch:311407
    Summary
    Numerical methods able to model high Rayleigh (Ra) and high Prandtl (Pr) number thermal convection are important to study large-scale geophysical phenomena occuring in very viscous fluids such as magma chamber dynamics (104 < Pr < 107 and 107 < Ra < 1011). The important variable to quantify the thermal state of a convective fluid is a generalized dimensionless heat transfer coefficient (the Nusselt number) whose measure indicates the relative efficiency of the thermal convection. In this paper we test the ability of Multi-distribution Function approach (MDF) Thermal Lattice Boltzmann method to study the well-established scaling result for the Nusselt number (Nu ∝ Ra 1/3) in Rayleigh Bénard convection for 104 ≤ Ra ≤ 109 and 101 ≤ Pr ≤ 104. We explore its main drawbacks in the range of Pr and Ra number under investigation: (1) high computational time N c required for the algorithm to converge and (2) high spatial accuracy needed to resolve the thickness of thermal plumes and both thermal and velocity boundary layer. We try to decrease the computational demands of the method using a multiscale approach based on the implicit dependence of the Pr number on the relaxation time, the spatial and temporal resolution characteristic of the MDF thermal model