Coriolis effect on thermal convection in a couple-stress fluid-saturated rotating rigid porous layer

Both linear and weakly nonlinear stability analyses are performed to study thermal convection in a rotating couple-stress fluid-saturated rigid porous layer. In the case of linear stability analysis, conditions for the occurrence of possible bifurcations are obtained. It is shown that Hopf bifurcation is possible due to Coriolis force, and it occurs at a lower value of the Rayleigh number at which the simple bifurcation occurs. In contrast to the nonrotating case, it is found that the couple-stress parameter plays a dual role in deciding the stability characteristics of the system, depending on the strength of rotation. Nonlinear stability analysis is carried out by constructing a set of coupled nonlinear ordinary differential equations using truncated representation of Fourier series. Sub-critical finite amplitude steady motions occur depending on the choice of physical parameters but at higher rotation rates oscillatory convection is found to be the preferred mode of instability. Besides, the stability of steady bifurcating equilibrium solution is discussed using modified perturbation theory. Heat transfer is calculated in terms of Nusselt number. Also, the transient behavior of the Nusselt number is investigated by solving the nonlinear differential equations numerically using the Runge–Kutta–Gill method. It is noted that increase in the value of Taylor number and the couple-stress parameter is to dampen the oscillations of Nusselt number and thereby to decrease the heat transfer.     

Thermal instability in a rotating anisotropic porous layer saturated by a viscoelastic fluid

Linear and non-linear thermal instability in a rotating anisotropic porous medium, saturated with viscoelastic fluid, has been investigated for free-free surfaces. The linear theory is being related to the normal mode method and non-linear analysis is based on minimal representation of the truncated Fourier series analysis containing only two terms. The extended Darcy model, which includes the time derivative and Coriolis terms has been employed in the momentum equation. The criteria for both stationary and oscillatory convection is derived analytically. The rotation inhibits the onset of convection in both stationary and oscillatory modes. A weak non-linear theory based on the truncated representation of Fourier series method is used to find the thermal Nusselt number. The transient behaviour of the Nusselt number is also investigated by solving the finite amplitude equations using a numerical method. The results obtained during the analysis have been presented graphically.     

Linear instability of the electrified free interface between two cylindrical shells of viscoelastic fluids through porous media

In this paper, we have discussed the linear stability analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interfacial insoluble surfactant and surface charge through porous media. The case of long waves interfacial stability has been studied. The dispersion relation is solved numerically and hence the effects of various parameters are illustrated graphically. Our results reveal that the influence of the physicochemical parameter β is to shrink the instability region of the surface and reduce the growth rate of the unstable normal modes. Such important effects of the surfactant on the shape of interfacial structures are more sensitive to the variation of the β corresponding to non-Newtonian fluids-model compared with the Newtonian fluids model. In the case of long wave limit, it is demonstrated that increasing β, has a dual role in-fluence (de-stabilizing effects) depending on the viscosity of the core fluid. It has a destabilizing effect at the large values of the core fluid viscosity coefficient, while this role is exchanged to a regularly stabilizing influence at small values of such coefficient.     

The convective instability of Maxwell fluid-saturated porous layer using a thermal non-equilibrium model

A linear stability analysis is carried out to study viscoelastic fluid convection in a horizontal porous layer heated from below and cooled from above when the solid and fluid phases are not in a local thermal equilibrium. The modified Darcy–Brinkman–Maxwell model is used for the momentum equation and two-field model is used for the energy equation each representing the solid and fluid phases separately. The conditions for the onset of stationary and oscillatory convection are obtained analytically. Linear stability analysis suggests that, there is a competition between the processes of viscoelasticity and thermal diffusion that causes the first convective instability to be oscillatory rather than stationary. Elasticity is found to destabilize the system. Besides, the effects of Darcy number, thermal non-equilibrium and the Darcy–Prandtl number on the stability of the system are analyzed in detail.     

Bounds on convective heat transport in a rotating porous layer

Using the background field variational method, bounds on convective heat transport in a rotating porous layer heated from below are derived from the primitive equations. The enhancement of heat transport beyond the minimal conduction value (the Nusselt number Nu) is bounded in terms of the dimensionless temperature difference across the layer (the Rayleigh number Ra) according to

This rigorous upper bound shows that rotation has a retarding effect on convective heat transport.     

Entry flow behaviour of viscoelastic fluids in an annulus

Developing and fully developed velocity profiles in the entrance region of an abrupt 2-to-1 annular contraction were measured for a number of visco-elastic polymer solutions. Experimental results were obtained for Reynolds number and flow behaviour index in the range 9.8 ? Re ? 355 and 0.372 ? n ? 0.55 respectively. A momentum-energy integral technique was employed in the boundary layer analysis. The deviation from inelastic behaviour was indicated by the ratio of elastic to inertial forces, Ws/Re. Within the limits of confidence of the experimental results, good agreement with theoretical predictions was obtained and very little deviation from inelastic behaviour was observed for Ws/Re < 0.08. For the test fluids investigated, the entrance length was found to be longer than that predicted for the corresponding inelastic fluids of the same n.     

Thermohaline instability in a rotating anisotropic porous medium

The onset of thermohaline convection in an anisotropic rotating porous layer of infinite horizontal extent is investigated. Numerical computations are made assuming horizontal isotropy in permeability. It is observed that (i) for certain parameters a bottom-heavy arrangement destabilises a rotating anisotropic porous layer, (ii) the lower the anisotropy parameter, the higher the range of bottom-heavy solute gradient for which there is destabilisation, (iii) increase in the anisotropy parameter stabilises the system, and (iv) for some values of the parameters rotation destabilises the system, though in general, it has a stabilising effect.     

Flow of viscoelastic fluids between rotating disks

Few boundary-value problems in fluid mechanics can match the attention that has been accorded to the flow of fluids, Newtonian and non-Newtonian, between parallel rotating disks rotating about a common axis or about distinct axes. An interesting feature which has been recently observed is the existence of solutions that are not axially symmetric even in the case of flow due to the rotation of disks about a common axis. In this article we review the recent efforts that have been expended in the study of both symmetric and asymmetric solutions in the case of both the classical linearly viscous fluid and viscoelastic fluids.The support of the Air Force Office of Scientific Research is gratefully acknowledged.     

Dynamics of convective thermal waves in a porous continuum

A one-dimensional motion of an incompressible fluid (gas) through a fixed bed of a granular material under a short thermal pulse at the bed inlet is considered. An analytical dimensionless solution of the boundary-value problem is presented in quadratures. Specific features of the propagation of convective thermal waves due to the fluid (gas) flow are studied. The data on the propagation velocity, amplitude, temperature difference between the fluid (gas) and the bed, and the length of the zone of intense heat transfer are obtained. The examples of the analysis of a continuum flow in engineering apparatuses which use convective thermal waves are presented.     

Nonlinear convective stability problems of viscoelastic fluids in finite domains

A Chebyshev pseudospectral method is generalized to solve the nonlinear hydrodynamic stability problems of Rayleigh-Bénard convection of viscoelastic fluids in finite domains, which are compatible with the experimental situations, for the range of viscoelastic parameters where the exchange of stabilities is valid. The effects of box aspect ratio, the Deborah number 5 and the dimensionless retardation time ) on the critical Rayleigh number and convection intensity are investigated. The comparison of these results with the experimental data might be used to guide the selection of constitutive equations and to estimate viscoelastic parameter values. The present technique of hydrodynamic stability analysis is quite versatile and can be employed to solve other hydrodynamic stability problems in finite domains.     

Electrohydrodynamic instability of two superposed Walters B´ viscoelastic fluids in relative motion through porous medium

Summary  The electrohydrodynamic Kelvin–Helmholtz instability of the interface between two uniform superposed viscoelastic (B′ model) dielectric fluids streaming through a porous medium is investigated. The considered system is influenced by applied electric fields acting normally to the interface between the two media, at which there are no surface charges present. In the absence of surface tension, perturbations transverse to the direction of streaming are found to be unaffected by either streaming and applied electric fields for the potentially unstable configuration, or streaming only for the potentially stable configuration, as long as perturbations in the direction of streaming are ignored. For perturbations in all other directions, there exists instability for a certain wavenumber range. The instability of this system can be enhanced (increased) by normal electric fields. In the presence of surface tension, it is found also that the normal electric fields have destabilizing effects, and that the surface tension is able to suppress the Kelvin–Helmholtz instability for small wavelength perturbations, and the medium porosity reduces the stability range given in terms of the velocities difference and the electric fields effect. Finally, it is shown that the presence of surface tension enhances the stabilizing effect played by the fluid velocities, and that the kinematic viscoelasticity has a stabilizing as well as a destabilizing effect on the considered system under certain conditions. Graphics have been plotted by giving numerical values to the parameters, to depict the stability characteristics. Received 27 March 2000; accepted for publication 3 May 2001     

On the convective instability of a thermal boundary layer

When the surface temperature of a liquid is a harmonic function of time with a frequency, a temperature wave propagates into the liquid. The amplitude of this wave decreases exponentially with distance from the surface. The temperature oscillation is essentially concentrated in a layer of the order of (2/)^1/2, where x is the thermal conductivity of the liquid (thermal boundary layer). Depending on the phase, at certain positions below the surface the temperature gradient is directed downwards and if its magnitude is sufficiently large (the magnitude is a function of the amplitude and frequency of the surface oscillations) the liquid can become unstable with respect to the onset of convection. In that case the convective motion may spread beyond the initial unstable layer. For low frequencies the stability condition can be derived from the usual static Rayleigh criterion, on the basis of the Rayleigh number and the average temperature gradient of the unstable layer. This quasi-static approach, used by Sal'nikov [1], is appropriate to those cases in which the period of the temperature oscillations is much larger than the characteristic time of the perturbations. But when these times are of the same order, the problem must be analyzed in dynamic terms. The stability problem must then be formulated as a problem of parametricresonance excitation of velocity oscillations due to the action of a variable parameter-the temperature gradient.In an earlier work [2] we considered the problem of the stability of a horizontal layer of liquid with a periodically varying temperature gradient. It was assumed that the thickness of the layer was much smaller than the penetration depth of the thermal wave, so that the temperature gradient could be assumed to be independent of position. In the present work we consider the opposite case, in which the liquid layer is assumed to be much larger than the penetration depth, i. e., a thermal boundary layer can be defined. The temperature gradient at equilibrium, which is a parameter in the equations determining the onset of perturbations, is here a periodic function of time and a relatively complicated function of the depth coordinate z. The periodic oscillations are solved by the Fourier method; the equations for the amplitudes are solved by the approximate method of KarmanPohlhausen.The authors are grateful to L. G. Loitsyanskii for helpful criticism.     

On thermal effects in a special class of viscoelastic fluids

Summary A brief account of a non-isothermal theory of a restricted class of viscoelastic fluids with fading memory followed by its application to some special viscometric flows.     

Steady flow of a maxwell fluid in a porous cylindrical annulus with circular cross-section

Summary When a fluid with memory is injected into any flow region some assumptions regarding the initial state of stress have to be made in order to determine the state of stress at any subsequent instant. For a Maxwell fluid, it is assumed that the fluid near the surface of injection is suddenly stressed and responds by starting flow in accordance with the mechanical model chosen. The flow of a Maxwell fluid with a single relaxation time has been determined under the above assumption in the following two cases: (i) annulus between two porous concentric circular cylinders, and (ii) space between two porous and infinitely extending parallel plates. The nature of flow in the present case is similar to that of the Reiner-Rivlin fluids obtained by Narasimhan²).     

Analysis of the convective instability of a binary mixture in a porous medium

The problem of the stability of a binary mixture in a porous medium is investigated in the complete formulation with allowance for cross kinetic and gravitational effects. Boundary conditions of the first and second kinds for a plane horizontal layer of the porous medium are considered. The boundaries of the region of instability are determined. The region of the parameters corresponding to the stability paradox effect, i.e., the instability of a mixture that becomes heavier with depth, is described. It is established that the multicomponent nature of the mixture helps to stabilize the equilibrium state.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 110–119, January–February, 1993.     

Thermal instability in an anisotropic rotating porous medium

 Influenced by the article of Vadasz [1], an analysis has been carried out to investigate convective instability due to centrifugal acceleration in an anisotropic porous medium. Results reveal that anisotropy in thermal diffusivity destabilizes the system whereas that in permeability has the opposite effect. Received on 26 February 1999     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

Experimental and theoretical studies of convective heat transfer in a cylindrical porous medium

Convective heat transfer at constant heat flux through unconsolidated porous media has been studied both experimentally and theoretically. Heat transfer measurements have been performed for convective heat transfer over a wide range of operational parameters at constant heat fluxes. In addition to heat transfer coefficients, pressure drop and temperature profiles both in radial and axial direction have been recorded. The equations of motion and energy which account for the non-Darcian effect are used to describe the flow and convective heat transfer through the porous medium. Mathematical models for the prediction of heat transfer coefficients and temperature profiles are presented which predict the experimental data with good accuracy.     

Coriolis effect on responses of rotating thin piezoelectric hollow cylinder

Coriolis effect is considered in the analysis of a rotating piezoelectric hollow cylinder. An inhomogeneous Bessel equation governing the radial mechanical displacement is derived, which can be approximated as an Euler type differential equation when the cylinder is very thin. Numerical examples show that the Coriolis effect can be significant under certain conditions.