Modèle thermique · données · limites

Méthodo­logie

Cette page décrit exactement comment NADIR calcule — les données, les équations, les seuils — et ce que ça vaut : où le modèle est solide, où il se trompe, et dans quel sens il se trompe. Aucun secret : tous les paramètres sont listés ici, tels qu'ils sont dans le code de la version 1.0.

La question posée

Ouvrir les fenêtres rafraîchit si l'air dehors est plus frais que chez vous. Mais deux choses compliquent l'affaire. D'abord, ouvrir suppose en général de relever les volets : le soleil qui entre peut coûter plus que l'air frais ne rapporte. Ensuite, la fraîcheur gagnée doit se stocker quelque part — dans les murs — sinon elle repart en une heure.

C'est donc un problème de timing et d'arbitrage, pas une règle fixe. NADIR le traite de la seule façon honnête : en simulant votre logement heure par heure sur les 30 prochaines heures, et en cherchant les créneaux où ouvrir fait vraiment gagner.

Les données d'entrée

Et ce que NADIR n'utilise pas, autant le dire : pas de données de rayonnement solaire mesuré, pas de vent, pas de nébulosité. Le soleil est un modèle géométrique de ciel clair (détaillé plus bas) — c'est une vraie limite, examinée dans la section sur les biais.

Tant qu'aucun lieu réel n'est choisi, l'app affiche une sinusoïde d'exemple étiquetée « Exemple » — et les raccourcis Siri refusent de répondre dessus plutôt que d'inventer un verdict.

Le modèle : deux nœuds, l'air et la masse

L'observation de départ : quand vous aérez, l'air se rafraîchit en minutes ; les murs mettent des heures. Un seul « réservoir de chaleur » ne peut pas représenter ça. Deux suffisent : un nœud air (l'air de la pièce et le mobilier léger) et un nœud masse (murs, planchers, cloisons), couplés entre eux et vers l'extérieur.

DEHORS T_ext (météo) AIR T_a · C_a=150 kJ/K MASSE T_m · C_m=2,2–5,2 MJ/K H_env 50 W/K + H_vent si ouvert H_am 300 W/K soleil 0–650 W/façade internes 150 W volets fermés : 30 % du soleil passe · fenêtres ouvertes (volets relevés) : 100 %
Le modèle complet tient dans ce schéma. La ventilation agit sur l'air ; l'air, lentement, sur la masse.

À chaque heure de météo, deux bilans d'énergie sont intégrés (10 sous-pas de 6 minutes, méthode d'Euler explicite — dimensionnés pour rester sous la limite de stabilité du cas le plus ventilé) :

Ca·dTa/dt = Hext·(Text − Ta) + Ham·(Tm − Ta) + Φ(t)
Cm·dTm/dt = Ham·(Ta − Tm)
avec Hext = Henv (+ Hvent fenêtres ouvertes) · Φ = internes + soleil admis

Tous les paramètres, sans exception :

ParamètreValeurD'où ça vient
Pièce type20 m² · 50 m³ Pièce de vie type ; voir « biais » pour ce que ça implique
Capacité de l'air Ca150 kJ/K Air seul ≈ 60 kJ/K (ρcp = 1 206 J/m³K × 50 m³) + mobilier léger
Capacité de la masse Cm = κ·A2,2 / 3,3 / 5,2 MJ/K κ = 110 / 165 / 260 kJ/(m²K) selon murs légers / moyens / lourds — classes de la norme EN ISO 13790
Couplage air ↔ masse Ham300 W/K Échange convectif/radiatif sur ~40 m² de parois, ≈ 7,7 W/m²K
Enveloppe Henv50 W/K (×1,15 traversant) Pertes parois + infiltrations, logement ancien type
Ventilation n (fenêtres ouvertes)3 / 8 / 15 vol/h Renouvellements d'air par heure : entrebâillé / bien ouvert / grand courant d'air
Facteur traversant×1,5 ou ×1,8 Deux façades adjacentes / opposées (le courant traverse)
Hvent = ρcp·n·V·f / 360050 → 452 W/K Selon vos réglages ; à comparer aux 50 W/K de l'enveloppe
Apports internes150 W Occupants + appareils, constant
Soleil direct (pic, par façade)650 W ~2 m² de vitrage en plein soleil, facteur solaire compris
Soleil diffus (pic)120 W Ciel, toutes orientations
Volets fermés30 % passe Volets + rideaux imparfaits
Hystérésis d'ouverture−0,75 / −0,25 °C Voir « la décision d'ouvrir »
Pas de calcul1 h ÷ 10 Sous-pas de 6 min, Euler explicite

Ces valeurs donnent au logement simulé trois échelles de temps, et c'est elles qui portent tout le comportement :

≈ 3–5 minl'air, fenêtres grandes ouvertes : τ = Ca/Htotal. L'air « suit » le dehors presque instantanément — conforme au vécu.
≈ 2–5 hla masse : τ = Cm/Ham. C'est pour ça que 10 minutes d'aération rafraîchissent l'air mais pas le logement — et qu'une nuit entière, oui.
≈ 13–30 hle bâtiment fermé : τ = (Ca+Cm)/Henv. Fermé, le logement dérive lentement vers la température moyenne extérieure.

Ces ordres de grandeur correspondent aux constantes de temps mesurées sur de vrais bâtiments (typiquement 15 à 50 h pour un logement fermé). C'est le premier argument de confiance : la dynamique est la bonne, même si les valeurs absolues restent approximatives.

L'état de départ

La courbe démarre exactement à la température que vous saisissez : l'air y est fixé, et la température de la masse est déduite en résolvant l'équilibre de l'air à l'instant initial, fenêtres fermées (dTa/dt = 0) :

Tm = Ta − [Henv·(Text − Ta) + Φ0] / Ham

Concrètement : s'il fait chaud dehors et que le soleil tape, la masse est supposée un peu plus fraîche que l'air (elle absorbe) ; la nuit après une journée chaude, un peu plus chaude (elle restitue). Pas de saut artificiel au départ de la courbe, et votre thermomètre calibre réellement le modèle.

Le soleil, sans données de rayonnement

Faute d'utiliser un rayonnement mesuré, chaque façade reçoit une cloche de ciel clair en cos², avec un pic direct de 650 W par façade exposée :

FaçadePicFenêtre de soleil direct
Est09 h 0005 h → 13 h
Sud13 h 0008 h 30 → 17 h 30
Ouest17 h 3013 h → 22 h
Norddiffus uniquement

S'y ajoute un diffus de 120 W au pic (cloche large centrée à 13 h 30). En logement traversant, le vitrage se répartit (×0,6 par façade) et le diffus augmente un peu (×1,2 — logement d'angle). Volets fermés, 30 % du soleil passe ; fenêtres ouvertes — donc volets relevés — tout passe. Ce dernier point est le cœur de l'arbitrage jour/nuit.

Autocritique immédiate : c'est un soleil de ciel clair, tous les jours, toute l'année. Pas de nuages, pas de saison, pas de latitude, pas de masques (arbre, immeuble en face). C'est le morceau le plus schématique du modèle. Sa circonstance atténuante : l'essentiel du verdict se joue la nuit, où le soleil est nul et où cette approximation ne pèse plus rien.

La décision d'ouvrir

À chaque heure, la règle compare l'extérieur à l'air intérieur simulé, avec une hystérésis et une pénalité solaire :

ouvrir  si  Text < Ta − 0,75 °C − P
fermer  si  Text > Ta − 0,25 °C − P
P = 0,7 · Φsol / Hvent  (le soleil admis en plus volets relevés, converti en °C équivalents)

L'hystérésis (ouvrir sous −0,75, ne refermer qu'au-dessus de −0,25) évite un verdict qui clignoterait à chaque frôlement des courbes. La pénalité P dit simplement : ouvrir ne vaut le coup que si le refroidissement par l'air dépasse le surplus de soleil qui entrera volets relevés.

Deux exemples chiffrés, tirés du modèle tel quel :

Ce que veulent dire les chiffres affichés

Les widgets et Siri recalculent exactement la même simulation, avec les mêmes données partagées — jamais une version simplifiée.

Là où le modèle est solide

Les biais, honnêtement

Chaque simplification pousse le résultat dans un sens. Les voici, du plus important au plus anecdotique — optimiste = l'app promet plus de fraîcheur que la réalité, prudent = elle en promet moins.

La bonne lecture découle de tout ça : le créneau à ±1 h près, les tendances, les comparaisons — pas la décimale. Le chiffre affiché avec une décimale est une sortie de calcul, pas une promesse de précision.

Se calibrer chez soi

Le meilleur correctif ne coûte rien : un thermomètre. La température intérieure que vous saisissez est l'état initial exact du modèle — c'est le levier de précision numéro un, bien devant tous les réglages.

Sources

Les capacités thermiques κ viennent des classes de la norme EN ISO 13790. La démarche (ventilation nocturne, masse thermique, arbitrage solaire) s'appuie sur la littérature publique — MIT CoolVent, OMS/NCBI NBK143285, Kolokotroni et al. 1998, Yam et al. 2007, Building & Environment 2021, AIVC, entre autres. La liste complète est sur la page copyright. Les données météo proviennent d'Open-Meteo (CC BY 4.0).

Une erreur dans cette page, un biais oublié, un désaccord sur un ordre de grandeur ? Écrivez-moi : solal.gendrin@gmail.com.