EXPERIMENTAL cardiovascular simulation

Research and education only. No clinical validity established. All numerical defaults are exploratory. Baseline has elevated pressure and is not a healthy reference. Independent parameter scenarios are not a timed treatment-response simulation.

ScenarioCO (L/min)Forward SV (mL)Periodicity
baseline7.6993102.658periodic_within_declared_gate
reduced_contractility6.699389.325periodic_within_declared_gate
partial_contractility_restoration7.331097.746periodic_within_declared_gate
Computed final-cycle pressure, flow, volume and pressure-volume plots Computed initialization transients and beatwise cardiac output

Model assumptions and limits

Six compartments; synchronous prescribed elastance; linear vascular compliance; four ideal diode valves; fixed external pressure zero. No atria, baroreflex, renal regulation, septal interaction, respiration or drug pharmacokinetics. Pulmonary venous pressure is a reservoir variable, not measured left-atrial pressure. Metadata checks cannot detect falsely labeled numerical data.

Actual parameters, solver settings and mechanism changes

{
  "experimental": true,
  "clinical_validity": "not established",
  "baseline_interpretation": "Exploratory elevated-pressure reference; not a normal healthy adult",
  "scenario_design": "Independent runs from identical initial volumes; changes applied at simulation time zero",
  "config": {
    "cycles": 60,
    "samples_per_cycle": 400,
    "rtol": 1e-09,
    "atol": 1e-10,
    "max_step_fraction": 0.005,
    "periodicity_tolerance_ml": 0.1
  },
  "solver": "DOP853",
  "all_periodic": true,
  "periodicity_definition": "maximum absolute six-volume difference between last two sampled beat waveforms",
  "scenarios": {
    "baseline": {
      "parameters": {
        "heart_rate": 75.0,
        "emin_lv": 0.06,
        "emax_lv": 2.5,
        "emin_rv": 0.03,
        "emax_rv": 0.8,
        "systolic_fraction": 0.4,
        "compliance_sa": 1.5,
        "compliance_sv": 100.0,
        "compliance_pa": 5.0,
        "compliance_pv": 20.0,
        "resistance_aortic": 0.01,
        "resistance_systemic": 1.0,
        "resistance_tricuspid": 0.01,
        "resistance_pulmonary_valve": 0.01,
        "resistance_pulmonary": 0.1,
        "resistance_mitral": 0.01,
        "unstressed_volumes": [
          10.0,
          700.0,
          2500.0,
          10.0,
          100.0,
          400.0
        ],
        "initial_volumes": [
          140.0,
          850.0,
          3200.0,
          140.0,
          180.0,
          490.0
        ]
      },
      "changes": [],
      "metrics": {
        "heart_rate_bpm": 75.0,
        "stroke_volume_ml": 102.65750116541429,
        "cardiac_output_l_min": 7.699312587406099,
        "cardiac_output_m3_s": 0.0001283218764567683,
        "lv_edv_ml": 177.73484201674057,
        "lv_esv_ml": 75.13490812304825,
        "lv_ejection_fraction": 0.5772640453019819,
        "systemic_systolic_pressure_mmhg": 163.00977266014849,
        "systemic_diastolic_pressure_mmhg": 104.70270377958953,
        "systemic_mean_pressure_mmhg": 133.4991899351408,
        "pulmonary_mean_pressure_mmhg": 23.763094206390836,
        "pulmonary_venous_mean_pressure_mmhg": 10.938101773336614,
        "volume_drift_ml": 1.0913936421275139e-11,
        "final_beat_max_volume_difference_ml": 3.880984422721667e-11
      },
      "audit": {
        "volume_drift_ml": 1.0913936421275139e-11,
        "volume_tolerance_ml": 1e-06
      },
      "connected_CO_l_min": 7.699312587406098,
      "status": "periodic_within_declared_gate"
    },
    "reduced_contractility": {
      "parameters": {
        "heart_rate": 75.0,
        "emin_lv": 0.06,
        "emax_lv": 1.25,
        "emin_rv": 0.03,
        "emax_rv": 0.8,
        "systolic_fraction": 0.4,
        "compliance_sa": 1.5,
        "compliance_sv": 100.0,
        "compliance_pa": 5.0,
        "compliance_pv": 20.0,
        "resistance_aortic": 0.01,
        "resistance_systemic": 1.0,
        "resistance_tricuspid": 0.01,
        "resistance_pulmonary_valve": 0.01,
        "resistance_pulmonary": 0.1,
        "resistance_mitral": 0.01,
        "unstressed_volumes": [
          10.0,
          700.0,
          2500.0,
          10.0,
          100.0,
          400.0
        ],
        "initial_volumes": [
          140.0,
          850.0,
          3200.0,
          140.0,
          180.0,
          490.0
        ]
      },
      "changes": [
        {
          "name": "Reduced LV contractility",
          "target": "emax_lv",
          "factor": 0.5,
          "mechanism": "Maximum LV systolic elastance",
          "rationale": "Exploratory isolated contractility impairment; not a clinical heart-failure diagnosis",
          "limitations": "Exploratory parameter change; not a drug-dose or clinical outcome prediction.",
          "before": 2.5,
          "after": 1.25,
          "time_s": 0
        }
      ],
      "metrics": {
        "heart_rate_bpm": 75.0,
        "stroke_volume_ml": 89.32464638298117,
        "cardiac_output_l_min": 6.6993484787236115,
        "cardiac_output_m3_s": 0.00011165580797872685,
        "lv_edv_ml": 214.0823304749263,
        "lv_esv_ml": 124.70492796195904,
        "lv_ejection_fraction": 0.41749079578258463,
        "systemic_systolic_pressure_mmhg": 142.82899525198815,
        "systemic_diastolic_pressure_mmhg": 91.42413186862798,
        "systemic_mean_pressure_mmhg": 116.4336003039731,
        "pulmonary_mean_pressure_mmhg": 24.170951662333643,
        "pulmonary_venous_mean_pressure_mmhg": 12.998776276028545,
        "volume_drift_ml": 1.000444171950221e-11,
        "final_beat_max_volume_difference_ml": 1.48929757415317e-11
      },
      "audit": {
        "volume_drift_ml": 1.000444171950221e-11,
        "volume_tolerance_ml": 1e-06
      },
      "connected_CO_l_min": 6.6993484787236115,
      "status": "periodic_within_declared_gate"
    },
    "partial_contractility_restoration": {
      "parameters": {
        "heart_rate": 75.0,
        "emin_lv": 0.06,
        "emax_lv": 1.875,
        "emin_rv": 0.03,
        "emax_rv": 0.8,
        "systolic_fraction": 0.4,
        "compliance_sa": 1.5,
        "compliance_sv": 100.0,
        "compliance_pa": 5.0,
        "compliance_pv": 20.0,
        "resistance_aortic": 0.01,
        "resistance_systemic": 1.0,
        "resistance_tricuspid": 0.01,
        "resistance_pulmonary_valve": 0.01,
        "resistance_pulmonary": 0.1,
        "resistance_mitral": 0.01,
        "unstressed_volumes": [
          10.0,
          700.0,
          2500.0,
          10.0,
          100.0,
          400.0
        ],
        "initial_volumes": [
          140.0,
          850.0,
          3200.0,
          140.0,
          180.0,
          490.0
        ]
      },
      "changes": [
        {
          "name": "Reduced LV contractility",
          "target": "emax_lv",
          "factor": 0.5,
          "mechanism": "Maximum LV systolic elastance",
          "rationale": "Exploratory isolated contractility impairment; not a clinical heart-failure diagnosis",
          "limitations": "Exploratory parameter change; not a drug-dose or clinical outcome prediction.",
          "before": 2.5,
          "after": 1.25,
          "time_s": 0
        },
        {
          "name": "Partial contractility restoration",
          "target": "emax_lv",
          "factor": 1.5,
          "mechanism": "Maximum LV systolic elastance",
          "rationale": "Illustrative increase relative to impaired parameter; independent run, not a timed drug treatment",
          "limitations": "Exploratory parameter change; not a drug-dose or clinical outcome prediction.",
          "before": 1.25,
          "after": 1.875,
          "time_s": 0
        }
      ],
      "metrics": {
        "heart_rate_bpm": 75.0,
        "stroke_volume_ml": 97.74630508633064,
        "cardiac_output_l_min": 7.330972881474824,
        "cardiac_output_m3_s": 0.0001221828813579137,
        "lv_edv_ml": 190.91920978583036,
        "lv_esv_ml": 93.13592611139043,
        "lv_ejection_fraction": 0.5121710056527649,
        "systemic_systolic_pressure_mmhg": 155.73361779190523,
        "systemic_diastolic_pressure_mmhg": 99.85324668896926,
        "systemic_mean_pressure_mmhg": 127.28245313676504,
        "pulmonary_mean_pressure_mmhg": 23.908034224849835,
        "pulmonary_venous_mean_pressure_mmhg": 11.685124307873833,
        "volume_drift_ml": 4.3655745685100555e-11,
        "final_beat_max_volume_difference_ml": 2.9103830456733704e-11
      },
      "audit": {
        "volume_drift_ml": 4.3655745685100555e-11,
        "volume_tolerance_ml": 1e-06
      },
      "connected_CO_l_min": 7.330972881474824,
      "status": "periodic_within_declared_gate"
    }
  },
  "annotated_variable_count": 34
}

Selected biological annotations and equations

External concept labels are source-verified; aggregate compartments and several quantities are explicitly local extensions. Project relations are not imported ontology entailments.

VariableAnatomyPropertyUnitEquation
V_lvcavity of left ventricleFluid volumemLdV/dt = Q_in - Q_out
P_lvcavity of left ventricleFluid pressuremmHgP = E(t)*(V - V0)
E_lvheart left ventricleFluid volumetric elastancemmHg/mLE(t) = Emin + (Emax-Emin)*activation(t,HR)
V_rvcavity of right ventricleFluid volumemLdV/dt = Q_in - Q_out
P_rvcavity of right ventricleFluid pressuremmHgP = E(t)*(V - V0)
E_rvheart right ventricleFluid volumetric elastancemmHg/mLE(t) = Emin + (Emax-Emin)*activation(t,HR)
V_saLumped systemic arterial compartmentFluid volumemLdV/dt = Q_in - Q_out
P_saLumped systemic arterial compartmentFluid pressuremmHgP = (V - V0)/C
C_saLumped systemic arterial compartmentFluid volume compliancemL/mmHgC = dV/dP; constant linear reservoir compliance
V_svLumped systemic venous compartmentFluid volumemLdV/dt = Q_in - Q_out
P_svLumped systemic venous compartmentFluid pressuremmHgP = (V - V0)/C
C_svLumped systemic venous compartmentFluid volume compliancemL/mmHgC = dV/dP; constant linear reservoir compliance
V_paLumped pulmonary arterial compartmentFluid volumemLdV/dt = Q_in - Q_out
P_paLumped pulmonary arterial compartmentFluid pressuremmHgP = (V - V0)/C
C_paLumped pulmonary arterial compartmentFluid volume compliancemL/mmHgC = dV/dP; constant linear reservoir compliance
V_pvLumped pulmonary venous compartmentFluid volumemLdV/dt = Q_in - Q_out
P_pvLumped pulmonary venous compartmentFluid pressuremmHgP = (V - V0)/C
C_pvLumped pulmonary venous compartmentFluid volume compliancemL/mmHgC = dV/dP; constant linear reservoir compliance
Q_aorticModel flow path lv to saFluid flow ratemL/sQ = max(P_up-P_down,0)/R
R_aorticModel flow path lv to saFluid flow resistancemmHg*s/mLR = pressure_drop / flow
Q_systemicModel flow path sa to svFluid flow ratemL/sQ = (P_up-P_down)/R
R_systemicModel flow path sa to svFluid flow resistancemmHg*s/mLR = pressure_drop / flow
Q_tricuspidModel flow path sv to rvFluid flow ratemL/sQ = max(P_up-P_down,0)/R
R_tricuspidModel flow path sv to rvFluid flow resistancemmHg*s/mLR = pressure_drop / flow
Q_pulmonary_valveModel flow path rv to paFluid flow ratemL/sQ = max(P_up-P_down,0)/R
R_pulmonary_valveModel flow path rv to paFluid flow resistancemmHg*s/mLR = pressure_drop / flow
Q_pulmonaryModel flow path pa to pvFluid flow ratemL/sQ = (P_up-P_down)/R
R_pulmonaryModel flow path pa to pvFluid flow resistancemmHg*s/mLR = pressure_drop / flow
Q_mitralModel flow path pv to lvFluid flow ratemL/sQ = max(P_up-P_down,0)/R
R_mitralModel flow path pv to lvFluid flow resistancemmHg*s/mLR = pressure_drop / flow
HRheartEvent frequency1/minT_seconds = 60 / HR_per_minute
COheart left ventricleCycle-mean aortic outflowmL/sCO = integral(Q_aortic, complete_beat)/T
SVheart left ventricleLeft ventricular stroke volume per complete beatmLSV = integral(Q_aortic, complete_beat)
V_totalClosed six-compartment blood circulationFluid volumemLV_total = sum(V_i)