Skip to content

Commit

Permalink
adding test networks from msx repo
Browse files Browse the repository at this point in the history
  • Loading branch information
@samhatchett
samhatchett committed Jul 14, 2016
1 parent 2ea0f45 commit b35f1f5
Show file tree
Hide file tree
Showing 35 changed files with 4,600 additions and 0 deletions.
59 changes: 59 additions & 0 deletions As5Adsorb/example-gc.msx
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
[TITLE]
Arsenic Oxidation/Adsorption Example

[OPTIONS]
AREA_UNITS M2 ;Surface concentration is mass/m2
RATE_UNITS HR ;Reaction rates are concentration/hour
SOLVER RK5 ;5-th order Runge-Kutta integrator
TIMESTEP 360 ;360 sec (5 min) solution time step
RTOL 0.001 ;Relative concentration tolerance
ATOL 0.0001 ;Absolute concentration tolerance
COMPILER GC

[SPECIES]
BULK AS3 UG ;Dissolved arsenite
BULK AS5 UG ;Dissolved arsenate
BULK AStot UG ;Total dissolved arsenic
WALL AS5s UG ;Adsorbed arsenate
BULK NH2CL MG ;Monochloramine

[COEFFICIENTS]
CONSTANT Ka 10.0 ;Arsenite oxidation rate coefficient
CONSTANT Kb 0.1 ;Monochloramine decay rate coefficient
CONSTANT K1 5.0 ;Arsenate adsorption coefficient
CONSTANT K2 1.0 ;Arsenate desorption coefficient
CONSTANT Smax 50 ;Arsenate adsorption saturation limit

[TERMS]
Ks K1/K2 ;Equil. adsorption coeff.

[PIPES]
;Arsenite oxidation
RATE AS3 -Ka*AS3*NH2CL
;Arsenate production
RATE AS5 Ka*AS3*NH2CL - Av*(K1*(Smax-AS5s)*AS5 - K2*AS5s)
;Monochloramine decay
RATE NH2CL -Kb*NH2CL
;Arsenate adsorption
EQUIL AS5s Ks*Smax*AS5/(1+Ks*AS5) - AS5s
;Total bulk arsenic
FORMULA AStot AS3 + AS5

[TANKS]
RATE AS3 -Ka*AS3*NH2CL
RATE AS5 Ka*AS3*NH2CL
RATE NH2CL -Kb*NH2CL
FORMULA AStot AS3 + AS5

[QUALITY]
;Initial conditions (= 0 if not specified here)
NODE Source AS3 10.0
NODE Source NH2CL 2.5

[REPORT]
NODES C D ;Report results for nodes C and D
LINKS 5 ;Report results for pipe 5
SPECIE AStot YES ;Report results for each specie
SPECIE AS5 YES
SPECIE AS5s YES
SPECIE NH2CL YES
59 changes: 59 additions & 0 deletions As5Adsorb/example-vc.msx
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
[TITLE]
Arsenic Oxidation/Adsorption Example

[OPTIONS]
AREA_UNITS M2 ;Surface concentration is mass/m2
RATE_UNITS HR ;Reaction rates are concentration/hour
SOLVER RK5 ;5-th order Runge-Kutta integrator
TIMESTEP 360 ;360 sec (5 min) solution time step
RTOL 0.001 ;Relative concentration tolerance
ATOL 0.0001 ;Absolute concentration tolerance
COMPILER VC

[SPECIES]
BULK AS3 UG ;Dissolved arsenite
BULK AS5 UG ;Dissolved arsenate
BULK AStot UG ;Total dissolved arsenic
WALL AS5s UG ;Adsorbed arsenate
BULK NH2CL MG ;Monochloramine

[COEFFICIENTS]
CONSTANT Ka 10.0 ;Arsenite oxidation rate coefficient
CONSTANT Kb 0.1 ;Monochloramine decay rate coefficient
CONSTANT K1 5.0 ;Arsenate adsorption coefficient
CONSTANT K2 1.0 ;Arsenate desorption coefficient
CONSTANT Smax 50 ;Arsenate adsorption saturation limit

[TERMS]
Ks K1/K2 ;Equil. adsorption coeff.

[PIPES]
;Arsenite oxidation
RATE AS3 -Ka*AS3*NH2CL
;Arsenate production
RATE AS5 Ka*AS3*NH2CL - Av*(K1*(Smax-AS5s)*AS5 - K2*AS5s)
;Monochloramine decay
RATE NH2CL -Kb*NH2CL
;Arsenate adsorption
EQUIL AS5s Ks*Smax*AS5/(1+Ks*AS5) - AS5s
;Total bulk arsenic
FORMULA AStot AS3 + AS5

[TANKS]
RATE AS3 -Ka*AS3*NH2CL
RATE AS5 Ka*AS3*NH2CL
RATE NH2CL -Kb*NH2CL
FORMULA AStot AS3 + AS5

[QUALITY]
;Initial conditions (= 0 if not specified here)
NODE Source AS3 10.0
NODE Source NH2CL 2.5

[REPORT]
NODES C D ;Report results for nodes C and D
LINKS 5 ;Report results for pipe 5
SPECIE AStot YES ;Report results for each specie
SPECIE AS5 YES
SPECIE AS5s YES
SPECIE NH2CL YES
36 changes: 36 additions & 0 deletions As5Adsorb/example.inp
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
[TITLE]
EPANET-MSX Example Network

[JUNCTIONS]
;ID Elev Demand Pattern
A 0 4.1
B 0 3.4
C 0 5.5
D 0 2.3

[RESERVOIRS]
;ID Head Pattern
Source 100

[PIPES]
;ID Node1 Node2 Length Diameter Roughness
1 Source A 1000 200 100
2 A B 800 150 100
3 A C 1200 200 100
4 B C 1000 150 100
5 C D 2000 150 100

[TIMES]
Duration 48
Hydraulic Timestep 1:00
Quality Timestep 0:05
Report Timestep 2
Report Start 0
Statistic NONE

[OPTIONS]
Units CMH
Headloss H-W
Quality NONE

[END]
59 changes: 59 additions & 0 deletions As5Adsorb/example.msx
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
[TITLE]
Arsenic Oxidation/Adsorption Example

[OPTIONS]
AREA_UNITS M2 ;Surface concentration is mass/m2
RATE_UNITS HR ;Reaction rates are concentration/hour
SOLVER RK5 ;5-th order Runge-Kutta integrator
TIMESTEP 360 ;360 sec (5 min) solution time step
RTOL 0.001 ;Relative concentration tolerance
ATOL 0.0001 ;Absolute concentration tolerance
COMPILER NONE

[SPECIES]
BULK AS3 UG ;Dissolved arsenite
BULK AS5 UG ;Dissolved arsenate
BULK AStot UG ;Total dissolved arsenic
WALL AS5s UG ;Adsorbed arsenate
BULK NH2CL MG ;Monochloramine

[COEFFICIENTS]
CONSTANT Ka 10.0 ;Arsenite oxidation rate coefficient
CONSTANT Kb 0.1 ;Monochloramine decay rate coefficient
CONSTANT K1 5.0 ;Arsenate adsorption coefficient
CONSTANT K2 1.0 ;Arsenate desorption coefficient
CONSTANT Smax 50 ;Arsenate adsorption saturation limit

[TERMS]
Ks K1/K2 ;Equil. adsorption coeff.

[PIPES]
;Arsenite oxidation
RATE AS3 -Ka*AS3*NH2CL
;Arsenate production
RATE AS5 Ka*AS3*NH2CL - Av*(K1*(Smax-AS5s)*AS5 - K2*AS5s)
;Monochloramine decay
RATE NH2CL -Kb*NH2CL
;Arsenate adsorption
EQUIL AS5s Ks*Smax*AS5/(1+Ks*AS5) - AS5s
;Total bulk arsenic
FORMULA AStot AS3 + AS5

[TANKS]
RATE AS3 -Ka*AS3*NH2CL
RATE AS5 Ka*AS3*NH2CL
RATE NH2CL -Kb*NH2CL
FORMULA AStot AS3 + AS5

[QUALITY]
;Initial conditions (= 0 if not specified here)
NODE Source AS3 10.0
NODE Source NH2CL 2.5

[REPORT]
NODES C D ;Report results for nodes C and D
LINKS 5 ;Report results for pipe 5
SPECIE AStot YES ;Report results for each specie
SPECIE AS5 YES
SPECIE AS5s YES
SPECIE NH2CL YES
146 changes: 146 additions & 0 deletions As5Adsorb/example.rpt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,146 @@
Page 1 Mon Aug 10 13:53:35 2009

******************************************************************
* E P A N E T *
* Hydraulic and Water Quality *
* Analysis for Pipe Networks *
* Version 2.00.12 *
******************************************************************

EPANET-MSX Example Network

Input Data File ................... As5Adsorb/example.inp
Number of Junctions................ 4
Number of Reservoirs............... 1
Number of Tanks ................... 0
Number of Pipes ................... 5
Number of Pumps ................... 0
Number of Valves .................. 0
Headloss Formula .................. Hazen-Williams
Hydraulic Timestep ................ 1.00 hrs
Hydraulic Accuracy ................ 0.001000
Status Check Frequency ............ 2
Maximum Trials Checked ............ 10
Damping Limit Threshold ........... 0.000000
Maximum Trials .................... 200
Quality Analysis .................. None
Specific Gravity .................. 1.00
Relative Kinematic Viscosity ...... 1.00
Relative Chemical Diffusivity ..... 1.00
Demand Multiplier ................. 1.00
Total Duration .................... 48.00 hrs
Reporting Criteria:
No Nodes
No Links

Analysis begun Mon Aug 10 13:53:36 2009

Processing MSX input file As5Adsorb/example.msx


Page 1 EPANET-MSX 1.1

******************************************************************
* E P A N E T - M S X *
* Multi-Species Water Quality *
* Analysis for Pipe Networks *
* Version 1.1 *
******************************************************************

Arsenic Oxidation/Adsorption Example


<<< Node C >>>

Time AS5 AStot NH2CL
hr:min UG/L UG/L MG/L
------- ---------- ---------- ----------
0:00 0.00 0.00 0.00
2:00 0.00 0.00 0.00
4:00 0.00 0.00 0.00
6:00 0.00 0.00 0.00
8:00 0.00 0.00 1.10
10:00 9.17 9.17 1.10
12:00 9.17 9.17 1.10
14:00 9.17 9.17 1.10
16:00 9.17 9.17 1.10
18:00 9.17 9.17 1.10
20:00 9.17 9.17 1.10
22:00 9.17 9.17 1.10
24:00 9.17 9.17 1.10
26:00 9.17 9.17 1.10
28:00 9.17 9.17 1.10
30:00 9.17 9.17 1.10
32:00 9.17 9.17 1.10
34:00 9.17 9.17 1.11
36:00 9.17 9.17 1.11
38:00 10.03 10.03 1.11
40:00 10.03 10.03 1.11
42:00 10.03 10.03 1.11
44:00 10.03 10.03 1.11
46:00 10.03 10.03 1.11
48:00 10.03 10.03 1.11

<<< Node D >>>

Time AS5 AStot NH2CL
hr:min UG/L UG/L MG/L
------- ---------- ---------- ----------
0:00 0.00 0.00 0.00
2:00 0.00 0.00 0.00
4:00 0.00 0.00 0.00
6:00 0.00 0.00 0.00
8:00 0.00 0.00 0.00
10:00 0.00 0.00 0.00
12:00 0.00 0.00 0.00
14:00 0.00 0.00 0.00
16:00 0.00 0.00 0.00
18:00 0.00 0.00 0.00
20:00 0.00 0.00 0.00
22:00 0.00 0.00 0.00
24:00 0.00 0.00 0.24
26:00 0.00 0.00 0.24
28:00 9.17 9.17 0.24
30:00 9.17 9.17 0.24
32:00 9.17 9.17 0.24
34:00 9.17 9.17 0.24
36:00 9.17 9.17 0.24
38:00 9.17 9.17 0.24
40:00 9.17 9.17 0.24
42:00 9.17 9.17 0.24
44:00 9.17 9.17 0.24
46:00 9.17 9.17 0.24
48:00 9.17 9.17 0.24

<<< Link 5 >>>

Time AS5 AStot AS5s NH2CL
hr:min UG/L UG/L UG/M2 MG/L
------- ---------- ---------- ---------- ----------
0:00 0.00 0.00 0.00 0.00
2:00 0.00 0.00 0.00 0.00
4:00 0.00 0.00 0.00 0.00
6:00 0.00 0.00 0.00 0.00
8:00 0.00 0.00 0.00 0.05
10:00 0.85 0.85 4.51 0.17
12:00 1.86 1.86 9.88 0.27
14:00 2.87 2.87 15.27 0.35
16:00 3.87 3.87 20.64 0.42
18:00 4.88 4.88 26.02 0.47
20:00 5.89 5.89 31.39 0.52
22:00 6.89 6.89 36.75 0.55
24:00 7.90 7.90 42.12 0.56
26:00 8.91 8.91 47.50 0.56
28:00 9.17 9.17 48.93 0.56
30:00 9.17 9.17 48.93 0.56
32:00 9.17 9.17 48.93 0.56
34:00 9.17 9.17 48.93 0.56
36:00 9.17 9.17 48.93 0.57
38:00 9.19 9.19 48.93 0.57
40:00 9.30 9.30 48.95 0.57
42:00 9.41 9.41 48.96 0.57
44:00 9.52 9.52 48.97 0.57
46:00 9.64 9.64 48.98 0.57
48:00 9.75 9.75 48.99 0.57

Analysis ended Mon Aug 10 13:53:37 2009
13 changes: 13 additions & 0 deletions As5Adsorb/readme.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
Arsenite/Arsenate adsorption/desorption example
==================================================

This example models arsenite/arsenate oxidation by
monochloramine, with adsorption of arsenate to the
pipe wall according to the local equilibrium
assumption. The example network is the simple 4-node
network described in section 3 of the EPANET-MSX user's manual.
The files included are:

example.inp -- the EPANET network input file.
example.msx -- the MSX reaction file
example.rpt -- the report file containing the MSX results
Loading

0 comments on commit b35f1f5

Please sign in to comment.