No Arabic abstract
Complex fluids flow in complex ways in complex structures. Transport of water and various organic and inorganic molecules in the central nervous system are important in a wide range of biological and medical processes [C. Nicholson, and S. Hrabv{e}tova, Biophysical Journal, 113(10), 2133(2017)]. However, the exact driving mechanisms are often not known. In this paper, we investigate flows induced by action potentials in an optic nerve as a prototype of the central nervous system (CNS). Different from traditional fluid dynamics problems, flows in biological tissues such as the CNS are coupled with ion transport. It is driven by osmosis created by concentration gradient of ionic solutions, which in term influence the transport of ions. Our mathematical model is based on the known structural and biophysical properties of the experimental system used by the Harvard group Orkand et al [R.K. Orkand, J.G. Nicholls, S.W. Kuffler, Journal of Neurophysiology, 29(4), 788(1966)]. Asymptotic analysis and numerical computation show the significant role of water in convective ion transport. The full model (including water) and the electrodiffusion model (excluding water) are compared in detail to reveal an interesting interplay between water and ion transport. In the full model, convection due to water flow dominates inside the glial domain. This water flow in the glia contributes significantly to the spatial buffering of potassium in the extracellular space. Convection in the extracellular domain does not contribute significantly to spatial buffering. Electrodiffusion is the dominant mechanism for flows confined to the extracellular domain.
It is suggested that the propagation of the action potential is accompanied by an axoplasmic pressure pulse propagating in the axoplasm along the axon length. The pressure pulse stretch-modulates voltage-gated Na (Nav) channels embedded in the axon membrane, causing their accelerated activation and inactivation and increasing peak channel conductance. As a result, the action potential propagates due to mechano-electrical activation of Nav channels by straggling ionic currents and the axoplasmic pressure pulse. The velocity of such propagation is higher than in the classical purely electrical Nav activation mechanism, and it may be close to the velocity of propagation of pressure pulses in the axoplasm. Extracellular Ca ions influxing during the voltage spike, or Ca ions released from intracellular stores, may trigger a mechanism that generates and augments the pressure pulse, thus opposing its viscous decay. The model can potentially explain a number of phenomena that are not contained within the purely electrical Hodgkin-Huxley-type framework: the Meyer-Overton rule for the effectiveness of anesthetics, as well as various mechanical, optical and thermodynamic phenomena accompanying the action potential. It is shown that the velocity of propagation of axoplasmic pressure pulses is close to the measured velocity of the nerve impulse, both in absolute magnitude and in dependence on axon diameter, degree of myelination and temperature.
The accumulation of potassium in the narrow space outside nerve cells is a classical subject of biophysics that has received much attention recently. It may be involved in potassium accumulation textcolor{black}{including} spreading depression, perhaps migraine and some kinds of epilepsy, even (speculatively) learning. Quantitative analysis is likely to help evaluate the role of potassium clearance from the extracellular space after a train of action potentials. Clearance involves three structures that extend down the length of the nerve: glia, extracellular space, and axon and so need to be described as systems distributed in space in the tradition used for electrical potential in the `cable equations of nerve since the work of Hodgkin in 1937. A three-compartment model is proposed here for the optic nerve and is used to study the accumulation of potassium and its clearance. The model allows the convection, diffusion, and electrical migration of water and ions. We depend on the data of Orkand et al to ensure the relevance of our model and align its parameters with the anatomy and properties of membranes, channels, and transporters: our model fits their experimental data quite well. The aligned model shows that glia has an important role in buffering potassium, as expected. The model shows that potassium is cleared mostly by convective flow through the syncytia of glia driven by osmotic pressure differences. A simplified model might be possible, but it must involve flow down the length of the optic nerve. It is easy for compartment models to neglect this flow. Our model can be used for structures quite different from the optic nerve that might have different distributions of channels and transporters in its three compartments. It can be generalized to include a fourth (distributed) compartment representing blood vessels to deal with the glymphatic flow into the circulatory system.
Despite their importance in many biological, ecological and physical processes, microorganismal fluid flows under tight confinement have not been investigated experimentally. Strong screening of Stokelets in this geometry suggests that the flow fields of different microorganisms should be universally dominated by the 2D source dipole from the swimmers finite-size body. Confinement therefore is poised to collapse differences across microorganisms, that are instead well-established in bulk. Here we combine experiments and theoretical modelling to show that, in general, this is not correct. Our results demonstrate that potentially minute details like microswimmers spinning and the physical arrangement of the propulsion appendages have in fact a leading role in setting qualitative topological properties of the hydrodynamic flow fields of micro-swimmers under confinement. This is well captured by an effective 2D model, even under relatively weak confinement. These results imply that active confined hydrodynamics is much richer than in bulk, and depends in a subtle manner on size, shape and propulsion mechanisms of the active components.
The most essential characteristic of any fluid is the velocity field v(r) and this is particularly true for macroscopic quantum fluids. Although rapid advances have occurred in quantum fluid v(r) imaging, the velocity field of a charged superfluid - a superconductor - has never been visualized. Here we use superconductive-tip scanning tunneling microscopy to image the electron-pair density r{ho}_S(r) and velocity v_S(r) fields of the flowing electron-pair fluid in superconducting NbSe2. Imaging v_S(r) surrounding a quantized vortex finds speeds reaching 10,000 km/hr. Together with independent imaging of r{ho}_S(r) via Josephson tunneling, we visualize the supercurrent density j_S(r)=r{ho}_S(r)v_S(r), which peaks above 3 x 10^7 A/cm^2. The spatial patterns in electronic fluid flow and magneto-hydrodynamics reveal hexagonal structures co-aligned to the crystal lattice and quasiparticle bound states, as long anticipated. These novel techniques pave the way for electronic fluid flow visualization in many other quantum fluids.
Surface effects become important in microfluidic setups because the surface to volume ratio becomes large. In such setups the surface roughness is not any longer small compared to the length scale of the system and the wetting properties of the wall have an important influence on the flow. However, the knowledge about the interplay of surface roughness and hydrophobic fluid-surface interaction is still very limited because these properties cannot be decoupled easily in experiments. We investigate the problem by means of lattice Boltzmann (LB) simulations of rough microchannels with a tunable fluid-wall interaction. We introduce an ``effective no-slip plane at an intermediate position between peaks and valleys of the surface and observe how the position of the wall may change due to surface roughness and hydrophobic interactions. We find that the position of the effective wall, in the case of a Gaussian distributed roughness depends linearly on the width of the distribution. Further we are able to show that roughness creates a non-linear effect on the slip length for hydrophobic boundaries.