Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 1. Experimental studies

A matrix of batch, column and two-dimensional (2-D) box experiments was conducted to investigate the coupled effects of rate-limited solubilization and layering on the entrapment and subsequent recovery of a representative dense NAPL, tetrachloroethylene (PCE), during surfactant flushing. Batch experiments were performed to determine the equilibrium solubilization capacity of the surfactant, polyoxyethylene (20) sorbitan monooleate (Tween 80), and to measure fluid viscosity, density and interfacial tension. Results of one-dimensional column studies indicated that micellar solubilization of residual PCE was rate-limited at Darcy velocities ranging from 0.8 to 8.2 cm/h and during periods of flow interruption. Effluent concentration data were used to develop effective mass transfer coefficient (Ke) expressions that were dependent upon the Darcy velocity and duration of flow interruption. To simulate subsurface heterogeneity, 2-D boxes were packed with layers of F-70 Ottawa sand and Wurtsmith aquifer material within 20-30 mesh Ottawa sand. A 4% Tween 80 solution was then flushed through PCE-contaminated boxes at several flow velocities, with periods of flow interruption. Effluent concentration data and visual observations indicated that both rate-limited solubilization and pooling of PCE above the fine layers reduced PCE recovery to levels below those anticipated from batch and column measurements. These experimental results demonstrate the potential impact of both mass transfer limitations and subsurface layering on the recovery of PCE during surfactant enhanced aquifer remediation.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 2. Numerical simulation

A numerical model of surfactant enhanced solubilization was developed and applied to the simulation of nonaqueous phase liquid recovery in two-dimensional heterogeneous laboratory sand tank systems. Model parameters were derived from independent, small-scale, batch and column experiments. These parameters included viscosity, density, solubilization capacity, surfactant sorption, interfacial tension, permeability, capillary retention functions, and interphase mass transfer correlations. Model predictive capability was assessed for the evaluation of the micellar solubilization of tetrachloroethylene (PCE) in the two-dimensional systems. Predicted effluent concentrations and mass recovery agreed reasonably well with measured values. Accurate prediction of enhanced solubilization behavior in the sand tanks was found to require the incorporation of pore-scale, system-dependent, interphase mass transfer limitations, including an explicit representation of specific interfacial contact area. Predicted effluent concentrations and mass recovery were also found to depend strongly upon the initial NAPL entrapment configuration. Numerical results collectively indicate that enhanced solubilization processes in heterogeneous, laboratory sand tank systems can be successfully simulated using independently measured soil parameters and column-measured mass transfer coefficients, provided that permeability and NAPL distributions are accurately known. This implies that the accuracy of model predictions at the field scale will be constrained by our ability to quantify soil heterogeneity and NAPL distribution.     

Effect of nonionic surfactant partitioning on the dissolution kinetics of residual perchloroethylene in a model porous medium

At concentrations above the critical micelle concentration, surfactants can significantly enhance the solubilization of residual nonaqueous phase liquids (NAPL) and, for this reason, are the focus of research on surfactant-enhanced aquifer remediation (SEAR). As a consequence of their amphiphilic nature, surfactants may also partition to various extents between the organic and aqueous phases, thereby affecting SEAR performance. We report here on the observation and analysis of the effect of surfactant partitioning on the dissolution kinetics of residual perchloroethylene (PCE) by aqueous solutions (1000 mg/L) of the non-ionic surfactant Triton X-100 in a model porous medium. For this fluid system, batch equilibration experiments showed that the surfactant partitions strongly into the NAPL (NAPL-water partition coefficient equal to 12.5). Dynamic interfacial tension (IFT) measurements were employed to study surfactant diffusion and interfacial adsorption. The dynamic IFT measurements were consistent with partitioning of the surfactant between the two liquid phases. PCE dissolution experiments, conducted in a transparent glass micromodel using an aqueous surfactant solution, were contrasted to experiments using clean water. Surfactant partitioning was observed to delay significantly the onset of micellar solubilization of PCE, an observation reproduced by a numerical model. This effect is attributed to the reduction of surfactant concentration in the immediate vicinity of the NAPL-water interface, which accompanies transport of the surfactant into the NAPL. Accordingly, it is suggested that both the rate and the extent of diffusion of the surfactant into the NAPL affect the onset of and the driving force for micellar solubilization. While many surfactants do not readily partition in NAPL, this possibility must be considered when selecting non-ionic surfactants for the enhanced solubilization of residual chlorinated solvents in porous media.     

Effect of scale and dimensionality on the surfactant-enhanced solubilization of a residual DNAPL contamination

The mass transfer rate from residual dense non-aqueous phase liquids (DNAPLs) to the mobile aqueous phase is an important parameter for the efficiency of surfactant-enhanced remediation through solubilization of this type of contamination. The mass transfer kinetics are highly dependent on the dimensionality of the system. In this study, irregularly shaped residual TCE saturations in two-dimensional saturated flow fields were flushed with a 2% polyoxyethylene sorbitan (20) monooleate (POESMO) solution until complete removal had been achieved. A numerical model was developed and used for the simulation of the various surfactant-flushing experiments with different initial saturation patterns and flow rates. Through optimization against in situ concentration and saturation data, a phenomenological power-law model for the relationship between the mass transfer rate from the DNAPL to the mobile aqueous phase on the one hand and the residual DNAPL saturation and the flow velocity on the other hand was derived. The obtained mass transfer rate parameters provide a reasonable fit to the experimental data, predicting the cleanup time and the general saturation and concentration pattern quite well but failing to predict the concentration curves at every individual sampling port. The obtained mass transfer rate model gives smaller values for the predicted mass transfer rate but shows a comparable dependence on water flow and saturation as in earlier published one-dimensional column experiments with identical characteristics for porous medium, DNAPL and surfactant. Mass transfer rate predictions were about one order of magnitude lower in the 2-D flow cell experiment than in 1-D column experiments. These results give an indication for the importance of dimensionality during surfactant remediation.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

A design of a vermiculite column adsorber for the removal of lead from water

Mini-column techniques were employed to determine the mass transfer coefficient for lead adsorption onto vermiculite. Variation of the mass transfer coefficient with flow rate, particle size of sorbent, and influent lead concentration were studied. Multiple linear regression (MLR) analysis showed that the mass transfer coefficient varied as the 0.43 power of the liquid flow rate and inversely as the 0.272 power of particle diameter of the vermiculite, but was independent of the influent concentration of lead. Different parameters of a fixed bed column design for the removal of lead by vermiculite were determined using the data from the batch sorption study. The performance of the liquid bed column in removing lead was in close agreement with predicted performance using the batch isotherm data.     

Biosurfactant-enhanced solubilization of NAPL mixtures

Remediation of nonaqueous phase liquids (NAPLs) by conventional pump-and-treat methods (i.e., water flushing) is generally considered to be ineffective due to low water solubilities of NAPLs and to mass-transfer constraints. Chemical flushing techniques, such as surfactant flushing, can greatly improve NAPL remediation primarily by increasing the apparent solubility of NAPL contaminants. NAPLs at hazardous waste sites are often complex mixtures. However, the equilibrium and nonequilibrium mass-transfer characteristics between NAPL mixtures and aqueous surfactant solutions are not well understood. This research investigates the equilibrium solubilization behavior of two- and three-component NAPL mixtures (containing akylbenzenes) in biosurfactant solutions. NAPL solubilization is found to be ideal in water (i.e., obeys Raoult's Law), while solubilization in biosurfactant solutions was observed to be nonideal. Specifically, the relatively hydrophobic compounds in the mixture experienced solubility enhancements that were greater than those predicted by ideal enhanced solubilization theory, while the solubility enhancements for the relatively hydrophilic compounds were less than predicted. The degree of nonideality is shown to be a nonlinear function of the NAPL-phase mole fraction. Empirical relationships based on the NAPL-phase mole fraction and/or micelle-aqueous partition coefficients measured in single-component NAPL systems are developed to estimate values for the multicomponent partition coefficients. Empirical relationships that incorporate both the NAPL-phase mole fraction and single-component partition coefficients yield much improved estimates for the multicomponent partition coefficient.     

Effects of ethanol addition on micellar solubilization and plume migration during surfactant enhanced recovery of tetrachloroethene

Alcohol addition has been suggested for use in combination with surfactant flushing to enhance solubilization kinetics and permit density control of dense non-aqueous phase liquid (DNAPL)-laden surfactant plumes. This study examined the effects of adding ethanol (EtOH) to a 4% Tween 80 (polyoxyethylene (20) sorbitan monooleate) solution used to flush tetrachloroethene (PCE)-contaminated porous media. The influence of EtOH concentration, subsurface layering and scale on flushing solution delivery and PCE recovery was investigated through a combination of experimental and mathematical modeling studies. Results of batch experiments demonstrated that the addition of 2.5%, 5% and 10% (wt.) EtOH incrementally increased the PCE solubilization capacity and viscosity of the surfactant solution, while reducing solution density from 1.002 to 0.986 g/cm3. Effluent concentration data obtained from one-dimensional (1-D) column experiments were used to characterize rate-limited micellar solubilization of residual PCE, which was strongly dependent upon flow velocity and weakly dependent upon EtOH concentration. Two-dimensional (2-D) box studies illustrated that minor differences (0.008 g/cm3) between flushing and resident solution density can strongly influence surfactant front propagation. A two-dimensional multiphase simulator, MISER, was used to model the influence of EtOH composition on the aqueous flow field and PCE mass recovery. The ability of the numerical simulator to predict effluent concentrations and front propagation was demonstrated for both 1-D columns and 2-D boxes flushed with EtOH-amended Tween 80 solutions. Results of this study quantify the potential influence of alcohol addition on surfactant solution properties and solubilization capacity, and demonstrate the importance of considering small density variations in remedial design.     

Investigation of surfactant-enhanced dissolution of entrapped nonaqueous phase liquid chemicals in a two-dimensional groundwater flow field.

Because of their low solubility, waste chemicals in the form of nonaqueous phase liquids (NAPLs) that are entrapped in subsurface formations act as long-term sources of groundwater contamination. In the design of remediation schemes that use surfactants, it is necessary to estimate the mass transfer rate coefficients under multi-dimensional flow fields that exit at field sites. In this study, we investigate mass transfer under a two-dimensional flow field to obtain an understanding of the basic mechanisms of surfactant-enhanced dissolution and to quantify the mass transfer rates. Enhanced dissolution experiments in a two-dimensional test cell were conducted to measure rates of mass depletion from entrapped NAPLs to a flowing aqueous phase containing a surfactant. In situ measurement of transient saturation changes using a gamma attenuation system revealed dissolution patterns that are affected by the dimensionality of the groundwater flow field. Numerical modeling of local flow fields that changed with time, due to depletion of NAPL sources, enabled the examination of the basic mechanisms of NAPL dissolution in complex groundwater systems. Through nonlinear regression analysis, mass transfer rates were correlated to porous media properties, NAPL saturation and aqueous phase velocity. Results from the experiments and numerical analyses were used to identify deficiencies in existing methods of analysis that uses assumptions of one-dimensional flow, homogeneity of aquifer properties, local equilibrium and idealized transient mass transfer.     

Anionic surfactant remediation of soil columns contaminated by nonaqueous phase liquids

A variety of column experiments have been completed for the purpose of selecting and evaluating suitable surfactants for remediation of nonaqueous phase liquids (NAPLs). The various NAPLs tested in the laboratory experiments were tetrachloroethylene (PCE), trichloroethylene (TCE), jet fuel (JP4) and a dense nonaqueous phase liquid from a site at Hill Air Force Base, UT. Both Ottawa sand and Hill field soil were used in these experiments. Surfactant candidates were first screened using phase behavior experiments and only the best ones were selected for the subsequent column experiments. Surfactants which showed high contaminant solubilization, fast coalescence times, and the absence of liquid crystal phases and gels during the phase behavior experiments were tested in soil column experiments. The primary objective of the soil column experiments was to identify surfactants that recovered at least 99% of the contaminant. The secondary objective was to identify surfactants that show low adsorption and little or no loss of hydraulic conductivity during the column experiments. Results demonstrated that up to 99.9% of the contaminants were removed as a result of surfactant flooding of the soil columns. The addition of xanthan gum polymer to the surfactant solution was shown to increase remediation efficiency as a lower volume of surfactant was required for recovering a given volume of NAPL. Based on these experimental results, guidelines for designing highly efficient and robust surfactant floods have been developed and applied to a field demonstration.     

Development of a continuously stirred flow cell for investigating sorption mass transfer

This paper introduces a new reversible-flow design for a continuously stirred reactor used to study sorption mass transfer in soil and solvent systems. The stirred reactor has potential advantages over conventional packed column or batch reactors because it isolates intraparticle sorption rate limitations from advective-dispersive transport, yet allows changes to flux through the reactor for analysis of sorption kinetics under dynamic conditions. Previously, stirred reactors have often failed due to clogging of sediment on the effluent frit. The reverse-flow backwashing design allows longer life and higher confidence in maintaining mixed conditions than previous designs. Mass transfer 'rate coefficients estimated from stirred and column experiments are compared; both techniques produced results consistent with a published correlation. The data also show that fitted sorption mass transfer coefficients can be strongly dependent on the choice of equilibrium partition coefficient (i.e. batch or first-moment derived values), and that the conventional two-site sorption kinetics model fails to accurately predict sorption mass transfer in the presence of changing solvent velocity through the reactor.     

Mass flux from a non-aqueous phase liquid pool considering spontaneous expansion of a discontinuous gas phase

The partitioning of non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase results in the repeated spontaneous expansion, snap-off, and vertical mobilization of the gas phase. This mechanism has the potential to significantly affect the mass transfer processes that control the dissolution of NAPL pools by increasing the vertical transport of NAPL mass and increasing the total mass transfer rate from the surface of the pool. The extent to which this mechanism affects mass transfer from a NAPL pool depends on the rate of expansion and the mass of NAPL compound in the gas phase. This study used well-controlled bench-scale experiments under no-flow conditions to quantify for the first time the expansion of a discontinuous gas phase in the presence of NAPL. Air bubbles placed in glass vials containing NAPL increased significantly in volume, from a radius of 1.0 mm to 2.0 mm over 215 days in the presence of tetrachloroethene (PCE), and from a radius of 1.2 mm to 2.3 mm over 22 days in the presence of trans-1,2-dichloroethene (tDCE). A one-dimensional mass transfer model, fit to the experimental data, showed that this expansion could result in a mass flux from the NAPL pool that was similar in magnitude to the mass flux expected for the dissolution of a NAPL pool in a two-fluid (NAPL and water) system. Conditions favouring the significant effect of a discontinuous gas phase on mass transfer were identified as groundwater velocities less than ~0.01 m/day, and a gas phase that covers greater than ~10% of the pool surface area and is located within ~0.01 m of the pool surface. Under these conditions the mass transfer via a discontinuous gas phase is expected to affect, for example, efforts to locate NAPL source zones using aqueous concentration data, and predict the lifetime and risk associated with NAPL source zones in a way that is not currently included in the common conceptual models used to assess NAPL-contaminated sites.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

Rate-limited mass transfer of octane, decane, and dodecane into nonionic surfactants solutions under laminar flow conditions

A key component to predicting the success of utilizing surfactants to enhance the removal of organic liquids from soil system is quantifying micellar solubilization kinetics. In this study, a flow reactor was employed to investigate the influence of surfactant ethoxylate chain length on the rates of solubilization of octane, decane, and dodecane in micellar solutions of a homologous series of purified dodecyl alcohol ethoxylates. Effluent concentration data were fit using a finite element model utilizing a linear-driving-force model to represent mass transfer at the interface. For flow rates between 0.1 and 2 ml min(-1), mass transfer coefficients ranged from 5 x 10(-8) to 7 x 10(-7)m s(-1) and did not vary in a systematic way with either solute structure or surfactant ethoxylate chain length and were lower than those found in pure water. Correlations developed for the Sherwood number based on diffusion coefficients of surfactant micelles containing organic material (organic-laden micelle) exhibit a velocity dependence similar to that found for systems based on aqueous diffusion. These results suggest that under gentle flowing conditions, the mass transfer is limited by diffusion of the organic-laden micelle. Although these trends are specific for this experimental system, the results demonstrate the importance of selecting the proper diffusion coefficient when modeling surfactant solubilization processes.     

Enhanced naphthalene solubility in the presence of sodium dodecyl sulfate: effect of critical micelle concentration

Surfactants can increase the solubility of non-polar compounds, and have been applied in areas such as soil washing and treatment of non-aqueous phase liquids (NAPLs). This investigation explored the feasibility of removing vapor phase polycyclic aromatic hydrocarbon (PAH) from gases using an anionic surfactant. The solubility of vapor phase naphthalene was measured herein using gas chromatograph (GC) with a photon ionization detector (PID). The measurement results indicated that surfactant molecules were not favorable to micelle formation when temperatures increased from 25 degrees C to 50 degrees C. Regardless of whether solutions were quiescent or agitated, equilibrium naphthalene apparent solubility increased linearly with surfactant concentrations exceeding critical micelle concentration (CMC). The pH effects on naphthalene apparent solubility were small. Agitation increased naphthalene apparent solubility and lumped mass transfer coefficients. Furthermore, lumped mass transfer coefficients decreased with increasing surfactant concentration owing to increase in interfacial resistance and viscosity and decreased spherical micelle diffusion coefficients. Finally, the net absorption rate increased because the solubilization effects of micelles exceeded the reduction effects of mass transfer coefficient above the CMC. The enhanced naphthalene apparent solubility from the addition of surfactant can be expressed by an enrichment factor (EF). The EF value of naphthalene for the surfactant solution at 0.1 M with agitation at 270 rpm relative to quiescent water could reach 18.6. This work confirms that anionic surfactant can improve the removal efficiency of hydrophobic organic compound (HOC) from the gas phase.     

Reaction-enhanced mass transfer and transport from non-aqueous phase liquid source zones

The objective of the following research is to theoretically quantify the enhancement of interphase mass transfer of dissolved non-aqueous phase liquid (NAPL) compounds from the non-aqueous phase to the aqueous phase and the enhancement of dispersive mass transport from a NAPL zone due to destruction of dissolved NAPL compounds. For relatively slow reaction rates, such as for permanganate and perchloroethene (PCE), local-scale mass transfer enhancement is expected to be small. Dispersive mass transport with reaction from a horizontal NAPL zone can be quantified using equations derived for a mathematically equivalent falling film reactor system. In contrast to local-scale interphase mass transfer, dispersive mass transport from NAPL zones may be significantly increased by reaction. Enhancement factors due to destruction of the NAPL compound(s) are mainly dependent on NAPL solubility and oxidant concentration and to a lesser extent on reaction rate, stoichiometry, and transverse dispersion coefficients. Higher NAPL solubility and/or lower oxidant concentration reduces the maximum expected enhancement factor. Reaction enhancement factors for mass transport from NAPL zones are expected to be in the range of 5-50 for permanganate and chlorinated solvents. Theoretical results suggest that assuming instantaneous reaction rates may be appropriate for dispersive mass transport from NAPL zones.     

Experimental Zn(II) retention in a sandy loam soil by very small columns

The use of column experiments, usually performed to better approximate field conditions, may provide information that is not available from batch experiments. In such experiments heavy metals are often adsorbed until saturation followed by desorption experiments. When the affinity of the metal to soil is high, the retention factor (R) could be greater than thousands and the duration of experiments can become impractically long. In order to use reasonable laboratory time, the flow rate should be increased or the column size decreased. The increase in flow rate produces undesirable kinetic and dispersion effects, so we used very small soil columns (pore volume = 0.31–0.70 ml) and relatively high flow rates (0.03–0.12 ml min⁻¹) in studies of Zn(II) adsorption and retention in soils. Conservative tracer flow column experiments under saturation conditions were carried out to determine flow parameters for different flow rates. Column pore volume (V_p), Peclet numbers (Pe) and longitudinal dispersion coefficients (D_L) were determined from breakthrough curves. The effect of type of electrolyte and ionic strength on the Zn(II) retention onto soil was determined. The influence of flow rate and bed height on the retention coefficient and on the mass transfer zone was also studied. The effect of different influent Zn(II) concentrations on the R values obtained was analyzed. Freundlich parameters from column experiments were compared with batch ones. The leaching efficiency of different electrolytes, salts of weak organic acids and EDTA was also studied.     

Modelling dissolution of multi-component non-aqueous phase liquids using Pareto distributions

A distributed-site non-aqueous phase liquid (NAPL) dissolution model that can capture the physical process of dissolution is presented. The model uses a Pareto distribution of effective interfacial area and a correlation for the mass transfer. Model parameters are obtained using a formal constrained non-linear optimisation technique. The model is applied on effluent data from two soil column experiments. The results of the Pareto distributed site model, with only two parameters to fit, were compared to two-site and multi-site models as well as the beta distributed site model. The Pareto distributed-site model produces equal or smaller error estimates than the other models.     

Parameters that control the cleanup of fractured permeable aquifers

This study develops a modeling approach for simulating and evaluating entrapped light nonaqueous-phase liquid (light NAPL-LNAPL) dissolution and transport of the solute in a fractured permeable aquifer (FPA). The term FPA refers to an aquifer made of porous blocks of high permeability that embed fractures. The fracture network is part of the domain characterized by high permeability and negligible storage. Previous studies show that sandstone aquifers often represent FPAs. The basic model developed in this study is a two-dimensional (2-D) model of permeable blocks that embed oblique equidistant fractures with constant aperture and orientation. According to this model, two major parameters govern NAPL dissolution and transport of the solute. These parameters are: 1) the dimensionless interphase mass transfer coefficient, K(f0), and 2) the mobility number, N(M0). These parameters represent measures of heterogeneity affecting flow, NAPL dissolution, and transport of the solute in the domain. The parameter K(f0) refers to the rate at which organic mass is transferred from the NAPL into the water phase. The parameter N(M0) represents the ratio of flow through the porous blocks to flow through the fracture network in regions free of entrapped NAPL. It also provides a measure of groundwater flow bypassing regions contaminated by entrapped NAPL. In regions contaminated by entrapped NAPL our simulations have often indicated very low permeability of the porous blocks, enabling a significant increase of the fracture flow at the expense of the permeable block flow. Two types of constitutive relationships also affect the rate of FPA cleanup: 1) the relationship between the saturation of the entrapped NAPL and the permeability of the porous blocks, and 2) the relationships representing effects of the entrapped NAPL saturation and the permeable block flow velocity on rates of interphase mass transfer. This study provides basic tools for evaluating the characteristics of pump-and-treat cleanup of FPAs by referring to sets of parameters and constitutive relationships typical of FPAs. The numerical simulations carried out in this study show that at high initial saturation of the entrapped NAPL, during initial stages of the FPA cleanup the contaminant concentration increases, but later it decreases. This phenomenon originates from significant groundwater bypassing the NAPL entrapped in the permeable blocks via the fracture network.