190 likes | 327 Views
This overview presents a comprehensive examination of Darcy's law, the transient mass balance equation, and steady-state conditions in groundwater flow. Emphasizing anisotropic and heterogeneous media, we analyze the characteristics of hydraulic conductivity (K) in various dimensional frameworks. The text details the implementation of the Gauss-Seidel iteration method to solve the Toth problem in 2D groundwater scenarios. Special attention is given to specified head boundary conditions, showcasing examples with initial guesses and iteration solutions to illustrate concepts effectively.
E N D
q = - Kgrad h Darcy’s law grad is a vector
Transient mass balance equation: q is a vector
Steady State Mass Balance Equation with W = 0 gradq = 0 q = 0 (del notation) The dot product of grad and q is div q = 0
Anisotropic medium: K at a point varies with direction Kx Ky Kz or Kx= Ky Kz Heterogeneous medium: K varies in space
Characteristics of K in two dimensions A Kx A B B Kz Kx Kz Kx =Kz = K Homogeneous, isotropic Homogeneous, anisotropic B A A B Heterogeneous, anisotropic Heterogeneous, isotropic
Characteristics of K in two dimensions A Kx A B B Kz Kx =Kz = K Kx Kz Homogeneous, isotropic Homogeneous, anisotropic B A A B Heterogeneous, anisotropic Heterogeneous, isotropic
homogeneous, isotropic porous medium: Kx = Ky = Kz = K steady state conditions: W = 0 Laplace equation h = 0 2
Toth Problem (2D) Groundwater divide Groundwater divide Impermeable Rock 2D, steady state
i i, j-1 The 5 point star in a regular grid j i, j i-1, j i+1, j y i, j+1 x = y x
Iteration indices Solution by iteration Gauss-Seidel Iteration
Gauss-Seidel Iteration solution m+3 Iteration planes m+2 m+1 initial guesses m
1 2 3 4 i 1 2 3 4 5 6 10 9 8 7 10 6 2D example with specified head boundary conditions 10 5 10 5 (Assume an initial guess of h =10 for all interior nodes) 6 10 9 8 10 7 j
1 2 3 4 i 1 2 3 4 5 6 10 9 8 7 10 6 9.75 Example with specified head boundary conditions 10 5 10 5 6 10 9 8 10 7 j
1 2 3 4 i 1 2 3 4 5 6 10 9 8 7 9.75 10 6 Example with specified head boundary conditions 8.44 10 5 10 5 6 10 9 8 10 7 j
1 2 3 4 i 1 2 3 4 5 6 10 9 8 7 8.44 9.75 10 6 Example with specified head boundary conditions 9.94 10 5 10 5 6 10 9 8 10 7 j
1 2 3 4 i 1 2 3 4 5 6 10 9 8 7 10 6 Example with specified head boundary conditions 9.94 10 5 10 5 6 10 9 8 10 7 j
i j i, j Gauss-Seidel Iteration
Example of spreadsheet formula Toth Problem