@open3cl/engine 1.0.0

package/9_chauffage.js

@@ -0,0 +1,291 @@
+import { requestInput, Tbase, tv, tvColumnIDs } from './utils.js';
+import { calc_emetteur_ch } from './9_emetteur_ch.js';
+import {
+  calc_generateur_ch,
+  checkForGeneratorType,
+  hasConsoForAuxDistribution
+} from './9_generateur_ch.js';
+import { tv_generateur_combustion } from './13.2_generateur_combustion.js';
+import { tv_temp_fonc_30_100 } from './13.2_generateur_combustion_ch.js';
+import enums from './enums.js';
+
+export default function calc_chauffage(
+  dpe,
+  ch,
+  ca_id,
+  zc_id,
+  inertie_id,
+  map_id,
+  bch,
+  bch_dep,
+  GV,
+  Sh,
+  hsp,
+  ac,
+  ilpa
+) {
+  const de = ch.donnee_entree;
+  const di = {};
+  const du = {};
+
+  const ca = enums.classe_altitude[ca_id];
+  const zc = enums.zone_climatique[zc_id];
+  const tbase = Tbase[ca][zc.slice(0, 2)];
+
+  di.besoin_ch = bch;
+  di.besoin_ch_depensier = bch_dep;
+
+  const em_ch = ch.emetteur_chauffage_collection.emetteur_chauffage;
+  em_ch.forEach((em_ch) => {
+    calc_emetteur_ch(em_ch, de, map_id, inertie_id);
+  });
+
+  if (de.surface_chauffee !== Sh) {
+    bch *= de.surface_chauffee / Sh;
+  }
+
+  const Fch = tv_ch_facteur_couverture_solaire(de, zc_id);
+  const cfg_id = requestInput(de, du, 'cfg_installation_ch');
+  const gen_ch = ch.generateur_chauffage_collection.generateur_chauffage;
+
+  gen_ch.forEach((gen) => {
+    const genChDe = gen.donnee_entree;
+    const genChDu = gen.donnee_utilisateur || {};
+    const genChDi = gen.donnee_intermediaire || {};
+
+    genChDe.ratio_virtualisation = de.ratio_virtualisation || 1;
+    genChDe.surface_chauffee = de.surface_chauffee || Sh;
+    genChDe.nombre_niveau_installation_ch = de.nombre_niveau_installation_ch || 1;
+    genChDe.fch = Fch || 0.5;
+
+    // Récupération du type de générateur
+    checkForGeneratorType(dpe, genChDe, genChDi, genChDu);
+
+    if (genChDu.isCombustionGenerator) {
+      const methodeSaisie = parseInt(genChDe.enum_methode_saisie_carac_sys_id);
+      tv_generateur_combustion(dpe, genChDi, genChDe, 'ch', GV, tbase, methodeSaisie);
+
+      const type_gen_ch_list = tvColumnIDs('temp_fonc_30', 'type_generateur_ch');
+      if (type_gen_ch_list.includes(genChDe.enum_type_generateur_ch_id)) {
+        /**
+         * Si la méthode de saisie n'est pas "Valeur forfaitaire" mais "caractéristiques saisies"
+         * Documentation 3CL : "Pour les installations récentes ou recommandées, les caractéristiques réelles des chaudières présentées sur les bases
+         * de données professionnelles peuvent être utilisées."
+         *
+         * 5 - caractéristiques saisies à partir de la plaque signalétique ou d'une documentation technique du système à combustion : pn, rpn,rpint,qp0,temp_fonc_30,temp_fonc_100
+         */
+        if (methodeSaisie !== 5 || !genChDi.temp_fonc_30 || !genChDi.temp_fonc_100) {
+          tv_temp_fonc_30_100(genChDi, genChDe, genChDu, em_ch, ac);
+        }
+      }
+    }
+
+    gen.donnee_intermediaire = genChDi;
+    gen.donnee_utilisateur = genChDu;
+    gen.donnee_entree = genChDe;
+  });
+
+  /**
+   * 13.2.1.3 Cascade de deux générateurs à combustion
+   * Une puissance relative pour chaque générateur est calculée et appliquée à la conso globale de chauffage
+   * Seuls les générateurs en cascade sont concernés
+   *
+   * 9.1.4.3 Les pompes à chaleur hybrides
+   * Cas particulier des PAC hybrides avec répartition forfaitaire du besoin
+   * @type {number|number}
+   */
+  const Pnominal = gen_ch.reduce((acc, gen) => acc + (gen.donnee_intermediaire.pn || 0), 0);
+  du.Pnominal = Pnominal;
+
+  const nbCascadeAndCombustion = gen_ch.filter(
+    (value) =>
+      value.donnee_utilisateur.isCombustionGenerator &&
+      Number.parseInt(de.enum_cfg_installation_ch_id) === 1
+  ).length;
+
+  // Nombre de générateurs avec une consommation des auxiliaires de distribution
+  const nbGenWithAuxConsoDistribution = gen_ch.reduce((acc, gen) => {
+    if (hasConsoForAuxDistribution(gen.donnee_entree.enum_type_generateur_ch_id)) {
+      acc++;
+    }
+    return acc;
+  }, 0);
+
+  gen_ch.forEach((gen, _pos) => {
+    // Nombre de générateurs en cascade pour un même émetteur si l'installation est une installation simple
+    const nbEmetteur = em_ch.filter(
+      (em) =>
+        em.donnee_entree.enum_lien_generateur_emetteur_id ===
+        gen.donnee_entree.enum_lien_generateur_emetteur_id
+    ).length;
+    const nbCascadeForSameEmetteur =
+      Number.parseInt(de.enum_cfg_installation_ch_id) === 1 && nbEmetteur === 1
+        ? gen_ch.filter(
+            (gen_ch) =>
+              gen_ch.donnee_entree.enum_lien_generateur_emetteur_id ===
+              gen.donnee_entree.enum_lien_generateur_emetteur_id
+          ).length
+        : 0;
+
+    const prorataGenerateur = getProrataGenerateur(
+      gen,
+      nbCascadeAndCombustion,
+      nbCascadeForSameEmetteur,
+      Pnominal,
+      zc
+    );
+    (gen.donnee_utilisateur = gen.donnee_utilisateur || {}).nbGenerateurCascade = gen_ch.length;
+
+    calc_generateur_ch(
+      dpe,
+      gen,
+      _pos,
+      em_ch,
+      cfg_id,
+      bch * prorataGenerateur,
+      bch_dep * prorataGenerateur,
+      GV,
+      Sh,
+      hsp,
+      ca_id,
+      zc_id,
+      ilpa
+    );
+
+    // Si plusieurs générateurs de chauffage, la consommation des auxiliaires est répartie sur chacun d'eux
+    if (
+      gen.donnee_intermediaire.conso_auxiliaire_distribution_ch &&
+      nbGenWithAuxConsoDistribution > 0
+    ) {
+      gen.donnee_intermediaire.conso_auxiliaire_distribution_ch /= nbGenWithAuxConsoDistribution;
+    }
+  });
+
+  di.conso_ch = gen_ch.reduce((acc, gen_ch) => acc + gen_ch.donnee_intermediaire.conso_ch, 0);
+  di.conso_ch_depensier = gen_ch.reduce(
+    (acc, gen_ch) => acc + gen_ch.donnee_intermediaire.conso_ch_depensier,
+    0
+  );
+
+  ch.donnee_intermediaire = di;
+  ch.donnee_utilisateur = du;
+}
+
+/**
+ * 13.2.1.3 Cascade de deux générateurs à combustion
+ * Une puissance relative pour chaque générateur est calculée et appliquée à la conso globale de chauffage
+ * Seuls les générateurs en cascade sont concernés
+ *
+ * 9.1.4.3 Les pompes à chaleur hybrides
+ * Cas particulier des PAC hybrides avec répartition forfaitaire du besoin
+ * @type {number|number}
+ */
+function getProrataGenerateur(genCh, nbCascadeAndCombustion, nbGenerateurCascade, Pnominal, zc) {
+  // IDs des pompes à chaleur hybrides
+  if (
+    genCh.donnee_entree.enum_type_generateur_ch_id >= 145 &&
+    genCh.donnee_entree.enum_type_generateur_ch_id <= 170
+  ) {
+    const zone = zc.slice(0, 2);
+    const hybrideProrata = {
+      h1: {
+        pac: 0.8,
+        chaudiere: 0.2
+      },
+      h2: {
+        pac: 0.83,
+        chaudiere: 0.17
+      },
+      h3: {
+        pac: 0.88,
+        chaudiere: 0.12
+      }
+    };
+
+    // Partie PAC du générateur
+    if (genCh.donnee_intermediaire.scop || genCh.donnee_intermediaire.cop) {
+      return hybrideProrata[zone]['pac'];
+    } else {
+      return hybrideProrata[zone]['chaudiere'];
+    }
+  }
+
+  return nbCascadeAndCombustion > 1
+    ? genCh.donnee_intermediaire.pn / Pnominal
+    : 1 / (nbGenerateurCascade || 1);
+}
+
+/**
+ * Retourne le facteur de couverture solaire pour les maisons avec chauffage solaire
+ * @param de {Donnee_entree}
+ * @param zc_id {int}
+ * @returns {number|*|null}
+ */
+function tv_ch_facteur_couverture_solaire(de, zc_id) {
+  /**
+   * 18.4 Facteur de couverture solaire
+   * Les facteurs de couverture solaire peuvent être saisi directement quand ils sont connus et peuvent être justifiés.
+   */
+  if (de.fch_saisi) {
+    return de.fch_saisi;
+  }
+
+  const matcher = {
+    enum_zone_climatique_id: zc_id
+  };
+
+  const row = tv('facteur_couverture_solaire', matcher);
+  if (row) {
+    return Number(row.facteur_couverture_solaire);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour facteur_couverture_solaire !!');
+    return null;
+  }
+}
+
+/**
+ * 13.2.1.2 Présence d’un ou plusieurs générateurs à combustion indépendants
+ * Calcul du taux de charge cdimref et cdimrefDep pour chacun des générateurs
+ *
+ * @param installationChauffage {InstallationChauffageItem[]}
+ * @param GV {number} déperdition de l'enveloppe
+ * @param zcId {string} id de la zone climatique du bien
+ * @param caId {string} id de la classe d'altitude du bien
+ */
+export function tauxChargeForGenerator(installationChauffage, GV, caId, zcId) {
+  // Récupération des installations de chauffage avec générateur à combustion
+  const installChauffageWithCombustion = [];
+  installationChauffage.forEach((ch) => {
+    const gen_ch = ch.generateur_chauffage_collection.generateur_chauffage;
+
+    const genCombustion = gen_ch.reduce((acc, gen) => {
+      if (gen.donnee_utilisateur.isCombustionGenerator) {
+        acc.push(gen);
+      }
+
+      return acc;
+    }, []);
+
+    if (genCombustion.length) {
+      ch.donnee_utilisateur.genCombustion = genCombustion;
+      installChauffageWithCombustion.push(ch);
+    }
+  });
+
+  const ca = enums.classe_altitude[caId];
+  const zc = enums.zone_climatique[zcId];
+  const tbase = Tbase[ca][zc.slice(0, 2)];
+
+  // Pour N générateurs à combustion, puissance totale de tous les générateurs
+  const Pn = installChauffageWithCombustion.reduce(
+    (acc, gen) => acc + gen.donnee_utilisateur.Pnominal,
+    0
+  );
+
+  installChauffageWithCombustion.forEach((installCh) => {
+    (installCh.donnee_utilisateur.genCombustion || []).forEach((gen) => {
+      gen.donnee_utilisateur.cdimref = Pn / (GV * (19 - tbase));
+      gen.donnee_utilisateur.cdimrefDep = Pn / (GV * (21 - tbase));
+    });
+  });
+}

package/9_chauffage.spec.js

@@ -0,0 +1,101 @@
+import calc_chauffage from './9_chauffage.js';
+
+describe('Recherche de bugs dans le calcul de la consommation', () => {
+  test('calcul de la consommation de chauffage pour 2475E2510509B', () => {
+    const Sh = 22.76; // Surface habitable
+    const zc_id = 1; // Zone climatique
+    const ca_id = 1; // Classe d'altitude
+    const inertie_id = 4; // Inertie du batiment
+    const map_id = 2; // Méthode d'application du DPE (enum_methode_application_dpe_log_id)
+    const bch = 7998.446638037776; // Besoin en chauffage
+    const bch_dep = 9714.925476692655; // Besoin en chauffage dépensier
+    const GV = 150.18558182166174; // Déperdition de l'enveloppe
+    const hsp = 2.5700000000000003; // Hauteur sous plafond
+    const ac = 1; // Période de construction
+
+    const ch = {
+      donnee_entree: {
+        description:
+          'Convecteur électrique NFC, NF** et NF*** avec programmateur pièce par pièce (système individuel)',
+        reference: '2024_07_11_09_30_44_722898200737262',
+        surface_chauffee: 22.76,
+        rdim: 1,
+        nombre_niveau_installation_ch: 1,
+        enum_cfg_installation_ch_id: '1',
+        enum_type_installation_id: '1',
+        enum_methode_calcul_conso_id: '1'
+      },
+      emetteur_chauffage_collection: {
+        emetteur_chauffage: [
+          {
+            donnee_entree: {
+              description: '',
+              reference: 'Emetteur:2024_07_11_09_30_44_722898200737262#1',
+              surface_chauffee: 22.76,
+              tv_rendement_emission_id: 1,
+              tv_rendement_distribution_ch_id: 1,
+              tv_rendement_regulation_id: 1,
+              enum_type_emission_distribution_id: '1',
+              tv_intermittence_id: 138,
+              reseau_distribution_isole: 0,
+              enum_equipement_intermittence_id: '4',
+              enum_type_regulation_id: '2',
+              enum_periode_installation_emetteur_id: '1',
+              enum_type_chauffage_id: '1',
+              enum_temp_distribution_ch_id: '1',
+              enum_lien_generateur_emetteur_id: '1'
+            }
+          }
+        ]
+      },
+      generateur_chauffage_collection: {
+        generateur_chauffage: [
+          {
+            donnee_entree: {
+              description: 'Electrique - Convecteur électrique NFC, NF** et NF***',
+              reference: 'Generateur:2024_07_11_09_30_44_722898200737262#1',
+              reference_generateur_mixte: '',
+              ref_produit_generateur_ch: 'Sans Objet',
+              enum_type_generateur_ch_id: '98',
+              enum_usage_generateur_id: '1',
+              enum_type_energie_id: '1',
+              position_volume_chauffe: 1,
+              tv_rendement_generation_id: 29,
+              identifiant_reseau_chaleur: '',
+              enum_methode_saisie_carac_sys_id: '1',
+              enum_lien_generateur_emetteur_id: '1'
+            }
+          }
+        ]
+      }
+    };
+
+    calc_chauffage(ch, ca_id, zc_id, inertie_id, map_id, bch, bch_dep, GV, Sh, hsp, ac);
+
+    expect(ch.donnee_intermediaire).toStrictEqual({
+      besoin_ch: 7998.446638037776,
+      besoin_ch_depensier: 9714.925476692655,
+      conso_ch: 6322.706769289927,
+      conso_ch_depensier: 7679.569278179351
+    });
+    expect(
+      ch.emetteur_chauffage_collection.emetteur_chauffage[0].donnee_intermediaire
+    ).toStrictEqual({
+      i0: 0.86,
+      rendement_distribution: 1,
+      rendement_emission: 0.95,
+      rendement_regulation: 0.99
+    });
+    expect(
+      ch.generateur_chauffage_collection.generateur_chauffage[0].donnee_intermediaire
+    ).toStrictEqual({
+      conso_auxiliaire_generation_ch: 0,
+      conso_auxiliaire_generation_ch_depensier: 0,
+      conso_ch: 6322.706769289927,
+      conso_ch_depensier: 7679.569278179351,
+      rendement_generation: 1,
+      rg: 1,
+      rg_dep: 1
+    });
+  });
+});

package/9_conso_ch.js

@@ -0,0 +1,95 @@
+import { requestInputID } from './utils.js';
+import { rendement_emission } from './9_emetteur_ch.js';
+import { calc_intermittence } from './8_intermittence.js';
+
+function coef_ch(Fch) {
+  return {
+    'installation de chauffage simple': {
+      0: 1,
+      1: 1
+    },
+    'installation de chauffage avec chauffage solaire': {
+      0: 1 - Fch
+    },
+    'installation de chauffage avec insert ou poêle bois en appoint': {
+      0: 0.75,
+      1: 0.25
+    },
+    // todo: 2 sdb
+    'installation de chauffage par insert, poêle bois (ou biomasse) avec un chauffage électrique dans la salle de bain':
+      {
+        0: 0.9,
+        1: 0.1
+      },
+    'installation de chauffage avec en appoint un insert ou poêle bois et un chauffage électrique dans la salle de bain (différent du chauffage principal)':
+      {
+        0: 0.75 * 0.9,
+        1: 0.25 * 0.9,
+        2: 0.1
+      },
+    'installation de chauffage avec une chaudière ou une pac en relève d’une chaudière bois': {
+      0: 0.75,
+      1: 0.25
+    },
+    'installation de chauffage avec chauffage solaire et insert ou poêle bois en appoint': {
+      0: 0.75 * (1 - Fch),
+      1: 0.25 * (1 - Fch)
+    },
+    'installation de chauffage avec chaudière en relève de pac': {
+      0: 0.8,
+      1: 0.2
+    },
+    'installation de chauffage avec chaudière en relève de pac avec insert ou poêle bois en appoint':
+      {
+        0: 0.8 * 0.75,
+        1: 0.2 * 0.75,
+        2: 0.25
+      },
+    'installation de chauffage collectif avec base + appoint': {
+      0: 0.75 * (1 - Fch),
+      1: 0.25 * (1 - Fch)
+    },
+    'convecteurs bi-jonction': {
+      0: 0.6,
+      1: 0.4
+    }
+  };
+}
+
+export function conso_ch(di, de, du, _pos, cfg_ch, em_list, GV, Sh, hsp, bch, bch_dep) {
+  const gen_lge_id = requestInputID(de, du, 'lien_generateur_emetteur');
+  const coef = coef_ch(de.fch || 0.5)[cfg_ch][_pos] || 1;
+  let em_filt;
+
+  if (em_list.length === 1) {
+    em_filt = em_list;
+  } else {
+    em_filt = em_list.filter(
+      (em) => em.donnee_entree.enum_lien_generateur_emetteur_id === gen_lge_id
+    );
+  }
+
+  const hasMultipleEmetteur = em_filt.length > 1;
+
+  const emetteur_eq = em_filt.reduce((acc, em) => {
+    const int = calc_intermittence(GV, Sh, hsp, em.donnee_intermediaire.i0);
+    const r_em = rendement_emission(em);
+
+    /**
+     * 9.1.3 Installation avec plusieurs émissions pour un même générateur
+     * La part de la consommation traitée par chaque émetteur est proratisé par le ratio des surfaces habitables.
+     * @type {number|number}
+     */
+    const ratio_s = hasMultipleEmetteur
+      ? em.donnee_entree.surface_chauffee / de.surface_chauffee
+      : 1;
+
+    const Ich = 1 / r_em;
+    return acc + ratio_s * int * Ich;
+  }, 0);
+
+  const Ich = emetteur_eq / di.rg;
+  const Ich_dep = emetteur_eq / di.rg_dep;
+  di.conso_ch = coef * Ich * bch;
+  di.conso_ch_depensier = coef * Ich_dep * bch_dep;
+}