In this paper, we develop a mathematical model of blood circulation in the liver lobule. We aim to find the pressure and flux distributions within a liver lobule. We also investigate the effects of changes in pressure that occur following a resection of part of the liver, which often leads to high pressure in the portal vein. The liver can be divided into functional units called lobules. Each lobule has a hexagonal cross-section, and we assume that its longitudinal extent is large compared with its width. We consider an infinite lattice of identical lobules and study the two-dimensional flow in the hexagonal cross-sections. We model the sinusoidal space as a porous medium, with blood entering from the portal tracts (located at each of the vertices of the cross-section of the lobule) and exiting via the centrilobular vein (located in the center of the cross-section). We first develop and solve an idealized mathematical model, treating the porous medium as rigid and isotropic and blood as a Newtonian fluid. The pressure drop across the lobule and the flux of blood through the lobule are proportional to one another. In spite of its simplicity, the model gives insight into the real pressure and velocity distribution in the lobule. We then consider three modifications of the model that are designed to make it more realistic. In the first modification, we account for the fact that the sinusoids tend to be preferentially aligned in the direction of the centrilobular vein by considering an anisotropic porous medium. In the second, we account more accurately for the true behavior of the blood by using a shear-thinning model. We show that both these modifications have a small quantitative effect on the behavior but no qualitative effect. The motivation for the final modification is to understand what happens either after a partial resection of the liver or after an implantation of a liver of small size. In these cases, the pressure is observed to rise significantly, which could cause deformation of the tissue. We show that including the effects of tissue compliance in the model means that the total blood flow increases more than linearly as the pressure rises.

We present a mathematical model of blood and interstitial flow in the liver. The liver is treated as a lattice of hexagonal ‘classic’ lobules, which are assumed to be long enough that end effects may be neglected and a two-dimensional problem considered. Since sinusoids and lymphatic vessels are numerous and small compared to the lobule, we use a homogenized approach, describing the sinusoidal and interstitial spaces as porous media. We model plasma filtration from sinusoids to the interstitium, lymph uptake by lymphatic ducts, and lymph outflow from the liver surface. Our results show that the effect of the liver surface only penetrates a depth of a few lobules’ thickness into the tissue. Thus, we separately consider a single lobule lying sufficiently far from all external boundaries that we may regard it as being in an infinite lattice, and also a model of the region near the liver surface. The model predicts that slightly more lymph is produced by interstitial fluid flowing through the liver surface than that taken up by the lymphatic vessels in the liver and that the non-peritonealized region of the surface of the liver results in the total lymph production (uptake by lymphatics plus fluid crossing surface) being about 5 % more than if the entire surface were covered by the Glisson–peritoneal membrane. Estimates of lymph outflow through the surface of the liver are in good agreement with experimental data. We also study the effect of non-physiological values of the controlling parameters, particularly focusing on the conditions of portal hypertension and ascites. To our knowledge, this is the first attempt to model lymph production in the liver. The model provides clinically relevant information about lymph outflow pathways and predicts the systemic response to pathological variations.

During eye rotations the vitreous humour moves with respect to the eye globe. This relative motion has been suggested to possibly have an important role in inducing degradation of the gel structure, which might lead to vitreous liquefaction and/or posterior vitreous detachment. Aim of the present work is to study the characteristics of vitreous motion induced by eye rotations. We use an experimental setup, consisting of a Perspex model of the vitreous chamber that, for simplicity, is taken to have a spherical shape. The model is filled with an artificial vitreous humour, prepared as a solution of agar powder and hyaluronic acid sodium salt in deionised water, which has viscoelastic mechanical properties similar to those of the real vitreous. The model rotates about an axis passing through the centre of the sphere and velocity measurements are taken on the equatorial plane orthogonal to the axis of rotation, using an optical technique. The results show that fluid viscoelasticity has a strong influence on flow characteristics. In particular, at certain frequencies of oscillation of the eye model, fluid motion can be resonantly excited. This means that fluid velocity within the domain can be significantly larger than that of the wall. The frequencies for which resonant excitation occurs are within the range of possible eye rotations frequencies. Therefore, the present results suggest that resonant excitation of vitreous motion is likely to occur in practice. This, in turn, implies that eye rotations produce large stresses on the retina and within the vitreous that may contribute to the disruption of the vitreous gel structure. The present results also have implications for the choice of the ideal properties for vitreous substitute fluids.

We aim to understand how mechanical causation influences retinal detachment and reattachment processes. In particular, myopes suffer retinal detachment more frequently than emmetropes, and following a retinal detachment, scleral buckling promotes retinal reattachment. We test the hypothesis that stresses arising from saccadic eye rotations are involved in the processes, and that the alteration in the stress due to the change in the vitreous chamber geometry is sufficient to explain the phenomena. The vitreous chamber of the eye has an approximately spherical shape and it is filled with vitreous humor. We developed a mathematical model, treating the vitreous chamber in emmetropic and myopic eyes as a spheroid and in eyes subjected to scleral buckling as a sphere with a circumferential indentation. We assume that the eye performs prescribed small-amplitude, periodic, torsional rotations and we solve semi-analytically for the fluid pressure, velocity, and stress distributions. The shape of the vitreous chamber has a large effect on the retinal stress. The vitreous and the retina of a highly myopic eye continuously experience shear stresses significantly higher than those of an emmetropic eye. An eye fitted with a scleral buckle experiences large stress levels localized around the buckle. Our results provide a mechanical explanation for the more frequent occurrence of posterior vitreous detachment and retinal detachment in myopic eyes. To understand how the stress distribution in a buckled eye facilitates reattachment, an additional model of the details of the reattachment process should be coupled to this model.

The retina is composed of two main layers—the neuroretina and the retinal pigment epithelium (RPE)—that are separated by a potential gap termed the sub-retinal space (SRS). Accumulation of fluid in the SRS may result in a retinal detachment. A key function of the RPE is to prevent fluid accumulation in the SRS by actively pumping fluid from this space to the choroid. We have developed a mathematical model of this process that incorporates the transport of seven chemical species: Na ⁺ , K ⁺ , Cl ⁻ , HCO 3 − , H ⁺ , CO ₂ and H ₂ CO ₃ . This allows us to estimate solute and water fluxes and to understand the role of the different membrane ion channels. We have performed a global sensitivity analysis using the extended Fourier amplitude sensitivity test to investigate the relative importance of parameters in generating the model outputs. The model predicts that flow across the RPE is driven by an osmotic gradient in the cleft gap between adjacent cells. Moreover, the model estimates how water flux is modified in response to inhibition of membrane ion channels and carbonic anhydrase (CA). It provides a possible explanation for how CA inhibitors, which are used clinically to prevent fluid accumulation in the SRS, may be acting.