google-labs-jules[bot]
ede2dee772
This commit refactors the organ simulation to enable dynamic interactions between organs, replacing the previous "faked" or hardcoded connections.
Key changes include:
- Major Refactoring: Changed the `Organ::update` method signature to `update(Patient& patient, double deltaTime_s)`, allowing organs to access the shared patient state and other organs. This was propagated to all organ classes.
- Blood Chemistry Model: Introduced a central `Blood` struct in the `Patient` model to track shared resources like oxygen, CO2, glucose, and toxins.
- Organ System Interconnections:
- Lungs & Brain: Lungs now perform gas exchange affecting the blood. The brain consumes O2, produces CO2, and its GCS is affected by hypoxia/hypercapnia.
- Liver-Gallbladder: Gallbladder now receives bile directly from the liver's production rate.
- Digestive System: Stomach passes chyme to the intestines, which absorb glucose into the blood. The pancreas responds to blood glucose changes.
- Renal System: Kidneys' GFR is now influenced by the heart's aortic pressure, and they produce urine that fills the bladder directly.
- Cardiovascular & Neurological: The heart rate responds to hypoxia, and the brain uses live aortic pressure from the heart.
- Comprehensive Test Scenario: Updated the main example to include a 60-second simulation with a meal and a lung injury event to verify the new interconnected system.
This creates a more realistic and scalable physiological simulation framework where organ behaviors are emergent from their interactions.
77 lines
2.6 KiB
C++
77 lines
2.6 KiB
C++
#include "MedicalLib/Esophagus.h"
|
|
#include "MedicalLib/Patient.h"
|
|
#include <random>
|
|
#include <algorithm>
|
|
#include <sstream>
|
|
#include <iomanip>
|
|
|
|
Esophagus::Esophagus(int id)
|
|
: Organ(id, "Esophagus"),
|
|
currentState(PeristalsisState::IDLE),
|
|
lowerEsophagealSphincterTone(20.0) {}
|
|
|
|
void Esophagus::update(Patient& patient, double deltaTime_s) {
|
|
// Simulate a swallow every 15 seconds for demonstration
|
|
static double timeSinceLastSwallow = 0.0;
|
|
timeSinceLastSwallow += deltaTime_s;
|
|
if (timeSinceLastSwallow > 15.0) {
|
|
initiateSwallow(15.0); // Swallow a 15mL bolus
|
|
timeSinceLastSwallow = 0.0;
|
|
}
|
|
|
|
if (!activeBoli.empty()) {
|
|
currentState = PeristalsisState::CONTRACTING;
|
|
const double peristalsisSpeed_cm_per_s = 3.0;
|
|
|
|
for (auto& bolus : activeBoli) {
|
|
bolus.position_cm += peristalsisSpeed_cm_per_s * deltaTime_s;
|
|
}
|
|
|
|
// Remove boli that have passed into the stomach
|
|
activeBoli.erase(std::remove_if(activeBoli.begin(), activeBoli.end(),
|
|
[this](const Bolus& b) { return b.position_cm >= this->length_cm; }),
|
|
activeBoli.end());
|
|
|
|
if (activeBoli.empty()) {
|
|
currentState = PeristalsisState::IDLE;
|
|
}
|
|
} else {
|
|
currentState = PeristalsisState::IDLE;
|
|
}
|
|
|
|
// Fluctuate LES tone slightly
|
|
static std::random_device rd;
|
|
static std::mt19937 gen(rd());
|
|
std::normal_distribution<> d(0, 0.1);
|
|
lowerEsophagealSphincterTone += d(gen) * deltaTime_s;
|
|
lowerEsophagealSphincterTone = std::clamp(lowerEsophagealSphincterTone, 18.0, 25.0);
|
|
}
|
|
|
|
void Esophagus::initiateSwallow(double bolusVolume_mL) {
|
|
activeBoli.push_back({bolusVolume_mL, 0.0});
|
|
}
|
|
|
|
std::string Esophagus::stateToString(PeristalsisState state) const {
|
|
switch (state) {
|
|
case PeristalsisState::IDLE: return "Idle";
|
|
case PeristalsisState::CONTRACTING: return "Contracting";
|
|
case PeristalsisState::RELAXING: return "Relaxing";
|
|
default: return "Unknown";
|
|
}
|
|
}
|
|
|
|
std::string Esophagus::getSummary() const {
|
|
std::stringstream ss;
|
|
ss << "--- Esophagus Summary ---\n"
|
|
<< "State: " << stateToString(currentState) << "\n"
|
|
<< "LES Tone: " << std::fixed << std::setprecision(1) << lowerEsophagealSphincterTone << " mmHg\n"
|
|
<< "Boluses in transit: " << activeBoli.size() << "\n";
|
|
if (!activeBoli.empty()) {
|
|
ss << " - Top bolus position: " << activeBoli.front().position_cm << " / " << length_cm << " cm\n";
|
|
}
|
|
return ss.str();
|
|
}
|
|
|
|
PeristalsisState Esophagus::getCurrentState() const { return currentState; }
|
|
bool Esophagus::isSwallowing() const { return !activeBoli.empty(); }
|