Introduction à la simulation du gel dans le domaine de la multiphysique COMSOL

Dans le domaine de l'ingénierie aérospatiale, le gel sur les surfaces des avions peut perturber le débit d'air et réduire l'ascenseur. Dans le domaine de la réfrigération, l'accumulation de gel sur les bobines d'évaporateur dégrade l'efficacité du transfert de chaleur et augmente la consommation d'énergie. Dans le domaine de la science climatique, l'accumulation de gel sur la végétation et le sol influence l'énergie locale et les bilans d'humidité. La simulation précise de ces processus nécessite un traitement couplé du transfert de chaleur, du transport de masse et du changement de phase.

Principes physiques de la formation et de la fusion du gel

Le gel se forme lorsque la vapeur d'eau dans l'air se sublime directement sur une surface solide dont la température est inférieure au point de gel. Le point de gel dépend de la pression et de la température de la vapeur locales; il est généralement inférieur au point de rosée parce que la pression de vapeur de saturation sur la glace est inférieure à celle de l'eau liquide. La fonte survient lorsque la température de surface s'élève au-dessus de 0 °C, ce qui provoque la transition de la glace vers l'eau liquide, qui peut alors s'évaporer ou s'écouler, selon les conditions de surface et de débit d'air.

Équations régissant la phase de fluidité

Dans le domaine du flux d'air, les équations de Navier-Stokes (RANS) moyennes de Reynolds sont généralement résolues pour modéliser le flux turbulent autour des surfaces de gel-prone. L'équation de continuité assure la conservation de la masse, tandis que les équations de moment intègrent la fermeture de turbulence, souvent à l'aide des modèles K-ε ou K---SST. L'équation d'énergie comprend le transfert convectif de chaleur et, lorsqu'elle est couplée au transport d'humidité, l'effet de la libération ou de l'absorption de chaleur latente. Pour le champ d'humidité, une équation convection-diffusion pour la fraction de masse de vapeur d'eau est résolue, avec un terme source ou évier représentant la sublimation ou le dépôt.

Modèle de changement de phase pour le gel

Dans le cadre du COMSOL, le changement de phase peut être mis en oeuvre par une fonction en échelon lisse qui définit la fraction de masse de glace comme étant fonction de la température et de la concentration de vapeur. La chaleur latente est incluse comme terme de source volumétrique dans l'équation d'énergie, proportionnelle à la dérivée temporelle de la fraction de masse de glace. Pour la croissance du gel, un sous-modèle suit la porosité et la conductivité thermique de la couche de gel, qui évoluent tous deux avec la densité et l'épaisseur.

Configuration du modèle COMSOL étape par étape

Création de géométrie et considérations de mesh

Un modèle typique comprend une plaque plate représentant une aile ou une aile d'échangeur de chaleur, placée dans un canal d'écoulement qui définit le domaine de l'air. Pour les géométries axisymétriques (par exemple un tube de refroidissement), une configuration axisymétrique 2D peut réduire le coût de calcul. Le maillage doit être affiné près de la surface de gelage pour résoudre les gradients aigus de température, de concentration de vapeur et de vitesse. Un maillage de la couche limite avec 5 à 10 couches de prisme et une première épaisseur cellulaire telle que y+ Φ 1 soit recommandé pour les modèles de turbulence qui résolvent la sous-couche visqueuse.

Définition des interfaces physiques

COMSOL fournit des interfaces de physique dédiées qui peuvent être combinées:

  • Transfert de chaleur dans les fluides (ht) – pour la distribution de la température dans l'air et dans le solide.
  • Transfert de chaleur dans les solides (ht) – pour la paroi ou le substrat.
  • Transport d'espèces diluées (tds) – pour la fraction de masse de vapeur d'eau dans l'air.
  • Flux de turbulence, k-ε ou k----- (spf) – pour le transport de l'élan.
  • Raccords multiphysiques[ – couplage de flux non isotherme (couplage de pseudo-couplage) et couplage personnalisé de dépôt de gel qui relie le flux de vapeur au bilan thermique par l'intermédiaire de la condition limite.

Pour incorporer la croissance du gel comme limite mobile, on peut utiliser la fonction Moving Mesh[ ou Deformed Mesh. Par ailleurs, pour plus de simplicité, on peut utiliser une méthode à grille fixe avec porosité efficace, où la couche de gel est traitée comme un milieu poreux avec une épaisseur et des propriétés thermiques dépendantes du temps.

Propriétés du matériau

Pour l'air, la densité suit la loi idéale sur les gaz, la conductivité thermique et la viscosité sont dépendantes de la température (utilisation des expressions intégrées). La diffusion de vapeur d'eau dans l'air peut être tirée de formules empiriques (p. ex., du manuel Cussler ou Perry). Pour la couche de gel, la densité du givre varie généralement de 50 à 300 kg/m3 selon les conditions de dépôt. Une expression telle que ρ frost = 650 − 2,0×106/(T frost+273,15) ou la corrélation de Hayashi (1977) peut être utilisée. La conductivité thermique du givre k frost (W/(m·K)) est souvent donnée par k frost = 0,024 + 7,52×10−4 ρ frost + 4,1×10−9 ρ frost3. La chaleur latente de la sous-limation est d'environ 2,83×106 J/kg, et pour la fusion elle est de 3,34×105 J/kg.

Limite et conditions initiales

Pour la simulation standard de gel:

  • Inlet: vitesse constante (1–5 m/s typique pour la réfrigération), température (p. ex., -5 à 10 °C) et humidité relative (50–90%).
  • Extrait: pression de jauge nulle ou état de sortie.
  • Surface de gel: paroi avec température prescrite (p. ex., -15°C pour une plaque froide) ou couplée à un bilan énergétique qui comprend convection et chaleur latente. La concentration de vapeur à la paroi est réglée à la valeur de saturation sur la glace à cette température de paroi (par une fonction ou une table de recherche).
  • Conditions d'entrée:[ domaine d'écoulement entier réglé à la température et à l'humidité d'entrée, fraction de masse de glace zéro. La simulation est fonction du temps; les durées d'exécution typiques s'étendent de minutes à heures en temps physique, avec des temps de 1 à 10 secondes.

Simulation de la formation de gel: flux de travail numérique

Une fois le modèle construit, exécutez une étude dépendante du temps. La séquence du résolveur COMSOL=S sépare automatiquement la physique — premier flux de résolution (stable ou transitoire) puis les équations scalaires — mais pour les cas fortement couplés, une approche entièrement couplée peut être plus stable. Surveillez les sorties suivantes:

  • Épaisseur du gel en fonction du temps (après traitement de l'interface mobile ou du champ de fraction de masse de glace).
  • Coefficient de transfert de chaleur de surface (calculé à partir du flux total de chaleur à travers la paroi divisé par la différence de température).
  • Masse totale de gel par unité de surface.
  • Profils de température et de vitesse près de la paroi.

Les principales conclusions d'un modèle de formation de gel comprennent la période de croissance rapide initiale (dominée par le transfert de masse), suivie d'une phase de croissance plus lente lorsque la couche de gel isole la surface. Le modèle peut également prédire la distribution spatiale du gel, souvent plus épaisse près du bord d'attaque d'une plaque plate en raison de taux de diffusion de vapeur plus élevés.

Modélisation des cycles de fusion et de dégivrage du gel

Dans le cadre de la COMSOL, une approche consiste à inclure un modèle de fusion à couche mince qui active lorsque la température de surface dépasse 0°C. Le front de fusion peut être suivi par une méthode de mise en niveau ou de suivi d'interface, mais pour de nombreuses applications techniques, une méthode intégrale simplifiée de bilan thermique suffit. Impose un flux thermique ou une condition de limite de température élevée sur la paroi (p. ex., 30°C pour un chauffage électrique au dégivrage) et résout pour le temps nécessaire pour fondre complètement une épaisseur de gel donnée. La simulation doit tenir compte de l'énergie consommée par la chaleur latente de fusion ainsi que du chauffage raisonnable du gel poreux et du substrat.

L'optimisation du cycle du dégivrage peut être effectuée en paramétrant la puissance thermique, la durée du chauffage et la période de drainage, et en utilisant le module d'optimisation de COMSOL, afin de minimiser l'utilisation de l'énergie tout en assurant l'élimination complète du gel.

Validation et pratiques exemplaires

La validation des modèles de gel par rapport aux données expérimentales est cruciale. Les points de repère publiés comprennent les travaux de Liu et Jacobi (2013) sur la croissance du gel sur des plaques plates en convection forcée et les études de Padhmanaban et Sherif (2000) sur les propriétés du gel.

  • Épaisseur du gel par rapport au temps
  • Température de surface du gel
  • Réduction globale du coefficient de transfert de chaleur due au gel

Les utilisateurs doivent étalonner les corrélations empiriques de conductivité thermique et de densité du gel à leurs plages de température et d'humidité spécifiques. Une analyse de sensibilité sur ces paramètres aide à quantifier l'incertitude. Les études de raffinement des mailles devraient cibler une incertitude d'épaisseur du gel de <5%. Additionally, for turbulent flows, ensure that the y+ value at the first cell is checked—the built‑in wall functions in the k-ε model may require the first cell to lie in the log‑law region (y+ > 30), alors qu'un modèle de faible turbulence (comme k----) nécessite y+ ~1.

Considérations avancées : Turbulence, couplage et demande numérique

L'interaction entre le gel et la turbulence est double : la rugosité du gel modifie les structures turbulentes proches de la paroi et les tourbillons turbulents améliorent le transport de vapeur. Certains modèles avancés intègrent une fonction de rugosité dans le traitement de la paroi de turbulence, augmentant la rugosité de surface à mesure que le gel s'accumule.

Pour les longs temps physiques (p. ex., heures de croissance du gel), il est recommandé de suivre le temps d'adaptation avec une tolérance de 0,001. Les utilisateurs peuvent profiter de balayages paramétriques efficaces de COMSOL pour évaluer les multiples conditions d'entrée. Pour les grands modèles 3D, un maillage grossier sur les régions non friables et un maillage fin seulement à proximité de la surface peuvent réduire le nombre de cellules. Le temps total de simulation pour un modèle 2D (plusieurs milliers d'éléments) avec 30 minutes de temps physique peut prendre 1 à 4 heures sur un poste de travail, tandis qu'un modèle 3D d'un échangeur de chaleur fin et tube peut prendre 24 heures ou plus.

Applications et pertinence industrielle

Systèmes anti-circulation aérospatiales

Les modèles COMSOL aident à concevoir des systèmes d'antigivrage électrothermique en prédisant la puissance nécessaire pour maintenir les surfaces au-dessus du point de gel. Des études paramétriques sur la disposition du chauffage et le cycle de travail peuvent être effectuées dans le même modèle, y compris l'effet de la fonte du gel et du ruissellement de l'eau.

Performances de réfrigération et de thermopompe

Dans les bobines d'évaporateur de pompes à chaleur à source d'air, l'accumulation de gel réduit le débit d'air et le transfert de chaleur, ce qui entraîne une dégradation de la performance. Les ingénieurs utilisent COMSOL pour déterminer les critères optimaux d'initiation au dégivrage en fonction de l'épaisseur de gel mesurée ou de la chute de pression.

Systèmes de construction économes en énergie

Les simulations aident à concevoir des stratégies de contournement ou de préchauffage qui réduisent l'accumulation de gel tout en maintenant le confort. Les modèles d'énergie de construction couplés peuvent être simplifiés en utilisant la co-simulation COMSOL avec des outils comme EnergyPlus (via interface externe).

Protection contre les glaces éoliennes

Les modèles multiphysiques qui comprennent la rotation des pales (avec des interfaces de machines rotatives) et les charges de vent stochastiques peuvent simuler l'accumulation de gel dans des conditions atmosphériques variables. Ces modèles soutiennent le développement de systèmes de dégivrage actifs, tels que le chauffage résistif ou les vibrations ultrasoniques, validés dans l'environnement COMSOL.

Défis et nouvelles orientations de la recherche

Malgré la puissance de COMSOL, la modélisation du gel reste difficile en raison de la complexité de la nucléation microéchelle et de la nature stochastique de la croissance du cristal de gel. De nombreux modèles actuels reposent sur des corrélations empiriques qui ne sont validées que dans des gammes de paramètres étroites. La recherche future se dirige vers la simulation numérique directe (DNS) des interactions de gouttelettes de congélation et de la formation de réseaux cristallins, bien qu'elles soient calculables coûteuses pour les problèmes macroéchelle.

Conclusion

La validation des données expérimentales demeure une étape critique pour assurer la crédibilité de la simulation. Grâce aux progrès continus des ressources informatiques et des sous-modèles physiques, COMSOL restera un outil clé pour concevoir des systèmes résistant au gel et optimiser les stratégies de dégivrage dans les secteurs de l'aérospatiale, de la réfrigération et de l'énergie. Pour plus de détails, consulter les bibliothèques d'application COMSOL pour des exemples sur le changement de phase et le transport d'humidité, et examiner les articles par COMSOL sur la modélisation du gel dans les échangeurs de chaleur et sur la simulation du dégivrage d'aéronef.