jeanflop/post_ocr_correction2 huggingface dataset:
problem of poorly conducting horizontal boundaries. One may just note that such a bifurcation does not agree with the linear theoretical work of Edwards (1988) on the effect of small aspect ratios, in the case of infinitely conducting horizontal walls, and, consequently, suggests that the effect of low conductivity of the horizontal walls is significant. The square pattern configuration of the flow field obtained for higher Rayleigh numbers also does not agree with the theoretical work of Edwards. In addition, it appears to be stable even far from the onset of convection (at least 16 Ra), illustrating the large effect of the finite conductivity of the walls on the cell stability. Further calculations are currently being carried out for a larger range of c values in order, first, to compare the evolution of the corresponding critical Rayleigh numbers with those predicted by the linear stability analysis and, second, to obtain an estimate of the bounding value of L, for roll pattern stability. Additionally, more information on the influence of the lateral walls on the flow pattern in the vicinity of the convective threshold would be obtained. 1. De nombreux travaux théoriques et numériques ont été réalisés pour caractériser la convection de Rayleigh-Bénard se développant dans une couche fluide située entre deux plans parallèles horizontaux à températures fixées. L'influence de la conductivité thermique des parois délimitant la couche fluide sur la valeur du nombre de Rayleigh et du nombre d'onde critiques a été notamment mise en évidence par Sparrow, Goldstein et Jonsson (1964). S'attachant plus particulièrement au rôle des parois horizontales, l'étude analytique de stabilité linéaire de Proctor (1981) pour une couche d'extension infinie prédit que la valeur du nombre de Rayleigh critique Rac diminue avec le rapport c entre la conductivité thermique des parois solides et celle du fluide. En particulier, la valeur asymptotique de Rac prédite lorsque c tend vers 0 est 720. De même, la valeur du nombre d'onde critique ac diminue avec c et tend vers 0 à la limite. L'étude analytique en convection faiblement non linéaire, plus récente, de Jenkins et Proctor (1984) montre en outre
de Rayleigh-Bénard se développant dans une couche fluide située entre deux plans parallèles horizontaux à températures fixées. L'influence de la conductivité thermique des parois délimitant la couche fluide sur la valeur du nombre de Rayleigh et du nombre d'onde critiques a été notamment mise en évidence par Sparrow, Goldstein et Jonsson (1964). S'attachant plus particulièrement au rôle des parois horizontales, l'étude analytique de stabilité linéaire de Proctor (1981) pour une couche d'extension infinie prédit que la valeur du nombre de Rayleigh critique Rac diminue avec le rapport c entre la conductivité thermique des parois solides et celle du fluide. En particulier, la valeur asymptotique de Rac prédite lorsque c tend vers 0 est 720. De même, la valeur du nombre d'onde critique ac diminue avec c et tend vers 0 à la limite. L'étude analytique en convection faiblement non linéaire, plus récente, de Jenkins et Proctor (1984) montre en outre qu'au voisinage du seuil de convection et pour des valeurs de c suffisamment faibles paramétrées par le nombre de Prandtl, la structure convective classique en rouleaux bidimensionnels est instable, l'écoulement se présente sous la forme de cellules carrées résultant de la présence de rouleaux orthogonaux identiques. Ainsi, pour des nombres de Prandtl supérieurs à l'unité, le changement de stabilité intervient lorsque c est inférieur à une valeur limite voisine de 1. À notre connaissance, de telles structures, pour lesquelles la stabilité serait reliée à la faible conductivité des parois horizontales, n'ont été observées que dans l'étude expérimentale effectuée par Le Gal et Croquette (1988). 2. La présente étude a pour but d'étudier la stabilité des cellules carrées par la simulation numérique. On considère pour cela une couche fluide d'épaisseur h comprise entre deux plaques parallèles horizontales conductrices identiques et d'épaisseur k"h/2. Les allongements dans les directions horizontales x et y du domaine étudié sont notés Ax h et Ay h respectivement. Par la suite, l'origine des coordonnées est prise au centre de la couche fluide. Les frontières verticales délimitant le domaine sont supposées solides et parfaitement isolées. Les caractéristiques physiques du fluide et celles du solide constituant les plaques horizontales
de Prandtl supérieurs à l'unité, le changement de stabilité intervient lorsque c est inférieur à une valeur limite voisine de 1. À notre connaissance, de telles structures, pour lesquelles la stabilité serait reliée à la faible conductivité des parois horizontales, n'ont été observées que dans l'étude expérimentale effectuée par Le Gal et Croquette (1988). 2. La présente étude a pour but d'étudier la stabilité des cellules carrées par la simulation numérique. On considère pour cela une couche fluide d'épaisseur h comprise entre deux plaques parallèles horizontales conductrices identiques et d'épaisseur k"h/2. Les allongements dans les directions horizontales x et y du domaine étudié sont notés Ax h et Ay h respectivement. Par la suite, l'origine des coordonnées est prise au centre de la couche fluide. Les frontières verticales délimitant le domaine sont supposées solides et parfaitement isolées. Les caractéristiques physiques du fluide et celles du solide constituant les plaques horizontales sont supposées constantes. En adoptant les approximations classiques de Boussinesq, on a alors, dans le milieu fluide, le système d'équations gouvernant la vitesse u, la pression p et la déviation de température par rapport à l'état de référence en l'absence de mouvement ef : tandis que dans le milieu solide, le problème se réduit à une équation de conduction sur la déviation de température es : où KS et Kf sont les diffusivités thermiques du milieu solide et du milieu fluide respectivement. On applique la condition d'adhérence sur l'ensemble des frontières limitant la couche fluide. Les températures sont supposées uniformes sur les frontières extérieures des plaques horizontales et l'adiabaticité est appliquée sur les frontières verticales. Sur les frontières internes des plaques solides, les conditions de continuité de la température et du flux de chaleur conduisent à : Le nombre de Rayleigh est ici défini à partir du gradient de température vertical traversant la couche fluide en l'absence de mouvement. La méthode numérique utilisée est une méthode d'éléments finis en espace couplée à une méthode de différences finies à pas fractionnaires en temps (voir Carrière et Jeandel, 1991) permettant le calcul simultané des solutions approchées dans le milieu fluide et le
le problème se réduit à une équation de conduction sur la déviation de température es : où KS et Kf sont les diffusivités thermiques du milieu solide et du milieu fluide respectivement. On applique la condition d'adhérence sur l'ensemble des frontières limitant la couche fluide. Les températures sont supposées uniformes sur les frontières extérieures des plaques horizontales et l'adiabaticité est appliquée sur les frontières verticales. Sur les frontières internes des plaques solides, les conditions de continuité de la température et du flux de chaleur conduisent à : Le nombre de Rayleigh est ici défini à partir du gradient de température vertical traversant la couche fluide en l'absence de mouvement. La méthode numérique utilisée est une méthode d'éléments finis en espace couplée à une méthode de différences finies à pas fractionnaires en temps (voir Carrière et Jeandel, 1991) permettant le calcul simultané des solutions approchées dans le milieu fluide et le milieu solide. Les éléments finis sont des tétraèdres P1 (approximation linéaire); le maillage utilisé contient 33 145 nœuds, la couche fluide proprement dite étant décrite par (41)2.17 nœuds où sont calculées les trois composantes de vitesse et la température. 3. La configuration présentée ici a été choisie de façon à se rapprocher du cas expérimental étudié par Le Gal et Croquette. Le rapport des conductivités est égal à 0,4, le nombre de Prandtl est égal à 7. L'extension verticale des murs conducteurs est λ = 4. Les principales différences entre le présent calcul et l'expérience (voir Le Gal et Croquette, 1988) proviennent de la nécessité de considérer un domaine de calcul de dimension finie. Compte tenu des limitations sur le plan numérique quand au nombre de nœuds du maillage, il n'a pas été possible de prendre en compte une grande extension horizontale tout en conservant une description fine de chacune des structures cellulaires de l'écoulement. La couche fluide choisie a une section horizontale carrée dont la longueur du côté est limitée à environ trois fois le diamètre des cellules déterminé par l'étude de stabilité linéaire, soit Ax = λy = 4,3. Le premier calcul a été effectué pour une valeur
se rapprocher du cas expérimental étudié par Le Gal et Croquette. Le rapport des conductivités est égal à 0,4, le nombre de Prandtl est égal à 7. L'extension verticale des murs conducteurs est λ = 4. Les principales différences entre le présent calcul et l'expérience (voir Le Gal et Croquette, 1988) proviennent de la nécessité de considérer un domaine de calcul de dimension finie. Compte tenu des limitations sur le plan numérique quand au nombre de nœuds du maillage, il n'a pas été possible de prendre en compte une grande extension horizontale tout en conservant une description fine de chacune des structures cellulaires de l'écoulement. La couche fluide choisie a une section horizontale carrée dont la longueur du côté est limitée à environ trois fois le diamètre des cellules déterminé par l'étude de stabilité linéaire, soit Ax = λy = 4,3. Le premier calcul a été effectué pour une valeur du nombre de Rayleigh de 2100 (correspondant à environ deux fois le Rayleigh critique obtenu par l'analyse de stabilité linéaire pour une extension horizontale infinie) à partir d'un champ initial correspondant à un champ statique auquel est superposée une faible (0,1 %) perturbation aléatoire du champ thermique. Le champ obtenu pour cette valeur du nombre de Rayleigh est stationnaire. La figure 1 montre deux projections du champ de vitesse dans les plans médians verticaux. En dehors des régions proches des parois latérales et indépendamment des interactions non linéaires entre structures, on peut interpréter le champ résultant comme la superposition de trois rouleaux dans chacune des directions horizontales. Une coupe du champ de température dans le plan médian horizontal est présentée sur la figure 2. Les maximums et minimums relatifs de température sont situés à environ 1/3 et 2/3 de la largeur ce qui correspond à des axes de rouleaux situés à 1/6, 1/2 et 5/6 de la largeur de la cavité. D'autre part, le champ thermique dans le plan médian horizontal présente une symétrie par rapport à la diagonale x=y et une antisymétrie par rapport à la diagonale x= - y. L'ensemble de ces résultats apparaît ainsi cohérent avec la
pour cette valeur du nombre de Rayleigh est stationnaire. La figure 1 montre deux projections du champ de vitesse dans les plans médians verticaux. En dehors des régions proches des parois latérales et indépendamment des interactions non linéaires entre structures, on peut interpréter le champ résultant comme la superposition de trois rouleaux dans chacune des directions horizontales. Une coupe du champ de température dans le plan médian horizontal est présentée sur la figure 2. Les maximums et minimums relatifs de température sont situés à environ 1/3 et 2/3 de la largeur ce qui correspond à des axes de rouleaux situés à 1/6, 1/2 et 5/6 de la largeur de la cavité. D'autre part, le champ thermique dans le plan médian horizontal présente une symétrie par rapport à la diagonale x=y et une antisymétrie par rapport à la diagonale x= - y. L'ensemble de ces résultats apparaît ainsi cohérent avec la description habituelle des cellules carrées, et notamment avec les observations expérimentales de Le Gal et Croquette (1988). 4. La recherche du seuil critique s'effectue par calculs successifs en imposant des valeurs décroissantes du nombre de Rayleigh jusqu'à obtenir une solution dont l'amplitude tende vers zéro. Une telle approche, pour une détermination directe du Rayleigh critique, est nécessairement limitée en précision car le taux de croissance ou de décroissance de l'amplitude de l'écoulement tend vers zéro au voisinage du seuil. La voie classique est d'interpoler la valeur critique de Ra à partir de la loi d'évolution d'une quantité caractéristique de l'amplitude de l'écoulement au voisinage du seuil; rappelons que les travaux expérimentaux de Dubois et Bergé (1978) ont montré que, pour une configuration classique (γx, γy, 1), l'amplitude de la vitesse suit une loi de puissance en (Ra-Rac)1/2. Pour notre part, compte tenu de la non unicité du type de structure rencontré qui est décrite par la suite, seule une quantité moyennée spatialement a pu être considérée. On présente figure 3 la variation de l'énergie cinétique moyenne sans dimension u dans l'écoulement en fonction du nombre de Rayleigh. On note qu'elle vérifie remarquablement une loi linéaire au-delà du seuil critique. La
détermination directe du Rayleigh critique, est nécessairement limitée en précision car le taux de croissance ou de décroissance de l'amplitude de l'écoulement tend vers zéro au voisinage du seuil. La voie classique est d'interpoler la valeur critique de Ra à partir de la loi d'évolution d'une quantité caractéristique de l'amplitude de l'écoulement au voisinage du seuil; rappelons que les travaux expérimentaux de Dubois et Bergé (1978) ont montré que, pour une configuration classique (γx, γy, 1), l'amplitude de la vitesse suit une loi de puissance en (Ra-Rac)1/2. Pour notre part, compte tenu de la non unicité du type de structure rencontré qui est décrite par la suite, seule une quantité moyennée spatialement a pu être considérée. On présente figure 3 la variation de l'énergie cinétique moyenne sans dimension u dans l'écoulement en fonction du nombre de Rayleigh. On note qu'elle vérifie remarquablement une loi linéaire au-delà du seuil critique. La valeur interpolée du nombre de Rayleigh critique est 1233 (ce qui correspond à un niveau de 15 supérieur à celui prédit par l'analyse de stabilité linéaire pour une extension infinie). Les résultats expérimentaux de Le Gal et Croquette (1988) se situent dans le même ordre de grandeur. Rappelons qu'il est par ailleurs connu que le fort confinement horizontal utilisé ici a pour effet d'augmenter la valeur du nombre de Rayleigh critique (voir Catton, 1972). Il faut également remarquer que, la largeur de la cavité étant imposée, la valeur du nombre d'onde critique est, bien évidemment, fortement contrainte : il apparaît néanmoins que le nombre de cellules obtenues reste cohérent avec celui prédit par l'analyse linéaire. 5. L'évolution des cellules carrées pour des nombres de Rayleigh décroissants ne montre pas d'évolution particulière excepté au voisinage du seuil critique où apparaît une bifurcation vers une structure différente, de tendance torique, qui est présentée sur les figures 4 et 5. Cette structure présente une parfaite symétrie par rapport aux plans médians verticaux, le fluide froid descendant par le centre de la cavité et le fluide chaud remontant par les bords. Le diamètre de cette cellule est 1,5 fois celui des cellules carrées. La
par ailleurs connu que le fort confinement horizontal utilisé ici a pour effet d'augmenter la valeur du nombre de Rayleigh critique (voir Catton, 1972). Il faut également remarquer que, la largeur de la cavité étant imposée, la valeur du nombre d'onde critique est, bien évidemment, fortement contrainte : il apparaît néanmoins que le nombre de cellules obtenues reste cohérent avec celui prédit par l'analyse linéaire. 5. L'évolution des cellules carrées pour des nombres de Rayleigh décroissants ne montre pas d'évolution particulière excepté au voisinage du seuil critique où apparaît une bifurcation vers une structure différente, de tendance torique, qui est présentée sur les figures 4 et 5. Cette structure présente une parfaite symétrie par rapport aux plans médians verticaux, le fluide froid descendant par le centre de la cavité et le fluide chaud remontant par les bords. Le diamètre de cette cellule est 1,5 fois celui des cellules carrées. La question de l'origine d'une telle bifurcation est à l'heure actuelle ouverte, mais son existence dans cette étude pose le problème de l'influence des parois verticales sur les types de structures rencontrées. Rappelons que des structures particulières peuvent être observées pour des cavités de section carrée comme l'on a montré Storck et Müller (1972), notamment lorsque le confinement est fort. Edwards (1988) a proposé une approche linéaire théorique de l'effet des parois latérales sur la forme des cellules au seuil de convection pour cela et a mis en évidence la présence de rouleaux croisés. On peut noter que la structure prédite par cette analyse correspondrait plutôt à des cellules non symétriques formées par la présence de trois rouleaux dans une direction et deux dans l'autre. Concernant la structure quasi-torique prédite ici, par la simulation, au voisinage du seuil, il semble donc naturel d'invoquer, en plus de l'influence propre du confinement, l'effet de la faible conductivité des parois horizontales. En s'éloignant du seuil critique, l'influence de la condition d'adhérence aux parois latérales sur la structuration de la partie interne du domaine fluide diminue, ce qui expliquerait l'émergence de cellules carrées dont la stabilité serait bien alors reliée de façon prépondérante à la
et Müller (1972), notamment lorsque le confinement est fort. Edwards (1988) a proposé une approche linéaire théorique de l'effet des parois latérales sur la forme des cellules au seuil de convection pour cela et a mis en évidence la présence de rouleaux croisés. On peut noter que la structure prédite par cette analyse correspondrait plutôt à des cellules non symétriques formées par la présence de trois rouleaux dans une direction et deux dans l'autre. Concernant la structure quasi-torique prédite ici, par la simulation, au voisinage du seuil, il semble donc naturel d'invoquer, en plus de l'influence propre du confinement, l'effet de la faible conductivité des parois horizontales. En s'éloignant du seuil critique, l'influence de la condition d'adhérence aux parois latérales sur la structuration de la partie interne du domaine fluide diminue, ce qui expliquerait l'émergence de cellules carrées dont la stabilité serait bien alors reliée de façon prépondérante à la faible conductivité des murs horizontaux en accord avec les analyses en convection faiblement non linéaires. Notons que les cellules carrées obtenues sont stables même loin du Rayleigh critique (au moins 16 Rac pour les calculs déjà effectués). 6. A l'heure actuelle, des calculs complémentaires sont effectués pour couvrir une gamme plus large de valeurs de ç. Ces calculs devraient permettre, en premier lieu, de caractériser l'évolution du nombre de Rayleigh critique en fonction de ç qui sera comparée à celle prédite par l'analyse de stabilité linéaire pour une extension infinie. De plus, il sera alors possible de distinguer les différentes structures stables au-delà du seuil critique et de rechercher ainsi un accord entre simulation numérique et analyse en convection faiblement non linéaire sur la valeur de ç limitant le domaine de stabilité des rouleaux bidimensionnels. D'autre part, cette étude permettra de mieux séparer les rôles respectifs du confinement latéral et celui de la conductivité des murs horizontaux sur le type de structures rencontrées dans le proche voisinage du seuil de convection; on pourra également compléter l'étude par des simulations pour des géométries de sections non carrées privilégiant la formation de structure en rouleaux de façon à évaluer la stabilité intrinsèque
plus large de valeurs de ç. Ces calculs devraient permettre, en premier lieu, de caractériser l'évolution du nombre de Rayleigh critique en fonction de ç qui sera comparée à celle prédite par l'analyse de stabilité linéaire pour une extension infinie. De plus, il sera alors possible de distinguer les différentes structures stables au-delà du seuil critique et de rechercher ainsi un accord entre simulation numérique et analyse en convection faiblement non linéaire sur la valeur de ç limitant le domaine de stabilité des rouleaux bidimensionnels. D'autre part, cette étude permettra de mieux séparer les rôles respectifs du confinement latéral et celui de la conductivité des murs horizontaux sur le type de structures rencontrées dans le proche voisinage du seuil de convection; on pourra également compléter l'étude par des simulations pour des géométries de sections non carrées privilégiant la formation de structure en rouleaux de façon à évaluer la stabilité intrinsèque des cellules particulières obtenues lors des présentes simulations. Note remise le 22 janvier 1993, acceptée le 8 février 1993. RÉFÉRENCES BIBLIOGRAPHIQUES Ph. CARRIERE et D. JEANDEL, A 3D finite element method for the simulation of thermoconvective flows and its performances on a vector-parallel computer, Int. J. numer. methodsfluids, 12, 1991, p. 929-946. I. CATTON, The effect of insulating vertical walls on the onset of motion in a fluid heated from below, Int. J. Heat Mass Transfer, 15, 1972, p. 665-672. M. DUBOIS et P. BERGE, Experimental study of the velocity field in Rayleigh-Bénard convection, J. Fluid Mech., 85, 1978, p. 641-653. B. F. EDWARDS, Crossed rolls at onset of convection in a rigid box, J. Fluid Mech., 191, 1988, p. 583-597. P. LE GAL et V. CROQUETTE, Appearance of a square pattern in a Rayleigh-Bénard experiment, Phys. Fluids, 31, (11), 1988, p. 3440-3442. D. R. JENKINS et M. R. E. PROCTOR, The transition from roll to square-cell solutions in Rayleigh-Bénard convection, J. Fluid Mech., 139, 1984, p. 461-471. M. R. E. PROCTOR, Planform selection by finite-amplitude thermal convection between poorly conducting slabs, J. Fluid Mech., 113, 1981, p. 469-485. E. M. SPARROW, R. J. GOLDSTEIN et V. K. JONSSON, Thermal
p. 929-946. I. CATTON, The effect of insulating vertical walls on the onset of motion in a fluid heated from below, Int. J. Heat Mass Transfer, 15, 1972, p. 665-672. M. DUBOIS et P. BERGE, Experimental study of the velocity field in Rayleigh-Bénard convection, J. Fluid Mech., 85, 1978, p. 641-653. B. F. EDWARDS, Crossed rolls at onset of convection in a rigid box, J. Fluid Mech., 191, 1988, p. 583-597. P. LE GAL et V. CROQUETTE, Appearance of a square pattern in a Rayleigh-Bénard experiment, Phys. Fluids, 31, (11), 1988, p. 3440-3442. D. R. JENKINS et M. R. E. PROCTOR, The transition from roll to square-cell solutions in Rayleigh-Bénard convection, J. Fluid Mech., 139, 1984, p. 461-471. M. R. E. PROCTOR, Planform selection by finite-amplitude thermal convection between poorly conducting slabs, J. Fluid Mech., 113, 1981, p. 469-485. E. M. SPARROW, R. J. GOLDSTEIN et V. K. JONSSON, Thermal instability in a horizontal fluid layer: effect of boundary conditions and non-linear temperature profile, J. Fluid Mech., 18, 1964, p. 513-528. K. STORCK et U. MÜLLER, Convection in boxes: experiments, J. Fluid Mech., 54, 1972, p. 599-611. P. C. et D. J. : LMFA, URA-CNRS n° 263, 36, avenue G.-de-Collongues, 69130 Ecully, France; A. C. P. B. : Université de Brasilia, Dpt. Mech. Eng., Brasilia, Brazil. In the gaussian case, if W is not positive definite, we introduce an auxiliary norm on C in order to define a norm on Z0 z. Such an auxiliary norm is not unique, but the associated topology on z is not ambiguous. When C is a Hilbert space of functions on ℝ, and Q a positive quadratic form on ℝ, the pair (C', C) makes a quadruplet with the Schwartz spaces S(ℝ) and Z0 (ℝ) which generalizes the Gel'fand triplet (S(ℝ) L2(ℝ) S(ℝ)). We prefer the quadruplet (7) to the triplet (8) because we need an operator D : C → Z0 more general than a map from C into itself: for instance, in quantum physics D is a differential operator. We prefer the pair (C', D) to the pair (S(ℝ), S(ℝ)) because the