Introduction à la modélisation de la chauffe-eau solaire

Les chauffe-eau solaires transforment le rayonnement solaire en énergie thermique pour fournir de l'eau chaude pour les applications résidentielles, commerciales et industrielles. Ces systèmes offrent une alternative durable aux chauffe-eau électriques ou au gaz classiques, réduisant à la fois les coûts énergétiques et les émissions de carbone. Cependant, les performances thermiques d'un chauffe-eau solaire dépendent d'un jeu complexe de mécanismes de transfert de chaleur – la conversion par des composants solides, la convection dans le fluide de travail et le rayonnement de la surface absorbante vers les environs.

Pour concevoir des chauffe-eau solaires à haute efficacité, les ingénieurs doivent prédire avec précision la distribution de la température, les vitesses d'écoulement et les taux de perte de chaleur dans des conditions de fonctionnement variables.Computational Fluid Dynamics (CFD) utilisant Ansys Fluent fournit un outil rigoureux pour simuler ces physique couplée.En créant un prototype virtuel du collecteur et de ses tuyauteries, les analystes peuvent explorer les modifications de conception, optimiser la géométrie et évaluer les performances avant de construire des prototypes physiques.

Les fondamentaux du transfert de chaleur dans les chauffe-eau solaires

Les chauffe-eau solaires sont généralement constitués d'une surface absorbante (souvent revêtue d'un matériau sélectif pour maximiser l'absorption solaire et minimiser l'émission thermique), d'un couvercle transparent (verre ou polymère), d'une isolation à l'arrière et aux côtés et d'un réseau de transport de fluides.

  • Conduction: La chaleur traverse la plaque d'absorbeur solide et les parois des tubes. La conductivité thermique des matériaux (cuivre, aluminium ou acier inoxydable) détermine la résistance à la conduction.
  • Convection: Le fluide de travail (eau ou mélange glycol-eau) absorbe la chaleur des parois du tube. Dans les systèmes de circulation forcée, une pompe entraîne le débit du fluide; dans les systèmes de circulation naturelle (thermosipon), le débit est déterminé par des différences de densité.
  • Radiation: L'absorbeur émet un rayonnement infrarouge à longue ondes vers le couvercle, et le couvercle re-radiage la partie arrière. Le rayonnement solaire passe par le couvercle et est absorbé par l'absorbeur. Le transfert de chaleur radiative est modélisé à l'aide de modèles surface-surface (S2S) ou d'ordonnées discrètes (DO) dans Fluent.

De plus, la chaleur est perdue du collecteur à l'environnement par le couvercle supérieur (convection plus rayonnement) et par l'isolation du dos et des côtés (convection et convection).Un modèle CFD détaillé doit tenir compte de tous ces chemins pour prédire avec précision le gain de chaleur utile et l'élévation de la température de l'eau qui en résulte.

Types de chauffe-eau solaires modélisés

L'approche de modélisation d'Ansys Fluent s'adapte à la configuration spécifique du chauffe-eau solaire. Les types les plus courants sont les suivants:

Collecteurs de plaques plates

Il s'agit d'une plaque d'absorbeur plate avec tubes de riser fixés ou d'un tube serpentin. Le fluide de travail coule à travers les tubes, extrayant la chaleur. Les collecteurs de plaques plates sont largement utilisés dans les applications résidentielles et sont relativement simples à modéliser: une simulation de transfert de chaleur conjugué 2D ou 3D avec un flux laminaire ou turbulent dans les tubes.

Collecteurs de tubes évacués

Les collecteurs de tubes évacués utilisent un vide entre l'absorbeur et le tube de verre extérieur pour réduire de façon significative les pertes de chaleur convectif. L'absorbeur peut être un tube à flux direct (fluide passant par un tube en U ou un tube concentrique) ou un type de tube de chaleur (un changement de phase se produit à l'intérieur d'un tube scellé).

Systèmes thermosiphoniques

Dans les chauffe-eau solaires thermosiphon, le réservoir de stockage est placé au-dessus du collecteur. La circulation naturelle entraînée par la flottabilité doit être capturée avec précision. Cela nécessite une simulation transitoire avec couplage entre le collecteur et l'hydrodynamique du réservoir.

Mise en place d'une simulation de fluide Ansys pour les chauffe-eau solaires

Pour construire une simulation crédible, il faut un flux de travail systématique. Les sous-sections suivantes détaillent chaque étape.

Création et simplification de la géométrie

Pour les capteurs à plaques plates, la géométrie comprend la couverture en verre, l'espace d'air, la plaque d'absorption, les tubes et l'isolation. Simplifiez si possible : retirez les petits filets, les boulons et les détails non essentiels. Pour les tubes périodiques, un seul tube représentatif avec des conditions de limite périodiques peut réduire le coût de calcul. Dans les capteurs de tubes évacués, modélisez l'absorbeur, le tube en verre et le canal fluide, mais traitez le vide comme une région solide avec une très faible conductivité thermique (p. ex. 0,005 W/m·K) ou utilisez la région vide avec un rayonnement seulement.

Génération de Mesh

Un maillage de haute qualité est crucial pour une DFC précise. Utilisez Ansys Meshing ou Fluent Meshing pour créer un maillage hybride :

  • Zones fluides: Utiliser des éléments hexaédriques ou polyédriques avec des couches de prisme près des parois pour capturer les couches de limite thermique et de vitesse. La première hauteur de cellule doit être telle que y+ ~1 pour les modèles de turbulences SST k----- ou bas-Re, ou y+ ~30 pour les fonctions de paroi avec k-ε.
  • Zones solides: Les éléments tétraédriques ou hexaédriques sont acceptables; assurer au moins quelques éléments à travers l'épaisseur des parois minces (plaques absorbantes, parois tubulaires) pour résoudre la conduction.
  • Indépendance de la maille:[ Effectuer une étude de convergence de grille en raffinant le maillage jusqu'à ce que les sorties clés (température de sortie, chute de pression) changent de moins de 1 à 2 %.

Configuration physique

Dans Ansys Fluent, activez l'équation d'énergie et choisissez un modèle de turbulence (si le débit est turbulent).Pour la plupart des débits de capteurs solaires, le nombre de Reynolds est faible (laminaire ou transitoire). Utilisez le modèle Transition SST pour une prévision de transition précise, ou k----SST pour un débit complètement turbulent.

Définir les propriétés du matériau:

  • Fluide de travail : mélange eau-eau, glycol-eau. Utiliser la densité dépendante de la température, la chaleur spécifique, la viscosité et la conductivité thermique (p. ex., à partir de courbes NIST REFPROP ou polynômes).
  • Matériaux solides: cuivre ou aluminium pour absorbeur, verre pour couverture, mousse de polyuréthane pour isolation.

Régler les conditions limites :

  • Inlet: Débit massique ou entrée de vitesse à température fixe (par exemple, 20°C).
  • Outlet: Sortie de pression, généralement à pression manométrique atmosphérique.
  • Flux de chaleur solaire: Appliquer un flux de chaleur uniforme ou variable dans le temps (p. ex. 1000 W/m2) sur la surface extérieure de l'absorbeur.
  • Convection et rayonnement externes:[ Utiliser une condition thermique mixte sur le revêtement de verre et l'isolation arrière. Spécifier la température ambiante (par exemple, 25°C) et le coefficient de transfert de chaleur (par exemple, 10 W/m2·K pour la convection naturelle, 30 W/m2·K pour les conditions venteuses).

Pour le transfert radiatif de chaleur au sein du collecteur, le modèle S2S fonctionne bien lorsque le support participant n'est pas participant (les vides d'air sont transparents). Si le fluide absorbe/émets les radiations (par exemple, dans les capteurs d'absorption directe), utiliser le modèle DO avec le coefficient d'absorption.

Paramètres du solvant et stratégie de solution

Pour les problèmes de circulation naturelle, il peut être nécessaire de réduire les facteurs de sous-relax pour l'impulsion et la densité (par exemple 0,3–0,5) pour stabiliser la convergence. Utilisez la discrétisation du vent ascendant de deuxième ordre pour l'impulsion et l'énergie. Initialiser avec un champ de température et d'écoulement uniforme, puis effectuer la simulation. Surveiller les résidus et suivre les quantités intégrées (température de sortie, taux de transfert de chaleur). La convergence est réalisée lorsque les résidus tombent sous 1e-5 pour la continuité et l'impulsion, 1e-6 pour l'énergie, et la température de sortie se stabilise dans un délai de 0,1 K sur 100 itérations.

Des simulations transitoires sont nécessaires pour les systèmes à rayonnement solaire variable (p. ex. profils quotidiens) ou dynamique thermosiphonienne. Régler une taille d'étape qui résout les échelles de temps et de débit thermique (généralement 0,1-1 s pour la réponse thermique du collecteur, 1-10 s pour la circulation naturelle).

Traitement et interprétation des résultats

Après la simulation, extraire et analyser les paramètres suivants:

  • Distribution de température:[ Les placettes de contour sur la plaque d'absorbeur révèlent des points chauds et un chauffage inégal.
  • Modèles de refoulement:[ Les vecteurs de vélocité ou les rationalisations dans les tubes montrent si le débit est complètement développé, recirculation ou contournant les zones chaudes.
  • Analyse du flux de chaleur:[ Calculer le gain de chaleur utile (intégral du flux de chaleur de paroi sur les surfaces du tube en contact avec le fluide) et comparer avec l'énergie solaire incidente pour calculer l'efficacité du collecteur.
  • Drop de pression:[ Évaluer la puissance de pompage requise.

Utilisez les outils de rapport d'Ansys Fluent pour calculer la température moyenne pondérée par zone à la sortie et la température moyenne par volume de l'absorbeur. Exportez les données pour un traitement ultérieur dans MATLAB ou Excel.

Validation et vérification

Pour faire confiance au modèle, valider par rapport aux données expérimentales de la littérature ou des essais en laboratoire. Comparer la température de sortie, la courbe d'efficacité (η vs. (T in – T amb)/G) et la chute de pression dans des conditions d'essai standard (par exemple, ISO 9806). Ajuster les paramètres de simulation tels que le coefficient de transfert de chaleur sur le couvercle de verre, l'émissivité de l'absorbeur ou les propriétés du fluide pour correspondre aux mesures dans un délai de 5 à 10 %.

Vérifier l'indépendance du maillage et la convergence temporelle. Documenter toutes les hypothèses (par exemple, irradiation solaire uniforme, absence de salissure, isolation idéale).

Exemple pratique : Diamètre optimal du tube de riser

Considérez un collecteur de plaque plate avec dix tubes de riser parallèles. Utilisez Ansys Fluent pour étudier l'effet du diamètre intérieur du tube sur l'efficacité thermique. Gardez la constante du flux solaire (800 W/m2), la température d'entrée 30°C, la température ambiante 25°C, le vent 2 m/s. Diamètre de la varie de 8 mm à 16 mm par paliers. Pour chaque cas, faites une simulation en état d'équilibre. Les résultats montrent que les diamètres plus petits donnent une vitesse de transfert de chaleur plus élevée et convectif, mais aussi une chute de pression plus élevée et un débit de masse plus faible pour une pompe donnée (ou circulation naturelle).

Défis de la modélisation des chauffe-eau solaires

Plusieurs facteurs compliquent la modélisation précise des CFD :

  • L'irradiation solaire non uniforme:[ L'ombrage du monde réel, le couvert nuageux et les variations de l'angle d'incidence nécessitent un profil d'irradiation transitoire.
  • Le débit en deux phases dans les collecteurs de tuyaux thermiques:[ La simulation de l'évaporation et de la condensation est coûteuse en calcul; une approche équivalente de source de chaleur peut être utilisée comme approximation.
  • Instabilité à la circulation naturelle: Le débit entraîné par la flottabilité peut osciller; un résolveur transitoire avec de petites étapes de temps est nécessaire.
  • Effets d'échelle:[ Modéliser un champ entier de collecteurs (ferme solaire) nécessite des simplifications comme des modèles de supports poreux ou de paramètres à boulets.

Malgré ces défis, Ansys Fluent fournit des capacités avancées – comme le modèle VOF pour les surfaces libres dans les systèmes à boucle ouverte ou le modèle multiphase eulérien pour le changement de phase – pour les résoudre.

Meilleures pratiques pour une simulation efficace

  • Commencez par un modèle axisymétrique 2D ou planaire 2D pour les études paramétriques, puis validez avec un modèle 3D pour la conception finale.
  • Utilisez des conditions de limites symétriques ou périodiques pour réduire la taille du modèle.
  • Employez le résolveur couplé pour des flux de convection naturelle fortement couplés.
  • Tirez parti de l'optimisation paramétrique d'Ansys Workbench pour balayer automatiquement les variables de conception.
  • Documenter la qualité du maillage (sauvage < 0.9, orthogonal quality > 0,1).

Tendances futures et méthodologies avancées

Le champ se dirige vers l'optimisation multiphysique en intégrant CFD avec des simulations structurales et optiques. La capacité d'Ansys Fluent à s'associer avec Ansys Mechanical[ pour l'analyse de la contrainte thermique ou avec Ansys Lumerical[ pour la modélisation optique détaillée deviendra plus courante. De plus, les modèles de substitution basés sur l'apprentissage de la machine formés sur les données Fluent peuvent accélérer le contrôle en temps réel des systèmes solaires de chauffage de l'eau.

Pour les ingénieurs et les chercheurs engagés à faire progresser la technologie solaire thermique, maîtriser CFD à Ansys Fluent est une compétence inestimable. Non seulement il réduit le coût et le temps du prototypage, mais permet également des conceptions innovantes qui repoussent les limites de l'efficacité.

Ressources extérieures

Conclusion

En modélisant le transfert de chaleur et le flux de fluides dans les chauffe-eau solaires avec Ansys Fluent, les ingénieurs peuvent optimiser la conception, prévoir les performances et réduire les risques de développement. En suivant un workflow de simulation structuré, de la simplification de la géométrie et de la génération de mailles par la configuration physique et le post-traitement, les analyseurs peuvent obtenir des informations quantitatives fiables.