Analytical and numerical stability analysis of Soret-driven convection in a horizontal porous layer

We present an analytical and numerical stability analysis of Soret-driven convection in a porous cavity saturated by a binary fluid. Both the mechanical equilibrium solution and the monocellular flow obtained for particular ranges of the physical parameters of the problem are considered. The porous cavity, bounded by horizontal infinite or finite boundaries, is heated from below or from above. The two horizontal plates are maintained at different constant temperatures while no mass flux is imposed. The influence of the governing parameters and more particularly the role of the separation ratio, characterizing the Soret effect and the normalized porosity, are investigated theoretically and numerically. From the linear stability analysis, we find that the equilibrium solution loses its stability via a stationary bifurcation or a Hopf bifurcation depending on the separation ratio and the normalized porosity of the medium. The role of the porosity is important, when it decreases, the stability of the equilibrium solution is reinforced. For a cell heated from below, the equilibrium solution loses its stability via a stationary bifurcation when the separation ratio >0(Le,), while for 0, while a stationary or an oscillatory bifurcation occurs if mono the monocellular flow loses stability via a Hopf bifurcation. As the Rayleigh number increases, the resulting oscillatory solution evolves to a stationary multicellular flow. For a cell heated from above and <0, the monocellular flow remains linearly stable. We verified numerically that this problem admits other stable multicellular stationary solutions for this range of parameters

Effect of a finite external heat transfer coefficient on the Darcy-Benard instability in a vertical porous cylinder

The onset of thermal convection in a vertical porous cylinder is studied by considering the heating from below and the cooling from above as caused by external forced convection processes. These processes are parametrised through a finite Biot number, and hence through third-kind, or Robin, temperature conditions imposed on the lower and upper boundaries of the cylinder. Both the horizontal plane boundaries and the cylindrical sidewall are assumed to be impermeable; the sidewall is modelled as a thermally insulated boundary. The linear stability analysis is carried out by studying separable normal modes, and the principle of exchange of stabilities is proved. It is shown that the Biot number does not affect the ordering of the instability modes that, when the radius-to-height aspect ratio increases, are displayed in sequence at the onset of convection. On the other hand, the Biot number plays a central role in determining the transition aspect ratios from one mode to its follower. In the limit of a vanishingly small Biot number, just the first (non-axisymmetric) mode is displayed at the onset of convection, for every value of the aspect ratio. (C) 2013 American Institute of Physic

Anomalous Scaling and Solitary Waves in Systems with Non-Linear Diffusion

We study a non-linear convective-diffusive equation, local in space and time, which has its background in the dynamics of the thickness of a wetting film. The presence of a non-linear diffusion predicts the existence of fronts as well as shock fronts. Despite the absence of memory effects, solutions in the case of pure non-linear diffusion exhibit an anomalous sub-diffusive scaling. Due to a balance between non-linear diffusion and convection we, in particular, show that solitary waves appear. For large times they merge into a single solitary wave exhibiting a topological stability. Even though our results concern a specific equation, numerical simulations supports the view that anomalous diffusion and the solitary waves disclosed will be general features in such non-linear convective-diffusive dynamics.Comment: Corrected typos, added 3 references and 2 figure

The effect of an embedded solid block on the onset of convection in a porous cavity

Purpose – The purpose of this paper is to determine how the presence of an embedded, centrally placed, solid but heat-conducting block affects the onset and development of Darcy-Bénard convection. Design/methodology/approach – Steady solutions are obtained using finite difference methods with SOR as the smoother. A detailed presentation is given of how the interface conditions are modelled, and how a continuity of pressure argument is used to determine the value of the streamfunction on the solid block. Findings – The presence of the block affects strongly both the onset of convection and the nonlinear properties such as the mean Nusselt number and the strength of the fluid circulation. The smallest possible critical Darcy-Rayleigh is found to be 22.0152, which is smaller than 4π2, the value when the block is absent. Research limitations/implications – The Darcy-Rayleigh number is restricted to values at or below 200, which is five times the critical value without a solid block, but the size and conductivity of the block vary between all admissible values. Originality/value – This is the first investigation of the effect of internal obstacles on Darcy-Bénard convection. </jats:sec

Effects of temperature-dependent viscosity variation on entropy generation, heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford, is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Péclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case

Incompressible flow in porous media with fractional diffusion

In this paper we study the heat transfer with a general fractional diffusion term of an incompressible fluid in a porous medium governed by Darcy's law. We show formation of singularities with infinite energy and for finite energy we obtain existence and uniqueness results of strong solutions for the sub-critical and critical cases. We prove global existence of weak solutions for different cases. Moreover, we obtain the decay of the solution in $L^p$, for any $p\geq2$, and the asymptotic behavior is shown. Finally, we prove the existence of an attractor in a weak sense and, for the sub-critical dissipative case with $\alpha\in (1,2]$, we obtain the existence of the global attractor for the solutions in the space $H^s$ for any $s > (N/2)+1-\alpha$