Compiling cpp

KratosMultiphysics · Dec 15, 2024 · 6009262 · 6009262
1 parent 402ea29
commit 6009262
Show file tree

Hide file tree

Showing 5 changed files with 1,235 additions and 1,159 deletions.
diff --git a/...tion/custom_elements/data_containers/low_mach_navier_stokes/low_mach_navier_stokes_data.h b/...tion/custom_elements/data_containers/low_mach_navier_stokes/low_mach_navier_stokes_data.h
@@ -58,10 +58,10 @@ class LowMachNavierStokesData : public FluidElementData<TDim,TNumNodes, true>
 
     NodalVectorData MeshVelocity;
 
-    NodalScalarData HeatFlux;
-
     NodalVectorData BodyForce;
 
+    NodalScalarData HeatFlux;
+
     NodalScalarData Pressure;
     NodalScalarData PressureOldStep1;
     NodalScalarData PressureOldStep2;
@@ -72,10 +72,16 @@ class LowMachNavierStokesData : public FluidElementData<TDim,TNumNodes, true>
 
     NodalScalarData Density;
 
+    double ThermodynamicPressure;
+
+    double ThermodynamicPressureDerivative;
+
     double Conductivity;
 
     double SpecificHeat;
 
+    double HeatCapacityRatio;
+
     double ThermalExpansionCoefficient;
 
     double DeltaTime; // Time increment
@@ -117,10 +123,18 @@ class LowMachNavierStokesData : public FluidElementData<TDim,TNumNodes, true>
         const auto& r_properties = rElement.GetProperties();
         this->FillFromProperties(Conductivity, CONDUCTIVITY, r_properties);
         this->FillFromProperties(SpecificHeat, SPECIFIC_HEAT, r_properties);
+        this->FillFromProperties(HeatCapacityRatio, HEAT_CAPACITY_RATIO, r_properties);
 
         // Fill data from ProcessInfo container
         this->FillFromProcessInfo(DeltaTime, DELTA_TIME, rProcessInfo);
         this->FillFromProcessInfo(DynamicTau, DYNAMIC_TAU, rProcessInfo);
+        if (rProcessInfo.Has(PRESSURE)) {
+            this->FillFromProcessInfo(ThermodynamicPressure, PRESSURE, rProcessInfo);
+            ThermodynamicPressureDerivative = 0.0;
+        } else {
+            ThermodynamicPressure = 0.0;
+            ThermodynamicPressureDerivative = 0.0;
+        }
 
         const auto& r_BDF_vector = rProcessInfo[BDF_COEFFICIENTS];
         bdf0 = r_BDF_vector[0];