Original Research ARTICLE
A Lumped-Parameter Subject-Specific Model of Blood Volume Response to Fluid Infusion
- 1Department of Mechanical Engineering, University of Maryland, College Park, MD, USA
- 2Department of Emergency Medicine, Massachusetts General Hospital, Boston, MA, USA
This paper presents a lumped-parameter model that can reproduce blood volume response to fluid infusion. The model represents the fluid shift between the intravascular and interstitial compartments as the output of a hypothetical feedback controller that regulates the ratio between the volume changes in the intravascular and interstitial fluid at a target value (called “target volume ratio”). The model is characterized by only three parameters: the target volume ratio, feedback gain (specifying the speed of fluid shift), and initial blood volume. This model can obviate the need to incorporate complex mechanisms involved in the fluid shift in reproducing blood volume response to fluid infusion. The ability of the model to reproduce real-world blood volume response to fluid infusion was evaluated by fitting it to a series of data reported in the literature. The model reproduced the data accurately with average error and root-mean-squared error (RMSE) of 0.6 and 9.5% across crystalloid and colloid fluids when normalized by the underlying responses. Further, the parameters derived for the model showed physiologically plausible behaviors. It was concluded that this simple model may accurately reproduce a variety of blood volume responses to fluid infusion throughout different physiological states by fitting three parameters to a given dataset. This offers a tool that can quantify the fluid shift in a dataset given the measured fractional blood volumes.
The fluid shift between the intravascular and interstitial compartments is an essential mechanism of homeostasis, and a key determinant of the physiological response to circulatory pathology and medical therapy. The net fluid shift is determined by the summary action occurring across the body's massive network of microvasculature, and the determinants of flow for each microscopic segment are complex, including the permeability of the vessels and the local Starling forces (i.e., hydrostatic and oncotic pressure gradients). Each of these determinants can be, in turn, altered by a wide range of other factors, including the vasomotion of upstream and downstream vessels, lymphatic flow, and a myriad of endocrine and exocrine signals that affect the preceding determinants.
Given the undeniable complexity that underlies fluid shift between the intravascular and interstitial compartments, it can be challenging to mathematically model the fluid shift in an individual subject or patient. Existing models of fluid shift can reproduce the volume changes in the blood and interstitial fluid (Cervera and Moss, 1974; Pirkle and Gann, 1975; Hedlund et al., 1988; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Gyenge et al., 2003; Tatara et al., 2007) and even intracellular fluid (Hedlund et al., 1988; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Ursino and Innocenti, 1997; Gyenge et al., 2003; Fernandez de Canete and Del Saz Huang, 2010; Siam et al., 2013). However, there are so many disparate factors that, without exhaustive and impractical measurements, it is not possible to comprehensively characterize the fluid shift associated with an individual. A reasonable alternative is to use “typical” values for certain parameters (i.e., values that represent average values in a population) (Champion et al., 1975; Lewis, 1986; Wears and Winton, 1990; Mardel et al., 1995; Simpson et al., 1996; Hirshberg et al., 2006; Saito et al., 2013) but this is no longer an individualized model. Recently reported volume kinetic models (Ståhle et al., 1997; Svensén and Hahn, 1997; Drobin and Hahn, 1999, 2002; Hahn, 2010) may offer value via a tradeoff between simplicity and transparency. However, these models still have weakness in terms of physiological transparency. As an example, volume kinetic models dictate that the ratio between the changes in the intravascular and interstitial volumes depends only on their baseline values (see Appendix for details), which is an oversimplification of the known physiology that the ratio actually varies with the level of intravascular and interstitial volumes as well as physiological conditions (Guyton et al., 1975). To overcome the challenges in modeling the fluid shift in an individual, we sought to develop a lumped-parameter model that could reproduce the fluid shift at the macroscopic level without microscopic details between the intravascular and interstitial compartments. Specifically, we sought a model which is simple enough to be fitted to individual subjects with only a rudimentary set of measurements, and at the same time accurate enough to characterize that subject's response to fluid infusion. The ability of the model to reproduce real-world blood volume response to fluid infusion was evaluated by fitting it to a series of data.
Materials and Methods
Development of a Lumped-Parameter Model of Blood Volume Response to Fluid Infusion
The proposed lumped-parameter model involves only three parameters to describe an individual's fluid shift after fluid infusion: target volume ratio (i.e., ratio between the volume changes in the intravascular and interstitial fluids), α; feedback gain (specifying the speed of fluid shift), K; and initial blood volume, VB0.
The first parameter of the proposed model is α, which varies depending on the overall physiological state of the subject. By way of background, the volume of fluid stored in the intravascular, interstitial and intracellular compartments is determined by the vessel permeability and the hydraulic and osmotic pressure gradients at the capillary walls and cell membranes. From microscopic standpoint, the kinetics of a number of ions and proteins as well as the pressure-volume relationships of the fluid compartments determine the hydraulic and osmotic pressure gradients, and therefore, the fluid volume stored in each compartment. However, from macroscopic standpoint, the consequence of the interaction among these complex mechanisms is that the ratio between the volume changes in the intravascular and interstitial fluid is summarized by a constant parameter value, denoted here as α (Guyton et al., 1975). Typically, interstitial fluid volume changes 2–3 times as much as intravascular volume changes (i.e., α = 2–3) when total blood volume is increased, up to a critical blood volume level. Beyond this level, however, majority of the fluid infused to the body is not stored in the intravascular compartment but is shifted to the interstitial compartment (i.e., α≫2–3) (Guyton et al., 1975; Bajwa and Kulshrestha, 2012). In sum, the fluid shift between intravascular and interstitial compartments, which is due to an array of complex physiological processes, may be viewed as the output of a hypothetical feedback controller that regulates the volume changes in the intravascular and interstitial fluid at a target ratio 1:α.
The second parameter of the model is K, which specifies the speed of fluid shift, i.e., the rate of fluid shift between the intravascular and interstitial compartments. K also varies depending on the physiological state of the subject. This parameter is predicated on a feedback control system analogy shown in Figure 1A. The left and right compartments represent the intravascular and interstitial compartments, respectively, while the valve represents a summary of all the physiological processes that produce the fluid shift between the two compartments. In this analogy, it is both fluid infusion (u) and loss (v; e.g., hemorrhage and urine) that act on the intravascular compartment to alter blood volume (VB). The magnitude of the valve opening is a function of the discrepancy between the target versus actual changes in VB: there is an increased rate of fluid shift,q, between the intravascular and interstitial compartments when the discrepancy between target versus actual VB grows larger.
Figure 1. Lumped-parameter model of blood volume response to fluid infusion. (A) Feedback control system analogy via connected compartments. The left and right compartments represent the intravascular and interstitial compartments, while the valve represents a summary of all the physiologic processes that produce the fluid shift between the two compartments. The fluid infusion (u) and loss (v; e.g., hemorrhage and urine) act on the intravascular compartment to alter the blood volume (VB), which in turn alters the interstitial fluid volume (VISF). The magnitude of valve opening is a function of the discrepancy between the target versus actual changes in blood volume. (B) Mathematical model. The fluid infusion (u) and loss (v; e.g., hemorrhage and urine) are the inputs to the model, while the change in blood volume (ΔVB) is the output. The objective of the feedback controller is to retain the 1/(1+α) fraction of the inputted fluid volume in the intravascular compartment while shifting the remaining α/(1 + α) fraction to the interstitial compartment in the steady state. The fluid shift from the intravascular to interstitial compartment (q) acts as feedback control to steer ΔVB to the target change in blood volume (rB).
The third parameter of the model is VB0. Its role is to normalize VB to yield its fractional change from the initial VB ((VB − VB0)/VB0) as the output of the model. Considering that fractional change in VB, rather than its absolute change, is reported in many experimental protocols (Ståhle et al., 1997; Svensén and Hahn, 1997; Drobin and Hahn, 1999, 2002; Hedin and Hahn, 2005; Hahn, 2010; Hahn et al., 2013), the above normalization allows the model to readily analyze experimental data obtained from such protocols.
The representation in Figure 1A can be formalized into the mathematical model (or block diagram) shown in Figure 1B. Both u and v are the inputs to the model, while the change in VB from VB0 (ΔVB = VB − VB0) is the output. The objective of the feedback controller is to retain the 1/(1+α) fraction of the inputted fluid volume in the intravascular compartment while shifting the remaining α/(1+α) fraction to the interstitial compartment in the steady state. q acts as feedback control to steer ΔVB to the target change in VB (rB). In this way, the steady-state volume changes associated with the intravascular and interstitial fluid achieve the 1:α ratio.
The governing equation associated with this mathematical model is derived as follows. The rate of change in VB is given by the following ordinary differential equation:
q is formulated to eliminate the discrepancy between the target versus actual changes in VB. As described above, rB in response to u and v is the 1/(1 + α) fraction of the accumulated resultant fluid intake u − v:
To capture the macroscopic behavior of q with a possibly simplest mathematical expression, q is modeled to be proportional to eB = rB − ΔVB:
This expression physically means that (1) q shifts from the intravascular to interstitial compartment if rB < ΔVB while from interstitial to intravascular compartment if rB > ΔVB, and that (2) q increases as the discrepancy between rB and ΔVB increases. Combining (Equations 1–3) yields the following ordinary differential equation as the lumped-parameter model that relates u − v to ΔVB:
where was used. To rewrite (Equation 4) in terms of the fractional VB response, ΔVB must be normalized by VB0. Dividing both sides of (Equation 4) by VB0 yields (where ):
The parameters in the model (Equation 5) can be uniquely determined if the data are informative (Ljung, 1999). To illustrate, (Equation 5) can be transformed to the following linear regression:
The vector θ in a linear regression is identifiable (Ljung, 1999) (meaning its elements can be uniquely determined provided the data are informative). Hence, all the parameters can be uniquely determined from θ once it is determined via, e.g., the standard least-squares technique (Ljung, 1999) or even a numerical optimization technique (Storn and Price, 1997): K from the first element of θ; VB0 from the second element of θ; and α from the third element of θ using K and VB0 thus determined.
Model Evaluation Method
The lumped-parameter model detailed above is intended to summarize VB response to u and v in an individual subject or patient using just three parameters specific to the individual. However, it is possible that the underlying physiology is too complex to be adequately reproduced by such a simple model. Hence, we sought to evaluate the ability of this lumped-parameter model to reproduce real-world blood volume response to fluid infusion by fitting it to a multitude of different experimental datasets, to assess whether or not such a simple model is capable of accurately reproducing the experimental data.
To analyze the model, a series of datasets reported in the literature were used that had response (versus time) to u and v. In all the datasets used in this study, was measured using hemoglobin as the tracer substance. Assuming that the mass of the tracer substance remains constant, an increase in VB results in a decrease in the concentration of the tracer substance in the blood. Hence, can be computed from the change in the blood hemoglobin concentration as follows (Drobin and Hahn, 1999, 2002; Hahn, 2010):
where Hct(t) is the hematocrit at time t, and Hgb(t) is hemoglobin concentration at time t. The datasets used in this study included (1) response to bolus infusion of fluid under different VB states (Dataset 1 Drobin and Hahn, 1999); (2) response to bolus infusion of crystalloid fluids (Dataset 2 Drobin and Hahn, 2002); and (3) response to bolus infusion of colloid fluids (Dataset 3 Hedin and Hahn, 2005).
From the aforementioned reports (Drobin and Hahn, 1999, 2002; Hedin and Hahn, 2005), the volume of fluid infused and response were extracted at 10 min intervals as the average of the maximum and minimum responses across all the subjects (although the model is intended for use in individual subjects, average responses were considered here because the response associated with each subject could not be extracted from the visually presented data in these reports). The urine excretion rate was estimated by dividing the total urine volume by the study duration (infusion time + post-infusion observation time) based on the simplifying assumption that it remained constant, because only the total urine volumes were available in these reports.
Overall, we studied a total of seven distinct protocols summarized in Table 1: in Dataset 1, there were three protocols (with 0, 450, and 900 ml pre-infusion hemorrhage); in Dataset 2, there were two protocols (with infusions of 0.9% saline and ringer's lactate); and in Dataset 3, there were two protocols (with infusions of 5% albumin and autologous plasma). The measured responses associated with all the datasets are shown in Figure 2.
Figure 2. Measured fractional blood volume responses to fluid used for evaluation of the lumped-parameter model. Ringer's AC, ringer's acetate (Dataset 1; Drobin and Hahn, 1999); Saline,0.9% saline (Dataset 2; Drobin and Hahn, 2002); Ringer's LA, ringer's lactate (Dataset 2; Drobin and Hahn, 2002); Albumin, 5% albumin (Dataset 3; Hedin and Hahn, 2005); Autologous Plasma, autologous plasma (Dataset 3; Hedin and Hahn, 2005).
Fitting the model to each dataset meant estimating the values of the three parameters α, K and VB0 associated with each dataset. First, Equation (5) was discretized into a difference equation using the Euler's method, where TS is sampling interval (Nise, 2010):
Second, an optimization problem was formulated to derive the set of optimal parameters that minimize the error between and its model-reproduced counterpart :
where Θ is a vector of (unknown) parameters, Θ* is the optimal Θ, is model-reproduced (i) at a discrete time instant i (i = 1⋯N) for a protocol. For a given set of Θ, was computed by inputting u and v data to (Equation 8). The optimization problem (Equation 9) was solved using MATLAB's Optimization Toolbox (MathWorks, MA, USA).
First, the ability of the model to reproduce data was assessed by analyzing the goodness of fit associated with each protocol, in terms of sample-by-sample error normalized by the average (e) and root-mean-squared error normalized by the average (RMSNE):
where is the average response over the study duration, and is reproduced by inputting u and v data to (Equation 8) characterized by Θ*.
Second, the intravascular and interstitial fluid volume responses to u and v were qualitatively compared across different protocols. associated with each protocol was estimated by , while the change in interstitial fluid volume was estimated as the numerical integration of q(i, Θ*):
Third, the following hypotheses on the parameters were generated to gauge if they were meaningfully derived by the model fitting: (1) α* decreases as pre-infusion hemorrhage increases; (2) decreases as pre-infusion hemorrhage increases; and (3) α* is larger in crystalloids than in colloids. These hypotheses were tested by examining the values of α* and .
Parametric Sensitivity Analysis
The parametric sensitivity analysis was performed in the frequency domain. The sensitivity functions associated with the model (Equation 5) to VB0, K, and α were derived as the normalized partial derivatives of G(jω), the frequency response function associated with (Equation 5), with respect to VB0, K, and α:
Then, the frequency responses of these sensitivity functions were examined to elucidate how VB0, K, and α are tuned to fit the model to the data on response to u and v.
Table 2 summarizes the results of data analysis, including the estimated parameters for each dataset and the associated error metrics. Overall, the model reproduced the data accurately with average normalized error and RMSNE of 0.6 and 9.5% across crystalloid and colloid fluids.
Figures 3–5 show responses reproduced by the model, together with the measured data, associated with Datasets 1, 2, and 3, respectively. The upper and lower panels present the measured versus model-reproduced responses as well as and .
Figure 3. Model evaluation for Dataset 1 (Drobin and Hahn, 1999). Upper panel, Measured fractional blood volume responses to ringer's acetate (“Data”) under different pre-infusion hemorrhage states, versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 0 ml pre-infusion hemorrhage; 2nd column, 450 ml pre-infusion hemorrhage; 3rd column, 900 ml pre-infusion hemorrhage.
Figure 4. Model evaluation for Dataset 2 (Drobin and Hahn, 2002). Upper panel, Measured fractional blood volume responses to crystalloid fluids (“Data”), versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 0.9% saline; 2nd column, ringer's lactate.
Figure 5. Model evaluation for Dataset 3 (Hedin and Hahn, 2005). Upper panel, Measured fractional blood volume responses to colloid fluids (“Data”), versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 5% albumin; 2nd column, autologous plasma.
Figure 6 illustrates (A) the frequency responses of the model's sensitivity functions (12) to VB0, K, and α and (B) the mechanism of how VB0, K, and α are tuned to fit the model to the data on response to u and v.
Figure 6. (A) Frequency responses of the model's sensitivity functions (Equation 12) to VB0 and α. (B) Mechanism of how VB0 and α are tuned in fitting the model to fractional blood volume response data.
This paper presented a lumped-parameter model that can reproduce VB response to u and v. This section considers the strengths and limitations of the model in terms of its parsimony, accuracy and physiological transparency.
The lumped-parameter model presented in this paper is essentially a minimal model (Equation 5) characterized by only three parameters: α, K, and VB0. That this simple model can be fully characterized just by these parameters and accurately reproduce subject-specific VB responses is its key benefit compared with the existing physiology-based first principles models (Pirkle and Gann, 1975; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Gyenge et al., 2003) whose complexity makes it extremely challenging to reproduce subject-specific VB responses to u and v without exhaustive and impractical measurements.
The results of data analysis (Table 2 and Figures 3–5) support the ability of the lumped-parameter model to accurately reproduce VB responses to u and v. The goodness of fit between measured versus model-reproduced responses is encouraging. These results, derived from five different fluids one of which was also associated with three different pre-infusion hemorrhage states (Dataset 1), suggest that the model may be suitable for summarizing VB responses to a wide range of fluids over diverse VB states.
The proposed lumped-parameter model is a simple and highly structured model, constrained by the physiological principle that the ratio between the volume changes in the intravascular and interstitial fluid is regulated at a target value (Guyton et al., 1975). Thus, it is not obvious if the parameters estimated from the model fitting (Equation 9) would retain the intended physiological meanings. Since the actual values of the model parameters (VB0, K, and α) were not directly measured, the physiological transparency of the estimated parameters could not be established in a strict sense (i.e., by comparing the actual versus estimated parameters). Yet, the physiological “plausibility” of the proposed model is supported by a few observations from the data analysis (Table 2), as detailed below.
First, the values of and α* estimated by fitting the lumped-parameter model to measured responses were comparable, at least in terms of the orders of magnitude, to the known nominal plasma volume (~3.2 l for an 70 kg male; Skillman et al., 1967; Gyenge et al., 1999; Hedin and Hahn, 2005) and the ratio between the volume changes in the intravascular and interstitial fluid due to a perturbation in the body fluid volume (~2.3; Guyton et al., 1975; Lewis, 1986; Wears and Winton, 1990; Simpson et al., 1996): the values of α* derived for the models associated with isotonic crystalloids (Datasets 1–2) ranged between 2.4 and 4.5, while the values of derived for the models (Datasets 1–3) ranged between 3.1 l and 5.6 l.
Second, α* decreased with an increase in the pre-infusion hemorrhage volume (Dataset 1), which is consistent with the physiological principle that more fluid is retained in the intravascular compartment as its volume decreases (Guyton et al., 1975). In fact, comparing the lower panels in Figure 3, the model estimated that q decreases as pre-infusion hemorrhage increases via a decrease in α*. At the same time, also showed physiologically plausible behavior, exhibiting a decreasing trend with an increase in the pre-infusion hemorrhage (Dataset 1), meeting the physiological expectation that pre-infusion hemorrhage would decrease VB0.
Third, α* estimated for colloids (Dataset 3) was smaller than α* estimated for crystalloids (Datasets 1–2), which is consistent with the fact that colloid is better retained in blood than crystalloid due to its relatively large molecule size (Hahn, 2013) (note that this property of colloid should decrease α* since less fluid is shifted to the interstitial compartment). In fact, comparing the lower panels in Figures 3–5, the model estimated that the fluid volume shifted from the intravascular to interstitial compartment relative to the change in VB is largely smaller for colloids than crystalloids (note that the colloid infusion dose is only 40% of the crystalloid infusion dose; Table 1).
In sum, the lumped-parameter model presented in this paper is (1) simple in that it can be characterized by just three parameters; (2) able to reproduce VB responses to u and v; and (3) showed physiologically plausible behaviors in terms of the estimated model parameters. We speculate that the lumped-parameter model may have clinical uses. For example, each time it is fitted to the response of a subject to bolus infusion of fluid, it may (1) characterize VB state of a subject (e.g., by way of ); (2) estimate unmeasured variables of interest, such as the volume of inter-compartmental fluid shift and the change in interstitial fluid volume; and (3) predict future physiological responses of a subject, such as predicting VB response to future u and v. These speculative applications will be evaluated in future work.
Parametric Sensitivity Analysis: Mechanism Underlying Model Fitting
Parametric sensitivity analysis provided in-depth insight as to how the model parameters are tuned to fit the model to the data. The Bode magnitude plots of the sensitivity functions (Equation 12) to VB0, K, and α suggest that the model's response is primarily affected by VB0 in the high frequency regime while comparably affected by VB0 and α in the low frequency regime, with the effect of K limited to mid-frequency regime (Figure 6A). Physically, this means that the model's transient response is largely affected by VB0 (and, to a modest extent, by K) while its steady-state response is affected comparably by both VB0 and α. Thus, in solving (Equation 9), the model tends to fit the transient portion of the response by primarily adjusting VB0, while it tends to fit the steady-state portion of the response by primarily adjusting α once a value of VB0 that best fits the transient response is determined (Figure 6B).
The above mechanistic insight on model fitting may elucidate why VB0 was often estimated to be larger than what is regarded as nominal as detailed below (see, e.g., VB0 estimated for Dataset 2 in Table 2). In response to a fluid bolus, VB first peaks and then decreases as the fluid is distributed to the extravascular compartment. Thus, the real α may first peak and then decrease to a steady-state value as well (Guyton et al., 1975). However, the model assumes that the value of α is fixed at its steady-state (i.e., the lowest) value. The model (Equation 5) dictates that this α, together with actual VB0, yields an increase in the model-reproduced response. To compensate for this discrepancy, VB0 is tuned to be larger than its actual value to decrease the transient response. This may be (at least) a mechanism underlying the large values of VB0 derived for the Datasets 1–2. Compared with crystalloid (Datasets 1–2), VB0 derived for colloid (Dataset 3) were physiologically more plausible, because the infusion dose of colloid was small (Table 1), resulting in less likelihood of volume overload, thereby making the assumption of constant α more realistic.
Despite promising initial results obtained for the proposed lumped-parameter model, there are a few limitations that must be investigated in more depth in future studies. First, the amount of data examined in this study was limited (seven datasets representing average subjects). Considering that the model is ultimately intended for use in individual subjects, the model must be assessed using datasets associated with individual (rather than average) subjects receiving fluid infusion under diverse physiological states. Noting that we only evaluated the model's ability to reproduce VB response to fluid infusion (i.e., testing the model against the same data used for fitting), particular focus must be given to evaluate the model's ability to predict VB response (i.e., testing the model against fresh data not used for fitting). Second, the physiological relevance of the model parameters was not fully established. Future work is required to compare estimated model parameters with their gold standard counterparts to assess the agreement between them. VB0 may be compared with its direct measurement, e.g., based on the spectrophotometric detection of indocyanine green (Henschen et al., 1993). α may be compared with the ratio between ΔVB and ΔVISF in the steady state, which may be derived from the measurements of VB0 and as well as u and v. It is noted that α varies with volume state and is also altered by many physiological factors (Guyton et al., 1975). Thus, α must be interpreted as a resultant (and perhaps) linearized ratio between ΔVB and ΔVISF within the physiological regime presented by the data, especially under the presence of diverse physiological changes. Third, blood hemoglobin concentration measurement can be subject to inaccuracy under large volume changes (Dasselaar et al., 2012). Future work must investigate the influence of the measurement inaccuracy on the quality of the estimated model parameters. Finally, the amount and attributes (e.g., frequency contents) of data required to accurately determine subject-specific model are yet to be examined. In this study, we have shown via the identifiability analysis that the model can be uniquely determined if the dataset is informative, and that model parameters estimated from dataset containing both transient and steady state exhibit physiologically plausible behaviors. Yet, for the model to be clinically useful, model parameters must be reliably determined from minimal data. Future analysis on this aspect will thus be necessary.
RB, AR, and JH contributed to the development and validation of the lumped-parameter model as well as the manuscript preparation.
This material is based on work supported by the Office of Naval Research (ONR) under Grant No. N000141410591 and N000141512018. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the ONR.
Conflict of Interest Statement
The authors have a pending patent related to the lumped-parameter model developed in this study.
Arturson, G., Groth, T., Hedlund, A., and Zaar, B. (1989). Computer simulation of fluid resuscitation in trauma: first pragmatic validation in thermal injury. J. Burn Care Res. 10, 292–299. doi: 10.1097/00004630-198907000-00002
Cervera, A. L., and Moss, G. (1974). Crystalloid distribution following hemorrhage and hemodilution: mathematical model and prediction of optimum volumes for equilibration at normovolemia. J. Trauma Acute Care Surg. 14, 506–520. doi: 10.1097/00005373-197406000-00007
Champion, H. R., Sturdivan, L., Nolan, J., Stega, M., Adams Cowley, R., Sacco, W., et al. (1975). A mathematical model for volume replacement in bleeding patients. J. Surg. Res. 19, 297–302. doi: 10.1016/0022-4804(75)90053-0
Fernandez de Canete, J., and Del Saz Huang, P. (2010). First-principles modeling of fluid and solute exchange in the human during normal and hemodialysis conditions. Comput. Biol. Med. 40, 740–750. doi: 10.1016/j.compbiomed.2010.07.001
Gyenge, C. C., Bowen, B. D., Reed, R. K., and Bert, J. L. (2003). Preliminary model of fluid and solute distribution and transport during hemorrhage. Ann. Biomed. Eng. 31, 823–839. doi: 10.1114/1.1581878
Hedin, A., and Hahn, R. G. (2005). Volume expansion and plasma protein clearance during intravenous infusion of 5% albumin and autologous plasma. Clin. Sci. Lond. 108, 217–224. doi: 10.1042/cs20040303
Hedlund, A., Zaar, B., Groth, T., and Arturson, G. (1988). Computer simulation of fluid resuscitation in trauma. I. Description of an extensive pathophysiological model and its first validation. Comput. Methods Programs Biomed. 27, 7–21. doi: 10.1016/0169-2607(88)90099-5
Hirshberg, A., Hoyt, D. B., and Mattox, K. L. (2006). Timing of fluid resuscitation shapes the hemodynamic response to uncontrolled hemorrhage: analysis using dynamic modeling. J. Trauma 60, 1221–1227. doi: 10.1097/01.ta.0000220392.36865.fa
Mardel, S. N., Simpson, S. H., Kelly, S., Wytch, R., Beattie, T. F., and Menezes, G. (1995). Validation of a computer model of haemorrhage and transcapillary refill. Med. Eng. Phys. 17, 215–218. doi: 10.1016/1350-4533(95)95712-J
Siam, J., Mandel, Y., and Barnea, O. (2013). Optimization of oxygen delivery in fluid resuscitation for hemorrhagic shock: a computer simulation study. Cardiovasc. Eng. Technol. 5, 82–95. doi: 10.1007/s13239-013-0169-z
Simpson, S. H., Menezes, G., Mardel, S. N., Kelly, S., White, R., and Beattie, T. (1996). A computer model of major haemorrhage and resuscitation. Med. Eng. Phys. 18, 339–343. doi: 10.1016/1350-4533(95)00044-5
Volume Kinetics Model
For illustration purposes, a two-compartment volume kinetic model is considered. However, the finding applies to general volume kinetic models. The governing equation is given by (Hahn, 2010):
where vc and vc0 are the blood volume and its baseline value, vt and vt0 the interstitial fluid volume and its baseline value, R0 fluid infusion rate, CL0, CL and CLd the clearances associated with the baseline fluid loss, dilution-dependent elimination and inter-compartmental distribution, respectively. Denoting the changes in the intravascular and interstitial fluid volumes as state variables x1 and x2, (Eq. A1) is reformulated into the following:
where K1 = CL/vc0, K2 = CLd/vc0, K3 = CLd/vt0 and u = R0 − CL0. So, taking the Laplace transform and solving for x1 and x2 yields:
This results in the following:
Thus, the steady-state ratio between x1 and x2 (equivalent to α in our lumped-parameter model) reduces to lims → 0x2(s)/x1(s) = K2/K3 = vt0/vc0, meaning that the ratio between the changes in the intravascular and interstitial fluid volumes is determined by the ratio of their initial values.
Keywords: fluid shifts, blood volume, fluid resuscitation, mathematical model, lumped-parameter model, subject-specific model
Citation: Bighamian R, Reisner AT and Hahn J-O (2016) A Lumped-Parameter Subject-Specific Model of Blood Volume Response to Fluid Infusion. Front. Physiol. 7:390. doi: 10.3389/fphys.2016.00390
Received: 16 May 2016; Accepted: 23 August 2016;
Published: 31 August 2016.
Edited by:Zbigniew R. Struzik, University of Tokyo, Japan
Reviewed by:Yossi Chait, University of Massachusetts, USA
Tsuneo Tatara, Hyogo College of Medicine, Japan
Ofer Barnea, Tel Aviv University, Israel
John K-J. Li, Rutgers University, USA
Copyright © 2016 Bighamian, Reisner and Hahn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jin-Oh Hahn, email@example.com