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.

Influence of the Methane-Zeolite A Interaction Potential on the Concentration Dependence of Self-diffusivity

Studies on the diffusion of methane in a zeolite structure type LTA (as per IZA nomenclature) have indicated that different types of methane zeolite potentials exist in the literature in which methane is treated within the united-atom model. One set of potentials, referred to as model A, has a methane oxygen diameter of 3.14 angstrom, while another set of potential parameters, model B, employs a larger value of 3.46 angstrom. Fritzsche and co-workers (1993) have shown that these two potentials lead to two distinctly different energetic barriers for the passage of methane through the eight-ring window in the cation-free form of zeolite A. Here, we compute the variation of the self-diffusivity (D) with loading (c) for these two types of potentials and show that this slight variation in the diameter changes the concentration dependence qualitatively: thus, D decreases monotonically with c for model A, while D increases and goes through a maximum before finally decreasing for model B. This effect and the surprising congruence of the diffusion coefficients for both models at high loadings is examined in detail at the molecular level. Simulations for different temperatures reveal the Arrhenius behaviour of the self-diffusion coefficient. The apparent activation energy is found to vary with the loading. We conclude that beside the cage-to-cage jumps, which are essential for the migration of the guest molecules, at high concentrations migration within the cage and guest guest interactions with other molecules become increasingly dominant influences on the diffusion coefficient and make the guest zeolite interaction less important for both model A and model B.

Review on atomization and sprays of biofuels for IC engine applications

Ever increasing energy requirements, environmental concerns and energy security needs are strongly influencing engine researchers to consider renewable biofuels as alternatives to fossil fuels. Spray process being important in IC engine combustion, existing literature on various biofuel sprays is reviewed and summarized. Both experimental and computational research findings are reviewed in a detailed manner for compression ignition (CI) engine sprays and briefly for spark ignition (SI) engine sprays. The physics of basic atomization process of sprays from various injectors is included to highlight the most recent research findings followed by discussion highlighting the effect of physico-chemical properties on spray atomization for both biofuels and fossil fuels. Biodiesel sprays are found to penetrate faster and haw narrow spray plume angle and larger droplet sizes compared to diesel. Results of analytical and computational models are shown to be useful in shedding light on the actual process of atomization. However, further studies on understanding primary atomization and the effect of fuel properties on primary atomization are required. As far as secondary atomization is concerned, changes in regimes are observed to occur at higher air-jet velocities for biodiesel compared to those of diesel. Evaporating sprays revealed that the liquid length is longer for biodiesel. Pure plant oil sprays with potential use in CI engines may require alternative injector technology due to slower breakup as compared to diesel. Application of ethanol to gasoline engines may be feasible without any modifications to port fuel injection (PFI) engines. More studies are required on the application of alternative fuels to high pressure sprays used in Gasoline Direct Injection (GDI) engines.

Vertical uplift resistance of two closely spaced horizontal strip anchors embedded in cohesive-frictional weightless medium

The vertical uplift resistance of two closely spaced horizontal strip plate anchors has been investigated by using lower and upper bound theorems of the limit analysis in combination with finite elements and linear optimization. The interference effect on uplift resistance of the two anchors is evaluated in terms of a nondimensional efficiency factor (eta(c)). The variation of eta(c) with changes in the clear spacing (S) between the two anchors has been established for different combinations of embedment ratio (H/B) and angle of internal friction of the soil (phi). An interference of the anchors leads to a continuous reduction in uplift resistance with a decrease in spacing between the anchors. The uplift resistance becomes a minimum when the two anchors are placed next to each other without any gap. The critical spacing (S-cr) between the two anchors required to eliminate the interference effect increases with an increase in the values of both H/B and phi. The value of S-cr was found to lie approximately in the range 0.65B-1.5B with H/B = 1 and 11B-14B with H/B = 7 for phi varying from 0 degrees to 30 degrees.

Pullout capacity of inclined plate anchors embedded in sand

The pullout capacity of an inclined strip plate anchor embedded in sand has been determined by using the lower bound theorem of the limit analysis in combination with finite elements and linear optimization. The numerical results in the form of pullout factors have been presented by changing gradually the inclination of the plate from horizontal to vertical. The pullout resistance increases significantly with an increase in the horizontal inclination (theta) of the plate especially for theta > 30 degrees. The effect of the anchor plate-soil interface friction angle (delta) on the pullout resistance becomes extensive for a vertical anchor but remains insignificant for a horizontal anchor. The development of the failure zone around the anchor plates was also studied by varying theta and delta. The results from the analysis match well with the theoretical and experimental results reported in literature.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

Bearing capacity factors for a conical footing using lower- and upper-bound finite elements limit analysis

Bearing capacity factors, N-c, N-q, and N-gamma, for a conical footing are determined by using the lower and upper bound axisymmetric formulation of the limit analysis in combination with finite elements and optimization. These factors are obtained in a bound form for a wide range of the values of cone apex angle (beta) and phi with delta = 0, 0.5 phi, and phi. The bearing capacity factors for a perfectly rough (delta = phi) conical footing generally increase with a decrease in beta. On the contrary, for delta = 0 degrees, the factors N-c and N-q reduce gradually with a decrease in beta. For delta = 0 degrees, the factor N-gamma for phi >= 35 degrees becomes a minimum for beta approximate to 90 degrees. For delta = 0 degrees, N-gamma for phi <= 30 degrees, as in the case of delta = phi, generally reduces with an increase in beta. The failure and nodal velocity patterns are also examined. The results compare well with different numerical solutions and centrifuge tests' data available from the literature.