Two simple mathematical models of advection and diffusion of hydrogen within the retina are discussed. The work is motivated by the hydrogen clearance technique, which is used to estimate blood flow in the retina. The first model assumes that the retina consists of three, well-mixed layers with different thickness, and the second is a two-dimensional model consisting of three regions that represent the layers in the retina. Diffusion between the layers and leakage through the outer edges are considered. Solutions to the governing equations are obtained by employing Fourier series and finite difference methods for the two models, respectively. The effect of important parameters on the hydrogen concentration is examined and discussed. The results contribute to understanding the dispersal of hydrogen in the retina and in particular the effect of flow in the vascular retina. It is shown that the predominant features of the process are captured by the simpler model.

Arterial diseases such as aneurysm and stenosis may result from the mechanical and/or morphological change of an arterial wall structure and correspondingly altered hemodynamics. The development of a 3D computational model of blood flow can be useful to study the hemodynamics in major blood vessels and may provide an insight into the noninvasive technique to detect arterial diseases in early stage. In this paper, we present a three-dimensional model of blood flow in the aorta, which is based on the immersed boundary method to describe the interaction of blood flow with the aortic wall. Our simulation results show that the hysteresis loop is evident in the pressure-diameter relationship of the normal aorta when the arterial wall is considered to be viscoelastic. In addition, it is shown that flow patterns and pressure distributions are altered in response to the change of aortic morphology.

We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by a finite difference scheme and solved iteratively using the Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Then, the algebraic multigrid method is applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed, which reduces computational complexity from O(N2) to O(NlogN), where N is the number of grid points. Integrals involving the Dirac delta function are evaluated directly by coordinate transformation, which yields more accurate results compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for lithiumion (Li-ion) batteries are shown to be in good agreement with the experimental data and the results from previous studies.

Experimental observations show that a strong magnetic field has a dramatic influence on the sedimentation of RBCs, which motivates us to model the sedimentation of red blood cell (RBC) under strong external magnetic body force. To model the sedimentation of a RBC in a square duct and a circular pipe, a recently developed technique derived from the lattice Boltzmann and the distributed Lagrange multiplier/fictitious domain methods (LBM-DLM/FD) is extended to employ the mesoscopic network model for simulations of the sedimentation of a RBC in flow. The flow is simulated by the LBM with a strong magnetic body force, while the network model is used for modeling RBC deformation. The fluid-RBC interactions are enforced by the Lagrange multiplier. The sedimentation of RBC in a square duct and a circular pipe is simulated, which demonstrates the developed method's capability to model the sedimentation of RBCs in various flows. Numerical results illustrate that the terminal settling velocity increases incrementally with the exerted body force. The deformation of RBC has a significant effect on the terminal settling velocity due to the change in the frontal area. The larger the exerted force, the smaller the frontal area and the larger the RBC deformation become. Additionally, the wall effect on the motion and deformation of RBC is also investigated.

Stent placement is considered as a promising and minimally invasive technique to prevent rupture of aneurysm and favor coagulation mechanism inside the aneurysm. Many scholars study the effect of the stent on blood flow in cerebral aneurysm by numerical simulations, and usually regard blood as the Newtonian fluid, blood, however, is a kind of non-Newtonian fluid in practice. The main purpose of the present paper is to investigate the effect of non-Newtonian behavior on the hemodynamic characteristics of blood flow in stented cerebral aneurysm with lattice Boltzmann method. The Casson model is used to describe the blood non-Newtonian character, which is one of the most popular models in depicting blood fluid. In particular, hemodynamic characteristics derived with Newtonian and non-Newtonian models are studied, and compared in detail. The results show that the non-Newtonian effect gives a great influence on hemodynamic characteristics of blood flow in stented cerebral aneurysm, especially in small necked ones.