File contents
package ASM2d "Component models for the Activated Sludge Model No.2d"
extends Modelica.Icons.Library;
annotation (
Coordsys(
extent=[0, 0; 442, 386],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.44,
height=0.65,
library=1,
autolayout=1),
Documentation(info="This library contains components to build models of biological municipal
wastewater treatment plants based on the Activated Sludge Model No.2d (ASM2d) by the
International Association on Water Quality (IAWQ) [1].
The library currently is structured in following sub-libraries.
- Interfaces (partial ASM2d models and connectors)
- PreClar (primary clarifier models)
- SecClar (several secondary settling tank models)
- Examples (wastewater treatment plant models)
Main Author:
Gerald Reichl
Technische Universitaet Ilmenau
Faculty of Informatics and Automation
Department Dynamics and Simulation of ecological Systems
P.O. Box 10 05 65
98684 Ilmenau
Germany
email: gerald.reichl@tu-ilmenau.de
References:
[1] M. Henze and W.Gujer and T. Mino and. M.v. Loosdrecht: Activated Sludge
Models ASM1, ASM2, ASM2d, and ASM3. IWA Task Group on Mathematical Modelling
for Design and Operation of Biological Wastewater Treatment, 2000.
This package is free software; it can be redistributed and/or modified under the terms of the Modelica license, see the license conditions and the accompanying
disclaimer in the documentation of package Modelica in file \"Modelica/package.mo\".
Copyright (C) 2002 - 2003, Gerald Reichl
"));
model deni "ASM2d denitrification tank"
//denitrification tank based on the ASM2d model
extends WasteWater.Icons.deni;
extends Interfaces.ASM2dbase;
// tank specific parameters
parameter Modelica.SIunits.Volume V=1000 "Volume of denitrification tank";
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -10; -90, 10]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, -10; 110, 10]);
Interfaces.WWFlowAsm2dout MeasurePort annotation (extent=[50, 40; 60, 50]);
Modelica.Blocks.Interfaces.InPort T(final n=1) annotation (extent=[-110, 30;
-90, 50]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component models the ASM2d processes and reactions taking place in an unaerated (denitrification) tank
of a wastewater treatment plant.
The InPort signal gives the tank temperature to the model.
Parameters:
- all conversion factors and stoichiometric and kinetic parameters of the activated sludge model No.2d (ASM2d)
V - volume of the tank [m3]
"));
equation
aeration = 0;
// no aeration in this tank //
// volume dependent dilution term of each concentration
inputSo = (In.So - So)*In.Q/V;
inputSf = (In.Sf - Sf)*In.Q/V;
inputSa = (In.Sa - Sa)*In.Q/V;
inputSnh = (In.Snh - Snh)*In.Q/V;
inputSno = (In.Sno - Sno)*In.Q/V;
inputSpo = (In.Spo - Spo)*In.Q/V;
inputSi = (In.Si - Si)*In.Q/V;
inputSalk = (In.Salk - Salk)*In.Q/V;
inputSn2 = (In.Sn2 - Sn2)*In.Q/V;
inputXi = (In.Xi - Xi)*In.Q/V;
inputXs = (In.Xs - Xs)*In.Q/V;
inputXh = (In.Xh - Xh)*In.Q/V;
inputXpao = (In.Xpao - Xpao)*In.Q/V;
inputXpp = (In.Xpp - Xpp)*In.Q/V;
inputXpha = (In.Xpha - Xpha)*In.Q/V;
inputXa = (In.Xa - Xa)*In.Q/V;
inputXtss = (In.Xtss - Xtss)*In.Q/V;
inputXmeoh = (In.Xmeoh - Xmeoh)*In.Q/V;
inputXmep = (In.Xmep - Xmep)*In.Q/V;
end deni;
model nitri "ASM2d nitrification tank"
// nitrification (aerated) tank, based on the ASM2d model
extends WasteWater.Icons.nitri;
extends Interfaces.ASM2dbase;
// aeration system dependend parameters
parameter Modelica.SIunits.Volume V=1000 "Volume of nitrification tank";
parameter Real alpha=0.7 "Oxygen transfer factor";
parameter Modelica.SIunits.Length de=4.5 "depth of aeration";
parameter Real R_air=23.5 "specific oxygen feed factor [gO2/(m^3*m)]";
WWU.MassConcentration So_sat "Dissolved oxygen saturation";
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -10; -90, 10]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, -10; 110, 10]);
Interfaces.WWFlowAsm2dout MeasurePort annotation (extent=[50, 40; 60, 50]);
Modelica.Blocks.Interfaces.InPort T(final n=1) annotation (extent=[-110, 30;
-90, 50]);
Interfaces.AirFlow AirIn annotation (extent=[-5, -103; 5, -93]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component models the ASM2d processes and reactions taking place in an aerated (nitrification) tank
of a wastewater treatment plant.
The InPort signal gives the tank temperature to the model.
Parameters:
- all conversion factors and soichiometric and kinetic parameters of the activated sludge model No.2d (ASM2d)
V - volume of the tank [m3]
alpha - oxygen transfer factor
de - depth of the aeration system [m]
R_air - specific oxygen feed factor [g O2/(m3*m)]
"));
equation
// Temperature dependend oxygen saturation by Simba
So_sat = 13.89 + (-0.3825 + (0.007311 - 0.00006588*T.signal[1])*T.signal[1])*
T.signal[1];
// extends the Oxygen differential equation by an aeration term
// aeration [mgO2/l]; AirIn.Q_air needs to be in
// Simulationtimeunit [m3/day^-1]
aeration = (alpha*(So_sat - So)/So_sat*AirIn.Q_air*R_air*de)/V;
// aeration = Kla * (So_sat -So);
// volume dependent dilution term of each concentration
inputSo = (In.So - So)*In.Q/V;
inputSf = (In.Sf - Sf)*In.Q/V;
inputSa = (In.Sa - Sa)*In.Q/V;
inputSnh = (In.Snh - Snh)*In.Q/V;
inputSno = (In.Sno - Sno)*In.Q/V;
inputSpo = (In.Spo - Spo)*In.Q/V;
inputSi = (In.Si - Si)*In.Q/V;
inputSalk = (In.Salk - Salk)*In.Q/V;
inputSn2 = (In.Sn2 - Sn2)*In.Q/V;
inputXi = (In.Xi - Xi)*In.Q/V;
inputXs = (In.Xs - Xs)*In.Q/V;
inputXh = (In.Xh - Xh)*In.Q/V;
inputXpao = (In.Xpao - Xpao)*In.Q/V;
inputXpp = (In.Xpp - Xpp)*In.Q/V;
inputXpha = (In.Xpha - Xpha)*In.Q/V;
inputXa = (In.Xa - Xa)*In.Q/V;
inputXtss = (In.Xtss - Xtss)*In.Q/V;
inputXmeoh = (In.Xmeoh - Xmeoh)*In.Q/V;
inputXmep = (In.Xmep - Xmep)*In.Q/V;
end nitri;
model precipitation "Phosphorus precipitation tank"
extends WasteWater.Icons.precipitation;
extends Interfaces.ASM2dbase;
parameter Modelica.SIunits.Volume V=50 "Volume of precipitation tank";
parameter Real Wsub=2.5 "effective substance in precipitant [mol/kg]";
parameter Real D=1.4 "density of precipitant [kg/l]";
parameter Real beta=5.0 "precipitant surplus";
parameter Real Qmin=5.0 "Minimum flow of precipitant [l/h]";
parameter Real MP=30.97 "Molar Mass of Phosphorus [g/mol]";
parameter Real Mpre=55.85 "Molar Mass of precipitant [g/mol]";
WWU.VolumeFlowRate Qpreci "Dosage flow of precipitant used";
Real Preci "Concentration of effective substance in precipitant flow
[g/m³]";
Real H;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -10; -90, 10]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, -10; 110, 10]);
Interfaces.WWFlowAsm2dout MeasurePort annotation (extent=[50, 40; 60, 50]);
Modelica.Blocks.Interfaces.InPort T(final n=1) annotation (extent=[-110, 30;
-90, 50]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This ASM2d component is used to model the chemical Phosphorus precipitation
with FeCl3 as precipitant.
Parameters:
- all ASM2d conversion factors and stoichiometric and kinetic parameters
V - Volume of precipitation tank [m3]
Wsub - effective substance in precipitant [mol/kg]
D - density of precipitant [kg/l]
beta - precipitant surplus
Qmin - minimum flow rate of precipitant [l/h]
MP - molar mass of phosphorus [g/mol]
Mpre - molar mass of precipitant [g/mol]
"));
equation
aeration = 0;
//no aerotation in this tank
H = (1/Wsub)*(1/(D*1000))*(1/Mpre)*Mpre/MP*beta*In.Spo*In.Q;
Qpreci = if H < (Qmin*24/1000) then Qmin*24/1000 else H;
// the factor 1.91 is a stoichiometric coefficient between Fe3+ and Xmeoh
Preci = 1.91*Wsub*D*Mpre*1000;
// volume dependent dilution term of each concentration
inputSo = (In.So - So)*In.Q/V;
inputSf = (In.Sf - Sf)*In.Q/V;
inputSa = (In.Sa - Sa)*In.Q/V;
inputSnh = (In.Snh - Snh)*In.Q/V;
inputSno = (In.Sno - Sno)*In.Q/V;
inputSpo = (In.Spo - Spo)*In.Q/V;
inputSi = (In.Si - Si)*In.Q/V;
inputSalk = (In.Salk - Salk)*In.Q/V;
inputSn2 = (In.Sn2 - Sn2)*In.Q/V;
inputXi = (In.Xi - Xi)*In.Q/V;
inputXs = (In.Xs - Xs)*In.Q/V;
inputXh = (In.Xh - Xh)*In.Q/V;
inputXpao = (In.Xpao - Xpao)*In.Q/V;
inputXpp = (In.Xpp - Xpp)*In.Q/V;
inputXpha = (In.Xpha - Xpha)*In.Q/V;
inputXa = (In.Xa - Xa)*In.Q/V;
inputXtss = (In.Xtss - Xtss)*In.Q/V;
inputXmeoh = (In.Xmeoh - Xmeoh)*In.Q/V + Qpreci*Preci/V;
inputXmep = (In.Xmep - Xmep)*In.Q/V;
end precipitation;
model SecClarModTakacs "ASM2d Secondary Clarifier Model based on Takacs"
extends WasteWater.Icons.SecClar;
extends ASM2d.SecClar.Takacs.Interfaces.ratios;
package SCP = ASM2d.SecClar.Takacs;
package SI = Modelica.SIunits;
package WI = ASM2d.Interfaces;
package WWU = WasteWater.WasteWaterUnits;
parameter SI.Length hsc=4.0 "height of secondary clarifier";
parameter Integer n=10 "number of layers of SC model";
parameter SI.Length zm=hsc/(1.0*n) "height of m-th secondary clarifier layer"
;
parameter SI.Area Asc=1500.0 "area of secondary clarifier";
parameter WWU.MassConcentration Xt=3000.0 "threshold for X";
// total sludge concentration in clarifier feed
WWU.MassConcentration Xf;
WI.WWFlowAsm2din Feed annotation (extent=[-110, 4; -90, 24]);
WI.WWFlowAsm2dout Effluent annotation (extent=[92, 47; 112, 67]);
WI.WWFlowAsm2dout Return annotation (extent=[-40, -106; -20, -86]);
WI.WWFlowAsm2dout Waste annotation (extent=[20, -106; 40, -86]);
// layers 1 to 10
SCP.bottom_layer S1(
zm=zm,
Asc=Asc,
Xf=Xf,
rXi=rXi,
rXs=rXs,
rXh=rXh,
rXpao=rXpao,
rXpp=rXpp,
rXpha=rXpha,
rXa=rXa,
rXmeoh=rXmeoh,
rXmep=rXmep) annotation (extent=[-35, -93; 35, -78]);
SCP.lower_layer S2(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-35, -74; 35, -59]);
SCP.lower_layer S3(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-35, -55; 35, -40]);
SCP.lower_layer S4(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-35, -36; 35, -21]);
SCP.lower_layer S5(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-35, -17; 35, -2]);
SCP.feed_layer S6(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-35, 2; 35, 17]);
SCP.upper_layer S7(
zm=zm,
Asc=Asc,
Xf=Xf,
Xt=Xt) annotation (extent=[-35, 21; 35, 36]);
SCP.upper_layer S8(
zm=zm,
Asc=Asc,
Xt=Xt,
Xf=Xf) annotation (extent=[-35, 40; 35, 55]);
SCP.upper_layer S9(
zm=zm,
Asc=Asc,
Xf=Xf,
Xt=Xt) annotation (extent=[-35, 59; 35, 74]);
SCP.top_layer S10(
zm=zm,
Asc=Asc,
Xf=Xf,
Xt=Xt,
rXi=rXi,
rXs=rXs,
rXh=rXh,
rXpao=rXpao,
rXpp=rXpp,
rXpha=rXpha,
rXa=rXa,
rXmeoh=rXmeoh,
rXmep=rXmep) annotation (extent=[-35, 78; 35, 93]);
equation
connect(S1.Up, S2.Dn) annotation (points=[-2.22045e-15, -78; -2.22045e-15, -
74]);
connect(S2.Up, S3.Dn) annotation (points=[-2.22045e-15, -59; -2.22045e-15, -
55]);
connect(S3.Up, S4.Dn) annotation (points=[-2.22045e-15, -40; -2.22045e-15, -
36]);
connect(S5.Up, S6.Dn) annotation (points=[-2.22045e-15, -2; -2.22045e-15, 2])
;
connect(S6.Up, S7.Dn) annotation (points=[-2.22045e-15, 17; -2.22045e-15, 21]
);
connect(S7.Up, S8.Dn) annotation (points=[-2.22045e-15, 36; -2.22045e-15, 40]
);
connect(S9.Up, S10.Dn) annotation (points=[-2.22045e-15, 74; -2.22045e-15, 78
]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.12,
y=0.03,
width=0.35,
height=0.49),
Documentation(info="This component models an ASM2d 10 - layer secondary clarifier model with 4 layers above the feed_layer (including top_layer)
and 5 layers below the feed_layer (including bottom_layer) based on Takac`s theory.
Parameters:
hsc - height of clarifier [m]
n - number of layers
Asc - surface area of sec. clar. [m2]
Xt - threshold value for Xtss [mg/l]
"));
connect(S4.Up, S5.Dn) annotation (points=[-2.22045e-15, -21; -2.22045e-15, -
17]);
connect(S8.Up, S9.Dn) annotation (points=[-2.22045e-15, 55; -2.22045e-15, 59]
);
connect(Feed, S6.In) annotation (points=[-100, 10; -67.5, 10; -67.5, 9.8; -35
, 9.8]);
connect(S1.PQw, Waste) annotation (points=[17.5, -93; 17.5, -100; 30, -100]);
connect(S10.Out, Effluent) annotation (points=[35, 85.5; 67.5, 85.5; 67.5, 57
; 100, 57]);
connect(S1.PQr, Return) annotation (points=[-21, -93; -21, -100; -30, -100]);
// total sludge concentration in clarifier feed
Xf = Feed.Xtss;
// ratios of solid components
rXi = Feed.Xi/Xf;
rXs = Feed.Xs/Xf;
rXh = Feed.Xh/Xf;
rXpao = Feed.Xpao/Xf;
rXpp = Feed.Xpp/Xf;
rXpha = Feed.Xpha/Xf;
rXa = Feed.Xa/Xf;
rXmeoh = Feed.Xmeoh/Xf;
rXmep = Feed.Xmep/Xf;
end SecClarModTakacs;
model blower "Blower for the aeration of the nitrification tanks"
package WWU = WasteWater.WasteWaterUnits;
extends WasteWater.Icons.blower;
parameter WWU.VolumeFlowRate Q_max=20000 "maximum blower capacity";
parameter WWU.VolumeFlowRate Q_min=0.0 "minimum blower capacity";
Real H;
// this is just a help variable to reduce expressions
Interfaces.AirFlow AirOut annotation (extent=[-20, 90; 0, 110]);
Modelica.Blocks.Interfaces.InPort u(final n=1) annotation (extent=[88, -40;
108, -20], rotation=180);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component models a blower of a wastewater treatment plant which generates an airflow that is needed
for the nitrification.
The blower is connected to the nitrification tank.
The airflow is controlled by a signal u (-1 <= u <= 1).
Parameter:
Qmax - maximum blower capacity [m3 Air/d], this is produced when the control signal u is 1 or greater.
Qmin - minimum blower capacity [m3 Air/d], this is produced when the control signal u is -1 or below.
"));
equation
H = 0.5*(-Q_min + Q_max) + u.signal[1]*0.5*(-Q_min + Q_max) + Q_min;
AirOut.Q_air = -(if H > Q_max then Q_max else if H < Q_min then Q_min else H)
;
end blower;
model pump "ASM2d wastewater pump"
package WWU = WasteWater.WasteWaterUnits;
extends WasteWater.Icons.pump;
parameter WWU.VolumeFlowRate Q_min=0.0 "minimum pump capacity";
parameter WWU.VolumeFlowRate Q_max=20000 "maximum pump capacity";
Real H;
//help variable to reduce expressions
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -43; -90, -23]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, 18; 110, 38]);
Modelica.Blocks.Interfaces.InPort u(final n=1) annotation (extent=[-98, 15; -
78, 35]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component models an ASM2d wastewater pump. It generates a wastewater flow
that is controlled by the signal u (-1 <= u <=1).
Parameter:
Qmax - maximum pump capacity [m3/d], this is produced when the control signal u is 1 or greater.
Qmin - minimum pump capacity [m3/d], this is produced when the control signal u is -1 or below.
"));
equation
H = 0.5*(-Q_min + Q_max) + u.signal[1]*0.5*(-Q_min + Q_max) + Q_min;
Out.Q = -(if H > Q_max then Q_max else if H < Q_min then Q_min else H);
Out.Q + In.Q = 0;
Out.So = In.So;
Out.Sf = In.Sf;
Out.Sa = In.Sa;
Out.Snh = In.Snh;
Out.Sno = In.Sno;
Out.Spo = In.Spo;
Out.Si = In.Si;
Out.Salk = In.Salk;
Out.Sn2 = In.Sn2;
Out.Xi = In.Xi;
Out.Xs = In.Xs;
Out.Xh = In.Xh;
Out.Xpao = In.Xpao;
Out.Xpp = In.Xpp;
Out.Xpha = In.Xpha;
Out.Xa = In.Xa;
Out.Xtss = In.Xtss;
Out.Xmeoh = In.Xmeoh;
Out.Xmep = In.Xmep;
end pump;
model FlowSource "Flowsource"
extends WasteWater.Icons.FlowSource;
Interfaces.WWFlowAsm2dout Out annotation (extent=[88, -80; 108, -60]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component is used to feed an ASM2d wwtp model with flow data from measurement
when e.g. concentration is measured after the primary clarifier.
The dimension of InPort is 1.
1 volumeflowrate Q of incoming wastewater [m3/d]"),
Diagram(
Ellipse(extent=[-54, 54; 56, -54], style(color=8, fillColor=8)),
Polygon(points=[-8, -54; -14, -52; -24, -48; -32, -44; -36, -40; -42, -34
; -48, -26; -50, -20; 52, -20; 50, -26; 46, -32; 42, -36; 38, -40;
34, -44; 30, -46; 26, -48; 22, -50; 16, -52; 10, -54; 4, -54; 0, -
54; -8, -54], style(pattern=0, fillColor=70)),
Ellipse(extent=[-54, 54; 56, -54], style(color=0, thickness=2)),
Rectangle(extent=[-4, -52; 4, -74], style(pattern=0, fillColor=70)),
Rectangle(extent=[4, -74; 88, -68], style(pattern=0, fillColor=70)),
Line(points=[-4, -54; -4, -74; 88, -74], style(color=0, thickness=2)),
Line(points=[4, -54; 4, -68; 88, -68], style(color=0, thickness=2))));
Modelica.Blocks.Interfaces.InPort data(final n=1) annotation (extent=[-100, -
10; -80, 10]);
equation
Out.Q = -data.signal[1];
end FlowSource;
model WWSource "Wastewater source"
extends WasteWater.Icons.WWSource;
extends ASM2d.Interfaces.conversion_factors;
Interfaces.WWFlowAsm2dout Out annotation (extent=[88, -80; 108, -60]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component provides all ASM2d data at the influent of a wastewater treatment plant.
The dimension of InPort is 19.
1 volumeflowrate Q of incoming wastewater [m3/d]
2 So [g O2/m3]
3 Sf [g COD/m3]
4 Sa [g COD/m3]
5 Snh [g N/m3]
6 Sno [g N/m3]
7 Spo [g P/m3]
8 Si [g COD/m3]
9 Salk [mmol/l]
10 Sn2 [g N/m3]
11 Xi [g COD/m3]
12 Xs [g COD/m3]
13 Xh [g COD/m3]
14 Xpao [g COD/m3]
15 Xpp [g P/m3]
16 Xpha [g COD/m3]
17 Xa [g COD/m3]
18 Xmeoh [g TSS/m3]
19 Xmep [g TSS/m3]
Parameters:
- all ASM2d conversion factors for the calculation of Xtss."));
Modelica.Blocks.Interfaces.InPort data(final n=19) annotation (extent=[-100,
-10; -80, 10]);
equation
Out.Q = -data.signal[1];
Out.So = data.signal[2];
Out.Sf = data.signal[3];
Out.Sa = data.signal[4];
Out.Snh = data.signal[5];
Out.Sno = data.signal[6];
Out.Spo = data.signal[7];
Out.Si = data.signal[8];
Out.Salk = data.signal[9];
Out.Sn2 = data.signal[10];
Out.Xi = data.signal[11];
Out.Xs = data.signal[12];
Out.Xh = data.signal[13];
Out.Xpao = data.signal[14];
Out.Xpp = data.signal[15];
Out.Xpha = data.signal[16];
Out.Xa = data.signal[17];
Out.Xtss = i_TSS_Xi*Out.Xi + i_TSS_Xs*Out.Xs + i_TSS_BM*(Out.Xh + Out.Xpao +
Out.Xa) + 3.23*Out.Xpp + 0.6*Out.Xpha + Out.Xmeoh + Out.Xmep;
Out.Xmeoh = data.signal[18];
Out.Xmep = data.signal[19];
end WWSource;
model EffluentSink "Receiving water (river)"
// only for graphical termination in diagram layer, no equations needed
extends WasteWater.Icons.EffluentSink;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, 10; -90, 30]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component terminates an ASM2d wastewater treatment plant model e.g. the wastewater flow to the receiving water.
"));
equation
end EffluentSink;
model SludgeSink "Wastesludge sink"
// only for graphical termination in diagram layer, no equations needed
extends WasteWater.Icons.SludgeSink;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -22; -90, -2]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component terminates the waste sludge stream of an ASM2d wastewater treatment plant model.
Storage or further sludge treatment is not jet considered."));
equation
end SludgeSink;
model ControlledDivider2 "Controlled flow divider"
// divides one flow of wastewater into 2 Flows controlled by the
// input signal u; u=1 means Out1.Q=In.Q and u=0 means Out2.Q=In.Q
extends WasteWater.Icons.ControlledDivider2;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -7; -90, 13]);
Interfaces.WWFlowAsm2dout Out1 annotation (extent=[90, 15; 110, 35]);
Interfaces.WWFlowAsm2dout Out2 annotation (extent=[90, -25; 110, -5]);
Modelica.Blocks.Interfaces.InPort u(final n=1) annotation (
extent=[-10, -70; 10, -50],
rotation=90,
layer="icon");
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.29,
y=0.31,
width=0.68,
height=0.6),
Documentation(info="This component divides one wastewater flow (ASM2d) into two flows which are controlled by the signal u (0...1).
Is u.signal=1, the flow goes to output 1 (Out1) and is u.signal=0, the flow goes to output 2 (Out2).
The concentrations of the outport-flows are equal to the concentration at inport."
));
equation
Out1.Q = -In.Q*u.signal[1];
Out2.Q = -In.Q*(1 - u.signal[1]);
Out1.So = In.So;
Out1.Sf = In.Sf;
Out1.Sa = In.Sa;
Out1.Snh = In.Snh;
Out1.Sno = In.Sno;
Out1.Spo = In.Spo;
Out1.Si = In.Si;
Out1.Salk = In.Salk;
Out1.Sn2 = In.Sn2;
Out1.Xi = In.Xi;
Out1.Xs = In.Xs;
Out1.Xh = In.Xh;
Out1.Xpao = In.Xpao;
Out1.Xpp = In.Xpp;
Out1.Xpha = In.Xpha;
Out1.Xa = In.Xa;
Out1.Xtss = In.Xtss;
Out1.Xmeoh = In.Xmeoh;
Out1.Xmep = In.Xmep;
Out2.So = In.So;
Out2.Sf = In.Sf;
Out2.Sa = In.Sa;
Out2.Snh = In.Snh;
Out2.Sno = In.Sno;
Out2.Spo = In.Spo;
Out2.Si = In.Si;
Out2.Salk = In.Salk;
Out2.Sn2 = In.Sn2;
Out2.Xi = In.Xi;
Out2.Xs = In.Xs;
Out2.Xh = In.Xh;
Out2.Xpao = In.Xpao;
Out2.Xpp = In.Xpp;
Out2.Xpha = In.Xpha;
Out2.Xa = In.Xa;
Out2.Xtss = In.Xtss;
Out2.Xmeoh = In.Xmeoh;
Out2.Xmep = In.Xmep;
end ControlledDivider2;
model divider2 "Flowdivider"
// divides one flow of wastewater into 2 Flows; one amout needs to be specified
extends WasteWater.Icons.divider2;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -7; -90, 13]);
Interfaces.WWFlowAsm2dout Out1 annotation (extent=[90, 15; 110, 35]);
Interfaces.WWFlowAsm2dout Out2 annotation (extent=[90, -25; 110, -5]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info=
"This component divides one ASM2d wastewater flow into two ASM2d wastewater flows."
));
equation
In.Q + Out1.Q + Out2.Q = 0;
Out1.So = In.So;
Out1.Sf = In.Sf;
Out1.Sa = In.Sa;
Out1.Snh = In.Snh;
Out1.Sno = In.Sno;
Out1.Spo = In.Spo;
Out1.Si = In.Si;
Out1.Salk = In.Salk;
Out1.Sn2 = In.Sn2;
Out1.Xi = In.Xi;
Out1.Xs = In.Xs;
Out1.Xh = In.Xh;
Out1.Xpao = In.Xpao;
Out1.Xpp = In.Xpp;
Out1.Xpha = In.Xpha;
Out1.Xa = In.Xa;
Out1.Xtss = In.Xtss;
Out1.Xmeoh = In.Xmeoh;
Out1.Xmep = In.Xmep;
Out2.So = In.So;
Out2.Sf = In.Sf;
Out2.Sa = In.Sa;
Out2.Snh = In.Snh;
Out2.Sno = In.Sno;
Out2.Spo = In.Spo;
Out2.Si = In.Si;
Out2.Salk = In.Salk;
Out2.Sn2 = In.Sn2;
Out2.Xi = In.Xi;
Out2.Xs = In.Xs;
Out2.Xh = In.Xh;
Out2.Xpao = In.Xpao;
Out2.Xpp = In.Xpp;
Out2.Xpha = In.Xpha;
Out2.Xa = In.Xa;
Out2.Xtss = In.Xtss;
Out2.Xmeoh = In.Xmeoh;
Out2.Xmep = In.Xmep;
end divider2;
model mixer2 "Mixer of two ASM2d characterised flows"
extends WasteWater.Icons.mixer2;
Interfaces.WWFlowAsm2din In1 annotation (extent=[-110, 15; -90, 35]);
Interfaces.WWFlowAsm2din In2 annotation (extent=[-110, -25; -90, -5]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, -6; 110, 14]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info=
"This component mixes two flows of wastewater (ASM2d) of different concentration and different amount."
));
equation
In1.Q + In2.Q + Out.Q = 0;
Out.So = (In1.So*In1.Q + In2.So*In2.Q)/(In1.Q + In2.Q);
Out.Sf = (In1.Sf*In1.Q + In2.Sf*In2.Q)/(In1.Q + In2.Q);
Out.Sa = (In1.Sa*In1.Q + In2.Sa*In2.Q)/(In1.Q + In2.Q);
Out.Snh = (In1.Snh*In1.Q + In2.Snh*In2.Q)/(In1.Q + In2.Q);
Out.Sno = (In1.Sno*In1.Q + In2.Sno*In2.Q)/(In1.Q + In2.Q);
Out.Spo = (In1.Spo*In1.Q + In2.Spo*In2.Q)/(In1.Q + In2.Q);
Out.Si = (In1.Si*In1.Q + In2.Si*In2.Q)/(In1.Q + In2.Q);
Out.Salk = (In1.Salk*In1.Q + In2.Salk*In2.Q)/(In1.Q + In2.Q);
Out.Sn2 = (In1.Sn2*In1.Q + In2.Sn2*In2.Q)/(In1.Q + In2.Q);
Out.Xi = (In1.Xi*In1.Q + In2.Xi*In2.Q)/(In1.Q + In2.Q);
Out.Xs = (In1.Xs*In1.Q + In2.Xs*In2.Q)/(In1.Q + In2.Q);
Out.Xh = (In1.Xh*In1.Q + In2.Xh*In2.Q)/(In1.Q + In2.Q);
Out.Xpao = (In1.Xpao*In1.Q + In2.Xpao*In2.Q)/(In1.Q + In2.Q);
Out.Xpp = (In1.Xpp*In1.Q + In2.Xpp*In2.Q)/(In1.Q + In2.Q);
Out.Xpha = (In1.Xpha*In1.Q + In2.Xpha*In2.Q)/(In1.Q + In2.Q);
Out.Xa = (In1.Xa*In1.Q + In2.Xa*In2.Q)/(In1.Q + In2.Q);
Out.Xtss = (In1.Xtss*In1.Q + In2.Xtss*In2.Q)/(In1.Q + In2.Q);
Out.Xmeoh = (In1.Xmeoh*In1.Q + In2.Xmeoh*In2.Q)/(In1.Q + In2.Q);
Out.Xmep = (In1.Xmep*In1.Q + In2.Xmep*In2.Q)/(In1.Q + In2.Q);
end mixer2;
model mixer3 "Mixer of 3 ASM2d characterised flows"
extends WasteWater.Icons.mixer3;
Interfaces.WWFlowAsm2din In1 annotation (extent=[-110, 25; -90, 45]);
Interfaces.WWFlowAsm2din In2 annotation (extent=[-110, -15; -90, 5]);
Interfaces.WWFlowAsm2din In3 annotation (extent=[-110, -55; -90, -35]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[90, -14; 110, 6]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info=
"This component mixes 3 flows of wastewater (ASM2d) of different concentration and different amount."
));
equation
In1.Q + In2.Q + In3.Q + Out.Q = 0;
Out.So = (In1.So*In1.Q + In2.So*In2.Q + In3.So*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Sf = (In1.Sf*In1.Q + In2.Sf*In2.Q + In3.Sf*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Sa = (In1.Sa*In1.Q + In2.Sa*In2.Q + In3.Sa*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Snh = (In1.Snh*In1.Q + In2.Snh*In2.Q + In3.Snh*In3.Q)/(In1.Q + In2.Q +
In3.Q);
Out.Sno = (In1.Sno*In1.Q + In2.Sno*In2.Q + In3.Sno*In3.Q)/(In1.Q + In2.Q +
In3.Q);
Out.Spo = (In1.Spo*In1.Q + In2.Spo*In2.Q + In3.Spo*In3.Q)/(In1.Q + In2.Q +
In3.Q);
Out.Si = (In1.Si*In1.Q + In2.Si*In2.Q + In3.Si*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Salk = (In1.Salk*In1.Q + In2.Salk*In2.Q + In3.Salk*In3.Q)/(In1.Q + In2.Q
+ In3.Q);
Out.Sn2 = (In1.Sn2*In1.Q + In2.Sn2*In2.Q + In3.Sn2*In3.Q)/(In1.Q + In2.Q +
In3.Q);
Out.Xi = (In1.Xi*In1.Q + In2.Xi*In2.Q + In3.Xi*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Xs = (In1.Xs*In1.Q + In2.Xs*In2.Q + In3.Xs*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Xh = (In1.Xh*In1.Q + In2.Xh*In2.Q + In3.Xh*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Xpao = (In1.Xpao*In1.Q + In2.Xpao*In2.Q + In3.Xpao*In3.Q)/(In1.Q + In2.Q
+ In3.Q);
Out.Xpp = (In1.Xpp*In1.Q + In2.Xpp*In2.Q + In3.Xpp*In3.Q)/(In1.Q + In2.Q +
In3.Q);
Out.Xpha = (In1.Xpha*In1.Q + In2.Xpha*In2.Q + In3.Xpha*In3.Q)/(In1.Q + In2.Q
+ In3.Q);
Out.Xa = (In1.Xa*In1.Q + In2.Xa*In2.Q + In3.Xa*In3.Q)/(In1.Q + In2.Q + In3.Q)
;
Out.Xtss = (In1.Xtss*In1.Q + In2.Xtss*In2.Q + In3.Xtss*In3.Q)/(In1.Q + In2.Q
+ In3.Q);
Out.Xmeoh = (In1.Xmeoh*In1.Q + In2.Xmeoh*In2.Q + In3.Xmeoh*In3.Q)/(In1.Q +
In2.Q + In3.Q);
Out.Xmep = (In1.Xmep*In1.Q + In2.Xmep*In2.Q + In3.Xmep*In3.Q)/(In1.Q + In2.Q
+ In3.Q);
end mixer3;
model sensor_COD "Ideal sensor to measure chemical oxygen demand (COD)"
extends WasteWater.Icons.sensor_COD;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort COD(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the chemical oxygen demand (COD) concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"));
equation
In.Q = 0.0;
COD.signal[1] = In.Sf + In.Sa + In.Si + In.Xi + In.Xs + In.Xh + In.Xpao + In.
Xpha + In.Xa;
end sensor_COD;
model sensor_NH "Ideal sensor to measure ammonium nitrogen"
extends WasteWater.Icons.sensor_NH;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort Snh(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.63,
y=0.04,
width=0.35,
height=0.49),
Documentation(info="This component measures the ammonium nitrogen concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"));
equation
In.Q = 0;
Snh.signal[1] = In.Snh;
end sensor_NH;
model sensor_NO "Ideal sensor to measure nitrate nitrogen"
extends WasteWater.Icons.sensor_NO;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort Sno(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the nitrate nitrogen concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"));
equation
In.Q = 0;
Sno.signal[1] = In.Sno;
end sensor_NO;
model sensor_O2 "Ideal sensor to measure dissolved oxygen concentration"
extends WasteWater.Icons.sensor_O2;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort So(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the dissolved oxygen concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"),
Diagram(
Ellipse(extent=[-50, 50; 50, -50], style(
color=0,
thickness=2,
fillColor=52)),
Text(extent=[-80, 100; 80, 60], string="%name"),
Line(points=[0, 50; 0, 38], style(color=0, thickness=2)),
Line(points=[-50, 0; 38, 0], style(color=0, thickness=2)),
Line(points=[50, 0; 38, 0], style(color=0, thickness=2)),
Line(points=[-36, 34; -28, 26], style(color=0, thickness=2)),
Line(points=[34, 36; 26, 28], style(color=0, thickness=2)),
Line(points=[0, 0; 26, 28], style(color=0, thickness=2)),
Polygon(points=[30, 32; 10, 24; 24, 12; 30, 32], style(color=0, fillColor
=0)),
Text(extent=[-36, -10; 36, -32], string="O2"),
Line(points=[0, -50; 0, -90], style(color=0, thickness=2)),
Line(points=[50, 0; 88, 0], style(color=0))));
equation
In.Q = 0;
So.signal[1] = In.So;
end sensor_O2;
model sensor_PO "Ideal sensor to measure dissolved phosphorus"
extends WasteWater.Icons.sensor_PO;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort Spo(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.57,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the dissolves phosphorus concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"));
equation
In.Q = 0;
Spo.signal[1] = In.Spo;
end sensor_PO;
model sensor_Q
"Ideal sensor to measure the flow rate of an ASM2d wastewater stream"
extends WasteWater.Icons.sensor_Q;
Interfaces.WWFlowAsm2din In annotation (extent=[-110, -10; -90, 10]);
Interfaces.WWFlowAsm2dout Out annotation (extent=[92, -10; 112, 10]);
Modelica.Blocks.Interfaces.OutPort Q(final n=1) annotation (extent=[-10, -110
; 10, -90], rotation=-90);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the flow of an ASM2d wastewater stream and provides
the result as output signal (to be further processed with blocks of
the Modelica.Blocks library).
"));
equation
In.Q + Out.Q = 0;
Q.signal[1] = In.Q;
// eventually abs(In.Q) to be shure to have pos. signal
In.So = Out.So;
In.Sf = Out.Sf;
In.Sa = Out.Sa;
In.Snh = Out.Snh;
In.Sno = Out.Sno;
In.Spo = Out.Spo;
In.Si = Out.Si;
In.Salk = Out.Salk;
In.Sn2 = Out.Sn2;
In.Xi = Out.Xi;
In.Xs = Out.Xs;
In.Xh = Out.Xh;
In.Xpao = Out.Xpao;
In.Xpp = Out.Xpp;
In.Xpha = Out.Xpha;
In.Xa = Out.Xa;
In.Xtss = Out.Xtss;
In.Xmeoh = Out.Xmeoh;
In.Xmep = Out.Xmep;
end sensor_Q;
model sensor_TKN "Ideal TKN and total nitrogen sensor"
extends WasteWater.Icons.sensor_TKN;
extends Interfaces.conversion_factors;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort TKN(final n=2) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.57,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the Total Kjeldal Nitrogen (TKN) and the
total nitrogen (N_total) concentration [g/m3] of ASM2d wastewater
and provides the result as output signal (to be further processed
with blocks of the Modelica.Blocks library).
signal[1] - TKN
signal[2] - N_total
"));
equation
In.Q = 0.0;
TKN.signal[1] = In.Snh + i_N_Si*In.Si + i_N_Sf*In.Sf + i_N_Xi*In.Xi + i_N_Xs*
In.Xs + i_N_BM*(In.Xh + In.Xpao + In.Xa);
TKN.signal[2] = TKN.signal[1] + In.Sno;
end sensor_TKN;
model sensor_TP
"Ideal sensor to measure the total phosphorus concentration in ASM2d wastewater"
extends WasteWater.Icons.sensor_TP;
extends Interfaces.conversion_factors;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
annotation (Documentation(info="This component measures the total phosphorus concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"), Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]));
Modelica.Blocks.Interfaces.OutPort TP(final n=1) annotation (extent=[88, -10
; 108, 10]);
equation
In.Q = 0.0;
TP.signal[1] = In.Spo + In.Xpp + i_P_Sf*In.Sf + i_P_Si*In.Si + i_P_Xi*In.Xi
+ i_P_Xs*In.Xs + i_P_BM*(In.Xh + In.Xa + In.Xpao) + In.Xmep/4.87;
end sensor_TP;
model sensor_TSS
"Ideal sensor to measure total suspended solids concentration (ASM2d)"
extends WasteWater.Icons.sensor_TSS;
Interfaces.WWFlowAsm2din In annotation (extent=[-10, -110; 10, -90]);
Modelica.Blocks.Interfaces.OutPort TSS(final n=1) annotation (extent=[88, -10
; 108, 10]);
annotation (
Coordsys(
extent=[-100, -100; 100, 100],
grid=[2, 2],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component measures the total suspended solids concentration [g/m3]
of ASM2d wastewater and provides the result as output signal (to be
further processed with blocks of the Modelica.Blocks library).
"));
equation
In.Q = 0;
TSS.signal[1] = In.Xtss;
end sensor_TSS;
end ASM2d;
Click here to get the file