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.

Incipient motion design of sand bed channels affected by bed suction

Seepage through sand bed channels in a downward direction (suction) reduces the stability of particles and initiates the sand movement. Incipient motion of sand bed channel with seepage cannot be designed by using the conventional approach. Metamodeling techniques, which employ a non-linear pattern analysis between input and output parameters and solely based on the experimental observations, can be used to model such phenomena. Traditional approach to find non-dimensional parameters has not been used in the present work. Parameters, which can influence the incipient motion with seepage, have been identified and non-dimensionalized in the present work. Non-dimensional stream power concept has been used to describe the process. By using these non-dimensional parameters; present work describes a radial basis function (RBF) metamodel for prediction of incipient motion condition affected by seepage. The coefficient of determination, R-2 of the model is 0.99. Thus, it can be said that model predicts the phenomena very well. With the help of the metamodel, design curves have been presented for designing the sand bed channel when it is affected by seepage. (C) 2010 Elsevier B.V. All rights reserved.

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.

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.

Chromium-manganese iron alloy system design cast in metal and sand moulds for erosion resistance: a positron lifetime study

Erosion characteristics of high chromium (Cr, 16-19%) alloy cast iron with 5% and 10% manganese (Mn) prepared in metal and sand moulds through induction melting are investigated using jet erosion test setup in both as-cast and heat-treated conditions. The samples were characterised for hardness and microstructural properties. A new and novel non-destructive evaluation technique namely positron lifetime spectroscopy has also been used for the first time to characterise the microstructure of the material in terms of defects and their concentration. We found that the hardness decreases irrespective of the sample condition when the mould type is changed from metal to sand, On the other hand, the erosion volume loss shows an increasing trend. Since the macroscopic properties have a bearing on the microstructure, good credence is obtained from the microstructural features as seen from light and scanning electron micrographs. Faster cooling in the metal mould yielded fine carbide precipitation on the surface. The defect size and their concentration derived from positron method are higher for sand mould compared to metal mould. Lower erosion loss corresponds to smaller size defects in metal mould are the results of quicker heat transfer in the metal mould compared to the sand mould. Heat treatment effects are clearly seen as the reduced concentration of defects and spherodisation of carbides points to this. The erosion loss with respect to the defects size and concentration correlate very well.

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.

Variation of stream power with seepage in sand-bed channels

Downward seepage (suction) increases the mobility of the channel. In this study, experimental investigations were carried out to analyse the suction effect on stream power along the downstream side of the flume. It was observed that stream power has a major influence on the stability and mobility of the bed particles, due to suction. Stream power is found to be greater at the upstream side and lower at the downstream side. This reduces the increment in the mobility of the sand particles due to suction at the downstream side. Thus, there is more erosion at the upstream side than the downstream side. It was also found that the amount of deposition of sand particles at the downstream side, because of the high stream power at the upstream side, is greater than the amount of erosion of sand particles from the downstream side.

Effect of anchor width on pullout capacity of strip anchors in sand

By incorporating the variation of peak soil friction angle (phi) with mean principal stress (sigma(m)), the effect of anchor width (B) on vertical uplift resistance of a strip anchor plate has been examined. The anchor was embedded horizontally in a granular medium. The analysis was performed using lower bound finite element limit analysis and linear programming. An iterative procedure, proposed recently by the authors, was implemented to incorporate the variation of phi with sigma(m). It is noted that for a given embedment ratio, with a decrease in anchor width (B), (i) the uplift factor (F-gamma) increases continuously and (ii) the average ultimate uplift pressure (q(u)) decreases quite significantly. The scale effect becomes more pronounced at greater embedment ratios.

Seepage effects on sand-bed channels

ExperimentS were conducted to find the effects of seepage on flow over a sand bed in a straight rectangular flume under two conditions: (1) When the channel bed is plane, horizontal, and nontransporting; and (2) when the bed is transporting at a constant sediment concentration. Effects of both injection and suction, caused by seepage flow into and out of the channel bed, are studied for condition 1; and only suction effects are studied for condition 2, Three sands, sizes 0.34 mm, 0.53 mm, and 0.80 mm, are used in the study. It is found that seepage can cause an increase or decrease in the bed shear stress relative to no seepage for the two conditions. The change in bed shear stress depends on the relative magnitudes of the bed shear stress and the critical shear stress of particles under the no-seepage condition, sediment concentration, and the seepage rate. Quantitative relationships giving the ratio of bed shear stresses with and without seepage are presented for both conditions of the bed. A procedure to estimate the changes in bed shear stress, friction factor, Manning's n, and stream power due to seepage for known initial conditions of the channel and the amount of applied seepage is presented

Ionic conductivity, mechanical strength and Li-ion battery performance of mono-functional and bi-functional (''Janus'') ``soggy sand'' electrolytes

Influence of dispersion of uniformly sized mono-functional and bi-functional (''Janus'') particles on ionic conductivity of novel ``soggy sand'' electrolytes and its implications on mechanical strength and lithium-ion battery performance are discussed here.

Mixed convection heat transfer from the bottom tip of a cylinder spinning about a vertical axis in a saturated porous medium

The present work is a numerical study of heat transfer characteristics from the bottom tip of a cylinder spinning about a vertical axis in an infinitely saturated porous medium. The problem is axisymmetric. The non-dimensionalized governing equations are solved using the SIMPLER algorithm on a staggered grid. The influence of rotational Reynolds numbers and Darcy numbers on the heat transfer for a Grashof number of 104 and Prandtl number of 7.0 is studied. It is found that for very high Darcy numbers, over a wide range of rotational Reynolds numbers, the heat transfer takes place mainly due to conduction. The convective heat transfer takes place for lower Darcy numbers and for higher rotational Reynolds numbers. Moreover, there is a rapid increase in the overall Nusselt number below a certain Darcy number with increase in the rotational Reynolds numbers. The effect of the Darcy number and the rotational Reynolds number on the heat transfer and fluid flow in the porous medium is depicted in the form of streamline and isotherm plots. The variation of the overall Nusselt number with respect to the Darcy number for various rotational Reynolds numbers is plotted. The variation of the local Nusselt number with respect to the radial coordinate at the heated tip of the vertical cylinder is plotted for various Darcy and rotational Reynolds numbers.

Geometry of sand-bed channels with seepage

Literature suggests that apart from a significant percentage of water loss from seepage, downward seepage from the channel causes an increase in the mobility of channel bed materials and thus changes the channel stability. Consequently, regime conditions (which provide the relationship among hydraulic parameters) should also be revised by incorporating downward seepage as an additional parameter. In the present work, regime conditions for the prediction of channel geometry in alluvial channels affected by downward seepage have been formulated on the basis of experimental observations. (C) 2011 Elsevier B.V. All rights reserved.

Saturated liquid viscosity correlations for alternative refrigerants

This article presents dimensionless equations for the temperature dependence of the saturated liquid viscosity of R32, R123, R124, R125, R134a, R141b, and R152a valid over a temperature range of engineering interest. The correlation has the form Phi(D)(n)=A+BTD where Phi(D) is the dimensionless fluidity (1/eta(D)) and T-D is a dimensionless temperature. n, A, and B are evaluated for each of the above refrigerants based on a least-squares fit to experimental data. This equation is found to provide an improved fit over those existing in the literature up to T-D=0.8.