Résistance thermiqueR
—m²·K/W
Par unité de surface
CFDapp
Yakhchal · pré-dimensionnement thermique
Ingénierie préliminaire · mur en terre
Comprendre les transferts
avant la CFD.
Estimez la conduction à travers la paroi et son effet sur le stockage de glace. Ajustez les paramètres pour explorer rapidement les ordres de grandeur d’un dôme simplifié.
Modèle stationnaire · paroi homogène
Résultats en direct
02 / SORTIESFlux surfaciqueQ
—W/m²
Flux signé Ti − Te
Surface latéraleS
—m²
Cône · base exclue
Puissance traversanteQ × S
—W
Valeur absolue · régime établi
Volume équivalent si le stock est entièrement en glace : — m³ pour — kg (ρglace = 917 kg/m³).
Sensibilité à l’épaisseur
Flux signé Q = λ/e · (Ti − Te), toutes choses égales par ailleurs.
Q(e)
Projection de l’état du stock
Bilan énergétique simplifié sur la durée choisie.
État après—
—
—
Fraction transformée en eau0 %
— kg/h estimés
Durée jusqu’à l’eau complète
—
Puissance nette reçue
—
La projection suppose une puissance constante et uniformément transmise au stock. Elle distingue le réchauffement de la glace, la fusion latente et le réchauffement éventuel de l’eau ; ce n’est pas une simulation transitoire CFD.
Préparation du cas CFD · OpenFOAM
Périmètre cible : air confiné + mur en terre + bassin cylindrique avec changement de phase. La simulation OpenFOAM se lance hors de ce site : exportez les paramètres d’entrée pour préparer le cas. Le calcul navigateur est un modèle réduit distinct d’une simulation OpenFOAM ; la projection de phase est un bilan énergétique indicatif, pas un solveur de fusion CFD.
Prêt · lancer le calcul
Champ de température dans la paroi
Maillage —— °C— °C
Convergence de l’itération
échelle logarithmiqueInitialisationItérationsRésidu —
—
—
—
—
—
—
—
—
Le bouton « Modèle web réduit » résout uniquement la conduction dans la paroi avec des corrélations de convection et de rayonnement ; Ra, Nu et Re* sont des estimations. Il ne résout ni Navier–Stokes ni un front de phase. « Exporter paramètres CFD » crée une fiche d’entrée JSON, pas un cas OpenFOAM maillé ou prêt à calculer. Pour exécuter le cas cible, il faudra produire et vérifier le maillage multi-régions, puis configurer et valider le solveur transitoire avec changement de phase.
Rappel méthodologique
Le tableau de référence utilise une conduction 1D stationnaire ; le solveur numérique ajoute une conduction 2D dans la paroi et des conditions convectives/radiatives estimées. La géométrie extérieure est un cône idéal et seule sa surface latérale échange. La projection de phase suppose une puissance constante et distingue chaleur sensible et chaleur latente. Infiltration, ponts thermiques, sol, pertes du bassin et inertie ne sont pas inclus. λ doit être mesurée sur la terre réelle du mur.
R = e / λ
U = 1 / R = λ / e
Q = U · (Ti − Te)
Vglace = M / ρglace
Slatérale = π · d/2 · √[(d/2)² + h²]
Ptraversante = |Q| · S
Ra = g β ΔT H³ / (ν α)
hi = Nu · kair / H
qrad = εσ (Tciel⁴ − Ts⁴)
Lfusion = 334 kJ/kg
Perspectives pour Star-CCM+
Faire évoluer le modèle simplifié à partir des mesures terrain et des essais matériaux. Les hypothèses et régimes devront être vérifiés avant de choisir les modèles numériques.
Boussinesq · convection naturelleRa · RayleighRe · ReynoldsNu · NusseltChangement de phase glace ↔ eauCompressibilité à évaluer
La projection de phase prend en compte le réchauffement sensible, L = 334 kJ/kg et le réchauffement de l’eau après fusion. La puissance nette analytique ou du solveur est supposée constante sur la durée choisie : c’est une borne haute indicative, pas une simulation transitoire avec front de phase.