@@ -484,9 +484,9 @@ int CDUMUXData::WriteInputForDuMux(CRFProcess* m_pcs,
484
484
// TF commented out since we want to use the improved PointDuMux
485
485
// class
486
486
// value =
487
- // pcs_vector[this->ProcessIndex_CO2inLiquid]->GetNodeValue(i,
488
- // indexConcentration_DIC) * (this->Molweight_CO2 / 1000) /
489
- // this->NodeData[i]->phase_density[0];
487
+ // pcs_vector[this->ProcessIndex_CO2inLiquid]->GetNodeValue(i,
488
+ // indexConcentration_DIC) * (this->Molweight_CO2 / 1000) /
489
+ // this->NodeData[i]->phase_density[0];
490
490
value = pcs_vector[this ->ProcessIndex_CO2inLiquid ]->GetNodeValue (
491
491
i, indexConcentration_DIC) *
492
492
(this ->Molweight_CO2 / 1000 ) /
@@ -583,10 +583,11 @@ int CDUMUXData::WriteInputForDuMux(CRFProcess* m_pcs,
583
583
tempstring = temp.str ();
584
584
// TF commented out since we want to use the improved PointDuMux class
585
585
// temp.str(""); temp.clear(); temp << this->NodeData[i]->x;
586
- // tempstring += "; " + temp.str(); temp.str(""); temp.clear(); temp <<
587
- // this->NodeData[i]->y; tempstring += "; " + temp.str(); temp.str("");
588
- // temp.clear(); temp << this->NodeData[i]->z; tempstring += "; " +
589
- // temp.str();
586
+ // tempstring += "; " + temp.str(); temp.str(""); temp.clear();
587
+ // temp
588
+ // << this->NodeData[i]->y; tempstring += "; " + temp.str();
589
+ // temp.str(""); temp.clear(); temp << this->NodeData[i]->z; tempstring
590
+ // += "; " + temp.str();
590
591
for (size_t k (0 ); k < 3 ; k++)
591
592
{
592
593
temp.str (" " );
@@ -597,29 +598,29 @@ int CDUMUXData::WriteInputForDuMux(CRFProcess* m_pcs,
597
598
598
599
// TF commented out since we want to use the improved PointDuMux class
599
600
// temp.str(""); temp.clear(); temp.precision(12); temp <<
600
- // this->NodeData[i]->phase_pressure[0]; tempstring
601
+ // this->NodeData[i]->phase_pressure[0]; tempstring
601
602
// += "; " + temp.str();
602
603
// temp.str(""); temp.clear(); temp.precision(12); temp <<
603
- // this->NodeData[i]->phase_pressure[1]; tempstring
604
+ // this->NodeData[i]->phase_pressure[1]; tempstring
604
605
// += "; " + temp.str();
605
606
// temp.str(""); temp.clear(); temp.precision(12); temp <<
606
- // this->NodeData[i]->phase_saturation[0];
607
+ // this->NodeData[i]->phase_saturation[0];
607
608
// tempstring
608
609
// += "; " + temp.str();
609
610
// temp.str(""); temp.clear(); temp.precision(12); temp <<
610
- // this->NodeData[i]->phase_saturation[1];
611
+ // this->NodeData[i]->phase_saturation[1];
611
612
// tempstring
612
613
// += "; " + temp.str();
613
614
// temp.str(""); temp.clear(); temp.precision(12); temp <<
614
- // this->NodeData[i]->phase_density[0]; tempstring
615
+ // this->NodeData[i]->phase_density[0]; tempstring
615
616
// +=
616
617
// "; " + temp.str();
617
618
// temp.str(""); temp.clear(); temp.precision(12); temp <<
618
- // this->NodeData[i]->phase_density[1]; tempstring
619
+ // this->NodeData[i]->phase_density[1]; tempstring
619
620
// +=
620
621
// "; " + temp.str();
621
622
// temp.str(""); temp.clear(); temp.precision(12); temp <<
622
- // this->NodeData[i]->CO2inLiquid; tempstring += ";
623
+ // this->NodeData[i]->CO2inLiquid; tempstring += ";
623
624
// "
624
625
// + temp.str();
625
626
temp.str (" " );
@@ -663,10 +664,10 @@ int CDUMUXData::WriteInputForDuMux(CRFProcess* m_pcs,
663
664
i, indexConcentration_DIC);
664
665
// TF commented out since we want to use the improved PointDuMux class
665
666
// mass_CO2_gas = node_volume *
666
- // this->NodeData[i]->phase_saturation[1] *
667
+ // this->NodeData[i]->phase_saturation[1] *
667
668
// this->NodeData[i]->phase_density[1];
668
669
// mass_CO2_water = node_volume *
669
- // this->NodeData[i]->phase_saturation[0] * concentration_CO2_water *
670
+ // this->NodeData[i]->phase_saturation[0] * concentration_CO2_water *
670
671
// this->Molweight_CO2 * 0.001;
671
672
mass_CO2_gas = node_volume *
672
673
this ->NodeData [i]->getPhaseSaturation ()[1 ] *
@@ -834,11 +835,11 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
834
835
// TF commented out since we want to use the improved PointDuMux
835
836
// class
836
837
// if ((atof(TextFile->Data_separated[i][1].data()) !=
837
- // this->NodeData[i]->x) ||
838
+ // this->NodeData[i]->x) ||
838
839
// (atof(TextFile->Data_separated[i][2].data()) !=
839
- // this->NodeData[i]->y) ||
840
+ // this->NodeData[i]->y) ||
840
841
// (atof(TextFile->Data_separated[i][3].data()) !=
841
- // this->NodeData[i]->z)) {
842
+ // this->NodeData[i]->z)) {
842
843
if ((atof (TextFile->Data_separated [i][1 ].data ()) !=
843
844
(*(this ->NodeData [i]))[0 ]) ||
844
845
(atof (TextFile->Data_separated [i][2 ].data ()) !=
@@ -865,7 +866,7 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
865
866
// TF commented out since we want to use the improved PointDuMux
866
867
// class
867
868
// this->NodeData[i]->phase_pressure[j] =
868
- // atof(TextFile->Data_separated[i][column].data()) *
869
+ // atof(TextFile->Data_separated[i][column].data()) *
869
870
// Multiplier;
870
871
this ->NodeData [i]->getPhasePressure ()[j] =
871
872
atof (TextFile->Data_separated [i][column].data ()) * Multiplier;
@@ -891,7 +892,7 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
891
892
// TF commented out since we want to use the improved PointDuMux
892
893
// class
893
894
// this->NodeData[i]->q[j][k] =
894
- // atof(TextFile->Data_separated[i][column].data()) *
895
+ // atof(TextFile->Data_separated[i][column].data()) *
895
896
// Multiplier;
896
897
this ->NodeData [i]->getQ ()[j][k] =
897
898
atof (TextFile->Data_separated [i][column].data ()) *
@@ -900,32 +901,32 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
900
901
column += 1 ;
901
902
// TF commented out since we want to use the improved PointDuMux class
902
903
// this->NodeData[i]->CO2inLiquid =
903
- // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
904
+ // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
904
905
this ->NodeData [i]->setCO2InLiquid (
905
906
atof (TextFile->Data_separated [i][column].data ()) * Multiplier);
906
907
// cout << i << " " << this->NodeData[i]->CO2inLiquid << "\n";
907
908
column += 1 ;
908
909
// TF commented out since we want to use the improved PointDuMux class
909
910
// this->NodeData[i]->NaClinLiquid =
910
- // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
911
+ // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
911
912
this ->NodeData [i]->setNaClInLiquid (
912
913
atof (TextFile->Data_separated [i][column].data ()) * Multiplier);
913
914
column += 1 ;
914
915
// TF commented out since we want to use the improved PointDuMux class
915
916
// this->NodeData[i]->temperature =
916
- // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
917
+ // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
917
918
this ->NodeData [i]->setTemperature (
918
919
atof (TextFile->Data_separated [i][column].data ()) * Multiplier);
919
920
column += 1 ;
920
921
// TF commented out since we want to use the improved PointDuMux class
921
922
// this->NodeData[i]->phase_density[0] =
922
- // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
923
+ // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
923
924
this ->NodeData [i]->getPhaseDensity ()[0 ] =
924
925
atof (TextFile->Data_separated [i][column].data ()) * Multiplier;
925
926
column += 1 ;
926
927
// TF commented out since we want to use the improved PointDuMux class
927
928
// this->NodeData[i]->phase_density[1] =
928
- // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
929
+ // atof(TextFile->Data_separated[i][column].data()) * Multiplier;
929
930
this ->NodeData [i]->getPhaseDensity ()[1 ] =
930
931
atof (TextFile->Data_separated [i][column].data ()) * Multiplier;
931
932
}
@@ -980,10 +981,11 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
980
981
tempstring = temp.str ();
981
982
// TF commented out since we want to use the improved PointDuMux class
982
983
// temp.str(""); temp.clear(); temp << this->NodeData[i]->x;
983
- // tempstring += "; " + temp.str(); temp.str(""); temp.clear(); temp <<
984
- // this->NodeData[i]->y; tempstring += "; " + temp.str(); temp.str("");
985
- // temp.clear(); temp << this->NodeData[i]->z; tempstring += "; " +
986
- // temp.str();
984
+ // tempstring += "; " + temp.str(); temp.str(""); temp.clear();
985
+ // temp
986
+ // << this->NodeData[i]->y; tempstring += "; " + temp.str();
987
+ // temp.str(""); temp.clear(); temp << this->NodeData[i]->z; tempstring
988
+ // += "; " + temp.str();
987
989
for (size_t k (0 ); k < 3 ; k++)
988
990
{
989
991
temp.str (" " );
@@ -994,29 +996,29 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
994
996
995
997
// TF commented out since we want to use the improved PointDuMux class
996
998
// temp.str(""); temp.clear(); temp.precision(12); temp <<
997
- // this->NodeData[i]->phase_pressure[0]; tempstring
999
+ // this->NodeData[i]->phase_pressure[0]; tempstring
998
1000
// += "; " + temp.str();
999
1001
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1000
- // this->NodeData[i]->phase_pressure[1]; tempstring
1002
+ // this->NodeData[i]->phase_pressure[1]; tempstring
1001
1003
// += "; " + temp.str();
1002
1004
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1003
- // this->NodeData[i]->phase_saturation[0];
1005
+ // this->NodeData[i]->phase_saturation[0];
1004
1006
// tempstring
1005
1007
// += "; " + temp.str();
1006
1008
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1007
- // this->NodeData[i]->phase_saturation[1];
1009
+ // this->NodeData[i]->phase_saturation[1];
1008
1010
// tempstring
1009
1011
// += "; " + temp.str();
1010
1012
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1011
- // this->NodeData[i]->phase_density[0]; tempstring
1013
+ // this->NodeData[i]->phase_density[0]; tempstring
1012
1014
// +=
1013
1015
// "; " + temp.str();
1014
1016
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1015
- // this->NodeData[i]->phase_density[1]; tempstring
1017
+ // this->NodeData[i]->phase_density[1]; tempstring
1016
1018
// +=
1017
1019
// "; " + temp.str();
1018
1020
// temp.str(""); temp.clear(); temp.precision(12); temp <<
1019
- // this->NodeData[i]->CO2inLiquid; tempstring += ";
1021
+ // this->NodeData[i]->CO2inLiquid; tempstring += ";
1020
1022
// "
1021
1023
// + temp.str();
1022
1024
temp.str (" " );
@@ -1057,12 +1059,12 @@ void CDUMUXData::ReadDuMuxData(CRFProcess* m_pcs,
1057
1059
1058
1060
// TF commented out since we want to use the improved PointDuMux class
1059
1061
// concentration_CO2_water = this->NodeData[i]->CO2inLiquid *
1060
- // this->NodeData[i]->phase_density[0] / (this->Molweight_CO2 * 1e-3);
1062
+ // this->NodeData[i]->phase_density[0] / (this->Molweight_CO2 * 1e-3);
1061
1063
// mass_CO2_gas = node_volume *
1062
- // this->NodeData[i]->phase_saturation[1] *
1064
+ // this->NodeData[i]->phase_saturation[1] *
1063
1065
// this->NodeData[i]->phase_density[1];
1064
1066
// mass_CO2_water = node_volume *
1065
- // this->NodeData[i]->phase_saturation[0] * concentration_CO2_water *
1067
+ // this->NodeData[i]->phase_saturation[0] * concentration_CO2_water *
1066
1068
// this->Molweight_CO2 * 0.001;
1067
1069
concentration_CO2_water = this ->NodeData [i]->getCO2InLiquid () *
1068
1070
this ->NodeData [i]->getPhaseDensity ()[0 ] /
@@ -1164,7 +1166,7 @@ void CDUMUXData::WriteDataToGeoSys(CRFProcess* m_pcs)
1164
1166
m_pcs->GetNodeValueIndex (" PRESSURE1" ) + 1 ; // +1... new time level
1165
1167
// TF commented out since we want to use the improved PointDuMux class
1166
1168
// value = this->NodeData[i]->getPhasePressure()[1] -
1167
- // this->NodeData[i]->getPhasePressure()[0];
1169
+ // this->NodeData[i]->getPhasePressure()[0];
1168
1170
value =
1169
1171
pnt_dumux->getPhasePressure ()[1 ] - pnt_dumux->getPhasePressure ()[0 ];
1170
1172
m_pcs->SetNodeValue (i, index, value);
@@ -1291,7 +1293,7 @@ void CDUMUXData::WriteDataToGeoSys(CRFProcess* m_pcs)
1291
1293
// * density_liq[kg/m�] / (Molweight_CO2 [g/mol] * 1e-3 [kg/g]) TF
1292
1294
// commented out since we want to use the improved PointDuMux class
1293
1295
// value = this->NodeData[i]->CO2inLiquid *
1294
- // this->NodeData[i]->phase_density[0] / (this->Molweight_CO2 *
1296
+ // this->NodeData[i]->phase_density[0] / (this->Molweight_CO2 *
1295
1297
// 1e-3);
1296
1298
value = pnt_dumux->getCO2InLiquid () * pnt_dumux->getPhaseDensity ()[0 ] /
1297
1299
(this ->Molweight_CO2 * 1e-3 );
@@ -1307,7 +1309,7 @@ void CDUMUXData::WriteDataToGeoSys(CRFProcess* m_pcs)
1307
1309
// TF commented out since we want to use the improved PointDuMux
1308
1310
// class
1309
1311
// cout << "Node: " << i << " Druck: " <<
1310
- // this->NodeData[i]->phase_pressure[1] << " X_CO2: " <<
1312
+ // this->NodeData[i]->phase_pressure[1] << " X_CO2: " <<
1311
1313
// this->NodeData[i]->CO2inLiquid << " Dichte: " <<
1312
1314
// this->NodeData[i]->phase_density[0] << " C(4): " << value
1313
1315
// << "\n";
@@ -1420,7 +1422,8 @@ void CDUMUXData::WriteDataToGeoSys(CRFProcess* m_pcs)
1420
1422
// position = path.find_last_of("\\");
1421
1423
// path = path.substr(0,position);
1422
1424
// temp.str(""); temp.clear(); temp << timestep; tempstring =
1423
- // temp.str(); string aus_file = path + "\\CheckVelocity_" + Phases[k] +
1425
+ // temp.str(); string aus_file = path + "\\CheckVelocity_" + Phases[k]
1426
+ // +
1424
1427
// "_" + tempstring + ".csv"; ofstream aus;
1425
1428
// aus.open(aus_file.data(),ios::out);
1426
1429
// for (int i = 0; i < vec_string.size(); i++) {
@@ -1655,4 +1658,4 @@ int CDUMUXData::RunDuMux(long Timestep, CRFProcess* m_pcs)
1655
1658
// ToDO:
1656
1659
// - Konstruktor nur einmal aufrufen
1657
1660
// - in st-file: Masse einlesen und mit aktueller Dichte in Volumenstrom
1658
- // umrechnen
1661
+ // umrechnen
0 commit comments