New mathematical model for fluid-glucose-albumin transport in peritoneal dialysis

A mathematical model for fluid transport in peritoneal dialysis is constructed. The model is based on a three-component nonlinear system of two-dimensional partial differential equations for fluid, glucose and albumin transport with the relevant boundary and initial conditions. Non-constant steady-state solutions of the model are studied. The restrictions on the parameters arising in the model are established with the aim to obtain exact formulae for the non-constant steady-state solutions. As the result, the exact formulae for the fluid fluxes from blood to tissue and across the tissue were constructed together with two linear autonomous ODEs for glucose and albumin concentrations. The analytical results were checked for their applicability for the description of fluid-glucose-albumin transport during peritoneal dialysis.

💡 Research Summary

The paper presents a comprehensive mathematical framework for describing the coupled transport of fluid, glucose, and albumin during peritoneal dialysis (PD). Starting from local mass balances, the authors formulate three coupled nonlinear partial differential equations (PDEs) governing the evolution of the interstitial fluid volume fraction (ν), glucose concentration (C_G), and albumin concentration (C_A). Fluid exchange across the tissue (j_U) and from blood to tissue (q_U) are modeled using linear non‑equilibrium thermodynamics, incorporating hydrostatic pressure gradients, osmotic pressure of glucose, and oncotic pressure of albumin. The solute fluxes include both diffusive and convective components, with permeability and reflection coefficients (σ, S) characterizing the capillary wall and tissue barriers.

A crucial closure relation links ν to the interstitial pressure P via an empirically derived monotonic function F(P), reflecting experimental observations of tissue compressibility. Boundary conditions are set at the peritoneal surface (Dirichlet values for pressure and solute concentrations) and at the far tissue edge (Neumann zero‑flux). Initial conditions describe the pre‑dialysis equilibrium state.

To obtain steady‑state solutions, the authors nondimensionalize the system, introduce scaled variables (p, u, w) for pressure and concentration deviations, and reduce the PDEs to ordinary differential equations (ODEs). By imposing the relations S_A = S_TA and S_G = S_TG, the system simplifies dramatically. Assuming a constant tissue porosity (ν = ν_m) leads to a linear second‑order ODE for q_U, whose solution is a combination of exponentials. The constants are fixed by the boundary conditions, yielding explicit formulas for the fluid flux q_U(x) and the associated pressure gradient j_U(x).

The glucose and albumin concentration profiles, u(x) and w(x), satisfy linear autonomous ODEs with source terms driven by the exponential fluid flux. While closed‑form solutions are not generally available, the authors note that these equations can be solved numerically (e.g., with Maple) or, in special cases, expressed via modified Bessel functions. When ν is allowed to vary linearly with position, the fluid‑flux ODE transforms into a modified Bessel equation, and the solution involves I_0 and K_0 functions.

Parameter values are drawn from clinical literature (blood pressure, blood glucose, plasma albumin, permeability coefficients, etc.). Using realistic parameter sets, the analytical steady‑state solutions reproduce key PD outcomes: ultrafiltration volume, glucose removal, and albumin loss. The authors compare these analytical predictions with earlier numerical simulations of a simplified model, demonstrating good agreement and thereby validating the new formulation.

The study concludes that the derived model captures the essential physics of PD, including the interplay of hydrostatic, osmotic, and oncotic forces, as well as the impact of tissue porosity on transport. Its analytical tractability under reasonable parameter constraints makes it a valuable tool for optimizing dialysis prescriptions, designing new dialysis solutions, and exploring patient‑specific treatment strategies.