Search

Nonlinear phenomena including multiple equilibria and spontaneous oscillations are common in fluid networks containing either multiple phases or constituents. In many systems, such behavior might be attributed to the complicated geometry of the network, the complex rheology of the constituent fluids, or, in the case of microvascular blood flow, biological control. In this paper we investigate two examples of a simple three-node fluid network containing two miscible Newtonian fluids of differing viscosities, the first modeling microvascular blood flow and the second modeling stratified laminar flow. We use a combination of analytic and numerical techniques to identify and track saddle-node and Hopf bifurcations through the large parameter space. In both models, we document sustained spontaneous oscillations and, for an experimentally relevant example of parameter analysis, investigate the sensitivity of these oscillations to changes in the viscosity contrast between the constituent fluids and the inlet flow rates. For the case of stratified laminar flow, we detail a physically realizable set of network parameters that exhibit rich dynamics. The tools and results developed here are general and could be applied to other physical systems. [ABSTRACT FROM AUTHOR]

We explore the stability of equilibrium solution(s) of a simple model of microvascular blood flow in a two-node network. The model takes the form of convection equations for red blood cell concentration, and contains two important rheological effects-the Fåhræus-Lindqvist effect, which governs viscosity of blood flow in a single vessel, and the plasma skimming effect, which describes the separation of red blood cells at diverging nodes. We show that stability is governed by a linear system of integral equations, and we study the roots of the associated characteristic equation in detail. We demonstrate using a combination of analytical and numerical techniques that it is the relative strength of the Fåhræus-Lindqvist effect and the plasma skimming effect which determines the existence of a set of network parameter values which lead to a Hopf bifurcation of the equilibrium solution. We confirm these predictions with direct numerical simulation and suggest several areas for future research and application. [ABSTRACT FROM AUTHOR]