10-05-2022 checking density model

Cristopher-Morales · Cristopher-Morales · commit f49c522175b0 · 2022-05-10T16:29:00.000+02:00
diff --git a/Common/src/CConfig.cpp b/Common/src/CConfig.cpp
@@ -3741,6 +3741,9 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i
     SU2_MPI::Error("The use of FLUID_MIXTURE requires the number of entries for MOLECULAR_WEIGHT and SPECIFIC_HEAT_CP,\n"
                    "to be equal to the number of entries of SCALAR_INIT + 1", CURRENT_FUNCTION);
     }
+    if (Kind_DensityModel != INC_DENSITYMODEL::VARIABLE) {
+    SU2_MPI::Error("The use of FLUID_MIXTURE requires the INC_DENSITY_MODEL option to be VARIABLE", CURRENT_FUNCTION);
+    }
   }
 
   /*--- Overrule the default values for viscosity if the US measurement system is used. ---*/
diff --git a/SU2_CFD/include/fluid/CFluidModel.hpp b/SU2_CFD/include/fluid/CFluidModel.hpp
@@ -132,13 +132,7 @@ class CFluidModel {
       OneOverMeanMolecularWeight += val_scalars[i_scalar]/(molar_masses[i_scalar]/1000);
       val_scalars_sum += val_scalars[i_scalar];
     }
-    /*if(val_scalars_sum>1){
-      cout<<val_scalars_sum<<endl;
-    }*/
     OneOverMeanMolecularWeight += (1 - val_scalars_sum)/(molar_masses[n_scalars]/1000);
-    /*if(OneOverMeanMolecularWeight<0){
-      cout<<OneOverMeanMolecularWeight<<endl;
-    }*/
     return OneOverMeanMolecularWeight;
   }
   /*!

Original file line number	Diff line number	Diff line change
`@@ -3741,6 +3741,9 @@ void CConfig::SetPostprocessing(SU2_COMPONENT val_software, unsigned short val_i`
`3741`	`3741`	`SU2_MPI::Error("The use of FLUID_MIXTURE requires the number of entries for MOLECULAR_WEIGHT and SPECIFIC_HEAT_CP,\n"`
`3742`	`3742`	`"to be equal to the number of entries of SCALAR_INIT + 1", CURRENT_FUNCTION);`
`3743`	`3743`	`}`
	`3744`	`+ if (Kind_DensityModel != INC_DENSITYMODEL::VARIABLE) {`
	`3745`	`+ SU2_MPI::Error("The use of FLUID_MIXTURE requires the INC_DENSITY_MODEL option to be VARIABLE", CURRENT_FUNCTION);`
	`3746`	`+ }`
`3744`	`3747`	`}`
`3745`	`3748`
`3746`	`3749`	`/--- Overrule the default values for viscosity if the US measurement system is used. ---/`
Original file line number	Diff line number	Diff line change
`@@ -132,13 +132,7 @@ class CFluidModel {`
`132`	`132`	`OneOverMeanMolecularWeight += val_scalars[i_scalar]/(molar_masses[i_scalar]/1000);`
`133`	`133`	`val_scalars_sum += val_scalars[i_scalar];`
`134`	`134`	`}`
`135`		`- /*if(val_scalars_sum>1){`
`136`		`- cout<<val_scalars_sum<<endl;`
`137`		`- }*/`
`138`	`135`	`OneOverMeanMolecularWeight += (1 - val_scalars_sum)/(molar_masses[n_scalars]/1000);`
`139`		`- /*if(OneOverMeanMolecularWeight<0){`
`140`		`- cout<<OneOverMeanMolecularWeight<<endl;`
`141`		`- }*/`
`142`	`136`	`return OneOverMeanMolecularWeight;`
`143`	`137`	`}`
`144`	`138`	`/*!`