Nervous Methodology



The BioGears Nervous System is a preliminary system representing the central nervous system (CNS) and peripheral nervous system (PNS). Currently, the system includes baroreceptor and chemoreceptor feedback models that are used to regulate arterial pressure and blood-gas levels. Additionally, the Nervous System provides for tracking of parameters related to traumatic brain injury events, pain stimulus response, and drug effects. In the future, other feedback mechanisms will be added for improved modeling of homeostatic and crisis states.


The Nervous System in the human body is comprised of the central nervous system (CNS) and the the peripheral nervous system (PNS). The CNS contains the brain and the spinal cord, while the PNS is subdivided into the somatic, autonomic, and enteric nervous systems. The somatic system controls voluntary movement. The autonomic nervous system is further subdivided into the parasympathetic and sympathetic nervous system. The sympathetic nervous system is activated in times of stress and crisis, commonly referred to as the "flight or fight" behaviors. The parasympathetic nervous system controls homeostasis functions during the relaxed state. The enteric nervous system controls the gastrointestinal system. Many of these behaviors are tightly connected to the Endocrine System. In many cases, the Nervous System receptors trigger hormone releases that cause systemic changes, such as dilation and constriction of vessels, heart contractility changes, and respiration changes [134].

BioGears will focus on a select few mechanisms within the Nervous System to provide feedback for both homeostasis and crisis behavior in the body. The list is shown below:

  • Baroreceptors
  • Chemoreceptors
  • Local autoregulation
  • Temperature Effects
  • Generic parasympathetic and sympathetic effects
  • Pain stimulus

Additionally, BioGears will model the following features related to brain function:

  • Simple model of traumatic brain injury (TBI)
  • Intracranial pressure
  • Pupillary response to drugs and TBI

System Design

Background and Scope

The BioGears Nervous System began as a low-fidelity model that tracked metrics associated with brain injury. Over time, feedback models local to other BioGears systems (i.e. baroreceptors, chemoreceptors) were relocated to the Nervous System with the goal of moving towards a centralized Nervous model that serves as a global controller. The primary difficulty with such a model is that most detailed Nervous models rely on differential equations that output nerve firing rates. If the frequencies of these firing rates are too high compared to the time-step used to model the equations, the solution will be unstable. Fortunately, as BioGears stability has been improved, the simulation time-step has been decreased to its current value of 0.02 s. This step size allowed us to recently implement a new chemoreceptor model that generates firing rates on a scale of 10 Hz. Future work will seek to implement a similar update to the Baroreceptor model and to create more synergy between the existing models.

Data Flow

An overview of the data flow in the BioGears Nervous system is shown in Figure 1.

Initialization and Stabilization

The BioGears initialization and stabilization is described in detail in the stabilization section of the System Methodology report. The mean arterial pressure set-point is updated after the Cardiovascular system reaches a homeostatic state.


Baroreceptor Feedback

The baroreceptor mechanism provides rapid negative feedback control of arterial pressure. A drop in arterial pressure is sensed by the baroreceptors and leads to an increase in heart rate, heart elastance, and vessel distensibility. These changes operate with the goal of maintaining arterial pressure at its healthy resting level. The baroreceptor mechanism may be divided into three parts: afferent, CNS, and efferent. The afferent part contains the receptors, which detect changes in the MAP. After that, the CNS portions the response into either sympathetic or parasympathetic. Lastly, the efferent part is used to provide feedback to the vasculature and organs within the Cardiovascular System [247]. The model chosen for BioGears only models the CNS and efferent portion of the baroreceptors. The fidelity of the model does require the modeling of the actual stretch receptors in the aorta or carotid arteries to provide the necessary feedback.

Chemoreceptor Feedback

The chemoreceptor mechanism provides feedback control to regulate the partial pressures of oxygen and carbon dioxide in the blood. The response is divided in to two components: peripheral and central. Peripheral chemoreceptors (which reside in the aortic arch and carotid arteries) are sensitive to deviations in both oxygen and carbon dioxide pressures. The central chemoreptors, so named because of their location in the central nervous system, respond only the carbon dioxide pressure changes. Assuming that each receptor class operates independently, the feedback from the central and peripheral chemoreceptors is summed to generate a respiratory and a cardiovascular response. The respiratory response sets the patient ventilation (which effects both the tidal volume and respiration rate) to restore blood gas levels to their set points. Likewise, a cardiovascular response is initiated by altering the patient heart rate.

Check Pain Stimulus

Pain can be initiated in BioGears via a Pain Stimulus action. The BioGears pain response takes into account both the magnitude of the pain action and the susceptibility of the virtual patient to pain (set in the Patient definition file). The result is an output corresponding to the pain Visual Analog Scale (VAS) on a 0-10 interval. The VAS score dictates the patient cardiovascular and respiratory response to the pain stimulus. Increased epinephrine production is also stimulated.


Check Nervous Status

Events associated with nervous system functionality are evaluated. First, intracranial pressure is checked to determine whether or not to throw the Intracranial Hypertension or Intracranial Hypotension events. Then, the presence of fasciculation (muscle spasms) is evaluated based on the activity of the compounds Sarin and Succinylcholine. The latter drug generally produces transient fasciculation which disappear when neuromuscular block is achieved [9]. For a description of the neurological effects of Sarin, see the Nerve Agent: Sarin section in the Drugs documentation. It should be noted that fasciculation can also be caused by an ion imbalance; however, this functionality is not yet active. When completed, the criteria for this type of event will be calcium deficiency. If the patient's arterial calcium concentration falls below one milligram per deciliter, then fasciculation will be triggered [128].


Pupil modifiers are retrieved from the Drugs System and combined with any modifiers from TBI before being applied to the eyes.

Post Process

There is no system specific functionality for the Nervous System in Post Process.


Assessments in BioGears are data collected and packaged to resemble a report or analysis that might be ordered by a physician. No BioGears assessments are associated with the Nervous System.

Figure 1. The data flow for the Nervous System consists of Preprocess, Process and Postprocess. In Preprocess, the baroreceptor feedback is calculated. In Process, the brain metrics are checked to see if event thresholds have been reached and pupil effects are determined and applied.

Features, Capabilities, and Dependencies


The baroreceptor model implemented in BioGears is adapted from the models described by Ottesen et al [247]. The model is used to drive the current MAP towards the resting set-point of the patient. This is accomplished by calculating the sympathetic and parasympathetic response. The fraction of the baroreceptor response that corresponds to the sympathetic effects is determined from Equation 1. The parsympathetic fraction is shown in Equation 2.

\[ \eta_{s} (P_{a}) = \left[ 1+ \left( \frac{P_{a}}{P_{a,setpoint}} \right)^{ \nu} \right]^{-1} \]

Equation 1.

\[ \eta_{p} (P_{a}) = \left[ 1+ \left( \frac{P_{a}}{P_{a,setpoint}} \right)^{- \nu} \right]^{-1} \]

Equation 2.

Where ν is a parameter that represents the response slope of the baroreceptors, pa is the current MAP, and pa,setpoint is the MAP set-point. An example of the sympathetic and parasympathetic responses as a function of MAP are shown in Figure 1. These were calculated with an assumed MAP set-point of 87 mmHg.

Figure 2. The sympathetic and parasympatheric response fractions are displayed as a function of mean arterial pressure (MAP). Both fractional forms show asymptotic behavior as divergence from the MAP set-point occurs. The response fractions are additive, always summing to a value of 1.0. At homeostasis (MAP equal to the set-point), the fractions are both equal to 0.5.

As described in the cardiovascular methodology report, the BioGears cardiovascular system is initialized according to patient definitions and the stabilized to a homeostatic state. The set-point is the resultant mean arterial pressure following the engine stabilization period. The set-point is adjusted dynamically with certain actions and insults, as shown in Equation 3.

\[ P^{n+1}_{a,setpoint} = P^{n}_{a,setpoint} + \Delta P_{a,drugs} + \Delta P_{a,exercise} \]

Equation 3.

Where Δpa,drugs and Δpa,exercise are changes in MAP due to drugs and exercise, respectively.

The sympathetic and parasympathetic fractional responses are used to determine the response of the following cardiovascular parameters:

  • heart rate
  • heart elastance
  • systemic vascular resistance
  • systemic vascular compliance

This is accomplished by tracking the time-dependent values of each parameter relative to their value at the MAP set-point. The time-dependent behavior of these parameters is presented via a set of ordinary differential equations, as shown in Equations 4-7.

\[ \frac{dx_{HR}}{dt} = \left(- \frac{1}{ \tau_{HR}} \right) \left( -x_{HR} + \alpha_{HR} \eta_{s}(P_{a}) + \beta_{HR} \eta_{p}(P_{a}) + \gamma_{HR} \right) \]

Equation 4.

\[ \frac{dx_{E}}{dt} = \left(- \frac{1}{ \tau_{E}} \right) \left( -x_{E} + \alpha_{E} \eta_{s}(P_{a}) + \gamma_{E} \right) \]

Equation 5.

\[ \frac{dx_{R}}{dt} = \left(- \frac{1}{ \tau_{R}} \right) \left( -x_{R} + \alpha_{R} \eta_{s}(P_{a}) + \gamma_{R} \right) \]

Equation 6.

\[ \frac{dx_{C}}{dt} = \left(- \frac{1}{ \tau_{C}} \right) \left( -x_{C} + \alpha_{C} \eta_{p}(P_{a}) + \gamma_{C} \right) \]

Equation 7.

Where xHR, xE, xR and xC are the relative values of heart rate, heart elastance, vascular resistance and vascular compliance, respectively. τHR, τE, τR and τC are the time constants for heart rate, heart elastance, vascular resistance and vascular compliance, respectively. The remaining α, β and γ parameters are a set of tuning variables used to achieve the correct responses in the Cardiovascular System during arterial pressure shifts. Note that the heart rate feedback is a function of both the sympathetic response and parasympathetic response, whereas the elastance feedback and vascular tone feedback depend on the sympathetic or parasympathetic responses individually. Figure 3 shows the normalized response curves.

Figure 3. The plot array demonstrates the normalized organ responses to sympathetic or parasympathetic activity, plotted against the normalized mean arterial pressure.


The chemoreceptors are cells that are sensitive to changes in blood gas concentration. Peripheral chemosensitive cells–located in the aortic arch–detect fluctuations in both arterial carbon dioxide and oxygen partial pressures. The central chemoreceptors, which reside in the central nervous system (CNS), respond to changes in cerebral blood pH caused by irregular carbon dioxide levels. It is assumed in this model that carbon dioxide transport across the blood-brain barrier and equilibration therein are rapid processes; that is, it assumed that the central chemoreceptor response is proportional to deviations in arterial carbon dioxide. As sympathetic activators, the chemoreceptors also increase the heart rate and contractility.

Respiratory Control

Control of respiration is based on the model developed by Albanese, Magosso, and Ursino in [347], [198], and [348]. It is assumed that the central and peripheral feedback define deviations from baseline patient ventilation (set in BioGears during system stabilization) and that their contributions are independent (and thus additive).

\[\dot{V_{target}} = \dot{V_{base}} + \delta\dot{V_{c}} + \delta\dot{V_{p}} \]

Equation 8.

Clearly, the blood gas levels are normal, central Vc and peripheral Vp changes in ventilation are 0, and target ventilation Vtarget equals patient baseline ventilation Vbase.

The evolution of the central chemoreceptor feedback is assumed to be a function of arterial carbon dioxide partial pressure (PCO2) deviation from its normal set point (PCO2,set = 40 mmHg).

\[ \frac{d\delta\dot{V_{c}}}{dt} = \left(\frac{1}{ \tau_{c}} \right) \left( -\delta\dot{V_{c}} + g_{c}\left(P_{CO_{2}}(t) - P_{CO_{2,set}}\right)\right) \]

Equation 9.

The parameters τc and τgc are the time constant and control gain associated with the central response, respectively. Equation 9 is identical to that reported in , except that a delay term is omitted. This omission is a result of being currently unable to model delay differential equations in BioGears.

The form of the peripheral feedback in Equation 10 is similar to Equation 9, but the response is dictated by a nerve firing rate that is a function of the combined action of carbon dioxide and oxygen. Equations 11-12 show how this firing rate (f(t)) is calculated.

\[ \frac{d\delta\dot{V_{p}}}{dt} = \left(\frac{1}{ \tau_{p}} \right) \left( -\delta\dot{V_{p}} + g_{p}\left(f(t) - f_{set}\right)\right) \]

Equation 10.

\[ \frac{df}{dt} = \frac{1}{\tau_{f}}\left(-f(t) + \psi\right) \]

Equation 11.

\[ \psi = \frac{f_{max} + f_{min}*exp\left(\frac{P_{O_{2}} - P_{O_{2,half}}}{k_{O_{2}}}\right)}{1 + exp\left(\frac{P_{O_{2}} - P_{O_{2,half}}}{k_{O_{2}}}\right)}*\left[K*ln\left(\frac{P_{CO_{2}}}{P_{CO_{2,set}}}\right) + \gamma\right] \]

Equation 12.

Equation 12 is a sigmoid that outputs a chemoreceptor firing rate (bounded by fmin and fmax) based on the current arterial oxygen pressure (PO2). The variables PO2,half and kO2 are the half-max and slope parameters that define the shape of the sigmoid, a shape which is increasingly adjusted as the arterial carbon dioxide pressure (PCO2) deviates from normal (PCO2,set. K and are tuning parameters required to maintain a steady output when oxygen and carbon dioxide pressures are at their normal values. Values for all parameters in Equations 9-12 are available in [198].

Cardiovascular Control

Only the heart rate effects of cardiovascular control are modeled in the current chemoreceptor implementation, but contractility modification will be included in a future release. Figure 4 shows the chemoreceptor effect on heart rate. In the figure, the abscissa represents the fractional deviation of gas concentration from baseline, while the ordinate shows the resultant change in heart rate as a fraction of patient baseline. The final heart rate modification due to chemoreceptors is the sum of the oxygen and carbon dioxide effects. For example, severe hypoxia and severe hypercapnia will result in a three-fold increase in heart rate (baseline + 2 * baseline).

Figure 4. The BioGears chemoreceptor model is a phenomenological model which elicits a tuned response to hypoxia and/or hypercapnia. A reverse effect is also present, but at a much lesser magnitude.


Traumatic brain injuries are relatively common, affecting millions annually. They are usually defined as a blunt or penetrating injury to the head that disrupts normal brain function. Furthermore, they are classified as either focal (e.g. cerebral contusions, lacerations, and hematomas) or diffuse (e.g. concussions and diffuse axonal injuries) [2]. The scope of the BioGears TBI model is somewhat limited by the low fidelity of the modeled brain. The brain is represented in the current BioGears Cardiovascular circuit with only two resistors and a compliance, all within one compartment. Because the circuit is modeled in this way, TBI is considered as an acute, non-localized, non-penetrating injury from a circuit perspective. Thus, the TBI model can account for intracranial pressure and cerebral blood flow on the basis of the whole brain, but not to specific areas of the brain. However, BioGears does provide for three injury inputs: diffuse, right focal, and left focal. Similarly to the Renal System, the BioGears brain circuit could be expanded to accommodate a higher-fidelity brain model.

Three important metrics are used to evaluate patients with traumatic brain injuries: intracranial pressure (ICP), cerebral blood flow (CBF), and cerebral perfusion pressure (CPP). Patients with brain injuries often exhibit increased intracranial pressure, decreased cerebral blood flow, and a cerebral perfusion pressure outside of a normal range [326]. Cerebral perfusion pressure is defined as

\[ CPP = MAP - ICP \]

Equation 8.

Where MAP is the mean arterial pressure. In order to model these behaviors, the Brain Injury action in BioGears will modify the resistors of the brain circuit, which is shown in Figure 5 below. The brain circuit is a section of the cardiovascular circuit.

Figure 5. The BioGears brain is represented by two resistors and a compliance. The upstream resistor, R1, is connected to the aorta, and the downstream resistor, R2, is connected to the vena cava.

By increasing R1 and R2, the ICP can be increased while CBF decreases. The resistors are tuned based on the severity (on a scale from 0 to 1) of TBI such that ICP is above 25 mmHg and CBF is near 8 mL per 100 grams of brain tissue per minute for the most severe injury.

Pupillary Response

The pupil's diameter is controlled by the autonomic nervous system in order to affect the amount of light entering the eye. This diameter control is carried out by two muscles: the pupilloconstrictor, which is parasympathetically controlled (via the neurotransmitter acetylcholine), and the pupillodilator, which is sympathetically controlled (via the neurotransmitter norepinephrine) [2]. The synaptic pathways for these two pupil effects differ, and thus, any deviation from normal pupil behavior can shed light on any interference, whether that be a brain injury causing increased pressure on a nerve or a drug effect interfering with synaptic transmission. Because pupil examination is informative yet minimally invasive, it is an excellent tool for helping to diagnose neurological problems. The PERRLA assessment is commonly used to classify pupillary behavior: Pupils Equal, Round, Reactive to Light, and Accomodating. Pupillary response in BioGears is modeled by tracking any modification to pupil size and light reactivity for each eye. Both TBI actions and drugs can apply modifiers to pupil size and reactivity.

Pupil size and reactivity can vary patient-to-patient based on many factors, including age, mental and emotional state, and ambient light conditions [381]. Furthermore, pupillary assessments like PERRLA are often qualitative rather than quantitative. Because of this, BioGears models pupillary modifiers rather than discrete pupil sizes and size changes over time. In this way, the pupillary effects can be assessed qualitatively without the need for baseline patient values for pupil size and reactivity.

There are two ways to affect the pupils in BioGears: drug pharmacodynamic effects and TBI. For a discussion of pharmacodynamic effects on pupillary response, see Drugs Methodology. The other method of influencing pupillary response is TBI. Increasing intracranial pressures constrict ocular nerves, causing pupillary disruptions [165]. Because of this relationship between ICP and pupil effects, ICP is the metric on which pupillary modifiers are based. For the pupil size modifier, Equation 9 was developed so that pupil size (ms) begins to see effects when ICP increases above around 20 mmHg, the pressure at which recoverable brain damage is often observed [326], hitting its maximum slope at 22.5 mmHg, and leveling off as ICP approaches 25 mmHg. For pupil light reactivity (mr), the curve seen in Equation 10 was developed so that it remains at 0 until ICP approaches 20 mmHg, then drops off sharply towards -1 as ICP approaches 25 mmHg. For both of these modifiers, a 1 represents "maximal" effect, -1 represents "minimal" effect, and 0 represents no effect. So a pupil size modifier of 1 would mean the pupil size is at its maximum possible diameter.


Equation 9.

\[m_{r} =-.001*10^{.3(ICP-15)} \]

Equation 10.

For diffuse injuries, these modifiers are applied to both eyes. For focal injuries, the modifiers are applied to the eye on the same side as the injury, as this reflects the most common observed behavior [165].

Patient Variability

The baroreceptor reflex gains and time constants are independent of patient configuration. However, set-points are computed after stabilization, and the baroreceptor reflex drives towards a configuration-specific homeostasis. Because TBI uses multiplicative factors to modify the brain section of the cardiovascular circuit, the patient variability of TBI is the same as that seen in the Cardiovascular System. There is no patient variability inherent in BioGears pupillary response model, as it uses modifiers instead of baseline values. A detailed discussion of patient configuration and variability in BioGears is available in the Patient Methodology report.


The baroreceptors interact directly with the cardiovascular system by modifying cardiovascular circuit and heart driver parameters. These include heart rate, left and right heart elastance, systemic vascular resistance and systemic vascular compliance. The cardiovascular parameters are modified by first calculating the normalized change in the parameter. This normalized change is passed into the common data model (CDM) where the cardiovascular system may extract the change and then apply it locally as a fractional adjustment. The resultant hemodynamic changes feedback directly to the MAP, and the updated MAP is used by the Nervous system in the next time slice to compute new normalized baroreceptor effects. TBI also directly modifies the cardiovascular circuit, specifically both the upstream and downstream brain resistors. Finally, the Drugs system can contribute additional pupillary effects, that combine with possible TBI pupil effects via simple summation.

Assumptions and Limitations

The current implementation of the baroreceptor model does not track adjustments in unstressed volume. Currently BioGears does not correctly represent the physiologic unstressed volume under resting conditions. Therefore, these changes to the unstressed volume cannot be reflected in the BioGears engine. This may be addressed in future releases of the engine.

The BioGears TBI model assumes an acute injury to the brain; chronic effects are not considered, and neither are consciousness assessments. The brain is considered to be one compartment and is not subdivided into discrete functional areas. This means that there is no circuit distinction between a diffuse TBI and a Left Focal TBI. The only difference between the TBI actions is the expected pupil effect. TBI metric validation assumes a supine adult. The current TBI model does not consider specific effects of CO2 saturation, O2 saturation, or blood pH on CBF.



Pain Stimulus

The Pain Stimulus action triggers changes to the cardiovascular, respiratory, and endocrine systems that are associated with sympathetic activation. Action initiation requires a severity (0-1 scale) and a location. If multiple locations are specified, the pain severities are summed until the upper boundary of 1 is reached. The severity is assumed to correspond directly to the Visual Analog Scale (VAS) used to quantify pain in clinical settings (e.g. Pain Severity = 0.5 –> VAS = 5). To account for the fact that individuals perceive pain differently, a Pain Susceptibility parameter was added to the BioGears patient definition. This modifier has a range of -1 to 1, with positive values implying heightened sensitivity to pain and negative values indicating blunted sensitivity. If not susceptibility is set in the patient file, then a value of 0 is assumed (i.e. "average susceptibility). The susceptibility parameter is added to the pain severity when a stimulus action is applied. For example, when exposed to a pain stimulus with a severity of 0.5, two patients with pain susceptibilities of -0.1 and 0.2 would report respective pain levels of 0.4 and 0.7 (4 and 7 on the Visual Analog Scale).

The pain scale calculated from the action severity and patient susceptibility is used to establish changes in physiology. Heart rate, respiration rate, and blood pressure all increase as the consolidated pain scale increases. Much of this response is achieved by increasing production of epinephrine in the BioGears Endocrine system. Administration of analgesic drugs (such as Morphine or Fentanyl) lowers the pain score and reverses these physiological trends.

Brain Injury

The Brain Injury action refers to an acute injury to the brain. Whereas in a real patient, a TBI could have manifold effects, both acute and chronic, BioGears only considers the effects of TBI on brain vascular resistance and pupillary response. The effects of the Brain Injury action depend on the severity value from 0 to 1 and a Type (Diffuse, Left Focal, or Right Focal) assigned to the injury. A severity of 0 will result in a multiplicative factor of 1 being applied to both the upstream and downstream resistors in the brain, which equates to no deviation from the normal, uninjured state. For an injury severity of 1, the upstream resistor is assigned a multiplicative factor of 4.775 and the downstream resistor is assigned a multiplicative factor of 30.409. Severity values between 0 and 1 are converted to multiplicative factors linearly. Increased flow resistance results in increased ICP and decreased CBF. A Diffuse-type injury affects both eyes equally, whereas a Left Focal injury affects only the left eye, and a Right Focal injury affects only the right eye.


Intracranial Hypertension

Intracranial Hypertension is triggered when the intracranial pressure exceeds 25 mmHg. The event is reversed when the pressure returns below 24 mmHg.

Intracranial Hypertension

Intracranial Hypotension is triggered when the intracranial pressure drops below 7 mmHg. The event is reversed when the pressure returns above 7.5 mmHg.

Results and Conclusions

Validation - Resting Physiologic State

No resting state physiology validation was completed, because the baroreceptors do not provide feedback during the initial resting and chronic stabilization phases. For more information on the stabilization phases, see System Methodology.

Validation - Actions and Conditions

Baroreceptor Reflex

The baroreceptor reflex is validated through simulation of an acute hemorrhage scenario. This scenario begins with the healthy male patient. After a short time a massive hemorrhage (bleeding rate of ~1000 mL/min) is initiated, and after 30 seconds the bleeding is stopped. The responses for cardiac output, heart rate, systemic vascular resistance and blood volume are shown in Figure 6. The responses shown in the plot are initially driven by decreasing blood volume, resulting in decreased preload and a subsequent reduction in cardiac output. The baroreceptors attempt to maintain the arterial pressure by increasing the heart rate and contractility, increasing the systemic vascular resistance, and reducing vascular compliance. For this specific hemorrhage scenario, there is an expected 30% increase in heart rate, 15-20% decrease in cardiac output, and 10-15% increase in systemic vascular resistance [136]. BioGears baroreceptor feedback does a fair job reaching these values, though it is a bit aggressive in doing so. The heart rate increase observed is much higher than expected, but the changes in cardiac output and systemic vascular resistance are in line with expectations. The validation results are shown in tabular form below. Additional validation related to the baraoreceptor reflex and hemodynamic stability can be found in the Cardiovascular Methodology report.

Figure 6. The hemodynamic response to hemorrhage with feedback from the baroreceptor reflex.

Table 1. The cardiac output, heart rate, systemic vascular resistance and blood volume are shown as a function of time for the acute hemorrhage scenario. There is a noticeable decrease in blood volume and cardiac output due to the fluid loss. The baroreceptor response in this situation leads to an increase in heart rate and systemic vascular resistance in an attempt to maintain arterial pressure.
Action Notes Sampled Scenario Time (s) Heart Rate (beats/min) Cardiac Output (mL/min) Systemic Vascular Resistance (mmHg s/mL)
Hemorrhage 10% blood loss in 30 s 325 Increase ~30% [136] [247] Decrease ~15-20% [136] [247] Increase 10-15% [136] [247]

Chemoreceptor Reflex

The performance of the Chemoreceptor model in BioGears was validated under both hypercanpic and hypoxic conditions. Hypercapnic validation was based on the study of Reynolds, Milhorn, and Holloman [271], in which healthy volunteers breathed air with supernormal carbon dioxide levels (3%, 5%, 6%, and 7% CO2). Likewise, hypoxic validation scenarios used data from Reynolds and Milhorn [270]. Volunteers in this study breathed 7%, 8%, and 9% O2 in nitrogen.

Brain Injury

The Brain Injury action is validated through repeated application and removal of increasing severities of TBI. The scenario begins with a healthy male patient. After a short time, a mild brain injury (Severity = 0.2, Type = Diffuse) is applied, and the patient is allowed to stabilize before the injury state is removed (only one TBI action can be in effect at a time, so adding a Diffuse Severity 0 TBI removes all TBI effects). This process is repeated for a more severe injury (Severity = 0.75, Type = Left Focal) and severe (Severity = 1, Type = Right Focal) brain injury. We expect to see increases in ICP, with the most severe injury resulting in an ICP greater than 25 mmHg and decreases in CBF, with the most severe case approaching 8 mL per 100 grams of brain per minute (108 mL per minute for the validated patient) [212] [326]. The scenario shows good agreement for these values, and the patient even dies after some time with the most severe injury. We expect CPP to either increase above its maximum normal value or decrease below its minimum normal value, but, though we see a drop, it isn't quite as pronounced as expected [326]. We can also see that for the low severity injury, ICP doesn't quite reach the threshold to strongly affect the pupils. For the Left Focal injury, only the left pupil is affected, and for the Right Focal injury, only the right pupil is affected.

Figure 7. Traumatic brain injury response at three different severity levels.

Figure 8. Pupillary response to the same TBI scenario as shown in Figure 5 where increasing severities are applied first as Diffuse, then as Left Focal, then as Right Focal.

Table 2. The validation data for the TBI scenario shows good agreement with expected results.
Action Notes Action Occurrence Time (s) Sampled Scenario Time (s) Intracranial Pressure (mmHg) Cerebral Blood Flow (mL/min) Cerebral Perfusion Pressure (mmHg) Heart Rate (1/min) Respiration Rate (1/min)
Brain Injury Severity 0.2, CPP=MAP-ICP, StandardMale brain mass=1450g 20 600 10% Increase [212] Decrease [326] Decrease [16] 0-10% Decrease [212] 0-10% Decrease [212]
Brain Injury Severity 0 620 900 7-15 mmHg [326] 50-65 mL/100g/min [134] 60-98 mmHg 72 [134] [12.0, 20.0], [13.0, 19.0] [316] [204]
Brain Injury Severity 0.75 920 1700 Increase [326] Decrease [326] Decrease [16] 0-15% Decrease [212] 0-20% Decrease [212]
Brain Injury Severity 0 1720 2000 7-15 mmHg [326] 50-65 mL/100g/min [134] 60-98 mmHg 72 [134] [12.0, 20.0], [13.0, 19.0] [316] [204]
Brain Injury Severity 1 2020 3220 >25 mmHg [326] <8 mL/100g/min [326] Decrease [16] Decrease [212] Decrease [212]


The Nervous System is currently in a preliminary state that contains only a baroreceptor feedback model, basic TBI, and pupillary response. The baroreceptor feedback is used to control rapid changes in arterial pressure by adjusting heart rate, heart elastance, and vascular resistance and compliance. The baroreceptor model has been validated by comparing the BioGears outputs to experimental data for hemorrhage. It currently shows good agreement with the expected trends, but the magnitude of the response is not a strong as the validation data. This is due to a large total vascular compliance, which allows for large changes in blood volume with small changes in pressure. Future adjustments to the cardiovascular circuit would correct the vascular compliance and improve the accuracy of the model. The TBI model shows good agreement for the most prominent TBI metrics, ICP and CBF, for acute brain injuries. Pupillary response behaves as expected and arms BioGears with yet another tool for matching output with clinical data.

Future Work

Coming Soon

  • Chemoreceptor modification of heart contractility
  • Local autoregulation
  • Improved sympathetic and parasympathetic control

Recommended Improvements

  • Including unstressed volume in the tone model
  • Patient variability incorporated into the baroreceptor reflex model constants
  • Positional variability
  • Higher-fidelity brain model with localized injuries
  • Interventions to treat TBI like drainage of cerebrospinal fluid, hypothermic treatment, hyperventilation, and mannitol therapy
  • Consciousness model, including reduced level of consciousness with carbon monoxide and other intoxications
  • Sleep
  • Heart rate and respiration rate effects from TBI



CBF - Cerebral Blood Flow

CNS - Central Nervous System

CPP - Cerebral Perfusion Pressure

ICP - Intracranial Pressure

PERRLA - Pupils Equal, Round, Reactive to Light, and Accommodating

PNS - Peripheral Nervous System

MAP - Mean Arterial Pressure

mmHg - Millimeters mercury

Data Model Implementation