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.

Nonsimilar boundary layers for non-Darcy mixed convection flow about a horizontal surface in a saturated porous medium

The nonsimilar non-Darcy mixed convection flow about a heated horizontal surface in a saturated porous medium has been studied when the surface temperature is a power function of distance (Tw = T∞ ± Axλ). The analysis is performed for the cases of parallel and stagnation flows with favourable induced pressure gradient. The partial differential equations governing the flow have been solved numerically using the Keller box method. The heat transfer is enhanced due to the buoyancy parameter and wall temperature, but the non-Darcy parameter reduces it. For non-Darcy flow, the similarity solution exists only for the case of parallel flow.

Engineering behavior of uplifted Smectite-rich Cochin and Mangalore marine clays

The present study aims to assess whether the smectite-rich Cochin and Mangalore clays, which were deposited in a marine medium and subsequently uplifted, exhibit consistency limits response typical of expanding lattice or nonexpanding (fixed) lattice-type clays on artificially changing the chemical environment. The chemical and engineering behaviors of Cochin and Mangalore marine clays are also compared with those of the smectite-rich Ariake Bay marine clay from Japan. Although Cochin, Mangalore, and Ariake clays contain comparable amounts of smectite (32-45%), Ariake clay exhibits lower consistency limits and much higher ranges of liquidity indices than the Indian marine clays. The lower consistency limits of the Ariake clay are attributed to the absence of well-developed, long-range, interparticle forces associated with the clay. Also, Ariake clay exhibits a significantly large (48-714 times) decrease in undrained strength on remolding in comparison to Cochin and Mangalore clays (sensitivity ranges between 1 and 4). A preponderance of long-range, interparticle forces reflected in the high consistency limits of Cochin and Mangalore clays (wL range from 75 to 180%) combined with low natural water contents yield low liquidity indices (typically <1) and high, remolded, undrained strengths and are considered to be responsible for the low sensitivity of the Indian marine clays.

Critical Reappraisal of Colloidal Activity of Clays

Soil properties and their behavior, apart from stress history, are influenced markedly by physicochemical characteristics of the constituent clay and nonclay minerals and their relative proportions. Atterberg limits and Skempton’s colloidal activity, which are simple quantitative parameters, reflect the composite effects of the soil constituents and their interactions with pore fluid. Micromechanistic interpretations of these parameters have been provided in this paper. It has been shown that, in general, the liquid limit of fine-grained soils reflects the physicochemical potential and that each of the factors of Skempton’s colloidal activity are interdependent. It has been illustrated that property correlations with colloidal activity, as well as with Atterberg limits, result in involved interrelationships due to the interdependence of the parameters.

Critical Reappraisal of Colloidal Activity of Clays

Soil properties and their behavior, apart from stress history, are influence markedly by physicochemical characteristics of the constituent clay and nonclay minerals and their relative proportions. Atterberg limits and Skempton's colloidal activity, which are simple quantitative parameters, reflect the composite effects of the soil constituents and their interactions with pore fluid. Micromechanistic interpretations of these parameters have been provided in this paper. It has been shown that, in general, the liquid limit of fine-grained soils reflects the physicochemical potential and that each of the factors of Skempton's colloidal activity are interdependent. It has been illustrated that property correlations with colloidal activity, as well as with Atterberg limits, result in involved interrelationships due to the interdependence of the parameters.

A reexamination of the permeability index of clays

The permeability index Ck, similar to the compression index, is the slope of the void ratio – coefficient of permeability relationship. Literature shows that, in general, for sensitive clays it can be related to initial void ratio by Ck = 0.5e0. The possibility of obtaining such a relationship for Cochin marine clays in terms of liquid limit void ratio is indicated in this paper. Analysis of permeability behaviour of Cochin marine clays and the test results available in published literature using generalized state parameter approach show that, in principle, these forms of equations for the permeability index are tenable, even though they were obtained based on experimental observation alone.

Non-Darcy mixed convection flow with thermal dispersion on a vertical cylinder in a saturated porous medium

The non-Darcy mixed convection flow on a vertical cylinder embedded in a saturated porous medium has been studied taking into account the effect of thermal dispersion. Both forced flow and buoyancy force dominated cases with constant wall temperature condition have been considered. The governing partial differential equations have been solved numerically using the Keller box method. The results are presented for the buoyancy parameter which cover the entire regime of mixed convection flow ranging from pure forced convection to pure free convection. The effect of thermal dispersion is found to be more pronounced on the heat transfer than on the skin friction and it enhances the heat transfer but reduces the skin friction.

Tropical Clays .1. Index Properties And Microstructural Aspects

Soils showing changes in plasticity characteristics upon driving form an important group in tropical soils. These changes are attributed to the grouping of particles into aggregates either due to mineralogy or presence of cementing agents and/or pore fluid characteristics. These changes are found to be permanent. In this paper, the effect of these changes leading to changes in index properties is discussed. The coefficient of permeability has been found to be comparable at liquid limit water content for different soils of varying liquid limit values. Permeability is an indirect reflection of microstructure and indicates the flow rate, which depends upon pore geometry. Other mechanical properties like compressibility and shear strength also depend upon pore geometry. These microstructural aspects of liquid limit as a reference state for the analysis of engineering behavior of tropical soils are examined in detail.

Generalized State Parameter for Partly Saturated Soils

Studies on compressibility and shear strength aspects are the concern of many investigators concerned with partly saturated soils. In soil engineering connected with partly saturated soils, there are no approaches connecting soil states and stress conditions. The present investigation is essentially a step in this direction. A generalized state parameter, identified with regard to material states is shown to be related to the compressibility and shear strength. The involved parameters are simple and normally determined in routine investigations. The advantage of this approach is that changes in soil states due to external stress conditions and the associated changes in strength can be examined particularly when different types of soils are involved.