File contents
package Examples "Demonstration examples of the components of the ASM1 library"
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 package contains example ASM1 wastewater treatment plant models to demonstrate the usage of
the WasteWater.ASM1 library.
Open the models and simulate them according to the description provided in the models.
The following demo models are present:
- SmallPlant
- BenchPlant
- ComplexPlant
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
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) 2003, Gerald Reichl
"));
class SmallPlant "Small WWTP Configuration"
import WasteWater;
extends Modelica.Icons.Example;
//Q_air=12100.99290780142 is equal to a Kla of 3.5 h^-1 from COST benchmark
//Q_air=34574.2654508612 is equal to a Kla of 10 h^-1 from COST benchmark
annotation (
Coordsys(
extent=[-120, -100; 120, 105],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.33,
y=0.04,
width=0.85,
height=0.87),
Diagram(Line(points=[-22, 58; -22, 58], style(color=3))),
Documentation(info="This fictitious plant provides an ASM1 example model with a small number of equations.
It consists of one denitrification and 2 nitrification tanks and a settler.
Change into the directory ../ASM1 and translate the model.
Before simulating the model load initial values from the script file small_asm1.mos
that is provided besides the model.
A 14 days dynamic influent data file is provided. So you may simulate up to 14 days.
But start with 1 day as it may take some time for simulation.
After simulation you may have a look at internal concentrations but most interesting
are the relevant concentrations at the effluent of a plant which can be viewed via the
sensors at the effluent of the secondary clarifier.
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
"));
ASM1.EffluentSink Effluent annotation (extent=[88, -28; 108, -8]);
ASM1.SludgeSink WasteSludge annotation (extent=[87, -51; 107, -31]);
ASM1.divider2 divider annotation (extent=[20, -6; 40, 14]);
ASM1.nitri tank3(V=1333) annotation (extent=[-6, -6; 14, 14]);
ASM1.nitri tank2(V=1333) annotation (extent=[-34, -6; -14, 14]);
ASM1.deni tank1(V=3000) annotation (extent=[-65, -6; -45, 14]);
ASM1.mixer3 mixer annotation (extent=[-104, 22; -84, 42]);
ModelicaAdditions.Tables.CombiTableTime CombiTableTime(
tableName="Inf_dry",
fileName="Inf_dry.txt",
icol={16,3,4,5,6,7,8,9,10,11,12,13,14,15}) annotation (extent=[-114, 78;
-94, 98]);
ASM1.WWSource WWSource annotation (extent=[-88, 78; -68, 98]);
ASM1.blower blower1(Q_max=34574.2654508612) annotation (extent=[-33, -62; -
13, -42]);
ASM1.blower blower2(Q_max=34574.2654508612) annotation (extent=[-6, -62; 14
, -42]);
ASM1.sensor_O2 sensor_O2 annotation (extent=[0, 24; 20, 44]);
Modelica.Blocks.Math.Feedback Feedback annotation (extent=[62, 40; 82, 60])
;
Modelica.Blocks.Continuous.PI PI1(k={500}, T={0.001}) annotation (extent=[
88, 40; 108, 60]);
Modelica.Blocks.Sources.Constant Constant1 annotation (extent=[-67, -87; -
47, -67]);
ASM1.pump RecyclePump(Q_max=46115) annotation (extent=[-94, -22; -74, -2],
rotation=180);
ASM1.pump ReturnPump(Q_max=9223) annotation (extent=[16, -36; 36, -16],
rotation=180);
ASM1.pump WastePump(Q_max=193) annotation (extent=[59, -55; 79, -35]);
Modelica.Blocks.Sources.Constant Constant2 annotation (extent=[22, -68; 42
, -48]);
Modelica.Blocks.Sources.Constant Temperature(k={15}) annotation (extent=[-
94, 50; -82, 62]);
sensor_NH sensor_NH1 annotation (extent=[64, 15; 80, 31]);
WasteWater.ASM1.sensor_NO sensor_NO1 annotation (extent=[81, 15; 97, 31]);
WasteWater.ASM1.sensor_TKN sensor_TKN1 annotation (extent=[97, 14; 113, 30]
);
WasteWater.ASM1.sensor_COD sensor_COD1 annotation (extent=[97, -5; 113, 11]
);
Modelica.Blocks.Sources.Step OxygenSetpoint(height={1.5}) annotation (
extent=[37, 40; 57, 60]);
WasteWater.ASM1.SecClar.Krebs.SecClarModKrebs Settler annotation (extent=[
48, -4; 68, 16]);
WasteWater.ASM1.sensor_TSS sensor_TSS1 annotation (extent=[32, 14; 49, 30])
;
equation
connect(tank3.Out, divider.In) annotation (points=[14, 4; 20, 4]
);
connect(mixer.Out, tank1.In) annotation (points=[-84, 32; -77, 32; -77, 4;
-65, 4]);
connect(mixer.In1, WWSource.Out) annotation (points=[-104, 36; -104, 74; -
68, 74; -68, 81; -68.2, 81]);
connect(CombiTableTime.outPort, WWSource.data) annotation (points=[-93, 88
; -88, 88]);
connect(blower2.AirOut, tank3.AirIn) annotation (points=[4, -42; 4, -5.8]);
connect(Feedback.outPort, PI1.inPort) annotation (points=[81, 50; 86, 50]);
connect(PI1.outPort, blower2.u) annotation (points=[109, 50; 114, 50; 114,
-84; 18, -84; 18, -55; 14, -55; 14, -56]);
connect(divider.Out2, RecyclePump.In) annotation (points=[40, 2; 44, 2; 44
, -8.7; -74, -8.7]);
connect(RecyclePump.Out, mixer.In3) annotation (points=[-94.2, -14.8; -104
, -14.8; -104, 28]);
connect(ReturnPump.Out, mixer.In2) annotation (points=[16, -28; 15.5, -28;
15.5, -30; -112, -30; -112, 32; -104, 32]);
connect(sensor_O2.So, Feedback.inPort2) annotation (points=[18.8, 34; 72,
34; 72, 42]);
connect(Temperature.outPort, tank1.T) annotation (points=[-81.5, 56; -71,
56; -71, 8; -65, 8], style(color=3));
connect(Temperature.outPort, tank2.T) annotation (points=[-81.5, 56; -39,
56; -39, 8; -34, 8], style(color=3));
connect(Temperature.outPort, tank3.T) annotation (points=[-81.5, 56; -39,
56; -39, 14; -5.9, 14; -5.9, 8; -5.8, 8], style(color=3));
connect(OxygenSetpoint.outPort, Feedback.inPort1) annotation (points=[58,
50; 64, 50], style(color=3));
connect(Constant1.outPort, blower1.u) annotation (points=[-46, -77; -7.2, -
77; -7.2, -55; -14.2, -55], style(color=3));
connect(WastePump.Out, WasteSludge.In) annotation (points=[79, -42; 87, -42
]);
connect(WastePump.u, Constant2.outPort) annotation (points=[60, -43; 46, -
43; 46, -58; 44, -58], style(color=3));
connect(tank2.Out, tank3.In) annotation (points=[-13, 4; -6, 4]);
connect(tank1.Out, tank2.In) annotation (points=[-44, 4; -34, 4]);
connect(blower1.AirOut, tank2.AirIn) annotation (points=[-23, -42; -23, -24
; -24, -24; -24, -6]);
connect(Constant1.outPort, RecyclePump.u) annotation (points=[-46, -77; -39
, -77; -39, -15; -75.1, -15], style(color=3));
connect(Settler.Effluent, Effluent.In) annotation (points=[69.2, 11.7; 78,
11.7; 78, -15; 88, -15]);
connect(Settler.Return, ReturnPump.In) annotation (points=[56, -4.6; 56, -
22.7; 36, -22.7]);
connect(WastePump.In, Settler.Waste) annotation (points=[59, -48; 52, -48;
52, -31; 62, -31; 62, -4.6]);
connect(sensor_NH1.In, Settler.Effluent) annotation (points=[72, 15; 72,
11.7; 69.2, 11.7]);
connect(sensor_NO1.In, Settler.Effluent) annotation (points=[89, 15; 89,
11.7; 69.2, 11.7]);
connect(sensor_TKN1.In, Settler.Effluent) annotation (points=[105, 14; 105
, 11.7; 69.2, 11.7]);
connect(sensor_COD1.In, Settler.Effluent) annotation (points=[105, -5; 105
, 11.7; 69.2, 11.7]);
connect(Constant2.outPort, ReturnPump.u) annotation (points=[43, -58; 48, -
58; 48, -28.5; 34.9, -28.5], style(color=3));
connect(divider.Out1, Settler.Feed) annotation (points=[40, 7; 44.5, 7;
44.5, 7.4; 49, 7.4]);
connect(tank3.MeasurePort, sensor_O2.In) annotation (points=[10, 9; 10, 25]
);
connect(sensor_TSS1.In, divider.Out1) annotation (points=[40, 14; 40, 6]);
end SmallPlant;
class BenchPlant "COST Benchmark WWTP Configuration"
import WasteWater;
//Q_air=34574.2654508612 is equal to a Kla of 10 h^-1 from COST benchmark
//Q_air=12100.99290780142 is equal to a Kla of 3.5 h^-1 from COST benchmark
extends Modelica.Icons.Example;
annotation (
Coordsys(
extent=[-120, -100; 120, 105],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.08,
y=0.03,
width=0.35,
height=0.49),
Diagram,
Documentation(info="This ASM1 plant consists of 2 denitrification tanks (tank1 and tank2),
3 nitrification tanks (tank3 - tank5) and a secondary clarifier by Takacs.
Furthermore there are 2 control loops modelled.
This configuration corresponds to the COST simulation benchmark [1].
Change into the directory ../ASM1 and translate the model.
Before simulating the model load initial values from the script file bench_asm1.mos
that is provided besides the model.
A 14 days dynamic influent data file is provided. So you may simulate up to 14 days.
But start with 1 day as it may take some time for simulation.
After simulation you may have a look at internal concentrations but most interesting
are the relevant concentrations at the effluent of a plant which can be viewed via the
sensors at the effluent of the secondary clarifier.
References:
[1] J.B. Copp: The COST Simulation Benchmark. 2000. http://www.ensic.u-nancy.fr/COSTWWTP/
PS: For those who want to reproduce the exact figures from the COST simulation benchmark some remarks:
The aeration system in this library is different from that in COST, so be sure to produce an airflow
corresponding to the desired Kla in COST. Furthermore in this library biological parameters are standard
parameters from the ASM1 distribution and implemented with temperature dependency which may vary a bit from
the parameter set used in COST.
But it is possible. During the validation phase of this library the steady state and dynamic results
from the COST simulation benchmark could exactly be reproduced.
"));
ASM1.EffluentSink Effluent annotation (extent=[88, -28; 108, -8]);
ASM1.SludgeSink WasteSludge annotation (extent=[87, -51; 107, -31]);
ASM1.SecClarModTakacs Settler annotation (extent=[48, -5; 68, 15]);
ASM1.divider2 divider annotation (extent=[20, -6; 40, 14]);
ASM1.nitri tank5(V=1333) annotation (extent=[-6, -6; 14, 14]);
ASM1.nitri tank4(V=1333) annotation (extent=[-32, -6; -12, 14]);
ASM1.nitri tank3(V=1333) annotation (extent=[-60, -6; -40, 14]);
ASM1.deni tank2 annotation (extent=[-48, 22; -28, 42]);
ASM1.deni tank1 annotation (extent=[-76, 22; -56, 42]);
ASM1.mixer3 mixer annotation (extent=[-104, 22; -84, 42]);
ModelicaAdditions.Tables.CombiTableTime CombiTableTime(
tableName="Inf_dry",
fileName="Inf_dry.txt",
icol={2,3,4,5,6,7,8,9,10,11,12,13,14,15}) annotation (extent=[-114, 78; -
94, 98]);
ASM1.WWSource WWSource annotation (extent=[-88, 78; -68, 98]);
ASM1.sensor_NO sensor_NO annotation (extent=[-42, 48; -22, 68]);
ASM1.blower blower1(Q_max=34574.2654508612) annotation (extent=[-60, -62; -
40, -42]);
ASM1.blower blower2(Q_max=34574.2654508612) annotation (extent=[-32, -62; -
12, -42]);
ASM1.blower blower3(Q_max=34574.2654508612) annotation (extent=[-6, -62; 14
, -42]);
ASM1.sensor_O2 sensor_O2 annotation (extent=[1, 25; 18, 42]);
Modelica.Blocks.Math.Feedback Feedback annotation (extent=[62, 40; 82, 60])
;
Modelica.Blocks.Continuous.PI PI1(k={500}, T={0.001}) annotation (extent=[
88, 40; 108, 60]);
Modelica.Blocks.Sources.Constant Constant1 annotation (extent=[-66, -89; -
46, -69]);
ASM1.pump RecyclePump(Q_max=92230) annotation (extent=[-94, -22; -74, -2],
rotation=180);
Modelica.Blocks.Continuous.PI PI2(k={15000}, T={0.05}) annotation (extent=[
58, 78; 78, 98]);
Modelica.Blocks.Sources.Constant NitrogenSetpoint annotation (extent=[-8,
78; 12, 98]);
Modelica.Blocks.Math.Feedback Feedback1 annotation (extent=[22, 78; 42, 98]
);
Modelica.Blocks.Nonlinear.Limiter Limiter1(uMax={10}, uMin={0.1})
annotation (extent=[8, 48; 28, 68]);
ASM1.pump ReturnPump(Q_max=18446) annotation (extent=[16, -36; 36, -16],
rotation=180);
ASM1.pump WastePump(Q_max=385) annotation (extent=[59, -55; 79, -35]);
Modelica.Blocks.Sources.Constant Constant2 annotation (extent=[22, -68; 42
, -48]);
Modelica.Blocks.Sources.Constant Temperature(k={15}) annotation (extent=[-
94, 50; -82, 62]);
ModelicaAdditions.Blocks.Nonlinear.FixedDelay FixedDelay1(delayTime=1/6/24)
annotation (extent=[-18, 48; 2, 68]);
sensor_NH sensor_NH1 annotation (extent=[64, 15; 80, 31]);
WasteWater.ASM1.sensor_NO sensor_NO1 annotation (extent=[81, 15; 97, 31]);
WasteWater.ASM1.sensor_TKN sensor_TKN1 annotation (extent=[97, 14; 113, 30]
);
WasteWater.ASM1.sensor_COD sensor_COD1 annotation (extent=[97, -5; 113, 11]
);
Modelica.Blocks.Sources.Step OxygenSetpoint(height={2}) annotation (extent=
[36, 40; 56, 60]);
WasteWater.ASM1.sensor_TSS sensor_TSS1 annotation (extent=[32, 15; 48, 30])
;
equation
connect(divider.Out1, Settler.Feed) annotation (points=[40, 6.6
; 44, 6.6; 44, 7.4; 48, 7.4]);
connect(tank5.Out, divider.In) annotation (points=[14, 4; 20, 4]);
connect(tank4.Out, tank5.In) annotation (points=[-12, 4; -6, 4]);
connect(tank3.Out, tank4.In) annotation (points=[-40, 4; -32, 4]);
connect(tank3.In, tank2.Out) annotation (points=[-60, 4; -70, 4; -70, 18; -
18, 18; -18, 32; -28, 32]);
connect(tank1.Out, tank2.In) annotation (points=[-56, 32; -48, 32]);
connect(mixer.Out, tank1.In) annotation (points=[-84, 32; -76, 32]);
connect(mixer.In1, WWSource.Out) annotation (points=[-104, 36; -104, 74; -
68, 74; -68, 81; -68.2, 81]);
connect(CombiTableTime.outPort, WWSource.data) annotation (points=[-93, 88
; -88, 88]);
connect(sensor_NO.In, tank2.MeasurePort) annotation (points=[-32, 48; -32,
36.5; -32.5, 36.5]);
connect(blower1.AirOut, tank3.AirIn) annotation (points=[-50, -42; -50, -
5.8]);
connect(blower2.AirOut, tank4.AirIn) annotation (points=[-22, -42; -22, -
5.8]);
connect(blower3.AirOut, tank5.AirIn) annotation (points=[4, -42; 4, -5.8]);
connect(Feedback.outPort, PI1.inPort) annotation (points=[81, 50; 86, 50]);
connect(PI1.outPort, blower3.u) annotation (points=[109, 50; 114, 50; 114,
-84; 18, -84; 18, -55; 14, -55; 14, -56]);
connect(divider.Out2, RecyclePump.In) annotation (points=[40, 2; 44, 2; 44
, -8.7; -74, -8.7]);
connect(RecyclePump.Out, mixer.In3) annotation (points=[-94.2, -14.8; -104
, -14.8; -104, 28]);
connect(Feedback1.outPort, PI2.inPort) annotation (points=[41, 88; 56, 88])
;
connect(NitrogenSetpoint.outPort, Feedback1.inPort1) annotation (points=[13
, 88; 24, 88]);
connect(PI2.outPort, RecyclePump.u) annotation (points=[79, 88; 118, 88;
118, -94; -70, -94; -70, -14.5; -75.1, -14.5]);
connect(Limiter1.outPort, Feedback1.inPort2) annotation (points=[29, 58; 32
, 58; 32, 80]);
connect(Settler.Return, ReturnPump.In) annotation (points=[55, -4.6; 55, -
22.7; 36, -22.7]);
connect(ReturnPump.Out, mixer.In2) annotation (points=[16, -28; 15.5, -28;
15.5, -30; -112, -30; -112, 32; -104, 32]);
connect(sensor_O2.So, Feedback.inPort2) annotation (points=[17.83, 34; 72,
34; 72, 42]);
connect(Temperature.outPort, tank1.T) annotation (points=[-81.5, 56; -78,
56; -78, 36; -75.8, 36], style(color=3));
connect(Temperature.outPort, tank2.T) annotation (points=[-81.5, 56; -52,
56; -52, 36; -48, 36], style(color=3));
connect(Temperature.outPort, tank3.T) annotation (points=[-81.5, 56; -52,
56; -52, 14; -60, 14; -60, 8; -59.5, 8.4; -59.5, 8.4; -60, 8], style(
color=3));
connect(Temperature.outPort, tank4.T) annotation (points=[-81.5, 56; -52,
56; -52, 14; -32, 14; -32, 8], style(color=3));
connect(Temperature.outPort, tank5.T) annotation (points=[-81.5, 56; -52,
56; -52, 14; -5.9, 14; -5.9, 8; -5.8, 8], style(color=3));
connect(sensor_NO.Sno, FixedDelay1.inPort) annotation (points=[-22, 58; -22
, 58], style(color=3));
connect(FixedDelay1.outPort, Limiter1.inPort) annotation (points=[3, 58; 4
, 58], style(color=3));
connect(OxygenSetpoint.outPort, Feedback.inPort1) annotation (points=[58,
50; 64, 50], style(color=3));
connect(Constant1.outPort, blower1.u) annotation (points=[-45, -79; -36, -
79; -36, -55; -40.2, -55], style(color=3));
connect(blower2.u, Constant1.outPort) annotation (points=[-12.2, -55; -9, -
55; -9, -79; -44, -79], style(color=3));
connect(WastePump.Out, WasteSludge.In) annotation (points=[79, -42; 87, -42
]);
connect(WastePump.In, Settler.Waste) annotation (points=[59, -48; 52, -48;
52, -31; 62, -31; 62, -3.6]);
connect(WastePump.u, Constant2.outPort) annotation (points=[60, -43; 46, -
43; 46, -58; 44, -58], style(color=3));
connect(sensor_NH1.In, Settler.Effluent) annotation (points=[72, 15; 72, 11
; 69, 11]);
connect(sensor_NO1.In, Settler.Effluent) annotation (points=[89, 15; 89, 11
; 68.2, 11]);
connect(sensor_TKN1.In, Settler.Effluent) annotation (points=[105, 14; 105
, 11; 68.2, 11]);
connect(sensor_COD1.In, Settler.Effluent) annotation (points=[105, -5; 105
, 11; 68.2, 11]);
connect(Effluent.In, Settler.Effluent) annotation (points=[88, -16; 78.5, -
16; 78.5, 11; 69, 11]);
connect(Constant2.outPort, ReturnPump.u) annotation (points=[43, -58; 46, -
58; 46, -29; 34.9, -29], style(color=3));
connect(tank5.MeasurePort, sensor_O2.In) annotation (points=[9.5, 8.5; 9.5
, 25; 9, 25]);
connect(sensor_TSS1.In, divider.Out1) annotation (points=[40, 15; 40, 7]);
end BenchPlant;
model ComplexPlant "Complex ASM1 WWTP"
import WasteWater;
extends Modelica.Icons.Example;
annotation (
Coordsys(
extent=[-200, -120; 200, 120],
grid=[1, 1],
component=[20, 20]),
Window(
x=0.45,
y=0.01,
width=0.35,
height=0.49),
Diagram,
Documentation(info="This ASM1 example plant configuration is from a real municipal wastewater treatment plant
with a size of 145.000 p.e. It is a cascade-type continuous flow plant for a mean dry
weather inflow of 28.500 m3/d. It consists of a preclarifier, an anaerobic tank,
3 denitrification and 3 nitrification tanks and a secondary settler.
This model is an example for the Wastewater library and is not adapted with its parameters
to the reality, therefore simulation results do not reflect the real plant behaviour.
Change into the directory ../ASM1 and translate the model.
Before simulating the model load initial values from the script file complex_asm1.mos
that is provided besides the model.
A 14 days dynamic influent data file is provided. So you may simulate up to 14 days.
But start with 1 day as it may take some time for simulation.
After simulation you may have a look at internal concentrations but most interesting
are the relevant concentrations at the effluent of a plant which can be viewed via the
sensors at the effluent of the secondary clarifier.
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
"));
ControlledDivider2 cdivider1 annotation (extent=[-168, 65; -148, 85]);
Modelica.Blocks.Sources.Constant Constant2(k={0.8}) annotation (extent=[-
178, 52; -168, 62]);
blower blower1(Q_max=162816)
"there exist 4 blowers of 4240 Nm3/h each, Q_max adusted according active aerated tanks"
annotation (extent=[145, -16; 165, 4]);
nitri nitri2(
V=2772,
alpha=0.305,
de=5.24,
R_air=20) annotation (extent=[110, 18; 130, 38]);
deni anaerob(V=1287) annotation (extent=[-138, 13; -118, 33]);
deni deni1(V=2772) annotation (extent=[-80, 14; -60, 34]);
deni deni3(V=2772) annotation (extent=[80, 18; 100, 38]);
deni deni2(V=2772) annotation (extent=[-20, 14; 0, 34]);
nitri nitri3(
V=5602,
alpha=0.305,
de=5.24,
R_air=21) annotation (extent=[144, 18; 164, 38]);
blower blower2(Q_max=81408)
"there exist 4 blowers of 4240 Nm3/h each, Q_max adjusted according active aerated tanks"
annotation (extent=[111, -13; 131, 7]);
pump ReturnPump(Q_max=60480) annotation (extent=[-54, -104; -34, -84],
rotation=180);
pump RecyclePump(Q_max=60480) annotation (extent=[0, -47; 20, -27],
rotation=180);
pump WastePump(Q_max=1920) annotation (extent=[128, -104; 148, -84]);
ControlledDivider2 cdivider2 annotation (extent=[-42, 68; -22, 88]);
EffluentSink Effluent annotation (extent=[170, -72; 190, -52]);
SludgeSink WasteSludge annotation (extent=[170, -100; 190, -80]);
mixer2 mixer2_1 annotation (extent=[-45, 16; -25, 36]);
mixer2 mixer2_2 annotation (extent=[50, 17; 70, 37]);
mixer3 mixer3_1 annotation (extent=[-107, 14; -87, 34]);
mixer2 mixer2_5 annotation (extent=[-140, -15; -120, 5]);
divider2 divider2_1 annotation (extent=[56, -42; 76, -22], rotation=180);
ControlledDivider2 cdivider3 annotation (extent=[-132, -97; -112, -77],
rotation=180);
Modelica.Blocks.Sources.Constant Constant4(k={0.5}) annotation (extent=[-
148, -80; -138, -70]);
nitri nitri1(
V=5602,
alpha=0.305,
de=5.24,
R_air=21) annotation (extent=[12, 14; 32, 34]);
blower blower3(Q_max=162816)
"there exist 4 blowers of max 4240 Nm3/h, Q_max adusted according active aerated tanks"
annotation (extent=[13, -23; 33, -3]);
Modelica.Blocks.Sources.Constant Constant7(k={0.56}) annotation (extent=[-
66, 98; -56, 108]);
mixer2 mixer2_3 annotation (extent=[-168, 18; -148, 38]);
PreClar.preclar3 Preclaryfier(V=1372, n_corr=2.138) annotation (extent=[-
136, 68; -116, 88]);
FlowSource FlowInput annotation (extent=[-176, 94; -156, 114]);
ControlledDivider2 ControlledDivider2_1 annotation (extent=[-51, -49; -31,
-29], rotation=-180);
Modelica.Blocks.Sources.Constant Constant6 annotation (extent=[-46, -10; -
36, 0]);
sensor_NO sensor_NO1 annotation (extent=[-14, 40; 6, 60]);
sensor_NO sensor_NO3 annotation (extent=[90, -56; 110, -36]);
sensor_NH sensor_NH2 annotation (extent=[112, -56; 132, -36]);
sensor_TSS sensor_TSS1 annotation (extent=[-6, -86; 14, -66]);
Modelica.Blocks.Sources.Constant Temperature(k={11.5}) annotation (extent=[
58, 44; 78, 64]);
WasteWater.Misc.RecycleController2 RecycleController1(NO3min=1.5)
annotation (extent=[15, -63; 29, -49]);
WasteWater.Misc.ReturnController ReturnController1 annotation (extent=[-28
, -115; -14, -101]);
sensor_Q sensor_Q1 annotation (extent=[-194, 66; -174, 86]);
WasteWater.Misc.TwoPoint TwoPoint1(
on=4.5,
off=4.0,
out_on=2.5,
out_off=1.5) annotation (extent=[-8, 100; 2, 110]);
WasteWater.Misc.TwoPoint TwoPoint2(
on=4.5,
off=4.0,
out_on=2.5,
out_off=1.5) annotation (extent=[110, 100; 120, 110]);
WasteWater.ASM1.sensor_NH sensor_NH1 annotation (extent=[142, 38; 162, 58])
;
Modelica.Blocks.Math.Feedback Feedback1 annotation (extent=[12, 100; 22,
110]);
Modelica.Blocks.Math.Feedback Feedback2 annotation (extent=[127, 99; 137,
109]);
WasteWater.ASM1.sensor_O2 sensor_O2_1 annotation (extent=[164, 38; 184, 58]
);
WasteWater.ASM1.sensor_O2 sensor_O2_2 annotation (extent=[14, 38; 34, 58]);
WasteWater.ASM1.sensor_COD sensor_COD1 annotation (extent=[134, -56; 154, -
36]);
WasteWater.ASM1.sensor_COD sensor_COD2 annotation (extent=[-118, 94; -98,
114]);
WasteWater.ASM1.Examples.JenaSecClarModTakacs Settler(hsc=3.46, Asc=3704)
"The depth is calculated based on V and A of the settler and not the true depth."
annotation (extent=[68, -84; 88, -64]);
WasteWater.ASM1.sensor_TKN sensor_TKN1 annotation (extent=[68, -56; 88, -36
]);
ModelicaAdditions.Tables.CombiTableTime CombiTableTime1(
tableName="drysim130303",
fileName="drysim130303.txt",
icol={2}) annotation (extent=[-198, 96; -184, 112]);
ModelicaAdditions.Tables.CombiTableTime CombiTableTime2(
tableName="drysim130303",
fileName="drysim130303.txt",
icol={3,7,4,6}) annotation (extent=[-134, 96; -118, 112]);
ModelicaAdditions.Blocks.Nonlinear.FixedDelay FixedDelay1(delayTime=1/24/6)
annotation (extent=[-16, -63; -2, -49]);
Modelica.Blocks.Math.Feedback Feedback3 annotation (extent=[89, -18; 79, -8
], rotation=-180);
WasteWater.ASM1.sensor_O2 sensor_O2_3 annotation (extent=[116, 40; 134, 58]
);
Modelica.Blocks.Math.Gain Gain1(k={500}) annotation (extent=[29, 95; 49,
115]);
Modelica.Blocks.Math.Gain Gain2(k={500}) annotation (extent=[144, 93; 164,
113]);
Modelica.Blocks.Math.Gain Gain3(k={500}) annotation (extent=[94, -28; 108,
-14]);
WasteWater.Misc.TwoPoint TwoPoint3(
on=4.5,
off=4.0,
out_on=2.0,
out_off=1.0) annotation (extent=[63, -19; 75, -7]);
Modelica.Blocks.Sources.Step Step1(
height={0.125},
offset={-0.5},
startTime={2.375}) annotation (extent=[100, -90; 110, -80]);
WasteWater.ASM1.sensor_TSS sensor_TSS2 annotation (extent=[166, 5; 184, 24]
);
equation
connect(deni3.Out, nitri2.In) annotation (points=[100, 28; 110,
28]);
connect(nitri2.Out, nitri3.In) annotation (points=[130, 28; 144, 28]);
connect(nitri3.Out, divider2_1.In) annotation (points=[164, 28; 190, 28;
190, -32.3; 76, -32.3]);
connect(anaerob.Out, mixer3_1.In2) annotation (points=[-118, 23; -110.5, 23
; -110.5, 24.5; -107, 24.5]);
connect(mixer2_5.Out, mixer3_1.In3) annotation (points=[-120, -4.6; -111, -
4.6; -111, 20.5; -107, 20.5]);
connect(Constant4.outPort, cdivider3.u) annotation (points=[-137.5, -75; -
122, -75; -122, -81], style(color=3));
connect(nitri1.Out, mixer2_2.In2) annotation (points=[33, 24; 39, 24; 39,
25.5; 51, 25.5]);
connect(deni2.Out, nitri1.In) annotation (points=[5.55112e-16, 24; 12, 24])
;
connect(cdivider2.Out1, mixer2_2.In1) annotation (points=[-22, 80.5; 40,
80.5; 40, 29.5; 50, 29.5]);
connect(cdivider3.Out2, mixer2_5.In1) annotation (points=[-132, -85.5; -172
, -85.5; -172, -2.5; -140, -2.5]);
connect(deni1.Out, mixer2_1.In1) annotation (points=[-60, 24; -56, 24; -56
, 29; -50.75, 29; -50.75, 28.5; -45, 28.5]);
connect(Constant7.outPort, cdivider2.u) annotation (points=[-55.5, 103; -
50.5, 103; -50.5, 54; -32, 54; -32, 72], style(color=3));
connect(cdivider3.In, ReturnPump.Out) annotation (points=[-112, -87.3; -
86.1, -87.3; -86.1, -96.8; -54.2, -96.8]);
connect(cdivider3.Out1, mixer2_3.In2) annotation (points=[-132, -89.5; -184
, -89.5; -184, 26.5; -168, 26.5]);
connect(cdivider1.Out2, mixer2_3.In1) annotation (points=[-147, 73.5; -143
, 73.5; -143, 40; -184, 40; -184, 30.5; -168, 30.5]);
connect(WastePump.Out, WasteSludge.In) annotation (points=[148.2, -91.2;
170, -91.2]);
connect(ControlledDivider2_1.Out2, mixer2_1.In2) annotation (points=[-51, -
38.5; -52, -38.5; -52, 24.5; -45, 24.5]);
connect(ControlledDivider2_1.Out1, mixer2_5.In2) annotation (points=[-51, -
42.5; -150, -42.5; -150, -6.5; -140, -6.5]);
connect(Constant6.outPort, ControlledDivider2_1.u) annotation (points=[-36
, -5; -31, -5; -31, -18; -41, -18; -41, -34], style(color=3));
connect(Preclaryfier.In, cdivider1.Out1) annotation (points=[-136, 78; -
146.5, 78; -146.5, 77.5; -147, 77.5]);
connect(mixer2_2.Out, deni3.In) annotation (points=[70, 27.4; 70, 28; 80,
28]);
connect(sensor_NO1.In, deni2.MeasurePort) annotation (points=[-4, 40; -4,
28.5; -4.5, 28.5]);
connect(Temperature.outPort, anaerob.T) annotation (points=[80, 54; 94, 54
; 94, 34; -138, 34; -138, 27], style(color=3));
connect(Temperature.outPort, deni1.T) annotation (points=[80, 54; 94, 54;
94, 34; -80, 34; -80, 28], style(color=3));
connect(Temperature.outPort, deni2.T) annotation (points=[80, 54; 94, 54;
94, 34; -20, 34; -20, 28], style(color=3));
connect(Temperature.outPort, nitri1.T) annotation (points=[80, 54; 94, 54;
94, 34; 12, 34; 12, 28; 12.2, 28], style(color=3));
connect(Temperature.outPort, deni3.T) annotation (points=[80, 54; 94, 54;
94, 34; 80, 34; 80, 32], style(color=3));
connect(Temperature.outPort, nitri2.T) annotation (points=[80, 54; 94, 54;
94, 34; 110.5, 34; 110.5, 32; 110, 32], style(color=3));
connect(Temperature.outPort, nitri3.T) annotation (points=[80, 54; 94, 54;
94, 34; 144, 34; 144, 32], style(color=3));
connect(RecycleController1.out, RecyclePump.u) annotation (points=[29.7, -
56; 34, -56; 34, -38.5; 18.8, -38.5], style(color=3));
connect(FlowInput.Out, sensor_Q1.In) annotation (points=[-156, 97; -148, 97
; -148, 90; -194, 90; -194, 76]);
connect(sensor_Q1.Q, ReturnController1.in1) annotation (points=[-184, 66; -
184, 65.5; -194, 65.5; -194, -108; -28.7, -108], style(color=3));
connect(sensor_NH1.In, nitri3.MeasurePort) annotation (points=[152, 38; 152
, 34; 160, 34; 160, 32.5]);
connect(sensor_O2_2.In, nitri1.MeasurePort) annotation (points=[24, 38; 24
, 28.5; 28, 28.5]);
connect(Feedback1.inPort2, sensor_O2_2.So) annotation (points=[17, 101; 17
, 83.5; 34, 83.5; 34, 48], style(color=3));
connect(sensor_O2_1.In, nitri3.MeasurePort) annotation (points=[174, 38;
174, 34; 160.25, 34; 160.25, 32; 160, 32]);
connect(Feedback2.inPort2, sensor_O2_1.So) annotation (points=[132, 101;
132, 81; 184, 81; 184, 48], style(color=3));
connect(sensor_NH1.Snh, TwoPoint1.e) annotation (points=[162, 48; 162, 72;
-14, 72; -14, 105; -8, 105], style(color=3));
connect(TwoPoint2.e, sensor_NH1.Snh) annotation (points=[110, 105; 100, 105
; 100, 72; 162, 72; 162, 48], style(color=3));
connect(sensor_COD2.In, Preclaryfier.In) annotation (points=[-108, 94; -108
, 90; -136, 90; -136, 78]);
connect(Settler.Effluent, Effluent.In) annotation (points=[88.2, -68.3;
108.65, -68.3; 108.65, -68; 128, -68; 128, -60; 170, -60]);
connect(sensor_NO3.In, Settler.Effluent) annotation (points=[100, -56; 100
, -68.3; 88.2, -68.3]);
connect(sensor_NH2.In, Settler.Effluent) annotation (points=[122, -56; 122
, -68.3; 88.2, -68.3]);
connect(sensor_COD1.In, Settler.Effluent) annotation (points=[144, -56; 144
, -60; 128, -60; 128, -68.3; 88.2, -68.3]);
connect(WastePump.In, Settler.Waste) annotation (points=[128, -97.3; 82, -
97.3; 82, -84]);
connect(Settler.Return, ReturnPump.In) annotation (points=[74, -84; 74, -
90.7; -34, -90.7]);
connect(sensor_TSS1.In, Settler.Return) annotation (points=[4, -86; 4, -91
; 74, -91; 74, -84]);
connect(sensor_TKN1.In, Settler.Effluent) annotation (points=[78, -56; 88,
-56; 88, -68]);
connect(CombiTableTime1.outPort, FlowInput.data) annotation (points=[-183.3
, 104; -176, 104], style(color=3));
connect(Preclaryfier.MeasurePort, CombiTableTime2.outPort) annotation (
points=[-122, 88; -122, 94; -116, 94; -116, 104; -117.2, 104], style(
color=3));
connect(FixedDelay1.inPort, sensor_NO1.Sno) annotation (points=[-18, -56; -
22, -56; -22, -8; 10, -8; 10, 50; 6, 50], style(color=3));
connect(sensor_O2_3.In, nitri2.MeasurePort) annotation (points=[125, 40;
125, 32.5; 126, 32.5]);
connect(sensor_Q1.Out, cdivider1.In) annotation (points=[-174, 76; -168, 76
]);
connect(cdivider2.Out2, mixer3_1.In1) annotation (points=[-22, 76; -18, 76
; -18, 47; -111, 47; -111, 28.5; -107, 28.5]);
connect(mixer3_1.Out, deni1.In) annotation (points=[-87, 24; -81, 24]);
connect(mixer2_1.Out, deni2.In) annotation (points=[-24, 27; -22.5, 27; -
22.5, 23; -21, 23]);
connect(FixedDelay1.outPort, RecycleController1.in1) annotation (points=[-1
, -56; 14.3, -56], style(color=3));
connect(Feedback2.inPort1, TwoPoint2.u) annotation (points=[128, 105; 120,
105], style(color=3));
connect(TwoPoint1.u, Feedback1.inPort1) annotation (points=[2.5, 105; 13,
105], style(color=3));
connect(mixer2_3.Out, anaerob.In) annotation (points=[-148, 29; -143, 29; -
143, 23; -138, 23]);
connect(Gain1.inPort, Feedback1.outPort) annotation (points=[27, 104; 24,
104; 24, 105; 21, 105], style(color=3));
connect(Gain2.inPort, Feedback2.outPort) annotation (points=[141, 103; 139
, 103; 139, 104; 137, 104], style(color=3));
connect(Feedback3.outPort, Gain3.inPort) annotation (points=[88.5, -13;
90.55, -13; 90.55, -21; 92.6, -21], style(color=3));
connect(sensor_O2_3.So, Feedback3.inPort2) annotation (points=[134, 49; 139
, 49; 139, 16; 83, 16; 83, -9], style(color=3));
connect(TwoPoint3.u, Feedback3.inPort1) annotation (points=[75, -13; 80, -
13], style(color=3));
connect(TwoPoint3.e, sensor_NH1.Snh) annotation (points=[63, -13; 56, -13;
56, 6; 104, 6; 104, 72; 162, 72; 162, 49], style(color=3));
connect(Step1.outPort, WastePump.u) annotation (points=[111, -85; 119, -85
; 119, -91; 129.2, -91], style(color=3));
connect(blower1.AirOut, nitri3.AirIn) annotation (points=[155, 4; 155, 5.5
; 154, 5.5; 154, 18]);
connect(Constant2.outPort, cdivider1.u) annotation (points=[-167.5, 57; -
158, 57; -158, 69], style(color=3));
connect(Gain1.outPort, blower3.u) annotation (points=[50, 105; 55, 105; 55
, 83; 49, 83; 49, -15; 32, -15], style(color=3));
connect(Gain3.outPort, blower2.u) annotation (points=[108.7, -21; 139, -21
; 139, -6; 130.8, -6], style(color=3));
connect(blower1.u, Gain2.outPort) annotation (points=[165.8, -9; 195, -9;
195, 103; 165, 103], style(color=3));
connect(Settler.Feed, divider2_1.Out1) annotation (points=[68, -72; 48, -72
; 48, -35; 55, -35]);
connect(ReturnPump.u, ReturnController1.out) annotation (points=[-35.1, -97
; -5, -97; -5, -108; -13, -108], style(color=3));
connect(RecyclePump.Out, ControlledDivider2_1.In) annotation (points=[-
8.88178e-16, -40; -30, -40]);
connect(RecyclePump.In, divider2_1.Out2) annotation (points=[20, -33.7; 38
, -33.7; 38, -31; 56, -31]);
connect(blower3.AirOut, nitri1.AirIn) annotation (points=[22, -3; 22, 14]);
connect(cdivider2.In, Preclaryfier.Out) annotation (points=[-43, 78; -115,
78]);
connect(blower2.AirOut, nitri2.AirIn) annotation (points=[120, 7; 120, 18])
;
connect(sensor_TSS2.In, nitri3.Out) annotation (points=[175, 5; 175, 28;
163, 28]);
end ComplexPlant;
model JenaSecClarModTakacs
"Secondary Clarifier Model based on Takacs prepared for a special plant"
extends WasteWater.Icons.SecClar;
extends ASM1.SecClar.Takacs.Interfaces.ratios;
package SCP = ASM1.SecClar.Takacs;
package SI = Modelica.SIunits;
package WI = ASM1.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;
// layers 1 to 10
WI.WWFlowAsm1in Feed annotation (extent=[-110, 4; -90, 24]);
WI.WWFlowAsm1out Effluent annotation (extent=[92, 47; 112, 67]);
WI.WWFlowAsm1out Return annotation (extent=[-40, -106; -20, -86]);
WI.WWFlowAsm1out Waste annotation (extent=[20, -106; 40, -86]);
SCP.bottom_layer S1(
zm=zm,
Asc=Asc,
Xf=Xf,
rXs=rXs,
rXbh=rXbh,
rXba=rXba,
rXp=rXp,
rXi=rXi,
rXnd=rXnd) 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.feed_layer S4(
zm=zm,
Asc=Asc,
Xf=Xf) annotation (extent=[-36, -36; 34, -22]);
SCP.upper_layer S5(
zm=zm,
Asc=Asc,
Xf=Xf,
Xt=Xt) annotation (extent=[-36, -16; 34, -4]);
SCP.upper_layer S6(
zm=zm,
Asc=Asc,
Xf=Xf,
Xt=Xt) annotation (extent=[-36, 2; 34, 16]);
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,
rXs=rXs,
rXbh=rXbh,
rXba=rXba,
rXp=rXp,
rXi=rXi,
rXnd=rXnd) 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(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.45,
y=0.01,
width=0.35,
height=0.49),
Documentation(info="This component models an ASM1 10 - layer secondary clarifier with 6 layers above the feed_layer (including top_layer)
and 3 layers below the feed_layer (including bottom_layer).
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(S8.Up, S9.Dn) annotation (points=[-2.22045e-15, 55; -2.22045e-15,
59]);
connect(S1.PQw, Waste) annotation (points=[17.5, -93; 30, -93; 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; -30, -93; -30, -100])
;
connect(S4.Dn, S3.Up) annotation (points=[0, -36; 0, -40]);
connect(S4.Up, S5.Dn) annotation (points=[-2, -22; -2, -16]);
connect(S5.Up, S6.Dn) annotation (points=[0, -4; 0, 2]);
connect(S6.Up, S7.Dn) annotation (points=[0, 16; 0, 21]);
connect(Feed, S4.In) annotation (points=[-100, 10; -67.5, 10; -67.5, -28.72
; -35, -28.72]);
// total sludge concentration in clarifier feed
Xf = 0.75*(Feed.Xs + Feed.Xbh + Feed.Xba + Feed.Xp + Feed.Xi);
// ratios of solid components
rXs = Feed.Xs/Xf;
rXbh = Feed.Xbh/Xf;
rXba = Feed.Xba/Xf;
rXp = Feed.Xp/Xf;
rXi = Feed.Xi/Xf;
rXnd = Feed.Xnd/Xf;
end JenaSecClarModTakacs;
end Examples;
Click here to get the file