# Computational assessment of the effects of a pulsatile pump on toxin removal in blood purification

- Ki Moo Lim
^{1}and - Eun Bo Shim
^{1}Email author

**9**:31

https://doi.org/10.1186/1475-925X-9-31

© Lim and Shim; licensee BioMed Central Ltd. 2010

**Received: **31 January 2010

**Accepted: **22 June 2010

**Published: **22 June 2010

## Abstract

### Background

For blood purification systems using a semipermeable membrane, the convective mass transfer by ultrafiltration plays an important role in toxin removal. The increase in the ultrafiltration rate can improve the toxin removal efficiency of the device, ultimately reducing treatment time and cost. In this study, we assessed the effects of pulsatile flow on the efficiency of the convective toxin removal in blood purification systems using theoretical methods.

### Methods

We devised a new mathematical lumped model to assess the toxin removal efficiency of blood purification systems in patients, integrating the mass transfer model for a human body with a dialyser. The human body model consists of a three-compartment model of body fluid dynamics and a two-compartment model of body solute kinetics. We simulated three types of blood purification therapy with the model, hemofiltration, hemodiafiltration, and high-flux dialysis, and compared the simulation results in terms of toxin (urea and beta-2 microglobulin) clearance and the treatment dose delivered under conditions of pulsatile and non-pulsatile pumping. *In vivo* experiments were also performed to verify the model results.

### Results

Simulation results revealed that pulsatile flow improved the convective clearance of the dialyser and delivered treatment dose for all three types of therapy. Compared with the non-pulsatile pumping method, the increases in the clearance of urea and beta-2 microglobulin with pulsatile pumping were highest with hemofiltration treatment (122.7% and 122.7%, respectively), followed by hemodiafiltration (3.6% and 8.3%, respectively), and high-flux dialysis (1.9% and 4.7%, respectively). EKRc and std Kt/V averaged 28% and 23% higher, respectively, in the pulsatile group than in the non-pulsatile group with hemofiltration treatment.

### Conclusions

The pulsatile effect was highly advantageous for all of the toxins in the hemofiltration treatment and for β_{2}-microglobulin in the hemodiafiltration and high-flux dialysis treatments.

## Background

Blood purification systems using semipermeable membranes

Type of blood purification system | Mechanism of toxin removal |
---|---|

Conventional hemodialysis | Dialysis |

HFx | Dialysis + ultrafiltration |

HF | Ultrafiltration |

HDF | Dialysis + ultrafiltration |

Several experimental studies have indicated that the UF rate in blood purification systems can be increased by using a pulsatile pump instead of a non-pulsatile pump [3, 4]. More specifically, our previous study demonstrated that pulsatile flow can generate a greater UF rate than non-pulsatile flow due to the increased transmembrane pressure (TMP) and water permeability of the membrane [3]. However, these studies were limited to basic *in vitro* experiments to observe the UF efficiency of the dialyser device itself. To our knowledge, *in vivo* evaluations of the device and its application to patients have not been performed due to difficulties in developing an appropriate animal model. Although animal experiments regarding urea removal during blood purification treatment are relatively easy to perform, serious difficulties exist in establishing an *in vivo* animal model with a high ? _{2}-microglobulin (B2M) concentration. As noted previously by Cameron [5], B2M, which is another significant toxin in blood purification, accumulates slowly in the body under CRF conditions, and therefore, a great deal of time is required for renal failure associated with a high B2M concentration. However, *in vivo* analysis is inevitable for a reasonable assessment of the toxin removal efficiency.

Here, we use a theoretical method as an alternative to *in vivo* experiment to assess the effects of pulsatile pump on toxin removal. For this purpose, we devised a new mathematical model to assess the toxin removal efficiency of blood purification systems in patients, integrating the mass transfer model for a human body [6] with a dialyser in a blood purification system [2]. The present study focused on whether pulsatile flow can improve the delivered treatment dose for patients during blood purification. We performed a simple *in vivo* clinical test to verify the present simulation model by comparing the simulated urea concentration profile with clinical observations.

Using the mathematical method, we then predicted the dialyser clearances of urea and B2M for various blood purification therapies using pulsatile and non-pulsatile systems. In addition, the treatment dose delivered to patients with CRF was predicted under conditions of pulsatile and non-pulsatile blood purification therapy with a long-term duration of 5 weeks using common blood purification systems, such as hemofiltration (HF), hemodiafiltration (HDF) and high-flux dialysis (HFx).

## Methods

To compare toxin removal efficacy between pulsatile and non-pulsatile systems for the three types of blood purification therapy (*i.e.*, HF, HDF and HFx), we used a computational model with *in vitro* and *in vivo* experiments. The *in vitro* experiment to analyze the contribution of pulsatile flow to the UF rate has already been reported [3] and is summarized in Appendix A1. We identified parameters used in the computational model from the results of the *in vitro* experiment, such as the pressures of the dialyser inlet and outlet of the blood circuit, and the UF coefficient (*K*
_{
uf
}) of the dialyser according to the blood pumping rate (*Q*
_{
b
}) and pumping type (pulsatile *vs*. non-pulsatile). From these parameters, we theoretically calculated the convective clearance for each type of therapy. By inserting the calculated convective clearance into the computational model, we can simulate the hemodynamic and physiological conditions of hemodialysis patients and assess the effects of pump type and pulsatility on toxin removal efficacy of the patients with a long-term duration. To validate the present computational model, the computation results of urea removal in patients with CRF were compared with the results observed *in vivo*.

### Computational model of the diffusion and convection in the dialyser system

The mass transfer inside the body that occurs during blood purification therapy can be explained roughly using body solutes and fluid exchange kinetics, as described by Ursino *et al*. [6]. To simulate the effect of pulsatile flow in a dialyser on the convective clearance of the three types of blood purification therapy, *i.e.*, HF, HDF, and HFx, we modified the dialyser model of Depner and Garred [2] for each therapy and combined this with the human body model of Ursino *et al*. The human body model is explained briefly below, and the dialyser model for each therapy is explained in the following two sections.

*et al*. [6] is presented in Appendix A2 and Table 2, respectively.

Parameter assignment of the mathematical model

Definition | Value | Units | Ref. | |
---|---|---|---|---|

| Diffusive mass transfer coefficient of the dialyser (urea) | 967 | mL/min | Specs. |

| Diffusive mass transfer coefficient of the dialyser (B2M) | 290 | mL/min | Specs. |

| Dialyser sieving coefficient (urea) | 1 | Specs. | |

| Dialyser sieving coefficient (B2M) | 0.8 | Specs. | |

| Transfer coefficient (urea) | 0.77 | L/min | 6 |

| Transfer coefficient (B2M) | 0.077 | L/min | 2 |

| Water-transfer coefficient | 0.24 | L | 6 |

| Equilibrium ratio (urea) | 1 | 6 | |

| Equilibrium ratio (B2M) | 1 | Est. | |

| Plasma water fraction | 0.94 | 6 | |

| Urea generation rate | 6.24 | mg/min | 2 |

### Convective toxin clearance for HF, HDF and HFx treatments

HF therapy works by passing the patient's blood through a dialyser that filters out waste products and water and then adds replacement fluid before returning the blood to the body. The replacement fluid maintains fluid volume in the blood and provides electrolytes. The only mechanism of toxin removal in HF therapy is convective transport through the dialyser membrane. HDF is a method that combines hemofiltration and counter-current dialysis for diffusive transport. The dialysate flows in the direction opposite the blood flow in order to maintain the maximum concentration gradient across the membrane, and the counter-current flow condition increases the efficiency of dialysis. Consequently, HDF uses both replacement fluid and dialysate solution, and the mechanism of toxin removal in HDF is both convective and diffusive transport. HFx therapy also provides both convective and diffusive transport simultaneously, but it uses only dialysate solution, not replacement fluid. The pressure distribution along the length of the dialyser in both the blood and dialysate compartment is an important factor in convective transport. At the dialyser inlet of the blood circuit is a region of filtration, as the blood pressure exceeds the dialysate pressure. Backfiltration occurs at the region of the dialyser outlet of the blood circuit, as the dialysate pressure exceeds the blood pressure.

where *P*
_{
b,in
}and *P*
_{
b,out
}are the mean pressures at the dialyser inlet and outlet of the blood circuit, respectively, *P*
_{
d,in
}and *P*
_{
d,out
}are those of the dialysate circuit, *K*
_{
uf
}is the UF coefficient of the dialyser and *Q*
_{
f
}is the UF rate.

For the HF treatment, *P*
_{
d,in
}and *P*
_{
d,out
}were set to zero and *P*
_{
b,in
}and *P*
_{
b,out
}were obtained from the *in vitro* reference data. Based on these pressures, TMPm was computed using Eq. (1). We calculated the UF rate (*Q*
_{
f
}) in Eq. (2) from the UF coefficient obtained from the *in vitro* experiment and the TMPm. Then, the convective toxin clearance for the HF treatment was calculated with Eq. (A11) shown in the Appendix using the UF rate.

In the HDF treatment, both diffusive and convective clearances exist simultaneously. We set the dialysate flow rate to 300 mL/min as in the *in vitro* experiments. Diffusive toxin clearance for the HDF treatment can be computed from Eq. (A8) in the Appendix using the specified values of the blood and dialysate flow rate through the dialyser. Based on the pressures and UF coefficients obtained from the *in vitro* experiments, UF rate was calculated using Eq. (1) and (2). Then, convective toxin clearance for the HDF treatment was calculated using Eqs. (A9) and (A10) shown in the Appendix.

HFx treatment also uses convection and diffusion simultaneously for toxin removal but does not use replacement fluid, in contrast to other blood purification systems. Therefore, to maintain a constant body fluid volume as in the other therapies, the net UF should be zero by controlling the UF and backfiltration. In this HFx model, we set the mean pressure at the dialyser inlet and outlet of the dialysate circuit equal to those in the blood circuit to assign a zero net UF. The procedure to compute the UF rate is presented in Eqs. (A12)-(A15) in the Appendix. Convective toxin clearance for the HFx treatment is also calculated in Eq. (A11) in the Appendix using the UF rate.

###
*In vivo* test for model validation

*In vivo*experiments were performed to validate the present computational model. First, we measured urea removal during blood purification in the dialysis unit of Kang's Clinic in Seoul, Republic of Korea, during ten hemodialysis sessions using the AK95 roller pump. Detailed conditions of the clinical measurements are presented in Table 3. Then, we performed model simulations under conditions almost identical to the clinical settings with respect to patient body weight, diffusive mass transfer coefficient of dialyser for urea (KoA

_{u}), UF rate and blood and dialysate flow rates. The computed results of urea removal were then compared with the

*in vivo*clinical data.

Patient conditions and system settings for *in vivo* dialysis experiments

Parameters | Value | Dimension |
---|---|---|

Number of patients | 18 | |

Patient body weight | 73 ± 9 | kg |

Predialysis urea concentration | 105 ± 12 | mg/dL |

KoA | 967 | mL/min |

Blood flow rate | 300 | mL/min |

Dialysate flow rate | 500 | mL/min |

UF rate | 0.5 | L/h |

Treatment time | 4 | h |

### Simulation models

The calculated convective toxin clearance was applied to the proposed mathematical model. Three blood purification therapies (*i.e.*, HF, HDF and HFx) were simulated for both pulsatile and non-pulsatile systems according to the blood flow rate. Then, we compared the total performances of urea and B2M for the pulsatile system with those for the non-pulsatile system.

Next, we performed long-term simulations to compare the delivered treatment dose of the pulsatile system with that of the non-pulsatile system. A treatment schedule of three times per week, 4 h per treatment, was applied to each system. The EKRc and the standard (std) *Kt/V* indices were used to quantify the treatment dose. Here, the EKRc and std *Kt/V* reflect the time-averaged concentration (TAC) and mean pre-treatment concentration (MPC), respectively [7–9]. According to the results of previous studies [9], the required minimum values of EKRc and std *Kt/V* for an anuric patient with a total body water volume of 40 L are 11 mL/min and 2, respectively. The indices are calculated using the time-varying urea concentration profiles in the plasma compartment during 5 weeks of simulation. The indices are defined in the Appendix.

## Results

###
*Parameter estimation* and *in vivo* clinical test

*in vitro*experimental result was used to determine the parameters of the computational model, such as the inlet and outlet pressures of the blood circuits, and the UF coefficient of the dialyser according to the blood pumping rate and pump type. Figure 2(a) shows the pressures at the dialyser inlet and outlet of the blood circuit according to the blood flow rate and pump type. Under the same mean flow rate, the dialyser inlet pressure of the pulsatile pump was higher than that of the non-pulsatile pump, whilst the dialyser outlet pressure was almost identical for both pulsatile and non-pulsatile pumps. Figure 2(b) shows the derived UF coefficient according to the blood flow rate and pump type. The UF coefficient of the pulsatile pump was higher than that of the non-pulsatile pump, and the difference in the UF coefficient between the two pumps decreased with increasing blood flow rate.

*in vivo*for short treatment times. We compared the measured data with computational results. As shown in Figure 3, the simulated plasma urea concentration profile was very similar to the clinical results; both showed a concentration rebound phenomenon after 4 h of treatment.

### Toxin removal dynamics

Here, HF, HDF and HFx were simulated according to the pump type and blood flow rate. Then, we assessed the effects of the pump type on the toxin clearance performances for urea and B2M.

Figure 4(b) shows the percent increase in B2M clearance by pulsatile pumping for the three types of therapy. The HF treatment showed the largest clearance increment followed by the HDF treatment over the entire range of flow rates examined, similar to the case of urea clearance. In HF, the percent increase of B2M clearance was identical to that of urea clearance throughout the entire flow rate. In HDF, the B2M clearance increased by a maximum of 8.3% at a flow rate of 300 mL/min and a minimum of 2.2% at a flow rate of 100 mL/min. In HFx, the B2M clearance showed maximum and minimum increases of 4.7% and 1.5% at flow rates of 300 mL/min and 100 mL/min, respectively. In the cases of the HDF and HFx treatments, the increase in B2M clearance by pulsatile pumping was larger than that of urea clearance.

### Comparison of treatment doses

*Kt/V*indices for 12-h weekly treatments in the three types of therapy over the entire range of flow rates examined (100-400 mL/min). For the HF treatment, EKRc and std

*Kt/V*were on average 28% and 23% higher in the pulsatile group than in the non-pulsatile group, respectively. However, the superiority of the pulsatile pump was not remarkable in the HDF or HFx treatments, with increments of <2% in both treatments.

To satisfy the adequacy line of EKRc ≥ 11.0 and std *Kt/V* ≥ 2.0, the range of the adequate blood flow rate was >200 mL/min for the HDF and HFx treatments under both pulsatile and non-pulsatile conditions. However, for the HF treatment, the ranges of blood flow rate to satisfy the adequacy line of EKRc ≥ 11.0 were >300 mL/min and 400 mL/min under pulsatile and non-pulsatile conditions, respectively. In addition, the ranges of the blood flow rate for the HF treatment to guarantee the adequacy line of std *Kt/V* ≥ 2.0 were 350 mL/min and 400 mL/min under pulsatile and non-pulsatile conditions, respectively.

## Discussion

In our previous *in vitro* study, we demonstrated the superiority of a pulsatile pump relative to a non-pulsatile or weak pulsatile pump in terms of the TMP and UF coefficients of the dialyser [3]. Compared to the non-pulsatile pump, the pulsatile pump induced a higher dialyser pressure, such as dialyser inlet pressure and TMP, and a higher UF coefficient at the same mean blood flow rate. The dialyser inlet pressure is influenced by the flow impedance. Therefore, it is clear that the flow wave-form through the dialyser produced by the pulsatile pump induced higher flow impedance than did that produced by the roller pump. As both blood pumps are flow generators and not pressure generators, the pump inducing the higher impedance develops a higher fluid power to maintain the target flow rate, and consequently produces a higher dialyser inlet pressure, which in turn results in a higher TMPm and UF rate. In addition, an instantaneous high pressure developed as the pulsatile pump increased the UF coefficient due to the reduction of membrane layering. However, the difference in the UF coefficient between pulsatile and roller pumps decreased as the pumping rate increased. This implies that an effective pumping frequency exists for reducing membrane layering. In other words, a greater pumping frequency makes the system less efficient in terms of the UF coefficient. These results were limited to basic *in vitro* experiments to determine the UF efficiency of the dialyser device itself.

In this study, theoretical evaluations of the apparatus and its application in patients were performed to assess the effects of pulsatile flow on the efficiency of convective toxin removal in blood purification systems. For this purpose, we devised a new mathematical model, integrating the mass transfer model for the human body [6] with that for the dialyser (Figure 1). For verification of the model, we simulated conventional hemodialysis therapy for 4 h, and the computed urea concentration profile inside the body was compared with clinical observations (Figure 3). Using the verified model, we theoretically predicted the effect of pulsatile pumping on convective toxin clearance (Figure 4) and delivered treatment dose to patients with CRF (Figure 5) for various blood purification therapies, *i.e.*, HF, HDF and HFx.

HF treatment provides only convective toxin removal by UF and injection of replacement fluid, whilst the HDF and HFx treatments provide diffusive toxin removal driven by concentration differences between the blood and dialysate as well as convective toxin removal. The convective toxin removal is dependent on the UF rate and sieving coefficient. The *in vitro* experiment demonstrated that pulsatile pumping improved the UF rate by increasing the UF coefficient and the TMP (Figure 2). The improved UF rate increased the convective toxin removal. Due to the effect of diffusive clearance, the contributions of convective toxin removal to the total clearance in the HDF and HFx treatments were lower than that in the HF treatment. This explains why the increased clearance by pulsatile pumping showed less benefit in the HDF and HFx treatments than in the HF treatment as shown in Figure 4.

The sieving coefficient of the dialyser for B2M was 0.8, whereas the value for urea was 1. These observations indicated that in the HF treatment, B2M removal was 20% less than that of urea under the same UF rate conditions. Thus, the convective clearance for B2M was 20% less than that for urea in HF during which convective toxin clearance occurred. However, HDF and HFx use the mechanisms of diffusive and convective toxin removal simultaneously. In these cases, B2M removal was 20% less than that of urea by the convective effect under the same UF rate conditions. In contrast, the diffusive mass transfer coefficient of the dialyser for B2M [*KoA*
_{
B2M
}in Eq. (A8) in the Appendix] was 70% less than that for urea, which exceeded the 20% decrease compared to urea by the convective effect. Therefore, in the case of HDF and HFx, the relative portion of convective toxin removal to diffusive toxin removal for B2M was more remarkable than that for urea. These observations explain why the percentage increase in total clearance by pulsatile pumping for B2M was more remarkable than that for urea, as shown in Figure 4.

To test the long-term performance of the three blood purification therapies and predict the effects of pump type on their performance, we simulated the therapies for 5 weeks and obtained the EKRc and std *Kt/V* indices for each case (Figure 5). In terms of the EKRc and std *Kt/V* indices, HDF therapy showed the best performance, followed by HFx therapy. HF therapy showed the poorest performance because it provided only convective toxin removal, whilst the others provided convective and diffusive toxin removals using the dialysis solution. However, the pulsatile effect on the delivered dose was best in HF therapy, with neither HDF nor HFx showing a remarkable effect. As shown in Figure 5, the increases in delivered treatment doses of about 25% for both EKRc and std *Kt/V* in HF therapy suggested that use of a pulsatile pump may become a practical alternative to the conventional method with a non-pulsatile pump for HF treatment.

Although we developed a mathematical model of a CRF patient to assess several blood purification methods, this model is not limited to chronic conditions. It can be used to model an acute renal failure patient by modifying the mass transfer coefficients and considering residual renal function.

Several animal studies have shown that a blood purification system using a pulsatile blood pump can effectively maintain a physiologically stable hemodynamic state, *i.e.*, the heart rate, arterial blood pressure, and hematocrit, and is a plausible alternative to systems using a conventional roller pump [10–12]. Furthermore, our *in vitro* experiment showed that the hematocrit was maintained within the normal range during pulsatile pumping. This means that no blood damage occurred due to the pulse peak pressure during pulsatile blood purification therapy.

The use of a pulsatile blood pump in a blood purification system with a semipermeable membrane was suggested to increase convective toxin removal, which could reduce both treatment time and associated costs.

Although we provided simulated results for the test of the effects of pulsatile pump on toxin removal, there are some limitations to our study. First, we have validated our numerical patient model by comparing the simulated results with that the clinical test only using a conventional non-pulsatile system. However, a comparison with the experiment using pulsatile pump was not attempted because of no available pulsatile pump that can be applied to patients among commercialized blood pumps. In addition, the input parameters of the model, such as the pressures at the dialyser inlet and outlet and UF coefficients, were estimated from *in vitro* rather than *in vivo* experiment.

## Conclusions

We proposed a new mathematical model to assess the toxin removal efficiency of blood purification systems in patients, integrating the mass transfer model for a human body with a dialyser. *In vivo* experiments were also performed to verify clinically the model results, respectively. Using the verified model, we simulated the effect of pulsatile pumping on convective toxin clearance and delivered treatment dose to patients with CRF for various blood purification therapies. Compared with non-pusatile pumping method, the increases of urea and B2M clearances by pulsatile pump were highest in HF treatment (122.7% and 122.7%, respectively), followed by HDF (3.6% and 8.3%, respectively) and HFx (1.9% and 4.7%, respectively). EKRc and std *Kt/V* were on average 28% and 23% higher, respectively, in the pulsatile group than in the non-pulsatile group in HF treatment.

The superiority of pulsatile pumping is significant in HF treatment, which relies only on convection for toxin removal. In both HDF and HFx treatments, the effectiveness of pulsatile pumping on B2M removal, which depends more on convection compared to urea removal, was more remarkable than that on urea removal. Finally, pulsatile system will help any blood purification treatments, whose mechanism of toxin removal is mainly convection, get better performance than conventional non-pulsatile system.

## Appendix

### A1. *In vitro* experiment for identification of model parameters

The goal of the *in vitro* experiment was to establish reference data for the variation in TMPm and UF coefficient according to pump type and flow rate. Thus, the parameters such as TMPm and the UF coefficient for a specific condition of flow rate and pump type can be obtained by interpolating the reference data.

### A2. Governing equations of the mathematical model

*V*

_{ ic }), interstitial fluid (

*V*

_{ is }) and plasma (

*V*

_{ pl }) can be expressed as

where *k*
_{
f
}is the osmotic filtration coefficient at the cellular membrane, *O*
_{
ic
}and *O*
_{
is
}are the osmotic concentrations of the intracellular and interstitial compartments, respectively, *F*
_{
a
}is fluid filtration rate at the arterial capillaries, *R*
_{
v
}is the fluid reabsorption rate at the venous capillaries and *Q*
_{
f
}is the UF rate.

*M*

_{ s,ic }) and extracellular (

*M*

_{ s,ex }) compartments are expressed as follows:

*η*

_{ s }is the mass transfer coefficient of solute

*s*at the cellular membrane,

*β*

_{ s }is the equilibrium ratio of solute

*s*,

*C*

_{ s,ic }and

*C*

_{ s,ex }are the intracellular and extracellular concentrations of solute

*s*, respectively,

*G*

_{ s }is the generation rate of solute

*s*and

*J*

_{ s }is the solute transfer rate through the dialyser membrane. We assumed

*G*

_{ urea }to be 6.24 mg/min [2], and

*J*

_{ s }is derived below:

*α*

_{ s }is the Gibbs-Donnan equilibrium ratio of solute

*s*,

*r*is the plasma water fraction and

*K*

_{ s }is the clearance of solute

*s.*HDF or HFx treatments use simultaneous diffusion and convection for toxin removal. Thus, toxin clearance is calculated as the sum of diffusive clearance (

*Kd*

_{ s }) and convective clearance (

*Kc*

_{ s }) as follows:

*KoA*

_{ s }is the diffusive mass transfer coefficient of the dialyser for solute

*s*and

*Q*

_{ b }and

*Q*

_{ d }are the blood and dialysate flow rates through the dialyser, respectively. The convective clearance is expressed as

*Si*

_{ s }is the membrane sieving coefficient of solute

*s. T*, termed the transmittance, represents the mL/min increase in clearance for each mL/min of filtration and can be calculated as follows:

More detailed equations were described previously by Ursino *et al*. [6] and Depner and Garred [2].

### A3. Derivation of the internal UF rate for the HFx treatment

*K*

_{ ufl }is the UF coefficient per unit length and

*L*is the dialyser length. As a simplified approximation, the local pressure distributions along the dialyser fibres in the blood and dialysate circuits are described as follows:

*x*indicates the local position along the dialyser,

*P*

_{ b }(

*x*) and

*P*

_{ d }(

*x*) are the mean pressures at local position

*x*in the blood circuit and dialysate circuit, respectively. The UF rate, which is the primary determinant of the convective clearance of the toxin, is calculated as follows:

where *M* is a local point where neither UF nor backfiltration is generated. The internal UF rate is directly related to the convective toxin clearance in the HFx treatment, and therefore, the convective clearance for the HFx treatment is calculated using Eq. (A9) in the Appendix.

### A4. Definition of criterion indices

*s*(

*TAC*

_{ s }) is calculated by integrating the concentration profile of the solute in the plasma compartments with respect to time as follows:

*EKRc*

_{ s }is the corrected value of

*EKR*

_{ s }for the normalised water volume, 40 L, as follows:

where *V* is the total body fluid.

*s*(

*MPC*

_{ s }) is acquired by averaging the pre-treatment concentrations in the plasma compartment at steady state, and the weekly std

*Kt/V*is then calculated as follows:

where *t* is the total time. All equations are adapted from Ursino *et al*. [6] and Depner and Garred [2].

## Authors' information

KML received Ph.D. from Dept of Biomedical Engineering of Seoul National University, South Korea at 2008 and thereafter worked for the Biosystems Engineering Lab of Kangwon National University, South Korea, as a post doc. And now he is working for the department of Biomedical Engineering of Johns Hopkins Univ. as post doc.

EBS finished undergraduate course at Seoul National University and received Ph.D. from the Mechanical Engineering Dept. of KAIST (Korea Advanced Institute of Science and Technology) at 1994. He received another Ph.D. at the Dept. of Physiology (Medical School) of Kyoto University, Japan, at 2008. He became a full professor of Kangwon National University from 2004.

## Declarations

### Acknowledgements

This work was supported by the NRL (National Research Lab) program of National Research Foundation of Korea (No. ROA-2008-000-20127-0), by the Korea Research Foundation Grant funded by the Korean Government (MEST) (The Regional Research Universities Program/Medical & Bio-Materials Research Center), and by the National Research Foundation of Korea (No. 313-2007-2-D00968).

## Authors’ Affiliations

## References

- Ronco C:
**Hemofiltration and hemodiafiltration.**In*Hemodialysis High-Efficiency Treatments*. Edited by: Bosch P. San Antonio, Churchill Livingstone; 1993:119–133.Google Scholar - Depner T, Garred L:
**Solute transport mechanisms in dialysis.**In*Replacement of Renal Function by Dialysis*. Edited by: W Horl W, Norwell. Kluwer Academic Publishers; 2004:73–93.View ArticleGoogle Scholar - Lim KM, Park JY, Lee JC, Kim JC, Min BG, Kang ET, Shim EB:
**Quantitative analysis of pulsatile flow contribution to ultrafiltration.***Artif Organs*2009,**33**(1):69–73. 10.1111/j.1525-1594.2008.00683.xView ArticleGoogle Scholar - Runge TM, Briceno JC, Sheller ME, Moritz CE, Sloan L, Bohls FO, Ottmers SE:
**Hemodialysis: evidence of enhanced molecular clearance and ultrafiltration volume by using pulsatile flow.***Int J Artif Organs*1993,**16**(9):645–652.Google Scholar - Cameron JS:
**Dialysis arthropathy, amyloidosis and beta-2 microglobulin.***Pediatr Nephrol*1987,**1**(2):224–229. 10.1007/BF00849296View ArticleGoogle Scholar - Ursino M, Coli L, Brighenti C, Chiari L, de Pascalis A, Avanzolini G:
**Prediction of solute kinetics, acid-base status, and blood volume changes during profiled hemodialysis.***Ann Biomed Eng*2000,**28**(2):204–216. 10.1114/1.245View ArticleGoogle Scholar - Jaber BL, Zimmerman DL, Teehan GS, Swedko P, Burns K, Meyer KB, Leypoldt JK:
**Daily hemofiltration for end-stage renal disease: a feasibility and efficacy trial.***Blood Purif*2004,**22**(6):481–489. 10.1159/000081722View ArticleGoogle Scholar - Yamashita AC, Kawanishi H:
**Kinetics and dose of daily hemofiltration.***Blood Purif*2004,**22:**14–19. 10.1159/000081868View ArticleGoogle Scholar - Lee JC, Park CY, Choi SW, Park JY, Kim JC, Yun YM, Lee KK, Sun K, Min BG:
**Computational dose predictions for combined treatment of hemofiltration with weekly hemodialysis.***Comput Methods Programs Biomed*2008,**89**(3):275–281. 10.1016/j.cmpb.2007.11.004View ArticleGoogle Scholar - Lopez-Herce J, Ruperez M, Sanchez C, Garcia C, Garcia E, Rodriguez D, Del Canizo JF:
**Continuous venovenous renal replacement therapy with a pulsatile tubular blood pump: analysis of efficacy parameters.***Artif Organs*2006,**30**(1):64–69. 10.1111/j.1525-1594.2006.00181.xView ArticleGoogle Scholar - Ruperez M, Lopez-Herce J, Sanchez C, Garcia C, Garcia E, Del Francisco CJ:
**Comparison of a tubular pulsatile pump and a volumetric pump for continuous venovenous renal replacement therapy in a pediatric animal model.***ASIO J*2005,**51**(4):372–375. 10.1097/01.mat.0000169112.40062.0eView ArticleGoogle Scholar - Ruperez M, Sanchez C, Garcia C, Garcia E, Lopez-Herce J, Del Canizo FJ, Vigil D:
**Continuous venovenous renal replacement therapy using a pulsatile blood pump.***Pediatr Nephrol*2003,**18**(1):29–32. 10.1007/s00467-002-0991-2View ArticleGoogle Scholar - Morti SM, Zydney AL:
**Protein-membrane interactions during hemodialysis: effects on solute transport.***ASAIO J*1998,**44:**319–326. 10.1097/00002480-199807000-00016View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.