Next Article in Journal
Assessment of Cyanobacterial Chlorophyll A as an Indicator of Water Quality in Two Wetlands Using Multi-Temporal Sentinel-2 Images
Previous Article in Journal
Flood Vulnerability Mapping Using MaxEnt Machine Learning and Analytical Hierarchy Process (AHP) of Kamrup Metropolitan District, Assam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Fuzzy Analytical Solution for the Case of a Semi-Infinite Unconfined Aquifer †

by
Christos Tzimopoulos
1,
Nikiforos Samarinas
1,*,
Kyriakos Papadopoulos
2 and
Christos Evangelides
1
1
Department of Rural and Surveying Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
2
Department of Mathematics, Kuwait University—Khaldiya Campus, Safat 13060, Kuwait
*
Author to whom correspondence should be addressed.
Presented at the 7th International Electronic Conference on Water Sciences, 15–30 March 2023; Available online: https://ecws-7.sciforum.net.
Environ. Sci. Proc. 2023, 25(1), 70; https://doi.org/10.3390/ECWS-7-14303
Published: 3 April 2023
(This article belongs to the Proceedings of The 7th International Electronic Conference on Water Sciences)

Abstract

:
The solution to the second-order fuzzy unsteady nonlinear partial differential one-dimensional Boussinesq equation is examined. The physical problem concerns unsteady flow in a semi-infinite, unconfined aquifer bordering a lake. There is a sudden rise and subsequent stabilization in the water level of the lake; thus, the aquifer is recharging from the lake. The fuzzy solution is presented by a simple algebraic equation transformed in a fourth-degree polynomial approximation for the head profiles. In order to solve this equation, the initial and boundary conditions, as well as the numerous soil properties, must be known. A fuzzy approach is used to solve the problem since the aforementioned auxiliary conditions are vulnerable to various types of uncertainty resulting from human and machine errors. The physical problem described by a partial differential equation and the generalized Hukuhara derivative and the application of this theory for the partial derivatives were chosen as solving methods. In order to evaluate the accuracy and effectiveness of the suggested fuzzy analytical method, this study compares the findings of fuzzy analysis to those obtained using the Runge–Kutta method. This comparison attests to the accuracy of the former. Additionally, this results in a fuzzy number for water level profiles as well as for the water volume variation, whose α-cuts, provide according to Possibility Theory, the water levels and the water volume confidence intervals with probability p = 1 − α.

1. Introduction

The horizontal water flow concerning unconfined aquifers without precipitation is described by the one-dimensional second-order unsteady nonlinear partial differential equation, called the Boussinesq equation:
h t = K S x ( h h x ) ,
where K = hydraulic conductivity (LT−1), S = effective porosity(L0T0), h = piezometric head (L), x = horizontal coordinate (L), t = time (T).
Boussinesq (1904) [1] first proposed the above equation according to the assumption that the horizontal component of velocity ux does not change with depth and is a function of x and t, while the inertial forces are negligible. A unique solution to this nonlinear equation was published by Boussinesq in the “Journal de Mathématiques Pures et Appliquées” in 1904. With boundary conditions such as those of soil drained by drains placed in the impermeable substratum, Boussinesq solution dealt with the case of an aquifer atop an impermeable layer. Using the small disturbance method, Polubarinova-Kochina (952, 1962) [2,3] published a solution to Boussinesq’s equation. By utilizing polynomial approximation and similarity transformation, Tolikas et al. (1984) [4] found an approximate solution. A weighted residual approach was used by Lockington (1997) [5] to provide an easily applicable analytical solution. Due to an abrupt change in the head at the origin, this approach was used for both the recharging and discharging of an unconfined aquifer. By using Adomian’s decomposition method, Moutsopoulos (2010) [6] arrived at a simple series solution with a limited number of terms while he also performed a benchmark test, which demonstrated the benefits of his solution. The problem of flow in a one-dimensional semi-infinite horizontal aquifer with an initially dry head and a power-law function of time at the origin was examined by Lockington et al. in 2000 [7]. In accordance with the numerical outcomes, an approximative quadratic solution was developed. Using the traveling wave method, Basha (2013) [8] obtained a nonlinear solution with an easily applicable logarithmic form. The solution allows for the results to be of practical value in hydrology and is adaptable to any flow situation, whether it be recharge or discharge condition. There are also algebraic formulae for the propagation front velocity, the location of the wetting front, and the linkage between the characteristics of the aquifer. The nonlinear Boussinesq equation was given by a series solution by Chor et al. (2013) [9] in terms of the Boltzmann transform in a semi-infinite domain. An approximate solution was recently obtained by Hayek (2019) [10], who introduced an empirical function with four parameters. Using Microsoft Excel Solver, a numerical fitting approach was used to acquire the parameters. An approximate analytical solution for the recharge and discharge of a homogeneous unconfined aquifer was published by Tzimopoulos et al. in 2022 [11]. Numerous other studies [12,13,14,15,16] offer helpful clarification on the solution, offering a way for testing and accuracy of the numerical methods.
The definition of the initial flow condition, the method of linearizing the Boussinesq equation, the definition of drain spacing and hydraulic conductivity, boundary conditions, etc., are just a few examples of the ambiguities and uncertainties that the physical problem described by the Boussinesq equation presents [17,18]. Without taking into account the ambiguities and uncertainties of the groundwater flow problems, wrong management decisions could be made, leading to a number of significant negative environmental, social, and economic effects. Fuzzy algorithms were used to solve this problem for all the aforementioned reasons.
The fuzzy logic theory is a useful tool for modeling ambiguity, developed by Lofti Zadeh (1965) [19]. Its development has had a significant impact on both theoretical problems [20,21,22,23] as well as engineering and hydraulic problems [24]. To solve fuzzy differential equations, some analytical and numerical approaches have recently been put forth. Chang and Zadeh (1972) [25] first introduced the concept of fuzzy derivative, while Dubois and Prade (1982) [26] followed by using the extension principle in their approach. Fuzzy differential functions were studied by Puri and Ralescu (1983) [27], who extended Hukuhara’s derivative (H-derivative) [28] of a set of values appearing in fuzzy sets. Kaleva (1987, 1990) [29,30] and Seikkala (1987) [31] developed the fuzzy initial value problem, while Abbasbandy and Allahviranloo (2002) [32] presented a numerical algorithm for solving fuzzy ordinary differential equations based on the second Taylor method. However, this method presented certain drawbacks, and in many cases, this solution was not a good generalization of the classic case. The generalized Hukuhara differentiability (gH—differentiability) was introduced by [33,34], overcoming this drawback. This new derivative is defined for a larger class of fuzzy functions than the Hukuhara derivative. Allahviranloo et al. (2015) [35] introduced the (gH-p) differentiability for partial derivatives as an extension of the above theory. The gH-p differentiability was used by Tzimopoulos et al. (2018, 2020) [17,36], providing a fuzzy linear analytical solution to a parabolic partial differential equation, and also Tzimopoulos et al. (2018) [18] obtained a fuzzy linear analytical solution to the Boussinesq equation in the case of an unconfined aquifer problem.
In this paper, a comparison between the fuzzy analytical nonlinear Boussinesq equation and the Runge–Kutta numerical method is presented for the proposed fuzzy analytical solution accuracy evaluation and effectiveness check. This comparison attests to the accuracy of the former. Additionally, an application of the Possibility Theory [37,38] to the α-cuts of the water level profiles, as well as of the water volume variation, provides the water levels and the water volume confidence intervals with probability p = 1 − α. Thus, a combination of fuzzy theory with the Possibility Theory allows managers and engineers to solve practical hydraulic problems, making the right decision.

2. Constructing the Fuzzy Model and Solution

2.1. Crisp Model

As mentioned in the introduction to the case of one-dimensional horizontal flow, the Equation (1) corresponds to Boussinesq equation
The initial and boundaries conditions are as follows:
t = 0 , h ( x , 0 ) = h 0 , t > 0 , h ( 0 , t ) = h 1 , h x ( x , t ) = h 0
The solution [11] of the above Equation (1) is as follows:
h = h 0 + ( h 0 h 1 ) Ω ( μ , ξ )
where
Ω ( μ , ξ ) = F ( ξ ) Φ ( ξ ) , Φ ( ξ ) = e r f c ( ξ ) , ξ = x 2 K h 1 t S , = ( h 0 h 1 ) / h 1 , μ = h 1 / h 0 , F ( ξ ) = 1 π + ( 1 2 + 1 π ) Φ 1 π ( 1 Φ ) ξ e ξ 2 + 1 π ( 1 e 2 ξ 2 ) 1 2 Φ 2 .

2.2. Fuzzy Model

We write Εquation (1), in its fuzzy form as follows:
h ˜ t = K S x ( h ˜ H ˜ ˜ x ) = K S { ( h ˜ x ) 2 + h ˜ 2 h ˜ x 2 } ,
with the new boundary and initial conditions:
Initial   conditions Boundary   conditions t = 0 , h ˜ ( x , 0 ) = h ˜ 0 , t > 0 , h ˜ ( 0 , t ) = h ˜ 1 , h ˜ x ( x , t ) = h ˜ 0
By converting the aforementioned fuzzy problem into a system of second order crisp boundary value problems—(referred as the corresponding system for the fuzzy problem)—and applying the theory of [18,35,39,40], we are able to solve the fuzzy problem (2), subject to the boundary and initial conditions (3). As a result, for the fuzzy problem, eight crisp BVPs systems are feasible under identical initial and boundary conditions.
Note: We shall now turn our attention to the first system’s solution as it offers a practical solution to the case of the lake refilling the aquifer.

2.3. Solution of the First System

The nonlinear one-dimensional horizontal flow equations related to the first case (left boundary) and to the second case (right boundary) are provided with the following expressions:
First   case Second   case h t = K S { h 2 h x 2 + ( h x ) 2 } h + t = K S { h + 2 h + x 2 + ( h + x ) 2 }

3. Results and Discussion

We have used the data of Lockington (1997) [5], which means K = hydraulic conductivity = 20 m/d, S = effective porosity = 0.27, h1 = 3m, h0 = 2m. For the first case, the solution is as follows [11]:
h ( μ , ξ ) = h 0 + ( h 0 h 1 ) Ω ( μ , ξ )
Figure 1 illustrates the water level profiles for t = 5d of the new analytical method vs Runge–Kutta method (reference solution) and for α = 0, 1 value. The two methods approach each other closely. In addition, the water table in a period of 5 days approaches 125 m in length. If we consider also the hydraulic conductivity value (K = 20 m/d), the results are absolutely reasonable regarding the physical problem. Figure 2 illustrates the stored volume variation vs. time, and Figure 3 illustrates the membership function of V ˜ ( x , t ) for t = 5d. According to possibility theory [37,38], every function [ V ˜ ] estimates the crisp function V, and the α-cut { V ˜ } a = [ V a ,   V a + ] should be interpreted as the confidence intervals of V with a probability p 1 a . In this regard, in Figure 3, it is seen that for α = 0.05, the value of water volume lies in the interval [5.226, 14.310] with a probability higher than 95%, according to the possibility theory.

4. Conclusions

Undoubtedly, groundwater flow problems involve a number of ambiguities and uncertainties, making the use of differential equations even more difficult to solve. However, nowadays, the opportunity is given through the fuzzy partial differential equations theories to include these uncertainties in the final calculations and to provide more accurate results supporting the sustainability of groundwater as well as the researchers and engineers to make better decisions and plannings.
This work presents an innovative analytical solution in the nonlinear Boussinesq equation, which describes the groundwater’s unsteady flow. According to the results, the proposed solution completely coincides with the Runge–Kutta method results used as a reference solution for comparison reasons in order to prove the accuracy and reliability of the proposed analytical solution. The volume membership function (Figure 3) could support decision-makers and planners with a higher degree of confidence than the previous years, thanks to the possibility theory.

Author Contributions

Conceptualization, C.T.; methodology, C.T. and N.S.; software, C.T.; validation, K.P. and C.E.; writing—original draft preparation, C.T.; writing—review and editing, N.S. and C.E.; visualization, C.T.; supervision, C.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boussinesq, J. Recherches théoriques sur l’écoulement des nappes d’eau infiltrées dans le sol et sur le débit des sources. J. De Mathématiques Pures Appliquées 1904, 10, 5–78. [Google Scholar]
  2. Polubarinova-Cochina, P.Y. On unsteady motions of groundwater during seepage from water reservoirs. PMM (Prinkladaya Matematica I Mekhanica) 1949, 13, 2. [Google Scholar]
  3. Polubarinova-Cochina, P.Y. Theory of Groundwater Movement; Moscow: Gostekhizdat, Old Soviet Union (Now Russia); de Wiest, R.J.M., Translator; Princeton University Press: Princeton, NJ, USA, 1962. (In Russia) [Google Scholar]
  4. Tolikas, P.K.; Sidiropoulos, E.; Tzimopoulos, C.D. A simple analytical solution for the Boussinesq one-dimensional groundwater flow equation. Water Resour. Res. 1984, 20, 24–28. [Google Scholar] [CrossRef]
  5. Lockington, D.A. Response of unconfined aquifer to sudden change in boundary head. J. Irrig. Drain. Eng. 1997, 123, 24–27. [Google Scholar] [CrossRef]
  6. Moutsopoulos, K.N. The analytical solution of the Boussinesq equation for flow induced by a step change of the water table elevation revisited. Transp. Porous Media 2010, 85, 919–940. [Google Scholar] [CrossRef]
  7. Lockington, D.A.; Parlange, J.Y.; Parlange, M.B.; Selker, J. Similarity solution of the Boussinesq equation. Adv. Water Resour. 2000, 23, 725–729. [Google Scholar] [CrossRef] [Green Version]
  8. Basha, H.A. Traveling wave solution of the Boussinesq equation for groundwater flow in horizontal aquifers. Water Resour. Res. 2013, 49, 1668–1679. [Google Scholar] [CrossRef]
  9. Chor, T.; Dias, N.L.; de Zarate, A.R. An exact series and improved numerical and approximate solutions for the Boussinesq equation. Water Resour. Res. 2013, 49, 7380–7387. [Google Scholar] [CrossRef]
  10. Hayek, M. Accurate approximate semi-analytical solutions to the Boussinesq groundwater flow equation for recharging and discharging of horizontal unconfined aquifers. J. Hydrol. 2019, 570, 411–422. [Google Scholar] [CrossRef]
  11. Tzimopoulos, C.; Papadopoulos, K.; Evangelides, C.; Spyrides, A. Analytical solution of nonlinear Boussinesq equation. Desalination Water Treat. 2022, 24, 1127–1147. [Google Scholar] [CrossRef]
  12. Chen, Z.-X.; Bodvarsson, G.S.; Witherspoon, E.A.; Yortsos, Y.C. An integral equation formulation for the unconfined flow of groundwater with variable inlet conditions. Transp. Porous Media 1995, 18, 15–36. [Google Scholar] [CrossRef]
  13. Parlange, J.-Y.; Hogarth, W.L.; Govindaraju, R.S.; Parlange, M.B.; Lockington, D. On an Exact Analytical Solution of the Boussinesq Equation. Transp. Porous Media 2000, 39, 339–345. [Google Scholar] [CrossRef] [Green Version]
  14. Pistiner, A. Similarity solution to unconfined flow in an aquifer. Transp. Porous Media 2008, 71, 265–272. [Google Scholar] [CrossRef]
  15. Olsen, J.S.; Telyakovskiy, A.S. Polynomial approximate solutions of a generalized Boussinesq equation. Water Resour. Res. 2013, 49, 3049–3053. [Google Scholar] [CrossRef]
  16. Bartlett, M.S.; Porporato, A. A Class of exact solutions of the Boussinesq equation for horizontal and sloping aquifers. Water Resour. Res. 2018, 54, 767–778. [Google Scholar] [CrossRef]
  17. Tzimopoulos, C.; Papadopoulos, K.; Evangelides, C.; Papadopoulos, B. Fuzzy solution to the unconfined aquifer problem. Water 2018, 11, 54. [Google Scholar] [CrossRef] [Green Version]
  18. Tzimopoulos, C.; Papadopoulos, K.; Evangelides, C. Fuzzy Solution to the Second Order Unsteady Partial Differential Equation. In Proceedings of the International Conference of Numerical Analysis and Applied Mathematics (ICNAAM 2018), Rhodes, Greece, 13–18 September 2018; Volume 2116, pp. 440004-1–440004-5. [Google Scholar] [CrossRef]
  19. Zadeh, L.A. Fuzzy Sets. Inf. Control 1965. [Google Scholar] [CrossRef] [Green Version]
  20. Bhaskar, V.; Lakshikantham, T.; Devi, V. Revisiting fuzzy differential equations. Nonlinear Anal. 2004, 58, 351–358. [Google Scholar]
  21. Nieto, J.J.; Rodríguez-López, R. Bounded solutions for fuzzy differential and integral equations. Chaos Solitons Fractals 2006, 27, 1376–1386. [Google Scholar] [CrossRef]
  22. Aminikhah, H. Approximate analytical solution for the one-dimensional nonlinear Boussinesq equation. Int. J. Numer. Methods Heat Fluid Flow 2015, 25, 831–840. [Google Scholar] [CrossRef]
  23. Chen, F.; Qiu, X.; Alattas, K.A.; Mohammadzadeh, A.; Ghaderpour, E. A New Fuzzy Robust Control for Linear Parameter-Varying Systems. Mathematics 2022, 10, 3319. [Google Scholar] [CrossRef]
  24. Guo, M.; Xue, X.; Li, R. The oscillation of delay differential inclusions and fuzzy biodynamics models. Math. Comput. Model. 2003. [Google Scholar] [CrossRef]
  25. Chang, S.S.L.; Zadeh, L.A. On Fuzzy Mapping and Control. IEEE Trans. Syst. Man Cybern. 1972. [CrossRef]
  26. Dubois, H.; Prade, D. Towards fuzzy differential calculus: Part 3, differentiation. Fuzzy Sets Syst. 1982, 8, 30–34. [Google Scholar] [CrossRef]
  27. Puri, M.L.; Ralescu, D.A. Differentials of fuzzy functions. J. Math. Anal. Appl. 1983, 91, 552–558. [Google Scholar] [CrossRef] [Green Version]
  28. Hukuhara, M. Integration des Applications Measurables dont la Valeur est un Compact Convexe. Funkc. Ekvacioj 1967, 10, 205–233. [Google Scholar]
  29. Kaleva, O. Fuzzy differential equations. Fuzzy Sets Syst. 1987, 24, 301–307. [Google Scholar] [CrossRef]
  30. Kaleva, O. The cauchy problem for fuzzy differential equations. Fuzzy Sets Syst. 1990, 24, 389–396. [Google Scholar] [CrossRef]
  31. Seikkala, S. On the fuzzy initial value problem. Fuzzy Sets Syst. 1987, 24, 319–330. [Google Scholar] [CrossRef]
  32. Abbasbandy, S.; Allahviranloo, T. Numerical solution of fuzzy differential equation. Math. Comput. Appl. 2002, 7, 41–52. [Google Scholar] [CrossRef] [Green Version]
  33. Bede, B.; Gal, S.G. Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar] [CrossRef]
  34. Stefanini, L.; Bede, B. Generalized Hukuhara differentiability of interval-valued functions and interval differential equations. Nonlinear Anal. Theory Methods Appl. 2009, 71, 1311–1328. [Google Scholar] [CrossRef] [Green Version]
  35. Allahviranloo, T.; Gouyandeh, Z.; Armand, A.; Hasanoglu, A. On fuzzy solutions for heat equation based on generalized Hukuhara differentiability. Fuzzy Sets Syst. 2015, 265, 1–23. [Google Scholar] [CrossRef]
  36. Tzimopoulos, C.; Papaevangelou, G.; Papadopoulos, K.; Evangelides, C. Fuzzy Analytical Solution to Vertical Infiltration. J. Softw. Eng. Appl. 2020, 13, 41–46. [Google Scholar] [CrossRef] [Green Version]
  37. Dubois, D.; Prade, H. When upper probabilities are possibility Measures. Fuzzy Sets Syst. 1992, 49, 65–74. [Google Scholar] [CrossRef]
  38. Mylonas, N. Applications in Fuzzy Statistic and Approximate Reasoning. Ph.D. Thesis, Dimokritos University of Thrace, Komotini, Greece, 2022. (In Greek). [Google Scholar]
  39. Khastan, A.; Nieto, J.J. A boundary value problem for second order fuzzy differential equations. Nonlinear Anal. Theory Methods Appl. 2010, 72, 3583–3593. [Google Scholar] [CrossRef]
  40. Bede, B.; Stefanini, L. Generalized differentiability of fuzzy-valued functions. Fuzzy Sets Syst. 2013, 230, 119–141. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Water level profiles for t = 5d.
Figure 1. Water level profiles for t = 5d.
Environsciproc 25 00070 g001
Figure 2. Stored volume variation vs. time.
Figure 2. Stored volume variation vs. time.
Environsciproc 25 00070 g002
Figure 3. Membership function of V ˜ ( x , t ) .
Figure 3. Membership function of V ˜ ( x , t ) .
Environsciproc 25 00070 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tzimopoulos, C.; Samarinas, N.; Papadopoulos, K.; Evangelides, C. Fuzzy Analytical Solution for the Case of a Semi-Infinite Unconfined Aquifer. Environ. Sci. Proc. 2023, 25, 70. https://doi.org/10.3390/ECWS-7-14303

AMA Style

Tzimopoulos C, Samarinas N, Papadopoulos K, Evangelides C. Fuzzy Analytical Solution for the Case of a Semi-Infinite Unconfined Aquifer. Environmental Sciences Proceedings. 2023; 25(1):70. https://doi.org/10.3390/ECWS-7-14303

Chicago/Turabian Style

Tzimopoulos, Christos, Nikiforos Samarinas, Kyriakos Papadopoulos, and Christos Evangelides. 2023. "Fuzzy Analytical Solution for the Case of a Semi-Infinite Unconfined Aquifer" Environmental Sciences Proceedings 25, no. 1: 70. https://doi.org/10.3390/ECWS-7-14303

APA Style

Tzimopoulos, C., Samarinas, N., Papadopoulos, K., & Evangelides, C. (2023). Fuzzy Analytical Solution for the Case of a Semi-Infinite Unconfined Aquifer. Environmental Sciences Proceedings, 25(1), 70. https://doi.org/10.3390/ECWS-7-14303

Article Metrics

Back to TopTop