JP5992590 | Medical picture system |
WO/2005/084733 | METHOD AND DEVICE FOR CALCULATING A BOLUS AMOUNT |
WO/2016/200970 | SYSTEM AND METHOD FOR TRACKING CHANGES IN AVERAGE GLYCEMIA IN DIABETICS |
GALUZIO PAULO (US)
KOTANKO PETER (US)
KLEWINGHAUS JUERGEN (DE)
THOMPSON DAVID (US)
FRESENIUS MEDICAL CARE DEUTSCHLAND GMBH (DE)
WO2020185829A1 | 2020-09-17 |
EP3842084A1 | 2021-06-30 | |||
US20200294676A1 | 2020-09-17 |
CHERIF ET AL.: "A mathematical model of the four cardinal acid-base disorders", MATHEMATICAL BIOSCIENCES AND ENGINEERING, vol. 17, no. 5, 2020, pages 4457 - 4476
MAHESHWARI ET AL.: "An In Silico Method to Predict Net Calcium Transfer During Hemodialysis", 39TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY (EMBC, 2017, pages 2740 - 2743, XP033152573, DOI: 10.1109/EMBC.2017.8037424
HABRAN ET AL.: "Mathematical Modeling of Extracorporeal CO2 Removal Therapy", MEDICAL & BIOLOGICAL ENGINEERING & COMPUTING, vol. 56, no. 3, 2018, pages 421 - 434
FORSTER ET AL.: "Low-flow CO2 removal integrated into a renal-replacement circuit can reduce acidosis and decrease vasopressor requirements", CRIT CARE, vol. 17, 2013, pages R154, XP021162868, Retrieved from the Internet
What is claimed is: 1. An apparatus, comprising: at least one processor; and a memory coupled to the at least one processor, the memory comprising instructions that, when executed by the at least one processor, cause the at least one processor to: access an acid-base model configured to model acid-base homeostasis of a patient, the acid-base model comprising a patient model, a dialyzer model, and an extracorporeal CO2 removal device (ECCO2RD) model, and determine predicted patient information using the acid-base model. 2. The apparatus of claim 1, the predicted patient information comprising at least one of a blood flow rate (Q), a serum pH level, a pCO2 level, or a HCO3 level. 3. The apparatus of claim 1, the instructions, when executed by the at least one processor, to cause the at least one processor to determine continuous renal replacement therapy (CRRT) parameters to control acid-base status based on the predicted patient information. 4. The apparatus of claim 1, the acid-base model configured to model the regulation of H+, CO2, and HCO3-. 5. The apparatus of claim 1, the patient model configured to model patient physiology having input of blood flow and output of hydrogen ion concentration, carbon dioxide concentration, and bicarbonate concentration. 6. The apparatus of claim 1, the dialyzer model configured to model continuous renal replacement therapy (CRRT). 7. The apparatus of claim 1, the ECCO2RD model configured to model a one-dimensional (1D) diffusion device between blood and air. 8. The apparatus of claim 1, the acid-base model comprising a blood flow circuit flowing from a patient, modeled by the patient model, to a dialyzer, modeled by the dialyzer mode, to an ECCO2RD, modeled by the ECCO2RD model, and back to the patient. 9. The apparatus of claim 1, the blood circuit comprising diffusion at any point in the blood circuit. 10. The apparatus of claim 1, the patient comprising a virtual patient. 11. The apparatus of claim 1, the predicted patient information comprising a treatment recommendation. 12. The apparatus of claim 11, the treatment recommendation comprising a treatment process for an acid-base disorder. 13. A computer-implemented method of acid-base homeostasis analysis, the method comprising, via a processor of a computing device: providing an acid-base model configured to model acid-base homeostasis of a patient, the acid-base model comprising a patient model, a dialyzer model, and an extracorporeal CO2 removal device (ECCO2RD) model, and determining predicted patient information using the acid-base model. 14. The computer-implemented method of claim 13, the predicted patient information comprising at least one of a blood flow rate (Q), a serum pH level, a pCO2 level, or a HCO3 level. 15. The computer-implemented method of claim 13, further comprising prescribing continuous renal replacement therapy (CRRT) parameters to control acid-base status based on the predicted patient information. 16. The computer-implemented method of claim 13, the acid-base model configured to model the regulation of H+, CO2 and HCO3-. 17. The computer-implemented method of claim 13, the patient model configured to model patient physiology having input of blood flow and output of hydrogen ion concentration, carbon dioxide concentration and bicarbonate concentration. 18. The computer-implemented method of claim 13, the dialyzer model configured to model continuous renal replacement therapy (CRRT). 19. The computer-implemented method of claim 13, the ECCO2RD model configured to model a one-dimensional (1D) diffusion device between blood and air. 20. The computer-implemented method of claim 13, the acid-base model comprising a blood flow circuit flowing from a patient, modeled by the patient model, to a dialyzer, modeled by the dialyzer mode, to an ECCO2RD, modeled by the ECCO2RD model, and back to the patient. 21. The computer-implemented method of claim 13, the blood circuit comprising diffusion at any point in the blood circuit. 22. The computer-implemented method of claim 13, the patient comprising a virtual patient. 23. The computer-implemented method of claim 13, the predicted patient information comprising a treatment recommendation. 24. The computer-implemented method of claim 23, the treatment recommendation comprising a treatment process for an acid-base disorder. |
HCO3 reabsorption rat reaction rates or pK a of the buffer system, production (PH) and removal or non-bicarbonate buffering of protons (γ H ). Therapies targeting these parameters may have a strong effect on correcting pH disturbances. Accordingly, these may be primary parameters for pH acid-base homeostasis for healthy individuals. For individuals with metabolic acidosis, primary parameters may include respiratory CO2 removal supplementation or therapy for example, NaHCO3 or HD), hydration reaction rate (K H + HC0 - ), and/or removal of excess protons (for example, through acid-binder supplementation) will be effective.
[0056] In some embodiments, acid-base models, such as implemented via acid-base model logic 132 (see, for example, FIG. 2) may operate to model acid-basis homeostasis for a patient and/or population of patients to determine predicted patient information (e.g., pH level, pCCE level, and/or HCO3 level). In various embodiments, acid-base models may model an acid-base disorder for a certain patient population to determine predicted patient information (e.g., pH level, pCO 2 level, and/or HCO3 level). The acid-base models may be used to determine predicted patient information for certain treatments (predicted treatment information), which may model treatment results for a patient. In some embodiments, a treatment recommendation may include maintaining/regulating acid-base homeostasis. Due to the accuracy of the acid-base models described according to some embodiments in predicting treatment outcomes, the acid- base models may be used to generate useful treatment recommendations for patients compared with conventional systems. In some embodiments, primary parameters may be determined based on predicted patient information (for instance, running an acid-base model to determine primary parameters affecting an acid-base disorder). Treatment recommendations 158 may be targeted to affect primary parameters in order to effectively address a condition. Embodiments are not limited in this context.
[0057] Some embodiments may provide an acid-base model in the form of a physiological acid- base model. In various embodiments, a physiological acid-base model may provide a physiologically-based model describing acid-base homeostasis under normal (or substantially normal) physiologic conditions that may be used, for example, to analyze the effects of pathophysiologic acid-base perturbations on the acid-base status of a patient (or virtual patient). In exemplary embodiments, a physiological acid-base model may operate to model, process, analyze, experiment, or otherwise simulate, among other things, the physiological regulation of HCO 3 – /CO 2 buffering system with Henderson-Hasselbalch mass-action kinetics, endogenous production of both CO2 and H + , non-bicarbonate buffering, and/or renal and respiratory regulations. [0058] Various embodiments may provide acid-base models in the form of intradialytic models. In some embodiments, the intradialytic models may be or may include a dialysis patient model (for instance, implemented via dialysis patient model logic 134) that may include an implementation of a physiological acid-base model for a dialysis patient, for example, characterized by impaired renal regulation that has been replaced by dialysis (for instance, hemodialysis (HD)). [0059] As described in more detail below, the acid-base models according to some embodiments may be used to predict treatment outcomes and/or to provide treatment recommendations for various acid-base disorders. In some embodiments, certain model parameters may be altered to determine their effect on acid-base homeostasis. For instance, changes to an acid secretion rate parameter and a renal filtration rate parameter (for example, reducing or eliminating these parameters (e.g., setting their model value to zero)) may simulate an acid-base disorder, such as renal failure, metabolic acidosis, metabolic alkalosis, respiratory acidosis, respiratory alkalosis, and/or the like. Based on the particular configuration and operating processes of acid-base models according to some embodiments, they may be used to accurately predict changes to parameters affecting acid-base homeostasis. [0060] In some embodiments, acid-base homeostasis logic 130 may operate to receive patient information 152 for a particular patient (for instance, gender, age, health (for example, renal failure, normal, and/or the like), HCO3 level, pH, and/or the like) and determine a treatment recommendation 156 for maintaining acid-base homeostasis and/or treating an acid-base disorder. For example, in a healthy person, correcting a deficient HCO 3 level may be different than for a patient experiencing metabolic acidosis. In some embodiments, acid-base models may be used to model, predict, or otherwise process various treatment models (for instance, acid- binder therapies, HCO 3 therapies, and/or the like) to determine, predict, or otherwise analyze treatment outcomes that may be used to determine actual patient therapy regimens. In some embodiments, acid-base homeostasis logic 130 may include, implement, or otherwise process a feedback loop (or iterative) function for ongoing use of the acid-base models to adjust acid-base conditions, including through dialysis modifications (for example, adjusting ultrafiltration rate and/or volume), or administration of drugs, bicarbonate, etc. based on continuous predicted patient information. In some embodiments, as a treatment is administered (for example, during dialysis, supplement dosages, and/or the like), the treatment information may be used as input into an acid-base model to generate updated predicted patient information, which may then be used to determine treatment results and/or to update a treatment recommendation (including in real or substantially real time, for example, during a dialysis treatment). [0061] FIG.2 illustrates an example of an operating environment 200 that may be representative of some embodiments. As shown in FIG.2, operating environment 200 may include an acid- base model 205 depicting an illustrative model structure in accordance with the present disclosure. For example, FIG.2 may represent the structure between different elements of an acid-base model according to some embodiments. In various embodiments, acid-base model 205 may include a patient model (or element) 210, a dialyzer model 220, an ECCO2RD model 230, and/or a dilution model 240. [0062] The configuration of acid-base model 205 depicted in FIG.2 is for non-limiting illustrative purposes. For example, acid-base model 205 may have different forms, configurations, and/or the like according to some embodiments. For instance, acid-base model 205 may have more or less elements, and/or elements arranged in different configurations. For example, acid-base model may not include dialyzer model 220 and/or ECCO 2 RD model 230. In another example, dialyzer model 220 and/or ECCO2RD model 230 may be in different positions (e.g., as indicated in FIG.2, dialyzer model 220 and ECCO 2 RD model 230 may have their positions switched). [0063] In a further example, dilution model 240 may be inserted anywhere in the circuit (e.g., between patient model 210 and dialyzer model 220; between dialyzer model 220 and ECCO2RD model 230; between ECCO 2 RD model 230 and patient model 210; and/or the like). In various embodiments, dilution model 240 may operate or be represented by the increase of the blood flow rate (Q) in between models 210, 220, and/or 230, and the appropriate decrease in concentrations. [0064] In some embodiments, models 210, 220, 230, and/or 240 of acid-base model 205 may exchange values of Q and of concentrations of H + , HCO − and CO 2 (e.g., C H , C HCO , C CO ). These input/outputs may be represented as boundary conditions within acid-base model 205. [0065] Concentrations of hydrogen and CO2 may be expressed in terms of the solution pH and carbon dioxide partial pressure (pCO 2 ). In some embodiments, the conversion between pH and may be made according to the following Equation (1): ^^^^ ^^^^ = − ^^^^ ^^^^ ^^^^ 10 ^^^^ ^^^^ + embodiments. In general, FIG.3 depicts input 310 and output 320 from a patient (or patient model) 210 perspective. In various embodiments, input 310 may be or may include Q in , C H , c CO2 , c HCO3 . In some embodiments, output 320 may be or may include Q out , C H , c CO2 , c HCO3 . [0067] In some embodiments, the homeostatic dynamics of the HCO 3 − /CO 2 acid-base system may be the same or similar to dynamics described in Cherif et al. In some embodiments, equations that describe the homeostatic dynamics of the HCO 3 − /CO 2 acid-base system may include the following Equations (2) – (4): ^ ^^^ ^^^^ ^^^^ = ^^^^ ^^ − ^^^^ ^^^^ + R (2) ^^^^ ^^^^ ^^ ^^^^ ( 3) (4), where R = KH,HCO 3 CHC HCO 3 +KCO 2 CCO 2 comes from the Henderson-Hasselbach mass- actio [006 n a blood circuit with a dialyzer (for instance, dialyzer model 220), changes in total blood volume (for example, through fluid removal) may be included. From a patient perspective, such as via patient model 210 such a change may be given by the following Equation (5): (5) [0069] In various embodiments, Equati may be re-written in terms of the number of moles, such as, for each concentration C i = n i /V, the derivative for n i , may be given by the following Equation (6): ^^^^ ^^^^ = ^^^^ ^^^^ ^^^^ + ^^^^ ^^^^ ^^^^ ^^^^ = ^^^^ ^^^^ ^^^^ + ^^^^ ^^^^ ^^^^ ^^^^ ^^^^, ^^^^ ^^^^ − ^^^^ ^^^^ ^^^^ ^^^^ ^^^^ ^^^^ [0070] In some embodiments, where it is included, the concentration of the incoming flow rate (C i,in ) may not be the same as the concentration of the outgoing flow; accordingly, the acid-base model 205 may includ h f ll i E i 7 E i 7 7d [0071] In various embodiments in which each concentration may be calculated as C i = n i /V , if Qin = Qout = 0, then V ˙ = 0 ⇒ V = constant, then Equations (7a)-(7d) can be rewritten back to their original form of Equations (2)-(4). [0072] FIG.4 illustrates an example Table 400 in accordance with the present disclosure. More specifically, Table 400 represents the values and units of the constants and variables in Cherif et al., as well as their respective values in the SI units. When applicable, a conversion factor is also present. The values in Table 400 were adjusted, for instance, to cause the point pH = 7.4, pCO 2 = 40 mmHg, and C HCO3 = 2.4mM to be fixed points of Equations (7a)-7(d). [0073] In some embodiments, dialyzer model 220 may be configured the same or very similar to a parallel plates heat exchanger or model thereof. In various embodiments, the equations for a dialyzer model may be the same as, similar to, based on, and/or adapted from equations described in the Acid-Base Dynamics Disclosure. [0074] In various embodiments, dialyzer model 220 may be configured to simulate intradialytic dynamics associated with pH and acid-base homeostasis. In some embodiments, dialyzer model 220 may be configured to model an HD dialyzer. In exemplary embodiments, dialyzer model 220 may be configured for quantitating the intradialytic dynamics of HCO3– and H+, which may be parameterize to model anuric patients receiving HD. In some embodiments, dialyzer model 220 may include, may be the same or substantially the same as, or may be an adaptation of the dialyzer model described in Maheshwari et al., “An In Silico Method to Predict Net Calcium Transfer During Hemodialysis,” 201739th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), pp.2740-2743 (2017). [0075] In various embodiments, acid-base model 205 and/or models thereof, may be the same or similar to a physiological acid-base model, for example, as described in the Acid-Base Dynamics Disclosure. [0076] For example, a physiological acid-base model according to some embodiments may operate to model the effect of systemic acid-base homeostasis. In some embodiments, a physiological acid-base model may be implemented using a system of coupled nonlinear ordinary differential equations, for example, to describe the acid-base buffering kinetics through the HCO 3 -CO 2 system and incorporating the relevant physiological regulatory mechanisms. In some embodiments, physiological acid-base models may focus on HCO3-CO2 buffering kinetics, which is the most effective buffer system that controls systemic pH. [0077] The pH of extracellular fluid is mainly regulated by the following three mechanisms, which act on different timescales: (i) chemical acid-base buffering, (ii) respiratory control, and (iii) renal filtration. The following Equation (8) provides an illustrative and non-restrictive HCO 3 buffering system according to some embodiments: (8) In some embodiments, in Equation (8), it may be assumed that carbonic anhydrase (CA) accelerates the carbonic acid reactions. [0078] Chemical acid-base buffering prevents excessive changes in pH, where the timescale of this process is usually in seconds. The ability of the lung to increase or decrease ventilation allows it to regulate CO 2 removal as a gas in the expired air from the extracellular fluid, thereby adjusting the pH. In particular, due to continuous production of CO2 as a by-product of cellular metabolism, the ventilation rate must be able to accommodate alterations in CO2 in order to equilibrate the pH of the extracellular fluid. Although the process is fast, for example, occurring in a matter of minutes, it is less effective than chemical buffering. If the acid-base imbalance persists, then kidneys excrete either excess acid or base, an adaptation process that takes hours to days. [0079] The kidneys, representing a very powerful regulatory system, have the ability to secrete large amounts of H + into the tubular lumen during metabolic acidosis. Also, excretion and reabsorption of HCO3 take place in the proximal tubule and distal tubule. In particular, in the proximal tubule, H + is secreted through a Na + /H + countertransport-facilitated process, while HCO 3 is reabsorbed by combining with H + to form carbonic acid(H 2 CO 3 ) which is converted into CO2 and H2O (via carbonic anhydrase enzymatic activity). In the intercalated cells, H + /Cl- cotransporter facilitates the secretion of H + . Although this process only accounts for 5% of the secreted H + , it provides a gradient for the further secretion of more H + into the tubule lumen. Most of the secreted H + is used to reclaim the filtered HCO3, and this rate of HCO3 reabsorption is related to the rate of acid excretion. [0080] To model the acid-base homeostatic process, physiological acid-base models according to some embodiments may track the concentrations of bicarbonate ( ^ ^^^ ^^^^ ^^^^3 − ), carbon dioxide 2 and free hydrogen protons ( ^ ^^^ +). Using Henderson-Hasselbalch mass-action kinetics with renal and pulmonary regulatory mechanisms, the homeostatic dynamics of the HCO 3 /CO 2 acid-base system, as modelled by physiological acid-base models according to some embodiments may be configured according to the following Equations (9)-(11) (i.e., the physiological acid-base model equations or bicarbonate buffer kinetic system): [0081] In some embodiments, in Equations ( the following initial conditions: ^^^^ (0) ^^^^ ^^^^ (0) ^^^^ d ^^^^ (0) ^^^^ may set a patient (or virtual patient) to a normal physiological state. The parameter P H + represents the cellular production of , denotes loss either due to renal clearance and/or non-bicarbonate buffering (for example, buffering with albumin, Ca 2+ , PO^-). The hydration and de-hydration reaction rates are given by the parameters , respectively, where the values may be adjusted to reflect the carbonic anhydrase activity. In addition denotes H therapy and/or supplementation, represents the acid secretion rate f s ^e body or cellular (mitochondrial) production of CO2, and D HC0 - represents the renal fdtration rate of The effective ventilation rate (D C02 V Q ) captures the pulmonary removal of CO2, where Vo is the minute volume ventilation, and D C02 is the ventilation rate. In some embodiments, the kinetics of H2O is not included in Equations (9)-(l 1) because H2O is assumed to be abundant as a solvent. In some embodiments, Equations (9)-( 11) may be formulated via general expressions (see, for example, Equations (12) and (13)), for instance, with non-linear ventilation
[0082] In various embodiments, the first two terms in Equation (9) may account for production of P H +(t) of H + from the body, either through consumption and/or through other processes (for example, cellular metabolism), and for removal of H + as either a titratable acid and/or ammonium or non-carbonate buffering (for instance, buffering with phosphate and/or calcium). These two terms collectively correspond to H + mobilization due to buffering with non- bicarbonate buffers and other processes. The third and last terms in Equation (9) correspond to buffering reaction kinetics.
[0083] For Equation (10), the first term, represents therapy, the second and third terms, may describe renal filtration processes, where it is assumed that the amount of HCO3 lost to kidney from the blood through filtration may be related to filtered load, and the equivalence of reabsorption and acid excretion is through the splitting of CO2 by intracellular carbonic anhydrase enzymatic activity. For example, an increase in CO2 concentration may increase the conversion of CO2 into H + and 3 in a normally functioning kidney, which may result in higher acid secretion into the urine and CO2 absorption into the bloodstream. Accordingly, the second term of Equation (3) may describe acid secretion, whereas the third term is th load to be filtered. [0084] Similarly, Equation (11) may be directed to buffering kinetics, where the first term may be production of CO2 in the body (for example, a cellular metabolic or mitochondrial process), and the second term, , is the removal of CO2 through respiratory ventilation by the lung, which may depend on blood volume, cardiac output, arteriovenous difference (for example, the concentration difference between arterial and venous blood) of CO2. To simplify the model, some embodiments may assume that ventilation rate V o is constant; however, E o may depend on pCCh, oxygen partial pressure (pCh), and/or pH.
[0085] In some embodiments, the term 2 2 , may be nonlinear and may represent the effective ventilation rate, 2 > where becomes the ventilation rate and V denotes the minute volume ventilation. The empirical ventilation function may take the form of one of the following Equations (12) or (13): where Vo and Vp represent baseline and slope parameters, and Ip is the cutoff threshold. In the ventilation term, the role of pH sensors in the medulla oblongata in regulating bicarbonate buffering homeostasis is disregarded as de minimus. The ventilation term above relates the interaction between partial pressure pCO2 and pO2 and their effects on ventilation. In uses an oxygen saturation function, So 2 which may be expressed in terms of are related by In some embodiments, we may also assume that pO2 is constant. At normal pO 2 , ventilation increases by 2.5 L/min for every 1 mm Hg increase of pCO 2 . From the above Equations (5) and (6) and FIG 2, lowering pO 2 increases ventilation for a given pCO 2 and the steepness of the net effective slope of ventilation. The term Vo can be replaced by either of the functional expressions above under the assumption that pO 2 is constant or exogenously provided. In some embodiments, the simplified version may be used to determine acid-base dynamics according to some embodiments, thereby reducing the number of model parameters that need to be identified. In other embodiments, all the effects of ventilation may be lumped into ventilatory rate parameters. In some embodiments, Equations (12) and (13) may provide general physiological acid-base models (for example, with Equations (9)-(11) being a specific implementation thereof). [0086] In various embodiments, acid-base model 205 and/or models thereof, may be the same or similar to an intradialytic acid-base model, for example, as described in the Acid-Base Dynamics Disclosure.
[0087] FIG. 5 illustrates a block diagram of an intradialytic acid-base model according to some embodiments. As shown in FIG. 5, an intradialytic acid-base model 502 may include a dialysis patient model or compartment 510 and a dialyzer model or compartment 512 having a dialyzer 514. In various embodiments, patient compartment 510 may include the distribution volume (V ex ) and concentrations of acid-base variables (for example, ). In some embodiments, Q p , Qd, and Q uf are the plasma flow rate, dialysate flow rate, and ultrafiltration rate, respectively. A single dialyzer fiber 516 may depict the counter-current flows interaction between blood and dialysate flows in dialyzer 514, in which there is an intradialytic transfer between blood and dialysate through the infinitesimal fiber segment Ax. The subscripts for q and C denote plasma (pl) or extracellular fluid (ex), dialyzer input (di) and output (do), dialyzer blood inlet (pi), and outlet (out).
[0088] For the dialysis patient model, some embodiments may use a physiologically-based dynamic model describing the regulation of buffering system with Henderson- Hasselbalch mass-action kinetics, in which we incorporate the endogenous production of both CO2 and H + , non-bicarbonate buffering, and the physiologic regulation of the buffering system through ventilation and renal excretion. Although there are other buffer systems, some embodiments may focus mainly on the buffering system, which is the most abundance and effective buffer system in the body. Under physiologic conditions, pH is regulated mainly by chemical acid-base buffering, respiratory control and renal glomerular filtration. The chemical acid-base buffering is modeled using Henderson-Hasselbalch mass- action kinetics, in which, in some embodiments, the action of carbonic anhydrase is fast, and association and disassociate rates, respectively.
[0089] During normal cellular metabolic activities, CO2 and H + are produced as by-products; this endogenous production is captured by , respectively. In addition, acid is mobilized or buffered with other non-bicarbonate buffers such as phosphate, and this is represented by is the mobilization or removal rate. In normal physiologic conditions, the kidney, the second regulatory organ, is responsible for regulating H CCU level by either secreting excess H+ into tubular lumen during metabolic acidosis, and/or excreting or reabsorbing HCOr in the proximal and distal segments of nephrons. Unlike the physiological acid-base model according to some embodiments, the dialysis patient model may not include (or may include lower functioning of) the renal regulation of and H + because the dialysis patient model is adapted to dialysis patients where the ability of excrete acid is impaired. As a result, in some embodiments, the dialysis patient model may assume that there is little or even no renal function. In some embodiments, the renal function is replaced by the use of the dialyzer (for example, by dialyzer fluxes) implemented via the dialyzer model. For example, the expression may account for bicarb flux from patient and post-dialyzer flux to patient, respectively, which is characterized by blood flow rate,
Q p , and ultrafiltrate rate, Q uf . For the third regulatory mechanism, in order to regulate CO2 removal, lung increases or decreases ventilation, which is triggered by the response of central and peripheral chemoreceptors to changes in describes removal of CO2 through ventilation by lung characterized by blood volume, cardiac output, arteriovenous difference of CO2. The parameters D C02 and V o are ventilation rate and minute ventilation, respectively. In some embodiments, the dialysis patient model describing acid-base homeostasis may be according to the following Equations ( 14)-(16):
[0090] During dialysis, a patient loses a significant amount of fluid. This is assumed to occur at a constant ultrafiltration rate, Q uf , and fluid removal by ultrafiltration occurs in proportion to the compartmental distribution volume. Accordingly, in some embodiments, extracellular fluid volume may be determined according to the following Equation (17):
[0091] In some embodiments, the dialyzer model may include two spatial temporal models describing both blood and dialysate sides using hyperbolic partial differential equations. The concentration of HCCh in the blood side may be determined according to the following
Equation (18): where N is the number of fibers, A denotes fiber cross-sectional area, represents the reflection coefficient - defines the Peclet, L is the fiber length, K 0 A is the effective membrane mass-transfer coefficient for H is Gibbs-Donnan correction factor
(i.e., a Gibbs-Donnan corrected dialyzer model). In some embodiments, β and may be set to a constant value of 1.05 which corresponds to 5% of HCO v. In various embodiments, for the dialyzer side, it may be assumed that the dialysate flow is uniform and equally shared by the N fibers present in the dialyzer housing. The hyperbolic partial differential equation describing the concentration within annulus dialysate flow boundary may have the form of the following Equation (19): where A d circular cross-sectional area of annulus space for dialysate flow around a fiber. In various embodiments, plasma flow rate (Q p ) may decrease along the fiber length in the dialyzer due to ultrafiltration, the decrease may be linearly along the fiber length. In some embodiments, the dialysate flow rate may increase by the amount of fluid removed by ultrafiltration from the blood side to the dialysate side, resulting in counter-current kinetics. The spatial aspects of plasma and dialysate flow rates may be determined according to the following equations (20) and (21): where Q pi and Q di are initial plasma and dialysate flow rates.
[0092] In some embodiments, Equations (14)-(21) may constitute components of intradialytic acid-base models describing, simulating, predicting, or otherwise processing intradialytic acidbase dynamics. [0093] FIG.6 illustrates table 600 of a first illustrative set of acid-base model parameters and values according to some embodiments. FIG.7 illustrates table 700 of a second illustrative set of acid-base model parameters and values according to some embodiments. In some embodiments, for the dialyzer models, the difference between the values of the hydration and dehydration rates (K H and K CO2 ), may differ from one order of magnitude between those described in the Acid-Base Dynamics Disclosure. [0094] FIG.8A and 8B illustrates results of a simulation of a physiological acid-base model according to some embodiments. More specifically, FIG.8A depicts results 810 from acid-base models from the Acid-Base Dynamics Disclosure. FIG.8B depicts results 820 from acid-base models according to some embodiments, for instance, acid-base model 205. The results in FIGS. 8A and 8B are obtained with only one ρ residual being zero, for instance, when no residual kidney function is included in the model. [0095] FIG.9 depicts an illustrative ECCO 2 RD model 230 according to some embodiments. In some embodiments, ECCO 2 RD model 230 may be approximately modeled as a one-dimensional (1D) diffusion device between blood and air. In various embodiments, ECCO2RD model 230 may receive input 910 to generate output 920. As shown in FIG.9, the dynamics of ECCO2RD model 230 may include the following Equations (22a) and (22b): [0096] In some embodiments, the equations that describe the dynamics of ECCO 2 RD model 230, such as Equations (22a) and (22b), may be the same, similar to, and/or adapted from the equations described in Habran et al., “Mathematical Modeling of Extracorporeal CO2 Removal Therapy,” Medical & Biological Engineering & Computing, 56(3):421–434 (2018) (“Habran et al.”), the contents of which are incorporated by reference as if fully set forth in the present disclosure. FIG.10 depicts Table 1000 providing illustrative values of parameters used in Equations (22a) and (22b): [0097] FIGS.11A-11D depict results for an acid-base model according to some embodiments. More specifically, FIG.11A depicts results for a patient with normal renal function; FIG.11B depicts results for a patient without renal function; FIG.11C depicts results for calculated PCO 2 as a function of time over measured points; and FIG.11D depicts results for time evolution of a model with validation data. [0098] In some embodiments, the results for the ECCO 2 RD model may be determined for Q d = 268 mL/min, the parameters for the patient may be Qp = 350 mL/min, V0 = 0 L/min and the patient has full residual kidney function (ρresidual = 1). The optimized parameters may be PH = 1.6962 × 10 −7 mol/L/min, P CO2 = 3.8074 × 10−6 mol/L/min, D CO2 = 9.9278 × 10 −5 L−1, γ = 2.9864 × 103 min−1, V blood = 0.5294 L, and V gas = 2.9694 L. EXAMPLE RESULTS: VALIDATION OF ACID-BASE MODELS [0099] An average of data from 10 ventilated critically ill patients with ARDS and AKI undergoing renal and respiratory replacement therapy according to Forster et al., “Low-flow CO2 removal integrated into a renal-replacement circuit can reduce acidosis and decrease vasopressor requirements,” Crit Care 17, R154 (2013). https://doi.org/10.1186/cc12833 was used to validate the model. The following parameter values were set at the beginning of treatment: blood flow 355±79 ml/min, dialysate flow 2.3±7 L/h, ultrafiltration rate 58±66 ml/h, gas flow 5.2±1.0 L/min, treatment time 95±68 h. The model was parameterized to obtain patient-specific parameters. FIG.12 depicts results for an acid-base model according to some embodiments, in which the results are depicted for model prediction for a patient with renal failure and low ventilation. In FIG.12, the solid line is the model prediction; the dots are patient data. As indicated by FIG.12, models according to some embodiments are capable of accurately predicting the pH, and pCO2 over time for example, for at least the first 24 hours of treatment. [0100] FIG.13 illustrates an embodiment of an exemplary computing architecture 1300 suitable for implementing various embodiments as previously described. In various embodiments, the computing architecture 1300 may comprise or be implemented as part of an electronic device. The embodiments are not limited in this context. \ [0101] As used in this application, the terms “system” and “component” and “module” are intended to refer to a computer-related entity, either hardware, a combination of hardware and software, software, or software in execution, examples of which are provided by the exemplary computing architecture 1300. For example, a component can be, but is not limited to being, a process running on a processor, a processor, a hard disk drive, multiple storage drives (of optical and/or magnetic storage medium), an object, an executable, a thread of execution, a program, and/or a computer. By way of illustration, both an application running on a server and the server can be a component. One or more components can reside within a process and/or thread of execution, and a component can be localized on one computer and/or distributed between two or more computers. [0102] Further, components may be communicatively coupled to each other by various types of communications media to coordinate operations. The coordination may involve the uni- directional or bi-directional exchange of information. For instance, the components may communicate information in the form of signals communicated over the communications media. The information can be implemented as signals allocated to various signal lines. In such allocations, each message is a signal. Further embodiments, however, may alternatively employ data messages. Such data messages may be sent across various connections. Exemplary connections include parallel interfaces, serial interfaces, and bus interfaces. [0103] The computing architecture 1300 includes various common computing elements, such as one or more processors, multi-core processors, co-processors, memory units, chipsets, controllers, peripherals, interfaces, oscillators, timing devices, video cards, audio cards, multimedia input/output (I/O) components, power supplies, and so forth. The embodiments, however, are not limited to implementation by the computing architecture 1300. [0104] As shown in FIG.13, the computing architecture 1300 comprises a processing unit 1304, a system memory 1306 and a system bus 1308. The processing unit 1304 may be a commercially available processor and may include dual microprocessors, multi-core processors, and other multi-processor architectures. [0105] The system bus 1308 provides an interface for system components including, but not limited to, the system memory 1306 to the processing unit 1304. The system bus 1308 can be any of several types of bus structure that may further interconnect to a memory bus (with or without a memory controller), a peripheral bus, and a local bus using any of a variety of commercially available bus architectures. Interface adapters may connect to the system bus 1308 via a slot architecture. Example slot architectures may include without limitation Accelerated Graphics Port (AGP), Card Bus, (Extended) Industry Standard Architecture ((E)ISA), Micro Channel Architecture (MCA), NuBus, Peripheral Component Interconnect (Extended) (PCI(X)), PCI Express, Personal Computer Memory Card International Association (PCMCIA), and the like. [0106] The system memory 1306 may include various types of computer-readable storage media in the form of one or more higher speed memory units, such as read-only memory (ROM), random-access memory (RAM), dynamic RAM (DRAM), Double-Data-Rate DRAM (DDRAM), synchronous DRAM (SDRAM), static RAM (SRAM), programmable ROM (PROM), erasable programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), flash memory, polymer memory such as ferroelectric polymer memory, ovonic memory, phase change or ferroelectric memory, silicon-oxide-nitride-oxide-silicon (SONOS) memory, magnetic or optical cards, an array of devices such as Redundant Array of Independent Disks (RAID) drives, solid state memory devices (for example, USB memory, solid state drives (SSD) and any other type of storage media suitable for storing information. In the illustrated embodiment shown in FIG.13, the system memory 1306 can include non-volatile memory 1310 and/or volatile memory 1312. A basic input/output system (BIOS) can be stored in the non- volatile memory 1310. The computer 1302 may include various types of computer-readable storage media in the form of one or more lower speed memory units, including an internal (or external) hard disk drive (HDD) 1314, a magnetic floppy disk drive (FDD) 1316 to read from or write to a removable magnetic disk 1311, and an optical disk drive 1320 to read from or write to a removable optical disk 1322 (for example, a CD-ROM or DVD). The HDD 1314, FDD 1316 and optical disk drive 1320 can be connected to the system bus 1308 by a HDD interface 1324, an FDD interface 1326 and an optical drive interface 1328, respectively. The HDD interface 1324 for external drive implementations can include at least one or both of Universal Serial Bus (USB) and IEEE 1114 interface technologies. [0107] The drives and associated computer-readable media provide volatile and/or nonvolatile storage of data, data structures, computer-executable instructions, and so forth. For example, a number of program modules can be stored in the drives and memory units 1310, 1312, including an operating system 1330, one or more application programs 1332, other program modules 1334, and program data 1336. In one embodiment, the one or more application programs 1332, other program modules 1334, and program data 1336 can include, for example, the various applications and/or components of a computing device. [0108] A user can enter commands and information into the computer 1302 through one or more wired/wireless input devices, for example, a keyboard 1338 and a pointing device, such as a mouse 1340. These and other input devices are often connected to the processing unit 1304 through an input device interface 1342 that is coupled to the system bus 1308, but can be connected by other interfaces. [0109] A monitor 1344 or other type of display device is also connected to the system bus 1308 via an interface, such as a video adaptor 1346. The monitor 1344 may be internal or external to the computer 1302. In addition to the monitor 1344, a computer typically includes other peripheral output devices, such as speakers, printers, and so forth. [0110] The computer 1302 may operate in a networked environment using logical connections via wired and/or wireless communications to one or more remote computers, such as a remote computer 1348. The remote computer 1348 can be a workstation, a server computer, a router, a personal computer, portable computer, microprocessor-based entertainment appliance, a peer device or other common network node, and typically includes many or all of the elements described relative to the computer 1302, although, for purposes of brevity, only a memory/storage device 1350 is illustrated. The logical connections depicted include wired/wireless connectivity to a local area network (LAN) 1352 and/or larger networks, for example, a wide area network (WAN) 1354. Such LAN and WAN networking environments are commonplace in offices and companies, and facilitate enterprise-wide computer networks, such as intranets, all of which may connect to a global communications network, for example, the Internet. [0111] The computer 1302 is operable to communicate with wired and wireless devices or entities using the IEEE 802 family of standards, such as wireless devices operatively disposed in wireless communication (for example, IEEE 802.16 over-the-air modulation techniques). This includes at least Wi-Fi (or Wireless Fidelity), WiMax, and Bluetooth™ wireless technologies, among others. Thus, the communication can be a predefined structure as with a conventional network or simply an ad hoc communication between at least two devices. Wi-Fi networks use radio technologies called IEEE 802.11x (a, b, g, n, etc.) to provide secure, reliable, fast wireless connectivity. A Wi-Fi network can be used to connect computers to each other, to the Internet, and to wire networks (which use IEEE 802.3-related media and functions). [0112] Numerous specific details have been set forth herein to provide a thorough understanding of the embodiments. It will be understood by those skilled in the art, however, that the embodiments may be practiced without these specific details. In other instances, well-known operations, components, and circuits have not been described in detail so as not to obscure the embodiments. It can be appreciated that the specific structural and functional details disclosed herein may be representative and do not necessarily limit the scope of the embodiments. [0113] Some embodiments may be described using the expression “coupled” and “connected” along with their derivatives. These terms are not intended as synonyms for each other. For example, some embodiments may be described using the terms “connected” and/or “coupled” to indicate that two or more elements are in direct physical or electrical contact with each other. The term “coupled,” however, may also mean that two or more elements are not in direct contact with each other, but yet still co-operate or interact with each other. [0114] Unless specifically stated otherwise, it may be appreciated that terms such as “processing,” “computing,” “calculating,” “determining,” or the like, refer to the action and/or processes of a computer or computing system, or similar electronic computing device, that manipulates and/or transforms data represented as physical quantities (for example, electronic) within the computing system’s registers and/or memories into other data similarly represented as physical quantities within the computing system’s memories, registers or other such information storage, transmission or display devices. The embodiments are not limited in this context. [0115] It should be noted that the methods described herein do not have to be executed in the order described, or in any particular order. Moreover, various activities described with respect to the methods identified herein can be executed in serial or parallel fashion. [0116] Although specific embodiments have been illustrated and described herein, it should be appreciated that any arrangement calculated to achieve the same purpose may be substituted for the specific embodiments shown. This disclosure is intended to cover any and all adaptations or variations of various embodiments. It is to be understood that the above description has been made in an illustrative fashion, and not a restrictive one. Combinations of the above embodiments, and other embodiments not specifically described herein will be apparent to those of skill in the art upon reviewing the above description. Thus, the scope of various embodiments includes any other applications in which the above compositions, structures, and methods are used. [0117] Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims. [0118] As used herein, an element or operation recited in the singular and proceeded with the word “a” or “an” should be understood as not excluding plural elements or operations, unless such exclusion is explicitly recited. Furthermore, references to “one embodiment” of the present disclosure are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. [0119] The present disclosure is not to be limited in scope by the specific embodiments described herein. Indeed, other various embodiments of and modifications to the present disclosure, in addition to those described herein, will be apparent to those of ordinary skill in the art from the foregoing description and accompanying drawings. Thus, such other embodiments and modifications are intended to fall within the scope of the present disclosure. Furthermore, although the present disclosure has been described herein in the context of a particular implementation in a particular environment for a particular purpose, those of ordinary skill in the art will recognize that its usefulness is not limited thereto and that the present disclosure may be beneficially implemented in any number of environments for any number of purposes. Accordingly, the claims set forth below should be construed in view of the full breadth and spirit of the present disclosure as described herein.