# HG changeset patch
# User Michiel Broek <mbroek@mbse.eu>
# Date 1650894741 -7200
# Node ID 19156b8b339f12d04eeee09517951369750ddddd
# Parent  b36d249512ccdfebd847a322a3e727bd13b90c62
Added calcSparge.

diff -r b36d249512cc -r 19156b8b339f src/EditRecipe.cpp
--- a/src/EditRecipe.cpp	Mon Apr 25 12:53:36 2022 +0200
+++ b/src/EditRecipe.cpp	Mon Apr 25 15:52:21 2022 +0200
@@ -712,6 +712,9 @@
     connect(ui->w1_nameEdit, QOverload<int>::of(&QComboBox::currentIndexChanged), this, &EditRecipe::w1_name_changed);
     connect(ui->w2_nameEdit, QOverload<int>::of(&QComboBox::currentIndexChanged), this, &EditRecipe::w2_name_changed);
     connect(ui->w2_volEdit, QOverload<double>::of(&QDoubleSpinBox::valueChanged), this, &EditRecipe::w2_volume_changed);
+    connect(ui->sp_sourceEdit, QOverload<int>::of(&QComboBox::currentIndexChanged), this, &EditRecipe::sp_source_changed);
+    connect(ui->sp_acidtypeEdit, QOverload<int>::of(&QComboBox::currentIndexChanged), this, &EditRecipe::sp_type_changed);
+    connect(ui->sp_phEdit, QOverload<double>::of(&QDoubleSpinBox::valueChanged), this, &EditRecipe::sp_ph_changed);
 
     ui->saveButton->setEnabled(false);
     ui->deleteButton->setEnabled((id >= 0) ? true:false);
diff -r b36d249512cc -r 19156b8b339f src/EditRecipe.h
--- a/src/EditRecipe.h	Mon Apr 25 12:53:36 2022 +0200
+++ b/src/EditRecipe.h	Mon Apr 25 15:52:21 2022 +0200
@@ -107,6 +107,9 @@
     void mw_ph_changed(double val);
     void mw_acid_changed(double val);
     void mw_type_changed(int val);
+    void sp_source_changed(int val);
+    void sp_type_changed(int val);
+    void sp_ph_changed(double val);
     void step_name_changed(QString val);
     void step_type_changed(int val);
     void step_temp_changed(double val);
@@ -168,6 +171,7 @@
     double ProtonDeficit(double pHZ);
     double MashpH();
     void calcWater();
+    void calcSparge();
     double GetBUGU();
     double GetOptSO4Clratio();
     void calcYeast();
diff -r b36d249512cc -r 19156b8b339f src/EditRecipeTab7.cpp
--- a/src/EditRecipeTab7.cpp	Mon Apr 25 12:53:36 2022 +0200
+++ b/src/EditRecipeTab7.cpp	Mon Apr 25 15:52:21 2022 +0200
@@ -396,7 +396,120 @@
     ui->wb_hco3Edit->setStyleSheet((bicarbonate > 250) ? "background-color: red":"background-color: green");
     ui->wb_caco3Edit->setStyleSheet((bicarbonate > 250) ? "background-color: red":"background-color: green");
 
-    // calcSparge();
+    calcSparge();
+}
+
+
+/*
+ * Based on the work of ajDeLange.
+ */
+void EditRecipe::calcSparge()
+{
+    double TargetpH = recipe->sparge_ph;
+    double Source_pH = recipe->w1_ph;
+    double Source_alkalinity = recipe->w1_total_alkalinity;
+
+    qDebug() << "calcSparge()";
+
+    const QSignalBlocker blocker1(ui->sp_sourceEdit);
+
+    // Select watersource or fallback to the first source.
+    if (recipe->sparge_source == 1) { // Source 2
+	if (recipe->w2_ph > 0.0 && recipe->w2_amount > 0) {
+	    Source_pH = recipe->w2_ph;
+	    Source_alkalinity = recipe->w2_total_alkalinity;
+	} else {
+	    recipe->sparge_source = 0; // Source 1
+	    ui->sp_sourceEdit->setCurrentIndex(0);
+	}
+    } else if (recipe->sparge_source == 2) { // Mixed
+	if (recipe->w2_ph > 0.0 && recipe->w2_amount > 0) {
+	    Source_pH = recipe->wg_ph;
+	    Source_alkalinity = recipe->wg_total_alkalinity;
+	} else {
+	    recipe->sparge_source = 0; // Source 1
+            ui->sp_sourceEdit->setCurrentIndex(0);
+	}
+    }
+
+    // Step 1: Compute the mole fractions of carbonic (f1o), bicarbonate (f2o) and carbonate(f3o) at the water pH
+    double r1 = pow(10, Source_pH - 6.35);
+    double r2 = pow(10, Source_pH - 10.33);
+    double d = 1 + r1 + r1 * r2;
+    double f1 = 1 / d;
+    double f3 = r1 * r2 / d;
+
+    // Step 2. Compute the mole fractions at pH = 4.3 (the pH which defines alkalinity)
+    double r143 = pow(10, 4.3 - 6.35);
+    double r243 = pow(10, 4.3 - 10.33);
+    double d43 = 1 + r143 + r143 * r243;
+    double f143 = 1 / d43;
+    double f343 = r143 * r243 / d43;
+
+    // Step 4. Solve
+    //double Ct = (Source_alkalinity - 1000 * (pow(10, -4.3) - pow(10, -Source_pH))) / ((f143 - f1) + (f3 - f343));
+    double Ct = Source_alkalinity / 50 / ((f143 - f1) + (f3 - f343));
+
+    // Step 5. Compute mole fractions at desired pH
+    double r1g = pow(10, TargetpH - 6.35);
+    double r2g = pow(10, TargetpH - 10.33);
+    double dg = 1 + r1g + r1g * r2g;
+    double f1g = 1 / dg;
+    double f3g = r1g * r2g / dg;
+
+    // Step 6. Use these to compute the milliequivalents acid required per liter (mEq/L)
+    double Acid = Ct * ((f1g - f1) + (f3 - f3g)) + pow(10, -TargetpH) - pow(10, -Source_pH);  //mEq/l
+    Acid += 0.01;   // Add acid that would be required for distilled water.
+
+    //Step 8. Get the acid data.
+    int AT = recipe->sparge_acid_type;
+    if (AT < 0 || AT >= my_acids.size()) {
+	AT = 0;
+	recipe->sparge_acid_type = 0;
+	ui->sp_acidtypeEdit->setCurrentIndex(0);
+	recipe->sparge_acid_perc = my_acids[0].AcidPrc;
+	ui->sp_acidpercEdit->setValue(recipe->sparge_acid_perc);
+    }
+    double fract = Utils::CalcFrac(TargetpH, my_acids[AT].pK1, my_acids[AT].pK2, my_acids[AT].pK3);
+
+    // Step 9. Now divide the mEq required by the "fraction". This is the required number of moles of acid.
+    Acid /= fract;
+
+    // Step 10. Multiply by molecular weight of the acid
+    Acid *= my_acids[AT].MolWt; //mg
+
+    // Step 11. Divide by Specific Gravity and Percentage to get the final ml.
+    Acid = Acid / my_acids[AT].AcidSG / (recipe->sparge_acid_perc / 100); //ml
+    Acid *= recipe->sparge_volume; //ml acid total
+    Acid = round(Acid * 100.0) / 100.0;
+    recipe->sparge_acid_amount = Acid / 1000;
+    ui->sp_acidvolEdit->setValue(Acid);
+}
+
+
+void EditRecipe::sp_source_changed(int val)
+{
+    recipe->sparge_source = val;
+    calcSparge();
+    is_changed();
+}
+
+
+void EditRecipe::sp_type_changed(int val)
+{
+    recipe->sparge_acid_type = val;
+    recipe->sparge_acid_perc = my_acids[val].AcidPrc;
+    ui->sp_acidpercEdit->setValue(recipe->sparge_acid_perc);
+    calcSparge();
+    is_changed();
+}
+
+
+void EditRecipe::sp_ph_changed(double val)
+{
+    recipe->sparge_ph = val;
+    calcSparge();
+    is_changed();
 }
 
 
diff -r b36d249512cc -r 19156b8b339f ui/EditRecipe.ui
--- a/ui/EditRecipe.ui	Mon Apr 25 12:53:36 2022 +0200
+++ b/ui/EditRecipe.ui	Mon Apr 25 15:52:21 2022 +0200
@@ -4253,6 +4253,9 @@
          <property name="accelerated">
           <bool>true</bool>
          </property>
+         <property name="suffix">
+          <string> ml</string>
+         </property>
          <property name="decimals">
           <number>2</number>
          </property>
@@ -4485,8 +4488,14 @@
          <property name="alignment">
           <set>Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter</set>
          </property>
+         <property name="readOnly">
+          <bool>true</bool>
+         </property>
+         <property name="buttonSymbols">
+          <enum>QAbstractSpinBox::NoButtons</enum>
+         </property>
          <property name="accelerated">
-          <bool>true</bool>
+          <bool>false</bool>
          </property>
          <property name="suffix">
           <string>%</string>
@@ -4513,8 +4522,17 @@
          <property name="alignment">
           <set>Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter</set>
          </property>
+         <property name="readOnly">
+          <bool>true</bool>
+         </property>
+         <property name="buttonSymbols">
+          <enum>QAbstractSpinBox::NoButtons</enum>
+         </property>
          <property name="accelerated">
-          <bool>true</bool>
+          <bool>false</bool>
+         </property>
+         <property name="suffix">
+          <string> ml</string>
          </property>
          <property name="decimals">
           <number>2</number>