Making adjustments and adding burn specificity based on review

Author: ntatum94

projects/biogears/libBiogears/include/biogears/engine/Systems/Energy.h

@@ -100,7 +100,7 @@ class BIOGEARS_API Energy : public SEEnergySystem, public BioGearsSystem {
   RunningAverage m_BicarbonateMolarity_mmol_Per_L;
   bool m_packOn;
   double m_previousWeightPack_kg;
-  double m_meanSkinTemperature;
+  double m_Test;
   // Stateless member variable (Set in SetUp())
   double m_dT_s;
   SEPatient* m_Patient;

projects/biogears/libBiogears/src/cdm/patient/actions/SEBurnWound.cpp

@@ -65,6 +65,8 @@ bool SEBurnWound::Load(const CDM::BurnWoundData& in)
 
   if (in.DegreeOfBurn().present()) {
     SetDegreeOfBurn(in.DegreeOfBurn().get());
+  } else {
+    SetDegreeOfBurn(CDM::enumBurnDegree::Third);
   }
 
   if (in.BurnInitiationTime().present()) {

projects/biogears/libBiogears/src/engine/Controller/BioGears.cpp

@@ -5602,25 +5602,25 @@ void BioGears::SetupInternalTemperature()
   double skinCapTuningFactor = 1.0;
 
   SEThermalCircuitNode& TorsoSkin = cIntemperature.CreateNode(BGE::InternalTemperatureNode::InternalTorsoSkin);
-  TorsoSkin.GetTemperature().SetValue(33.0, TemperatureUnit::C);
+  TorsoSkin.GetTemperature().SetValue(34.4, TemperatureUnit::C);
   SEThermalCircuitPath& TorsoSkinToTemperatureGround = cIntemperature.CreatePath(TorsoSkin, Ground, BGE::InternalTemperaturePath::InternalTorsoSkinToTemperatureGround);
   TorsoSkinToTemperatureGround.GetCapacitanceBaseline().SetValue(skinCapTuningFactor * 1.21, HeatCapacitanceUnit::kcal_Per_C);
   TorsoSkin.GetHeatBaseline().SetValue(TorsoSkinToTemperatureGround.GetCapacitanceBaseline().GetValue(HeatCapacitanceUnit::J_Per_K) * TorsoSkin.GetTemperature().GetValue(TemperatureUnit::K), EnergyUnit::J);
 
   SEThermalCircuitNode& HeadSkin = cIntemperature.CreateNode(BGE::InternalTemperatureNode::InternalHeadSkin);
-  HeadSkin.GetTemperature().SetValue(33.0, TemperatureUnit::C);
+  HeadSkin.GetTemperature().SetValue(35.2, TemperatureUnit::C);
   SEThermalCircuitPath& HeadSkinToTemperatureGround = cIntemperature.CreatePath(HeadSkin, Ground, BGE::InternalTemperaturePath::InternalHeadSkinToTemperatureGround);
   HeadSkinToTemperatureGround.GetCapacitanceBaseline().SetValue(skinCapTuningFactor * 0.24, HeatCapacitanceUnit::kcal_Per_C);
   HeadSkin.GetHeatBaseline().SetValue(HeadSkinToTemperatureGround.GetCapacitanceBaseline().GetValue(HeatCapacitanceUnit::J_Per_K) * HeadSkin.GetTemperature().GetValue(TemperatureUnit::K), EnergyUnit::J);
 
   SEThermalCircuitNode& LeftArmSkin = cIntemperature.CreateNode(BGE::InternalTemperatureNode::InternalLeftArmSkin);
-  LeftArmSkin.GetTemperature().SetValue(33.0, TemperatureUnit::C);
+  LeftArmSkin.GetTemperature().SetValue(33.6, TemperatureUnit::C);
   SEThermalCircuitPath& LeftArmSkinToTemperatureGround = cIntemperature.CreatePath(LeftArmSkin, Ground, BGE::InternalTemperaturePath::InternalLeftArmSkinToTemperatureGround);
   LeftArmSkinToTemperatureGround.GetCapacitanceBaseline().SetValue(skinCapTuningFactor * 0.30, HeatCapacitanceUnit::kcal_Per_C);
   LeftArmSkin.GetHeatBaseline().SetValue(LeftArmSkinToTemperatureGround.GetCapacitanceBaseline().GetValue(HeatCapacitanceUnit::J_Per_K) * LeftArmSkin.GetTemperature().GetValue(TemperatureUnit::K), EnergyUnit::J);
 
   SEThermalCircuitNode& RightArmSkin = cIntemperature.CreateNode(BGE::InternalTemperatureNode::InternalRightArmSkin);
-  RightArmSkin.GetTemperature().SetValue(33.0, TemperatureUnit::C);
+  RightArmSkin.GetTemperature().SetValue(33.6, TemperatureUnit::C);
   SEThermalCircuitPath& RightArmSkinToTemperatureGround = cIntemperature.CreatePath(RightArmSkin, Ground, BGE::InternalTemperaturePath::InternalRightArmSkinToTemperatureGround);
   RightArmSkinToTemperatureGround.GetCapacitanceBaseline().SetValue(skinCapTuningFactor * 0.30, HeatCapacitanceUnit::kcal_Per_C);
   RightArmSkin.GetHeatBaseline().SetValue(RightArmSkinToTemperatureGround.GetCapacitanceBaseline().GetValue(HeatCapacitanceUnit::J_Per_K) * RightArmSkin.GetTemperature().GetValue(TemperatureUnit::K), EnergyUnit::J);

projects/biogears/libBiogears/src/engine/Systems/Energy.cpp

@@ -88,6 +88,7 @@ void Energy::Clear()
   m_BloodpH.Reset();
   m_BicarbonateMolarity_mmol_Per_L.Reset();
   m_previousWeightPack_kg = 0.0;
+  m_Test = 0.0;
 
   m_packOn = false;
 }
@@ -104,14 +105,14 @@ void Energy::Initialize()
   //Initialization of other system variables
   /// \cite herman2008physics
   GetCoreTemperature().SetValue(37.0, TemperatureUnit::C);
-  GetSkinTemperatureTorso().SetValue(33.0, TemperatureUnit::C);
-  GetSkinTemperatureHead().SetValue(33.0, TemperatureUnit::C);
-  GetSkinTemperatureLeftArm().SetValue(33.0, TemperatureUnit::C);
-  GetSkinTemperatureRightArm().SetValue(33.0, TemperatureUnit::C);
+  GetSkinTemperatureTorso().SetValue(34.4, TemperatureUnit::C);
+  GetSkinTemperatureHead().SetValue(35.2, TemperatureUnit::C);
+  GetSkinTemperatureLeftArm().SetValue(33.6, TemperatureUnit::C);
+  GetSkinTemperatureRightArm().SetValue(33.6, TemperatureUnit::C);
   GetSkinTemperatureLeftLeg().SetValue(33.0, TemperatureUnit::C);
   GetSkinTemperatureRightLeg().SetValue(33.0, TemperatureUnit::C);
-  GetBurnSkinTemperature().SetValue(33.0, TemperatureUnit::C);
-  GetSkinTemperature().SetValue(33.0, TemperatureUnit::C);
+  GetBurnSkinTemperature().SetValue(34.4, TemperatureUnit::C);
+  GetSkinTemperature().SetValue(34.4, TemperatureUnit::C);
   /// \cite phypers2006lactate
   GetLactateProductionRate().SetValue(1.3, AmountPerTimeUnit::mol_Per_day);
   /// \cite guyton2006medical
@@ -225,6 +226,7 @@ void Energy::AtSteadyState()
 void Energy::PreProcess()
 {
   CalculateMetabolicHeatGeneration();
+  m_data.GetDataTrack().Probe("CTSresistance", m_Test);
   CalculateSweatRate();
   UpdateHeatResistance();
   //ManageEnergyDeficit();
@@ -610,7 +612,7 @@ void Energy::CalculateSweatRate()
   double dWaterVaporPressureInAmbientAir_mmHg = GeneralMath::AntoineEquation(dAirTemperature_C);
   double m_dWaterVaporPressureInAmbientAir_Pa = Convert(dWaterVaporPressureInAmbientAir_mmHg, PressureUnit::mmHg, PressureUnit::Pa);
   // double ambientAtmosphericPressure_Pa = m_data.GetEnvironment().GetConditions().GetAtmosphericPressure().GetValue(PressureUnit::Pa);
-  double maximumEvaporativeCapacity_W = 14.21 * (m_Patient->GetSkinSurfaceArea().GetValue(AreaUnit::m2)) * effectiveClothingEvaporation_im_Per_clo * (133.322 * (std::pow(10, (8.1076 - (1750.286 / (235.0 + (m_skinNodes[0]->GetTemperature(TemperatureUnit::C))))))) - (m_dWaterVaporPressureInAmbientAir_Pa)); // Still needs effective clothing evaporation
+  double maximumEvaporativeCapacity_W = 14.21 * (m_Patient->GetSkinSurfaceArea().GetValue(AreaUnit::m2)) * effectiveClothingEvaporation_im_Per_clo * (133.322 * (std::pow(10, (8.1076 - (1750.286 / (235.0 + (GetSkinTemperature(TemperatureUnit::C))))))) - (m_dWaterVaporPressureInAmbientAir_Pa)); // Still needs effective clothing evaporation
 
   double vaporizationEnergy_J_Per_kg = m_data.GetConfiguration().GetVaporizationEnergy(EnergyPerMassUnit::J_Per_kg);
   double sweatSodiumConcentration_mM = 51.0; /// \cite shirreffs1997whole
@@ -620,7 +622,6 @@ void Energy::CalculateSweatRate()
 
   double currentEvaporativeCapacity_W = sweatHeatTranferCoefficient_W_Per_K * (m_skinNodes[0]->GetTemperature(TemperatureUnit::K));
   if (currentEvaporativeCapacity_W > maximumEvaporativeCapacity_W) {
-
     sweatHeatTranferCoefficient_W_Per_K = maximumEvaporativeCapacity_W / (m_skinNodes[0]->GetTemperature(TemperatureUnit::K));
   }
 
@@ -698,27 +699,32 @@ void Energy::UpdateHeatResistance()
     double bloodDensity_kg_Per_m3 = m_data.GetBloodChemistry().GetBloodDensity().GetValue(MassPerVolumeUnit::kg_Per_m3);
     double bloodSpecificHeat_J_Per_K_kg = m_data.GetBloodChemistry().GetBloodSpecificHeat().GetValue(HeatCapacitancePerMassUnit::J_Per_K_kg);
 
-    double alphaScale = 0.5; // Scaling factor for convective heat transfer from core to skin (35 seems to be near the upper limit before non-stabilization)
+    double alphaScale = 0.42; // Scaling factor for convective heat transfer from core to skin (35 seems to be near the upper limit before non-stabilization)
     bool isBurnWound = false;
+    double burnSurfaceAreaFraction = 0.0;
     if (m_data.GetBloodChemistry().GetInflammatoryResponse().HasInflammationSource(CDM::enumInflammationSource::Burn)) {
       SEBurnWound* burnAction = m_data.GetActions().GetPatientActions().GetBurnWound();
       std::vector<std::string> burnComptVector = burnAction->GetCompartments();
       // Check if burn is on specific compartment. Skip head since burns cannot currently be initialized on the head
       if (index == 0 && burnAction->HasCompartment("Trunk")) {
         isBurnWound = true;
+        burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->getTrunkBurnIntensity();
       } else if (index == 2 && burnAction->HasCompartment("LeftArm")) {
         isBurnWound = true;
+        burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->getLeftArmBurnIntensity();
       } else if (index == 3 && burnAction->HasCompartment("RightArm")) {
         isBurnWound = true;
+        burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->getRightArmBurnIntensity();
       } else if (index == 4 && burnAction->HasCompartment("LeftLeg")) {
         isBurnWound = true;
+        burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->getLeftLegBurnIntensity();
       } else if (index == 5 && burnAction->HasCompartment("RightLeg")) {
         isBurnWound = true;
+        burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->getRightLegBurnIntensity();
       }
     }
     if (isBurnWound) {
       isAnySegmentBurned = true;
-      const double burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->GetBurnIntensity();
       const double resInput = std::min(2.0 * burnSurfaceAreaFraction, 1.0); // Make >50% burn the worse case scenario
       const double targetAlpha = GeneralMath::LinearInterpolator(0.0, resInput, alphaScale, 20.0, resInput);
       const double lastAlpha = 1.0 / (coreToSkinPath->GetResistance(HeatResistanceUnit::K_Per_W) * bloodDensity_kg_Per_m3 * bloodSpecificHeat_J_Per_K_kg * segmentedSkinBloodFlows[index]);
@@ -730,9 +736,9 @@ void Energy::UpdateHeatResistance()
     // When skin blood flow increases, then heat transfer resistance decreases leading to more heat transfer from core to skin.
     // The opposite occurs for skin blood flow decrease.
     double coreToSkinResistance_K_Per_W = 1.0 / (alphaScale * bloodDensity_kg_Per_m3 * bloodSpecificHeat_J_Per_K_kg * segmentedSkinBloodFlows[index]);
-
+    
     coreToSkinResistance_K_Per_W = BLIM(coreToSkinResistance_K_Per_W, 0.0001, 20.0);
-
+    m_Test = coreToSkinResistance_K_Per_W;
     coreToSkinPath->GetNextResistance().SetValue(coreToSkinResistance_K_Per_W, HeatResistanceUnit::K_Per_W);
     index += 1;
   }