Pumping SolutionGeology 724Problem Set #2The Island Recharge Problem-Poisson EquationProblem set #2 consists of variations of the following problems from Wang and Anderson: 3.1, 3.2,3.4, and 3.5. The changes are described below. We will consider four conceptual models: 1Daquifer with constant thickness (i.e., 1D confined aquifer), 2D aquifer with constant thickness (i.e.,2D confined aquifer), 1D unconfined, and 2D unconfined.3.1. 1-D Confined Aquifer: Refer to Fig. 3.2 (p. 43) and follow the directions for problem 3.1 onp. 65 of Wang and Anderson. Then, use the analytical solution you derived to solve the inverseproblem by calculating a value of R for a system with: T = 10,000 ft2/day; l = 12,000 ft; h(0) = 20 ft(e.g., at the center of the island). Also calculate a value for h(0) if R = 0.00305 ft/day.3.2 & 3.5. 2-D Confined Aquifer: Construct a spreadsheet model to solve for the distribution ofheads using Gauss/Seidel iteration for the upper right hand quadrant of the island pictured in Figure3.3. Assume a mesh-centered grid. Print out the array of heads to the fourth decimal place. Includea water balance in your program as described in Problem 3.5 in Wang and Anderson. Use thefollowing values for the aquifer parameters: T = 10,000 ft2/day, l = 12,000 ft and R =0.00305 ft/day.Construct two worksheets in the spreadsheet as follows: 4 x 7 grid cells (x = y = 4000 ft) with error tolerance of 1E-4 13 x 25 grid (x = y = 1000 ft) with error tolerance of 1E-4 Compare your solution with the analytical solution on p. 4 of this handout, which was calculatedfrom the equation given in Problem 3.2 of Wang and Anderson. 3.4. Unconfined Aquifer: Follow the directions in Wang and Anderson. Solve the 1D problemanalytically and the two-dimensional problem using Gauss/Seidel iteration for the upper righthand quadrant of the island with R = 0.00305 ft/day. Be sure to include a water balance in yourspreadsheet. Produce a solution using a 13 x 25 node grid (nodal spacing = 1000 ft). Use anerror tolerance of 1E-4 and set K = 100 ft/day.Add a pumping well to your spreadsheet model of the 2D unconfined aquifer (x = y = 1000 ft). Assume that the well is located in the node at the center of the island (i.e., at the intersection of the groundwater divides). Assume that the well pumps 10% of the water that enters the model via recharge.Modify the finite difference equation in the cell in your spreadsheet in which the pumping well is located. Also modify the outflow term in your water balance to include the volume of water pumped out of the model by the well.In your results, report the value of head in the cell that contains the pumping well, i.e. the cell in the center of the island. Also use equation 5.7 on p. 149 of Anderson and Woessner to estimate the head in the well, if we assume the radius of the well is 1 foot.WRITE-UP. Please type your write-up following the format given below. It is okay to write inequations and figures by hand. Please be brief and concise. Also turn in a disk containing yourspreadsheet(s).Problem Description. Present the modeling objective: to compare the results of all the solutionsby focusing on the head at the center of the island. Describe the conceptual model. Present themathematical model (governing equation and boundary conditions) for each of the following versionsof the problem: the 1-D confined aquifer, the 2-D confined aquifer, the 1-D unconfined aquifer andthe 2-D unconfined aquifer. Show the grids used for the numerical solutions and identify thenumerical solution method. Results. You have 3 analytical solutions for this problem. Use the template below to tabulate theresults from each solution by listing the head in the center of the island when R = 0.00305 ft/day.Also give the results from problem 3.1 for the 1-D confined aquifer for which you assumed that thehead was 20 ft and then calculated the recharge rate needed to maintain that head. You also have 4 numerical solutions. Fill in the tables (below) showing the head in the center of theisland (reported to the 4th decimal place), and the water balance fluxes and the water balance errorfor each numerical solution.Discussion. Discuss and explain the differences in the solutions as reported in the tables presented inyour results section. Conclusions. Make some conclusions about your modeling work relevant to the modelingobjective.Appendix. Include the derivations of the analytical solutions for the 1-D confined and 1-Dunconfined cases. Also include a printout of each numerical model run. Label each printout. Besure your printout includes the distribution of head in the upper right hand quadrant of the island,convergence criterion (error tolerance), water balance inflow and outflow, difference in the waterbalance and percent error in the water balance. Analytical SolutionsSolution Type Head at center of island (ft)1d confined21d unconfined2d confinedNumerical SolutionsConceptual Modelx, y(ft)Error Tolerance (ft)Head at centerof island (ft)Water BalanceError (%)Inflow, Outflow (ft3/day)2D confined 4000 1E-42D confined 1000 1E-42D unconfined 1000 1E-4Pumping SolutionConceptual Modelx, y(ft)Error Tolerance (ft)Head in node at centerof island Head in well (ft)Water BalanceError (%)Inflowft3/dayQ(well)ft3/dayOutflow tothe ocean(ft3/day)2D unconfined1000 1E-4Analytical solution of the Island Recharge Problem using the equation given in Problem 3.2 (p. 65) inWang and Anderson. Note that the solution is for the lower right hand quadrant of the problem domain.20.00 19.86 19.46 18.78 17.82 16.59 15.08 13.29 11.22 8.86 6.20 3.25 .0019.98 19.85 19.44 18.76 17.81 16.58 15.07 13.28 11.21 8.85 6.20 3.25 .00319.93 19.80 19.39 18.71 17.76 16.54 15.04 13.25 11.19 8.83 6.18 3.24 .0019.85 19.71 19.31 18.64 17.69 16.47 14.98 13.20 11.14 8.80 6.16 3.23 .0019.73 19.59 19.19 18.52 17.58 16.37 14.89 13.13 11.08 8.75 6.13 3.22 .0019.57 19.43 19.04 18.38 17.45 16.25 14.78 13.03 11.00 8.69 6.09 3.19 .0019.37 19.24 18.84 18.19 17.27 16.09 14.63 12.91 10.90 8.61 6.04 3.17 .0019.12 18.99 18.61 17.96 17.06 15.89 14.46 12.76 10.78 8.52 5.97 3.14 .0018.83 18.70 18.32 17.69 16.81 15.66 14.25 12.58 10.63 8.41 5.90 3.10 .0018.48 18.36 17.99 17.37 16.51 15.38 14.01 12.37 10.46 8.27 5.81 3.05 .0018.08 17.96 17.60 17.00 16.15 15.06 13.72 12.12 10.25 8.12 5.70