Stability of Poiseuille flow of a Bingham fluid overlying an anisotropic and inhomogeneous porous layer

2019 ◽  
Vol 874 ◽  
pp. 573-605 ◽  
Author(s):  
Sourav Sengupta ◽  
Sirshendu De

Modal and non-modal stability analyses are performed for Poiseuille flow of a Bingham fluid overlying an anisotropic and inhomogeneous porous layer saturated with the same fluid. In the case of modal analysis, the resultant Orr–Sommerfeld type eigenvalue problem is formulated and solved via the Chebyshev collocation method, using QZ decomposition. It is found that no unstable eigenvalues are present for the problem, indicating that the flow is linearly stable. Therefore, non-modal analysis is attempted in order to observe the short-time response. For non-modal analysis, the initial value problem is solved, and the response of the system to initial conditions is assessed. The aim is to evaluate the effects on the flow stability of porous layer parameters in terms of depth ratio (ratio of the fluid layer thickness $d$ to the porous layer thickness $d_{m}$), Bingham number, Darcy number and slip coefficient. The effects of anisotropy and inhomogeneity of the porous layer on flow transition are also investigated. In addition, the shapes of the optimal perturbations are constructed. The mechanism of transient growth is explored to comprehend the complex interplay of various factors that lead to intermediate amplifications. The present analysis is perhaps the first attempt at analysing flow stability of viscoplastic fluids over a porous medium, and would possibly lead to better and efficient designing of flow environments involving such flow.

1994 ◽  
Vol 263 ◽  
pp. 133-150 ◽  
Author(s):  
I. A. Frigaard ◽  
S. D. Howison ◽  
I. J. Sobey

The stability to linearized two-dimensional disturbances of plane Poiseuille flow of a Bingham fluid is considered. Bingham fluids exhibit a yield stress in addition to a plastic viscosity and this description is typically applied to drilling muds. A non-zero yield stress results in an additional parameter, a Bingham number, and it is found that the minimum Reynolds number for linear instability increases almost linearly with increasing Bingham number.


2017 ◽  
Vol 826 ◽  
pp. 376-395 ◽  
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang

This paper performs a linear stability analysis to investigate the stability of plane Poiseuille–Couette flow in a fluid layer overlying a porous medium saturated with the same fluid. The effect of superimposed Couette flow on the associated Poiseuille flow in such a two-layer system is explored carefully. The result shows that the presence of Couette flow may destabilize the Poiseuille flow at small depth ratio $\hat{d}$, defined by the ratio of the depth of the fluid layer to the depth of the porous layer, and induce a tri-modal structure to the neutral curves. At moderate $\hat{d}$, the Couette component generally produces a stabilization effect on the flow. When the velocity of the upper moving wall is large enough, a bi-modal behaviour of the neutral curves appears and a shift of instability mode occurs from the long-wave fluid-layer mode to the porous-layer mode with higher wavenumber. These stability characteristics are remarkably different from those of the plane Poiseuille–Couette flow in a single fluid layer in that the flow becomes absolutely stable when the wall velocity is over 70 % of the maximum velocity of the Poiseuille component of flow. The stability of pure Couette flow in such a two-layer system is also studied. It is found that the flow is still absolutely stable with respect to infinitesimal disturbances, which is the same as the stability characteristic of a single-layer plane Couette flow.


1989 ◽  
Vol 111 (2) ◽  
pp. 357-362 ◽  
Author(s):  
M. E. Taslim ◽  
U. Narusawa

The results of stability analyses for the onset of convective motion are reported for the following three horizontally superposed systems of porous and fluid layers: (a) a porous layer sandwiched between two fluid layers with rigid top and bottom boundaries, (b) a fluid layer overlying a layer of porous medium, and (c) a fluid layer sandwiched between two porous layers. By changing the depth ratio dˆ from zero to infinity, a set of stability criteria (i.e., the critical Rayleigh number Rac and the critical wave number ac) is obtained, ranging from the case of a fluid layer between two rigid boundaries to the case of a porous layer between two impermeable boundaries. The effects of k/km (the thermal conductivity ratio), δ (the square root of the Darcy number), and α (the nondimensional proportionality constant in the slip condition) on Rac and ac are also examined in detail. The results in this paper, combined with those reported previously for Case (a) (Pillatsis et al., 1987), will provide a comprehensive picture of the interaction between a porous and a fluid layer.


Author(s):  
M. Subotic ◽  
F. C. Lai

Flow and temperature fields in an annulus between two rotating cylinders have been examined in this study. While the outer cylinder is stationary, the inner cylinder is rotating with a constant angular speed. A homogeneous and isotropic porous layer is press-fit to the inner surface of the outer cylinder. The porous sleeve is saturated with the fluid that fills the annulus. The Brinkman-extended Darcy equations are used to model the flow in the porous layer while Navier-Stokes equations are used for the fluid layer. The conditions applied at the interface between the porous and fluid layers are the continuity of temperature, heat flux, tangential velocity and shear stress. Analytical solutions have been attempted. Through these solutions, the effects of Darcy number, Brinkman number, and porous sleeve thickness on the velocity profile and temperature distribution are studied.


2012 ◽  
Vol 710 ◽  
pp. 505-544 ◽  
Author(s):  
R. Liu ◽  
Q. S. Liu

AbstractThe problem of an incompressible flow between two coaxial cylinders with radii$a$and$b$subjected to a sliding motion of the inner cylinder in the axial direction is considered. The energy stability and the non-modal stability have been investigated for both axisymmetric and non-axisymmetric disturbances. For the non-modal stability, we focus on two problems: response to external excitations and response to initial conditions. The former is studied by examining the$\epsilon $-pseudospectrum, and the latter by examining the energy growth function$G(t)$. Unlike the results of the modal analysis in which the stability of sliding Couette flow is determined by axisymmetric disturbances, the energy analysis shows that a non-axisymmetric disturbance has a critical energy Reynolds number for all radius ratios$\eta = a/ b$. The results for non-modal stability show that rather large transient growth occurs over a wide range of azimuthal wavenumber$n$and streamwise wavenumber$\ensuremath{\alpha} $, even though the Reynolds number is far below its critical value. For the problem of response to external excitations, the response is sensitive to low-frequency external excitations. For all values of the radius ratio, the maximum response is achieved by non-axisymmetric and streamwise-independent disturbances when the frequency of external forcing$\omega = 0$. For the problem of response to initial conditions, the optimal disturbance is in the form of helical streaks at low Reynolds numbers. With the increase of$\mathit{Re}$, the optimal disturbance becomes very close to straight streaks. For each$\eta $, the maximum energy growth of streamwise-independent disturbances is of the order of${\mathit{Re}}^{2} $, and the optimal time is of the order of$\mathit{Re}$. This relation is qualitatively similar to that for plane Couette flow, plane Poiseuille flow and pipe Poiseuille flow. Direct numerical simulations are applied to investigate the transition of the streamwise vortex (SV) scenario at$\mathit{Re}= 1000$and 1500 for various$\eta $. The initial disturbances are the optimal streamwise vortices predicted by the non-modal analysis. We studied the streak breakdown phase of the SV scenarios by examining the instability of streaks. Moreover, we have investigated the sustainment of the energy of disturbances in the SV scenario.


Author(s):  
C. C. Ngo ◽  
F. C. Lai

Heat transfer induced by buoyancy from a pipe buried in a semi-infinite porous medium with a superimposed fluid layer has been numerically examined in this study. Due to the complexity involved, finite difference method along with body-fitted coordinate systems has been employed. The Brinkman-extended Darcy equations are used to model flow in the porous medium while Navier-Stokes equations are used for the fluid layer. The conditions applied at the interface between the fluid and porous layers are the continuity of temperature, heat flux, normal and tangential velocity, shear stress and pressure. A parametric study has been performed to investigate the effects of Rayleigh number, Prandtl number, Darcy number, and fluid layer thickness on the flow patterns and heat transfer rates. The results show that heat transfer increases with the Rayleigh number, but the convective strength decreases with the Darcy number. The heat transfer rate is smaller when the superimposed fluid is air instead of water. For a porous layer with Da ≤ 0.0005 and an overlaying fluid layer thickness of L/ri ≥ 1, convection is initiated in the fluid layer and it may develop into multiple recirculating cells at a moderate Rayleigh number (i.e., Ra ≤ 104), and may further develop into a single cell at a higher Rayleigh number of 105.


2007 ◽  
Vol 577 ◽  
pp. 211-239 ◽  
Author(s):  
C. NOUAR ◽  
N. KABOUYA ◽  
J. DUSEK ◽  
M. MAMOU

The receptivity problem of plane Bingham–Poiseuille flow with respect to weak perturbations is addressed. The relevance of this study is highlighted by the linear stability analysis results (spectra and pseudospectra). The first part of the present paper thus deals with the classical normal-mode approach in which the resulting eigenvalue problem is solved using the Chebychev collocation method. Within the range of parameters considered, the Poiseuille flow of Bingham fluid is found to be linearly stable. The second part investigates the most amplified perturbations using the non-modal approach. At a very low Bingham number (B ≪ 1), the optimal disturbance consists of almost streamwise vortices, whereas at moderate or large B the optimal disturbance becomes oblique. The evolution of the obliqueness as function of B is determined. The linear analysis presented also indicates, as a first stage of a theoretical investigation, the principal challenges of a more complete nonlinear study.


2008 ◽  
Vol 130 (10) ◽  
Author(s):  
M. Subotic ◽  
F. C. Lai

Flow and temperature fields in an annulus between two rotating cylinders have been examined in this study. While the outer cylinder is stationary, the inner cylinder is rotating with a constant angular speed. A homogeneous and isotropic porous layer is press fit to the inner surface of the outer cylinder. The porous sleeve is saturated with the fluid that fills the annulus. The Brinkman-extended Darcy equations are used to model the flow in the porous layer while the Navier–Stokes equations are used for the fluid layer. The conditions applied at the interface between the porous and fluid layers are the continuity of temperature, heat flux, tangential velocity, and shear stress. Analytical solutions have been attempted. Through these solutions, the effects of Darcy number, Brinkman number, and porous sleeve thickness on the velocity profile and temperature distribution are studied.


Author(s):  
Joseph J. Webber ◽  
Herbert E. Huppert

AbstractMotivated by shallow ocean waves propagating over coral reefs, we investigate the drift velocities due to surface wave motion in an effectively inviscid fluid that overlies a saturated porous bed of finite depth. Previous work in this area either neglects the large-scale flow between layers (Phillips in Flow and reactions in permeable rocks, Cambridge University Press, Cambridge, 1991) or only considers the drift above the porous layer (Monismith in Ann Rev Fluid Mech 39:37–55, 2007). Overcoming these limitations, we propose a model where flow is described by a velocity potential above the porous layer and by Darcy’s law in the porous bed, with derived matching conditions at the interface between the two layers. Both a horizontal and a novel vertical drift effect arise from the damping of the porous bed, which requires the use of a complex wavenumber k. This is in contrast to the purely horizontal second-order drift first derived by Stokes (Trans Camb Philos Soc 8:441–455, 1847) when working with solely a pure fluid layer. Our work provides a physical model for coral reefs in shallow seas, where fluid drift both above and within the reef is vitally important for maintaining a healthy reef ecosystem (Koehl et al. In: Proceedings of the 8th International Coral Reef Symposium, vol 2, pp 1087–1092, 1997; Monismith in Ann Rev Fluid Mech 39:37–55, 2007). We compare our model with field measurements by Koehl and Hadfield (J Mar Syst 49:75–88, 2004) and also explain the vertical drift effects as documented by Koehl et al. (Mar Ecol Prog Ser 335:1–18, 2007), who measured the exchange between a coral reef layer and the (relatively shallow) sea above.


Sign in / Sign up

Export Citation Format

Share Document