This model extends the original autoregulatory model presented by Feldberg et al., adding an explicit glomerulus and calculating model SNGFR in addition to afferent flow.

The source code is available under the BSD 3-Clause license.

The original model comprised only an afferent arteriole (AA) and an efferent
arteriole (EA), with the hydrostatic pressure defined at the start of the AA
(the *inflow* pressure) and the end of the EA (the *venous return* pressure).
The glomerulus was treated only as a single spatial point, at the junction of
the AA and EA segments.
This meant that once the resistance of each AA segments was known, the
pressure drop along the entire model could be directly obtained (the EA was
assumed to have a constant resistance).

We added an explicit glomerular capillary bed (GC), of unit length and constant resistance, whose membrane is permeable to water and solutes. There is a pressure gradient across the membrane which drives filtration. The single-nephron glomerular filtration rate (SNGFR) is the rate at which plasma is filtered from the GC into Bowmanâ€™s capsule.

The flow rate in the GC decreases as a function of distance from the AA, in response to this filtration. This change in flow rate affects the pressure drop along the GC, and the filtration pressure gradient.

Accordingly, once the resistance of each AA segment is known we need to
*solve* for the hydrostatic pressure at the start of the GC, such that the
pressure drop across the system equals the difference between the *inflow*
and *venous return* pressures.
Once this pressure is known, we can then calculate the flow rate and pressure
*everywhere* in the model.

In the updated model, there is an explicit pressure profile along the glomerular capillary bed and glomerular filtration is accounted for. The initial hydrostatic pressure , blood flow , and plasma flow depend upon the afferent blood flow and resistance and the plasma hematocrit :

From these glomerular inflows we can calculate the SNGFR and efferent blood flow , given the glomerular resistance , filtration coefficient , and oncotic pressure :

Oncotic pressure is defined as a cubic function of plasma protein concentration :

In order to obtain a solution to these equations, we need to find the value of
for which *conservation of energy* is preserved:

That is, for the correct value of the pressure drop across the EA as determined by the pressure profile will equal the pressure drop across the EA as determined by the efferent blood flow rate .

Given the net afferent resistance , a solution for is obtained using a bisection method. This determines the hydrostatic pressure and flow rate at every spatial location in the model. The muscle activation of each AA segment can then be updated and a new value for the inner radius of each AA segment is obtained. The net afferent resistance is updated accordingly, and can again be solved.

This iterative process begins with the normalised inner radius of each segment set to 1, and is repeated until the relative change in inner radius is less than .

The following figures demonstrate that this implementation can reproduce the results presented in the original manuscript.

**Figure 1:** The relationship between stress \(\sigma\) and
strain \(\varepsilon\), for varying levels of muscle activation \(\Psi\).
Solid lines show net stress \(\sigma_e + \sigma_a\) and dashed lines show
active stress \(\sigma_a\).

**Figure 2:** The relationship between transmural pressure
\(P\) and normalised inner radius \(r_P\), for varying levels of muscle
activation \(\Psi\).

**Figure 3:** The relationship between transmural pressure
\(P\) and normalised inner radius \(r_P\), for three myogenic responses
\(\Psi_{MYO}(P)\) with different slopes.