Open access
Technical Papers
Aug 1, 2014

Study of Flow Transitions during Air Sparging Using the Geotechnical Centrifuge

Publication: Journal of Environmental Engineering
Volume 141, Issue 1

Abstract

Air sparging (AS) is an in-situ groundwater remediation technique for volatile organic compounds (VOCs) in contaminated sites. Knowledge of effect of the injected air on air flow patterns and extent of the zone known as zone of influence (ZOI) is crucial for design of in situ air sparging. The two main objectives of this research are to study transition of air flow patterns and characteristics of ZOI under different air injection pressures during centrifugal tests and to find the scaling factors for centrifugal tests under different gravitational levels (N times gravitational force). In this research, the modeling of models test was first performed to validate centrifugal technique for simulating air sparging. Then, a set of two-dimensional (2D) centrifugal modeling tests were performed under various gravitational levels to study the transition of air flow patterns and characteristics of ZOI during air sparging. The centrifugal model tests showed that the observed 2D ZOI was truncated-cone under various gravitational levels. It was also observed that flow pattern for the same soil transits from channel flow under lower gravitational levels to bubble flow under higher gravitational level. If there are same air flow regimes during air sparging tests under varied gravitational levels, then the scaling factor for air sparging is N, and a formula was proposed to check the flow regime. In this research, a formula was proposed to explain the transition of air flow patterns during air sparging and determine the boundary of flow regime transition based on grain size and unit-sparging rate.

Introduction

Contaminated soils and groundwater is currently one of the major environmental problems facing the earth. There are many technologies to remediate contaminants in soil and groundwater (U.S. EPA 1995). Air sparging is an in-situ remediation technology where air is injected under pressure into a well installed in the saturated zone beneath the contaminated zone. Air injected below the water table during air sparging volatilizes contaminants dissolved in groundwater. Volatilized contaminants removed during air sparging migrate upward into the vadose zone and are removed using a soil vapor extraction technique. In addition to air stripping, air sparging also promotes biodegradation of organic contaminants by increasing oxygen concentrations in the subsurface. Air sparging systems must be designed with adequate air flow rates and pressures for the effective removal of contaminants from the sites.
The nature and extent of the air pathways determine the zone of influence (ZOI) during air sparging. Two-phase flow in porous media is one of the theoretical models for describing the air sparging process. McCray and Falta (1997) used a multiphase flow model (T2VOC) based on the integrated finite difference method to simulate the air sparging process. McCray (2000) presented a comprehensive review on numerical modeling of air sparging. Several numerical models were employed for air sparging modeling and compared with the experiments data (Mei and Cheng 2002; Jang and Aral 2009; Wang and Hu 2012). Gao et al. (2011, 2012, 2013a, b) proposed a microscopic model to describe air migration behavior during air sparging.
Several laboratory investigations have been conducted to determine the contributing factors and efficiency of air sparging (Ji et al. 1993; Nyer and Suthersan 1993; Johnson et al. 1993; Elder and Benson 1999; Reddy and Adams 2001; Adams and Reddy 2003; Hu and Liu 2008). The efficiency of an in-situ sparging system is controlled by the extent of contact between injected air and contaminated soil or pore fluid. However, most of the experiments were conducted under 1g using scaled models, and the in-situ stress state was not simulated. Also only small injection pressure can be applied with 1g models to avoid soil failure. The geotechnical centrifugal modeling can be used to simulate the movement and remediation of contaminants under in-situ stress state in both soil mass and interstitial fluid of porous media (Culligan-Hensley and Savvidou 1995). The scaled model is typically constructed in the laboratory and then loaded onto the end of the centrifuge. The purpose of spinning the models on the geotechnical centrifuge is to increase the centrifugal force on the model so that stresses in the model are made equal to stresses in the prototype. A 1/N scale model tested at N times the gravity because of centrifugal acceleration experiences stress conditions identical to that in the prototype. Soils are compressible, and the hydrogeological properties are affected by in-situ stress. Hence the behavior of a prototype subjected to the earth’s gravity g, can be reproduced in the scaled model subjected to a centrifugal force. The centrifuge models provide the unique advantage of testing field soils under exact in-situ stresses and fluid pressures, and can simulate the in-situ pore fluid flow as well as contaminant transport, which one-dimensional (1D) laboratory scale models are incapable of (Culligan-Hensley and Savvidou 1995). Arulanandan et al. (1988) developed scaling laws to be used with the centrifuge in simulating geoenvironmental problems. With the use of scaling laws, the phenomena being observed in the centrifuge can be correlated to the prototype, and centrifuge modeling can be applied in environmental engineering, including contaminants transport, multiphase flow, site remediation, and so on (Soga et al. 2003; Lo and Hu 2004a, b; Lo et al. 2005; Hu et al. 2006a, b; Meegoda and Hu 2011). Some of the scaling factors used during centrifugal tests are listed in Table 1. Recently the centrifuge technique was employed to investigate the mechanism of air sparging (Marulanda et al. 2000; Hu et al. 2010, 2011).
Table 1. Basic Scaling Laws Used in Centrifuge Modeling
ParameterScaling law
Gravity (LT2)N
Length (L)1/N
Mass (M)1/N3
Density (ML3)1
Pressure (MLT2)1
Intrinsic permeability (L2)1
Permeability (LT1)N
Viscosity (ML1T1)1
Time (T)1/N2
With centrifugal tests, in-situ stress state is always maintained while applying a wide range of air pressures during air sparging. Marulanda et al. (2000) investigated the air flow characteristics through pure fluid and saturated porous media by means of centrifugal tests under a wide range of gravitational accelerations, and the effect of injector geometry on the breakthrough velocity and air plume shape was studied under a specific injection pressure. Hu et al. (2010) investigated the effect of particle size and air injection pressure on size and shape of ZOI using a point sparger, and showed that ZOI can be expressed by two components: the horizontal expansion from preferential intrusion around the injection point and the angle of ZOI, which is the angle between the vertical direction and the boundary of ZOI. The centrifuge investigation by Hu et al. (2011) showed that air mass flow rate increases linearly with the effective sparging pressure ratio, which is the difference between sparging pressure and hydrostatic pressure normalized with respect to the effective overburden pressure at the sparging point. With increasing sparging pressure, the zone of influence will increase in both lateral intrusion and the cone angle until a limited ZOI, which will not expand further.
During all the above centrifugal studies (Marulanda et al. 2000; Hu et al. 2010, 2011), as the flow regimes under different g levels changed, however, the scaling factors for the ZOI was not defined. Hence, it is not clear if the test results under different g levels are comparable. In order to develop a scaling factor for ZOI, the applicability of centrifugal models to simulate the field air sparging should be carefully checked. The best way to very this would be to simulate actual field test in the centrifuge. This is rather difficult to accomplish unless air sparging test data are available from rather uniform field test sites. The modeling of models technique is an alternative technique to validate centrifuge test results. Arulanandan et al. (1988) showed that if several centrifugal model tests simulating the same prototype can predict similar prototype behavior, then centrifugal tests can be used to simulate that prototype geoenvironmental problem. In addition, during all these centrifugal studies mentioned above, point spargers, which simulate three-dimensional (3D) profiles, were employed; however, while performing result analysis, only two-dimensional (2D) images were captured and analyzed in using narrow or thin centrifugal models. The use of point spargers to observe 2D images introduces errors. Hence 2D spargers should be used to observe 2D images in order to obtain air migration mechanisms during air sparging.
There are two main research objectives of this research: one is to investigate the flow pattern transition during centrifugal tests under different g levels, and the other is to study the scaling factor for ZOI during centrifugal tests by considering flow pattern transition both using a 2D or line sparger. To obtain the scaling factor for ZOI and to validate the applicability of using the centrifuge to model air sparging, the modeling of models technique was employed in this research.

Experimental Procedure

Experimental Setup

A centrifuge strongbox was used as a soil container for the centrifugal experiments. The internal dimensions of the strongbox were 600mm(length)×800mm(height)×43mm(width). The strong box was fabricated using Plexiglas to view airflow patterns in soil during testing. A linear air sparger with a length of 43 mm, which was the same as the width of the centrifugal model, was installed at the middle of the strongbox to simulate 2D prototype tests. The 43-mm long line sparger had an aperture width of the sparger that was either 2.0 mm or 1.5 mm depending on the acceleration of the centrifugal test. The sparger was connected to a compressed air supply, which included a compressor, a regulator, a pressure gauge, a flow meter, and a valve. The regulator can alter the air pressure, and a constant air flow rate can be measured using the flow meter. The mass flow rate was computed using the known volume flow rate and the corresponding pressure. Three pore-water pressure transducers (PPTs) were installed at the bottom of the strongbox to monitor water pressure during air sparging. A schematic diagraph for centrifugal testing setup is shown in Fig. 1.
Fig. 1. Testing setup for centrifugal modeling of air sparging
Centrifugal tests were performed at the Tsinghua University geotechnical centrifugal facility. This centrifuge has a radius of 2.5 m and is capable of spinning a maximum payload of 250 kg (model mass). The maximum acceleration that can be applied is 200 times gravitational acceleration. 2D digital images were collected during centrifuge flight by a high resolution charge coupled device (CCD) camera. The data collecting system was mounted in the centrifuge and was transferred to the control room via slip rings.

Testing Materials

The opacity of natural soil prevents the direct observation of air flow patterns during air sparging. Hence in this study, transparent fused silica glass beads, having almost the similar physical properties as sand, were used as a replacement for sand. Water was used as pore fluid. Hu et al. (2010, 2011) showed that air movement in glass beads saturated with water could be adequately captured in flight by a video recorder.
Two types of glass beads, named coarse beads and fine beads, were used for the centrifugal models. The physical properties of the coarse glass beads are as follows: particle size 5.0–6.0 mm; effective size (d10) 5.1 mm; coefficient of uniformity (Cu) 1.10; and coefficient curvature (Cc) 0.98. The physical properties of the fine glass beads are as follows: particle size 0.8–1.0 mm; effective size (d10) 0.82 mm; coefficient of uniformity (Cu) 1.12; and coefficient curvature (Cc) 0.98.

Centrifugal Models

The dry pluviation method was adopted to obtain uniform samples (Hu et al. 2010). Then the sample was connected to a water tank allowing the water level to gradually rise and saturate the soil by soaking in water for 12 h before the centrifugal testing. During saturation, the water level was kept at 2 cm above the surface of the soil sample to simulate the constant groundwater level. The height of the glass beads model was 45 cm with a porosity of 0.34 for coarse glass beads and 0.41 for fine glass beads, and test samples were subjected to 5, 10, 15, 20, 30, 40, and 50 gravitational levels.
Additional two centrifugal tests were conducted (20 and 40 gravitational levels) with coarse glass beads to simulate the same prototype obtained from 15 to 30 gravitational levels tests with 45-cm high coarse glass beads models. For those two tests, the aperture width of linear air sparger and sample height were reduced to 1.5 mm and 33.7 cm respectively to produce the same prototype dimensions as the 15 and 30 gravitational level tests of the first test series.

Experiment Procedure

Centrifugal tests were performed under a wide range of sparging pressures and centrifugal gravitational levels for two soil types. The soil model was spun at a desired gravitational level for 30 min for further compaction and homogenization before air sparging tests. The water table was maintained at a constant level during all tests, and pore water pressures were measured using a miniature transducer near the sparger. Then the air pressure was gradually increased in small 5 kPa steps until air was able to be injected into the soil sample. This threshold air pressure was referred to as the minimum sparging pressure. Beyond this pressure level, the sparging pressure increased in 10 kPa steps. For each step, the pressure level was maintained for at least 10 min to stabilize the air flow and to facilitate full development of ZOI. The air flow rate was measured, and the full development of ZOI using a video camera was monitored. The centrifugal tests were terminated when there was no further expansion of the ZOI. The rise of the water level was observed in the centrifugal model during air sparging, which is from the air movement within ZOI. Typical image of the zone of influence in centrifugal model was shown in Fig. 2.
Fig. 2. Typical image of the zone of influence during air sparging in centrifugal model

Experimental Results

Modeling of Model Tests

The modeling of model tests was first performed to validate the applicability of a centrifugal testing technique for air sparging. Arulanandan et al. (1988) provided the scaling laws to be used in predicting prototype behavior. For example depth has a scaling factor of N. Hence, a 45-cm high model at a 15 gravitational level would produce a 6.75-m high prototype (15×45cm=6.75m). In this research, two prototypes were modeled using centrifugal tests at various gravitational levels.
1.
The first prototype was a 6.75-m high coarse glass beads layer with a 30-mm wide sparger. Two different centrifuge models were run to simulate the same prototype: one was a 45-cm high model with 2-mm wide sparger aperture at a 15 gravitational level, and the other was a 33.7-cm high model with 1.5-mm wide sparger at a 20 gravitational level.
2.
The second prototype was a 13.5-m high coarse glass beads layer with a 60-mm wide linear sparger, which was simulated using two centrifuge models: one is a 45-cm high model with 2-mm wide sparger at a gravitational level of 30 and the other a 33.7-cm high model with 1.5-mm wide sparger with a gravitational level of 40.
As per the test results shown in Table 2, both modeling of model tests produced prototypes with similar sized ZOI and bubble flow, confirming that air sparging tests could be reproduced in a centrifuge with a scaling law of N for a centrifugal model simulated at N gravitational level.
Table 2. Zone of Influence for Modeling of Model Tests
Model test numberCentrifugal modelPrototypeEffective sparging pressure (kPa)Cone angle from sparging point (°)Cone angle (°)Prototype lateral intrusion (cm)
Sparger width (mm)Sample height (cm)g levelSparger width (mm)Soil mass depth (m)
Test 12.04515306.751607.257.022.5
Test 21.533.720306.751607.257.028.0
Test 32.045306013.52005.755.530.0
Test 41.533.7406013.52006.055.836.0
The effective sparging pressure, which is defined as the air injection pressure subtracting the corresponding hydrostatic pressure, is used to study the correlation between air flow rate and air injection pressure under different gravitational levels. Table 3 lists the minimum effective entry air pressures. Although the entry pressure increased with the increase in centrifugal gravitational level, the minimum effective entry air pressure did not change with the centrifugal gravitational level. Hence it can be stated that minimum effective entry air pressure only depends on soil properties, such as the grain size distribution and the pore network structure. The relationship between air mass flow rate per unit width and effective sparging pressure for four modeling of model tests is shown in Fig. 3. From Fig. 3 one can observe that for the same prototype, the variation of air mass flow rate per unit width with effective sparging pressure for different model tests are similar. In the following sections, the air flow rate, characteristics of ZOI, and air flow patterns and air flow pattern transition are further investigated and discussed.
Table 3. Minimum Injection Pressure and Zone of Influence for Air Sparging with Various g Levels
Particle sizeg levelSample height (cm)Minimum effective AS pressure Pmin (kPa)Cone angle (°)Length of lateral intrusion (cm)Cone angle from sparging point (°)
Coarse1450.112.02.014.4
Coarse5456.38.51.910.9
Coarse10457.37.51.59.4
Coarse15456.07.01.58.9
Coarse20456.46.51.28.0
Coarse30458.55.51.06.8
Coarse40457.65.00.55.6
Coarse504510.55.00.35.4
Fine1450.917.04.422.0
Fine5457.516.52.519.4
Fine10457.610.02.112.6
Fine15456.89.51.511.3
Fine20457.08.01.49.7
Fine30458.26.01.37.6
Fine40459.05.50.76.4
Fine504510.84.00.74.9
Fig. 3. Relationship between air mass flow rate per unit width and effective sparging pressure for modeling of model tests

Air Flow Rate during Centrifugal Tests

The stable air sparging was obtained when a constant air flow rate occurred under a given sparging pressure. The effective sparging pressure and corresponding air mass flow rate Q were monitored during stable air sparging, and the relationship was shown in Fig. 4. The air mass flow rate Q increased with the increase in injection pressure for all tests. Initially with the increase of air pressure, increase of flow rate was due to the pressure gradient and the apparent expansion of ZOI. Then the ZOI reached its maximum size, and the pressure gradient became the only contributing factor for air flow rate. Therefore, the slope of air mass flow rate to effective sparging pressure was larger at the early stage during the increase of injection pressure and leveled off after full development of the ZOI.
Fig. 4. Relationship between air mass flow rate and effective sparging pressure for coarse glass beads

Development of the ZOI with Increase of Sparging Pressure

The captured camera images show that the observed shape of 2D ZOI was a truncated cone. However, the camera images obtained during this 2D test are much sharper than those reported by Hu et al. (2010, 2011) using a point sparger. Fig. 5 shows a selected stable ZOI boundary for different gravitational levels under similar normalized effective sparging pressure. For the data shown in Fig. 5, the normalized effective sparging pressure was chosen as 0.5 for all gravitational level tests. The test results showed that for the same normalized effective sparging pressure, the shapes of ZOI are quite similar for different gravitational levels. Fig. 6 shows that the cone angle increased with the increase of air sparging pressure and reached a stable value of around 7° for the model with coarse glass beads. The ZOI expanded with the increase in injection air pressure when the sparging pressure was low. Hu et al. (2010) reported that ZOI contains two components: the lateral expansion due to pneumatic fracturing or preferential intrusion around the injection port and the cone angle between the vertical direction and boundary of ZOI, and this conclusion was confirmed by results from this research. The increase in air pressure produced well-developed air pathways. However, with the increase of the sparging pressure, the cone angle rapidly reached a limiting value, and thereafter the lateral expansion was the main contributor to the expansion of the ZOI.
Fig. 5. Selected stable ZOI boundary under normalized effective sparging pressure of 0.5 for different gravitational levels tests
Fig. 6. Variation of ZOI angle with increase of effective sparging pressure ratio for different gravitational levels for coarse beads model
Table 3 summarizes the ultimate boundaries of ZOI under various gravitational levels of centrifuge tests for the fine glass beads for both coarse beads and fine beads. According to Table 3, the ZOI decreased as the centrifugal gravitational levels increased. It can be further observed that for 1 and 5 gravitational levels tests for fine glass beads, the cone angle values were similar but much higher than that for 10g tests because of different air flow patterns. This observation will be discussed later in this manuscript. Please note that previous researchers have defined the ZOI as a cone, and the slope of the sparging angle as the slope of the line connecting the sparging point to the edge of the ZOI based on 1g tests (Reddy and Adams 2008). Hence the reported angles of the ZOI values in Table 3 are much smaller than previously reported values for similar soil types.

Air Flow Pattern Transition

Bubble flow patterns were observed in all centrifugal tests with coarse glass beads. In centrifugal tests with fine glass beads, both channel and bubble flows were observed. Air flowed through microchannels in 5 gravitational level centrifugal tests with fine glass beads, whereas those for 10 and higher gravitational levels tests bubble flow was observed. In the following part, the air flow pattern transition is discussed.
When the air sparging is initiated, the air bubbles are formed at the sparging point, and they enter the soil matrix. The minimum sparging or injection pressure Pmin is calculated by the balancing force from injection pressure by that due to surface tension as given below (Marulanda et al. 2000)
Pmin=Phyd+Pcap=ρghw+4σcosθD
(1)
Where Phyd is hydrostatic pressure at the sparging point, Pcap is the capillary pressure for the soil, σ is the interfacial tension between air and the pore water that is 0.074N/m at 20°C, ρ is the density of the pore water, hw is the height of the water table above the sparge point, θ is the contact angle that was assumed equal to zero for complete wetting, and D is the average pore size diameter assumed as d10/5. Once an air bubble is formed, it is moved into the path of least resistance from the pressure gradient created by the sparging pressure. Hence, depending on the inclination of injection port and the hydraulic conductivity tensor of the pore structure, air bubbles can move both vertically as well as laterally due to anisotropy (Meegoda et al. 1989).
Selker et al. (2007) proposed that during air sparging tests, the zone of influence can be separated into two regions: one is around the air sparger called near source region whereas the other one is far from the air sparger called far source region. The main characteristic of air driven mode in the near source region distinguishes from that in far source region is the driving force of air flow. Within near source region, air is primarily driven by gradients of air pressure whereas in far source region, the buoyant force is the main driving force for the air phase. In the far source region, the bubble movement is governed by the following three forces, including the surface tension force (πDσcosθ), buoyant force [πD3g(ρρa)/6], and frictional or Stokes force (3πDμv), where ρa is the air density, which is approximately zero, μ is absolute viscosity which in 0.001Pa·s, and v is the bubble velocity that varied from 30 to 200cm/s for all tests. Please note that the two pore diameters (D10/5) for fine and coarse glass beads are 0.164 mm and 1.02 mm, respectively. When the numerical values of three forces are calcuated for fine glass beads, it was found that buoyant force was at least two orders of magnitude smaller than other two forces for 1 gravitational level test. However, with higer gravitational levels, magnitude of buoyant force becomes comparable with the other two forces. Under such situations, bubble flow occurs. Under low gravitational levels, bubbles collapse into the channel flow, so that a buoyant force of combined tube shape bubbles in the channel becomes comparable with the other two forces that cause the tube shape bubbles to rise. The low cone angles observed in this research for high gravitational levels suggest that there was bubble flow in both fine and coarse glass bead samples. This observation also suggests that in order to apply the scaling law of N, particle size should also be modeled. However, this requirement could be avoided if flow regimes of both prototype and the model are similar. This can be achieved by equating the buoyant force to surface tension force and computing a test diameter [Dtest=(6σ/ρNg)0.5]. All the test diameter values in all tests conducted in this research are listed in Table 4.
Table 4. Test Diameters of All Tests
Particle sizeg levelDtest=(6σ/ρNg)0.5(mm)Pore diameter (d10/5) (mm)Ratio of test diameter to pore diameter
Coarse16.731.027
Coarse53.011.023
Coarse102.131.022
Coarse151.741.022
Coarse201.511.021
Coarse301.231.021
Coarse401.061.021
Coarse500.951.021
Fine16.730.1642
Fine53.010.1619
Fine102.130.1613
Fine151.740.1611
Fine201.510.169
Fine301.230.168
Fine401.060.167
Fine500.950.166
The relationship between the ratio of test diameter and the pore diameter and cone angle is drawn in Fig. 7. From Fig. 7 and Table 4, it can be observed that for the fine glass beads, the ratio of test diameter to pore diameter sharply changed from 19 (for 5g centrifugal test) to 13 (10g centrifugal test), which indicates an air pattern transition from bubble flow under higher gravitation levels to channel flow of the ZOI under lower gravitational level. Hence, one could estimate the critical value of the ratio of test diameter to pore diameter (Dtest/D) where a transition from bubble to channel flow occurs. Based on Fig. 7, it can be concluded that if the Dtest/D value is greater than about 19, channel flow would occur. Also a Dtest/D value less than 4 would gurarantee stable bubble flow pattern irrespective of the soil particle size. Hence, with the above information, centrifuge modeling can be used as a design tool for remediating contaminated soils used in-situ sparging eliminating the need for expensive pilot tests. However, the soil heterogeneity plays a significant role in ZOI dimensions, so the pilot tests are still recommended for air sparging system design sites with heterogeneous soil deposits. Based on the above discussion, the observation from Fig. 7 and Table 3 is consisent with the conclusion that the cone angle of the ZOI is much higher for channel flow than that for bubble flow.
Fig. 7. Relationship between ratio of test diameter to pore diameter and cone angle
In this research, the centrifugal tests under different gravitational levels were performed to investigate air flow pattern transition during centrifugal tests. It was found that in the centrifugal tests, air migrates in the form of discrete bubble instead of air channel at high gravitational levels. A test diameter [Dtest=(6σ/ρNg)0.5] was computed by equating the buoyant force to surface tension force. The ratio of test diameter to pore diameter can help inferring the flow regimes.
Similarly, Geistlinger et al. (2006) proposed a flow rate and grain-size dependent stability criterion based on constricted capillary tube model. They used carefully conducted experimental data from Brooks et al. (1999) with a wide range of grain sizes and air flow rates. The air flow rate depends on the size of the sparger. Hence, the variation of air flow rate per unit area values with grain-size from this research as well as from Brooks et al. (1999) were plotted in Fig. 8. For the centrifugal test data, the grain size is multiplied by the gravitational level to theoretically account for the increased gravity in the centrifuge. Please note that only confirmed bubble and channel flow data from the centrifugal tests were included in Fig. 8. Fig. 8 shows there is a clear boundary from bubble flow to channel flow and that boundary is a straight line, which can be described as
q=0.2·D4
(2)
where q = air flow rate per unit area in ml/min/mm2, and D = grain size in mm.
Fig. 8. Classification of channel flow and bubble flow
Based on the reported results shown in Fig. 8, if the Poiseuille equation is used to back calculate the air flow rate assuming an average air channel diameter of around D/4, one would obtain Eq. (2). Hence, based on Fig. 8 and Eq. (2), it can be concluded that the initiation of channel flow in a soil mass subjected to air sparging may be described using Poiseuille equation. The above finding is very preliminary and should be confirmed by carefully conducted micro mechanics research program (Gao et al. 2011, 2012, 2013a, b). Also, previous researchers were not able to obtain channel flow patterns for grain diameters larger than 4 mm. However, based on Fig. 8 with centrifugal tests, there should not be such limitation. Channel flow requires higher pressure-gradients than that for bubble flows for any grain size. Based on results from Brooks et al. (1999), the maximum gradient that could be applied to soil is around 3.0 without causing soil liquefaction. Based on Eq. (2), this would occur around grain diameter of 4 mm, making it impossible to have channeled flows for grain diameters larger than 4 mm for small-scale laboratory tests conducted under regular gravitational levels. However, with the centrifugal tests, there is no such limitation, and Eq. (2) is valid for all grain sizes.

Summary and Conclusions

The modeling of model tests was first performed to validate the applicability of centrifugal tests to model air migration through saturated soils during air sparging, and confirmed that air sparging tests could be simulated using the centrifuge. Then a series of 2D centrifugal model tests were performed under various gravitational levels to study air flow patterns transition under prototype stresses and with a wide range of air injection pressures. The following conclusions can be drawn from this research:
The observed 2D ZOI was a truncated cone under various gravitational levels, and ZOI under low gravitational level is always larger than that under high gravitational level. In addition, the observed cone angle of the ZOI is much higher for channel flow, when compared with that for bubble flow.
For the same normalized effective air sparging pressure at the injection port, sizes of the 2D ZOI under different gravitational levels are quite similar.
The research confirmed that the scaling factor for air mass flow during air sparging is N, if there are similar air flow regimes during air sparging tests under varied gravitational levels.
It was observed that flow pattern for the same soil changed from channel flow under lower gravitational levels to bubble flow under higher gravitational levels. In order to make the centrifugal tests results under varied gravitational levels comparable, soil should be selected such that its size is scaled by N. This requirement of scaling of particles could be avoided by ensuring similar air flow regimes during air sparging modeling tests under varied gravitational levels. A test diameter [Dtest=(6σ/ρNg)0.5] was computed by equating the buoyant force to surface tension force. The ratio of test diameter to pore diameter can help inferring the flow regimes. Based on the test results, if the ratio of test diameter to pore diameter was greater than about 15, one would expect channel flow. Stable bubble flow patterns would occur irrespective of the soil particle size if the ratio of test diameter to pore diameter (Dtest/D) is less than 4.
The flow rate per unit area versus grain size data were plotted using centrifugal test results from this research as well as those from 1g test with either clear bubble flow or channel flows. There is a clear demarcation between bubble flow and channel flow. The air flow pattern transition can be predicted using Poiseuille equation.

Acknowledgments

The financial support from State Key Laboratory of Hydro Science and Engineering (SKLHSE-2012-KY-01, SKLHSE-2013-D-01), National Key Basic Research Program (2012CB719804), and National Natural Science Foundation of China (Project No. 51323014, 41372352, 50879038) are gratefully acknowledged.

References

Adams, J. A., and Reddy, K. R. (2003). “Extent of benzene biodegradation in saturated soil columns during air sparging.” Ground Water Monit. Rem., 23(3), 85–94.
Arulanandan, K., Thompson, P. Y., Kutter, B. L., Meegoda, N. J., Muraleetharan, K. K., and Yogachandran, C. (1988). “Centrifuge modeling of transport processes for pollutants in soils.” J. Geotech. Geoenviron. Eng., 185–205.
Brooks, M. C., Wise, W. R., and Annable, M. D. (1999). “Fundamental changes in in situ air sparging flow patterns.” Ground Water Monit. Rem., 19(2), 105–113.
Culligan-Hensley, P. J., and Savvidou, C. (1995). “Environmental geomechanics and transport processes.” Geotechnical centrifuge technology, R. N. Taylor, ed., Blackie Academic and Professional, Glasgow, U.K.
Elder, C., and Benson, C. (1999). “Air channel formation, size, spacing and tortuosity during air sparging.” Ground Water Monit. Rem., 19(3), 171–181.
Gao, S., Meegoda, J. N., and Hu, L. (2011). “Microscopic modeling of air migration during air sparging.” J. Hazard., Toxic Radioactive Waste Manage., 70–79.
Gao, S., Meegoda, J. N., and Hu, L. (2013a). “A dynamic two-phase flow model for air sparging.” Int. J. Numer. Anal. Methods Geomech., 37(12), 1801–1821.
Gao, S. Y., Meegoda, J. N., and Hu, L. (2012). “Two methods for pore-network of porous media.” Int. J. Numer. Anal. Methods Geomech., 36(18), 1954–1970.
Gao, S. Y., Meegoda, J. N., and Hu, L. (2013b). “Simulation of dynamic two-phase flow during multistep air sparging.” Transp. Porous Media, 96(1), 173–192.
Geistlinger, H., Krauss, G., Lazik, D., and Luckner, D. L. (2006). “Direct gas injection into saturated glass beads: Transition from incoherent to coherent gas flow pattern.” Water Resour. Res., 42(7), W07403.
Hu, L., and Liu, Y. (2008). “Physical modeling of air-sparging technique for groundwater remediation.” Chin. J. Geotech. Eng., 30(6), 835–839.
Hu, L., Lo, I. M. C., and Meegoda, J. N. (2006a). “Centrifuge testing of LNAPL migration and soil vapor extraction for soil remediation.” Pract. Period. Hazard. Toxic Radioact. Waste Manage., 33–40.
Hu, L., Lo, I. M. C., and Meegoda, J. N. (2006b). “Numerical analysis and centrifugal modeling of LNAPLs transport in subsurface system.” Prog. Nat. Sci., 16(4), 416–424.
Hu, L., Meegoda, J., Du, J., Gao, S., and Wu, X. (2011). “Centrifugal study of zone of influence during air-sparging.” J. Environ. Monit. RSC, 13(9), 2443–2449.
Hu, L., Wu, X., Liu, Y., Meegoda, J. N., and Gao, S. (2010). “Physical modeling of air flow during air sparging remediation.” Environ. Sci. Technol., 44(10), 3883–3888.
Jang, W., and Aral, M. M. (2009). “Multiphase flow fields in in-situ air sparging and its effect on remediation.” Transp. Porous Media, 76(1), 99–119.
Ji, W., Ahefeld, A., Lin, D. P., and Hill, E. (1993). “Laboratory study of air sparging: Air flow visualization.” Ground Water Monit. Rem., 13(4), 115–126.
Johnson, R., Johnson, P., McWorter, D., Hinchee, R., and Goodman, I. (1993). “An overview of in situ air sparging.” Ground Water Monit. Rem., 13(4), 127–135.
Lo, I. M. C., and Hu, L. (2004a). “Centrifuge modeling of LNAPLs transport in unsaturated soils.” J. Geotech. Geoenviron. Eng., 535–539.
Lo, I. M. C., and Hu, L. (2004b) “Long-term migration of LNAPL in unsaturated soils: Clayey silt and fine sand.” Pract. Period. Hazard. Toxic Radioact. Waste Manage., 228–237.
Lo, I. M. C., Hu, L., and Meegoda, N. J. (2005). “Feasibility study of using centrifuge for investigation of LNAPL migration in unsaturated soils.” Soil Sediment Contam., 14(1), 85–103.
Marulanda, C., Culligan, P. J., and Germaine, J. T. (2000). “Centrifuge modeling of air sparging—A study of air flow through saturated porous media.” J. Hazard. Mater., 72(2–3), 179–215.
McCray, J. (2000). “Mathematical modeling of air sparging for subsurface remediation: State of the art.” J. Hazard. Mater., 72(2–3), 237–263.
McCray, J. E., and Falta, R. W. (1997). “Numerical simulation of air sparging for remediation of NAPL contamination.” Ground Water, 35(1), 99–110.
Meegoda, J. N., and Hu, L. (2011). “A review of centrifugal testing of gasoline contamination and remediation.” Int. J. Environ. Res. Publ. Health, 8(12), 3496–3513.
Meegoda, N. J., King, I. P., and Arulanandan, K. (1989). “An expression for permeability of anisotropic granular media.” Int. J. Numer. Anal. Meth. Geomech., 13(6), 575–598.
Mei, C. C., and Cheng, Z. (2002). “A model for flow induced by steady air venting and air sparging.” Appl. Math. Model., 26(7), 727–750.
Nyer, E. K., and Suthersan, S. S. (1993). “Air sparging: Savior of ground water remediation or just blowing bubbles in the bathtub.” Ground Water Monit. Rem., 13(4), 87–91.
Reddy, K. R., and Adams, J. (2008). “Conceptual modeling of air sparing for groundwater remediation.” Proc., 9th Int. Symp. on Environmental Geotechnology and Global Sustainable Development, International Society for Environmental Geotechnology (ISEG).
Reddy, K. R., and Adams, J. A. (2001). “Effect of soil heterogeneity on airflow patterns and hydrocarbon removal during in situ air sparging.” J. Geotech. Geoenviron. Eng., 234–247.
Selker, J. S., Niemet, M., Mcduffie, N. G., Gorelick, G. M., and Parlange, J. Y. (2007). “The local geometry of gas injection into saturated homogeneous porous media.” Transp. Porous Media, 68(1), 107–127.
Soga, K., Kawabata, J., Kechavarzi, C., Coumoulos, H., and Waduge, W. A. P. (2003). “Centrifuge modeling of nonaqueous phase liquid movement and entrapment in unsaturated layered soils.” J. Geotech. Geoenviron. Eng., 173–182.
U.S. EPA. (1995). “How to evaluate alternative cleanup technologies for underground storage tank sites.”, Washington, DC.
Wang, J., and Hu, L. (2012). “Numerical simulation of air flow during air sparging remediation.” Appl. Mech. Mater., 138–189, 27–32.

Information & Authors

Information

Published In

Go to Journal of Environmental Engineering
Journal of Environmental Engineering
Volume 141Issue 1January 2015

History

Received: Jul 27, 2013
Accepted: Jun 18, 2014
Published online: Aug 1, 2014
Published in print: Jan 1, 2015
Discussion open until: Jan 1, 2015

Authors

Affiliations

Liming Hu, A.M.ASCE [email protected]
Associate Professor, State Key Laboratory of Hydro-Science and Engineering, Dept. of Hydraulic Engineering, Tsinghua Univ., Haidian District, Beijing 100084, P.R. China (corresponding author). E-mail: [email protected]
Jay N. Meegoda, F.ASCE [email protected]
Visiting Professor, Dept. of Hydraulic Engineering, Tsinghua Univ., Haidian District, Beijing 100084, P.R. China; and Professor, Dept. of Civil and Environmental Engineering, New Jersey Institute of Technology, Univ. Heights, Newark, NJ 07102. E-mail: [email protected]
Hengzhen Li
Research Assistant, State Key Laboratory of Hydro-Science and Engineering, Dept. of Hydraulic Engineering, Tsinghua Univ., Haidian District, Beijing 100084, P.R. China.
Jianting Du
Research Assistant, State Key Laboratory of Hydro-Science and Engineering, Dept. of Hydraulic Engineering, Tsinghua Univ., Haidian District, Beijing 100084, P.R. China.
Shengyan Gao
Research Associate, Dept. of Hydraulic Engineering, Tsinghua Univ., Haidian District, Beijing 100084, P.R. China; and Postdoctoral Fellow, Dept. of Civil and Environmental Engineering, New Jersey Institute of Technology, Univ. Heights, Newark, NJ 07102.

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share