Skip to content

Energy Surrogate Model #52

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
jakechurchi wants to merge 12 commits into
base: main
Choose a base branch
from
Draft

Energy Surrogate Model #52

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
7085849
Initial commit with parameter sweep components
May 6, 2026
bd312fe
Working param sweep function. Issues with finding solutions remain
May 7, 2026
11f7696
Initial results from working version of param sweep. Unclear if recov…
May 8, 2026
0d32366
Added uf pump energy sweep
May 8, 2026
212c452
Now ro_sweep actual works
May 8, 2026
42f2960
Results with different operational limits. Min brine flow instead of …
May 8, 2026
18b4092
Created sample energy surrogate
May 8, 2026
f22e56d
Updated surrogates
May 11, 2026
019d4fa
black and deleting unneeded files
May 11, 2026
011442c
Param sweeps is working and pump efficiency has been updated accordin…
May 13, 2026
1291413
Created new surrogate
May 14, 2026
f6a21a4
black .
May 19, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Binary file added ps_demo.h5
View file
Open in desktop
Binary file not shown.
70 changes: 70 additions & 0 deletions ro_sweep_brine_lim(1).csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
# Recovery,Feed Flow,Stage1 RR,Stage2 RR,Stage3 RR,Stage1 Qin (m3/s),Stage2 Qin (m3/s),Stage3 Qin (m3/s),Stage1 Brine (m3/s),Stage2 Brine (m3/s),Stage3 Brine (m3/s),Stage1 Pin (Pa),Stage2 Pin (Pa),Stage3 Pin (Pa),Stage1 Pout (Pa),Stage2 Pout (Pa),Stage3 Pout (Pa),Total Power (W),Stage1 Power (W),Stage2 Power (W),Stage3 Power (W),Stage1 First Elem. Prod. (m3/s),Stage2 First Elem. Prod. (m3/s),Stage3 First Elem. Prod. (m3/s),Stage1 Perm. Conc. (kg/m3),Stage2 Perm. Conc. (kg/m3),Stage3 Perm. Conc. (kg/m3)
8.80E-01,1.45E-01,5.84E-01,5.84E-01,3.07E-01,1.45E-01,6.03E-02,2.51E-02,6.03E-02,2.51E-02,1.74E-02,2.44E+05,7.89E+05,6.33E+05,8.63E+05,9.63E+05,6.57E+05,1.57E+05,1.42E+05,1.40E+04,1.23E+03,1.83E-04,1.86E-04,9.08E-05,1.39E-02,2.09E-02,2.20E-01
8.80E-01,1.51E-01,5.84E-01,5.84E-01,3.05E-01,1.51E-01,6.28E-02,2.61E-02,6.28E-02,2.61E-02,1.81E-02,2.44E+05,8.16E+05,6.61E+05,8.90E+05,9.91E+05,6.66E+05,1.69E+05,1.54E+05,1.47E+04,2.88E+02,1.90E-04,1.93E-04,9.33E-05,1.34E-02,2.02E-02,2.14E-01
8.80E-01,1.57E-01,5.85E-01,5.85E-01,3.03E-01,1.57E-01,6.53E-02,2.71E-02,6.53E-02,2.71E-02,1.89E-02,2.44E+05,8.44E+05,6.89E+05,9.17E+05,1.02E+06,6.75E+05,1.81E+05,1.66E+05,1.54E+04,-7.47E+02,1.98E-04,2.01E-04,9.57E-05,1.29E-02,1.95E-02,2.09E-01
8.80E-01,1.64E-01,5.86E-01,5.86E-01,3.01E-01,1.64E-01,6.78E-02,2.81E-02,6.78E-02,2.81E-02,1.96E-02,2.44E+05,8.71E+05,7.17E+05,9.45E+05,1.05E+06,6.83E+05,1.93E+05,1.79E+05,1.61E+04,-1.88E+03,2.05E-04,2.08E-04,9.80E-05,1.24E-02,1.88E-02,2.04E-01
8.80E-01,1.70E-01,5.86E-01,5.86E-01,2.99E-01,1.70E-01,7.03E-02,2.91E-02,7.03E-02,2.91E-02,2.04E-02,2.44E+05,8.99E+05,7.45E+05,9.73E+05,1.08E+06,6.91E+05,2.06E+05,1.92E+05,1.69E+04,-3.11E+03,2.13E-04,2.15E-04,1.00E-04,1.20E-02,1.82E-02,1.99E-01
8.80E-01,1.76E-01,5.87E-01,5.87E-01,2.96E-01,1.76E-01,7.27E-02,3.00E-02,7.27E-02,3.00E-02,2.11E-02,2.44E+05,9.27E+05,7.73E+05,1.00E+06,1.10E+06,6.99E+05,2.20E+05,2.06E+05,1.77E+04,-4.43E+03,2.20E-04,2.22E-04,1.02E-04,1.16E-02,1.77E-02,1.95E-01
8.89E-01,1.45E-01,5.85E-01,5.85E-01,3.56E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.61E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,7.09E+05,1.60E+05,1.42E+05,1.39E+04,3.83E+03,1.83E-04,1.86E-04,1.06E-04,1.39E-02,2.10E-02,2.02E-01
8.89E-01,1.51E-01,5.88E-01,5.88E-01,3.46E-01,1.51E-01,6.23E-02,2.57E-02,6.23E-02,2.57E-02,1.68E-02,2.44E+05,8.20E+05,6.61E+05,8.94E+05,9.91E+05,7.11E+05,1.72E+05,1.55E+05,1.43E+04,2.56E+03,1.91E-04,1.93E-04,1.05E-04,1.34E-02,2.05E-02,2.04E-01
8.89E-01,1.57E-01,5.89E-01,5.89E-01,3.44E-01,1.57E-01,6.48E-02,2.67E-02,6.48E-02,2.67E-02,1.75E-02,2.44E+05,8.48E+05,6.89E+05,9.22E+05,1.02E+06,7.21E+05,1.84E+05,1.67E+05,1.50E+04,1.70E+03,1.99E-04,2.00E-04,1.08E-04,1.29E-02,1.98E-02,1.99E-01
8.89E-01,1.64E-01,5.89E-01,5.89E-01,3.42E-01,1.64E-01,6.72E-02,2.76E-02,6.72E-02,2.76E-02,1.82E-02,2.44E+05,8.76E+05,7.17E+05,9.50E+05,1.05E+06,7.31E+05,1.96E+05,1.80E+05,1.56E+04,7.41E+02,2.06E-04,2.08E-04,1.10E-04,1.24E-02,1.91E-02,1.94E-01
8.89E-01,1.70E-01,5.90E-01,5.90E-01,3.40E-01,1.70E-01,6.97E-02,2.86E-02,6.97E-02,2.86E-02,1.89E-02,2.44E+05,9.04E+05,7.45E+05,9.77E+05,1.08E+06,7.40E+05,2.10E+05,1.94E+05,1.64E+04,-3.05E+02,2.14E-04,2.15E-04,1.13E-04,1.20E-02,1.85E-02,1.89E-01
8.89E-01,1.76E-01,5.91E-01,5.91E-01,3.38E-01,1.76E-01,7.21E-02,2.95E-02,7.21E-02,2.95E-02,1.95E-02,2.44E+05,9.32E+05,7.73E+05,1.01E+06,1.10E+06,7.49E+05,2.23E+05,2.08E+05,1.71E+04,-1.44E+03,2.21E-04,2.22E-04,1.15E-04,1.16E-02,1.80E-02,1.85E-01
8.90E-01,1.45E-01,5.85E-01,5.85E-01,3.62E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.60E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,7.16E+05,1.60E+05,1.42E+05,1.39E+04,4.14E+03,1.83E-04,1.86E-04,1.07E-04,1.39E-02,2.10E-02,2.00E-01
8.90E-01,1.51E-01,5.88E-01,5.88E-01,3.51E-01,1.51E-01,6.23E-02,2.56E-02,6.23E-02,2.56E-02,1.66E-02,2.44E+05,8.21E+05,6.61E+05,8.95E+05,9.91E+05,7.16E+05,1.72E+05,1.55E+05,1.42E+04,2.81E+03,1.91E-04,1.93E-04,1.06E-04,1.34E-02,2.05E-02,2.03E-01
8.90E-01,1.57E-01,5.89E-01,5.89E-01,3.49E-01,1.57E-01,6.47E-02,2.66E-02,6.47E-02,2.66E-02,1.73E-02,2.44E+05,8.48E+05,6.89E+05,9.22E+05,1.02E+06,7.26E+05,1.84E+05,1.67E+05,1.49E+04,1.97E+03,1.99E-04,2.00E-04,1.09E-04,1.29E-02,1.98E-02,1.98E-01
8.90E-01,1.64E-01,5.90E-01,5.90E-01,3.47E-01,1.64E-01,6.72E-02,2.76E-02,6.72E-02,2.76E-02,1.80E-02,2.44E+05,8.76E+05,7.17E+05,9.50E+05,1.05E+06,7.36E+05,1.97E+05,1.80E+05,1.56E+04,1.03E+03,2.06E-04,2.08E-04,1.12E-04,1.24E-02,1.92E-02,1.93E-01
8.90E-01,1.70E-01,5.90E-01,5.90E-01,3.45E-01,1.70E-01,6.96E-02,2.85E-02,6.96E-02,2.85E-02,1.87E-02,2.44E+05,9.04E+05,7.45E+05,9.78E+05,1.08E+06,7.45E+05,2.10E+05,1.94E+05,1.63E+04,3.49E+00,2.14E-04,2.15E-04,1.14E-04,1.20E-02,1.86E-02,1.89E-01
8.90E-01,1.76E-01,5.91E-01,5.91E-01,3.42E-01,1.76E-01,7.20E-02,2.95E-02,7.20E-02,2.95E-02,1.94E-02,2.44E+05,9.32E+05,7.73E+05,1.01E+06,1.10E+06,7.54E+05,2.24E+05,2.08E+05,1.70E+04,-1.11E+03,2.22E-04,2.22E-04,1.17E-04,1.16E-02,1.80E-02,1.85E-01
8.98E-01,1.45E-01,5.85E-01,5.85E-01,4.08E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.48E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,7.66E+05,1.63E+05,1.42E+05,1.39E+04,6.67E+03,1.83E-04,1.86E-04,1.22E-04,1.39E-02,2.10E-02,1.88E-01
8.98E-01,1.51E-01,5.92E-01,5.92E-01,3.88E-01,1.51E-01,6.18E-02,2.52E-02,6.18E-02,2.52E-02,1.54E-02,2.44E+05,8.25E+05,6.62E+05,8.99E+05,9.92E+05,7.57E+05,1.75E+05,1.56E+05,1.38E+04,4.81E+03,1.93E-04,1.93E-04,1.17E-04,1.34E-02,2.08E-02,1.98E-01
8.98E-01,1.57E-01,5.92E-01,5.92E-01,3.87E-01,1.57E-01,6.42E-02,2.62E-02,6.42E-02,2.62E-02,1.61E-02,2.44E+05,8.53E+05,6.90E+05,9.26E+05,1.02E+06,7.68E+05,1.87E+05,1.68E+05,1.45E+04,4.11E+03,2.00E-04,2.00E-04,1.20E-04,1.29E-02,2.01E-02,1.92E-01
8.98E-01,1.64E-01,5.93E-01,5.93E-01,3.85E-01,1.64E-01,6.66E-02,2.71E-02,6.66E-02,2.71E-02,1.67E-02,2.44E+05,8.81E+05,7.17E+05,9.54E+05,1.05E+06,7.79E+05,2.00E+05,1.81E+05,1.51E+04,3.33E+03,2.08E-04,2.07E-04,1.23E-04,1.24E-02,1.94E-02,1.87E-01
8.98E-01,1.70E-01,5.94E-01,5.94E-01,3.82E-01,1.70E-01,6.90E-02,2.81E-02,6.90E-02,2.81E-02,1.73E-02,2.44E+05,9.09E+05,7.45E+05,9.82E+05,1.08E+06,7.89E+05,2.13E+05,1.95E+05,1.58E+04,2.47E+03,2.15E-04,2.15E-04,1.26E-04,1.20E-02,1.88E-02,1.83E-01
8.98E-01,1.76E-01,5.94E-01,5.94E-01,3.80E-01,1.76E-01,7.14E-02,2.90E-02,7.14E-02,2.90E-02,1.80E-02,2.44E+05,9.37E+05,7.73E+05,1.01E+06,1.10E+06,7.99E+05,2.27E+05,2.09E+05,1.64E+04,1.51E+03,2.23E-04,2.22E-04,1.29E-04,1.16E-02,1.83E-02,1.79E-01
9.00E-01,1.45E-01,5.85E-01,5.85E-01,4.20E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.45E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,7.79E+05,1.63E+05,1.42E+05,1.39E+04,7.31E+03,1.83E-04,1.86E-04,1.26E-04,1.39E-02,2.10E-02,1.85E-01
9.00E-01,1.51E-01,5.92E-01,5.92E-01,3.98E-01,1.51E-01,6.16E-02,2.51E-02,6.16E-02,2.51E-02,1.51E-02,2.44E+05,8.26E+05,6.62E+05,9.00E+05,9.92E+05,7.67E+05,1.75E+05,1.56E+05,1.37E+04,5.31E+03,1.93E-04,1.93E-04,1.20E-04,1.33E-02,2.09E-02,1.97E-01
9.00E-01,1.57E-01,5.93E-01,5.93E-01,3.96E-01,1.57E-01,6.41E-02,2.61E-02,6.41E-02,2.61E-02,1.57E-02,2.44E+05,8.54E+05,6.90E+05,9.27E+05,1.02E+06,7.79E+05,1.88E+05,1.69E+05,1.43E+04,4.65E+03,2.00E-04,2.00E-04,1.23E-04,1.29E-02,2.02E-02,1.91E-01
9.00E-01,1.64E-01,5.94E-01,5.94E-01,3.94E-01,1.64E-01,6.65E-02,2.70E-02,6.65E-02,2.70E-02,1.64E-02,2.44E+05,8.82E+05,7.18E+05,9.55E+05,1.05E+06,7.90E+05,2.01E+05,1.82E+05,1.50E+04,3.91E+03,2.08E-04,2.07E-04,1.26E-04,1.24E-02,1.95E-02,1.86E-01
9.00E-01,1.70E-01,5.94E-01,5.94E-01,3.92E-01,1.70E-01,6.89E-02,2.79E-02,6.89E-02,2.79E-02,1.70E-02,2.44E+05,9.10E+05,7.45E+05,9.84E+05,1.08E+06,8.01E+05,2.14E+05,1.95E+05,1.56E+04,3.08E+03,2.16E-04,2.14E-04,1.29E-04,1.20E-02,1.89E-02,1.82E-01
9.00E-01,1.76E-01,5.95E-01,5.95E-01,3.90E-01,1.76E-01,7.13E-02,2.89E-02,7.13E-02,2.89E-02,1.76E-02,2.44E+05,9.38E+05,7.73E+05,1.01E+06,1.10E+06,8.11E+05,2.28E+05,2.10E+05,1.63E+04,2.17E+03,2.23E-04,2.22E-04,1.32E-04,1.16E-02,1.84E-02,1.78E-01
9.07E-01,1.45E-01,5.85E-01,5.85E-01,4.61E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.35E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,8.25E+05,1.66E+05,1.42E+05,1.39E+04,9.61E+03,1.83E-04,1.86E-04,1.39E-04,1.39E-02,2.10E-02,1.77E-01
9.07E-01,1.51E-01,5.93E-01,5.93E-01,4.37E-01,1.51E-01,6.15E-02,2.50E-02,6.15E-02,2.50E-02,1.41E-02,2.44E+05,8.27E+05,6.62E+05,9.01E+05,9.92E+05,8.12E+05,1.78E+05,1.57E+05,1.36E+04,7.51E+03,1.93E-04,1.93E-04,1.32E-04,1.33E-02,2.10E-02,1.89E-01
9.07E-01,1.57E-01,5.96E-01,5.96E-01,4.30E-01,1.57E-01,6.36E-02,2.57E-02,6.36E-02,2.57E-02,1.46E-02,2.44E+05,8.58E+05,6.90E+05,9.31E+05,1.02E+06,8.17E+05,1.90E+05,1.70E+05,1.39E+04,6.53E+03,2.01E-04,2.00E-04,1.33E-04,1.28E-02,2.05E-02,1.89E-01
9.07E-01,1.64E-01,5.97E-01,5.97E-01,4.28E-01,1.64E-01,6.60E-02,2.66E-02,6.60E-02,2.66E-02,1.52E-02,2.44E+05,8.86E+05,7.18E+05,9.59E+05,1.05E+06,8.29E+05,2.03E+05,1.83E+05,1.45E+04,5.93E+03,2.09E-04,2.07E-04,1.36E-04,1.24E-02,1.98E-02,1.84E-01
9.07E-01,1.70E-01,5.98E-01,5.98E-01,4.26E-01,1.70E-01,6.84E-02,2.75E-02,6.84E-02,2.75E-02,1.58E-02,2.44E+05,9.14E+05,7.46E+05,9.88E+05,1.08E+06,8.41E+05,2.17E+05,1.97E+05,1.52E+04,5.24E+03,2.17E-04,2.14E-04,1.39E-04,1.20E-02,1.92E-02,1.80E-01
9.07E-01,1.76E-01,5.98E-01,5.98E-01,4.24E-01,1.76E-01,7.08E-02,2.84E-02,7.08E-02,2.84E-02,1.64E-02,2.44E+05,9.42E+05,7.73E+05,1.02E+06,1.10E+06,8.52E+05,2.31E+05,2.11E+05,1.58E+04,4.47E+03,2.24E-04,2.21E-04,1.42E-04,1.16E-02,1.86E-02,1.76E-01
9.10E-01,1.45E-01,5.85E-01,5.85E-01,4.78E-01,1.45E-01,6.02E-02,2.50E-02,6.02E-02,2.50E-02,1.31E-02,2.44E+05,7.90E+05,6.33E+05,8.64E+05,9.63E+05,8.45E+05,1.67E+05,1.42E+05,1.39E+04,1.06E+04,1.83E-04,1.86E-04,1.45E-04,1.39E-02,2.10E-02,1.75E-01
9.10E-01,1.51E-01,5.93E-01,5.93E-01,4.56E-01,1.51E-01,6.15E-02,2.50E-02,6.15E-02,2.50E-02,1.36E-02,2.44E+05,8.27E+05,6.62E+05,9.01E+05,9.92E+05,8.33E+05,1.79E+05,1.57E+05,1.36E+04,8.55E+03,1.93E-04,1.93E-04,1.38E-04,1.33E-02,2.10E-02,1.86E-01
9.10E-01,1.57E-01,5.98E-01,5.98E-01,4.44E-01,1.57E-01,6.34E-02,2.55E-02,6.34E-02,2.55E-02,1.42E-02,2.44E+05,8.59E+05,6.90E+05,9.33E+05,1.02E+06,8.34E+05,1.91E+05,1.70E+05,1.38E+04,7.35E+03,2.02E-04,2.00E-04,1.37E-04,1.28E-02,2.06E-02,1.89E-01
9.10E-01,1.64E-01,5.98E-01,5.98E-01,4.43E-01,1.64E-01,6.58E-02,2.64E-02,6.58E-02,2.64E-02,1.47E-02,2.44E+05,8.87E+05,7.18E+05,9.61E+05,1.05E+06,8.46E+05,2.05E+05,1.83E+05,1.44E+04,6.80E+03,2.09E-04,2.07E-04,1.41E-04,1.24E-02,1.99E-02,1.84E-01
9.10E-01,1.70E-01,5.99E-01,5.99E-01,4.41E-01,1.70E-01,6.81E-02,2.73E-02,6.81E-02,2.73E-02,1.53E-02,2.44E+05,9.16E+05,7.46E+05,9.89E+05,1.08E+06,8.58E+05,2.18E+05,1.97E+05,1.50E+04,6.17E+03,2.17E-04,2.14E-04,1.44E-04,1.20E-02,1.93E-02,1.79E-01
9.10E-01,1.76E-01,6.00E-01,6.00E-01,4.39E-01,1.76E-01,7.05E-02,2.82E-02,7.05E-02,2.82E-02,1.59E-02,2.44E+05,9.44E+05,7.73E+05,1.02E+06,1.10E+06,8.70E+05,2.33E+05,2.11E+05,1.56E+04,5.46E+03,2.25E-04,2.21E-04,1.47E-04,1.16E-02,1.87E-02,1.75E-01
9.16E-01,1.45E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.16E-01,1.51E-01,5.93E-01,5.93E-01,4.92E-01,1.51E-01,6.15E-02,2.50E-02,6.15E-02,2.50E-02,1.27E-02,2.44E+05,8.27E+05,6.62E+05,9.01E+05,9.92E+05,8.76E+05,1.81E+05,1.57E+05,1.36E+04,1.07E+04,1.93E-04,1.93E-04,1.51E-04,1.33E-02,2.10E-02,1.81E-01
9.16E-01,1.57E-01,6.00E-01,6.00E-01,4.74E-01,1.57E-01,6.29E-02,2.51E-02,6.29E-02,2.51E-02,1.32E-02,2.44E+05,8.63E+05,6.90E+05,9.36E+05,1.02E+06,8.69E+05,1.94E+05,1.71E+05,1.34E+04,8.99E+03,2.03E-04,2.00E-04,1.46E-04,1.28E-02,2.09E-02,1.89E-01
9.16E-01,1.64E-01,6.01E-01,6.01E-01,4.72E-01,1.64E-01,6.53E-02,2.61E-02,6.53E-02,2.61E-02,1.37E-02,2.44E+05,8.91E+05,7.18E+05,9.65E+05,1.05E+06,8.82E+05,2.07E+05,1.84E+05,1.40E+04,8.56E+03,2.10E-04,2.07E-04,1.50E-04,1.24E-02,2.02E-02,1.84E-01
9.16E-01,1.70E-01,6.02E-01,6.02E-01,4.71E-01,1.70E-01,6.77E-02,2.70E-02,6.77E-02,2.70E-02,1.43E-02,2.44E+05,9.19E+05,7.46E+05,9.93E+05,1.08E+06,8.95E+05,2.21E+05,1.98E+05,1.45E+04,8.05E+03,2.18E-04,2.14E-04,1.53E-04,1.20E-02,1.96E-02,1.79E-01
9.16E-01,1.76E-01,6.02E-01,6.02E-01,4.69E-01,1.76E-01,7.00E-02,2.78E-02,7.00E-02,2.78E-02,1.48E-02,2.44E+05,9.48E+05,7.73E+05,1.02E+06,1.10E+06,9.07E+05,2.35E+05,2.13E+05,1.51E+04,7.46E+03,2.26E-04,2.21E-04,1.57E-04,1.16E-02,1.90E-02,1.75E-01
9.20E-01,1.45E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.20E-01,1.51E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.20E-01,1.57E-01,6.02E-01,6.02E-01,4.96E-01,1.57E-01,6.27E-02,2.50E-02,6.27E-02,2.50E-02,1.26E-02,2.44E+05,8.64E+05,6.90E+05,9.38E+05,1.02E+06,8.96E+05,1.95E+05,1.72E+05,1.32E+04,1.03E+04,2.03E-04,2.00E-04,1.53E-04,1.28E-02,2.10E-02,1.88E-01
9.20E-01,1.64E-01,6.03E-01,6.03E-01,4.93E-01,1.64E-01,6.50E-02,2.58E-02,6.50E-02,2.58E-02,1.31E-02,2.44E+05,8.94E+05,7.18E+05,9.67E+05,1.05E+06,9.07E+05,2.09E+05,1.85E+05,1.37E+04,9.75E+03,2.11E-04,2.07E-04,1.56E-04,1.24E-02,2.04E-02,1.84E-01
9.20E-01,1.70E-01,6.04E-01,6.04E-01,4.91E-01,1.70E-01,6.73E-02,2.67E-02,6.73E-02,2.67E-02,1.36E-02,2.44E+05,9.22E+05,7.46E+05,9.96E+05,1.08E+06,9.20E+05,2.23E+05,1.99E+05,1.42E+04,9.32E+03,2.19E-04,2.14E-04,1.60E-04,1.19E-02,1.97E-02,1.79E-01
9.20E-01,1.76E-01,6.04E-01,6.04E-01,4.89E-01,1.76E-01,6.97E-02,2.76E-02,6.97E-02,2.76E-02,1.41E-02,2.44E+05,9.51E+05,7.73E+05,1.02E+06,1.10E+06,9.33E+05,2.37E+05,2.13E+05,1.48E+04,8.82E+03,2.27E-04,2.21E-04,1.63E-04,1.16E-02,1.92E-02,1.75E-01
9.25E-01,1.45E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.25E-01,1.51E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.25E-01,1.57E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.25E-01,1.64E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.25E-01,1.70E-01,6.06E-01,6.06E-01,5.17E-01,1.70E-01,6.69E-02,2.64E-02,6.69E-02,2.64E-02,1.27E-02,2.44E+05,9.25E+05,7.46E+05,9.99E+05,1.08E+06,9.53E+05,2.25E+05,2.00E+05,1.39E+04,1.09E+04,2.20E-04,2.14E-04,1.68E-04,1.19E-02,2.00E-02,1.81E-01
9.25E-01,1.76E-01,6.07E-01,6.07E-01,5.15E-01,1.76E-01,6.92E-02,2.72E-02,6.92E-02,2.72E-02,1.32E-02,2.44E+05,9.54E+05,7.73E+05,1.03E+06,1.10E+06,9.67E+05,2.39E+05,2.14E+05,1.44E+04,1.05E+04,2.28E-04,2.21E-04,1.72E-04,1.15E-02,1.94E-02,1.77E-01
9.30E-01,1.45E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.30E-01,1.51E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.30E-01,1.57E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.30E-01,1.64E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.30E-01,1.70E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.30E-01,1.76E-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
,,,,,,,,,,,,,,,,,,,,,,,,,,
,,,,,,,,,,,,,,,,,,,,,,,,,,
,,,,,,,,,,,,,,,,,,,,,,,,,,
Binary file added solution.pickle
View file
Open in desktop
Binary file not shown.
42 changes: 33 additions & 9 deletions src/wrd/components/ro_stage.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,6 +3,7 @@
assert_optimal_termination,
units as pyunits,
TransformationFactory,
value,
)
from pyomo.network import Arc

Expand All @@ -12,6 +13,7 @@
from idaes.core import FlowsheetBlock
from idaes.models.unit_models import StateJunction
from idaes.core.util.scaling import calculate_scaling_factors
from idaes.core.util.exceptions import InitializationError

from watertap.costing import WaterTAPCosting
from watertap.core.util.model_diagnostics.infeasible import *
Expand Down Expand Up @@ -131,9 +133,9 @@ def build_ro_stage(
build_ro(blk.ro, stage_num=stage_num, file=file, prop_package=prop_package)

blk.product = StateJunction(property_package=prop_package)
touch_flow_and_conc(blk.feed)
touch_flow_and_conc(blk.product)
blk.disposal = StateJunction(property_package=prop_package)
touch_flow_and_conc(blk.feed)
touch_flow_and_conc(blk.disposal)

# Arcs to connect the unit models
blk.feed_to_pump = Arc(source=blk.feed.outlet, destination=blk.pump.feed.inlet)
Expand All @@ -159,10 +161,10 @@ def initialize_system(m):


def initialize_ro_stage(blk):

blk.feed.initialize()
propagate_state(blk.feed_to_pump)

touch_flow_and_conc(blk.pump.unit)
initialize_pump(blk.pump)

propagate_state(blk.pump_to_ro)
Expand All @@ -171,7 +173,29 @@ def initialize_ro_stage(blk):
propagate_state(blk.ro_to_product)
blk.product.initialize()
propagate_state(blk.ro_to_disposal)
blk.disposal.initialize()
try:
blk.disposal.initialize()
except InitializationError:
ro_disposal_props = blk.ro.disposal.properties[0]
try:
blk.disposal.initialize(
state_args={
"flow_mass_phase_comp": {
("Liq", "H2O"): value(
ro_disposal_props.flow_mass_phase_comp["Liq", "H2O"]
),
("Liq", "NaCl"): value(
ro_disposal_props.flow_mass_phase_comp["Liq", "NaCl"]
),
},
"temperature": value(ro_disposal_props.temperature),
"pressure": value(ro_disposal_props.pressure),
}
)
except InitializationError:
# This state block is a pass-through junction; if initialization fails,
# continue with propagated state so the full train can still initialize.
pass


def report_ro_stage(blk, w=30, add_costing=True):
Expand Down Expand Up @@ -294,12 +318,12 @@ def run_march_stages():
# run_august_stages()
# run_march_stages()
m = main(
Qin=2800, # gpm?
Cin=0.528,
Qin=396 * 0.29 / 0.25, # gpm?
Cin=2.58,
Tin=295,
Pin=35.4 * pyunits.psi,
stage_num=1,
file="wrd_inputs_2800_gpm.yaml",
Pin=100 * pyunits.psi,
stage_num=3,
file="wrd_inputs_8_19_21.yaml",
)

m.fs.ro_stage.ro.unit.recovery_vol_phase[0, "Liq"].fix(0.5)
Expand Down
2 changes: 2 additions & 0 deletions src/wrd/surrogate_models/S1_S2_RR_not_equal.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,2 @@
# Recovery,Feed Flow,Stage1 RR,Stage2 RR,Stage3 RR,Stage1 Qin (m3/s),Stage2 Qin (m3/s),Stage3 Qin (m3/s),Stage1 Brine (m3/s),Stage2 Brine (m3/s),Stage3 Brine (m3/s),Stage1 Pin (Pa),Stage2 Pin (Pa),Stage3 Pin (Pa),Stage1 Pout (Pa),Stage2 Pout (Pa),Stage3 Pout (Pa),Total Power (W),Stage1 Power (W),Stage2 Power (W),Stage3 Power (W),Stage1 First Elem. Prod. (m3/s),Stage2 First Elem. Prod. (m3/s),Stage3 First Elem. Prod. (m3/s),Stage1 Perm. Conc. (kg/m3),Stage2 Perm. Conc. (kg/m3),Stage3 Perm. Conc. (kg/m3),Stage1 Pump Speed (-),Stage2 Pump Speed (-)
8.80E-01,1.45E-01,5.58E-01,5.96E-01,3.28E-01,1.45E-01,6.41E-02,2.59E-02,6.41E-02,2.59E-02,1.74E-02,2.44E+05,7.59E+05,6.80E+05,8.33E+05,1.01E+06,6.80E+05,1.56E+05,1.34E+05,2.15E+04,1.25E-08,1.75E-04,2.00E-04,9.95E-05,1.41E-02,1.85E-02,1.99E-01,8.50E-01,1.00E+00
5 changes: 5 additions & 0 deletions src/wrd/surrogate_models/S1_S2_TEST_POINT.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
# Recovery,Feed Flow,Stage1 RR,Stage2 RR,Stage3 RR,Stage1 Qin (m3/s),Stage2 Qin (m3/s),Stage3 Qin (m3/s),Stage1 Brine (m3/s),Stage2 Brine (m3/s),Stage3 Brine (m3/s),Stage1 Pin (Pa),Stage2 Pin (Pa),Stage3 Pin (Pa),Stage1 Pout (Pa),Stage2 Pout (Pa),Stage3 Pout (Pa),Total Power (W),Stage1 Power (W),Stage2 Power (W),Stage3 Power (W),Stage1 First Elem. Prod. (m3/s),Stage2 First Elem. Prod. (m3/s),Stage3 First Elem. Prod. (m3/s),Stage1 Perm. Conc. (kg/m3),Stage2 Perm. Conc. (kg/m3),Stage3 Perm. Conc. (kg/m3),Stage1 Pump Speed (-),Stage2 Pump Speed (-)
8.800000e-01,1.450000e-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
8.800000e-01,1.761111e-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.250000e-01,1.450000e-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
9.250000e-01,1.761111e-01,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan
Loading
Loading