Numerical simulation for the 3D seepage flow with fractional derivatives in porous media

In this paper, the numerical simulation of the 3D seepage flow with fractional derivatives in porous media is considered under two special cases: non-continued seepage flow in uniform media (NCSFUM) and continued seepage flow in non-uniform media (CSF-NUM). A fractional alternating direction implicit scheme (FADIS) for the NCSF-UM and a modified Douglas scheme (MDS) for the CSF-NUM are proposed. The stability, consistency and convergence of both FADIS and MDS in a bounded domain are discussed. A method for improving the speed of convergence by Richardson extrapolation for the MDS is also presented. Finally, numerical results are presented to support our theoretical analysis.

Modified alternating direction methods for solving a two-dimensional non-continuous seepage flow with fractional derivatives

In this paper, a two-dimensional non-continuous seepage flow with fractional derivatives (2D-NCSF-FD) in uniform media is considered, which has modified the well known Darcy law. Using the relationship between Riemann-Liouville and Grunwald-Letnikov fractional derivatives, two modified alternating direction methods: a modified alternating direction implicit Euler method and a modified Peaceman-Rachford method, are proposed for solving the 2D-NCSF-FD in uniform media. The stability and consistency, thus convergence of the two methods in a bounded domain are discussed. Finally, numerical results are given.

Influence of spatial variability of permeability property on steady state seepage flow and slope stability analysis

In recent years, spatial variability modeling of soil parameters using random field theory has gained distinct importance in geotechnical analysis. In the present Study, commercially available finite difference numerical code FLAC 5.0 is used for modeling the permeability parameter as spatially correlated log-normally distributed random variable and its influence on the steady state seepage flow and on the slope stability analysis are studied. Considering the case of a 5.0 m high cohesive-frictional soil slope of 30 degrees, a range of coefficients of variation (CoV%) from 60 to 90% in the permeability Values, and taking different values of correlation distance in the range of 0.5-15 m, parametric studies, using Monte Carlo simulations, are performed to study the following three aspects, i.e., (i) effect ostochastic soil permeability on the statistics of seepage flow in comparison to the analytic (Dupuit's) solution available for the uniformly constant permeability property; (ii) strain and deformation pattern, and (iii) stability of the given slope assessed in terms of factor of safety (FS). The results obtained in this study are useful to understand the role of permeability variations in slope stability analysis under different slope conditions and material properties. (C) 2009 Elsevier B.V. All rights reserved.

Design of plane sand-bed channels affected by seepage

The importance of seepage in the design of channels is discussed. Experimental investigations reveal that seepage, either in the downward direction (suction) or in the upward direction (injection), can significantly change the resistance as well as the mobility of the sand-bed particles. A resistance equation relating 'particle Reynolds number' and 'shear Reynolds number' under seepage conditions is developed for plane sediment beds. Finally, a detailed design procedure of the plane sediment beds affected by seepage is presented.

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.

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

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.

Simplified seepage consolidation test for soft sediments

The consistency of very soft sediments prevents the conventional oedometer test from being applied to study their compressibility and permeability characteristics. The hydraulic consolidation test in existence requires sophisticated instrumentation and testing procedures. The present paper proposes a seepage-force-induced consolidation testing procedure for studying the compressibility and permeability behavior of soft sediments at low effective stress levels. The good agreement that has been observed between the results obtained from the proposed method and the conventional oedometer test at overlapping effective stress levels indicates that the proposed method can be used to study the compressibility and permeability characteristics of soft sediments at low effective stress levels satisfactorily.

Stability and mobility of sand-bed channels affected by seepage

Seepage effects on the stability, mobility, and incipient motion of sand-bed particles are experimentally investigated. Seepage through a sand bed in a downward direction (suction) reduces the stability of particles, and it can even initiate their movement. The bed erosion is increased with the increased rates of suction. Whereas the seepage in an upward direction (injection) increases the stability of bed particles, it does not aid initiating their movement. The rate of bed erosion is reduced or even stopped by the increased infection rates. Hydrodynamic conditions leading to the so-called "pseudoincipient motion'' with suction (for the initiation of particles movement that are otherwise at rest under no-seepage conditions), and with injection (for only arresting the particles movement that are otherwise moving initially) are evaluated. The conventional Shields curve cannot be used to predict such pseudoincipient motion conditions with seepage. The concepts thus developed are useful for a better understanding of the sediment transport mechanics and in the design of stable alluvial channels affected by seepage.

Modeling of seepage flow through layered solis

A two-dimensional finite difference model, which solves mixed type of Richards' equation, whose non-linearity is dealt with modified Picard's iteration and strongly implicit procedure to solve the resulting equations, is presented. Modeling of seepage flow through heterogeneous soils, which is common in the field is addressed in the present study. The present model can be applied to both unsaturated and saturated soils and can handle very dry initial condition and steep wetting fronts. The model is validated by comparing experimental results reported in the literature. Newness of this two dimensional model is its application on layered soils with transient seepage face development, which has not been reported in the literature. Application of the two dimensional model for studying unconfined drainage due to sudden drop of water table at seepage face in layered soils is demonstrated. In the present work different sizes of rectangular flow domain with different types of layering are chosen. Sensitivity of seepage height due to problem dimension of layered system is studied. The effect of aspect ratio on seepage face development in case of the flow through layered soil media is demonstrated. The model is also applied to random heterogeneous soils in which the randomness of the model parameters is generated using the turning band technique. The results are discussed in terms of phreatic surface and seepage height development and also flux across the seepage face. Such accurate modeling of seepage face development and quantification of flux moving across the seepage face becomes important while modeling transport problems in variably saturated media.

Bearing Capacity of Foundations with Inclined Groundwater Seepage

The ultimate bearing capacity of strip foundations in the presence of inclined groundwater flow, considering both upward and downward flow directions, has been determined by using the lower bound finite-element limit analysis. A numerical solution has been generated for both smooth and rough footings placed on frictional soils. A correction factor (f gamma), which needs to be multiplied with the N gamma-term, has been computed to account for groundwater seepage. The variation of f gamma has been obtained as a function of the hydraulic gradient (i) for various inclinations of groundwater flow. For a given magnitude of i, there exists a certain critical inclination of the flow for which the value of f gamma is minimized. With an upward flow, for all flow inclinations, the magnitude of f gamma always reduces with an increase in the value of i. An example has also been provided to illustrate the application of the obtained results when designing foundations in the presence of groundwater seepage.

Effect of Groundwater Seepage on Uplift Resistance of Buried Pipelines

The uplift resistance of pipelines buried in sands, in the presence of inclined groundwater flow, considering both upward and downward flow directions, has been determined by using the lower bound finite elements limit analysis in conjunction with nonlinear optimization. A correction factor (f (gamma) ), which needs to be multiplied with the uplift factor (F (gamma) ), has been computed to account for groundwater seepage. The variation of f (gamma) has been obtained as a function of i(gamma (w) /gamma (sub) ) for different horizontal inclinations (theta) of groundwater flow; where i = absolute magnitude of hydraulic gradient along the direction of flow, gamma (w) is the unit weight of water and gamma (sub) is the submerged unit weight of soil mass. For a given magnitude of i, there exists a certain critical value of theta for which the magnitude of f (gamma) becomes the minimum. An example has also been presented to illustrate the application of the results obtained for designing pipelines in presence of groundwater seepage.

Seepage beneath Hoover Dike southern shore of Lake Okeechobee, Florida

Future water needs in southern Florida call for an increase in the storage capacity of Lake Okeechobee. Seepage from the lake is expected to increase as a result of raising the lake level. Data concerning the occurrence and amounts of seepage are needed for the design and operation of flood-control works which will remove excess water from the rich agricultural lands along the southern shore. Intensive studies at five sites along the southern shore of Lake Okeechobee between the Caloosahatchee Canal and the St. Lucie Canal indicate that seepage occurs chiefly through beds of shell and limestone which underlie the Hoover Dike at shallow depth. Seepage rates at the five sites range from about 0.1 to 0.9 cfs per mile per foot of head across the dike. Seepage beneath the 50-mile length of dike should increase from about 22 to 50 cfs if the average stage of the lake is raised from 14 to 16.5 feet. Seepage is greatest between Moore Haven and Clewiston, where deep borrows have been excavated on the landward and lakeward sides of the dike. Most of the seepage from the lake can be controlled by properly spaced toe ditches which would intercept the seepage and return it to the lake. (PDF contains 108 pages.)