Davide Picchi

Characteristics of stratified flows of Newtonian/non-Newtonian shear-thinning fluids

Abstract Exact solutions for laminar stratified flows of Newtonian/non-Newtonian shear-thinning fluids in horizontal and inclined channels are presented. An iterative algorithm is proposed to compute the laminar solution for the general case of a Carreau non-Newtonian fluid. The exact solution is used to study the effect of the rheology of the shear-thinning liquid on two-phase flow characteristics considering both gas/liquid and liquid/liquid systems. Concurrent and counter-current inclined systems are investigated, including the mapping of multiple solution boundaries. Aspects relevant to practical applications are discussed, such as the insitu hold-up, or lubrication effects achieved by adding a less viscous phase. A characteristic of this family of systems is that, even if the liquid has a complex rheology (Carreau fluid), the two-phase stratified flow can behave like the liquid is Newtonian for a wide range of operational conditions. The capability of the two-fluid model to yield satisfactory predictions in the presence of shear-thinning liquids is tested, and an algorithm is proposed to a priori predict if the Newtonian (zero shear rate viscosity) behaviour arises for a given operational conditions in order to avoid large errors in the predictions of flow characteristics when the power-law is considered for modelling the shear-thinning behaviour. Two-fluid model closures implied by the exact solution and the effect of a turbulent gas layer are also addressed.

Stability of stratified two-phase channel flows of Newtonian/non-Newtonian shear-thinning fluids

Abstract Linear stability of horizontal and inclined stratified channel flows of Newtonian/non-Newtonian shear-thinning fluids is investigated with respect to all-wavelength perturbations. The Carreau model has been chosen for the modeling of the rheology of a shear-thinning fluid, owing to its capability to describe properly the constant viscosity limits (Newtonian behavior) at low and high shear rates. The results are presented in the form of stability boundaries on flow pattern maps (with the phases’ superficial velocities as coordinates) for several practically important gas–liquid and liquid–liquid systems. The stability maps are accompanied by spatial profiles of the critical perturbations, along with the distributions of the effective and tangent viscosities in the non-Newtonian layer, to show the influence of the complex rheological behavior of shear-thinning liquids on the mechanisms responsible for triggering instability. Due to the complexity of the considered problem, a working methodology is proposed to alleviate the search for the stability boundary. Implementation of the proposed methodology helps to reveal that in many practical cases the investigation of the simpler Newtonian problem is sufficient for the prediction of the exact (non-Newtonian) stability boundary of smooth stratified flow (i.e., in case of horizontal gas–liquid flow). Therefore, the knowledge gained from the stability analysis of Newtonian fluids is applicable to those (usually highly viscous) non-Newtonian systems. Since the stability of stratified flow involving highly viscous Newtonian liquids has not been researched in the literature, interesting findings on the viscosity effects are also obtained. The study highlights the limitations of applying the simpler and widely used power-law model for characterizing the shear-thinning behavior of the liquid. That model would predict a rigid layer (infinite viscosity) at the interface, where the shear rates in the viscous liquid are low, and thereby unphysical representation of the interaction between the phases and its impact on the stability boundary.

Bistability of buoyancy-driven exchange flows in vertical tubes

Buoyancy-driven exchange flows are common to a variety of natural and engineering systems ranging from persistently active volcanoes to counterflows in oceanic straits. Experiments of exchange flows in closed vertical tubes have been used as surrogates to elucidate the basic features of such flows. The resulting data have historically been analyzed and interpreted through core-annular flow solutions, the most common flow configuration at finite viscosity contrasts. These models have been successful in fitting experimental data, but less effective at explaining the variability observed in natural systems. In this paper, we formulate a core-annular solution to the classical problem of buoyancy-driven exchange flows in vertical tubes. The model posits the existence of two mathematically valid solutions, i.e. thin- and thick-core solutions. The theoretical existence of two solutions, however, does not necessarily imply that the system is bistable in the sense that flow switching may occur. Using direct numerical simulations, we test the hypothesis that core-annular flow in vertical tubes is bistable, which implies that the realized flow field is not uniquely defined by the material parameters of the flow. Our numerical experiments, which fully predict experimental data without fitting parameters, demonstrate that buoyancy-driven exchange flows are indeed inherently bistable systems. This finding is consistent with previous experimental data, but in contrast to the underlying hypothesis of previous analytical models that the solution is unique and can be identified by maximizing the flux or extremizing the dissipation in the system. These results have important implications for data interpretation by analytical models, and may also have relevant ramifications for understanding volcanic degassing.