/// <summary> /// Measure TS content inside a digester /// /// type 7 /// </summary> /// <param name="x">ADM state vector</param> /// <param name="myPlant"></param> /// <param name="mySubstrates">list of substrates</param> /// <param name="mySensors"></param> /// <param name="Q"> /// substrate feed and recirculation sludge going into the digester /// first values are Q for substrates, then pumped sludge going into digester /// dimension: always number of substrates + number of digesters /// </param> /// <param name="par">not used</param> /// <returns></returns> override protected physValue[] doMeasurement(double[] x, biogas.plant myPlant, biogas.substrates mySubstrates, biogas.sensors mySensors, double[] Q, params double[] par) { // TODO // so erweitern, dass auch TS in Fermenter Input gemessen werden kann // evtl. substrate_sensor nutzen, kann aber auch direkt über mySubstrates gemessen werden physValue[] values = new physValue[1]; // number of substrates int n_substrate = mySubstrates.getNumSubstrates(); if (Q.Length < n_substrate) { throw new exception(String.Format( "Q.Length < n_substrate: {0} < {1}!", Q.Length, n_substrate)); } double[] Qsubstrates = new double[n_substrate]; // volumeflow for substrates for (int isubstrate = 0; isubstrate < n_substrate; isubstrate++) { Qsubstrates[isubstrate] = Q[isubstrate]; } // biogas.substrates substrates_or_sludge; // List <double> Q_s_or_s = new List <double>(); // if no substrate is going into the fermenter we have to take the // TS from the digesters sludge going into this digester to calculate a // sludge. If there is substrate going into the digester we ignore // recirculation sludge, because the COD content inside the digester is // already influenced by the recirculated COD, so a high recirculation leads // to a reduction of the COD content inside the digester and then also to a // TS content reduction. if (math.sum(Qsubstrates) == 0) { substrates_or_sludge = new biogas.substrates(); int ifermenter = 0; // for (int iflux = n_substrate; iflux < Q.Length; iflux++) { string digester_id = myPlant.getDigesterID(ifermenter + 1); double TS_digester = 0; try { mySensors.getCurrentMeasurementD("TS_" + digester_id + "_3", out TS_digester); if (TS_digester < double.Epsilon) // vermutlich wurde noch nicht gemessen? { TS_digester = 11; } } catch // TODO: wann passiert das??, sollte eigentlich nicht passieren { TS_digester = 11; } substrates_or_sludge.addSubstrate( new biogas.sludge(mySubstrates, math.ones(n_substrate), TS_digester)); ifermenter = ifermenter + 1; } // for (int isubstrate = n_substrate; isubstrate < Q.Length; isubstrate++) { Q_s_or_s.Add(Q[isubstrate]); } } else { substrates_or_sludge = mySubstrates; for (int isubstrate = 0; isubstrate < Qsubstrates.Length; isubstrate++) { Q_s_or_s.Add(Qsubstrates[isubstrate]); } } if (id.EndsWith("3")) // out sensor { // calc TS inside digester due to the substrates or the sludge if digester is not fed values[0] = biogas.digester.calcTS(x, substrates_or_sludge, Q_s_or_s.ToArray()); } else if (id.EndsWith("2")) // in sensor { // TODO // ändern, da momentan entweder substrate oder sludge genommen wird und nicht beides // selbe Problem wie bei OLR_sensor denke ich substrates_or_sludge.get_weighted_mean_of(Q_s_or_s.ToArray(), "TS", out values[0]); values[0].Symbol = "TS"; } else { throw new exception(String.Format("id of TS sensor not valid: {0}", id)); } // values[0].Label = "total solids"; return(values); }
/// <summary> /// Measure OLR of a digester /// /// type 7 /// /// TODO: größtenteils identisch mit TS_sensor, man könnte was zusammen legen /// </summary> /// <param name="x">ADM state vector</param> /// <param name="myPlant"></param> /// <param name="mySubstrates">list of substrates</param> /// <param name="mySensors"></param> /// <param name="Q"> /// substrate feed and recirculation sludge going into fermenter in m³/d /// first values are Q for substrates, then pumped sludge going into digester /// dimension: always number of substrates + number of digesters /// </param> /// <param name="par"></param> /// <returns></returns> override protected physValue[] doMeasurement(double[] x, biogas.plant myPlant, biogas.substrates mySubstrates, biogas.sensors mySensors, double[] Q, params double[] par) { physValue[] values = new physValue[1]; // TODO: so abändern, dass die aufgerufene methode calcOLR immer // feed und sludge übergeben werden. nicht oder // number of substrates int n_substrate = mySubstrates.getNumSubstrates(); double Qsum = math.sum(Q); if (Q.Length < n_substrate) { throw new exception(String.Format( "Q.Length < n_substrate: {0} < {1}!", Q.Length, n_substrate)); } double[] Qsubstrates = new double[n_substrate]; // volumeflow for substrates for (int isubstrate = 0; isubstrate < n_substrate; isubstrate++) { Qsubstrates[isubstrate] = Q[isubstrate]; } // biogas.substrates substrates_or_sludge; // List <double> Q_s_or_s = new List <double>(); // if no substrate is going into the fermenter we have to take the // TS from the digesters sludge going into this digester to calculate a // sludge. If there is substrate going into the digester we ignore // recirculation sludge, because the COD content inside the digester is // already influenced by the recirculated COD, so a high recirculation leads // to a reduction of the COD content inside the digester and then also to a // TS content reduction. if (math.sum(Qsubstrates) == 0) { substrates_or_sludge = new biogas.substrates(); int ifermenter = 0; // for (int iflux = n_substrate; iflux < Q.Length; iflux++) { string digester_id = myPlant.getDigesterID(ifermenter + 1); double TS_digester, VS_digester = 0; try { mySensors.getCurrentMeasurementD("TS_" + digester_id + "_3", out TS_digester); mySensors.getCurrentMeasurementD("VS_" + digester_id + "_3", out VS_digester); // TODO - warum 4, wenn Fermenter abgestürzt ist, dann ist TS < 2 // was zu doofen Fehlermeldungen führt mit calcNfE und boundNDF, wenn man // VS unten in biogas.sludge setzt. deshalb hier abfrage if (TS_digester < 4 /*double.Epsilon*/) { TS_digester = 11; } // TODO - herausfinden warum 15 if (VS_digester < 20 /*double.Epsilon*/) { VS_digester = 85; // 85 % TS } } catch { TS_digester = 11; VS_digester = 85; } try { substrates_or_sludge.addSubstrate( new biogas.sludge(mySubstrates, math.ones(n_substrate), TS_digester, VS_digester)); } catch (exception e) { LogError.Log_Err(String.Format("OLR_sensor.doMeasurement, VS: {0}, TS: {1}", VS_digester, TS_digester), e); throw (e); } ifermenter = ifermenter + 1; } // for (int isubstrate = n_substrate; isubstrate < Q.Length; isubstrate++) { Q_s_or_s.Add(Q[isubstrate]); } } else { substrates_or_sludge = mySubstrates; for (int isubstrate = 0; isubstrate < Qsubstrates.Length; isubstrate++) { Q_s_or_s.Add(Qsubstrates[isubstrate]); } } // digester myDigester = myPlant.getDigesterByID(id_suffix); try { values[0] = myDigester.calcOLR(x, substrates_or_sludge, Q_s_or_s.ToArray(), Qsum); } catch (exception e) { LogError.Log_Err("OLR_sensor.doMeasurement2", e); throw (e); } // return(values); }