@open3cl/engine 1.0.11 → 1.0.13

package/features/engine/domain/apport_et_besoin/ch/besoin-ch.service.js

@@ -0,0 +1,150 @@
+import { mois_liste } from '../../../../../utils.js';
+import { inject } from 'dioma';
+import { FrTvStore } from '../../../../dpe/infrastructure/froid/frTv.store.js';
+
+/**
+ * Calcul du besoin en chauffage
+ * Chapitre 2 - Expression du besoin de chauffage
+ *
+ * Methode_de_calcul_3CL_DPE_2021 - Page 7
+ * Octobre 2021
+ * @see consolide_anne…arrete_du_31_03_2021_relatif_aux_methodes_et_procedures_applicables.pdf
+ */
+export class BesoinChService {
+  /**
+   * @type {FrTvStore}
+   */
+  #tvStore;
+
+  /**
+   * @param tvStore {FrTvStore}
+   */
+  constructor(tvStore = inject(FrTvStore)) {
+    this.#tvStore = tvStore;
+  }
+
+  /**
+   * Besoin de chauffage hors pertes récupérées
+   * 9.1.1 - Consommation de chauffage
+   *
+   * @param ctx {Contexte}
+   * @param logement {Logement}
+   * @return {{besoin_ch_hp: number, besoin_ch_depensier_hp: number, fraction_apport_gratuit_ch: number, fraction_apport_gratuit_depensier_ch: number}}
+   */
+  execute(ctx, logement) {
+    const besoinCh = mois_liste.reduce(
+      (acc, mois) => {
+        const dh19 = this.#tvStore.getData(
+          'dh19',
+          ctx.altitude.value,
+          ctx.zoneClimatique.value,
+          mois,
+          ctx.inertie.ilpa
+        );
+        const dh21 = this.#tvStore.getData(
+          'dh21',
+          ctx.altitude.value,
+          ctx.zoneClimatique.value,
+          mois,
+          ctx.inertie.ilpa
+        );
+
+        const besoinCh = this.besoinChHorsPertesRecuperees(ctx, logement, mois, dh19);
+        acc.besoin_ch_hp += besoinCh.besoinChMoisHP;
+        acc.fraction_apport_gratuit_ch += besoinCh.fractionApportGratuitMois;
+        logement.donnees_de_calcul.besoinChauffageHP[mois] = besoinCh.besoinChMoisHP;
+
+        const besoinChDepensier = this.besoinChHorsPertesRecuperees(ctx, logement, mois, dh21);
+        acc.besoin_ch_depensier_hp += besoinChDepensier.besoinChMoisHP;
+        acc.fraction_apport_gratuit_depensier_ch += besoinChDepensier.fractionApportGratuitMois;
+        logement.donnees_de_calcul.besoinChauffageDepensierHP[mois] =
+          besoinChDepensier.besoinChMoisHP;
+
+        acc.dh19 += dh19;
+        acc.dh21 += dh21;
+        return acc;
+      },
+      {
+        besoin_ch_hp: 0,
+        besoin_ch_depensier_hp: 0,
+        fraction_apport_gratuit_ch: 0,
+        fraction_apport_gratuit_depensier_ch: 0,
+        dh19: 0,
+        dh21: 0
+      }
+    );
+
+    besoinCh.fraction_apport_gratuit_ch /= besoinCh.dh19;
+    besoinCh.fraction_apport_gratuit_depensier_ch /= besoinCh.dh21;
+
+    delete besoinCh.dh19;
+    delete besoinCh.dh19;
+
+    return besoinCh;
+  }
+
+  /**
+   * Fraction des besoins de chauffage couverts par les apports gratuits pour un mois donné
+   * 6.1 Calcul de F
+   *
+   * @param ctx {Contexte}
+   * @param logement {Logement}
+   * @param mois {string}
+   * @param dh {number} - degrés-heures de chauffage sur le mois j (°Ch)
+   * @returns {number}
+   */
+  fractionBesoinCh(ctx, logement, mois, dh) {
+    // Apports internes dans le logement sur le mois
+    const Ai = logement.donnees_de_calcul.apportsInterneCh[mois];
+    // Apports solaires dans le logement sur le mois durant la période de chauffe
+    const As = logement.donnees_de_calcul.apportsSolaire[mois];
+
+    if (dh === 0) return 0;
+
+    let pow;
+    switch (ctx.inertie.id) {
+      // Inertie très lourde ou lourde
+      case 1:
+      case 2:
+        pow = 3.6;
+        break;
+      // Inertie moyenne
+      case 3:
+        pow = 2.9;
+        break;
+      // Inertie légère
+      case 4:
+        pow = 2.5;
+        break;
+    }
+
+    const Xj = (As + Ai) / (logement.sortie.deperdition.deperdition_enveloppe * dh);
+    return (Xj - Xj ** pow) / (1 - Xj ** pow);
+  }
+
+  /**
+   * Besoin de chauffage hors pertes récupérées sur le mois (kWh) :
+   * 9.1.1 - Consommation de chauffage
+   *
+   * Fraction des besoins de chauffage couverts par les apports gratuits sur le mois
+   * 6.1 - Calcul de F
+   *
+   * @param ctx {Contexte}
+   * @param logement {Logement}
+   * @param mois {string}
+   * @param dh {number} - degrés-heures de chauffage sur le mois j (°Ch)
+   * @returns {{besoinChMoisHP: number, fractionApportGratuitMois: number}}
+   */
+  besoinChHorsPertesRecuperees(ctx, logement, mois, dh) {
+    // Fraction des besoins de chauffage couverts par les apports gratuits pour un mois donné
+    const F = this.fractionBesoinCh(ctx, logement, mois, dh);
+
+    // Besoin de chauffage d’un logement par kelvin sur le mois j (W/K)
+    const BV = logement.sortie.deperdition.deperdition_enveloppe * (1 - F);
+
+    return {
+      besoinChMoisHP: (BV * dh) / 1000,
+      fractionApportGratuitMois: F * dh
+    };
+  }
+}

package/features/engine/domain/apport_et_besoin/ch/besoin-ch.service.spec.js

@@ -0,0 +1,145 @@
+import { beforeEach, describe, expect, test, vi } from 'vitest';
+import { FrTvStore } from '../../../../dpe/infrastructure/froid/frTv.store.js';
+import { BesoinChService } from './besoin-ch.service.js';
+import { mois_liste } from '../../../../../utils.js';
+
+/** @type {BesoinChService} **/
+let service;
+
+/** @type {FrTvStore} **/
+let tvStore;
+
+describe('Calcul des besoins en chauffage du logement', () => {
+  beforeEach(() => {
+    tvStore = new FrTvStore();
+    service = new BesoinChService(tvStore);
+  });
+
+  test.each([
+    {
+      inertie: 1,
+      expected: 0.24
+    },
+    {
+      inertie: 2,
+      expected: 0.24
+    },
+    {
+      inertie: 3,
+      expected: 0.23
+    },
+    {
+      inertie: 4,
+      expected: 0.22
+    }
+  ])(
+    'Determination de la fraction des besoins de chauffage couverts par les apports gratuits pour un mois donné avec inertie = $inertie',
+    ({ inertie, expected }) => {
+      /** @type {Contexte} */
+      const ctx = { inertie: { id: inertie } };
+
+      /** @type { Logement } **/
+      const logement = {
+        donnees_de_calcul: {
+          apportsInterneCh: { Janvier: 102.5 },
+          apportsSolaire: { Janvier: 12.5 }
+        },
+        sortie: { deperdition: { deperdition_enveloppe: 25.9 } }
+      };
+      expect(service.fractionBesoinCh(ctx, logement, 'Janvier', 18)).toBeCloseTo(expected, 2);
+    }
+  );
+
+  test('Besoin de chauffage hors pertes récupérées et fraction des apports gratuits pour un mois', () => {
+    /** @type {Contexte} */
+    const ctx = { inertie: { id: 1 } };
+
+    /** @type { Logement } **/
+    const logement = {
+      donnees_de_calcul: {
+        apportsInterneCh: { Janvier: 102.5 },
+        apportsSolaire: { Janvier: 12.5 }
+      },
+      sortie: { deperdition: { deperdition_enveloppe: 25.9 } }
+    };
+
+    const besoinHP = service.besoinChHorsPertesRecuperees(ctx, logement, 'Janvier', 12.5);
+    expect(besoinHP.besoinChMoisHP).toBeCloseTo(0.2139, 4);
+    expect(besoinHP.fractionApportGratuitMois).toBeCloseTo(4.2412, 4);
+  });
+
+  test('Besoin total de chauffage hors pertes récupérées et fraction des apports gratuits', () => {
+    vi.spyOn(tvStore, 'getData').mockReturnValue(10);
+
+    /** @type {Contexte} */
+    const ctx = {
+      inertie: { id: 1, ilpa: 1 },
+      altitude: { value: '400-800m' },
+      zoneClimatique: { value: 'h3c' }
+    };
+
+    /** @type { Logement } **/
+    const logement = {
+      donnees_de_calcul: {
+        apportsInterneCh: {
+          Janvier: 102.5,
+          Février: 102.5,
+          Mars: 102.5,
+          Avril: 102.5,
+          Mai: 102.5,
+          Juin: 102.5,
+          Juillet: 102.5,
+          Aout: 102.5,
+          Septembre: 102.5,
+          Octobre: 102.5,
+          Novembre: 102.5,
+          Décembre: 102.5
+        },
+        apportsSolaire: {
+          Janvier: 102.5,
+          Février: 102.5,
+          Mars: 102.5,
+          Avril: 102.5,
+          Mai: 102.5,
+          Juin: 102.5,
+          Juillet: 102.5,
+          Aout: 102.5,
+          Septembre: 102.5,
+          Octobre: 102.5,
+          Novembre: 102.5,
+          Décembre: 102.5
+        },
+        besoinChauffageHP: {},
+        besoinChauffageDepensierHP: {}
+      },
+      sortie: { deperdition: { deperdition_enveloppe: 25.9 } }
+    };
+
+    const besoinTotalHP = service.execute(ctx, logement);
+    expect(besoinTotalHP.besoin_ch_hp).toBeCloseTo(1.1388, 4);
+    expect(besoinTotalHP.besoin_ch_depensier_hp).toBeCloseTo(1.1388, 4);
+    expect(besoinTotalHP.fraction_apport_gratuit_ch).toBeCloseTo(0.6336, 4);
+    expect(besoinTotalHP.fraction_apport_gratuit_depensier_ch).toBeCloseTo(0.6336, 4);
+
+    for (const mois of mois_liste) {
+      expect(tvStore.getData).toHaveBeenCalledWith(
+        'dh19',
+        ctx.altitude.value,
+        ctx.zoneClimatique.value,
+        mois,
+        ctx.inertie.ilpa
+      );
+
+      expect(tvStore.getData).toHaveBeenCalledWith(
+        'dh21',
+        ctx.altitude.value,
+        ctx.zoneClimatique.value,
+        mois,
+        ctx.inertie.ilpa
+      );
+
+      expect(logement.donnees_de_calcul.besoinChauffageHP[mois]).toBeCloseTo(0.0949, 4);
+      expect(logement.donnees_de_calcul.besoinChauffageDepensierHP[mois]).toBeCloseTo(0.0949, 4);
+    }
+  });
+});

package/features/engine/domain/apport_et_besoin/ch/perte-ch-recup.service.js

@@ -0,0 +1,168 @@
+import { mois_liste } from '../../../../../utils.js';
+import { inject } from 'dioma';
+import { FrTvStore } from '../../../../dpe/infrastructure/froid/frTv.store.js';
+
+/**
+ * Calcul des pertes récupérées de génération de chauffage
+ * Chapitre 9.1.1 Consommation de chauffage
+ *
+ * Données calculées
+ *  — pertes_generateur_ch_recup : pertes récupérées de génération pour le chauffage (kWh)
+ *  — pertes_generateur_ch_recup_depensier : pertes récupérées de génération pour le chauffage (kWh) pour le scénario dépensier
+ *
+ * Methode_de_calcul_3CL_DPE_2021 - Page 57
+ * Octobre 2021
+ * @see consolide_anne…arrete_du_31_03_2021_relatif_aux_methodes_et_procedures_applicables.pdf
+ */
+export class PerteChRecupService {
+  /**
+   * @type {FrTvStore}
+   */
+  #tvStore;
+
+  /**
+   * @param tvStore {FrTvStore}
+   */
+  constructor(tvStore = inject(FrTvStore)) {
+    this.#tvStore = tvStore;
+  }
+
+  /**
+   * Pertes récupérées des générateurs de chauffage pour le chauffage (kWh)
+   *
+   * @param ctx {Contexte}
+   * @param logement {Logement}
+   * @return {{pertes_generateur_ch_recup: number, pertes_generateur_ch_recup_depensier: number}}
+   */
+  execute(ctx, logement) {
+    return {
+      pertes_generateur_ch_recup: this.pertesGenerateurChRecup(ctx, logement, false) / 1000,
+      pertes_generateur_ch_recup_depensier: this.pertesGenerateurChRecup(ctx, logement, true) / 1000
+    };
+  }
+
+  /**
+   * Pertes récupérées des générateurs de chauffage pour le chauffage (Wh)
+   * 9.1.1 Consommation de chauffage
+   *
+   * @param ctx {Contexte}
+   * @param logement {Logement}
+   * @param depensier {boolean}
+   * @returns {number}
+   */
+  pertesGenerateurChRecup(ctx, logement, depensier) {
+    const generateursWithPertesGeneration = this.generateursWithPertesGeneration(logement);
+
+    return mois_liste.reduce((acc, mois) => {
+      return (
+        acc +
+        this.Qrec(
+          generateursWithPertesGeneration,
+          this.#tvStore.getData(
+            depensier ? 'nref21' : 'nref19',
+            ctx.altitude.value,
+            ctx.zoneClimatique.value,
+            mois,
+            ctx.inertie.ilpa
+          ),
+          depensier
+            ? logement.donnees_de_calcul.besoinChauffageDepensierHP[mois]
+            : logement.donnees_de_calcul.besoinChauffageHP[mois]
+        )
+      );
+    }, 0);
+  }
+
+  /**
+   * Pertes récupérées des générateurs de chauffage pour le chauffage
+   * 9.1.1 Consommation de chauffage
+   *
+   * Seules les pertes des générateurs en volume chauffé sont récupérables.
+   * Les pertes récupérées des générateurs d’air chaud sont nulles.
+   *
+   * @param logement {Logement}
+   */
+  generateursWithPertesGeneration(logement) {
+    const installationsCh =
+      logement.installation_chauffage_collection?.installation_chauffage || [];
+
+    return installationsCh.reduce((acc, installation) => {
+      // Liste des générateurs de chauffage pour lesquels il y a une récupération d'énergie
+      return acc.concat(
+        (installation.generateur_chauffage_collection?.generateur_chauffage || []).filter(
+          (generateurChauffage) => {
+            const generateurChauffageDE = generateurChauffage.donnee_entree;
+
+            /**
+             * 50 - générateur à air chaud à combustion avant 2006
+             * 51 - générateur à air chaud à combustion standard a partir de 2006
+             * 52 - générateur à air chaud à combustion à condensation a partir de 2006
+             */
+            return (
+              (generateurChauffageDE.position_volume_chauffe ?? 0) === 1 &&
+              ![50, 51, 52].includes(parseInt(generateurChauffageDE.enum_type_generateur_ch_id))
+            );
+          }
+        )
+      );
+    }, []);
+  }
+
+  /**
+   * Pertes récupérées de génération pour le chauffage sur le mois j (Wh)
+   *
+   * @param generateurs {GenerateurChauffage[]}
+   * @param nref {number}
+   * @param besoinChauffageMois {number}
+   * @constructor
+   */
+  Qrec(generateurs, nref, besoinChauffageMois) {
+    return generateurs.reduce((acc, generateur) => {
+      const generateurChauffageDE = generateur.donnee_entree;
+      const generateurChauffageDI = generateur.donnee_intermediaire;
+
+      /**
+       * Part des pertes par les parois égale à 0,75 pour les équipements à ventouse ou assistés par ventilateur
+       * et 0,5 pour les autres
+       * @type {number}
+       */
+      const Cper = parseInt(generateurChauffageDE.presence_ventouse || 0) === 1 ? 0.75 : 0.5;
+
+      /**
+       * Durée pendant laquelle les pertes sont récupérées sur le mois (h)
+       */
+      let Dperj;
+      switch (parseInt(generateurChauffageDE.enum_usage_generateur_id)) {
+        case 1:
+          // Pour les générateurs assurant le chauffage uniquement
+          Dperj = Math.min(
+            nref,
+            this.#pertes_gen_ch(besoinChauffageMois, generateurChauffageDI.pn)
+          );
+          break;
+        case 2:
+          // Pour les générateurs assurant l’ECS uniquement
+          Dperj = this.#pertes_gen_ecs(nref);
+          break;
+        case 3:
+          // Pour les générateurs assurant le chauffage et l’ECS
+          Dperj = Math.min(
+            nref,
+            this.#pertes_gen_ch(besoinChauffageMois, generateurChauffageDI.pn) +
+              this.#pertes_gen_ecs(nref)
+          );
+          break;
+      }
+
+      return acc + 0.48 * Cper * generateurChauffageDI.qp0 * Dperj || 0;
+    }, 0);
+  }
+
+  #pertes_gen_ch(Bch, pn) {
+    return (1.3 * Bch * 1000) / (0.3 * pn);
+  }
+
+  #pertes_gen_ecs(nref) {
+    return (nref * 1790) / 8760;
+  }
+}