Hemodialysis (HD) is one type of procedure for eliminating toxic chemicals and infusing bicarbonate in patients with end-stage renal disease (ESRD). Research and development in the hemodialyzer industry have, hitherto, depended mostly on empirical evidence to optimize HD therapy. This is often costly and involves numerous clinical trials. Developing a comprehensive time-dependent mathematical model to examine the dynamic exchange of solutes ( − 3