1 / 25

General governing equation for steady-state , heterogeneous , anisotropic conditions

2D Laplace Eqn. --Homogeneous and isotropic aquifer without a sink/source term. --2D flow in a profile; (Unconfined aquifer with a water table boundary condition; recharge occurs as a result of the boundary condition.). General governing equation

ludlow
Télécharger la présentation

General governing equation for steady-state , heterogeneous , anisotropic conditions

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. 2D Laplace Eqn. --Homogeneous and isotropic aquifer without a sink/source term. --2D flow in a profile; (Unconfined aquifer with a water table boundary condition; recharge occurs as a result of the boundary condition.) General governing equation for steady-state, heterogeneous, anisotropic conditions

  2. Mathematical Model of the Toth Problem Unconfined aquifer h = c x + zo

  3. Aquifer b z x x z b z x Toth problem b = 1 m

  4. Aquifer b 2D horizontal flow in an aquifer with constant thickness, b. y x

  5. unconfined aquifer b is not constant b confined aquifer Figure from Hornberger et al. 1998

  6. Aquifer b 2D horizontal flow in an aquifer with constant thickness, b. with recharge y x

  7. 2D horizontal flow 2D horizontal flow; homogeneous and isotropic aquifer with constant aquifer thickness, b, so that T=Kb. with a source/sink term Poisson Equation

  8. South Fork Map of Long Island, N.Y. Charles Edward Jacob (1914-1970) Consultant to the Town of Southampton, NY December 1968

  9. h datum C.E. Jacob’s Conceptual Model of the South Fork of Long Island water table R groundwater divide ocean ocean b x = - L x = 0 x = L We can simulate this system assuming horizontal flow in a “confined” aquifer if we assume that T= Kb.

  10. at x =0 h(L) = 0 1D approximation used by C.E. Jacob R ocean ocean x = - L x = 0 x = L Boundary Conditions Governing Eqn.

  11. C.E. Jacob’s Model Governing Eqn. at x =0 Boundary conditions h(L) = 0 Analytical solution for 1D “confined” version of the problem h(x) = R (L2 – x2) / 2T Forward solution R = (2 T) h(x) / ( L2 – x2) Inverse solution for R

  12. Inverse solution for R R = (2 T) h(x) / ( L2 – x2) Solve for R with h(x) = h(0) = 20 ft. Inverse solution for T Rearrange eqn to solve for T, given value for R and h(0) = 20 ft. Observation well on the groundwater divide L

  13. Island Recharge Problem L y 2L ocean ocean well T= 10,000 ft2/day L = 12,000 ft x ocean

  14. Targets used in Model Calibration • Head measured in an observation well is known as a head target. • The simulated head at the node representing the observation well is compared with the measured head. • During model calibration, parameter values (e.g., R and T) are adjusted until the simulated head matches the observed value. • Model calibration solves the inverse problem.

  15. Island Recharge Problem L = 12,000 ft y 2L ocean Solve the forward problem: Given R= 0.00305 ft/d T= 10,000 ft2/day Solve for h at each nodal point ocean well x ocean

  16. Write the finite difference approximation: Gauss-Seidel Iteration Formula for 2D Poisson Equation with x = y = a

  17. Island Recharge Problem 4 X 7 Grid

  18. Water Balance IN = Out IN = R x AREA Out = outflow to the ocean

  19. Top 4 rows Head at a node is the average head in the area surrounding the node. Red dots represent specified head cells, which are treated as inactive nodes. Black dots are active nodes. (Note that the nodes along the groundwater divides are active nodes.)

  20. L 2L Top 4 rows IN =R x Area = R (L-x/2) (2L - y/2) Also: IN = R (2.5)(5.5)(a2)

  21. x x/2 x Top 4 rows

  22. Qy = (Th) /2 Qy x/2 Qx Top 4 rows OUT = Qy +  Qx Qy = K (x b) (h/y) Note: x = y  Qy = T h Qx = K (y b)(h/x) or Qx = T h y

  23. Bottom 4 rows y/2 Well Qx = (T h)/2

  24. Island Recharge Problem 4 X 7 Grid

  25. Water Budget Error

More Related