Pore-scale modeling of viscoelastic flow in porous media using a Bautista–Manero fluid

This article examines the extensional flow and viscosity and the converging–diverging geometry as the basis of the peculiar viscoelastic behavior in porous media. The modified Bautista–Manero model, which successfully describes elasticity, thixotropic time dependency and shear-thinning, was used for modeling the flow of viscoelastic materials which also show thixotropic attributes. An algorithm, originally proposed by Philippe Tardy, that employs this model to simulate steady-state time-dependent flow was implemented in a non-Newtonian flow simulation code using pore-scale modeling. The simulation results using two topologically-complex networks confirmed the importance of the extensional flow and converging–diverging geometry on the behavior of non-Newtonian fluids in porous media. The analysis also identified a number of correct trends (qualitative and quantitative) and revealed the effect of various fluid and flow parameters on the flow process. The impact of some numerical parameters was also assessed and verified.     

Instability of a liquid–vapor phase transition front in inhomogeneous wettable porous media

The stability of vertical flows through a horizontally extended two-dimensional region of a porous medium is considered in the case of presence of a phase transition front. It is shown that the plane steady-state phase transition front may have several steady-state positions in the wettable porous medium and the necessary condition of their existence is obtained. The spectral stability of the plane phase transition interface is investigated. It is found that in the presence of capillary forces exerted on the phase transition front in the wettable medium the plane front can be destabilized on the mode with both infinite and zero wavenumbers (short- and long-wave instabilities); the short-wave instability can then exist even in the case of the sole steady-state position of the front.     

Stokes＇ first problem for the fourth order fluid in a porous half space

In this study, the flow of a fourth order fluid in a porous half space is modeled. By using the modified Darcy’s law, the flow over a suddenly moving flat plate is studied numerically. The influence of various parameters of interest on the velocity profile is revealed. The English text was polished by Keren Wang.     

Onset of Darcy–Brinkman convection in a binary viscoelastic fluid saturated porous layer

The onset of Darcy–Brinkman double-diffusive convection in a binary viscoelastic fluid-saturated porous layer is studied using both linear and weakly nonlinear stability analyses. The Oldroyd-B model is employed to describe the rheological behavior of the fluid. An extended form of Darcy–Oldroyd law incorporating the Brinkman’s correction and time derivative is used to describe the fluid flow and the Oberbeck–Boussinesq approximation is invoked. The onset criterion for stationary and oscillatory convection is derived analytically. The effects of rheological parameters, Darcy number, normalized porosity, Lewis number, solute Rayleigh number, and Darcy–Prandtl number on the stability of the system is investigated. The results indicated that there is a competition among the processes of thermal, solute diffusions and viscoelasticity that causes the convection to set in through the oscillatory modes rather than the stationary. The Darcy–Prandtl number has a dual effect on the threshold of oscillatory convection. The nonlinear theory based on the method of truncated representation of Fourier series is used to find the transient heat and mass transfer. Some existing results are reproduced as the particular cases of present study.     

On a Darcy–Stefan Problem arising in freezing and thawing of saturated porous media

A model with phase change for material convection in a saturated porous medium with a frozen region is formulated as a Darcy-Stefan problem. We propose a new generalized formulation for this Stefan-type problem with convection governed by Darcy's law. This approach, which is valid for irregular geometries with irregular subregions, has the advantage of not requiring the smoothness of the temperature, that restricted previous mathematical works to two-dimensional particular cases. We show existence of generalized solutions, passing to the limit in suitable approximated problems, which in principle can be solved numerically by the finite element method. Received October 20, 1998     

Homotopy Analysis Method to Walter’s B fluid in a vertical channel with porous wall

In this research the steady three-dimensional flow of a Walter’s B fluid in a vertical channel with porous wall, through which the fluid is injected uniformly into the channel through one side of the channel, is studied analytically using Homotopy Analysis Method (HAM). The channel is assumed to be infinite and uniform. The effects of the elasticity of the fluid on the flow and heat transfer on the walls of the channel are discussed.     

Bio-Marangoni convection ?ow of Casson nanoliquid througha porous medium in the presence of chemicallyreactive activation energy?

正J. K. MADHUKESH1, G. K. RAMESH2, B. C. PRASANNAKUMARA1,S. A. SHEHZAD3,, F. M. ABBASI4     

MHD mixed convection–radiation interaction along a permeable surface immersed in a porous medium in the presence of Soret and Dufour’s Effects

This work is focused on the numerical modeling of steady, laminar, heat and mass transfer by MHD mixed convection from a semi-infinite, isothermal, vertical and permeable surface immersed in a uniform porous medium in the presence of thermal radiation and Dufour and Soret effects. A mixed convection parameter for the entire range of free-forced-mixed convection is employed and the governing equations are transformed into non-similar equations. These equations are solved numerically by an efficient, implicit, iterative, finite-difference scheme. The obtained results are checked against previously published work on special cases of the problem and are found to be in excellent agreement. A parametric study illustrating the influence of the thermal radiation coefficient, magnetic field, porous medium inertia parameter, concentration to thermal buoyancy ratio, and the Dufour and Soret numbers on the fluid velocity, temperature and concentration as well as the local Nusselt and the Sherwood numbers is conducted. The obtained results are shown graphically and the physical aspects of the problem are discussed.     

The influence of thermal radiation on hydromagnetic Darcy-Forchheimer mixed convection flow past a stretching sheet embedded in a porous medium

A boundary layer analysis is performed to study the influence of thermal radiation and buoyancy force on two-dimensional magnetohydrodynamic flow of an incompressible viscous and electrically conducting fluid over a vertical stretching sheet embedded in a porous medium in the presence of inertia effect. The governing system of partial differential equations is first transformed into system of ordinary differential equations using self-similarity transformation. A special form for magnetic field is chosen to obtain the similarity solution. The transformed boundary layer equations are solved numerically for some important values of the physical parameters. The present results are compared with the previously published papers and the results are found to be in excellent agreement. The important features of the flow, heat and mass transfer characteristics for different values of thermal radiation, porous permeability, magnetic field and buoyancy parameters are analyzed and discussed. The effects of various physical parameters on the skin friction coefficient, local Nusselt number and local Sherwood number are also presented. It is found that increase in the value of thermal radiation parameter R ₁ increases the skin friction coefficient and Sherwood number whereas reverse trend is seen for the local Nusselt number.     

Linear and nonlinear stability of double diffusive convection in a couple stress fluid–saturated porous layer

The linear and nonlinear stability of double diffusive convection in a layer of couple stress fluid–saturated porous medium is theoretically investigated in this work. Applying the linear stability theory, the criterion for the onset of steady and oscillatory convection is obtained. Emphasizing the presence of couple stresses, it is shown that their effect is to delay the onset of convection and oscillatory convection always occurs at a lower value of the Rayleigh number at which steady convection sets in. The nonlinear stability analysis is carried out by constructing a system of nonlinear autonomous ordinary differential equations using a truncated representation of Fourier series method and also employing modified perturbation theory with the help of self-adjoint operator technique. The results obtained from these two methods are found to complement each other. Besides, heat and mass transport are calculated in terms of Nusselt numbers. In addition, the transient behavior of Nusselt numbers is analyzed by solving the nonlinear system of ordinary differential equations numerically using the Runge–Kutta–Gill method. Streamlines, isotherms, and isohalines are also displayed.     

Linear stability of a fluid channel with a porous layer in the center简

We perform a Poiseuille flow in a channel linear stability analysis of a inserted with one porous layer in the centre, and focus mainly on the effect of porous filling ratio. The spectral collocation technique is adopted to solve the coupled linear stability problem. We investigate the effect of permeability, σ, with fixed porous filling ratio ψ = 1/3 and then the effect of change in porous filling ratio. As shown in the paper, with increasing σ, almost each eigenvalue on the upper left branch has two subbranches at ψ = 1/3. The channel flow with one porous layer inserted at its middle （ψ = 1/3） is more stable than the structure of two porous layers at upper and bottom walls with the same parameters. By decreasing the filling ratio ψ, the modes on the upper left branch are almost in pairs and move in opposite directions, especially one of the two unstable modes moves back to a stable mode, while the other becomes more instable. It is concluded that there are at most two unstable modes with decreasing filling ratio ψ. By analyzing the relation between ψ and the maximum imaginary part of the streamwise phase speed, Cimax, we find that increasing Re has a destabilizing effect and the effect is more obvious for small Re, where ψ a remarkable drop in Cimax can be observed. The most unstable mode is more sensitive at small filling ratio ψ, and decreasing ψ can not always increase the linear stability. There is a maximum value of Cimax which appears at a small porous filling ratio when Re is larger than 2 000. And the value of filling ratio 0 corresponding to the maximum value of Cimax in the most unstable state is increased with in- creasing Re. There is a critical value of porous filling ratio （= 0.24） for Re = 500; the structure will become stable as ψ grows to surpass the threshold of 0.24; When porous filling ratio ψ increases from 0.4 to 0.6, there is hardly any changes in the values of Cimax. We have also observed that the critical Reynolds number is especially sensitive for small ψ where the fastest drop is observed, and there may be a wide range in which the porous filling ratio has less effect on the stability （ψ ranges from 0.2 to 0.6 at σ = 0.002）. At larger permeability, σ, the critical Reynolds number tends to converge no matter what the value of porous filling ratio is.     

“Cold convection” in porous layers salted from above

A multicomponent fluid mixture saturating a porous rotating horizontal layer, heated from below and salted partly from below and partly from above, in the Darcy–Boussinesq scheme, is investigated. Conditions guaranteeing the “cold convection” i.e. the instability of the thermal conduction solution irrespective of the temperature gradient, are furnished.     

Simulation of conjugate radiation–forced convection heat transfer in a porous medium using the lattice Boltzmann method

In this paper, a lattice Boltzmann method is employed to simulate the conjugate radiation–forced convection heat transfer in a porous medium. The absorbing, emitting, and scattering phenomena are fully included in the model. The effects of different parameters of a silicon carbide porous medium including porosity, pore size, conduction–radiation ratio, extinction coefficient and kinematic viscosity ratio on the temperature and velocity distributions are investigated. The convergence times of modified and regular LBMs for this problem are 15 s and 94 s, respectively, indicating a considerable reduction in the solution time through using the modified LBM. Further, the thermal plume formed behind the porous cylinder elongates as the porosity and pore size increase. This result reveals that the thermal penetration of the porous cylinder increases with increasing the porosity and pore size. Finally, the mean temperature at the channel output increases by about 22% as the extinction coefficient of fluid increases in the range of 0–0.03.

     

Rayleigh–Taylor instability of immiscible fluids in porous media

The time development of an interface separating two immiscible fluids of different densities in heterogeneous two-dimensional porous media is studied. The governing equations are simplified with the help of approximate Green’s functions which allow computation of the shape of the interface directly without resolving the fluid flow in the entire domain. The new formulation is amenable to numerical approximation, and the reduction in dimension leads to a significant gain in efficiency in the numerical simulation of the interfacial dynamics. Several test cases are investigated, and the numerical solutions are compared to known exact solutions and experimental data.     

Steady crack growth in elastic–plastic fluid-saturated porous media

An asymptotic solution is obtained for stress and pore pressure fields near the tip of a crack steadily propagating in an elastic–plastic fluid-saturated porous material displaying linear isotropic hardening. Quasi-static crack growth is considered under plane strain and Mode I loading conditions. In particular, the effective stress is assumed to obey the Drucker–Prager yield condition with associative or non-associative flow-rule and linear isotropic hardening is adopted. Both permeable and impermeable crack faces are considered. As for the problem of crack propagation in poroelastic media, the behavior is asymptotically drained at the crack-tip. Plastic dilatancy is observed to have a strong effect on the distribution and intensity of pore water pressure and to increase its flux towards the crack-tip.     

Subgrid modeling of filtration in porous self‐similar media

Subgrid modeling of a filtration flow of a fluid in an inhomogeneous porous medium is considered. An expression for the effective permeability coefficient for the largescale component of the flow is derived using the scaleinvariance hypothesis. The model obtained is verified by numerical simulation of the complete problem.     

Chemical reaction and variable viscosity effects on flow and mass transfer of a non-Newtonian visco-elastic fluid past a stretching surface embedded in a porous medium

This work deals with the study of the boundary layer flow and mass transfer of a visco-elastic fluid immersed in a porous medium over a stretching surface in the presence of surface slip, chemical reaction and variable viscosity. The partial differential equations governing the flow have been transformed by similarity transformation into a system of coupled nonlinear ordinary differential equations which is solved numerically by means of the fourth order Runge-Kutta integration scheme coupled with the shooting technique. The effects of various involved interesting parameters on the velocity fields and concentration fields are shown graphically and investigated. In addition, tabulated results for the local skin-friction coefficient and the local Sherwood number are presented and discussed.     

Modeling of multiple noncatalytic gas–solid reactions in a moving bed of porous pellets based on finite volume method

A mathematical model is presented to simulate the multiple heterogeneous reactions with complex set of physicochemical and thermal phenomena in a moving bed of porous pellets. This model is based on both heat and mass transfer phenomena of gaseous species in a porous medium including chemical reactions at interfaces whose areas vary during the conversion. This model accounts for both the exothermic and endothermic reactions which can be equimolar or nonequimolar. Furthermore it considers simultaneously the reactions in the nonisothermal transient condition. A powerful technique based upon finite volume fully implicit approach has been implemented to solve the complicated governing equations numerically. The model has been validated by comparing with various experimental and analytical results in two cases: the single pellet scale as well as the counter current moving bed reactor.

Effects of bottom shear stresses on the wave-induced dynamic response in a porous seabed：PORO-WSSI （shear） model

When ocean waves propagate over the sea floor,dynamic wave pressures and bottom shear stresses exert on the surface of seabed.The bottom shear stresses provide a horizontal loading in the wave-seabed interaction system,while dynamic wave pressures provide a vertical loading in the system.However,the bottom shear stresses have been ignored in most previous studies in the past.In this study,the effects of the bottom shear stresses on the dynamic response in a seabed of finite thickness under wave loading will be examined,based on Biot’s dynamic poro-elastic theory.In the model,an "u-p" approximation will be adopted instead of quasi-static model that have been used in most previous studies.Numerical results indicate that the bottom shear stresses has certain influences on the wave-induced seabed dynamic response.Furthermore,wave and soil characteristics have considerable influences on the relative difference of seabed response between the previous model(without shear stresses) and the present model(with shear stresses).As shown in the parametric study,the relative differences between two models could up to 10% of p0,depending on the amplitude of bottom shear stresses.     

The râle of microscopic cross flow in idealized porous media

In this paper, results from a combined network/averaging study are presented. The emphasis is placed on understanding the flow phenomena, rather than predicting results for real porous media. Idealized porous media, consisting of networks of tubes, are used to interpret two of the terms in the averaged momentum equation. In particular, it is demonstrated that the pressure term accounts for microscopic cross flow, and that the magnitude of this term is proportional to the variation of the cross-sectional areas of the tubes in the macroscopic flow direction. For one-dimensional macroscopic flow in these idealized porous media, the agreement of network theory and averaging theory permeabilities depends on areosity (a term related to the area open to flow in a direction) remaining constant in the macroscopic flow direction; it may vary in other directions.