Discretizing burn affect on discretized skin nodes

Author: ntatum94

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

@@ -697,10 +697,25 @@ void Energy::UpdateHeatResistance()
     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)
-    // Rough implementation to correlate burn location to skin location
-    // ?????????????????????????? Wait until we can merge the burn updates to get better segmentation
+    bool isBurnWound = false;
     if (m_data.GetBloodChemistry().GetInflammatoryResponse().HasInflammationSource(CDM::enumInflammationSource::Burn)) {
-      const double burnSurfaceAreaFraction = m_data.GetActions().GetPatientActions().GetBurnWound()->GetTotalBodySurfaceArea().GetValue();
+      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;
+      } else if (index == 2 && burnAction->HasCompartment("LeftArm")) {
+        isBurnWound = true;
+      } else if (index == 3 && burnAction->HasCompartment("RightArm")) {
+        isBurnWound = true;
+      } else if (index == 4 && burnAction->HasCompartment("LeftLeg")) {
+        isBurnWound = true;
+      } else if (index == 5 && burnAction->HasCompartment("RightLeg")) {
+        isBurnWound = true;
+      }
+    }
+    if (isBurnWound) {
+      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]);

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

@@ -374,8 +374,8 @@ void Environment::ProcessActions()
     const double lastConvectiveResistance = m_ClothingToEnvironmentPath->GetResistance(HeatResistanceUnit::K_Per_W);
 
     //Target resistance are based on burn surface area (using "next" values since they were just set to basal values earlier in PreProcess)
-    const double targetRadiativeResistance = nextRadiativeResistance * std::pow(1.0 - burnSurfaceAreaFraction, 4.0);
-    const double targetConvectiveResistance = nextConvectiveResistance * std::pow(1.0 - burnSurfaceAreaFraction, 4.0);
+    const double targetRadiativeResistance = nextRadiativeResistance * std::pow(1.0 - burnSurfaceAreaIntensityFraction, 4.0);
+    const double targetConvectiveResistance = nextConvectiveResistance * std::pow(1.0 - burnSurfaceAreaIntensityFraction, 4.0);
 
     const double rampGain = 1.0e-5; //Smooth the response so that we reach targets over several minutes
 
@@ -389,7 +389,7 @@ void Environment::ProcessActions()
     for (SEThermalCircuitPath* skinToClothing : m_SkinToClothingPaths) {
       double nextSkinToClothingResistance = skinToClothing->GetNextResistance(HeatResistanceUnit::K_Per_W);
       const double lastSkinToClothingResistance = skinToClothing->GetResistance(HeatResistanceUnit::K_Per_W);
-      const double targetSkinToClothingResistance = nextSkinToClothingResistance * std::pow(1.0 - burnSurfaceAreaFraction, 4.0);
+      const double targetSkinToClothingResistance = nextSkinToClothingResistance * std::pow(1.0 - burnSurfaceAreaIntensityFraction, 4.0);
       nextSkinToClothingResistance = lastSkinToClothingResistance + rampGain * (targetSkinToClothingResistance - lastSkinToClothingResistance);
       skinToClothing->GetNextResistance().SetValue(nextSkinToClothingResistance, HeatResistanceUnit::K_Per_W);
     }
@@ -735,38 +735,57 @@ void Environment::CalculateEvaporation()
       auto burn_inflamation = std::find(inflamationSources.begin(), inflamationSources.end(), CDM::enumInflammationSource::Burn);
       if (burn_inflamation != inflamationSources.end()) {
         if (m_data.GetActions().GetPatientActions().HasBurnWound()) {
-          skinWettednessDiffusion = m_data.GetActions().GetPatientActions().GetBurnWound()->GetTotalBodySurfaceArea().GetValue();
-          dClothingResistance_m2_kPa_Per_W = 0.0;
-          fCl = 1.0;
+          bool isBurnWoundLocal = false;
+          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")) {
+              isBurnWoundLocal = true;
+            } else if (index == 2 && burnAction->HasCompartment("LeftArm")) {
+              isBurnWoundLocal = true;
+            } else if (index == 3 && burnAction->HasCompartment("RightArm")) {
+              isBurnWoundLocal = true;
+            } else if (index == 4 && burnAction->HasCompartment("LeftLeg")) {
+              isBurnWoundLocal = true;
+            } else if (index == 5 && burnAction->HasCompartment("RightLeg")) {
+              isBurnWoundLocal = true;
+            }
+          }
+          if (isBurnWoundLocal) {
+            skinWettednessDiffusion = m_data.GetActions().GetPatientActions().GetBurnWound()->GetBurnIntensity();
+            dClothingResistance_m2_kPa_Per_W = 0.0;
+            fCl = 1.0;
+          }
         } else {
           // Ok, so we likely want to evaluate the vector of inflamations
           // This really slows down optimization of adding/removing and accessing InflamationSources
           // If we can not have duplicate InflamationSources of a type then why store it as a vector.
           inflamationSources.erase(burn_inflamation);
         }
-      }
 
-      ///\ToDo:  The units of this constant are incorrect on the CDM--they should be W_Per_m2_Pa (no support for this unit currently)
-      GetEvaporativeHeatTranferCoefficient().SetValue(dEvaporativeHeatTransferCoefficient_W_Per_m2_kPa, HeatConductancePerAreaUnit::W_Per_m2_K);
+        ///\ToDo:  The units of this constant are incorrect on the CDM--they should be W_Per_m2_Pa (no support for this unit currently)
+        GetEvaporativeHeatTranferCoefficient().SetValue(dEvaporativeHeatTransferCoefficient_W_Per_m2_kPa, HeatConductancePerAreaUnit::W_Per_m2_K);
 
-      double dMaxEvaporativePotential = (1.0 / 1000.0) * (m_dWaterVaporPressureAtSkin_Pa - m_dWaterVaporPressureInAmbientAir_Pa) / (dClothingResistance_m2_kPa_Per_W + 1.0 / (fCl * dEvaporativeHeatTransferCoefficient_W_Per_m2_kPa));
-      double dSurfaceArea_m2 = m_Patient->GetSkinSurfaceArea(AreaUnit::m2) * segmentedSkinSurfaceAreaPercents[index];
+        double dMaxEvaporativePotential = (1.0 / 1000.0) * (m_dWaterVaporPressureAtSkin_Pa - m_dWaterVaporPressureInAmbientAir_Pa) / (dClothingResistance_m2_kPa_Per_W + 1.0 / (fCl * dEvaporativeHeatTransferCoefficient_W_Per_m2_kPa));
+        double dSurfaceArea_m2 = m_Patient->GetSkinSurfaceArea(AreaUnit::m2) * segmentedSkinSurfaceAreaPercents[index];
 
-      double dSweatRate_kgPers = 0.0;
-      if (m_data.GetEnergy().HasSweatRate()) {
-        dSweatRate_kgPers = m_data.GetEnergy().GetSweatRate(MassPerTimeUnit::kg_Per_s) / dSurfaceArea_m2;
-      }
-      double dSweatingControlMechanisms = dSweatRate_kgPers * m_dHeatOfVaporizationOfWater_J_Per_kg;
-      double dWettedPortion = 0.0;
-      if (dMaxEvaporativePotential != 0) {
-        dWettedPortion = dSweatingControlMechanisms / dMaxEvaporativePotential;
-      }
-      double dDiffusionOfWater = (1.0 - dWettedPortion) * skinWettednessDiffusion * dMaxEvaporativePotential;
-      double EvaporativeHeatLossFromSkin_W = dSweatingControlMechanisms + dDiffusionOfWater;
+        double dSweatRate_kgPers = 0.0;
+        if (m_data.GetEnergy().HasSweatRate()) {
+          dSweatRate_kgPers = m_data.GetEnergy().GetSweatRate(MassPerTimeUnit::kg_Per_s) / dSurfaceArea_m2;
+        }
+        double dSweatingControlMechanisms = dSweatRate_kgPers * m_dHeatOfVaporizationOfWater_J_Per_kg;
+        double dWettedPortion = 0.0;
+        if (dMaxEvaporativePotential != 0) {
+          dWettedPortion = dSweatingControlMechanisms / dMaxEvaporativePotential;
+        }
+        double dDiffusionOfWater = (1.0 - dWettedPortion) * skinWettednessDiffusion * dMaxEvaporativePotential;
+        double EvaporativeHeatLossFromSkin_W = dSweatingControlMechanisms + dDiffusionOfWater;
 
-      // Set the source
-      envSkinToGround->GetNextHeatSource().SetValue(dSurfaceArea_m2 * EvaporativeHeatLossFromSkin_W, PowerUnit::W);
-      index += 1;
+        // Set the source
+        envSkinToGround->GetNextHeatSource().SetValue(dSurfaceArea_m2 * EvaporativeHeatLossFromSkin_W, PowerUnit::W);
+        index += 1;
+      }
     }
   }
 }

share/data/config/VerificationScenarios.config

@@ -300,5 +300,4 @@ ScenarioTest( BurnWoundTreatment, Scenario=Burncare/BurnWoundTreatment.xml  )
 ScenarioTest( BurnWound10, Scenario=Burncare/BurnWound10.xml  )
 ScenarioTest( BurnWound20, Scenario=Burncare/BurnWound20.xml  )
 ScenarioTest( BurnWound25, Scenario=Burncare/BurnWound25.xml  )
-ScenarioTest( BurnWound40, Scenario=Burncare/BurnWound40.xml  )
-
+ScenarioTest( BurnWound40, Scenario=Burncare/BurnWound40.xml  )