### Abstract

Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

Original language | English |
---|---|

Pages (from-to) | 278-282 |

Number of pages | 5 |

Journal | Journal of Hydrology |

Volume | 523 |

DOIs | |

Publication status | Published - 1 Apr 2015 |

### Fingerprint

### Keywords

- Exponential integral
- Latin hypercube sampling
- Orthogonal array
- Orthogonal array-based latin hypercube sampling
- Well function

### ASJC Scopus subject areas

- Water Science and Technology

### Cite this

**Approximation of the exponential integral (well function) using sampling methods.** / Baalousha, Husam.

Research output: Contribution to journal › Article

}

TY - JOUR

T1 - Approximation of the exponential integral (well function) using sampling methods

AU - Baalousha, Husam

PY - 2015/4/1

Y1 - 2015/4/1

N2 - Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

AB - Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

KW - Exponential integral

KW - Latin hypercube sampling

KW - Orthogonal array

KW - Orthogonal array-based latin hypercube sampling

KW - Well function

UR - http://www.scopus.com/inward/record.url?scp=84922636311&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84922636311&partnerID=8YFLogxK

U2 - 10.1016/j.jhydrol.2015.01.062

DO - 10.1016/j.jhydrol.2015.01.062

M3 - Article

AN - SCOPUS:84922636311

VL - 523

SP - 278

EP - 282

JO - Journal of Hydrology

JF - Journal of Hydrology

SN - 0022-1694

ER -