Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Modeling the co-transport of viruses and colloids in unsaturated porous media

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations,are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Cluysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWL A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AIM, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWL The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SW!, maximum adsorption capacity of the AWI for viruses and inlet virus concentration. (C) 2015 Elsevier B.V. All rights reserved.

Non-darcy double-diffusive mixed convection from heated vertical and horizontal plates in saturated porous media

The non-darcy mixed convection flows from heated vertical and horizontal plates in saturated porous media have been considered using boundary layer approximations. The flows are considered to be driven by multiple buoyancy forces. The similarity solutions for both vertical and horizontal plates have been obtained. The governing equations have been solved numerically using a shooting method. The heat transfer, mass transfer and skin friction are reduced due to inertial forces. Also, they increase with the buoyancy parameter for aiding flow and decrease for the opposing flow. For aiding flow, the heat and mass transfer coefficients are found to approach asymptotically the forced or free convection values as the buoyancy parameter approaches zero or infinity.

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

We present the linear stability analysis of horizontal Poiseuille flow in a fluid overlying a porous medium with anisotropic and inhomogeneous permeability. The generalized Darcy model is used to describe the flow in the porous medium with the Beavers-Joseph condition at the interface of the two layers and the eigenvalue problem is solved numerically. The effect of major system parameters on the stability characteristics is addressed in detail. It is shown that the anisotropic and inhomogeneous modulation of the permeability of the underlying porous layer provides an effective means for passive control of the flow stability.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Unsteady free convection flow in the stagnation-point region of a three-dimensional body

The unsteady free convection flow in the stagnation-point region of a heated three-dimensional body placed in an ambient fluid is studied under boundary layer approximations. We have considered the case where there is an initial steady state that is perturbed by a step-change in the wall temperature. The non-linear coupled partial differential equations governing the free convection flow are solved numerically using a finite difference scheme. The presented results show the temporal development of the momentum and thermal boundary layer characteristics.

The effects of slip velocity at a membrane surface on blood flow in the microcirculation

Closed-form solutions are presented for blood flow in the microcirculation by taking into account the influence of slip velocity at the membrane surface. In this study, the convective inertia force is neglected in comparison with that of blood viscosity on the basis of the smallness of the Reynolds number of the flow in microcirculation. The permeability property of the blood vessel is based on the well known Starling's hypothesis [11]. The effects of slip coefficient on the velocity and pressure fields are clearly depicted.

A couple stress model of blood flow in the microcirculation

A simple mathematical model depicting blood flow in the capillary is developed with an emphasis on the permeability property of the blood vessel based on Starling's hypothesis. In this study the effect of inertia has been neglected in comparison with the viscosity on the basis of the smallness of the Reynolds number of the flow in the capillary. The capillary blood vessel is approximated by a circular cylindrical tube with a permeable wall. The blood is represented by a couple stress fluid. With such an ideal model the velocity and pressure fields are determined. It is shown that an increase in the couple stress parameter increases the resistance to the flow and thereby decreases the volume rate flow. A comparison of the results with those of the Newtonian case has also been made.

Compressible axially symmetric laminar boundary layer flow in the stagnation region of a blunt body in the presence of magnetic field

In this paper, the steady laminar viscous hypersonic flow of an electrically conducting fluid in the region of the stagnation point of an insulating blunt body in the presence of a radial magnetic field is studied by similarity solution approach, taking into account the variation of the product of density and viscosity across the boundary layer. The two coupled non-linear ordinary differential equations are solved simultaneously using Runge-Kutta-Gill method. It has been found that the effect of the variation of the product of density and viscosity on skin friction coefficient and Nusselt number is appreciable. The skin friction coefficient increases but Nusselt number decreases as the magnetic field or the total enthalpy at the wall increases

Unsteady compressible boundary layer flow in the stagnation region of a sphere with a magnetic field

Abstract: An analysis is performed to study the unsteady compressible laminar boundary layer flow in the forward stagnation-point region of a sphere with a magnetic field applied normal, to the surface. We have considered the case where there is an initial steady state that is perturbed by the step change in the total enthalpy at the wall. The nonlinear coupled parabolic partial differential equations governing the flow and heat transfer have been solved numerically using a finite-difference scheme. The numerical results are presented, which show the temporal development of the boundary layer. The magnetic field in the presence of variable electrical conductivity causes an overshoot in the velocity profile. Also, when the total enthalpy at the wall is suddenly increased, there is a change in the direction of transfer of heat in a small interval of time.

Effect of capillarity and heterogeneity in the numerical modelling of multiphase flow of fluids in unsaturated porous medium

The multiphase flow of fluids in the unsaturated porous medium is considered as a three phase flow of water, NAPL, and air simultaneously in the porous medium. The adaptive solution fully implicit modified sequential method is used for the numerical modelling. The effect of capillarity and heterogeneity effect at the interface between the media is studied and it is observed that the interface criteria has to be taken into account for the correct prediction of NAPL migration especially in heterogeneous media. The modified Newton Raphson method is used for the linearization and Hestines and Steifel Conjugate Gradient method is used as the solver.

Pulsatile blood flow in a rigid pulmonary alveolar sheet with porous walls

This paper deals with the pulsatile blood flow in the lung alveolar sheets by idealizing each of them as a channel covered by porous media. As the blood flow in the lung is of low Reynolds number, a creeping flow is assumed in the channel. The analytical and numerical results for the velocity and pressure distribution in the porous medium are presented. The effect of an imposed slip condition is also studied. Comparisons with the corresponding results for the steady-state case are made at the end.

TRANSIENT NATURAL CONVECTION FLOW IN A RECTANGULAR CAVITY FILLED WITH A POROUS MATERIAL WITH LOCALIZED HEATING FROM BELOW AND THERMAL STRATIFICATION

The transient natural convection flow with thermal stratification in a rectangular cavity filled with fluid saturated porous medium obeying Darcy's law has been studied. Prior to the time t* = 0, the flow in the cavity is assumed to be motionless and all four walls of the cavity are at the same constant temperature. At time t* = 0, the temperatures of the vertical walls are suddenly increased which vary linearly with the distance y and at the same time on the bottom wall an isothermal heat source is placed centrally. This sudden change in the wall temperatures gives rise to unsteadiness in the problem. The horizontal temperature difference induces and sustains a buoyancy driven flow in the cavity which is then controlled by the vertical temperature difference. The partial differential equations governing the transient natural convection flow have been solved numerically. The local and average Nusselt numbers decrease rapidly in a small time interval after the start of the impulsive change in the wall temperatures and the steady state is reached quickly. The time required to reach the steady state depends on the Rayleigh number and the thermal stratification parameter.

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

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 Narasimhan2).