Linear and non-linear analyses of convection in a micropolar fluid occupying a porous medium

Linear and weakly non-linear analyses of convection in a micropolar fluid occupying a high-porosity medium are performed. The Brinkman–Eringen momentum equation is considered. The linear and non-linear analyses are, respectively, based on the normal mode technique and truncated representation of Fourier series. The linear theory for a two-phase system reiterates that the preferred mode of convection is stationary as in the case of a single-phase system. An autonomous system of differential equations representing cellular convection arising in the study is considered to analyse the critical points. The Nusselt number is obtained as a function of micropolar and porous medium parameters.     

An analytical study of linear and non-linear double diffusive convection with Soret and Dufour effects in couple stress fluid

The onset of double diffusive convection in a two component couple stress fluid layer with Soret and Dufour effects has been studied using both linear and non-linear stability analysis. The linear theory depends on normal mode technique and non-linear analysis depends on a minimal representation of double Fourier series. The effect of couple stress parameter, the Soret and Dufour parameters, and the Prandtl number on the stationary and oscillatory convection are presented graphically. The Dufour parameter enhances the stability of the couple stress fluid system in case of both stationary and oscillatory mode. The effect of positive Soret parameter is to destabilize the system in case of stationary mode while it stabilizes the system in case of oscillatory mode. The negative Soret parameter enhances the stability in both stationary and oscillatory mode. The couple stress parameter enhances the stability of the system in both stationary and oscillatory modes. The Dufour parameter increases the heat transfer while the couple stress parameter has reverse effect. The Soret parameter has negligible influence on heat transfer. Both Dufour and Soret parameters increases the mass transfer while the couple stress parameter has dual effect depending on the value of the Rayleigh number.     

Three dimensional finite element analysis for Rayleight-Benard convection in rectangular enclosures

The pattern of Rayleigh-Benard convection of air in a rectangular box heated-from-below is studied by numerically solving the three-dimensional time-dependent Navier-Stokes equations under the Boussinesq approximation. Slightly supercritical Rayleigh number was adopted to track the evolutions of flow structure as a function of enclosure's aspect ratio (A=L/H). The flow will asymptotically evolve to different patterns, among which, two possible types of flow pattern are found. One consists of the pair of straight vortex rolls and the other appears as closed vortex rings. The transition between the flow patterns indicates that there exists a flow bifurcation with the variation of container's aspect ratio. In addition, both steady and oscillatory flows have been observed, corresponding to the pair of straight vortex rolls and the vortex ring, respectively. The complexity of flow structure tends to increase with the increasing aspect ratio of the rectangular enclosure. The project supported by the National Natural Science Foundation of China (19889210), the National Distinguished Young Fund (10125210), the Hundred Talents Program of CAS, and the Training Program for the Trans-Century Outstanding Young of MOE     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

双流体Rayleigh—Benard对流中的局部行波

Rayleigh-Benard对流是研究非平衡对流的斑图（Pattem）及非线性动力学特性的典型模型之一。据此,通过流体力学基本方程组的数值模拟,探讨了分离比沙-0．4时双流体局部行波的形成过程;讨论了双流体局部行波流速场、温度场、浓度场、平均浓度流的结构和特性,分析了局部行波被局部化的原因。研究结果表明：局部行波始终在背离端壁方向传播并被限定在此端,它存在于腔体左端还是右端取决于经过过渡后对传波中最终控制腔体的一支行波的传播方向;在局部行波的存在范围内（1．519≤r≤1．604）,随相对瑞利数增加,衡量对流振幅的最大垂直流速、特征通过流体层的垂直热流量的努塞尔数、局部行波宽度等都在增加,但反映浓度特性的混合参数减小。     

A unified asymptotic theory of the anelastic approximation in geophysical gases and liquids

An asymptotic theory of the anelastic approximation is developed for fluids having arbitrary equations of state under two assumptions: weak compressibility and small Brunt–Väisälä frequency. We show that both Boussinesq approximation (BA) and anelastic approximation (AA) may be included in a unique quasi-incompressible approximation (QIA) already constructed by Durran for polytropic gases. The only difference between AA and BA is that, in the BA, the equations are with slowly varying coefficients, while in the AA the coefficients are fast varying. Applications are made to atmospheric air and to sea-water.     

Concentration dependence of the linear viscoelastic properties of particle suspensions

The response under small amplitude oscillatory deformations of a suspension of non-Brownian spheres dispersed in a viscoelastic fluid is investigated. The correspondence principle of linear viscoelasticity is used to derive a simple constitutive model from a model for a suspension in a Newtonian liquid. The theory predicts that for a specific particulate system the concentration dependence of the viscoelastic properties should collapse to a single master curve when the values are normalized with those of the carrier fluid alone. Measurements with the micro-Fourier rheometer using oscillatory squeeze flow are carried out on two suspensions of 60 and 80 μm sized particles dispersed in polymeric fluid and in silicon oil, and the master curve is verified. Received: 27 April 1999/Accepted: 15 October 1999     

Quasiperiodic structures in the problem of fluid convection between horizontal planes

A steady flow induced by quasiperiodic temperature distributions at the boundaries is considered with reference to the example of fluid convection in a plane horizontal layer heated from below. A spectral-finite-difference method is employed to obtain convective vortex structures varying quasiperiodically along the horizontal coordinate. The properties of the quasiperiodic solutions, together with their spectra and integral characteristics (Nusselt number) are studied as functions of the Grashof number and the boundary conditions of the problem.     

Rayleigh-Bénard convection heat transfer in nanoparticle suspensions

Natural convection heat transfer of nanofluids in horizontal enclosures heated from below is investigated theoretically. The main idea upon which the present work is based is that nanofluids behave more like a single-phase fluid rather than like a conventional solid-liquid mixture, which implies that all the convective heat transfer correlations available for single-phase flows can be extended to nanoparticle suspensions, provided that the thermophysical properties appearing in them are the nanofluid effective properties calculated at the reference temperature. In this connection, two empirical equations, based on a wide variety of experimental data reported in the literature, are developed for the evaluation of the nanofluid effective thermal conductivity and dynamic viscosity, whereas the other effective properties are evaluated by the traditional mixing theory. The heat transfer enhancement that derives from the dispersion of nano-sized solid particles into the base liquid is calculated for different operating conditions, nanoparticle diameters, and combinations of solid and liquid phases. One of the fundamental results is the existence of an optimal particle loading for maximum heat transfer across the bottom-heated enclosure. In particular, for any assigned combination of suspended nanoparticles and base liquid, it is found that the optimal volume fraction increases as the nanofluid average temperature increases, and may either increase or decrease with increasing the nanoparticle size according as the flow is laminar or turbulent. Moreover, the optimal volume fraction has a peak at a definite value of the Rayleigh number of the base fluid, that depends on both the average temperature of the nanofluid and the diameter of the suspended nanoparticles.     

Non-linear approximations for natural convection in a horizontal annulus

Two different problems are proposed as approximations of the usual system modelling natural convection under the Oberbeck-Boussinesq assumptions. The error is evaluated by means of the norm of its gradient in the Hilbert space. The average Nusselt number is also estimated.     

Thermohydrodynamic characteristics of combined double-diffusive radiation convection heat transfer in a cavity

This work deals with the numerical analysis of a radiating gas flow caused by both temperature and buoyancy concentration gradients in a square cavity; in this regard, the set of governing equations, including conservation of mass, momentum, species, and energy are solved by a numerical technique. In terms of radiation, since the fluid is considered as a semitransparent medium, the radiative term in the energy equation appears and is calculated by numerical solving of the radiative transfer equation (RTE). Furthermore, all of the surrounding cavity walls are considered to be opaque, gray, and diffuse with constant emissivity. All of the flow equations are solved by the finite difference method (FDM) and the RTE by the discrete ordinate one (DOM). In the present study, an attempt is made to verify the optical thickness effects on flow, thermal behavior, and mass transform in a cavity flow, such that reciprocating trends were seen in this manner. Our numerical results show that the thermal field in double-diffusive convection flow reaches very fast its steady-state situation in comparison to the concentration distribution. Besides, it is found that the thermohydrodynamic characteristics of a double-diffusive convection flow of a radiating gas are much affected by optical thickness.     

An analytical study of the non-linear vibrations of cylindrical shells

The primary resonance response of simply supported circular cylindrical shells is investigated using the perturbation method. Donnell's non-linear shallow-shell theory is used to derive the governing partial differential equations of motion. The Galerkin technique is then employed to transform the equations of motion into a set of temporal ordinary differential equations. Considering only the primary resonance case, the method of multiple scales is used to study the periodic solutions and their stability. The necessary and sufficient conditions for appearance of the so-called companion mode are also discussed. To this end, a range of the possible multi-mode solution is obtained for response and excitation amplitudes and also excitation frequency as a function of damping, geometry and material properties of the shell. Parametric studies are performed to illustrate the effect of different values of thickness, length and material composition on the possibility of the companion mode participation in primary resonance response.     

Natural convection heat transfer in the annular region between porous confocal ellipses

The Darcy–Boussinesq equations are solved in two dimensions and in elliptical cylindrical co‐ordinates using a second‐order‐accurate finite difference code and a very fine grid. For the limiting case of a circular geometry, the results show that a hysteresis loop is possible for some values of the radius ratio, in agreement both with previous calculations using cylindrical co‐ordinates and with the available experimental data. For the general case of an annulus of elliptical cross‐section, two configurations, blunt or slender, are considered. When the major axes are horizontal (blunt case) a hysteresis loop appears for a certain range of Raleigh numbers. For the slender configuration, when the major axes are vertical, a transition from a steady to a periodic regime (Hopf bifurcation) has been evidenced. In all cases, the heat transfer rate from the slender geometry is greater than that obtained in the blunt case. Copyright © 1999 John Wiley & Sons, Ltd.     

Diffuse approximation method for solving natural convection in porous media

The diffuse approximation is presented and applied to natural convection problems in porous media. A comparison with the control volume-based finite-element method shows that, overall, the diffuse approximation appears to be fairly attractive.Nomenclature H height of the cavities - I functional - K permeability - p(M _i ,M) line vector of monomials - p ^T p-transpose - M current point - Nu Nusselt number - Ri inner radius - Ro outer radius - Ra Rayleigh number - x, y cartesian coordinates - u, v velocity components - T temperature - M vector of estimated derivatives - _t thermal diffusivity - coefficient of thermal expansion - practical aperture of the weighting function - scalar field - (M, M _i ) weighting function - streamfunction - kinematic viscosity     

Non-Newtonian end effects in falling ball viscometry of concentrated suspensions

In a Newtonian fluid contained in a cylinder, a small ball initially at rest released just below the surface would accelerate to achieve a steady-state velocity within one cylinder diameter. After traversing the center section of the cylinder, the ball would begin slowing down within one cylinder diameter of the bottom. This behavior is also observed in suspensions where the size of the suspended particles is small relative to the containing cylinder. However, in concentrated suspensions of larger suspended particles, balls released near the upper surface travel faster than the steady state velocity. In addition, the length of the upper surface end effect, where the falling ball decelerates to the steady state velocity, and the lower end effect zone, where the ball decelerates to rest at the bottom, is many times longer than in a Newtonian single-phase liquid. These non-Newtonian end effects are reduced if the suspended particles are polydisperse in their size distribution.     

An example of entropy balance in natural convection,Part 1: the usual Boussinesq equations

Numerical simulations of natural convection in cavities performed with the usual Boussinesq equations result in an unbalanced irreversibility budget. Thermodynamic analysis shows that these equations represent a system that exchanges with the surroundings, not only two heat fluxes, but also two fluxes of mechanical energy: an input, that generates the fluid motion, and an output, due to viscous friction. After this analysis, the thermodynamic discrepancies can be explained. To cite this article: M. Pons, P. Le Quéré, C. R. Mecanique 333 (2005).     

Application of a non-linear local analysis method for the problem of mixed convection instability

We consider the problem of laminar mixed convection flow between parallel, vertical and uniformly heated plates where the governing dimensionless parameters are the Prandtl, Rayleigh and Reynolds numbers. Using the method based on the centre manifold theorem which was derived from the general theory of dynamical systems, we reduce a three-dimensional simplified model of ordinary differential amplitude equations emanating from the original Navier-Stokes system of the problem in the vicinity of a trivial stationary solution. We have found that when the forcing parameter, the Rayleigh number, increases beyond the critical value Ra_s, the stationary solution is a pitchfork bifurcation point of the system.     

An example of entropy balance in natural convection,Part 2: the thermodynamic Boussinesq equations

Numerical simulations of natural convection performed with the usual Boussinesq equations result in unbalanced irreversibility budget. The thermodynamic Boussinesq equations solve this problem, especially because they simulate production of kinetic energy within the fluid through its expansion and contraction. These fluid volume changes, without which natural convection would not occur, also induce heat transfer by piston effect. The piston effect, which appears then as an intrinsic component of buoyancy-induced natural convection, introduces the non-dimensional adiabatic temperature gradient as a control parameter of natural convection. To cite this article: M. Pons, P. Le Quéré, C. R. Mecanique 333 (2005).     

The Effect of Thermal Expansion on Porous Media Convection Part 1: Thermal Expansion Solution

The effect of thermal expansion on porous media convection is investigated by isolating first the solution of thermal expansion in the absence of convection which allows to evaluate the leading order effects that need to be included in the convection problem that is solved later. A relaxation of the Boussinesq approximation is applied and the relevant time scales for the formulated problem are identified from the equations as well as from the derived analytical solutions. Particular attention is paid to the problem of waves propagation in porous media and a significant conceptual difference between the isothermal compression problem in flows in porous media and its non-isothermal counterpart is established. The contrast between these two distinct problems, in terms of the different time scales involved, is evident from the results. While the thermal expansion is identified as a transient phenomenon, its impact on the post-transient solutions is found to be sensitive to the symmetry of the particular temperature initial conditions that are applied.     

Generalization of the Lorenz model to the two-dimensional convection of second-grade fluids

We apply the truncation of the Navier-Stokes-Fourier equations which leads to the Lorenz model, to the investigation of second-grade fluids. The new set of equations proves to work as an approximated approach to 2D-convective dynamics, under the same restrictions as for Newtonian fluids. The different behaviour depends only on α₁ and consists in a “modified” Prandtl number.