Simulation of seawater intrusion in coastal aquifers: Some typical responses
|
|
- Arron Small
- 5 years ago
- Views:
Transcription
1 Sādhanā, Vol. 26, Part 4, August 2001, pp Printed in India Simulation of seawater intrusion in coastal aquifers: Some typical responses 1. Introduction AMLAN DAS 1 and BITHIN DATTA 2 1 Department of Civil Engineering, Regional Engineering College, Durgapur , India. 2 Department of Civil Engineering, Indian Institute of Technology, Kanpur , India. adas wrpm@yahoo.com; bithin@iitk.ac.in MS received 1 March 1999; revised 9 April 2001 Abstract. Seawater intrusion in coastal aquifers is generally three dimensional (3-D) in nature. In the literature, there is a general lack of reported results on 3-D simulations. This paper presents some typical example simulations of 3-D seawater intrusion process for a specified hypothetical study area. The simulation results presented in this paper are based on the density-dependent miscible flow and transport modelling approach for simulation of seawater intrusion in coastal aquifers. A nonlinear optimization-based simulation methodology was used in this study. Various steady state simulations are performed for a specified study area. Response evaluations consider the effects of vertical recharge on seawater intrusion, effects of boundary conditions, and effects of spatially varying pumping from the aquifer. The 3-D simulations demonstrate the viability of using a planned strategy of spatially varying withdrawals from the aquifer to manage seawater intrusion. It is demonstrated that series of pumps near the ocean-face boundary induce a hydraulic head distribution that can be effectively used for controlling seawater intrusion. Keywords. Seawater intrusion; coastal aquifers; density-dependent flow and transport; 3-D aquifer responses. The seawater intrusion mechanism in coastal aquifers generally causes the occurrence of a freshwater zone over a seawater zone. In between these two zones, there occurs a transition zone. Depending upon the relative thickness of this transition zone, two approaches are used to model the phenomenon of seawater intrusion in coastal aquifers. When the transition zone is narrow in thickness and areal extent, the seawater intrusion phenomenon is modelled as a two-phase fluid flow separated by a sharp interface. Mercer et al (1980) and Essaid (1990) present some of the recently developed simulation models based on A list of symbols is given at the end of the paper 317
2 318 Amlan Das and Bithin Datta this sharp interface approach. On the other hand, when the transition zone stretches areally as well as vertically to a considerable extent, the seawater intrusion phenomenon is modelled using the density-dependent miscible flow and transport approach. Huyakorn et al (1997), Galeati et al (1992), Putti & Paniconi (1995), and Das & Datta (1995, 2000) present some of the recently developed simulation models that are based on this second approach. These simulation models establish their applicability mostly using a standard two-dimensional example case of seawater intrusion in coastal aquifers, viz. the Henry problem (Henry 1964). However, seawater intrusion in coastal aquifers is generally three-dimensional (3-D) in nature. In the literature, there is a general lack of reported results on 3-D simulations. Till date, Huyakorn et al (1987) reported a 3-D steady state simulation example. This paper, therefore, presents some typical example simulations of the 3-D seawater intrusion process for a specified study area. The simulation results presented in this paper are based on the density-dependent miscible flow and transport modelling approach for simulation of the seawater intrusion in coastal aquifers. Seawater intrusion can be caused by over-pumping in sensitive portions of coastal aquifers. It often results in the degradation of water quality within the aquifer and may ultimately require costly remedial measures. Development of a well-planned pumping strategy can arrest further degradation of the aquifer water quality, remediate an already contaminated aquifer, and contain the contamination within a certain region of the aquifer. Mathematical modelling of seawater intrusion plays a key role in the development of such optimal operating strategies. Ideally, these mathematical models should be capable of simulating 3-D seawater intrusion cases. However, the mathematical models reported in the literature mostly consider two-dimensional example cases. In order to develop any simulation and also optimization model that necessarily seeks to optimally exploit a coastal aquifer, it is necessary to investigate the coastal aquifer responses to the plausible stress scenarios. Theoretically, there can be an infinite number of plausible scenarios of aquifer stresses. However, a few typical stress scenarios might prove to be sufficient to (i) analyse the nature of aquifer responses, (ii) demonstrate the adequacy and viability of incorporating a mathematical model for simulation of the seawater intrusion phenomenon in coastal aquifers within a management model, and (iii) motivate the future formulation of optimization models. This paper presents the responses of coastal aquifers to some typical scenarios of aquifer stresses. An optimization-based simulation model was used in this particular study. The simulation model essentially incorporates the finite difference form of the governing partial differential equations of the seawater intrusion in coastal aquifers as equality type constraints. The simulation model was applied to an aquifer system and the aquifer responses were analysed. These aquifer responses are presented in this paper. 2. Simulation model The simulation model adopted in this study, describing the seawater intrusion process is based on the following assumptions (Huyakorn et al 1987). (1) The physical process is described by density-dependent miscible transport approach. (2) Reference hydraulic head, referred to as freshwater head, and solute concentration are the independent variables. (3) In the three-dimensional Cartesian coordinate system, the x 2 axis points in the upward vertical direction.
3 Simulation of seawater intrusion in coastal aquifers 319 (4) The dynamic viscosity of the fluid does not change with respect to salinity and remains equal to that of freshwater. (5) The density of the mixed fluid is a linear function of the solute concentration. (6) The actual hydraulic conductivity is a linear function of the solute concentration of mixed fluid (seawater and freshwater). (7) Solute concentration is zero for freshwater (minimum density) and one for seawater condition (maximum density). (8) In order to define the Darcy velocity vector, the actual hydraulic conductivity term is replaced by freshwater hydraulic conductivity. (9) Intrinsic permeability is independent of solute concentration. The three-dimensional flow equation may be written as (Huyakorn et al 1987): ( ) h h K x ij + ηce i x j = S S j t + φη c t ρ ( ) qr q ρ p, (1) o where, K ij is the actual hydraulic conductivity tensor (LT 1 ), h is the reference hydraulic head (L), x j (j = 1, 2, 3) are Cartesian coordinates (L), η is the density coupling coefficient, c is the solute concentration (0 c 1) (ML 3 ), e j is the jth component of the gravitational unit vector (e 1 = 0,e 2 = 1,e 3 = 0), S S is the specific storage (L 1 ), t is time (T), φ is porosity, q r and q p are the volumetric flow rate of sources and sinks respectively, per unit volume of the porous medium (L 3 T 1 L 3 ), and, ρ and ρ o are the density of mixed fluid (freshwater and seawater) and reference (freshwater) density (ML 3 ). Note that Huyakorn et al (1987) used a single term q to define the volumetric flow rate of sources (or sinks) per unit volume of the porous medium (L 3 T 1 L 3 ). The q of Huyakorn et al (1987) is replaced by q r q p in (1). This ensures the nonnegativity of source and sink terms. Reference head is defined as: h = ( p/ρ o g ) + Y, (2) where, p is the fluid pressure (ML 1 T 2 ), g is the gravitational acceleration (LT 2 ), Y is the elevation above datum (L). Density coupling coefficient is expressed as (Huyakorn et al 1987): η = ε/c s, (3) where, ε is the density difference ratio defined as: ε = ( ) ρ s ρ o /ρ o, (4) and, c s is the solute concentration (ML 3 ) corresponding to the maximum density ρ s (ML 3 ). Actual hydraulic conductivity is defined as: K ij = k ij ρg/µ (5) where, k ij is the intrinsic permeability tensor (L 2 ), µ is the dynamic viscosity of fluid (µ = µ o )(ML 1 T 1 ), µ o is the viscosity of freshwater (ML 1 T 1 ).
4 320 Amlan Das and Bithin Datta Density of the actual fluid is expressed as: ρ = ρ o (1 + εc/c s ) (6) The advective dispersive equation to describe the salt transport is written as (Huyakorn et al 1987): ( ) c c D x ij V i x i = φ c j x i t + q rc, (7) where, D ij = φ D ij with D ij being the dispersion tensor (L 2 T 1 ) defined by Bear (1979, p. 234), V i is the Darcy velocity vector (LT 1 ). The Darcy velocity vector is expressed as: V i = Kij o h + ηce x j, (8) j where, Kij o is the hydraulic conductivity at the reference (freshwater) condition (LT 1 ). Note that Huyakorn et al (1987) used a term q(c c ) where c is the solute concentration in injected (or withdrawn) fluid (ML 3 ). The following simple algebraic manipulation is accomplished to replace the q(c c ) term with q r c: q(c c ) = (q r q p )(c c ) (a) = q r (c cr ) q p (c c p ), (b) where, cr is the solute concentration in the injected fluid (ML 3 ), and cp is the solute concentration in the withdrawn fluid (ML 3 ). As the concentration in the withdrawn fluid is the same as the existing concentration at that location, cp = c and if only freshwater is assumed to be recharged, cr = 0. Therefore, q(c c ) = q r c, (c) is the reduced form. The initial and boundary conditions for the flow equation may be expressed as: and h(x i, 0) = h o (x i ), (9) h(x i,t) = h (for some parts of the boundary), (10) V i n i = V n (for the remaining part of the boundary), (11) where h o is the initial head (L), h is the prescribed head (L), n i is the outward unit normal vector, V n is the prescribed fluid flux (LT 1 ) which is positive for inflow and negative for outflow. The initial and boundary conditions for the transport equation may be written as: c(x i, 0) = c o, (12) c(x i,t)=c (for some part of the boundary), (13) c D ij n x i = qc D (for some other part of the boundary), (14) j
5 Simulation of seawater intrusion in coastal aquifers 321 and D ij c x j n i V i n i c = q T c (on the remaining part of the boundary), (15) where, c o is the initial concentration (ML 3 ), c is the prescribed concentration in specified parts of the boundary (ML 3 ), qc D is the dispersive mass flux of solute (ML 2 T 1 ), and qc T is the total mass flux of solute (ML 2 T 1 ) in which qc D and qc T are positive for inward mass flux and negative for outward mass flux. 3. Simulation methodology Numerical simulation of seawater intrusion with a considerable transition zone requires simultaneous solution of the discretized flow and transport equations. A nonlinear optimization method for solving the finite difference approximated governing equations for simulation of seawater intrusion is used here. The main advantage of this method is that sequential iterative solutions of the flow and transport equations are not necessary. Also, simultaneous solution of the flow and transport equations using this optimization approach uniquely suits the transformation of this simulation approach to a management model. Finite difference approximation of the flow and transport equations, i.e. (1) and (7) constitute the simulation model of density-dependent seawater intrusion in coastal aquifers. The discretized equations are treated as binding equality constraints of an optimization model in the present simulation approach. The finite difference approximated simulation equations are given in Appendix A. The discretized nonlinear equations may be expressed in a functional form as: f I (u J ) = 0; I =1,..., N t ; J = 1,..., 2xN t. (16) Here, u J is a vector consisting of values of h, c, q r and q p s in the three-dimensional space and time domain of the problem; N t is the total number of equality type simulation constraints. The set of equations (16) represents all the finite difference approximated equations required to simulate the intrusion process. These equations can be solved using an appropriate numerical method. Alternatively, an optimization-based method is used here. In this solution method the system of nonlinear equations (16) are modified as f I (u J ) = ε I ; I =1,..., N t, (17) where, ε I represents error terms which are measures of the deviations of the left hand side of (16) from zero. The actual solution of the seawater intrusion problem has to identify a solution vector for which all ε I = 0 or are close to zero, i.e. within acceptable tolerance. Optimization based simulation approach finds the solution vectors by solving the following optimization model. This model seeks to minimize the sum of the absolute values of ε I for all I s. Formally, the objective function is expressed as: N t Minimize: ε I, (18) I=1 Subject to: f I (u J ) = ε I ; I =1,..., N t. (19) The objective function (18) together with the constraint set (19) represents a nonlinear programming problem. The finite difference approximated density-dependent miscible
6 322 Amlan Das and Bithin Datta flow and transport equations are the constraints in this mathematical programming problem. The density coupling term causes the nonlinearity in the optimization model. The h, c, q r, and q p s defined in the discretized space time domain of the physical problem are the decision variables of the optimization model. In this study, the aforesaid optimization model is solved by the optimization system solver MINOS. The MINOS algorithm uses the projected augmented Lagrangian method of nonlinear programming (Murtagh & Saunders 1993). The detailed solution method is given in Das (1995), Das & Datta (1995, 1998a,b, 2000). 4. Steady state simulations A number of simulations are performed for 3-D steady state flow and transport conditions. These simulation results are helpful in studying the simulated or predicted response of the aquifer to 3-D seawater intrusion. The aquifer specified is a hypothetical unconfined aquifer with uniformly distributed vertical recharge. The aquifer and the boundary conditions are described schematically in figure 1. The geometric dimensions of the aquifer are 1000m (length) 1000m (width) 50 m (depth). The right hand face represents the ocean face. Opposite to the ocean face is the inland face, which is a source of freshwater. The front, back and bottom faces of the aquifer are assumed to be impermeable. The top face is the idealized phreatic surface. The idealization of the phreatic surface is identical to that by Huyakorn et al (1987). In order to define the transport boundary conditions, the ocean face is divided into two regions (Huyakorn et al 1987). The top narrow outflow region allows outflow of water from the aquifer to the ocean. In this part, the dispersive flux is zero at the boundary. The remaining part, the influx region of the ocean face, has a constant concentration boundary (with c = 1). The seawater intrudes into the coastal aquifer through the influx region of the boundary. Throughout this study, identical boundary conditions are prescribed for the ocean face, in which the top 20% portion of the aquifer allows outflow of water into the ocean. For the phreatic surface of the aquifer the boundary condition is described by zero concentration of the injected fluid (Huyakorn et al 1987). Assuming a freshwater vertical recharge, the advective influx component becomes zero (Huyakorn et al 1987). Boundary conditions for the front, back, and bottom faces are specified as zero dispersive flux conditions. Boundary conditions for the inland face are specified as zero concentration condition, by assuming a freshwater source. The flow boundary condition for the ocean face is specified as a hydrostatic boundary condition assuming that the water surface elevation does not vary along the length of the coast line. Boundary conditions in the inland face are specified as a constant reference hydraulic head along the vertical direction but linearly varying along the length of the inland face. This ensures a perfect 3-D flow condition. Because of the assumption of impermeable boundaries, zero flux boundary conditions are specified for the front, back, and bottom faces. For the phreatic surface, the boundary conditions are specified in a manner similar to that by Galeati et al (1992). Various parameter values specified for the illustrative simulations are: Kxx o = Ko zz = 20 m/day; Kyy o = 0.2 m/day; V r = 0.02 m/year; h f = 0.95 m; h = 0.2 m; y = 10 m; z = x = 100 m; φ = 0.26; ε = 0.025; α L = 25 m; α T = 5m;d o =0.66 m 2 /day = m 2 /s. Here, Kxx o,ko yy and Ko zz are the freshwater hydraulic conductivities (L/T) in x,y and z directions respectively; V r is the uniformly distributed vertical recharge intensity (LT 1 ), h f is the reference hydraulic head at location (0, 0) in plan (L); h
7 Simulation of seawater intrusion in coastal aquifers 323 Figure 1. Schematic description of the hypothetical unconfined aquifer and the boundary conditions: (a) plan view, (b) elevation. is the difference in reference hydraulic heads between locations (0, 0) and (0, 1000 m) (L); x, y and z are the grid sizes in x,y and z directions respectively (L); φ is the porosity; ε is the density difference ratio; α L and α T are the longitudinal and transverse dispersivities (L); and d o is the molecular diffusion (L 2 T 1 ). In these simulations, pumping from a finite difference cell is considered synonymous to pumping from a single well. Also, the grid sizes chosen are not large enough to allow installation of more than one well in a single cell. The screen location for a specified well denotes the vertical layer of pumping in the finite difference grid system. Note that the simulation constraints for steady state simulations use the fully implicit scheme, and a one-time period simulation. Further, the S S ( h/ t) and φ ( c/ t) terms in flow and transport equations are also set equal to zero. Initially, the effect of grid size on the solution results was studied. This study consists of obtaining simulation results for various specifications of the grid sizes and comparing
8 324 Amlan Das and Bithin Datta Figure 2. Simulation results for the base example: (a) isoheads (m), (b) isochlors. the solution results for the steady state 3-D seawater intrusion in the unstressed aquifer. The solution results show that the obtained solutions were acceptable for grid sizes of x = z = 100 m, and y = 10 m, and were specified for further evaluations. Figures 2a and b show the simulated isoheads and isochlors for this case (henceforth, termed as base example solutions). The isoheads represent contours of equal reference hydraulic heads, and isochlors represent contours of equal salt concentrations. 4.1 Effect of vertical recharge on seawater intrusion The distribution of salt concentration in the aquifer is affected by the amount of freshwater recharge entering the saturated zone vertically. The following examples demonstrate the
9 Simulation of seawater intrusion in coastal aquifers 325 effects of varying amounts of distributed vertical recharge, V r, on the 3-D seawater intrusion. In addition to the base example, three other example cases were studied (viz. A1, A2, and A3 respectively). The vertical recharge values as specified for these three additional example cases were 0.1, 0.2 and 0.3 m/year respectively. The aquifer remains unstressed as no pumping is specified. The other parameter values in these examples remain the same as that of the base example. Simulation results in terms of the isoheads and isochlors for these three examples are shown in figures 3a and b, 4a, 4b, 5a, and 5b respectively. It is interesting to observe the changes in isohead patterns with the increase in V r. Steady state seawater intrusion is observed to reduce with increase in V r, as evident from the relative decrease in salt concentration in the aquifer. It is seen from these figures that due to increase in distributed vertical recharge, the reference hydraulic head increases. This increase in the reference head is responsible for reducing Figure 3. Simulation results for example A1. (a) Isoheads (m), (b) isochlors.
10 326 Amlan Das and Bithin Datta Figure 4. Simulation results for example A2. (a) Isoheads (m), (b) isochlors. the amount of seawater intrusion in the aquifer. From these examples, it is observed that the amount of seawater intrusion is inversely related to the amount of freshwater vertical recharge, as intuitively expected. 4.2 Effect of constant head boundary condition on the seawater intrusion A constant head boundary condition signifies the potential for unlimited recharge through the constant head boundary. The intrusion of seawater in the coastal aquifers will be certainly affected by the presence of constant head boundary conditions. The effect of the magnitude of the head in the constant head boundary is investigated in the following examples. In addition to the base example case, three more example cases were studied
11 Simulation of seawater intrusion in coastal aquifers 327 Figure 5. Simulation results for example A3. (a) Isoheads (m), (b) isochlors. (viz. B1, B2, and B3 respectively). Three values of h f = 1.0, 0.975, and m respectively, were specified for these particular example cases, to represent three different constant head magnitudes. Figures 6a, 6b, 7a, 7b, 8a, and 8b show the relative positions of isoheads and isochlors in these three examples respectively. In these examples, a nonzero value of h is maintained to represent 3-D flow conditions. It may be observed that isohead patterns gradually change with increase in head at the constant head boundary. Seawater intrusion increases with a decrease in h f. In these examples, it is observed that the amount of seawater intrusion as noted from the simulated isochlors is inversely proportional to the specified constant heads in the boundary. Volker & Rushton (1982) show similar decrease in amount of seawater intrusion with increase in the amount of constant freshwater flux through the boundary, for 2-D simulation. Since an increase in
12 328 Amlan Das and Bithin Datta Figure 6. Simulation results for example B1. (a) Isoheads (m), (b) isochlors. constant head at the boundary is also likely to increase the flux through it, similarity in resultant effects is expected. 4.3 Effect of pumping on seawater intrusion One important aim of seawater intrusion management in coastal aquifers should be to prescribe optimum level of exploitation of the aquifer. In order to assess the viability of exploiting a coastal aquifer for various purposes, it is required to study the effects of spatially varying pumping on the seawater intrusion. Infinite number of feasible but nonoptimal exploitation scenarios are possible in aquifer management problems. The following examples simulate few simple cases of spatially varying pumping and the resulting effects on seawater intrusion.
13 Simulation of seawater intrusion in coastal aquifers 329 Figure 7. Simulation results for example B2. (a) Isoheads (m), (b) isochlors. 4.3a Variation of pumping at a given location: The following four examples show the effect of varying the steady state pumping on seawater intrusion. In all these examples the location of the well as well as that of the screen are prespecified. Figure 9 shows the locations of the well and the screen. The coordinate of the well location in plan is (550 m, 550 m). In the vertical direction the central 10 m portion of the length of the well is screened. The parameter values and the boundary conditions identical to that in the base example are used in these examples. Pumping rates of 100 m 3 /day, 200 m 3 /day, 300 m 3 /day and 400 m 3 /day are specified in examples C1, C2, C3, and C4 respectively. The simulated isoheads and isochlors are shown in figures 10a, 10b, 11a, 11b, 12a, 12b, 13a, and 13b respectively. It may be observed from these figures that as the magnitude
14 330 Amlan Das and Bithin Datta Figure 8. Simulation results for example B3. (a) Isoheads (m), (b) isochlors. of pumping from the well increases, the steady state positions of 0.01, 0.05, 0.1, 0.3 and 0.5 isochlors move towards the inland face of the aquifer. This signifies an increased amount of seawater intrusion into the aquifer. As compared to the base example, these simulated isochlors represent increased intrusion of seawater into the aquifer, particularly in the vicinity of the pumping well. The isoheads also show a decrease in head near the well as pumping rate is increased. These examples demonstrate that the steady state pumping causes a reduction of reference hydraulic heads near the well. This reduction of the hydraulic head causes the seawater to intrude more towards the inland face of the aquifer. It is further observed that as pumping increases, the pumped water gets more contaminated with saline. These simulation examples establish the fact that stressing an aquifer, by increasing pumping
15 Simulation of seawater intrusion in coastal aquifers 331 Figure 9. Locations of the well and screen for examples C1 to C4. from it, will definitely affect the salt concentration distribution in the aquifer. These simple illustrative examples only show the adverse effects of pumping on salinity in the aquifer, as expected. However, a planned spatial and temporal variation in pumping may have a resultant beneficial effect in controlling aquifer contamination at required locations and times (Das, 1995). 4.3b Variation of screen position at a given well: The location of the screen in a pumping well determines the vertical layer through which water is entering the well. The effect of changing the screen location in a pumping well on steady state seawater intrusion is investigated in the following examples (D1 and D2). The location of the well in plan, boundary conditions, other parameters and quantity of pumping for these two examples are specified to be the same as that in example C4. The 10 m long screen is located near the top of the aquifer in example D1 and near the bottom of the aquifer in example D2. The isoheads and isochlors, as obtained from the simulation results are shown in figures 14a and b, 15a and 15b, for examples D1 and D2 respectively. The isohead patterns in examples C4, D1 and D2 are different from each other. It may be
16 332 Amlan Das and Bithin Datta Figure 10. Simulation results for example C1. (a) Isoheads (m), (b) isochlors. observed that the values of salinity of pumped water (as evident from the isochlors) in example D1 and D2 are less than that in example C4. As expected, pumping from different vertical locations induces different hydraulic head distributions, which results in different salt concentration distributions in the aquifer. 4.3c Variation of well location: Spatial locations of pumping wells significantly influence the hydraulic head distribution in an aquifer. It is therefore expected that the spatial distribution of pumping wells will significantly influence the salt concentration distribution in the aquifer also. When localized demand for water is to be satisfied and, simultaneously, a certain upper limit on salinity of the pumped water is to be maintained, the location of the well becomes an important aspect in coastal aquifer management. The effect of
17 Simulation of seawater intrusion in coastal aquifers 333 Figure 11. Simulation results for example C2. (a) Isoheads (m), (b) isochlors. changing the well location on the steady state seawater intrusion is investigated in the following examples. The location of the screen, boundary conditions, other parameters and quantity of pumping in these four examples are specified to be the same as that in example C4. The coordinates of well locations in the plan view are (550 m, 250 m), (550 m, 750 m), (250 m, 550 m) and (750 m, 550 m) for examples E1, E2, E3, and E4 respectively. Computed results in terms of isoheads and isochlors are shown in figures 16a and b, 17a and b, 18a and b, and 19a and b for examples E1, E2, E3 and E4 respectively. It is observed that contamination of pumped water (which is same as the salt concentration at the location of pumping) is highly dependent on the location of the well. In all these examples the contamination of the aquifer is higher in regions in the immediate
18 334 Amlan Das and Bithin Datta Figure 12. Simulation results for example C3. (a) Isoheads (m), (b) isochlors. vicinity of the pumping well compared to that in nearby locations. In examples E1 and E2, the contamination of the pumped water is higher than that in example C4. In example E4, due to the constant pumping at location (750 m, 550 m), the steady state 0.01 isochlor occupies a new position compared to the unstressed situation in the base example, indicating remediation to a certain extent. Examples C4, E1, E2, E3 and E4 show the effect of identical pumping with spatial variation in location. As spatial variations in location induce different hydraulic head configurations, these effects on concentration distribution are expected. 4.3d Series of wells close to ocean face: Todd (1980) prescribed heavy pumping in a row along the sea shore to form an extraction barrier as one of the useful methods of
19 Simulation of seawater intrusion in coastal aquifers 335 Figure 13. Simulation results for example C4. (a) Isoheads (m), (b) isochlors. controlling seawater intrusion. In the following examples, the effect of pumping from a series of wells located close to the ocean face is investigated. The boundary conditions and other parameters for these two examples are the same as in the base example. The locations of the wells are shown in figure 20. The wells are assumed to be screened throughout the depth of the aquifer. There are 10 wells located along the sea face in both the examples. Each well pumps at uniform rates of 200 m 3 /day and 100 m 3 /day in examples F1 and F2 respectively. Figures 21a and 21b, and 22a and 22b show the isoheads, and the isochlors as simulated from these two examples respectively. In example F2, the pumping rate is half of that in example F1, and the locations of the wells are closer to the ocean face.
20 336 Amlan Das and Bithin Datta Figure 14. Simulation results for example D1. (a) Isoheads (m), (b) isochlors. Simulation results show that the intrusion of seawater as evident from the salt concentration distribution is less in example F2, compared to the base example (unstressed) and example F1. This shows the efficiency of an extraction barrier closer to the sea face in controlling seawater intrusion in the aquifer. The extraction barrier is essentially a meeting point of hydraulic gradients in two different directions developed due to pumping. A series of pumps closer to the sea face is therefore more efficient for controlling the entry of seawater. A comparison of simulations for the base example, and F1 and F2 demonstrate these expected results. The illustrative simulations of the 3-D steady state cases of seawater intrusion presented in this study were limited to a small hypothetical study area. For these examples, the
21 Simulation of seawater intrusion in coastal aquifers 337 Figure 15. Simulation results for example D2. (a) Isoheads (m), (b) isochlors. total number of binding nonlinear constraints were All the simulations reported here were performed by using a DEC-ALPHA 3000/600 AXP Workstation (175 MHz, 40 MFLOPS) with DEC OSF/1 V2.0 operating system. CPU time required for a single steady state simulation, on the average, is about 15 minutes. 5. Conclusions Responses for 3-D cases of seawater intrusion are studied and presented. Various steadystate simulations are performed for a specified study area. These simulated responses
22 338 Amlan Das and Bithin Datta Figure 16. Simulation results for example E1. (a) Isoheads (m), (b) isochlors. will undoubtedly enhance the future development of simulation and also management models for optimal exploitation of coastal aquifers. The response evaluations consider the effects of vertical recharge on seawater intrusion, effects of boundary conditions, and effects of spatially varying pumping from the aquifer. The solution results are intuitively acceptable. The 3-D simulations also demonstrate the viability of using a planned strategy of spatially varying withdrawals from the aquifer to manage seawater intrusion. It is also shown that the modifications in the salt concentration distributions, due to an imposed pattern of pumping, can be beneficially used for controlling seawater intrusion. In particular, it is noted that a series of pumps near the ocean face boundary induce a hydraulic head distribution that can be effectively used for controlling seawater intrusion.
23 Simulation of seawater intrusion in coastal aquifers 339 Figure 17. Simulation results for example E2. (a) Isoheads (m), (b) isochlors. List of symbols c solute concentration (0 c 1)(ML 3 ); c solute concentration of the injected fluid (0 c 1)(ML 3 ); c s solute concentration (ML 3 ) corresponding to the maximum density ρ s (ML 3 ); c o initial concentration (ML 3 ); c prescribed concentration in specified parts of the boundary (ML 3 ); D ij φ D ij ;
24 340 Amlan Das and Bithin Datta Figure 18. Simulation results for example E3. (a) Isoheads (m), (b) isochlors. D ij dispersion tensor L 2 T 1 ; d 0 molecular diffusion coefficient; e j jth component of gravitational unit vector (e 1 = 0,e 2 =1,e 3 =0); g gravitational acceleration (LT 2 ); h reference hydraulic head (L); h o initial head (L); h prescribed head (L); K ij actual hydraulic conductivity tensor (LT 1 ); Kxx o,ko yy,ko zz freshwater hydraulic conductivities corresponding to x,y and z directions respectively;
25 Simulation of seawater intrusion in coastal aquifers 341 Figure 19. Simulation results for example E4. (a) Isoheads (m), (b) isochlors. Kij o hydraulic conductivity at the reference (freshwater) condition (LT 1 ); k ij intrinsic permeability tensor (L 2 ); outward unit normal vector; n i p fluid pressure (ML 1 T 2 ); q r,q p volumetric flow rate of sources and sinks respectively, per unit volume of the porous medium (L 3 T 1 L 3 ); qc D dispersive mass flux of solute (ML 2 T 1 ); q c total mass flux of solute (ML 2 T 1 ); S S specific storage (L 1 ); t time (T );
26 342 Amlan Das and Bithin Datta Figure 20. Location of wells for examples F1 and F2. U fresh water influx through the inland face; V i Darcy velocity vector (LT 1 ); V n prescribed fluid flux (LT 1 ); x j (j = 1, 2, 3) Cartesian coordinates (L); Y elevation above datum (L); α L,α T longitudinal and lateral dispersivities respectively; η density coupling coefficient; φ porosity; ρ,ρ o density of mixed fluid (freshwater and seawater) and reference (freshwater) density (ML 3 ) respectively; ε density difference ratio; µ dynamic viscosity of fluid (µ = µ o )(ML 1 T 1 ); µ o viscosity of the freshwater (ML 1 T 1 ); x, y, z finite difference grid sizes along x,y, and z directions respectively.
27 Simulation of seawater intrusion in coastal aquifers 343 Figure 21. Simulation results for example F1. (a) Isoheads (m), (b) isochlors. Appendix A. Finite difference approximated simulation equations. Finite difference approximation of the flow and transport equations are based on the point centred central difference scheme (Bear 1979). In order to approximate the flow equation, the left hand side of (1) is expressed as: ( ) h K x ij + ηce i x j j = K ij x i ( ) h + ηce x j j ( ) h + K ij + ηce x i x j. (A1) j
28 344 Amlan Das and Bithin Datta Figure 22. Simulation results for example F2. (a) Isoheads (m), (b) isochlors. As an illustration, the finite difference form of K xx / x is expressed as ( ) Kxx )/ = θ (K n K n 2 x xxi,j+1,k xxi,j 1,k x ( +(1 θ) Kxx n 1 i,j+1,k Kxx n 1 i,j 1,k )/2 x, (A2) where, θ is the time weighting factor (for an explicit scheme θ = 0, for an implicit scheme θ = 1, and for a Crank Nicolson scheme θ = 0.5), x is the grid size along the x direction, Kxx n is the hydraulic conductivity along the x direction for the mixed fluid at node () at time level n(lt 1 ), and is expressed as: Kxx n = Kxx o ( ) 1 + ηc n, (A3)
29 Simulation of seawater intrusion in coastal aquifers 345 where Kxx o is the hydraulic conductivity along the x direction for the reference fluid (freshwater) at node ()(LT 1 ) which does not vary with respect to time, c n is the solute concentration at node () at time level n(ml 3 ). In order to approximate the transport equation, the terms containing the dispersion tensor in (7) are approximated as ( ) c 2 c D x ij D i x ij. (A4) j x i x j / Here the ( D ij xi )( c / ) x j term is assumed to be very small in comparison to (D ij ( 2 c / )) x i x j. This assumption is, however, strictly valid only when the longitudinal dispersivity, α L (L), and transverse dispersivity, α T (L) are zero or close to zero. In the finite difference form of the transport equation, the dispersion tensor is computed at every node for all time periods. Thus, the dispersion coefficients are functions of hydraulic heads and solute concentration, and also varies spatially and temporally. The hydrodynamic dispersion i.e., dispersion tensor is expressed as (Huyakorn et al 1987; Bear 1979), D xx = φ D xx, ( ) / D xx = α T V 2 y + V 2 z + α L V 2 x V +d o, ( ) / D yy = α T V 2 x + V 2 z + α L V 2 y V +d o, ( ) / D zz = α T V 2 x + V 2 y + α L V 2 z V +d o, (A5) (A6) (A7) (A8) D xy = ( α L α T ) V x V y / V = D yx, D xz = ( α L α T ) V x V z / V = D zx, D yz = ( α L α T ) V y V z / V = D yz, (A9) (A10) (A11) where V x, V y and V z are the absolute values of Darcy velocities V x,v y and V z along x,y and z directions respectively (LT 1 ), d o is the molecular diffusion (L 2 T 1 ), and V is the absolute value of resultant Darcy velocity V(LT 1 ) which is expressed as: V = ( Vx 2 + V y 2 + V z 2 ) 1/2. (A12) The flow and transport equations, i.e. (1) and (7) are approximated by using the point centred central difference scheme (Bear 1979). The flow equation is now expressed as: A 1 h n i,j+1,k + A 2 hn i,j 1,k + A 3 hn i+1,j,k + A 4 hn i 1,j,k + A 5 hn +1 +A 6 h n 1 + A 7 hn + A 8 cn i+1,j,k + A 9 cn i 1,j,k + A 10 cn +A 11 qr n + A 12 qp n + B 1 h n 1 i,j+1,k + B 2 hn 1 i,j 1,k + B 3 hn 1 i+1,j,k +B 4 h n 1 i 1,j,k + B 5 hn B 6 hn B 7 hn 1 + B 8 cn 1 i+1,j,k +B 9 c n 1 i 1,j,k + B 10 cn 1 + B 11 qn 1 r + B 12 qp n 1 = 0, (A13)
30 346 Amlan Das and Bithin Datta where, A 1 = θ ) ( ) θ (K n K n + (1 θ) K n 1 4 x 2 xxi,j+1,k xxi,j 1,k xx i,j+1,k Kxx n 1 i,j 1,k A 2 = θ θ 4 x 2 A 3 = + θ θk n x 2 xx + (1 θ) Kxx n 1, (A14) ) ( ) (K n K n + (1 θ) K n 1 xxi,j+1,k xxi,j 1,k xx i,j+1,k Kxx n 1 i,j 1,k + θ θk n x 2 xx + (1 θ) Kxx n 1, (A15) θ ) ( ) θ (K n K n + (1 θ) K n 1 4 y 2 yyi+1,j,k yyi 1,j,k yy i+1,j,k Kyy n 1 i 1,j,k + θ θk n y 2 yy + (1 θ) Kyy n 1, (A16) A 4 = θ θ 4 y 2 A 5 = ) ( ) (K n K n + (1 θ) K n 1 yyi+1,j,k yyi 1,j,k yy i+1,j,k Kyy n 1 i 1,j,k + θ θk n y 2 yy + (1 θ) Kyy n 1, (A17) θ ) ( ) θ (K n K n + (1 θ) K n 1 4 z 2 zz+1 zz 1 zz +1 Kzz n 1 1 A 6 = θ θ 4 z 2 + θ θk n z 2 zz + (1 θ) Kzz n 1, (A18) ) ( ) (K n K n + (1 θ) K n 1 zz+1 zz 1 zz +1 Kzz n θ θk n z 2 zz + (1 θ) Kzz n 1, (A19) A 7 = 2θ θk n x 2 xx + (1 θ) Kxx n 1 2θ θk n y 2 yy + (1 θ) Kyy n 1 2θ θk n z 2 zz + (1 θ) Kzz n 1 S S t, (A20) A 8 = θη θkyy n 2 y + (1 θ) Kyy n 1, (A21) A 9 = θη θkyy n 2 y + (1 θ) Kyy n 1, (A22) A 10 = θη ( ) ( ) θ Kyy n 2 y i+1,j,k Kyy n 1 i 1,j,k + (1 θ) Kyy n 1 i+1,j,k Kyy n 1 i 1,j,k φη t, (A23)
31 Simulation of seawater intrusion in coastal aquifers 347 { } A 11 = θ 1 + η θc n + (1 θ) cn 1, (A24) A 12 = θ B 1 = { 1+η θc n } + (1 θ) cn 1, (A25) (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 x 2 xxi,j+1,k xxi,j 1,k xx i,j+1,k Kxx n 1 i,j 1,k (1 θ) + θk n x 2 xx + (1 θ) Kxx n 1, (A26) B 2 = (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 x 2 xxi,j+1,k xxi,j 1,k xx i,j+1,k Kxx n 1 i,j 1,k (1 θ) + θk n x 2 xx + (1 θ) Kxx n 1, (A27) B 3 = (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 y 2 yyi+1,j,k yyi 1,j,k yy i+1,j,k Kyy n 1 i 1,j,k (1 θ) + θk n y 2 yy + (1 θ) Kyy n 1, (A28) B 4 = (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 y 2 yyi+1,j,k yyi 1,j,k yy i+1,j,k Kyy n 1 i 1,j,k (1 θ) + θk n y 2 yy + (1 θ) Kyy n 1, (A29) B 5 = (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 z 2 zz+1 zz 1 zz +1 Kzz n 1 1 (1 θ) + θk n z 2 zz + (1 θ) Kzz n 1, (A30) B 6 = (1 θ) ) ( ) θ (K n K n + (1 θ) K n 1 4 z 2 zz+1 zz 1 zz +1 Kzz n 1 1 (1 θ) + θk n z 2 zz + (1 θ) Kzz n 1, (A31) B 7 = 2(1 θ) θk n x 2 xx + (1 θ) Kxx n 1 2 (1 θ) θk n y 2 yy + (1 θ) Kyy n 1 2 (1 θ) θk n z 2 zz + (1 θ) Kzz n 1 + S S t, (A32) B 8 = (1 θ)η 2 y B 9 =+ (1 θ)η 2 y θkyy n + (1 θ) Kyy n 1, (A33) θk n yy + (1 θ) K n 1 yy, (A34)
32 348 Amlan Das and Bithin Datta B 10 = (1 θ)η 2 y + (1 θ) B 11 = (1 θ) 1 + η ) θ (K n K n yyi+1,j,k yyi 1,j,k ( ) Kyy n 1 i+1,j,k Kyy n 1 i 1,j,k + φη t, { θc n { B 12 = (1 θ) 1+η θc n + (1 θ) cn 1 + (1 θ) cn 1 (A35) }, (A36) }, (A37) K n xx = K o xx ( 1 + ηc n ), (A38) K n yy = K o yy ( 1 + ηc n ), (A39) K n zz = K o zz ( 1 + ηc n ). (A40) The transport equation can be similarly expressed as: where E 1 ci,j+1,k n + E 2 cn i,j 1,k + E 3 cn i+1,j,k + E 4 cn i 1,j,k + E 5 cn +1 +E 6 c 1 n + E 7 cn i+1,j+1,k + E 8 cn i 1,j 1,k + E 9 cn i+1,j 1,k + E 10 cn i 1,j+1,k +E 11 ci,j+1,k+1 n + E 12 cn i,j 1,k 1 + E 13 cn i,j+1,k 1 + E 14 cn i,j 1,k+1 + E 15 cn i+1,j,k+1 +E 16 ci 1,j,k 1 n + E 17 cn i+1,j,k 1 + E 18 cn i 1,j,k+1 + E 19 cn i,j,j,k F 1 c n 1 i,j+1,k + F 2 cn 1 i,j 1,k + F 3 cn 1 i+1,j,k + F 4 cn 1 i 1,j,k + F 5 cn F 6 c n F 7 cn 1 i+1,j+1,k + F 8 cn 1 i 1,j 1,k + F 9 cn 1 i+1,j 1,k + F 10 cn 1 i 1,j+1,k +F 11 c n 1 i,j+1,k+1 + F 12 cn 1 i,j 1,k 1 + F 13 cn 1 i,j+1,k 1 + F 14 cn 1 i,j 1,k+1 + F 15 cn 1 i+1,j,k+1 +F 16 c n 1 i 1,j,k 1 + F 17 cn 1 i+1,j,k 1 + F 18 cn 1 i 1,j,k+1 + F 19 cn 1 i,j,j,k = 0, (A41) E 1 = E 2 = E 3 = E 4 = x 2 2 x x x Dxx n + (1 ) D n 1 Vx n + (1 ) Vx n 1, (A42) xx Dxx n + (1 ) D n 1 Vx n + (1 ) Vx n 1, (A43) xx D n y 2 yy + (1 ) D n 1 2 y yy Vy n + (1 ) Vy n 1, (A44) D n y 2 yy + (1 ) D n y yy V n y + (1 ) V n 1 y, (A45)
33 Simulation of seawater intrusion in coastal aquifers 349 E 5 = E 6 = E 7 = D n z 2 zz + (1 ) D n 1 2 z zz Vz n + (1 ) Vz n 1, (A46) D n z 2 zz + (1 ) D n z 2 x y zz Vz n + (1 ) Vz n 1, (A47) Dxy n + (1 ) Dxy n 1, (A48) E 8 = E 7, E 9 = E 7, (A49) (A50) E 10 = E 7, E 11 = (A51) Dyz n 2 y z + (1 ) Dyz n 1, (A52) E 12 = E 11, E 13 = E 11, (A53) (A54) E 14 = E 11, (A55) E 15 = Dxz n 2 z x + (1 ) Dxz n 1, (A56) E 16 = E 15, (A57) E 17 = E 15, E 18 = E 15, (A58) (A59) E 19 = 2 D n x 2 xx + (1 ) Dxx n 1 2 D n y 2 yy + (1 ) Dyy n 1 F 1 = F 2 = F 3 = 2 D n z 2 zz + (1 ) Dzz n 1 φ t qr n + (1 ) qr n 1, (A60) (1 ) D n x 2 xx + (1 ) Dxx n 1 (1 ) Vx n 2 x + (1 ) Vx n 1, (A61) (1 ) D n x 2 xx + (1 ) Dxx n 1 (1 ) + Vx n 2 x + (1 ) Vx n 1, (A62) (1 ) D n y 2 yy + (1 ) Dyy n 1 (1 ) 2 y V n y + (1 ) V n 1 y, (A63)
34 350 Amlan Das and Bithin Datta F 4 = F 5 = F 6 = (1 ) D n y 2 yy + (1 ) Dyy n 1 (1 ) + Vy n 2 y + (1 ) Vy n 1, (A64) (1 ) D n z 2 zz + (1 ) Dzz n 1 (1 ) Vz n 2 z + (1 ) Vz n 1, (A65) (1 ) D n z 2 zz + (1 ) Dzz n 1 (1 ) + Vz n 2 z + (1 ) Vz n 1, (A66) (1 ) F 7 = Dxy n 2 x y + (1 ) Dxy n 1, (A67) F 8 = F 7, (A68) F 9 = F 7, (A69) F 10 = F 7, (A70) (1 ) F 11 = Dyz n 2 y z + (1 ) Dyz n 1, (A71) F 12 = F 11, (A72) F 13 = F 11, (A73) F 14 = F 11, (1 ) F 15 = 2 z x F 16 = F 15, (A74) Dxz n + (1 ) Dxz n 1, (A75) (A76) F 17 = F 15, F 18 = F 15, (A77) (A78) F 19 = 2(1 ) D n x 2 xx + (1 ) Dxx n 1 2 (1 ) D n y 2 yy + (1 ) Dyy n 1 2 (1 ) D n z 2 zz + (1 ) Dzz n 1 + φ t (1 ) qr n + (1 ) qr n 1, (A79) / (h Vx n = (K 2 x) o n xx i,j+1,k h n i,j 1,k), (A80) {( (h Vy n = Kyy o n i+1,j,k h n i 1,j,k) / ) } 2 y ηc n, (A81) / (h Vz n = (K 2 z) o n zz +1 h n 1), (A82)
35 Simulation of seawater intrusion in coastal aquifers 351 V n = V n2 x + V n2 y + V n2 z 1/2, D n xx = D n yy = D n zz = D n xy = D n yz = D n xz = (A83) {α L V n2 x + α T (V n2 y + V n2 z )} /V n + d o φ, (A84) {α L V n2 y + α T (V n2 z + V n2 x )} /V n + d o φ, (A85) {α L V n2 z + α T (V n2 x + V n2 y )} /V n + d o φ, (A86) (αl ) α T V n V n x y/ V n φ, (A87) (αl ) α T V n V n y z/ V n φ, (A88) (αl ) α T V n V n x z/ V n φ, (A89) where V n is the absolute value of V n, V n x is the absolute value of V n x, V n y is the absolute value of Vy n, V n z is the absolute value of Vz n. For a steady state solution (performed over one time step), the storage coefficients for flow and transport, S S and φ are set to zero and a value of θ = = 1 is used. This amounts to using the fully implicit scheme. Generally, θ = is used throughout this study. References Bear J 1979 Hydraulics of groundwater (New York: McGraw Hill) Das A 1995 Development of nonlinear optimization based single and multiple objective management models for transient 3-D density dependent seawater intrusion in coastal aquifers, Ph D thesis, Civil Engineering Department, Indian Institute of Technology, Kanpur Das A, Datta B 1995 Simulation of density dependent 2-d seawater intrusion in coastal aquifers using nonlinear optimization algorithm. Proc. Am. Water Res. Assoc. Annu. Summer Symp. on Water resources and environmental hazards: Emphasis on hydrologic and cultural insight in the Pacific rim, Honolulu, Hawaii, pp Das A, Datta B 1998a Development of multiobjective management models for coastal aquifers. J. Water Resour. Planning Manage., Am. Soc. Civil Eng. 125: Das A, Datta B 1998b Development of management models for sustainable use of coastal aquifers. J. Irrig. Drainage Eng., Am. Soc. Civil Eng. 125: Das A, Datta B 2000 Optimization based solution of density dependent seawater intrusion in coastal aquifers. J. Hydrol. Eng., Am. Soc. Civil Eng. 5: Essaid H I 1990 A multilayered sharp interface model of coupled freshwater and saltwater in coastal systems: model development and application. Water Resour. Res. 27: Galeati G, Gambolati G, Neuman S P 1992 Coupled and partially coupled Eulerian Lagrangian model of freshwater-seawater mixing. Water Resour. Res. 28: Henry H R 1964 Effects of dispersion of salt encroachment in coastal aquifers, sea water in coastal aquifers. US Geol. Surv. Water Supply Pap C, C70 C84
36 352 Amlan Das and Bithin Datta Huyakorn P S, Anderson P F, Mercer J W, White W O Jr 1987 Saltwater intrusion in aquifers: development and testing of a three dimensional finite element model. Water Resour. Res. 23: Mercer J W, Larson S P, Faust C R 1980 Simulation of saltwater interface motion. Ground Water 18: Murtagh B A, Saunders M A 1993 Minos 5.4 user s guide. Technical report, SOL 83 20R, Systems Optimization Laboratory, Department of Operations Research, Stanford University, Stanford, CA Putti M, Paniconi C 1995 Picard and Newton linearization for the coupled model of saltwater intrusion in aquifers. Adv. Water Resour. 18: Todd D K 1980 Groundwater hydrology (New York: John Wiley) Volker R, Rushton K 1982 An assessment of the importance of some parameters for seawater intrusion in aquifers and comparison of dispersive and sharp-interface modeling approaches. J. Hydrol. 56:
Numerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
Evaluation of the hydraulic gradient at an island for low-level nuclear waste disposal
A New Focus on Groundwater Seawater Interactions (Proceedings of Symposium HS1001 at IUGG2007, Perugia, July 2007). IAHS Publ. 312, 2007. 237 Evaluation of the hydraulic gradient at an island for low-level
The lens of freshwater in a tropical island: the two interface case
ANZIAM J. 51 (EMAC2009) pp.c1 C15, 2010 C1 The lens of freshwater in a tropical island: the two interface case S. A. Chen 1 G. C. Hocking 2 (Received 29 October 2009; revised 10 February 2010) Abstract
MODELING THE SALTWATER INTRUSION PHENOMENON IN COASTAL AQUIFERS - A CASE STUDY IN THE INDUSTRIAL ZONE OF HERAKLEIO IN CRETE
Global NEST Journal, Vol. 7, No 2, pp 197-203, 2005 Copyright 2005 Global NEST Printed in Greece. All rights reserved MODELING THE SALTWATER INTRUSION PHENOMENON IN COASTAL AQUIFERS - A CASE STUDY IN THE
Improving the worthiness of the Henry problem as a benchmark for density-dependent groundwater flow models
WATER RESOURCES RESEARCH, VOL. 40, W01504, doi:10.1029/2003wr002199, 2004 Improving the worthiness of the Henry problem as a benchmark for density-dependent groundwater flow models Matthew J. Simpson 1
Two-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
A 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
Pollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
Darcy 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 ]
Geophysical Surveys for Groundwater Modelling of Coastal Golf Courses
1 Geophysical Surveys for Groundwater Modelling of Coastal Golf Courses C. RICHARD BATES and RUTH ROBINSON Sedimentary Systems Research Group, University of St. Andrews, St. Andrews, Scotland Abstract
Computer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design
CVEN 302-501 Computer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design Date distributed : 12/2/2015 Date due : 12/9/2015 at
dynamics 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
J. Environ. Res. Develop. Journal of Environmental Research And Development Vol. 8 No. 1, July-September 2013
SENSITIVITY ANALYSIS ON INPUT PARAMETERS OF 1-D GROUNDWATER FLOW GOVERNING EQUATION : SOLVED BY FINITE-DIFFERENCE AND THOMAS ALGORITHM IN MICROSOFT EXCEL Goh E.G.* 1 and Noborio K. 2 1. Department of Engineering
GG655/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
EFFECTIVE DISPERSION IN SEAWATER INTRUSION THROUGH HETEROGENEOUS AQUIFERS
18 SWIM. Cartagena 2004, Spain. (Ed. Araguás, Custodio and Manzano). IGME EFFECTIVE DISPERSION IN SEAWATER INTRUSION THROUGH HETEROGENEOUS AQUIFERS E. ABARCA, 1, J. CARRERA 1, R. HELD 2, X. SANCHEZ-VILA
THREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
EVALUATION 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
Darcy s law in 3-D. K * xx K * yy K * zz
PART 7 Equations of flow Darcy s law in 3-D Specific discarge (vector) is calculated by multiplying te ydraulic conductivity (second-order tensor) by te ydraulic gradient (vector). We obtain a general
Analytical solutions for water flow and solute transport in the unsaturated zone
Models for Assessing and Monitoring Groundwater Quality (Procsedines of a Boulder Symposium July 1995). IAHS Publ. no. 227, 1995. 125 Analytical solutions for water flow and solute transport in the unsaturated
16 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
LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort
- 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [
In 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
ESTIMATION OF TRANSVERSE DISPERSTVITY IN THE MIXING ZONE OF FRESH-SALT GROUNDWATER
ModelCARE 90: Calibration and Reliability in Groundwater Modelling (Proceedings of the conference held in The Hague, September 1990). IAHS Publ. no. 195, 1990. ESTIMATION OF TRANSVERSE DISPERSTVITY IN
Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
Australian Journal of Basic and Applied Sciences, 5(6): 536-543, 0 ISSN 99-878 Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
PLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
The Saltwater Intrusion Problem: Concepts, Challenges, Modeling, Experiments and Case Studies
The Saltwater Intrusion Problem: Concepts, Challenges, Modeling, Experiments and Case Studies Manfred Koch Department of Geohydraulics and Engineering Hydrology, University of Kassel, Germany Email: kochm@uni-kassel.de
Ping Cheng Department of Mechanical Engineering University of Hawaii Honolulu, Hawaii 96822
NUMERICAL AND ANALYTICAL STUDIES ON HEAT AND MASS TRANSFER IN VOLCANIC ISLAND GEOTHERMAL RESERVOIRS Ping Cheng Department of Mechanical Engineering University of Hawaii Honolulu, Hawaii 96822 The Hawaii
Frozen Ground Containment Barrier
Frozen Ground Containment Barrier GEO-SLOPE International Ltd. www.geo-slope.com 1200, 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: +1 403 269 2002 Fax: +1 888 463 2239 Introduction Frozen soil barriers
Getting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
WQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
RT3D Rate-Limited Sorption Reaction
GMS TUTORIALS RT3D Rate-Limited Sorption Reaction This tutorial illustrates the steps involved in using GMS and RT3D to model sorption reactions under mass-transfer limited conditions. The flow model used
Temperature 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
Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams
Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty
Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
Semi-Analytical 3D solution for assessing radial collector well pumping. impacts on groundwater-surface water interaction
1 2 Semi-Analytical 3D solution for assessing radial collector well pumping impacts on groundwater-surface water interaction 3 Ali A. Ameli 1 and James R. Craig 2 4 5 1 Global Institute for Water Security,
Chapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
Deflated CG Method for Modelling Groundwater Flow Near Faults. L.A. Ros
Deflated CG Method for Modelling Groundwater Flow Near Faults Interim Report L.A. Ros Delft, Utrecht January 2008 Delft University of Technology Deltares Supervisors: Prof.dr.ir. C. Vuik (TU Delft) Dr.
Hydraulic 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
CH.3. COMPATIBILITY EQUATIONS. Multimedia Course on Continuum Mechanics
CH.3. COMPATIBILITY EQUATIONS Multimedia Course on Continuum Mechanics Overview Introduction Lecture 1 Compatibility Conditions Lecture Compatibility Equations of a Potential Vector Field Lecture 3 Compatibility
Chapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
The similarity solution of the longitudinal dispersion phenomenon of miscible fluids in porous media
Int. J. Adv. Appl. Math. and Mech. 1(3) (2014) 77-82 ISSN: 2347-2529 Available online at www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics The similarity solution of
ABSTRACT GOVERNING EQUATIONS
A three dimensional finite element model for fluid flow and transport in confined or unconfined aquifer F. Jacob*, J.M. Crolef, P. Lesaint*, J. Mania* "Laboratoire de calcul scientifique, ^Laboratoire
Finite Element Solution to Groundwater Transport of Solutes Undergoing Decay and Non-Linear Sorption
Hydrology Days 2010 1 Finite Element Solution to Groundwater Transport of Solutes Undergoing Decay and Non-Linear Sorption Domenico Baú Department of Civil & Environmental Engineering Colorado State University,
Advanced 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.
Chapter 14: Groundwater. Fig 14.5b
Chapter 14: Groundwater Fig 14.5b OBJECTIVES Recognize that groundwater is a vital source of accessible freshwater. Describe how groundwater forms below the water table. Explain the origin of aquifers,
Contents Page TABLE OF FIGURES... VII GLOSSARY...XI INTRODUCTION... 1
ii Abstract Despite the long history of work on saltwater wedge intrusion in coastal aquifers, very little attention has been paid to the effect of basement irregularities on the wedge. This work combines
Earth dam steady state seepage analysis
Engineering manual No. 32 Updated 3/2018 Earth dam steady state seepage analysis Program: FEM Water Flow File: Demo_manual_32.gmk Introduction This example illustrates an application of the GEO5 FEM module
Simulation of salt migrations in density dependent groundwater flow
Simulation of salt migrations in density dependent groundwater flow E.S. van Baaren Master s Thesis Applied Mathematics June 27 Thesis Committee: Dr.ir. F.J. Vermolen (Delft University of Technology) Prof.dr.ir.
1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
Continuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets
1.7, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets Equation of Continuity Our equations of hydrogeology are a combination of o Conservation of mass o Some empirical
Todd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
Building ground level
TMA4195 MATHEMATICAL MODELLING PROJECT 212: AQUIFER THERMAL ENERGY STORAGE 1. Introduction In the project we will study a so-called Aquifer Thermal Energy Storage (ATES) system with the aim of climitizing
The Use of COMSOL for Integrated Hydrological Modeling
The Use of for Integrated Hydrological Modeling Ting Fong May Chui *, David L. Freyerg Department of Civil and Environmental Engineering, Stanford University *Corresponding author: 380 Panama Mall, Terman
A THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 40, 1623 1637 (1997) A THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION KANGLE HUANG, JIR[ I
Safety assessment for disposal of hazardous waste in abandoned underground mines
Safety assessment for disposal of hazardous waste in abandoned underground mines A. Peratta & V. Popov Wessex Institute of Technology, Southampton, UK Abstract Disposal of hazardous chemical waste in abandoned
Code Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
Chapter 8 - pg. 1 CHAPTER 8 ESTUARIES. To paraphrase Pritchard, a pioneer in studies of estuarine circulation,
Chapter 8 - pg 1 CHAPTER 8 ESTUARIES Estuaries are semi-closed basins in which a rather complex interaction between river inputs, tidal currents and wind leads to the turbulent mixing of salt from the
Measure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling
Measure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling Steven Mattis and Lindley Graham RMSWUQ 2015 July 17, 2015 Collaborators 2/45 The University of
Computational modelling of reactive transport in hydrogeological systems
Water Resources Management III 239 Computational modelling of reactive transport in hydrogeological systems N. J. Kiani, M. K. Patel & C.-H. Lai School of Computing and Mathematical Sciences, University
V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
OCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
Reservoir Oscillations with Through Flow
American Journal of Environmental Sciences 3 (): 37-42, 27 ISSN 553-345X 27 Science Publications Reservoir Oscillations with Through Flow A. A. Khan 28 Lowry Hall, epartment of Civil Engineering, Clemson
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
H4: Steady Flow through a Channel Network
August 9, Chapter 7 H4: Steady Flow through a Channel Network Contents 7. Problem Specification............................. 7-7. Background................................... 7-3 7.3 Contra Costa Water
THE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL
Tenth International Water Technology Conference, IWTC10 2006, Alexandria, Egypt 281 THE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL Karima Attia 1 and Gamal El Saied 2 1
Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions
Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions David R. Wilkins Copyright c David R. Wilkins 2000 2010 Contents 4 Vectors and Quaternions 47 4.1 Vectors...............................
1906 Dilip kumar Jaiswal et al./ Elixir Pollution 31 (2011) Available online at (Elixir International Journal)
196 Dilip kumar Jaiswal et al./ Eliir Pollution 31 (211) 196-191 ARTICLE INF O Article history: Received: 21 December 21; Received in revised form: 16 January 211; Accepted: 1 February 211; Keywords Advection,
1.061 / 1.61 Transport Processes in the Environment
MIT OpenCourseWare http://ocw.mit.edu 1.061 / 1.61 Transport Processes in the Environment Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 6.
Modular 3-D Transport model. for accommodating add-on reaction packages. Downloads & a wealth of information at: and.
MT3DMS Modular 3-D Transport model MS denotes the Multi-Species i structure t for accommodating add-on reaction packages Downloads & a wealth of information at: http://hydro.geo.ua.edu/mt3d/ edu/mt3d/
ψ 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.........................
The Power of Spreadsheet Models. Mary P. Anderson 1, E. Scott Bair 2 ABSTRACT
MODFLOW 00 and Other Modeling Odysseys Proceedings, 00, International Ground Water Modeling Center, Colorado School of Mines, Golden, CO, p. 85-8 The Power of Spreadsheet Models Mary P. Anderson, E. Scott
Dynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
The Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
v. 8.0 GMS 8.0 Tutorial RT3D Double Monod Model Prerequisite Tutorials None Time minutes Required Components Grid MODFLOW RT3D
v. 8.0 GMS 8.0 Tutorial Objectives Use GMS and RT3D to model the reaction between an electron donor and an electron acceptor, mediated by an actively growing microbial population that exists in both soil
NUMERICAL AND EXPERIMENTAL ANALYSIS OF SEEPAGE BENEATH A MODEL OF A GRAVITY DAM
Engineering Review Vol. 33, Issue 2, 75-84, 2013. 75 NUMERICAL AND EXPERIMENTAL ANALYSIS OF SEEPAGE BENEATH A MODEL OF A GRAVITY DAM T. Jelenkovi V. Travaš * Chair of Hydraulic Engineering, Faculty of
Updating 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
Thermohaline 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
A hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
Contaminant Modeling
Contaminant Modeling with CTRAN/W An Engineering Methodology February 2012 Edition GEO-SLOPE International Ltd. Copyright 2004-2012 by GEO-SLOPE International, Ltd. All rights reserved. No part of this
6. Continuous Release - Point Source
1 6. Continuous Release - Point Source A scalar released continuously into a moving fluid and from a discrete point will form a plume that grows in the lateral dimension through diffusion and extends downstream
Chapter 4: Fluid Kinematics
Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to
Chapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
Contents. Foreword... (vii) Preface...(ix) About the Author...(xi) List of Tables...(xxi) List of Figures... (xxvii)
(xiii) Contents Foreword... (vii) Preface...(ix) About the Author...(xi) List of Tables...(xxi) List of Figures... (xxvii) CHAPTER 1: INTRODUCTION TO WATER RESOURCES 1.1 Global Water Resources... 3 1.2
(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
Determination of Locally Varying Directions through Mass Moment of Inertia Tensor
Determination of Locally Varying Directions through Mass Moment of Inertia Tensor R. M. Hassanpour and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering
STEADY SOLUTE DISPERSION IN COMPOSITE POROUS MEDIUM BETWEEN TWO PARALLEL PLATES
Journal of Porous Media, 6 : 087 05 03 STEADY SOLUTE DISPERSION IN COMPOSITE POROUS MEDIUM BETWEEN TWO PARALLEL PLATES J. Prathap Kumar, J. C. Umavathi, & Ali J. Chamkha, Department of Mathematics, Gulbarga
Chapter 3 Permeability
3.2 Darcy s Law In 1856, Darcy investigated the flow of water through sand filters for water purification. His experimental apparatus is shown in Figure 3.11. By empirical observation Figure 3.11 Schematic
On the advective-diffusive transport in porous media in the presence of time-dependent velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 3, L355, doi:.29/24gl9646, 24 On the advective-diffusive transport in porous media in the presence of time-dependent velocities A. P. S. Selvadurai Department of Civil
Three-dimensional Modelling of Reactive Solutes Transport in Porous Media
151 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 41, 214 Guest Editors: Simonetta Palmas, Michele Mascia, Annalisa Vacca Copyright 214, AIDIC Servizi S.r.l., ISBN 978-88-9568-32-7; ISSN 2283-9216
1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009)
EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) Dr. Mohamad Hekarl Uzir-chhekarl@eng.usm.my School of Chemical Engineering Engineering Campus, Universiti
Seawater circulation in sediments driven by interactions between seabed topography and fluid density
WATER RESOURCES RESEARCH, VOL. 49, 1386 1399, doi:10.1002/wrcr.20121, 2013 Seawater circulation in sediments driven by interactions between seabed topography and fluid density L. F. Konikow, 1 M. Akhavan,
Chapter 2 Linear Optimization
Chapter 2 Linear Optimization Abstract The purpose of this chapter is to cover the basic concept of linear optimization analysis and its applications in water resources engineering. The graphical and simplex
Development of a Combined Quantity and Quality Model for Optimal Unsteady Groundwater Management
Utah State University DigitalCommons@USU Civil and Environmental Engineering Faculty Publications Civil and Environmental Engineering 10-1988 Development of a Combined Quantity and Quality Model for Optimal
3.3 Classification Diagrams Estuarine Zone Coastal Lagoons References Physical Properties and Experiments in
Contents 1 Introduction to Estuary Studies... 1 1.1 Why to Study Estuaries?.... 1 1.2 Origin and Geological Age... 4 1.3 Definition and Terminology... 7 1.4 Policy and Actions to Estuary Preservation....
Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.
FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.Bhosale 1 Research scholar, 2 Head of department & Asst professor Department
RATE 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