M S Mohan Kumar Mini Mathew & Shibani jha Dept. of Civil Engineering Indian Institute of Science
|
|
- Gavin Garrison
- 6 years ago
- Views:
Transcription
1 Multiphase Flow in Porous /Fracture Media M S Mohan Kumar Mini Mathew & Shibani jha Dept. of Civil Engineering Indian Institute of Science Bangalore, India
2 2 Multiphase fluids ----Fluids which are immiscible and is slightly soluble. Represent as individual phase in the subsurface, flow behavior is described as multiphase problem. Common trait of these substance are NAPLs 2 NAPLs LNAPLs DNAPLs 2 The phase does not fill the pore space completely. 2 Multiphase flow--- Simultaneous flow of more than two fluids, does not take place as a piston like process. 2 Flow experiments show that Darcy s law can be extended for multiphase flows 2 Darcy s velocity for each phase in a porous medium is described d by the generalized Darcy s law.
3 Source of NAPLs 2 Subsurface leakage of hydrocarbon fuels and other immiscible organic liquids due to leaky storage tanks or pipelines. 2 Coal tar from illuminating gas production, wastes from steel industry and wood treating operation 2 Organic substances used din industries i as Mineral lfuels ( ex : Petrol, fuel oil etc), Solvents Detergents ( ex : Chlorinated hydrocarbons) Goal 2 Estimate the potential danger of NAPL infiltration & to plan eventual remediation techniques. 2 To reduce the field investigation effort and cost.
4 2 Vertical migration in the vadose zone predominantly by gravity 2 Some lateral spreading due to capillar forces and media properties 2 Migration occurs when enough pressure is available to over come the displacement pressure 2 In saturated zone the movement is by displacement of water 2 In saturated medium phase system 2 In unsaturated medium phase system or phase system with static air pressure General migration pattern and process of NAPLs (after Helmig)
5 Processes of Multiphase Multicomponent System Water Water Air Organic compound Air Air Water vapour Organic compound Air NAPL Water NAPL
6 Multiphase Multicomponent Flow in the Subsurface - Governing Equations (φ α ρ α) L(φ α,ρα, v α) = + div(φ αρα v α) ρα qα ] dg= 0 t G k (φ S ρ X ) k α α α L(vX α, α ) = + divρ α α α α pmα α α α = t G k k k k { ρ v X φs D (ρ gradx ) } r ] dg 0 φs α D k pmα vi vj = δij αlv+ (α L α T) δijϕ v i j + 3/4 S 10/3 α D N T (u ρ S ) P L(uα, Pα, T) = s s α α α pm α = t t ρ α= 1 G [ { } { α α α α (1 φ)ρ c + ϕ + divv ρ (u + ) + div λ gradt } r ] 0 α
7 Interior Conditions Interior Conditions 1 S N 1 1 α S α N = = 1 1 α X α N k = = ) (, β α β α β α S k k P P C =P = ) ( α α α S k k r r
8 Present study consists of 2 Modelling and analysis of NAPLs migration in saturated porous medium----two phase NAPL-Water system 2 Study the influence of air phase on the infiltration of water in unsaturated porous medium ---- two phase Air-Water system. 2 Modelling and analysis of NAPL migration in unsaturated porous medium. Three phase Air-NAPL-Water system Two phase NAPL-Water system with constant air pressure 2 Modelling and analysis of NAPL migration in combined saturated-unsaturated porous medium 2 Parallel computation of multiphase flow in saturated porous medium
9 NAPL Infiltration Two phase region Air - Water Three phase region Water- NAPL-Air Two phase region Air - Water Interface region Interface region Interface region Single phase region of water Two phase region Water-NAPL Interface region Single phase region of water NAPL infiltration in the subsurface
10 Two Phase NAPL-Water Simulation in Saturated Porous Media 2 Develop a robust model and numerical method for the simulation of NAPL-Water in saturated porous media 2 Models of all the combinations of pressures and saturations of the phases have been developed 2 Comparative study is made between 6 models with simultaneous and modified sequential methods. 2 Comparative study between conventional simultaneous, modified sequential, and adaptive solution fully implicit modified sequential methods using pressure and saturation of wetting fluid are made 2 Effect of different types of approximations of nodal coefficient
11 2 Effect of different types of linearization methods in numerical modelling 2 Effect of different types of iterative methods in numerical modelling 2 Influence of capillarity and heterogeneity in ht heterogeneous media 2 Effect of Peclet, Courant numbers, and convergence criteria in two phase systems 2 Effect of different types of constitutive relations in the Effect of different types of constitutive relations in the numerical simulation
12 Solution Methodology used h Conventional simultaneous method h Modified sequential method h Adaptive solution fully implicit modified sequential method Simultaneous method Primary variables of the phases(pw & Sw) are solved together [ ] m [ ] m +1 A X = [ R] m
13 Sequential method Two step implicit technique Solve for pressure ( Saturation at previous iteration) [ m m + 1/2 A ] [ P ] = [ ] m W R1 1 = Solve for saturation ( Pressure at current level) [ ] m+ 1/2 [ ] m+ 1 [ ] m+ 1/ 2 A S = R 2 W 2
14 Adaptive solution fully implicit sequential method Identification of active and inactive nodes in each iteration for each primary variable. Difference in pressures between two consecutive iteration is greater than permissible level in any node or any of its 4 surrounding nodes. If the relative permeability of NAPL and water is greater than 0 and less than 1 in any node and any of its neighboring g4 nodes Difference in saturation between two consecutive iteration is greater than the permissible level in any node or any of its 4 surrounding nodes..
15 For active nodes Two step modified sequential method Solve for pressure ( Saturation at previous iteration) [ ] m [ ] m + 1/2 P [ R ] m A1 W = 1 Solve for saturation ( Pressure at current level) [ m+ 1/2 m+ 1 m+ 1/ 2 A ] [ S ] = [ R ] 2 W 2 For inactive nodes Pressure and saturation from the previous iteration
16 Verification of models Analytical solution based on McWhorter and Sunada(1990) Porosity of the medium [m 2 ] 0.35 Permeability 5E-11 Pore size distribution index 2 Displacement pressure [Pa] 2000 Residual wetting phase saturation 0.05 Wetting phase viscosity[pa-sec] 1E-03 Nonwetting phase viscosity[pa-sec] 0.5E-04 P -1/2 0t)=S 0 W (0,t) - P W (L,t) = A t - C t <= t 0 S W 0,t) W t>0
17 Verification of two dimensional models Laboratory experimental results(helmig, 1998) sand 1 sand 2 Intrinsic permeability [m 2 ] : 6.64E E-12 Porosity : Wetting Phase residual saturation : Displacement pressure : Pore size distribution index : 2.7 2
18 NAPL saturation distribution at different times T = 1000sec Experimental result T = 3000sec Present Model results T = 5000sec
19 Comparison of different types of models Between models Between numerical methods
20 Nonwetting fluid distribution profile after 2400 sec Density of NAPL : 1460kg/m 3 Viscosity of NAPL : pa-sec Experimental results Present model
21 Sand 1 Sand 2 Sand 3 Sand 4 Φ k 5.04E E E E S wr λ Pd
22 Sensitivity study for SIM, MSEM and ASMSEM
23 Nonwetting fluid distribution Host permeability : 7.08E-12m 2 Density of TCE : 1460kg/m 3 Host displacement pressure : Pa Porosity : 0.34 Residual wetting phase saturation : Pore size distribution index : 2.48 Lens displacement pressure : Pa
24 Host permeability : 9.28E-12 m 2 Density of TCE : 1400kg/m 3 Lens permeabliity : 7.53E-12 m 2 Viscosity of TCE : 1E-03 Pa-sec Host displacement pressure : 1667Pa Residual saturation : 0.12 Lens displacement pressure : 2353Pa Porosity : 0.38 Pore size distribution index : 2
25 Sand 1 Sand 2 k 6.64E E-12 Φ S Wr Pd λ ρ NW 1460kg/m 3 μ NW 0.9e-03pa-s
26 Effect of different types of approximations at the cell faces a Arithmetic mean or Harmonic mean b Harmonic mean with upstream weighing of relative permeabilities c Harmonic mean with upstream weighing of relative permeability and capillarity diffusivities d. Fully upwind method
27 Effect of different types of linearizations in the two phase system Classical Newton Raphson method Modified Newton Raphson method Picard s method
28 Effect of classical Newton Raphson method
29 Influence of capillarity and heterogeneity in two phase simulation J( s P C P P + * e ) = Pe = σ φ k 1/2 [ ] J ( S ) L ( W Pc S* - + P d - P d + 0 Wetting fluid saturation 1
30 Effect of without and with heterogeneity and capillary pressure effect
31 Identification of robust iterative method for the two phase simulation Explicit and alternating direction implicit methods are unstable (Peaceman, 1967) Incomplete Cholesky conjugate gradient method(iccg) Hesteness and Stiefel s conjugate gradient method(cghs) Diagonal scaling conjugate gradient(dscg) Krylov subsurface solver BiCGSTAB Strongly Implicit Procedure(SIP)
32 Influence of Peclet, Courant numbers and convergence criteria
33 Effect of different Constitutive Relations used Brooks and Corey Van-Genuchten
34 Effect of NAPL migration in Random Heterogeneous Media
35 Two Phase Air Water Flow Simulation in Unsaturated Porous Media 2 To study the effect of air flow for infiltration and distribution of water in unsaturated zone 2 Development of Numerical models of one dimensional and two dimensional using SIM, MSEM, and ASMSEM 2 Validation of one phase, quasi two phase and two phase models 2 Air water simulation in one dimensional homogeneous porous medium 2 Air water simulation in one dimensional heterogeneous porous media 2 Air water simulation in two dimensional heterogeneous porous media One phase model - Richard s equation Quasi two phase model - Two phase air water system with constant air pressure Two phase air water model - Air and water moves simultaneously
36 Verification of quasi two phase model using one phase model One phase model - Richard s equation Quasi two phase model - Arithmetic approximation - Harmonic approximation with upstream nodal relative permeabilities Flux = 3.29m/day Initial pressure = m
37
38 Validation of two phase Air-Water model with experimental results(touma and Vaucline, 1986) Soil parameters are 93.5cm Saturated water content = Residual moisture content = van Genuchten parameters α = 0.044cm -1 β = 2.2 Gas constant of air = Molecular weight = Impermeable boundary
39 Flux greater than hydraulic conductivity - 20cm/h Bounded column Two phase model One phase model
40 Open column Two phase model One phase model
41 Flux less than saturated hydraulic conductivity - 8cm/hr Two phase model One phase model
42 Soil Properties Residual moisture content : Constant head van-genuchten parameters : α= 0.044cm 044cm -1, β = Saturated moisture content : Porosity : 0.37 Intrinsic permeability :0.0145cm 2 P W =0.906 inch 3.07 ft
43 Two phase Air-Water in heterogeneous porous media Initial pressure = -100cm Flux at the inflow end = 9.5E-05cm/sec 150cm lens 20 cm 40cm Impermeable layer
44 Sand and sandy loam Open column One phase model Two phase model
45 Bounded column One phase model Two phase model
46 Sand and Silt Open column One phase model Two phase model
47 Bounded column One phase model Two phase model
48 Sand and Clay
49
50 Two dimensional two phase air water simulation in heterogeneous media Initial pressure = -150cm Flux = 1.13x cm/sec 10cm sand 50cm Low permeability soil 60cm 15cm 20cm 100 cm Test One phase Two phase Two phase Two phase problem model MSEM SIM ASMSEM
51 Sand and Sandy loam
52 Modelling and Simulation of NAPL Migration in Saturated - unsaturated us u zone oe 2 Development of three phase model of water, NAPL, and air 2 Quasi ithree phase model - Air - NAPL - Wt Water with static air pressure 2 Fully three phase model - Air- NAPL - Water with air pressure 2 Validation of the present models 2 Effect of capillary pressure and heterogeneity in the numerical modelling of three phase systems 2 Comparative study between simultaneous, sequential and adaptive modified sequential methods for the saturated unsaturated systems 2 Effect of air pressure on the numerical modelling of three phase systems
53 Model Validation Oil infiltration 3 cm 15 cm 20 cm 30 cm Watertable Infiltration : 30ml of NAPL, 20,5,5----1hr interval -----case a : 30ml ml for 3hrs case b
54 NAPL profile after 3hrs
55
56 Air-NAPL-Water simulation in variably saturated media Infiltartion rate- 5.07E-05m/sec
57
58 Comparative study between simultaneous, modified sequential and adaptive modified sequential methods. SIM MSEM ASMSEM Quasi three phase Full three phase Quasi three phase Full three phase
59 Effect of capillarity and heterogeneity in the numerical modelling of three phase systems LNAPL Flux = 7.1x10-4 m/sec
60 DNAPL
61 Effect of with and without air pressure LNAPL Open column
62 LNAPL Bounded column
63 Open column DNAPL
64 DNAPL Bounded column
65 Parallel computation of NAPL migartion in saturated porous media Need of Parallel Computation To reduce the computational time To simulate bigger domain multiphase flow simulations which is not possible by conventional computation General Parallel Architecture - Interaction between the processors Shared memory programming g Distributed memory with Message Passing Interface programming Data parallel programming
66 Parallel machine used for computing 32 nodes Scalable ab Parallel(SP2) a system Nodes 1-8 : IBM RS 66MHz with 512MB RAM Nodes 9-24 : IBM RS 77MHz with 256MB RAM Nodes : IBM RS 133MHz with 256MB RAM Distributed memory machine with Message Passing Interface Programming(MPI) Parallelization Methodology Divide the domain into sub-domains Identify the initial and boundary conditions of each sub-domain Compute the primary variables(p W &S W ) Communicate primary variables between sub-domains using Message Passing Interface (MPI) programming
67 Domain Decomposition methods used Overlapping region Row wise Column wise
68 1.0m Low permeability soil Source 1.5m Impermeable boundary Flux at source = 5.16E-06 m/s m m
69
70 Properties of media and fluid Properties Sand 1 Sand 2 Porosity Displacement pressure[pa] Intrinsic permeability [m 2 ] 5.04E E-12 Residual saturation pore size distribution index Density of NAPL[kg/m 3 ] 1460 Viscosity[Pa-s ] 0.90E-03
71 Nonwetting fluid saturation distribution after 1day after 3days
72 Modelling of Multiphase Flow in a Fracture Media
73 Purpose of the Study A numerical model to simulate multiphase flow within a fracture zone To study the conditions under which a DNAPL can enter a rough walled, initially water saturated fracture To study the subsequent behaviour of DNAPL within the fracture A rough walled fracture defined in terms of aperture distribution A single pair of parallel plates (homogeneous) A set of parallel plate pairs (spatially correlated) Model study for 1D variable aperture fracture Extended for 2D homogeneous and variable aperture fracture The model has been studied for various aperture e distribution with different correlation length Model study for sensitivity to fluid and fracture properties on the migration rate of DNAPL through fractures.
74 Model Conceptualization Capillary yp pressure (Bear, 1972) P c =P nw -P w Parallel plate P e =(2σcosθ)/e Circular P e =(4σcosθ)/e Entry for DNAPL P >P P c > P e H d =(2σ)/(Δρge)
75 Height of Pool versus Aperture Invaded 10 2 σ = N/m ( μm ) FRACTURE APE ERTURE INVADED ( σ = N/m σ = N/m σ = N/m σ = 005 N/m HEIGHT OF DNAPL POOL (m)
76 Mathematical Equations Mass conservation - (ρ w q wi e) x i =( (φs w ρ w )/t) e - (ρ nw q nwi e) x i = ( (φ S nw ρ nw ) t)e Darcy s law q wi = - (k ij k rw /μ w )( P w / x j +ρ w g y/ x j ) q nwi = - (k ij k rnw /μ nw ) ( P nw / x j +ρ nw g y/ x j ) Continuity equations ( / x i )[(ek ij k rw /μ w )( P w / x j + (ρ w g) ( y/ x j ))] = eφ( S w / t) ( / x i )[(ek ij k rnw /μ nw )( P nw / x j + (ρ nw g) ( y/ x j ))] = eφ( S nw / t)
77 Continuity Equations in Terms of P w and S w as Dependent Variables Saturation constraint S w +S nw = 1.0 For wetting phase ( / x)[(ekk rw /μ w )( P w / x)] + ( / y)[(ekk rw /μ w ) ( P w / y+ρ w g)] = eφ( S w / t), For Non-Aqueous phase ( / x)[(ekk rnw /μ nw )( (P w+ P c )/ x)] + ( / y)[(ekk rnw /μ nw ) ( (P ( w + P c )/ y + ρ nw g)] = - eφ( S w / t),
78 Boundary Conditions Constant mass flow rate or constant pressure-saturation Neumann boundary with zero flux Dirichlet boundary Dirichlet boundary Neumann boundary with zero flux a) Dirichlet boundary condition: u 1 = f 1 and u 2 = f 2 b) Neumann boundary condition: (k k /μ )( P / x +ρ g y/ x )=f - (k ij k rw /μ w ) ( P w / x j +ρ w g y/ x j ) = f 3 -(k ij k rnw /μ nw ) ( P nw / x j +ρ nw g y/ x j ) = f 4
79 Relative and Fracture Permeabilities Brooks and Corey model k rw =S (2+3λ)/λ e k rnw =(1-S e ) 2 (1-S e (2+3λ)/λ ) Flow between parallel plates F wi = - (e 3 /12μ w )( P w / x j + ρ w g y/ x j ) S e =(S w -S r )/(1-S r ), 0<=S e <=1 F nwi = - (e 3 /12μ nw ) ( P nw / x j + ρ nw g y/ x j ) S e =(P c /P d ) -λ Fracture permeability k = e 2 /12 P d= P e
80 One Dimensional Model Equivalent porous media approach k=(e 2 )/12( R 1.5 r ) (Marshily, 1986) For moderately rough fracture plane R r = 0.1 Solution domain
81 Fluid and Fracture Properties Properties Values Units Wetting phase viscosity Pa.s Nonwetting phase viscosity Pa.s Wetting phase density Kg/m 3 Nonwetting phase density Kg/m 3 Interfacial tension N/m Porosity Pore size distribution index Wetting phase residual saturation 0.1 -
82 Model Verification Pool Height (meters) present model for density = 1460 kg/m 3 Kueper s model for density = 1460 kg/m 3 present model for density = 1200 kg/m Kueper s model for density = 1200 kg/m Aperture (microns) present model for density = 1460 kg/m 3 Kueper s model for density = 1460 kg/m 3 present model for density = 1200 kg/m 3 Kueper s model for density = 1200kg/m Time (hours) DNAPL pool vs time Time (hours) Aperture vs time Kueper s Result Present Model es) Fracture Dip (degree Time (hours) Fracture dip vs time
83 Sensitivity with DNAPL Pooled Above Fracture Opening DNAPL pool height = 0.5 m DNAPL pool height = 1.0 m DNAPL pool height = 1.5 m DNAPL pool height = 2.0 m DNAPL pool = 0.50 m DNAPL pool = 0.40 m DNAPL pool = 0.30 m DNAPL pool = 0.35 m Fracture Ap perture (microns) DNAP PL Saturation Time (hours) Length of Fracture (m) Fracture aperture vs travel time DNAPL saturation at t = 15000s
84 Sensitivity with aperture of Fracture Opening Fracture Aperture = 100 μm 0.45 Fracture Aperture = 75 μm Fracture Aperture = 50 μm Fracture Aperture = 25 μm DNAPL Poo ol Height (meters) DNAP PL Saturation e = 25 μm e = 50 μm e = 75 μm e = 100 μm Time (hours) Length of Fracture (m) DNAPL Pool vs travel time DNAPL saturation at t = 15000s
85 Sensitivity With Fracture Dip Dip = 45 o Dip = 30 o Dip = 15 o Dip = 0 o 0.35 Dip = 60 o APL Saturation DN Dip = 90 o Length of Fracture (m) DNAPL saturation at t = 15000s
86 Effect of Viscosity and Density Fracture Dip = 90 degree Aperture = 75 μ m DNAPL Pool = 0.5 m Viscosity = 0.57E 03 Time = 1000 seconds Time = 2000 seconds Time = 5000 seconds Time = seconds Fracture Dip = 90 degree Aperture = 75 μ m DNAPL Pool = 0.5 m DNAPL Density = 1460 kg/m 3 Time = 1000 seconds Time = 2000 seconds Time = 5000 seconds Time = seconds DNAPL Saturation DNAPL Saturation Length of Fracture (meters) Length of Fracture (meters) Fracture Dip = 90 degree Aperture = 75 μ m DNAPL Pool = 0.5 m Viscosity = 0.9E 03 Time = 1000 seconds Time = 2000 seconds Time = 5000 seconds Time = seconds Fracture Dip = 90 degree Aperture = 75 μ m DNAPL Pool = 0.5 m DNAPL Density = 1600 kg/m 3 Time = 1000 seconds Time = 2000 seconds Time = 5000 seconds Time = seconds uration DNAPL Satu aturation DNAPL Sa Length of Fracture (meters) Length of Fracture (meters) DNAPL saturation plots for viscosity effect DNAPL saturation plots for density effect
87 Spatially Correlated Aperture Field for One Dimensional Fracture Aperture (mm) Length of Fracture (meters) Mean aperture e 75 μm Std.dev σ Correlation length 0.2m Maximum e μm Minimum e 19.1 μm Aperture distribution along the length of fracture
88 DNAPL and Water Migration T=5000s T=10000s T=25000s T=50000s DNAPL Velocity (m/s) DNAP PL Saturation Length of Fracture (meters) DNAPL velocity T=5000s T=10000s T=25000s T=50000s Time = 5000 s 0.3 Time = s Time = s Time = s Length of Fracture (meters) s) Water Velocity (m/s DNAPL distribution ib ti at various times Length of Fracture (meters) Water velocity
89 Effect of DNAPL Pool and Fracture Dip in a Correlated Aperture Field L Saturation DNAP DNAPL Pool = 0.15m DNAPL Pool = 0.20m DNAPL Pool = 0.25m DNAPL Pool = 0.50m DNAPL Pool = 0.75m DNAPL Saturation Fracture Dip = 90 o Fracture Dip = 60 o Fracture Dip = 45 o Fracture Dip = 30 o Fracture Dip = 15 o Fracture Dip = 1 o Length of Fracture (meters) Length of Fracture (meters) DNAPL distribution ib ti at t = 50000s DNAPL distribution ib ti at t = 50000s for various pool height for various fracture dip
90 Effect of Correlation Length on DNAPL and Water Migration Fracture Aperture (mm) Nodal Discretization = 250 mm correlation length = 100 mm correlation length = 200 mm correlation length = 300 mm correlation length = 400 mm DNAPL Velocity (m/s) correlation length =200mm correlation length =300mm correlation length=400mm correlation length =100mm Length of Fracture (meters) Aperture distribution Length of Fracture (meters) DNAPL velocity at t = 50000s on DNAPL Saturatio correlation length = 200mm correlation length = 300mm Water Velocity (m m/s) correlation length = 400mm correlation length = 300mm correlation length = 200mm correlation length = 100mm 0.3 correlation length = 400mm correlation length = 100mm Length of Fracture (meters) DNAPL distribution at t = 50000s Length of Fracture (meters) Water velocity at t = 50000s
91 Range of Apertures for Various Fracture Aperture fields generated Correlation Maximum Minimum Aperture at Aperture just lengths mm aperture μm aperture μm the top of fracture μm below the top cell μm
92 Two Dimensional Model DNAPL Saturation Time = 8000 s Time = s Time = s Time = s Length of fracture along y axis (m) DNAPL saturation along vertical centre line Saturation of DNAPL 0.5 Solution domain for homogeneous fracture plane Time = 8000 s Time = s Time = s Time = s n DNAPL Saturatio Length of fracture along x axis (m) DNAPL saturation along top boundary
93 ture along y axis (m) Length of Fract 0.35 Time = 8000 s Length of Fracture along x axis (m) Length of Fracture along y axis (m) DNAPL Velocity vector plot Time = 8000 s Length of Fracture along x axis (m) Length of Fractu ure along y axis (m) 0.35 Time = s Length of Fracture along x axis (m) along y axis (m) Length of Fracture DNAPL Velocity vector Time = s Length of Fracture along x axis (m) Length of Fracture along y axis (m) 0.35 Time = s Length of Fracture along x axis (m) g y axis (m) Length of Fracture alon DNAPL Velocity vector plot Time = s Length of Fracture along x axis (m) DNAPL distribution at various times DNAPL velocity vector at various times
94 Length of Frac cture along y axis (m) 0.35 Fracture Aperture = 25 μm Length of Fracture along x axis (m) along y axis (m) Fracture Aperture = 35 μm Base case Length of Fracture Length of Fracture along x axis (m) Le ength of Fracture alon ng y axis (m) Horizontal Fracture Effect of aperture Length of Fracture along x axis (m) Effect of dip
95 DNAPL Courant Number
96 Effect of Heterogeneity Solution domain for heterogeneous fracture plane DNAPL Saturation T = 8000s T = 30000s T = 40000s turation DNAPL Sat Time = 8000 s Time = s Time = s uration DNAPL Satu Length of Fracture along y axis (m) DNAPL saturation along vertical centre line Length of Fracture along x axis (m) DNAPL saturation along top boundary
97 Effect of Heterogeneity on DNAPL Migration
98 Grid Convergence
99 Rough Walled Fracture Plane in Terms of Spatially Correlated Aperture Distribution Solution domain for a rough walled fracture plane Aperture in mm Aperture distribution for a rough walled fracture plane
100 Fluid Properties and Source Condition Properties Units Value Wetting phase density Kg/m DNAPL density Kg/m Wetting phase viscosity Pa.s DNAPL viscosity Pa.s Interfacial tension N/m Wetting phase pressure at source Pa 0.0 Wetting phase saturation at source - 0.5
101 Range of Apertures in the Field Generated Along the top Along the bottom 4.5 x Along the top Along the bottom 3.5 riation of Aperture (mm) Var Perm meability of Fracture (m 2 ) Length of Fracture along x axis (m) Length of Fracture along x axis (m) Aperture distribution along the top and bottom boundary Fracture permeability along the top and bottom boundary Range of apertures within a fracture Range Unit Domain Top boundary Bottom boundary Maximum μm minimum μm
102 DNAPL Migration Length of Fracture along y axis (m) Time = 5000s Length of Fracture along x axis (m) ength of Fracture along y axis (m) L Time = seconds Length of Fracture along x axis (m) DNAPL distribution at t = 5000s DNAPL distribution at t = 50000s Length of Fractu ure along y axis (m) Time = seconds Length of Fracture along x axis (m) Length of Fractu ure along y axis (m) Time = seconds Length of Fracture along x axis (m) DNAPL distribution at t = 10000s DNAPL distribution at t = s
103 Isotropic and Anisotropic Aperture Field Aperture (mm) 0.11 Fracture along y axis (mm) Length of Length of Fracture along x axis (mm) 0.02 Solution domain Aperture (mm) Correlation length (30,30)mm30) Aperture (mm) Length of Fracture along y axis (m mm) s (mm) Length of Fracture along y axis Length of Fracture along x axis (mm) Correlation length (30,20)mm Length of Fracture along x axis (mm) Correlation length (20,20)mm 0.02
104 Anisotropic Aperture Field 0.7 Length of Frac cture along y axis (mm m) Length of Fracture along x axis (mm) 0 Isotropic Aperture Field 0.7 Isotropic Aperture Field Le ength of Fracture along y axis (mm) Le ength of Fracture along y axis (mm) Length of Fracture along x axis (mm) Length of Fracture along x axis (mm) 0
105 Buoyancy flow and Capillary Trapping Local capillary trapping Gravity fingering Buoyancy flow -most dense fluid- downward, lateral and sliding movement Gravity fingering - once the residual DNAPL pools up on the top of the wedge
106 Observations Fractures provides preferential and faster pathways DNAPL enters the fracture at the points of largest aperture and continue to migrate through the larger aperture regions The certain regions of the fracture may remain void of DNAPL at all times The ability of DNAPL to enter smaller aperture regions of fracture increases as a function of fdepth of penetration Traverse time for DNAPL is inversely proportional to the fracture aperture, fracture dip and DNAPL pooled above the fracture Fracture aperture is most sensitive parameter Shallow fractures (30 o to0 o ) shows significant change in migration For aperture field with correlation length close to grid size, DNAPL migration shows homogeneous nature Anisotropic distribution of aperture field provides higher rate of DNAPL movement
107 Immiscible-miscible miscible flow in coastal reservoirs Seawater intrusion is an ideal problem to study the buoyancy flow Natural intrusion of heavier fluid with matched viscosity
108 Ghyben-Herzberg approximation for seawater intrusion Groundwater flow patterns in an idealized, homogeneous coastal aquifer. (Source: USGS) Analytical solution on the basis of Ghyben-herzberg Relationship (Bear,1972) and Single Potential Theory (strack,1976). The assumption of a sharp interface between freshwater and saltwater. Variable density flow in both time and space dimension (Frind,1982; Voss and Souza; Kolditz et al.,1998). 2/10/
109 Immiscible-miscible flow of Seawater - freshwater A freshwater confined aquifer of size 2m by 1m, Impermeable boundary conditions on the top and bottom. The left open boundary, a constant freshwater influx (q f =6.6x10-5 m 2 /s, ρ f =1000 kg/ m 3, μ f =10-3 m 2 /s) The open right hand boundary, a hydrostatic pressure distribution (ρ s =1025 kg/ m 3, μ s =10-3 m 2 /s) 50 % saltwater saturation (S S ) is assumed 25 % of the vertical dimension of the aquifer is assumed to be freshwater mouth. The interface between the 2 fluids is sharp without any smearing Seawater saturation contours at 6000 seconds
110 Miscible and immiscible models y dimension in (m) Immiscible model for linear retention laws x dimension in (m) Mixing accelerates the front movement Immiscible-Miscible model for linear retention laws 110
111 Miscible and immiscible models Immiscible model for linear retention ti laws Mixing accelerates the front movement Immiscible-Miscible model for linear retention laws 111
112 With buoyancy Effect of freshwater flux Without buoyancy Effect of buoyancy Buoyancy in a multiphase model acts as a dispersive mechanism Freshwater flux reduced to half
113 Effect of relative permeability on the seawater-freshwater transition zone Brooks - Corey Linear relation Quadratic relation Cubic relation shows nearest resemblance to Brooks - Corey Cubic relation
114 Effect of relative permeability on the seawater-freshwater transition zone
115 3-Phase Immiscible Flow in Coastal Reservoir (Buoyancy Flow and Capillary Trapping) Seawater intrusion and oil migration in a coastal aquifer This study is done to demonstrate buoyancy flow and capillary trapping in immiscible flow
116 Conceptual Model The pore space is completely filled by all the three existing phases S w + Sn + Ss =1 P P = n P n P s P w P s P cnw = w P = P = P P cns cnw cns csw P csw Interface conditions Capillary equilibrium This leads to capillary trapping 1. Possibility of interface curvature: DNAPL most nonwetting among three 2. Possibility of capillary trapping (fingering) w - Freshwater (most wetting phase) n - DNAPL (most non-wetting phase) s- Seawater (intermediate phase)
117 P P cns cnw = P = P d d ( Ss + Sn Srn) ( 1.0 Srn Srw Sw + Sn Srw S (1.0 S rn S ) ( rn) ( rw ) Linear relation for capillary pressure (conditions for the interfaces) Non-linear relation (Brooks-Corey) for relative permeability of phases (conditions for the interfaces) k rs = k ( ) 4 rw = S ew ( S ) 3 es (2.0 Ses ) ( S ) 3 (1.0 S S ) k = + rn en es ew S ew = ( Sw Srw ) ( 1.0 S S ) S es = rw S rn S ( 1.0 S S ) S en = s rw S n rn S ( 1.0 S S ) rw rn rn
118 3-Phase buoyancy flow and capillary trapping y dimension in (m) Saltwater wedge moved seaward 2. Some local blobs (trapping) and fingering developed in three-phase region x dimension in (m) y dire ection in (m) Initial condition of saline aquifer x direction in (m)
119 Buoyancy flow and Capillary Trapping Local capillary trapping Gravity fingering Buoyancy flow -most dense fluid- downward, lateral and sliding movement Gravity fingering - once the residual DNAPL pools up on the top of the wedge
120 Study of miscible and non-isothermal flow models Fracture plane developed due to external stresses Idealized fracture plane P e =(2σcosθ)/e k = 2 e 12 δ n 1 = k n 2 2 [ σ cos β + σ β ] 1 3 sin Brady and Brown, 1993 Normal deformation components of a joint
121 Coupling Deformation with Fluid Pressures Assumption δ n 1 = 1 3 sin β k n [ σ β + σ P ] 2 2 cos Modified normal df deformation If P >P If P w >P n w If P w <P n δ 1 [ σ β + σ 2 2 cos β ] n = 1 β 3 sin β P n kn Coupling Deformation with Immiscible flow Model Aperture of fracture, e t e t = e 0 ± δ 121 n
122 Influence of deformation on phase distribution Deformation of the medium enhances the diffusion process Phase distribution is more diffusive
123 Influence of deformation on energy transfer Deformation of the medium enhances the diffusion process
124 Influence of deformation on mass transfer Deformation of the medium enhances the diffusion process Mass transfer diffuses very fast In the storage system this kind of diffusive process secures better storage
125 Immiscible flow with miscible, non-isothermal and medium deformation This test problem shows that immiscible flow system should be coupled with miscible and non-isothermal flow which can represent the complete multiphase system Also the medium deformation through hydromechanical coupling should be considered if the multiphase flow system has to beconsidered d in any geological environment, shallow or deep
126 Thank You
DNAPL migration through interbedded clay-sand sequences
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield. UK. June 2001). IAHS Publ. no. 275. 2002. 455 DNAPL
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationMixed-hybrid finite element method for modelling two-phase flow in porous media
Journal of Math-for-Industry, Vol. 3 (211C-2), pp. 9 19 Mixed-hybrid finite element method for modelling two-phase flow in porous media Radek Fučík and Jiří Mikyška Revised on September 23, 211 Abstract.
More informationAnalysis of Multiphase Flow under the Ground Water
Analysis of Multiphase Flow under the Ground Water Pramod Kumar Pant Department of Mathematics, Bhagwant University, Ajmer, Rajasthan, India Abstract The single-phase fluid flow through a porous medium
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationEvaporation-driven soil salinization
Evaporation-driven soil salinization Vishal Jambhekar 1 Karen Schmid 1, Rainer Helmig 1,Sorin Pop 2 and Nima Shokri 3 1 Department of Hydromechanics and Hydrosystemmodeling, University of Stuttgart 2 Department
More informationNon-equilibrium alcohol flooding model for immiscible phase remediation: 2. Model development and application
Advances in Water Resources Vol. 1, No. 8, pp. 663 678, 1998 1998 Elsevier Science Ltd Printed in Great Britain. All rights reserved PII: S 0 3 0 9-1 7 0 8 ( 9 7 ) 0 0 05-0 0309-1708/98/$19.00 + 0.00 Non-equilibrium
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationAppendix D Fractured Rock Appendix
Appendix D Fractured Rock Appendix 1.0 Introduction The behavior of LNAPL in fractured bedrock is not necessarily intuitive and is not as easily described using the principles and techniques adopted for
More informationField Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer. Bo Ren, Larry Lake, Steven Bryant
Field Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer Bo Ren, Larry Lake, Steven Bryant 2 nd Biennial CO 2 for EOR as CCUS Conference Houston, TX October 4-6, 2015
More informationGG655/CEE623 Groundwater Modeling. Aly I. El-Kadi
GG655/CEE63 Groundwater Modeling Model Theory Water Flow Aly I. El-Kadi Hydrogeology 1 Saline water in oceans = 97.% Ice caps and glaciers =.14% Groundwater = 0.61% Surface water = 0.009% Soil moisture
More informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationDNAPL MIGRATION IN SINGLE FRACTURES:
DNAPL MIGRATION IN SINGLE FRACTURES: ISSUES OF SCALE, APERTURE VARIABILITY & MATRIX DIFFUSION By Katherine I. Hill A thesis submitted to the, Faculty of Engineering, Computing and Mathematics, in conformity
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
More informationCHAPTER III BEHAVIOR OF HYDROCARBONS IN THE SUBSURFACE
CHAPTER III BEHAVIOR OF HYDROCARBONS IN THE SUBSURFACE CHAPTER III BEHAVIOR OF HYDROCARBONS IN THE SUBSURFACE The purpose of this chapter is to supplement your knowledge of hydrocarbon behavior in the
More informationψ ae is equal to the height of the capillary rise in the soil. Ranges from about 10mm for gravel to 1.5m for silt to several meters for clay.
Contents 1 Infiltration 1 1a Hydrologic soil horizons...................... 1 1b Infiltration Process......................... 2 1c Measurement............................ 2 1d Richard s Equation.........................
More informationDarcy s Law, Richards Equation, and Green-Ampt Equation
Darcy s Law, Richards Equation, and Green-Ampt Equation 1. Darcy s Law Fluid potential: in classic hydraulics, the fluid potential M is stated in terms of Bernoulli Equation (1.1) P, pressure, [F L!2 ]
More informationSimulation of Unsaturated Flow Using Richards Equation
Simulation of Unsaturated Flow Using Richards Equation Rowan Cockett Department of Earth and Ocean Science University of British Columbia rcockett@eos.ubc.ca Abstract Groundwater flow in the unsaturated
More informationChapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:
Chapter Seven Horizontal, steady-state flow of an ideal gas This case is presented for compressible gases, and their properties, especially density, vary appreciably with pressure. The conditions of the
More informationTwo-phase gravity currents in porous media
J. Fluid Mech. (211), vol. 678, pp. 248 27. c Cambridge University Press 211 doi:1.117/jfm.211.11 Two-phase gravity currents in porous media MADELEINE J. GOLDING, JEROME A. NEUFELD, MARC A. HESSE AND HERBERT
More informationUpdating the Coupling Algorithm in HYDRUS Package for MODFLOW
Updating the Coupling Algorithm in HYDRUS Package for MODFLOW SAHILA BEEGUM Guided by Dr. K P Sudheer, Dr. Indumathi M Nambi & Dr. Jirka Šimunek Department of Civil Engineering, Indian Institute of Technology
More informationINL Capabilities and Approach to CO 2 Sequestration. 4 th U.S.-China CO2 Emissions Control Science & Technology Symposium
www.inl.gov INL Capabilities and Approach to CO 2 Sequestration 4 th U.S.-China CO2 Emissions Control Science & Technology Symposium Travis McLing, Ph.D. Energy Resources Recovery and Sustainability September
More information1. Water in Soils: Infiltration and Redistribution
Contents 1 Water in Soils: Infiltration and Redistribution 1 1a Material Properties of Soil..................... 2 1b Soil Water Flow........................... 4 i Incorporating K - θ and ψ - θ Relations
More informationI. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010
: an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS
More informationDavid de Courcy-Bower and Samuel Mohr
Applicability and Limitations of LNAPL Transmissivity as a Metric within Bedrock Formations Insert then choose Picture select your picture. Right click your picture and Send to back. David de Courcy-Bower
More informationAnalysis of oil displacement by water in oil reservoirs with horizontal wells
Analysis of oil displacement by water in oil reservoirs with horizontal wells Paulo Dore Fernandes, Thiago Judson L. de Oliveira and Rodrigo A. C. Dias Problem Description This work involves near-well
More informationAvailable online at Energy Procedia 1 (2009) (2008) GHGT-9
Available online at www.sciencedirect.com Energy Procedia 1 (2009) (2008) 3331 3338 000 000 Energy Procedia www.elsevier.com/locate/procedia www.elsevier.com/locate/xxx GHGT-9 Application of gravity currents
More informationNumerical Simulations of Two-Phase Flow in Rigid Porous Media
University of Colorado, Boulder CU Scholar Civil Engineering Graduate Theses & Dissertations Civil, Environmental, and Architectural Engineering Spring 1-1-2017 Numerical Simulations of Two-Phase Flow
More informationAn Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD
An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface
More informationDistribution of pore water pressure in an earthen dam considering unsaturated-saturated seepage analysis
E3S Web of Conferences 9, 194 (16) DOI: 1.11/ e3sconf/169194 E-UNSAT 16 Distribution of pore water in an earthen dam considering unsaturated-saturated seepage analysis 1a Kumar Venkatesh, Siva Ram Karumanchi
More informationReservoir Geomechanics and Faults
Reservoir Geomechanics and Faults Dr David McNamara National University of Ireland, Galway david.d.mcnamara@nuigalway.ie @mcnamadd What is a Geological Structure? Geological structures include fractures
More informationPore-scale modeling extension of constitutive relationships in the range of residual saturations
WATER RESOURCES RESEARCH, VOL. 37, NO. 1, PAGES 165 170, JANUARY 2001 Pore-scale modeling extension of constitutive relationships in the range of residual saturations Rudolf J. Held and Michael A. Celia
More information6. THE BOREHOLE ENVIRONMENT. 6.1 Introduction. 6.2 Overburden Pressures
6. THE BOREHOLE ENVIRONMENT 6.1 Introduction Wireline logging has a single clearly defined purpose: to give accurate and representative data on the physical properties of the rock formations and fluids
More informationpressure-saturations relationships for two-phase flow in 3D porous media: implications micro-heterogeneities
Loughborough University Institutional Repository Dynamic effects in capillary pressure-saturations relationships for two-phase flow in 3D porous media: implications of micro-heterogeneities This item was
More informationModelling of pumping from heterogeneous unsaturated-saturated porous media M. Mavroulidou & R.I. Woods
Modelling of pumping from heterogeneous unsaturated-saturated porous media M. Mavroulidou & R.I. Woods Email: M.Mavroulidou@surrey.ac.uk; R. Woods@surrey.ac.uk Abstract Practising civil engineers often
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationRATE OF FLUID FLOW THROUGH POROUS MEDIA
RATE OF FLUID FLOW THROUGH POROUS MEDIA Submitted by Xu Ming Xin Kiong Min Yi Kimberly Yip Juen Chen Nicole A project presented to the Singapore Mathematical Society Essay Competition 2013 1 Abstract Fluid
More informationInstructor : Dr. Jehad Hamad. Chapter (7)
Instructor : Dr. Jehad Hamad Chapter (7) 2017-2016 Soil Properties Physical Properties Mechanical Properties Gradation and Structure Compressibility Soil-Water Relationships Shear Strength Bearing Capacity
More informationAvailable online at ScienceDirect
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (17) 3343 3352 www.elsevier.com/locate/procedia 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13,
More informationInstitute of Fluid Mechanics. and Computer Applications in Civil Engineering University of Hannover. Prof. Dr.-Ing. W. Zielke.
Institute of Fluid Mechanics and Computer Applications in Civil Engineering University of Hannover Prof. Dr.-Ing. W. Zielke ROCKFLOW Tutorial E Version 0.1.0 for RockFlow 3.8.68 Software developers: René
More informationEnabling Technologies
Enabling Technologies Mechanical Modelling 1 Key Parameter Mineral System Exploration is reflected in scale-dependent translation A. Gradient in hydraulic potential B. Permeability C. Solubility sensitivity
More informationMathematical modelling on transport of petroleum hydrocarbons in saturated fractured rocks
Sādhanā Vol. 39, Part 5, October 2014, pp. 1119 1139. c Indian Academy of Sciences Mathematical modelling on transport of petroleum hydrocarbons in saturated fractured rocks G SURESH KUMAR Petroleum Engineering
More informationCHAPTER 2. SOIL-WATER POTENTIAL: CONCEPTS AND MEASUREMENT
SSC107 Fall 2000 Chapter 2, Page - 1 - CHAPTER 2. SOIL-WATER POTENTIAL: CONCEPTS AND MEASUREMENT Contents: Transport mechanisms Water properties Definition of soil-water potential Measurement of soil-water
More informationAvailable online at ScienceDirect. Energy Procedia 114 (2017 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 114 (2017 ) 3312 3321 13th International Conference on Greenhouse Gas Control Technologies, GHGT-13, 14-18 November 2016, Lausanne,
More informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationMacro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media
Transp Porous Med (2005) 58:121 145 Springer 2005 DOI 10.1007/s11242-004-5472-6 Macro-Scale Dynamic Effects in Homogeneous and Heterogeneous Porous Media SABINE MANTHEY 1,, S. MAJID HASSANIZADEH 2 and
More informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
More informationA note on benchmarking of numerical models for density dependent flow in porous media
Advances in Water Resources 29 (2006) 1918 1923 www.elsevier.com/locate/advwatres A note on benchmarking of numerical models for density dependent flow in porous media B. Ataie-Ashtiani *, M.M. Aghayi
More informationRelative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure
GRC Transactions, Vol. 36, 2012 Relative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure Noriaki Watanabe, Keisuke
More information*** ***! " " ) * % )!( & ' % # $. 0 1 %./ +, - 7 : %8% 9 ) 7 / ( * 7 : %8% 9 < ;14. " > /' ;-,=. / ١
١ ******!" #$ % & '!( ) % * ") +,-./ % 01. 3 ( 4 56 7/4 ) 8%9 % : 7 ;14 < 8%9 % : *7./ = ;-, >/'." Soil Permeability & Seepage ٢ Soil Permeability- Definition ٣ What is Permeability? Permeability is the
More informationAbstract. 1 Introduction
Correlations for mass transfer coefficients applicable to NAPL pool dissolution in subsurface formations C.V. Chrysikopoulos & T.-J. Kim Department of Civil and Environmental Engineering, University of
More informationDepartment of Ocean Engineering, Indian Institute of Technology-Madras, Chennai , India. *Corresponding author.
J. Earth Syst. Sci. (2018) 127:53 c Indian Academy of Sciences https://doi.org/10.1007/s12040-018-0950-3 Interaction of dissolution, sorption and biodegradation on transport of BTEX in a saturated groundwater
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationTemperature dependent multiphase flow and transport
Temperature dependent multiphase flow and transport J.F. Sykes, A.G. Merry and J. Zhu Department of Civil Engineering, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1 E-mail: sykesj@uwaterloo.ca
More informationII. Dealing with Physical Laws: A Stochastic Flowpath Analysis Example
II. Dealing with Physical Laws: A Stochastic Flowpath Analysis Example In this chapter we will examine a physical problem where the significance of observing the underlying principles will be displayed.
More informationEvaporation-driven transport and precipitation of salt in porous media: A multi-domain approach
Evaporation-driven transport and precipitation of salt in porous media: A multi-domain approach Vishal Jambhekar Karen Schmid, Rainer Helmig Department of Hydromechanics and Hydrosystemmodeling EGU General
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationChapter 2 Overview of Processes Occurring During CO 2 Geological Storage and Their Relevance to Key Questions of Performance
Chapter 2 Overview of Processes Occurring During CO 2 Geological Storage and Their Relevance to Key Questions of Performance Chin-Fu Tsang and Auli Niemi Abstract The objective of this chapter is to provide
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationInfluence of Soil Characteristic and Rainfall Intensity on Matric Suction of Unsaturated Weathered Soil Slope
Journal of the Korean Society of Civil Engineers Kim, Yong Min*, Lee, Kwang Woo**, Kim, Jung Hwan*** Influence of Soil Characteristic and Rainfall Intensity on Matric Suction of Unsaturated Weathered Soil
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationFully implicit higher-order schemes applied to polymer flooding. K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
Fully implicit higher-order schemes applied to polymer flooding K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF) Motivation Challenging and complex to simulate water-based EOR methods, unresolved
More informationAbility of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium Situations
Research Article imedpub Journals http://www.imedpub.com Resources, Recycling and Waste Management Ability of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium
More informationSENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN SALINE AQUIFERS USING THE PROBABILISTIC COLLOCATION APPROACH
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 SENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN
More informationChapter 1 INTRODUCTION
Chapter 1 INTRODUCTION 1-1 The Fluid. 1-2 Dimensions. 1-3 Units. 1-4 Fluid Properties. 1 1-1 The Fluid: It is the substance that deforms continuously when subjected to a shear stress. Matter Solid Fluid
More informationDarcy's Law. Laboratory 2 HWR 531/431
Darcy's Law Laboratory HWR 531/431-1 Introduction In 1856, Henry Darcy, a French hydraulic engineer, published a report in which he described a series of experiments he had performed in an attempt to quantify
More informationA Generalized Numerical Approach for Modeling Multiphase Flow and Transport in Fractured Porous Media
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 1, pp. 85-108 Commun. Comput. Phys. July 2009 A Generalized Numerical Approach for Modeling Multiphase Flow and Transport in Fractured Porous Media Yu-Shu
More informationNumerical fractional flow modelling of inhomogeneous
Numerical fractional flow modelling of inhomogeneous air sparging E.F. Kaasschieter\ G.J. Mulder^ & J.D. van der Werff ten Bosch^ Department of Mathematics and Computing Science, hoven University of Technology,
More informationA Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations
A Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations Investigators Hamdi Tchelepi, Associate Professor of Energy Resources Engineering;
More informationWESTCARB Regional Partnership
WESTCARB Regional Partnership Subsurface Flow Modeling at King Island Christine Doughty, CADoughty@lbl.gov Curtis Oldenburg, CMOldenburg@lbl.gov Staff Scientists Lawrence Berkeley National Laboratory WESTCARB
More information1 Modeling Immiscible Fluid Flow in Porous Media
Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.
More informationRADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA
GEOPHYSICS RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA C. BUCUR 1, M. OLTEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, Pitesti, Romania, crina.bucur@scn.ro 2 Academy of Scientists Bucharest,
More informationAdvanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur. Lecture 6
Advanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur Lecture 6 Good morning and welcome to the next lecture of this video course on Advanced Hydrology.
More informationMulti-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media Jan Tecklenburg a,, Insa Neuweiler a, Jesus Carrera b, Marco Dentz b a Institute of Fluid Mechanics
More informationABSTRACT INTRODUCTION
Transport of contaminants from non-aqueous phase liquid pool dissolution in subsurface formations C.V. Chrysikopoulos Department of Civil Engineering, University of California, ABSTRACT The transient contaminant
More informationAPPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell
APPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell AaronJ.Mango* Mark W. Schmeeckle* David Jon Furbish* Department of Geological Sciences, Florida
More informationSoil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND
Soil Water Atmosphere Plant (SWAP) Model: I. INTRODUCTION AND THEORETICAL BACKGROUND Reinder A.Feddes Jos van Dam Joop Kroes Angel Utset, Main processes Rain fall / irrigation Transpiration Soil evaporation
More informationTable 17 1 Some general field equation terms. Heat Power. Current Source. 0 0 Boundary Current Porous Media Flow. Flow Source
17 Related Analogies 17.1 Basic Concepts The differential equation used in a finite element study in one discipline often appears in a different discipline, but with a different physical meaning for the
More informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
More informationPore-Level Bénard Marangoni Convection in Microgravity
Pore-Level Bénard Marangoni Convection in Microgravity Peyman Mohammadmoradi, and Apostolos Kantzas * Chemical and Petroleum Engineering Department, University of Calgary *Corresponding author: 2500 University
More informationFracture relative permeability revisited
Fracture relative permeability revisited NOROLLAH KASIRI and GHASEM BASHIRI, Iran University of Science and Technology Relative permeability is one of the most uncertain terms in multiphase flow through
More informationHydrocarbon Fluid Wave and it's resonance G.K. Batta*, D.S. Manral, G.V.J. Rao, and S.K. Basha, Oil India limited, India
Hydrocarbon Fluid Wave and it's resonance G.K. Batta*, D.S. Manral, G.V.J. Rao, and S.K. Basha, Oil India limited, India gnangp@gmail.com Keywords Hydrocarbon fluid wave, seismic acoustic impedance, law
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationThermohaline Flow and Reactive Solute Transport in Fractured Rock
Thermohaline Flow and Reactive Solute Transport in Fractured Rock Thomas GRAF Institute of Fluid Mechanics in Civil Engineering Alberta Geological Survey, Edmonton 2008-2009 Basin Hydrogeologist Sasdf
More informationRelative permeability characteristics necessary for simulating DNAPL infiltration, redistribution, and immobilization in saturated porous media
WATER RESOURCES RESEARCH, VOL. 39, NO. 8, 1213, doi:10.1029/2002wr001490, 2003 Relative permeability characteristics necessary for simulating DNAPL infiltration, redistribution, and immobilization in saturated
More informationNumerical Treatment of Two-Phase Flow in Porous Media Including Specific Interfacial Area
Procedia Computer cience Volume 51, 2015, Pages 1249 1258 ICC 2015 International Conference On Computational cience Numerical Treatment of To-Phase Flo in Porous Media Including pecific Interfacial Area
More informationNumerical modelling of coupled TH processes to facilitate analysis of geothermal reservoirs - Groß Schönebeck as an example-
Numerical modelling of coupled TH processes to facilitate analysis of geothermal reservoirs - Groß Schönebeck as an example- N Watanabe 1, G Blöcher 2, M Cacace 2, H Shao 1, O Kolditz 1 1 Helmholtz Centre
More informationSimulating Fluid-Fluid Interfacial Area
Simulating Fluid-Fluid Interfacial Area revealed by a pore-network model V. Joekar-Niasar S. M. Hassanizadeh Utrecht University, The Netherlands July 22, 2009 Outline 1 What s a Porous medium 2 Intro to
More informationMODELING OF GAS MIGRATION THROUGH LOW-PERMEABILITY CLAY USING INFORMATION ON PRESSURE AND DEFORMATION FROM FAST AIR INJECTION TESTS
PROCEEDINGS, TOUGH Symposium 2015 Lawrence Berkeley National Laboratory, Berkeley, California, September 28-30, 2015 MODELING OF GAS MIGRATION THROUGH LOW-PERMEABILITY CLAY USING INFORMATION ON PRESSURE
More information(Refer Slide Time: 02:10)
Soil Mechanics Prof. B.V.S. Viswanathan Department of Civil Engineering Indian Institute of Technology, Bombay Lecture 24 Flow of water through soils-v Welcome to lecture five of flow of water through
More informationWater in Soil Sections in Craig
Water in Soil Sections 2.1-2.6 in Craig Outlines Introduction Darcy s Law Volume of water flowing per unit time Measuring K in laboratory Seepage Theory Flow Net Introduction All soils are permeable materials,
More informationPiping Systems and Flow Analysis (Chapter 3)
Piping Systems and Flow Analysis (Chapter 3) 2 Learning Outcomes (Chapter 3) Losses in Piping Systems Major losses Minor losses Pipe Networks Pipes in series Pipes in parallel Manifolds and Distribution
More informationA BENCHMARK CALCULATION OF 3D HORIZONTAL WELL SIMULATIONS
INTERNATINAL JURNAL F NUMERICAL ANALYSIS AND MDELING Volume 1, Number 2, Pages 189 201 c 2004 Institute for Scientific Computing and Information A BENCHMARK CALCULATIN F 3D HRIZNTAL WELL SIMULATINS ZHANGIN
More informationPHYSICS OF FLUID SPREADING ON ROUGH SURFACES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Supp, Pages 85 92 c 2008 Institute for Scientific Computing and Information PHYSICS OF FLUID SPREADING ON ROUGH SURFACES K. M. HAY AND
More informationTable of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials
Table of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials 1.3 Scope of This Book 1.4 Historical Development of Geotechnical
More information1 FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE
FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE 1-1 1 FLUID-MECHANICAL INTERACTION SINGLE FLUID PHASE 1.1 Introduction FLAC models the flow of fluid (e.g., groundwater) through a permeable solid, such
More informationCO2 Storage Trapping Mechanisms Quantification
CO2 Storage Trapping Mechanisms Quantification CO2 storage The capture and storage of in deep geological formations is one of the proposed solutions to reduce emissions to the atmosphere. CO2 storage is
More informationAgain we will consider the following one dimensional slab of porous material:
page 1 of 7 REVIEW OF BASIC STEPS IN DERIVATION OF FLOW EQUATIONS Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations,
More information