@open3cl/engine 1.8.24 → 1.8.26

package/13.2_generateur_combustion_ch.js

@@ -272,6 +272,27 @@ function Tch_xfinal(x, Cdimref) {
   return Math.min(1, x / Cdimref);
 }
 
+/**
+ * Pertes à l'arrêt QP0 (kW, sur PCS) d'un générateur à combustion.
+ *
+ * Pour les générateurs d'air chaud, la valeur forfaitaire de la table 13.2
+ * (Pn × (1,75 − 0,55 × log Pn)) est exprimée en % de Pn : elle est ramenée en W (÷ 100).
+ * Cette conversion ne doit pas modifier di.qp0 : QPx() étant appelée pour chaque taux de charge,
+ * une division en place ramenait QP0 à ~0 (÷ 100 à chaque appel), et le terme 0,45 × QP0 du
+ * rendement de génération en était absent.
+ *
+ * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.1
+ * @param de {Donnee_entree}
+ * @param di {Donnee_intermediaire}
+ * @param k {number} rapport PCI/PCS de l'énergie
+ * @returns {number}
+ */
+export function qp0PcsKw(de, di, k) {
+  const type_gen_ch = enums.type_generateur_ch[de.enum_type_generateur_ch_id] || '';
+  const qp0 = type_gen_ch.includes('générateur à air chaud') ? di.qp0 / 100 : di.qp0;
+  return (qp0 * k) / 1000;
+}
+
 function QPx(x, de, di) {
   const type_gen_ch = enums.type_generateur_ch[de.enum_type_generateur_ch_id];
   const type_energie = enums.type_energie[de.enum_type_energie_id];
@@ -318,9 +339,7 @@ function QPx(x, de, di) {
       else QPx = QP30 + ((QP100 - QP30) * (x - 0.3)) / 0.7;
     }
   } else if (type_gen_ch.includes('générateur à air chaud')) {
-    // Change QP0 to be in decimal
-    di.qp0 = di.qp0 / 100;
-    qp0 = (di.qp0 * k) / 1000;
+    qp0 = qp0PcsKw(de, di, k);
 
     const QP50 = (0.5 * pn * (100 - rpint)) / rpint;
     const QP100 = (pn * (100 - rpn)) / rpn;
@@ -376,10 +395,11 @@ export function calc_generateur_combustion_ch(dpe, di, de, du) {
   const k = K[type_energie];
 
   // Pveil and QP0 are in kW
-  const Pveil = di.pveil / 1000;
-  // rg est d'abord calculé sur PCS (puis × k pour revenir en PCI) : QP0 doit donc être exprimé
-  // sur PCS, comme dans QPx() (qp0 × k). Cf. issue #205.
-  const QP0 = (di.qp0 * k) / 1000;
+  // rg est d'abord calculé sur PCS (puis × k pour revenir en PCI) : QP0 et Pveil (pertes de la
+  // veilleuse, consommation de combustible) doivent donc être exprimés sur PCS, comme dans QPx()
+  // (qp0 × k). Cf. issue #205 et autotests État MI4-0-x (Tribu : Pveil × k).
+  const Pveil = (di.pveil * k) / 1000;
+  const QP0 = qp0PcsKw(de, di, k);
 
   const rg_pcs = Pmfou / (Pmcons + 0.45 * QP0 + Pveil);
   const rg_pcs_dep = Pmfou_dep / (Pmcons_dep + 0.45 * QP0 + Pveil);

package/13.2_generateur_combustion_ch.spec.js

@@ -14,6 +14,7 @@ vi.mock('./enums.js', () => ({
   default: {
     type_generateur_ch: {
       89: 'chaudière gaz standard 2001-2015',
+      86: 'chaudière gaz classique 1981-1985',
       55: 'chaudière bois bûche avant 1978',
       96: 'chaudière gaz à condensation 2001-2015',
       70: 'chaudière gaz basse température 2001-2015',
@@ -50,7 +51,8 @@ const {
   periodeEmetteursDeduiteDuDpe,
   periodeEmetteurAnneeConstruction,
   PERIODES_INSTALLATION_EMETTEUR,
-  K
+  K,
+  qp0PcsKw
 } = await import('./13.2_generateur_combustion_ch.js');
 const { tv } = await import('./utils.js');
 
@@ -59,6 +61,37 @@ beforeEach(() => {
   state.bug = false;
 });
 
+// valeurs de référence de régression (générateur à air chaud, pn 20 kW, qp0 200 % de Pn)
+const RG_AIR_CHAUD = 0.8654909941701177;
+const RG_DEP_AIR_CHAUD = 0.874622228362935;
+
+/**
+ * Pertes à l'arrêt QP0 (kW sur PCS) : pour les générateurs d'air chaud, la valeur de la table 13.2
+ * est en % de Pn et doit être ramenée en W.
+ * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.1
+ */
+describe('qp0PcsKw - pertes à l’arrêt sur PCS', () => {
+  test('générateur à air chaud : qp0 ramené en W (÷ 100) puis en kW PCS', () => {
+    const di = { qp0: 1824.5 };
+    expect(qp0PcsKw({ enum_type_generateur_ch_id: '33' }, di, 1.11)).toBeCloseTo(
+      (18.245 * 1.11) / 1000,
+      12
+    );
+    expect(di.qp0).toBe(1824.5);
+  });
+
+  test('autre générateur (chaudière) : qp0 déjà en W', () => {
+    expect(qp0PcsKw({ enum_type_generateur_ch_id: '89' }, { qp0: 120 }, 1.11)).toBeCloseTo(
+      (120 * 1.11) / 1000,
+      12
+    );
+  });
+
+  test('type de générateur inconnu des enums : qp0 pris tel quel', () => {
+    expect(qp0PcsKw({ enum_type_generateur_ch_id: '0' }, { qp0: 50 }, 1)).toBeCloseTo(0.05, 12);
+  });
+});
+
 /**
  * 13.2.1.1 - Températures de fonctionnement à 30 % et 100 % de charge.
  * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.1.1
@@ -575,8 +608,40 @@ describe('calc_generateur_combustion_ch - rendement de génération', () => {
     };
     calc_generateur_combustion_ch({}, di, de, { cdimref: 0.8, cdimrefDep: 0.6 });
     // valeurs de référence de régression
-    expect(di.rg).toBeCloseTo(0.8655833354015402, 9);
-    expect(di.rg_dep).toBeCloseTo(0.8746954076850985, 9);
+    expect(di.rg).toBeCloseTo(RG_AIR_CHAUD, 9);
+    expect(di.rg_dep).toBeCloseTo(RG_DEP_AIR_CHAUD, 9);
+    // qp0 (% de Pn) n'est pas modifié par le calcul (il était divisé par 100 à chaque QPx)
+    expect(di.qp0).toBe(200);
+  });
+
+  test('générateur à air chaud : rendement indépendant du nombre de calculs (pas d’effet de bord)', () => {
+    const di = { pn: 20000, rpn: 0.9, rpint: 0.85, qp0: 200, temp_fonc_30: 40, temp_fonc_100: 70 };
+    const de = {
+      enum_type_generateur_ch_id: '33',
+      enum_type_energie_id: '2',
+      type_energie: 'gaz naturel',
+      presence_regulation_combustion: false,
+      description: 't'
+    };
+    calc_generateur_combustion_ch({}, di, de, { cdimref: 0.8, cdimrefDep: 0.6 });
+    calc_generateur_combustion_ch({}, di, de, { cdimref: 0.8, cdimrefDep: 0.6 });
+    expect(di.rg).toBeCloseTo(RG_AIR_CHAUD, 9);
+    expect(di.qp0).toBe(200);
+  });
+
+  test('générateur à air chaud : les pertes à l’arrêt réduisent le rendement', () => {
+    const base = { pn: 20000, rpn: 0.9, rpint: 0.85, temp_fonc_30: 40, temp_fonc_100: 70 };
+    const de = {
+      enum_type_generateur_ch_id: '33',
+      enum_type_energie_id: '2',
+      type_energie: 'gaz naturel',
+      description: 't'
+    };
+    const sansPertes = { ...base, qp0: 0 };
+    const avecPertes = { ...base, qp0: 200 };
+    calc_generateur_combustion_ch({}, sansPertes, de, { cdimref: 0.8, cdimrefDep: 0.6 });
+    calc_generateur_combustion_ch({}, avecPertes, de, { cdimref: 0.8, cdimrefDep: 0.6 });
+    expect(avecPertes.rg).toBeLessThan(sansPertes.rg);
   });
 
   test('chaudière basse température : coefficients dédiés (a = 0.1, b = 40)', () => {
@@ -710,6 +775,88 @@ describe('calc_generateur_combustion_ch - rendement de génération', () => {
     expect(di.rg).toBeCloseTo(0.8511886946001722, 12);
   });
 
+  /**
+   * Référence externe : autotest État CSTB `MI4-0-1`, moteur Moteur_DPE.dll (Tribu). Chaudière gaz
+   * classique 1981-1985 avec veilleuse, Pn = 24 kW, valeurs par défaut (Rpn, Rpint, Qp0 = 480 W,
+   * Pveil = 150 W, Tfonc30/100 = 56/70, Cdimref = 0.05 / TchFinal(Ch_5)).
+   * Comme QP0, la puissance de veilleuse est exprimée sur PCS (Pveil × k) dans le dénominateur.
+   * Avant correctif (Pveil sur PCI) : rg = 0.740607866698453 (Cch −0,42 %).
+   * @see : Methode_de_calcul_3CL_DPE_2021-338.pdf - §13.2.1.2
+   */
+  test('autotest État MI4-0-1 : rg identique à Tribu (Pveil exprimée sur PCS)', () => {
+    const di = {
+      pn: 24000,
+      rpn: 0.8676042248342321,
+      rpint: 0.8414063372513482,
+      qp0: 480,
+      pveilleuse: 150,
+      pveil: 150,
+      temp_fonc_30: 56,
+      temp_fonc_100: 70
+    };
+    const de = {
+      enum_type_generateur_ch_id: '86',
+      enum_type_energie_id: '2',
+      type_energie: 'gaz naturel',
+      presence_regulation_combustion: true,
+      description: 'Chaudière classique 81-85'
+    };
+    const cdimref = 0.05 / 0.01710986374076308;
+    calc_generateur_combustion_ch({}, di, de, { cdimref, cdimrefDep: cdimref });
+    // Rg Tribu (sortie MI4-0-1_Sortie.xml) : 0.7374941696172325
+    expect(di.rg).toBeCloseTo(0.7374941696172325, 12);
+  });
+
+  test('la veilleuse est convertie sur PCS : rg_pcs = Pmfou / (Pmcons + 0.45 × QP0 × k + Pveil × k)', () => {
+    const base = () => ({
+      pn: 24000,
+      rpn: 0.87,
+      rpint: 0.84,
+      qp0: 0,
+      pveilleuse: 150,
+      pveil: 150,
+      temp_fonc_30: 56,
+      temp_fonc_100: 70
+    });
+    const de = {
+      enum_type_generateur_ch_id: '86',
+      enum_type_energie_id: '2',
+      presence_regulation_combustion: true,
+      description: 't'
+    };
+    const avecVeilleuse = base();
+    calc_generateur_combustion_ch(
+      {},
+      avecVeilleuse,
+      { ...de, type_energie: 'gaz naturel' },
+      {
+        cdimref: 3
+      }
+    );
+    const sansVeilleuse = { ...base(), pveilleuse: 0, pveil: 0 };
+    calc_generateur_combustion_ch(
+      {},
+      sansVeilleuse,
+      { ...de, type_energie: 'gaz naturel' },
+      {
+        cdimref: 3
+      }
+    );
+    // rg = k × Pmfou / Pmcons sans veilleuse ; avec veilleuse : k × Pmfou / (Pmcons + 0.15 × k)
+    const k = 1.11;
+    const coefPond = { 0.05: 0.1, 0.15: 0.25, 0.25: 0.2, 0.35: 0.15, 0.45: 0.1, 0.55: 0.1 };
+    Object.assign(coefPond, { 0.65: 0.05, 0.75: 0.025, 0.85: 0.025, 0.95: 0 });
+    const pmfou = Object.entries(coefPond).reduce(
+      (acc, [x, c]) => acc + 24 * Math.min(1, Number(x) / 3) * c,
+      0
+    );
+    const pmcons = (k * pmfou) / sansVeilleuse.rg;
+    const attendu = (k * pmfou) / (pmcons + 0.15 * k);
+    expect(avecVeilleuse.rg).toBeCloseTo(attendu, 12);
+    // sur PCI (avant correctif), le rendement serait plus élevé
+    expect(avecVeilleuse.rg).toBeLessThan((k * pmfou) / (pmcons + 0.15));
+  });
+
   test('les pertes à charge nulle QP0 réduisent le rendement de génération', () => {
     const di = () => ({
       pn: 23000,

package/README.md

@@ -370,57 +370,57 @@ Ces résultats sont **générés automatiquement** à la fin de `npm run test:co
 
 <!-- CORPUS:START -->
 
-> **Version `1.8.23`** · branche `main` · généré le 2026-10-04
+> **Version `1.8.25`** · branche `main` · généré le 2026-10-08
 > Seuil de tolérance **5%**
 
 <table>
 <tr>
 <td align="center"><strong>9</strong><br/><sub>corpus</sub></td>
-<td align="center"><strong>89 990</strong><br/><sub>DPE analysés</sub></td>
-<td align="center"><strong>60 717</strong><br/><sub>DPE conformes</sub></td>
-<td align="center"><strong>67,47 %</strong><br/><sub>réussite globale</sub></td>
+<td align="center"><strong>89 989</strong><br/><sub>DPE analysés</sub></td>
+<td align="center"><strong>60 904</strong><br/><sub>DPE conformes</sub></td>
+<td align="center"><strong>67,68 %</strong><br/><sub>réussite globale</sub></td>
 </tr>
 </table>
 
 |     | Corpus                                                                                                                  |    Réussite |                        |  DPE conformes |
 | :-: | :---------------------------------------------------------------------------------------------------------------------- | ----------: | :--------------------- | -------------: |
-| 🔴  | **Généraliste**<br/><sub>`corpus_dpe.csv`</sub>                                                                         | **46,02 %** | `█████████░░░░░░░░░░░` | 4 602 / 10 000 |
+| 🔴  | **Généraliste**<br/><sub>`corpus_dpe.csv`</sub>                                                                         | **46,03 %** | `█████████░░░░░░░░░░░` | 4 603 / 10 000 |
 | 🟢  | **Appartement · chauffage individuel (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_individuel_2025.csv`</sub> | **92,53 %** | `███████████████████░` | 9 253 / 10 000 |
-| 🟢  | **Logement individuel (2025)**<br/><sub>`dpe_logement_individuel_2025.csv`</sub>                                        | **89,49 %** | `██████████████████░░` |  8 942 / 9 992 |
-| 🟢  | **Maison individuelle (2025)**<br/><sub>`dpe_maison_individuelle_2025.csv`</sub>                                        | **88,58 %** | `██████████████████░░` | 8 858 / 10 000 |
-| 🟡  | **Immeuble · chauffage individuel**<br/><sub>`dpe_immeuble_chauffage_individuel.csv`</sub>                              | **76,89 %** | `███████████████░░░░░` |  7 688 / 9 999 |
-| 🟡  | **Appartement · chauffage collectif (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_collectif_2025.csv`</sub>   | **71,24 %** | `██████████████░░░░░░` | 7 124 / 10 000 |
-| 🟡  | **Immeuble · chauffage collectif**<br/><sub>`dpe_immeuble_chauffage_collectif.csv`</sub>                                | **63,93 %** | `█████████████░░░░░░░` |  6 392 / 9 999 |
-| 🔴  | **Immeuble · chauffage mixte**<br/><sub>`dpe_immeuble_chauffage_mixte.csv`</sub>                                        | **51,64 %** | `██████████░░░░░░░░░░` | 5 164 / 10 000 |
+| 🟢  | **Logement individuel (2025)**<br/><sub>`dpe_logement_individuel_2025.csv`</sub>                                        | **89,54 %** | `██████████████████░░` |  8 947 / 9 992 |
+| 🟢  | **Maison individuelle (2025)**<br/><sub>`dpe_maison_individuelle_2025.csv`</sub>                                        | **88,59 %** | `██████████████████░░` | 8 859 / 10 000 |
+| 🟡  | **Immeuble · chauffage individuel**<br/><sub>`dpe_immeuble_chauffage_individuel.csv`</sub>                              | **76,90 %** | `███████████████░░░░░` |  7 688 / 9 998 |
+| 🟡  | **Appartement · chauffage collectif (2025)**<br/><sub>`dpe_appartement_individuel_chauffage_collectif_2025.csv`</sub>   | **71,59 %** | `██████████████░░░░░░` | 7 159 / 10 000 |
+| 🟡  | **Immeuble · chauffage collectif**<br/><sub>`dpe_immeuble_chauffage_collectif.csv`</sub>                                | **65,11 %** | `█████████████░░░░░░░` |  6 510 / 9 999 |
+| 🔴  | **Immeuble · chauffage mixte**<br/><sub>`dpe_immeuble_chauffage_mixte.csv`</sub>                                        | **51,91 %** | `██████████░░░░░░░░░░` | 5 191 / 10 000 |
 | 🔴  | **Individuel généré depuis l'immeuble (2026)**<br/><sub>`dpe_individuel_a_partir_dpe_immeuble_2026.csv`</sub>           | **26,94 %** | `█████░░░░░░░░░░░░░░░` | 2 694 / 10 000 |
 
 <sub>🟢 ≥ 85 % · 🟡 ≥ 60 % · 🔴 < 60 %</sub>
 
 #### Temps d’exécution
 
-> Durée de l’appel à `calcul_3cl` par DPE, sur 89 956 calculs.
+> Durée de l’appel à `calcul_3cl` par DPE, sur 89 955 calculs.
 > La copie défensive de l’entrée et la lecture du fichier sont exclues de la mesure.
 
 <table>
 <tr>
-<td align="center"><strong>5,32 ms</strong><br/><sub>moyenne</sub></td>
-<td align="center"><strong>4,66 ms</strong><br/><sub>médiane</sub></td>
-<td align="center"><strong>0,12 ms</strong><br/><sub>min</sub></td>
-<td align="center"><strong>933,6 ms</strong><br/><sub>max</sub></td>
+<td align="center"><strong>5,91 ms</strong><br/><sub>moyenne</sub></td>
+<td align="center"><strong>5,23 ms</strong><br/><sub>médiane</sub></td>
+<td align="center"><strong>0,17 ms</strong><br/><sub>min</sub></td>
+<td align="center"><strong>900,0 ms</strong><br/><sub>max</sub></td>
 </tr>
 </table>
 
 | Corpus                                         | Moyenne | Médiane |     Min |      Max |     p95 |     p99 |
 | :--------------------------------------------- | ------: | ------: | ------: | -------: | ------: | ------: |
-| **Généraliste**                                | 4,90 ms | 4,46 ms | 0,22 ms |  38,9 ms | 8,82 ms | 12,0 ms |
-| **Appartement · chauffage individuel (2025)**  | 3,82 ms | 3,51 ms | 1,34 ms |  30,2 ms | 6,52 ms | 8,52 ms |
-| **Logement individuel (2025)**                 | 5,25 ms | 4,67 ms | 1,66 ms |  38,9 ms | 9,78 ms | 13,2 ms |
-| **Maison individuelle (2025)**                 | 6,36 ms | 5,77 ms | 2,05 ms |  50,5 ms | 11,3 ms | 14,5 ms |
-| **Immeuble · chauffage individuel**            | 6,23 ms | 5,25 ms | 0,12 ms | 122,1 ms | 12,5 ms | 20,5 ms |
-| **Appartement · chauffage collectif (2025)**   | 4,33 ms | 3,99 ms | 1,46 ms |  35,9 ms | 7,01 ms | 9,22 ms |
-| **Immeuble · chauffage collectif**             | 5,18 ms | 4,35 ms | 0,12 ms | 129,8 ms | 9,51 ms | 17,0 ms |
-| **Immeuble · chauffage mixte**                 | 5,66 ms | 4,59 ms | 0,22 ms | 933,6 ms | 11,8 ms | 18,9 ms |
-| **Individuel généré depuis l'immeuble (2026)** | 6,18 ms | 5,35 ms | 1,62 ms |  44,3 ms | 11,6 ms | 17,6 ms |
+| **Généraliste**                                | 6,25 ms | 5,23 ms | 0,22 ms |  96,9 ms | 13,6 ms | 22,4 ms |
+| **Appartement · chauffage individuel (2025)**  | 4,81 ms | 4,51 ms | 1,55 ms |  33,1 ms | 8,15 ms | 11,3 ms |
+| **Logement individuel (2025)**                 | 5,17 ms | 4,63 ms | 1,64 ms |  36,8 ms | 9,69 ms | 12,9 ms |
+| **Maison individuelle (2025)**                 | 6,35 ms | 5,77 ms | 1,92 ms |  56,4 ms | 11,3 ms | 14,4 ms |
+| **Immeuble · chauffage individuel**            | 7,51 ms | 6,50 ms | 0,17 ms | 108,3 ms | 14,9 ms | 23,8 ms |
+| **Appartement · chauffage collectif (2025)**   | 5,06 ms | 4,89 ms | 1,65 ms |  36,1 ms | 8,09 ms | 11,2 ms |
+| **Immeuble · chauffage collectif**             | 6,54 ms | 5,84 ms | 0,22 ms | 114,6 ms | 11,9 ms | 20,0 ms |
+| **Immeuble · chauffage mixte**                 | 5,63 ms | 4,63 ms | 0,18 ms | 900,0 ms | 11,6 ms | 18,6 ms |
+| **Individuel généré depuis l'immeuble (2026)** | 5,87 ms | 5,07 ms | 1,71 ms |  46,4 ms | 11,1 ms | 17,6 ms |
 
 <sub>La moyenne est tirée vers le haut par les DPE collectifs, dont le coût atteint plusieurs
 dizaines de fois la médiane : c’est la médiane qui décrit le cas courant, et p95/p99 le cas

package/engine.js

@@ -33,7 +33,7 @@ import getFicheTechnique from './ficheTechnique.js';
 import { ProductionENR } from './16.2_production_enr.js';
 import DpeSanitizerService from './dpe-sanitizer.service.js';
 
-const LIB_VERSION = '1.8.24';
+const LIB_VERSION = '1.8.26';
 
 const dpeSanitizerService = new DpeSanitizerService();
 

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@open3cl/engine",
-  "version": "1.8.24",
+  "version": "1.8.26",
   "description": "Open Source 3CL-DPE engine",
   "main": "index.js",
   "directories": {