Commit 7922b4c2 authored by hazrmard's avatar hazrmard
Browse files

created cooling tower Env, made preprocessing a sub-package, re-factored...

created cooling tower Env, made preprocessing a sub-package, re-factored utils, new plotting function, ICAART notebook
parent 81d5d209
Loading
Loading
Loading
Loading
+1 −0
Original line number Diff line number Diff line
@@ -6,6 +6,7 @@ questions.txt
desktop.ini

.vscode/
src/bin/*
*.pyc
__pycache__
*.ipynb_checkpoints
 No newline at end of file

TODO_libraries.md

0 → 100644
+5 −0
Original line number Diff line number Diff line
* [Platypus - Multi-objective optimization][1]
* [Distributed evolutionary algorithms in Pyhon - ][2]

[1]: https://platypus.readthedocs.io/en/latest/
[2]: https://deap.readthedocs.io/en/master/
 No newline at end of file
+1 −1
Original line number Diff line number Diff line
@@ -52,7 +52,7 @@ Each cooling tower/chiller system has the following parameters:

15. `PerHumidity`: Relative ambient humidity.

16. `TempWetbulb`: Wet-bulb temperature.
16. `TempWetBulb`: Wet-bulb temperature.

17. `PerChiLoad`: Cooling load of the chiller as a fraction of maximum electrical capacity. The maximum cooling capacity in tons is 800 tons. The ratio `Tons / 800` should give roughly the same value as `PerChilLoad`.

+6 −0
Original line number Diff line number Diff line
@@ -18,9 +18,15 @@ Normalized [mutual information][3] was calculated for each pair of variables. A

![Mutual information (MI)](img/7-mutual-information.png)

###


## Causality

### Vector auto-regression

TODO

### Granger causality test

TODO
+4 −4
Original line number Diff line number Diff line
@@ -51,7 +51,7 @@ The total cooling depends on:
A first-order approximation would be:

$$
\frac{d}{dt} E_{evap}(T) \propto \frac{T(t) T_a v_{air} R}{T_w}
\frac{d}{dt} E_{evap}(T) \propto T(t) (T_a - T_w) v_{air} R
$$

Where:
@@ -67,9 +67,9 @@ Combining, and adding constant of proportionality $k$:

$$
\begin{align*}
-c_m m \frac{d}{dt} T           &= k \frac{T(t) T_a v_{air} R}{T_w} \\
\frac{1}{T(t)} \frac{d}{dt} T   &= -\frac{k T_a v_{air} R}{T_w c_m m} \\
T(t) &= T(0) e^{-\frac{k T_a v_{air} R}{T_w c_m m} t}
-c_m m \frac{d}{dt} T           &= k T(t) (T_a - T_w) v_{air} R \\
\frac{1}{T(t)} \frac{d}{dt} T   &= -\frac{k (T_a - T_w) v_{air} R}{c_m m} \\
T(t) &= T(0) e^{-\frac{k (T_a - T_w) v_{air} R}{c_m m} t}
\end{align*}
$$

Loading