1 files changed, 31 insertions, 0 deletions
diff --git a/src/main/Simulator/Simulator/UnitOperations/PFR/Integral.mo b/src/main/Simulator/Simulator/UnitOperations/PFR/Integral.mo
new file mode 100644
index 0000000..bd03e89
--- /dev/null
+++ b/src/main/Simulator/Simulator/UnitOperations/PFR/Integral.mo
@@ -0,0 +1,31 @@
+within Simulator.UnitOperations.PFR;
+
+  function Integral
+    extends Modelica.Math.Nonlinear.Interfaces.partialScalarFunction;
+
+        input Integer Nc;
+        input Integer Nr;
+        input Integer Base_comp;
+        input Real Co_dummy[Nc - 1];
+        input Real DO_dummy[Nc - 1, Nr];
+        input Real X_dummy[Nc - 1];
+        input Real X;
+        input Integer Order;
+        input Real DO[Nc, Nr];
+        input Real Co[Nc];
+        input Real Sc[Nc, Nr];
+        input Real Bc[Nr];
+        input Real Fao;
+        input Real k;
+        Real Rate;
+      algorithm
+        Rate := 1;
+        for i in 2:Nc loop
+          if DO[Base_comp, 1] == 0 then
+            Rate := Rate * product((Co[i] + Sc[i, 1] / Bc[1] * Co[Base_comp] * u) ^ DO[i, 1]);
+          else
+            Rate := Rate * product((Co_dummy[i - 1] * (1 - X_dummy[i - 1])) ^ DO_dummy[i - 1, 1]);
+          end if;
+        end for;
+        y := Fao / (k * (Co[Base_comp] * (1 - u)) ^ DO[Base_comp, 1] * Rate);
+      end Integral;