Пример #1
0
  private Hashtable<Spezies, Double> calc_dmj_dp() {
    Hashtable<Spezies, Double> massenBruchHash_Zone = gg_Zone.get_speziesMassenBruecheDetail();
    Hashtable<Spezies, Double> dmj_dp = new Hashtable<Spezies, Double>();

    Enumeration<Spezies> e = massenBruchHash_Zone.keys();
    Spezies spez;

    if (loeseChemGleichgwicht == true && burns == true) {
      double delta_p = 5000;
      GasGemisch rauchgas =
          new GasGemisch(
              GG_SOLVER.get_GG_molenBrueche(p_Zone + delta_p, T_Zone, gg_Zone),
              "rauchgas_temp_" + ID);
      Hashtable<Spezies, Double> rgSpezMassenBrueche = rauchgas.get_speziesMassenBruecheDetail();

      while (e.hasMoreElements()) {
        spez = e.nextElement();
        dmj_dp.put(
            spez,
            (m_Zone * (rgSpezMassenBrueche.get(spez) - massenBruchHash_Zone.get(spez))) / delta_p);
      }

    } else {

      while (e.hasMoreElements()) {
        spez = e.nextElement();
        dmj_dp.put(spez, 0.0);
      }
    }

    return dmj_dp;
  }
Пример #2
0
 /**
  * Mit dieser Methode kann ein Massenelement eines Gases mit der Temperatur T_m der Zone mit der
  * Temperatur T_Zone zugemischt werden. Die sich ergebende Mischungstemperatur wird für
  * VERÄNDERLICHE CVs und ein adiabates System nach dem ersten HS berechnet. Das Volumen der Zone
  * bleibt Konstant --> das zugemischte Massenelement leistet Verschiebearbeit!
  *
  * @param m_Zu --> zugemischtes Massenelement [kg]
  * @param T_Zu --> Temperatur des zugemischten Massenelements [K]
  * @param s_Zu --> Zusammensetzung des Massenelements
  */
 public void massenElementZumischenKonstVol(double m_Zu, double T_Zu, Spezies s_Zu) {
   // adiabate Mischungstemperatur nach dem ersten Hauptsatz
   double Tm = 278.15, U1, U_zu, U1m, U_zum, Cvm1, Cv_zum, F, dF, Tm_buffer;
   U1 = gg_Zone.get_u_mass(T_Zone) * m_Zone;
   U_zu = s_Zu.get_u_mass(T_Zu) * m_Zu;
   int idx = 0;
   // Newtonverfahren für F(Tm)= U1(T1)+U2(T2)-U1(Tm)-U2(Tm)=0
   // mit dF/dT=-m1*Cv1(Tm)-m2*Cv2(Tm)
   do {
     Tm_buffer = Tm;
     U1m = gg_Zone.get_u_mass(Tm) * m_Zone;
     U_zum = s_Zu.get_u_mass(Tm) * m_Zu;
     Cvm1 = gg_Zone.get_cv_mass(Tm) * m_Zone;
     Cv_zum = s_Zu.get_cv_mass(Tm) * m_Zu;
     F = U1 + U_zu - U1m - U_zum;
     dF = -Cvm1 - Cv_zum;
     Tm = Tm - F / dF;
     idx++;
   } while (idx < 1000 && Math.abs(Tm - Tm_buffer) > 0.1);
   if (idx >= 1000) {
     try {
       throw new MiscException("t");
     } catch (MiscException me) {
       me.log_Warning("Mangelnde Konvergenz bei der Berechnung der Mischungstemperatur");
     }
   }
   T_Zone = Tm;
   // dieser Aufruf muss nach der Temperaturberechnung erfolgen da sich hier die Masse der Zone
   // ändert
   massenElementZumischen(m_Zu, s_Zu);
 }
Пример #3
0
 /**
  * Mit dieser Methode kann ein Massenelement eines Gases mit der Temperatur T_m der Zone mit der
  * Temperatur T_Zone zugemischt werden. Die sich ergebende Mischungstemperatur wird für
  * VERÄNDERLICHE CVs und ein adiabates System nach dem ersten HS berechnet. Das Volumen der Zone
  * bleibt Konstant --> das zugemischte Massenelement leistet Verschiebearbeit!
  *
  * @param m_Zu --> zugemischtes Massenelement [kg]
  * @param T_Zu --> Temperatur des zugemischten Massenelements [K]
  * @param s_Zu --> Zusammensetzung des Massenelements
  */
 public void massenElementZumischenKonstVolEinfacher(double m_Zu, double T_Zu, Spezies s_Zu) {
   // Nach dem 1HS gilt: dm*h=dU
   double dH = m_Zu * s_Zu.get_h_mass(T_Zu);
   double U_Zone_mass_temp = m_Zone * gg_Zone.get_u_mass(T_Zone) + dH;
   // dieser Aufruf muss vor der Temperaturberechnung erfolgen da sich hier die Masse und die
   // Zusammensetzung
   // der Zone ändert
   massenElementZumischen(m_Zu, s_Zu);
   // Berechnen der Mischungstemperatur
   U_Zone_mass_temp = U_Zone_mass_temp / m_Zone;
   T_Zone = gg_Zone.get_T4u_mass(U_Zone_mass_temp);
 }
Пример #4
0
  private double[] calc_dudT_dRdT() {
    double[] du_dt_dR_dT = new double[2];
    double delta_T = 1;
    // if-Abfrage nur dann sinnvoll wenn Einfriertemperatur unter 1000K
    // oder wenn Luft als Rauchgas mit sehr gro§em Lambda modelliert wird
    if (T_Zone > 1000 && T_Zone <= 1001) {
      delta_T = -1;
    }
    // Berechnen der kalorischen Daten für Rauchgas
    if (loeseChemGleichgwicht == true && T_Zone >= GG_SOLVER.get_T_Freeze()) {
      GasGemisch zoneGG =
          new GasGemisch(GG_SOLVER.get_GG_molenBrueche(p_Zone, T_Zone, gg_Zone), "zoneGG" + ID);
      GasGemisch rauchgas =
          new GasGemisch(
              GG_SOLVER.get_GG_molenBrueche(p_Zone, T_Zone - delta_T, gg_Zone),
              "rauchgasZone" + ID);

      // ist zwar umstaendlich aber ich hab schon alles mit Java erlebt :-(
      double a = (zoneGG.get_u_mass(T_Zone) - rauchgas.get_u_mass(T_Zone - delta_T));
      du_dt_dR_dT[0] = a / delta_T;

      a = (zoneGG.get_R() - rauchgas.get_R());
      du_dt_dR_dT[1] = a / delta_T;
    } else {
      du_dt_dR_dT[0] =
          gg_Zone.get_cv_mass(T_Zone); // hier ändert sich die Zusammensetzung des Gases nicht mehr
      du_dt_dR_dT[1] = 0;
    }
    return du_dt_dR_dT;
  }
Пример #5
0
  /**
   * Mit dieser Methode kann eine Zone mit einer anderen gemischt werden. Das Volumen der beiden
   * Zonen wird addiert entsprechent wird keine Verschiebearbeit geleistet. Temperaur, Masse und
   * Zusammensetzung werden entsprechend {@link massenElementZumischen} berechnet
   *
   * @param Zone zone2 --> die Zone die zugemischt werden soll
   */
  public void zoneZumischen(Zone zone2) {

    double m_Zu = zone2.get_m();
    double T_Zu = zone2.get_T();
    double V_Zu = zone2.get_V();
    Spezies s_Zu = zone2.get_ggZone();

    // Nach dem 1HS gilt: -pdV+dm*h=dU --> -p*deltaV+deltaH=deltaU=U2-U1 --> U2=U1+deltaH-p*deltaV
    double H_zu = m_Zu * s_Zu.get_h_mass(T_Zu);
    // deltaV=V_zu  außerdem ist der Druck in beiden Zonen identisch
    double U2_Zone_mass = m_Zone * gg_Zone.get_u_mass(T_Zone) + H_zu - p_Zone * V_Zu;
    // dieser Aufruf muss vor der Temperaturberechnung erfolgen da sich hier die Masse und die
    // Zusammensetzung
    // der Zone ändert
    massenElementZumischen(m_Zu, s_Zu);
    // Berechnen der Mischungstemperatur
    T_Zone = gg_Zone.get_T4u_mass(U2_Zone_mass / m_Zone);
  }
Пример #6
0
 public double[] get_MoleFractions() {
   Hashtable<Spezies, Double> mf = gg_Zone.get_speziesMolenBruecheDetailToIntegrate();
   double x[] = new double[CP.SPEZIES_FABRIK.get_NmbrOfAllSpez()];
   Enumeration<Spezies> e = mf.keys();
   Spezies spez;
   while (e.hasMoreElements()) {
     spez = e.nextElement();
     x[CP.SPEZIES_FABRIK.get_indexOf(spez)] = mf.get(spez);
   }
   return x;
 }
Пример #7
0
 public double[] get_mi() {
   // TODO mach mich mi[] zur Klassenvariablen, dann duerfte es etwas schneller rechnen
   Hashtable<Spezies, Double> mb = gg_Zone.get_speziesMassenBruecheDetailToIntegrate();
   double mi[] = new double[CP.SPEZIES_FABRIK.get_NmbrOfAllSpez()];
   Enumeration<Spezies> e = mb.keys();
   Spezies spez;
   while (e.hasMoreElements()) {
     spez = e.nextElement();
     mi[CP.SPEZIES_FABRIK.get_indexOf(spez)] = mb.get(spez) * m_Zone;
   }
   return mi;
 }
Пример #8
0
  private double[] calc_dudp_dRdp() {
    double[] du_dp_dR_dp = new double[2];

    double delta_p = 5000;
    // Berechnen der kalorischen Daten für Rauchgas
    if (loeseChemGleichgwicht == true && T_Zone >= GG_SOLVER.get_T_Freeze()) {
      GasGemisch zoneGG =
          new GasGemisch(GG_SOLVER.get_GG_molenBrueche(p_Zone, T_Zone, gg_Zone), "zoneGG" + ID);
      GasGemisch rauchgas =
          new GasGemisch(
              GG_SOLVER.get_GG_molenBrueche(p_Zone + delta_p, T_Zone, gg_Zone),
              "rauchgasZone" + ID);

      double a = (rauchgas.get_u_mass(T_Zone) - zoneGG.get_u_mass(T_Zone));
      du_dp_dR_dp[0] = a / delta_p;

      a = (rauchgas.get_R() - zoneGG.get_R());
      du_dp_dR_dp[1] = a / delta_p;

    } else {
      du_dp_dR_dp[0] = 0;
      du_dp_dR_dp[1] = 0;
    }
    return du_dp_dR_dp;
  }
Пример #9
0
  private Hashtable<Spezies, Double> calc_dmj_dT() {
    Hashtable<Spezies, Double> massenBruchHash_Zone = gg_Zone.get_speziesMassenBruecheDetail();
    Hashtable<Spezies, Double> dmj_dT = new Hashtable<Spezies, Double>();
    Enumeration<Spezies> e = massenBruchHash_Zone.keys();
    Spezies spez;

    // Berechnen der kalorischen Daten für Rauchgas
    if (loeseChemGleichgwicht == true && burns == true) {
      double delta_T = 1;
      // if-Abfrage nur dann sinnvoll wenn Einfriertemperatur unter 1000K
      // oder wenn Luft als Rauchgas mit sehr gro§em Lambda modelliert wird
      if (T_Zone > 1000 && T_Zone <= 1001) {
        delta_T = -1;
      }

      GasGemisch rauchgas =
          new GasGemisch(
              GG_SOLVER.get_GG_molenBrueche(p_Zone, T_Zone - delta_T, gg_Zone), "rauchgas");
      Hashtable<Spezies, Double> rgSpezMassenBrueche = rauchgas.get_speziesMassenBruecheDetail();

      while (e.hasMoreElements()) {
        spez = e.nextElement();
        dmj_dT.put(
            spez,
            (m_Zone * (massenBruchHash_Zone.get(spez) - rgSpezMassenBrueche.get(spez)) / delta_T));
      }

    } else {

      while (e.hasMoreElements()) {
        spez = e.nextElement();
        dmj_dT.put(spez, 0.0);
      }
    }
    return dmj_dT;
  }
Пример #10
0
  public static Zone zonenMischen(CasePara cp, Zone z1, Zone z2, boolean burns, int ID) {
    double p = z1.get_p(); // fuer beide Zonen gleich
    double m1 = z1.get_m();
    double T1 = z1.get_T();
    double V1 = z1.get_V();
    Spezies s1 = z1.get_ggZone();

    double m2 = z2.get_m();
    double T2 = z2.get_T();
    double V2 = z2.get_V();
    Spezies s2 = z2.get_ggZone();

    // Erzeugen einer Spezies die durch die Mischung der beiden ZonenSpezies entsteht
    GasGemisch s0 = new GasGemisch("zonenGemisch");
    Hashtable<Spezies, Double> ht0 = new Hashtable<Spezies, Double>(3);
    double mGes = m1 + m2;
    ht0.put(s1, m1 / mGes);
    ht0.put(s2, m2 / mGes);
    s0.set_Gasmischung_massenBruch(ht0);
    // Berechnung der Mischungstemperatur
    // 1.HS --> Kontrollraumgrenze wird um beide Zonen gelegt
    // Dann gilt: Umisch=U1+U2;
    double Umisch = m1 * s1.get_u_mass(T1) + m2 * s2.get_u_mass(T2);
    // Berechnung der Temperatur bei der die innere Energie von Spezies s0=Umisch ist.
    // Dies entspricht der Mischungstemperatur
    double T0 = s0.get_T4u_mass(Umisch / mGes);
    // Das Volumen beider Zonen wird addiert
    double V0 = V1 + V2;

    // Test ob alles stimmt
    double pV = z1.get_p() * V0;
    double mRT = mGes * T0 * s0.get_R();
    double T =
        (m1 * s1.get_cv_mass(T1) * T1 + m2 * s2.get_cv_mass(T2) * T2)
            / (m1 * s1.get_cv_mass(T1) + m2 * s2.get_cv_mass(T2));
    double deltaT = (T - T0) / T0 * 100;
    double T3 = pV / mGes / s0.get_R();
    double mRT2 = mGes * T * s0.get_R();
    double mRT3 = mGes * T3 * s0.get_R();
    // Fazit ueber die Gasgleichung geht es auch und zwar viel einfacher und genauer!

    Zone z0 = new Zone(cp, p, V0, T0, mGes, s0, burns, ID);
    return z0;
  }
Пример #11
0
  /**
   * Mit dieser Methode kann ein Massenelement eines Gases der Zone zugemischt werden. --> es wird
   * nur die Zusammensetzung und die Gesamtmasse der Zone geändert nicht aber die Temperatur
   *
   * @param idx --> Rechenindex
   * @param mZu --> zugemischtes Massenelement [kg]
   * @param sZu --> SpeziesObjekt des Massenelements
   */
  public void massenElementZumischen(double mZu, Spezies sZu) {

    Hashtable<Spezies, Double> massenBruchHash_Zone = gg_Zone.get_speziesMassenBruecheDetail();

    Hashtable<Spezies, Double> massenBruchHashNeu = new Hashtable<Spezies, Double>(4);
    double m_Zone_neu = mZu + m_Zone;
    if (m_Zone_neu < 0) {
      try {
        throw new MiscException(
            "ERROR: Aus Zone " + ID + " wird mehr Masse entnommen als vorhanden ist");

      } catch (MiscException miscE) {
        miscE.stopBremo();
      }
    }
    Hashtable<Spezies, Double> einzelMassenHashIn = new Hashtable<Spezies, Double>();
    // einzelmassen der zugefuehrten Spezies bestimmen
    if (mZu < 0) {
      if (sZu.isGasGemisch()) {
        einzelMassenHashIn = // enthaelt negative Werte!!
            this.get_einzelMassenHash(((GasGemisch) sZu).get_speziesMassenBruecheDetail(), mZu);
      } else {
        einzelMassenHashIn.put(sZu, mZu);
      }
      // einzelmassen der Zone bestimmen
      Hashtable<Spezies, Double> einzelMassenHashZone =
          this.get_einzelMassenHash(massenBruchHash_Zone, m_Zone);

      // Addieren der einzelnen Massen von Zone und sZu
      Enumeration<Spezies> e = einzelMassenHashZone.keys();
      Spezies spez;
      double mTemp;

      while (e.hasMoreElements()) {
        spez = e.nextElement();
        mTemp = einzelMassenHashZone.get(spez);

        if (einzelMassenHashIn.containsKey(spez))
          mTemp = mTemp + einzelMassenHashIn.get(spez); // Wert in Hashtable ist negativ

        double massenBruch = mTemp / m_Zone_neu;
        massenBruchHashNeu.put(spez, massenBruch);

        if (massenBruch
            < 0) { // kommt vor wenn eine Masse entnommen werden soll die gar nicht in der Zone ist!
          try {
            throw new MiscException(
                "ERROR: In einer Zone tritt ein negativer Massenbruch auf \n"
                    + spez.get_name()
                    + ": "
                    + massenBruch);

          } catch (MiscException miscE) {
            massenBruchHashNeu.put(spez, 0d);
            //							miscE.stopBremo(); //TODO check mich
          }
        }
      }
    } else {
      massenBruchHashNeu.put(gg_Zone, m_Zone / m_Zone_neu);
      if (massenBruchHashNeu.containsKey(
          sZu)) { // ==True wenn der Zone eine Masse von sich selbst zugefuehrt wird
        massenBruchHashNeu.put(gg_Zone, 1D);
      } else massenBruchHashNeu.put(sZu, mZu / m_Zone_neu);
    }

    Hashtable<Spezies, Double> molenBruchHash = new Hashtable<Spezies, Double>(4);
    molenBruchHash = GasGemisch.Gasmischer.massenBruch2molenBruch(massenBruchHashNeu);

    GasGemisch gasGemisch = new GasGemisch(molenBruchHash, "");
    // würde man hier die
    gg_Zone.set_Gasmischung_molenBruch(gasGemisch.get_speziesMolenBruecheDetail());
    m_Zone = m_Zone_neu;
  }
Пример #12
0
 public double get_bT() {
   double bT = -m_Zone * gg_Zone.get_R() - m_Zone * T_Zone * dRdT;
   return bT;
 }
Пример #13
0
  /**
   * Berechnet den Massenverlust bei entnahme einer BESTIMMTEN Spezies Dieser Aufruf sollte nur vom
   * einzonigen Modell benoetigt werden
   *
   * @param dm_a [kg] muss groeser Null sein
   * @throws NegativeMassException
   */
  public void set_dm_aus(double dm_a, Spezies spez) throws NegativeMassException {
    if (dm_a < 0)
      throw new IllegalArgumentException(
          "\"set_dm_aus()\"  wurde fuer Zone "
              + ID
              + " mit einem negativen Massenstrom aufgerufen");
    if (dm_a > 0) {
      // Hinzufügen der Einzelmassen der entnommenen Spezies zur Hashtable dmj
      Hashtable<Spezies, Double> einzelMassenHash_aus = new Hashtable<Spezies, Double>();
      if (spez.isGasGemisch() && CP.SPEZIES_FABRIK.isToIntegrate(spez) == false) {
        Hashtable<Spezies, Double> massenBruchHash_aus =
            ((GasGemisch) spez).get_speziesMassenBruecheDetailToIntegrate();

        // Erstellen der Hashtable mit den Einzelmassen der Grundspezies (CO2, O2, usw....
        einzelMassenHash_aus = this.get_einzelMassenHash(massenBruchHash_aus, dm_a);

      } else {
        einzelMassenHash_aus.put(spez, dm_a);
      }
      // Hinzufügen der Änderung der Grundspeziesmassen
      Enumeration<Spezies> e = einzelMassenHash_aus.keys();
      Spezies spez1;
      while (e.hasMoreElements()) {
        spez1 = e.nextElement();
        if (!CP.SPEZIES_FABRIK.isToIntegrate(spez1)) {
          try {
            throw new BirdBrainedProgrammerException(
                "Aus Zone "
                    + this.ID
                    + " wurde eine Spezies "
                    + "("
                    + spez1.get_name()
                    + ") entnommen die nicht integriert werden soll"
                    + " (\"isToIntegrate()\" liefert false)");
          } catch (BirdBrainedProgrammerException bbpe) {
            bbpe.stopBremo();
          }
        }
        if (!gg_Zone.get_speziesMassenBruecheDetailToIntegrate().containsKey(spez1)
            && einzelMassenHash_aus.get(spez1)
                > 0) // Es kann sein, dass Spezies mit einem Massenbruch von 0 in der Hashtable
                     // stehen
        throw new NegativeMassException(
              "Aus Zone "
                  + this.ID
                  + " soll eine Spezies"
                  + "("
                  + spez1.get_name()
                  + ") entnommen werden "
                  + "die in der Zone nicht vorhanden ist");
      }

      e = null;
      e = einzelMassenHash_aus.keys();
      spez1 = null;
      while (e.hasMoreElements()) {
        spez1 = e.nextElement();
        if (this.dm_aus.containsKey(spez1)) {
          this.dm_aus.put(spez1, einzelMassenHash_aus.get(spez1) + dm_aus.get(spez1));
        } else {
          this.dm_aus.put(spez1, einzelMassenHash_aus.get(spez1));
        }
      }

      //			e=null;
      //			e=einzelMassenHash_aus.keys();
      //			spez1=null;
      //			double dmCheck=0;
      //			while(e.hasMoreElements()){
      //				spez1=e.nextElement();
      //				dmCheck+=einzelMassenHash_aus.get(spez1);
      //			}
      //			if(dmCheck!=dm_a)
      //				System.out.println("Zone.set_dm_aus: dm_a !=dmCheck " + ((dmCheck-dm_a)) );

      // das hier muss nach der NegativeMassException stehen!! Wenn die Exception geworfen
      // wird und dieser Code ist schon ausgefuehrt wuerde die Energie aber nicht dei Masse
      // beruecksichtigt
      double dU = spez.get_u_mass(T_Zone) * dm_a;
      sumdU = sumdU - dU;
      double dH = spez.get_h_mass(T_Zone) * dm_a;
      sumdH = sumdH - dH;
    }
  }
Пример #14
0
  public void set_p_V_T_mi(
      double p_init, double V_init, double T_init, Hashtable<Spezies, Double> mi) {
    this.p_Zone = p_init;
    this.V_Zone = V_init;
    this.T_Zone = T_init;
    if (p_Zone <= 0
        || V_Zone <= 0
        || T_Zone <= 0
        || (((Double) p_Zone).isNaN())
        || (((Double) T_Zone).isNaN())
        || (((Double) V_Zone).isNaN())) {
      try {
        throw new NegativeMassException(
            "Falsche Werte in Zone "
                + this.ID
                + "\n V= "
                + V_Zone
                + "\n p= "
                + p_Zone
                + "\n T= "
                + T_Zone);
      } catch (NegativeMassException nmE) {
        nmE.log_Warning();
      }
    }
    if (CP.get_aktuelle_Rechenzeit() > CP.SYS.WRITE_INTERVAL_SEC) {
      // System.out.println("xsi: "+this.get_xsi());
      //	System.out.println("eta: " +this.get_eta());
    }

    if (T_Zone > T_max_Zone) T_max_Zone = T_Zone;

    // Überprüfen ob die Berechnung der Dissoziation Sinn macht
    if (T_max_Zone >= GG_SOLVER.get_T_Freeze() && burns) loeseChemGleichgwicht = true;

    // Berechnen der Gesamtmasse
    double mGes = 0.0D;
    double mi_;
    Enumeration<Spezies> e = mi.keys();
    Spezies spez;
    while (e.hasMoreElements()) {
      spez = e.nextElement();
      mi_ = mi.get(spez);
      if (mi_ < 0) {
        if (Math.abs(mi_) >= 0.00001 * CP.SYS.MINIMALE_ZONENMASSE) {
          try {
            throw new NegativeMassException(
                spez.get_name()
                    + " in Zone "
                    + this.ID
                    + " hatte eine negative Masse ("
                    + mi_
                    + "kg). Die Masse wurde auf null gesetzt!");
          } catch (NegativeMassException nmE) {
            nmE.log_Warning();
          }
        }
        mi_ = 0;
        mi.put(spez, mi_);
      }
      mGes = mGes + mi_;
    }
    spez = null;
    e = null;
    this.m_Zone = mGes;

    if (this.m_Zone > 0) {
      // Berechnen der Massenbrueche
      Hashtable<Spezies, Double> massenBrueche = new Hashtable<Spezies, Double>();
      e = mi.keys();
      double xi;
      while (e.hasMoreElements()) {
        spez = e.nextElement();
        xi = mi.get(spez) / mGes;
        massenBrueche.put(spez, xi);
      }
      gg_Zone.set_Gasmischung_massenBruch(massenBrueche);
      // Berechnen der Dissoziation
      if (loeseChemGleichgwicht == true)
        gg_Zone.set_Gasmischung_molenBruch(GG_SOLVER.get_GG_molenBrueche(p_Zone, T_Zone, gg_Zone));

    } else {
      Hashtable<Spezies, Double> massenBrueche = new Hashtable<Spezies, Double>();
      Spezies co2 = CP.SPEZIES_FABRIK.get_spezCO2();
      massenBrueche.put(co2, 1D); // TODO irgendwie uncool
      gg_Zone.set_Gasmischung_massenBruch(massenBrueche);
    }

    double pV = p_Zone * V_Zone;
    double mRT = m_Zone * gg_Zone.get_R() * T_Zone;
    // This accounts for inaccuracy during the integration process.
    // The differences in temperature should be very very small!!
    //			if(CP.get_time()>CP.convert_KW2SEC(5.2))
    //				System.out.println();
    if (pV != mRT && m_Zone > 0 && V_Zone > 0) {
      //	System.out.println("(pV-mRT)/pV: "+(pV-mRT)/pV*100);
      double T_temp = pV / m_Zone / gg_Zone.get_R();
      this.T_Zone = T_temp;
      double chkPv = pV - m_Zone * gg_Zone.get_R() * T_temp;
      double chkT = Math.abs((T_Zone - T_temp) / T_temp);
      if (chkT > 7.555e-3 && m_Zone > CP.SYS.MINIMALE_ZONENMASSE) {
        try {
          throw new MiscException(
              "The temperature of  zone  "
                  + this.ID
                  + " had to be adjusted too much! "
                  + "chkT= "
                  + chkT);
        } catch (MiscException me) {
          me.log_Warning();
        }
      }
    }

    dmj_ein = new Hashtable<Spezies, Double>();
    dm_aus = new Hashtable<Spezies, Double>();
    sumdQ = 0;
    sumdH = 0;
    sumdU = 0;
    dudp = this.calc_dudp_dRdp()[0];
    dRdp = this.calc_dudp_dRdp()[1];
    dudT = this.calc_dudT_dRdT()[0];
    dRdT = this.calc_dudT_dRdT()[1];
  }