Next Article in Journal
Assessing the Influences of Land Use Change on Groundwater Hydrochemistry in an Oasis-Desert Region of Central Asia
Next Article in Special Issue
Discrimination between Pore and Throat Resistances against Single-Phase Flow in Porous Media
Previous Article in Journal
Assessing the Groundwater Reserves of the Udaipur District, Aravalli Range, India, Using Geospatial Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flow Discharge Prediction Study Using a CFD-Based Numerical Model and Gene Expression Programming

1
Department of Water Engineering, Urmia University, Urmia 57561, Iran
2
School of Civil Engineering, College of Engineering, University of Tehran, Tehran 13145, Iran
3
Center for Artificial Intelligence Research and Optimization, Torrens University Australia, Brisbane 4006, Australia
*
Author to whom correspondence should be addressed.
Water 2022, 14(4), 650; https://doi.org/10.3390/w14040650
Submission received: 17 December 2021 / Revised: 5 February 2022 / Accepted: 14 February 2022 / Published: 19 February 2022
(This article belongs to the Special Issue Numerical Modelling of Single and Multi-Phase Flow)

Abstract

:
The significance of spillways is to allow the flood to be safely discharged from downstream. There is a strong correlation between the poor design of spillways and the failures of dams. In order to address this concern, the present study investigates the flow over the Nazloo-ogee spillway using the CFD 3D numerical model and an artificial intelligence method called Gene Expression Programming (GEP). In a physical model, discharge and flow depths were calculated for 21 different total heads. Among different turbulence models, the R N G turbulence model achieved the maximum compatibility in computational fluid dynamic simulation. In addition, GEP was used to estimate Q, in which 70% of collected data was dedicated to training and 30% to testing. R 2 , R M S E , and M A E were obtained as performance criteria, and the new mathematical equation for the prediction of discharge was obtained using this model. Finally, the numerical model and GEP outputs were compared with the experimental data. According to the results, the numerical model and GEP exhibited a high level of correspondence in simulating flow over an ogee-crested spillway.

1. Introduction

Flood control, reliable water supply, navigation, recreation, and hydroelectric power generation are the most important reasons for dam development around the world [1,2,3,4]. Spillways are hydraulic structures built on dams to convey excess flood flow beyond the dam’s capacity. In order to avoid serious damage, spillways should be strong, reliable, and high-structures. Accordingly, design and construction play a role in dam spillways. Flood risks can be prevented if a dam spillway is properly designed and constructed [5,6,7]. Given their excellent hydraulic features, ogee-crested spillways are among the most investigated hydraulic structures. The U.S. Bureau of Reclamation, USBR, (1987) performed detailed laboratory experiments to investigate the behavior of water flow over a spillway, which resulted in the development and publication of spillway design manuals [8]. Physical studies have a number of problems, including high construction costs, the considerable time required for development and testing, and the presence of errors in results due to scaling effects [9,10,11]. The available computing capacity and algorithm improvements helped in finding solutions for a variety of problems, such as flow over the spillway. Numerical models provide a tool the fundamental design of spillways to identify operational concerns at a lower cost and in a shorter time [12,13,14]. Savage and Johnson investigated both pressure and discharge in the physical and numerical models in an ogee spillway, which indicated a satisfactory performance [15]. Kim and Parkn investigated the pressure and velocity distribution over an ogee spillway using the R N G K ϵ turbulence model in the CFD model, taking into account the surface roughness effects [16]. Peltier et al. used experimental modeling of an ogee spillway to estimate velocity and pressure for heads that were much higher than the design head and validated the data with experimental data [17]. Nowadays, numerous researchers have utilized the K ϵ turbulence model in their studies, including (Jahad et al. [18]; Aydin et al. [19]; Wan et al. [20]). Moreover, R N G K ϵ model has been applied in (Ghanbari and Heidarnejad [21]; Bayon et al. [22]; Valero et al. [23]), and remarkable accuracy in stimulating the overflow spillway. Artificial intelligence (AI) is a field of computer science that focuses on developing machines that can engage in intelligent behaviors. In recent years, machine learning techniques (e.g., gene expression programming (GEP), adaptive neuro-fuzzy inference system (ANFIS) and artificial neural network ( A N N s ), have been used in a wide range of publications. These methods are outstanding forecasting tools that have been used in a variety of civil engineering, hydraulic, and hydrological research in the last decade. Ferreira proposed gene-expression programming as a new adaptive approach for solving issues [24]. Yildiz et al. investigated discharge and flow depth over an ogee spillway using the ANFIS andCFD models in Flow-3D software. According to their results, therewas a reasonable agreement between the physical, numerical, and ANFIS artificial neural networks models [25]. The GEP technique has been used to represent a variety of water resource system components. Ebtehaj et al. utilized this method to predict the discharge coefficient in rectangular side weirs [26].
Using soft computing techniques, Kisi et al. examined the prediction of lateral outflow over triangular labyrinth side weirs under subcritical conditions [27]. Salmasi found a new equation for predicting discharge coefficients in an ogee weir, using gene expression programming and multiple regression techniques. The results demonstrated that the GEP technique was more successful than the regression equations [28]. Khan et al. employed a GEP method, to establish a functional relationship for bridge pier scour. The performanceof GEP was compared to that of other Al-based techniques, such as artificial neural networks (ANNs) and conventional reqression-based techniques [29]. Roushangar et al. used GEP and ANN techniques to estimate energy dissipation over stepped spillways, and the findings showed that GEP and ANN were extremely beneficial and encouraging in these circumstances [30]. Using GEP, Baxgatur and Onen created flood routing prediction models [31]. Bertonse et al. used GEP to simulate the concentration of dissolved oxygen (Do) in lakes [32]. New algorithms and models, particularly those based on soft computing, enable researchers to address the most complex systems in a number of ways [33,34]. In different engineering fields, the forecast methods that are not dependent on physics equations, including remote sensing methods and AI-based approaches likethe GEP method, are becoming more common [35,36]. The GEP has the advantage of being explicit in its formulation. This provides some insight into the nature of the phenomenon under investigation. It is simple to apply in reality [37]. The GEP is an artificial intelligence technique that utilizes key principles of genetic algorithms (GA) and genetic programming (GP) to create a calculation algorithm for forecasting a certain phenomenon. It mimics biological evolution [38]. Genetic programming is based on evolutionary principles developed for mathematical modeling. For solving regression and classification problems, GP provides several computer programs. The optimal values of some predefined parameters are obtained by GA, while finding both the best models and best parameters for a set of variables produced by GP is based on the Darwinian evolution theory [39]. There are also several other applications of such techniques in the literature that have been investigated insightfully in recent years [40,41].
A review of the literature indicates a lack of exclusive comprehensive research on the use of CFD and GEP models to estimate discharge values over an ogee-crested spillway with a pier. In the present study, the statistical performance of the model of an ogee-crested spillway with a pier was evaluated using error criteria, such as the root-mean-square error ( R M S E ) , which is a significant criterion, as well as the determination coefficient ( R 2 ) and mean absolute error ( M A E ) and a new mathematical equation was developed for predicting the spillway’s discharge. The results of the GEP were then compared to those obtained from the physical model. In addition to this analysis, knowledge of discharge and flow depth can be valuable for hydraulic engineering research as well as determining the GEP approach’s capability of forecasting a new output variable. The CFD model was validated with the experimental and Flow-3D data. The results showed that all of the model’s outputs overlap and are accurate when it comes to tackling fluid problems.
This paper consists of four sections. Section 2 contains a brief discussion of the methods and components utilized to simulate the ogee-crested spillway as well as the provided gene expression programming and a Flow-3D numerical model. Section 3 discusses the findings of the investigation regarding numerical and experimental data comparisons, and mathematical equations. Finally, Section 4 summarizes the paper’s main findings.

2. Materials and Methods Methods

The head, the inclination of the upstream face of the overflow section, and the height of the overflow section above the entrance channel’s floor determine the shape of the nappe profile over an ogee crest, which influences the velocity of approach to the crest. The general discharge equation for an ogee-crested spillway is as follows:
Q = C o L H o 3 2
where Q is the design discharge ( m 3 s ), C o is the variable coefficient of discharge for free-flow condition ( m 0.5 s ), L is the effective length of the crest (m), H o is the design head on the crest(m), including the velocity of approach, h a (m). Figure 1 demonstrates the ogee profile.
In this study, a physical model of a standard ogee spillway was developed at the Water Research Institute of the Ministry of Energy, Iran. On the west side of lake Urmia, on the Nazloo river, the Nazloo dam spillway is under construction with a width of 42 m, a crest height of 6 m, and a crest level of 1492 m above sea level. The hydraulic performance of the Nazloo spillway was studied using a physical model with a scale of 1:40 made of Plexiglas Figure 2. The effective highest measured at a distance of 4 m from the top of the crest of the water surface. The pier over the crest had a height of 3.51 m, a width of 2 m, and a length of 10 m.
Crest piers cause the flow to be constructed, reduce the effective length of the crest, and decrease the discharge compared to an uncontrolled crest. The shape and location of the pier nose are shown in Figure 3, as well as the thickness of the pier, the head relative to the design head, and the approach velocity, which all influence the pier construction coefficient, K p [42].
L e = L 2 ( N K p + K a ) H o
where L e is the effective length of the crest for calculatingthe discharge, L is the net clear length of the spillway crest, N is the number of piers, K p is the pier construction coefficient, K a is the abutment contraction coefficient and H o is the total design head on the crest including velocity head. For round-nosed piers, the K p and K a are set to 0.01 and 0.05, respectively. Dimensional analysis is a fundamental tool used in experimental research to determine dimensionless quantities. For this purpose, characteristics of the effect of the ogee spillway discharge should be identified at first. Then, Buckingham’s theory should be utilized to estimate the dimensionless parameters, which can be used to analyze their impact on the spillway’s discharge and identify the logical correlation between them [21]. The geometric parameters and flow characteristics might be used to calculate the discharge capacity of an ogee spillway. The quantities influencing ogee crest discharge can be summarised as follows:
f ( H o , L e , H d , P , g , μ , ρ , σ ) = 0
where f is a functional symbol; H e and H o are the total upstream head and the design head, respectively; P is the ogee crest upstream spillway height; μ and ρ are dynamic viscosity and density, respectively; and σ is the surface tension. The possible effects of surface tension on discharge were small in all experiments H > 30 mm. The Weber number was excluded from the analysis. Furthermore, since the flow was turbulent, the viscosity impact was minor compared to gravity. As a result, the Reynolds Number effect was also removed from the analysis. In this study, the ultimate relationship between the dimensionless parameters impacting the ogee spillway discharge would be stated as:
Q = f ( H e H o , L e H o , H e p )
The experiments were performed with 21 total heads over the spillway, ranging from H e H o = 0.123 to H e H o = 1.22.

2.1. Flow-3D Numerical Modeling

FLOW-3D software is a useful tool for analyzing complex fluid issues such as free-surface transient three-dimensional flow problems with complex geometries. This software uses the finite volume method along with regular rectangular grids. In order to simulate turbulence in hydraulic problems, two-equation models are commonly utilized. To derive time-averaged Reynolds equations, the Renormalization Group R N G turbulence model was used in this study. Furthermore, the numerical problems were solved using the FLOW-3D software, and the transient governing equations were numerically solved using the finite volume model. The geometry of the finite volume approach was defined in this software utilizing the Fractional Area-Volume Obstacle Representation (FAVOR) algorithm, as shown in Figure 4. The in-field obstacles in the computational cells were treated as a fractional value between 0 and 1. Therefore, if the entire cell was filled with obstacles, the fractional value of the area-value was equal to 1. The volume-of-fluid (VOF) algorithm was used to determine the flow’s free surface [43].
For the considered field, the FLOW-3D numerical model provided a three-dimensional structural grid made up of cuboid cells as shown in Figure 5. As a result, a three-dimensional model was created in AutoCAD software to build up the geometry of the models with a stereolithography (STL) file, based on the laboratory models’ specifications. The data was then loaded into FLOW-3D, where the gird was generated using VOF and FAVOR, and the boundaries and computational network were determined. The considered area was constructed using VOF and FAVOR methods after loading the geometric data into the software.
The boundary conditions’ most critical aspect is to create a flow situation similar to the physical status. Each type of boundary condition in the Flow-3D software can be used for the unique condition of the models. The X-minimum boundary condition, specified pressure, was applied in this study. An outflow was chosen for the X-maximum boundary condition since there was nothing to calculate at the flume’s end. Both (Y-min) and (Y-max) were computed as symmetry, with the bottom (Z-min) as a wall boundary condition and the top (Z-max) as a symmetry boundary condition. The mesh sensitivity evaluation comes second. In order to identify and ensure the representation of all the phenomena involved, a computational mesh sensitivity analysis must be done. The mesh sizes were examined from the largest to the smallest in a trial-and-error approach. The model’s precision was improved by fine-tuning the mesh size. A mesh size of 0.007 was used for these simulations. The determination of the correct mesh plays an indispensable role in any numerical model simulation. It is important to minimize the number of cells due to having sufficient resolution as well as this mesh and cell size affect simulation time. The mesh and cell size were tested in a sensitivity analysis, and the 0.007 cell size performed best in terms of results and time spent solving equations.A further decrease in this value in cell size does not affect the accuracy of results for the Q tests. Different cell sizes were used in this study (e.g., 0.06, 0.02, 0.01 and 0.006). These were then reduced to 0.007 to provide more accurate results for the flow tests, as shown in Figure 6. After simulating several different models, the numerical model was set to 40 s, which is a sufficient time to obtain a stable result. The SI system was used to determine the length unit, and degrees Celsius was chosen as the temperature unit. The fluid database was utilized to choose water with a temperature of 20 C and a viscosity of 10 pa s in the model. A review of the previous research indicated that a Renormalized group (RNG) turbulence model was appropriate for the numerical model. The numerical model functioned as a no-slip condition over the ogee’s surface, No specific material characteristics were defined for the ogee spillway, and the roughness was 0.002.

Estimation of Uncertainty in a CFD Application

Computational Fluid Dynamics (CFD) analysis provided a more detailed evaluation of flow characteristics than an experiment. The complexity of the phenomenon makes computer modeling difficult, and its reliability must be evaluated [44]. In this study, error estimates such as an approximate relative error ( e a 21 ) , an extrapolated relative error ( e e x t 21 ) and the fine-grid convergence index ( G C I f i n e 21 ) were used to evaluate the CFD model.
e a 21 = | Q 1 Q 2 Q 1 |
e e x t 21 = | Q e x t 21 Q 1 Q e x t 21 |
G C I f i n e 21 = 1.25 e a 21 r p 21 1
The calculation procedure for three selected grids showed in Table 1. Where Q represents the discharge, r indicates the grid refinement factor, which is greater than 1.3 [45].

2.2. Gene Expression Programming Approaches

The gene expression programming method, which is based on genetic programming and genetic algorithms, was created by Ferreira. GEP is a heuristic search and optimization technique that uses biological evolution to create computer software that can forecast a certain event [24].
One of the GEP’s advantages is that it is multigenic, allowing for the creation of increasingly complex programs with many subprograms. From genetic algorithms, it inherited fixed-length linear chromosomes, as well as expressive parse trees of various sizes and shapes via genetic programming. The G A is a probabilistic search approach that is based on evolution in nature. The following provides the G A ’s general procedure [46]: Step 1: a population is created by randomly selecting n chromosomes (potential problem solutions). Step 2: determine the compatibility of each chromosome ( x ) using the fitness function [ f ( x ) ] . Step 3: The following stages are iterated to create new ones. (i) (Selection): Based on their compatibility, two parent chromosomes are chosen from the population. (ii) (Crossover): Step 1’s parent chromosomes intercross with a particular probability and create two children. (iii) (Mutation): The chromosomes of the created children mutate at a random rate with a given probability. (iv) The children generated as a result of the genetic operators (selection, crossover, and mutation) are added to the new population. Step 4: in the new generation, the parents of the created population are replaced. Step 5: the algorithm terminates and the existing population shows the desired response if favorable circumstances, such as the desired accuracy or the number of iterations stated in the problem, are reached. Step 6: if the algorithm does not end in Step 5, Steps 2–6 are repeated until the desired results are reached. The chromosomes and Expression Trees ( E T s ) are two primary components of GEP. One or more genes representing a mathematical expression may be found on the chromosomes. A gene’s mathematical code can be expressed in two languages: the language of genes and the language of Expression Trees (ET). The GEP genes are divided into two sections: the head and the tail. The head contains mathematical operators, variables, and constants ( + , , * , / , ln , e x p , 1 , a , b , c ) used to represent a mathematical expression. Variables and constants ( 1 , a , b , c ) are included as terminal symbols in the tail. Additional symbols are used if the terminal symbols in the head are insufficient to explain a mathematical equation. Starting at the top line of the tree and reading left to right, top to bottom, the translation of the expression tree is performed. This method’s gene sequences are similar to biological gene sequences [31].
Three input parameters were used in the GEP model (flow head ratio, crest height, crest length). Table 2 shows the parameters of the GEP models. The fitness function used in this study was the root mean square error ( R M S E ) of the training set. The powerful soft computing software package GeneXproTools (Version 5.0.3926) was applied to generate GEP-based discharge prediction models. The program was run for a number of generations before being stopped since the fitness function value did not change. In the GEP model, five options were considered using the various operators listed in Table 3, as well as the output results.

Performance Criteria

A number of evaluation indicators can be used to assess the created models’ forecasting performance using statistical metrics of the goodness of fit [47]. In this study, Performance measures such as root mean square error ( R M S E ) , determination coefficient ( R 2 ) and mean absolute error ( M A E ) were used to evaluate the GEP model.
R M S E = i = 1 N ( Q p Q o ) 2 N R 2 = i = 1 N [ ( Q o Q ¯ o ) ( Q p Q ¯ p ) ] 2 i = 1 N ( Q o Q ¯ o ) 2 ( Q p Q ¯ p ) 2 M A E = 1 N i = 1 N | Q o Q p |
where N is the number of observations, Q o indicates the observed data, Q p indicates the predicted data, Q ¯ o refers to observed data, and Q ¯ p refers to the mean of the predicted data. The mean of the errors is represented by the RMSE, which varies from 0 to , with lower values indicating better model performance. The simulation accuracy of the model is described by R 2 , ranging from 0 to 1. The MAE is a comparison of two continuous variables. The average vertical distance between each point and the y = x line, commonly known as the one-to-one line, is called MAE.

3. Discussion

The numerical model had to be calibrated with experimental data in the first stage. The experimental investigation conducted at the Ministry of Energy’s Water Research Institute, provided discharge values and flow depth that were compared with those of the numerical model evaluated by Flow-3D and GEP models.
To reach steady-state conditions, it is crucial to extract the exact values from the data of a physical or numerical model. In this research, after running numerous different models with the existing numerical model, the acceptable time to extract the results was determined to be 40 s. There is no separation between the flow and the crest bottom, demonstrating the solid boundary condition’s true simulation.
Baffles are two-dimensional surfaces that come in a variety of shapes, including planes, cylinders, cones, and spheres. They come in both porous and non-porous varieties, and they may be used to quantify mass flows, heat streams, and applied forces. The number of elements in the computational mesh depends on the transit flow, which is directly acquired from the section where the baffle is positioned [48]. The flow rate was measured using baffle, and the X-minimum boundary condition was used with a specified pressure. The design head was set to 0.19 m, and the design discharge coefficient and discharge were calculated at 2.15 and 0.0175 m 3 /s, respectively, based on the scale of 1:40. For all data, the physical and numerical models are in good agreement. To make a more objective comparison between models, the results are non-dimensionalized. H e H o and Q Q d values are calculated and shown on graphs in Figure 7 for all models. The experiments conducted in water institute studies were simulated in three dimensions by Flow-3D, as shown in Figure 8.
In order to verify the numerical simulation results the grid size and boundary condition are critical, and it should be remembered that a wrong boundary condition results in a completely different result. Darwin’s theory of evolution is the foundation for gene expression programming [49]. This technique can select the input variables that have the greatest impact on the model automatically. The GEP model formulations were used to extract the most effective indicators. Ferreira [49] provided more information on the technical formulation of the GEP approach.
The GEP is the newest evolutionary algorithm approach, and it is becoming more prominent due to its high accuracy [50]. The main advantage of gene expression programming is that it may express the link between variables explicitly. It should be noted that each run of the GEP method produces a different formula. Accordingly, each model was run multiple times with various GEP parameter settings, and the best models (with the lowest amount of error) were chosen. 70% and 30% of the data were used to train and test discharge estimating models, respectively. The GEP approach demonstratedthat the predicted discharge and experimental discharge were in good agreement. The head over the spillway crest was  H o = h + V a 2 /2 g, in which V a was the approach velocity. The crest height and crest length were utilized as inputs, while discharge was used as an output.
In the GEP model, five options were explored utilizing different operators listed in Table 3, in addition to the ultimate results. The four basic arithmetic operators (+, −,*, / ), were used as well as several basic mathematical functions (x 2 , exp, ln, cube root, Atan, Tanh, min, max ). The GEP’s mathematical equation for option 5 is as follows (9). The parameter d in this equation represents the input parameters, while the parameter c donates constant values determined by GEP. For R 2 = 0.972, RMSE = 0.85, and MAE = 0.64, option 5 was the best alternative. Before 5 was chosen as the best option, a variety of performances were carried out, some of which are listed in Table 4. As well as this, different equations of these runs are shown based on changing the test and train stages. Moreover, the scatter plot in Figure 9 indicates that the training and testing data points are close to the y = x line and, providingthat, GEP can properly predict relative discharge values. The GEP model’s output is typically displayed as mathematical equations, and decision trees Figure 10. The following is the optimal model obtained through GEP modeling for discharge prediction:
Y = 0.0 Y = ( d ( 2 ) ( 1.0 / ( ( ( ( d ( 1 ) + d ( 1 ) ) / G 1 C 1 ) ( d ( 0 ) + d ( 0 ) 2 ) ) ) ) ) Y = Y + a t a n ( ( ( ( ( 1.0 / d ( 0 ) ) ) + ( 1.0 / d ( 0 ) ) ) ) / 2.0 ) + ( ( 1.0 / ( d ( 0 ) ) ) d ( 0 ) ) ) ) Y = Y + ( ( a t a n ( ( a t a n ( G 3 C 6 ) d ( 2 ) ) ) + t a n h ( d ( 0 ) d ( 2 ) ) ) ) / 2.0 ) r e s u l t = Y
where d 0 , d 1 and d 2 correspond to the model input variables and G1C1 and G3C6 are constant values. The numerical constants are shown in Table 5 and the equation derived from the expression trees is shown in Figure 10. The following GEP prediction approach is highlighted in the Materials and Methods section.

4. Conclusions

An ogee spillway is an essential hydraulic structure that can be developed on a variety of sites where concurrent flooding is a concern. Numerous studies based on numerical and analytical relationships have been undertaken to evaluate the discharge values in ogee spillways. The analysis of the behavior and hydraulic parameters of flow overa spillway dam is a difficult and time-consuming task. Flow-3D was utilized for modeling 21 various water levels. In computational fluid dynamics modeling, R N G turbulence model has the highest level of compatibility. Artificial intelligence approaches were also applied to estimate the flow over the spillway. The input (independent) variables in the proposed GEP technique included H e H o , L e H e and H e p , whereas the output (dependent) variable was the discharge of the ogee spillway. R 2 , RMSE and MAE were employed as evaluation indicators. Using these tools aids a designer in establishing the best condition for a hydraulic construction in which the flow pattern is critical. The modeling findings suggest that the CFD methods and GEP models are useful tools for simulating the flow pattern in the ogee spillways. They can be utilized to estimate Q with no need for complicated and time-consuming laboratory techniques. The proposed method is simple to use and accurate enough for practical usage.

Author Contributions

Conceptualization, S.M. and E.A.; data preparation, S.M., E.A. and M.N.; formal analysis, S.M. and E.A.; investigation, E.A. and M.N.; methodology, S.M.; resources, H.M.; supervision, M.N.; validation, S.M., E.A.; visualization, S.M., H.M. and E.A.; writing—original draft, S.M., E.A. and M.N.; writing—review and editing, M.N. and H.M.; correspondence, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GEPGene expression Programming
RMSERoot Mean Square Error
R 2 Coefficient of Determination
MAEMean Absolute Error

Appendix A

Appendix Mathematical equations from GEP model for estimating Q showed in Table A1, performance results of GEP model for option 5 described in Table 4. where d 1 , d 2 and d 3 correspond to the model input variables and numbers related to coefficients are also included in the equations.
Table A1. Each case’s formulation for option 5, described in Table 4.
Table A1. Each case’s formulation for option 5, described in Table 4.
Case NumberEquation
Case 1
Y = 3.87 Y = Y + ( d ( 3 ) / a t a n ( ( ( ( d ( 1 ) 2 ) d ( 2 ) ) + ( 3.4 ) Y = Y + g e p 3 R t ( ( ( ( ( d ( 3 ) + d ( 3 ) ) 2 ) + ( ( 8.61 8.61 ) + d ( 3 ) ) ) / 2.0 ) ) r e s u l t = Y
Case 2
Y = m i n ( ( ( m i n ( ( 3.95 + 9.42 ) , m i n ( d ( 3 ) , 3.95 ) ) + ( 4.6 d ( 2 ) ) ) / 2.0 ) , ( 2.4 ) ) Y = Y + m a x ( a t a n ( e x p ( a t a n ( d ( 1 ) ) ) ) , ( ( ( d ( 3 ) + d ( 2 ) ) + ( d ( 2 ) / d ( 1 ) ) ) / 2.0 ) ) Y = Y + t a n h ( ( ( ( m i n ( 4.18 , d ( 2 ) ) + ( d ( 1 ) d ( 1 ) ) ) / 2.0 ) ( ( d ( 3 ) 3.3 ) + 2.52 ) ) ) r e s u l t = Y
Case 3
Y = ( 1.0 / >( ( ( ( ( ( d ( 1 ) d ( 2 ) ) ( d ( 1 ) 2 ) + r e a l l o g ( 0.51 ) ) / 2.0 ) 2 ) ) ) Y = Y + ( ( ( ( ( 7.2 + d ( 3 ) ) / 2.0 ) + ( d ( 2 ) d ( 2 ) ) ) / 2.0 ) ( ( 7.57 + 2.09 ) ( ( d ( 2 ) + 7.57 ) / 2.0 ) ) ) Y = Y + r e a l l o g ( ( ( ( r e a l l o g ( d ( 2 ) ) + ( 4.157 + d ( 2 ) ) ) + ( ( ( 3.74 + d ( 3 ) ) / 2.0 ) ) / 2.0 ) ) r e s u l t = Y
Case 4
Y = ( 1.0 / ( a t a n ( g e p 3 R t ( ( ( r e a l l o g ( d ( 2 ) + ( 6.03 ) ) + ( ( d ( 2 ) + ( 6.03 ) ) / 2.0 ) ) ) ) ) ) Y = Y + ( ( 1.0 m a x ( r e a l l o g ( d ( 2 ) ) , ( 5.78 ) ) ) ( ( d ( 2 ) d ( 3 ) ) + ( 1.0 / ( 2.28 ) ) ) ) Y = Y + d ( 3 ) r e s u l t = Y
Case 5
Y = ( d ( 3 ) ( ( a t a n ( 2.39 ) ( d ( 1 ) d ( 1 ) ) ) ( 1.0 / ( d ( 1 ) ) ) ) ) Y = Y + t a n h ( ( ( m a x ( 2.86 , d ( 1 ) ) + ( 7.93 + d ( 1 ) ) ) + ( ( d ( 1 ) + d ( 1 ) ) d ( 2 ) ) ) ) Y = Y + t a n h ( ( ( ( ( 6.7 ( 1.05 ) + ( 5.52 ( 7.54 ) ) ) / 2.0 ) ( ( d ( 1 ) 2 ) d ( 3 ) ) ) ) r e s u l t = Y
Case 6
Y = ( ( 1.0 / ( d ( 2 ) ) ) 1.58 ) Y = Y + d ( 3 ) Y = Y + ( 1.0 / ( ( ( ( ( d ( 1 ) d ( 3 ) ) ( 8.1 + ( 1.55 ) ) + ( 1.0 ( ( 3.74 + 1.30 ) / 2.0 ) ) ) r e s u l t = Y
Case 7
Y = g e p 3 R t ( ( ( m i n ( d ( 3 ) , d ( 1 ) ) / t a n h ( d ( 1 ) ) ) + ( m i n ( 2.68 , 2.34 ) d ( 1 ) ) ) ) Y = Y + t a n h ( t a n h ( ( ( ( 3.98 d ( 3 ) ) ( d ( 3 ) ( 7.59 ) ) + ( d ( 1 ) d ( 1 ) ) ) ) ) Y = Y + ( ( ( ( ( d ( 2 ) / d ( 1 ) ) + ( ( 7.26 + d ( 1 ) ) / 2.0 ) ) / 2.0 ) + ( ( d ( 3 ) / 0.97 ) a t a n ( d ( 3 ) ) ) ) / 2.0 ) r e s u l t = Y
Case 8
Y = ( ( g e p 3 R t ( ( 1.0 t a n h ( d ( 2 ) ) ) ) + ( ( d ( 2 ) + r e a l l o g ( d ( 3 ) ) ) / 2.0 ) ) / 2.0 ) Y = Y + ( ( m i n ( g e p 3 R t ( ( d ( 1 ) 6.6 ) ) , ( ( d ( 3 ) + ( 8.98 ) / 2.0 ) ) + d ( 3 ) ) / 2.0 ) Y = Y + ( ( g e p 3 R t ( d ( 3 ) ) / ( d ( 1 ) ) ) + a t a n ( t a n h ( 3.3 ) ) ) r e s u l t = Y
Case 9
Y = ( ( ( ( d ( 3 ) 7.21 ) ( 4.07 d ( 3 ) ) ) + ( ( 1.0 / ( d ( 1 ) ) ) + ( 1.0 / ( d ( 2 ) ) ) ) ) / 2.0 ) Y = Y + a t a n ( ( ( ( t a n h ( 1.7 ) d ( 2 ) ) + d ( 1 ) ) * d ( 1 ) ) ) Y = Y + a t a n ( ( ( 9 t a n h ( 1.706 ) ( d ( 1 ) d ( 3 ) ) ) + a t a n ( ( d ( 3 ) d ( 1 ) ) ) ) ) r e s u l t = Y
Case 10
Y = g e p 3 R t ( 6.92 ) Y = Y + ( 1.0 ( ( ( ( ( d ( 2 ) ( 0.327 ) + r e a l l o g ( d ( 2 ) ) / 2.0 ) ( d ( 3 ) d ( 2 ) ) ) 2 ) ) Y = Y + ( ( d ( 2 ) + d ( 2 ) ) / 2.0 ) r e s u l t = Y

References

  1. Demeke, G.K.; Asfaw, D.H.; Shiferaw, Y.S. 3D hydrodynamic modelling enhances the design of Tendaho Dam Spillway, Ethiopia. Water 2019, 11, 82. [Google Scholar] [CrossRef] [Green Version]
  2. Heravi, G.; Salehi, M.M.; Rostami, M. Identifying cost-optimal options for a typical residential nearly zero energy building’s design in developing countries. Clean Technol. Environ. Policy 2020, 22, 2107–2128. [Google Scholar] [CrossRef]
  3. Korzeniowski, A.; Ghorbani, N. Put Options with Linear Investment for Hull-White Interest Rates. J. Math. Financ. 2021, 11, 152. [Google Scholar] [CrossRef]
  4. Amini, E.; Golbaz, D.; Asadi, R.; Nasiri, M.; Ceylan, O.; Majidi Nezhad, M.; Neshat, M. A Comparative Study of Metaheuristic Algorithms for Wave Energy Converter Power Take-Off Optimisation: A Case Study for Eastern Australia. J. Mar. Sci. Eng. 2021, 9, 490. [Google Scholar] [CrossRef]
  5. Boiten, W. Flow measurement structures. Flow Meas. Instrum. 2002, 13, 203–207. [Google Scholar] [CrossRef]
  6. Tafazzoli, S.; Shafaghat, R.; Alamian, R. Numerical investigation on the multi-body hydrodynamic interactions under Caspian Sea environmental conditions. Ocean. Eng. 2021, 232, 109048. [Google Scholar] [CrossRef]
  7. Mazaheri, P.; Asgarian, B.; Gholami, H. Assessment of strengthening, modification, and repair techniques for aging fixed offshore steel platforms. Appl. Ocean. Res. 2021, 110, 102612. [Google Scholar] [CrossRef]
  8. US Department of the Interior, Bureau of Reclamation. Design of Small Dams; US Department of the Interior: Washington, DC, USA, 1987. [Google Scholar]
  9. Ho, D.; Riddette, K. Application of computational fluid dynamics to evaluate hydraulic performance of spillways in Australia. Aust. J. Civ. Eng. 2010, 6, 81–104. [Google Scholar] [CrossRef]
  10. Tafazzoli, S.; Shafaghat, R.; Alamian, R. Optimization study of a catenary mooring system for a spar floating wind turbine based on its hydrodynamic responses. J. Eng. Marit. Environ. 2021, 235, 657–674. [Google Scholar] [CrossRef]
  11. Choupin, O.; Henriksen, M.; Etemad-Shahidi, A.; Tomlinson, R. Breaking-down and parameterising wave energy converter costs using the CapEx and similitude methods. Energies 2021, 14, 902. [Google Scholar] [CrossRef]
  12. Salazar, F.; San-Mauro, J.; Celigueta, M.Á.; Oñate, E. Shockwaves in spillways with the particle finite element method. Comput. Part. Mech. 2020, 7, 87–99. [Google Scholar] [CrossRef]
  13. Ghorbani, N.; Korzeniowski, A. Adaptive Risk Hedging for Call Options under Cox-Ingersoll- Ross Interest Rates. J. Math. Financ. 2020, 10, 697–704. [Google Scholar] [CrossRef]
  14. Choupin, O.; Andutta, F.P.; Etemad-Shahidi, A.; Tomlinson, R. A decision-making process for wave energy converter and location pairing. Renew. Sustain. Energy Rev. 2021, 147, 111225. [Google Scholar] [CrossRef]
  15. Savage, B.M.; Johnson, M.C. Flow over ogee spillway: Physical and numerical model case study. J. Hydraul. Eng. 2001, 127, 640–649. [Google Scholar] [CrossRef]
  16. Kim, D.G.; Park, J.H. Analysis of flow structure over ogee-spillway in consideration of scale and roughness effects by using CFD model. KSCE J. Civ. Eng. 2005, 9, 161–169. [Google Scholar] [CrossRef]
  17. Peltier, Y.; Dewals, B.; Archambeau, P.; Pirotton, M.; Erpicum, S. Pressure and velocity on an ogee spillway crest operating at high head ratio: Experimental measurements and validation. J. Hydro-Environ. Res. 2018, 19, 128–136. [Google Scholar] [CrossRef]
  18. Jahad, U.; Al-Ameri, R.; Chua, L.; Das, S. Investigating the Effects of Geometry on the Flow Characteristics and Energy Dissipation of Stepped Spillway Using Two-Dimensional Flow Modelling. In Proceedings of the International Association for Hydro-Environment Engineering and Research. Asia Pacific Division Congress, Yogyakarta, Indonesia, 2–5 September 2018; pp. 289–296. [Google Scholar]
  19. Aydin, M.C.; Isik, E.; Ulu, A.E. Numerical modeling of spillway aerators in high-head dams. Appl. Water Sci. 2020, 10, 1–9. [Google Scholar] [CrossRef] [Green Version]
  20. Wan, W.; Liu, B.; Raza, A. Numerical prediction and risk analysis of hydraulic cavitation damage in a High-Speed-Flow spillway. Shock Vib. 2018, 2018, 1817307. [Google Scholar] [CrossRef]
  21. Ghanbari, R.; Heidarnejad, M. Experimental and numerical analysis of flow hydraulics in triangular and rectangular piano key weirs. Water Sci. 2020, 34, 32–38. [Google Scholar] [CrossRef] [Green Version]
  22. Bayon, A.; Valero, D.; García-Bartual, R.; López-Jiménez, P.A. Performance assessment of OpenFOAM and FLOW-3D in the numerical modeling of a low Reynolds number hydraulic jump. Environ. Model. Softw. 2016, 80, 322–335. [Google Scholar] [CrossRef]
  23. Valero, D.; Bung, D.B.; Crookston, B.M. Energy dissipation of a Type III basin under design and adverse conditions for stepped and smooth spillways. J. Hydraul. Eng. 2018, 144, 04018036. [Google Scholar] [CrossRef]
  24. Ferreira, C. Gene expression programming: A new adaptive algorithm for solving problems. arXiv 2001, arXiv:cs/0102027. [Google Scholar]
  25. Yildiz, A.; Yarar, A.; Kumcu, S.Y.; Marti, A.I. Numerical and ANFIS modeling of flow over an ogee-crested spillway. Appl. Water Sci. 2020, 10, 1–10. [Google Scholar] [CrossRef] [Green Version]
  26. Ebtehaj, I.; Bonakdari, H.; Zaji, A.H.; Azimi, H.; Sharifi, A. Gene expression programming to predict the discharge coefficient in rectangular side weirs. Appl. Soft Comput. 2015, 35, 618–628. [Google Scholar] [CrossRef]
  27. Kisi, O.; Emiroglu, M.E.; Bilhan, O.; Guven, A. Prediction of lateral outflow over triangular labyrinth side weirs under subcritical conditions using soft computing approaches. Expert Syst. Appl. 2012, 39, 3454–3460. [Google Scholar] [CrossRef]
  28. Salmasi, F. Effect of downstream apron elevation and downstream submergence in discharge coefficient of ogee weir. ISH J. Hydraul. Eng. 2018, 27, 375–384. [Google Scholar] [CrossRef]
  29. Khan, M.; Azamathulla, H.M.; Tufail, M. Gene-expression programming to predict pier scour depth using laboratory data. J. Hydroinform. 2012, 14, 628–645. [Google Scholar] [CrossRef] [Green Version]
  30. Roushangar, K.; Akhgar, S.; Salmasi, F.; Shiri, J. Modeling energy dissipation over stepped spillways using machine learning approaches. J. Hydrol. 2014, 508, 254–265. [Google Scholar] [CrossRef]
  31. Bagatur, T.; Onen, F. Computation of design coefficients in ogee-crested spillway structure using GEP and regression models. KSCE J. Civ. Eng. 2016, 20, 951–959. [Google Scholar] [CrossRef]
  32. Bertone, E.; Stewart, R.A.; Zhang, H.; Veal, C. Data-driven recursive input–output multivariate statistical forecasting model: Case of DO concentration prediction in Advancetown Lake, Australia. J. Hydroinform. 2015, 17, 817–833. [Google Scholar] [CrossRef]
  33. Esmaeilpour, M.; Ghanbarian, B.; Liang, F.; Liu, H.H. Scale-dependent permeability and formation factor in porous media: Applications of percolation theory. Fuel 2021, 301, 121090. [Google Scholar] [CrossRef]
  34. Azita, M.; Mohammad, I. A Blockchain Information Management Framework for Construction Safety. In Proceedings of the International Conference on Computing in Civil Engineering, Orlando, FL, USA, 12–14 September 2021. [Google Scholar]
  35. Rahimi, Z.; Othman, F.; Shariff, A.R.M. Improving geometric construction of high resolution SAR images using Kriging-based surrogate modelling in mountainous terrain of Malaysia. Int. J. Remote Sens. 2021, 42, 8624–8639. [Google Scholar] [CrossRef]
  36. A GNSS-based weather forecasting approach using Nonlinear Auto Regressive Approach with Exogenous Input (NARX). J. Atmos. Sol. Terr. Phys. 2018, 178, 74–84. [CrossRef]
  37. Mohammadrezapour, O.; Piri, J.; Kisi, O. Comparison of SVM, ANFIS and GEP in modeling monthly potential evapotranspiration in an arid region (Case study: Sistan and Baluchestan Province, Iran). Water Supply 2019, 19, 392–403. [Google Scholar] [CrossRef] [Green Version]
  38. Roushangar, K.; Alami, M.T.; Saghebian, S.M. Modeling open channel flow resistance with dune bedform via heuristic and nonlinear approaches. J. Hydroinform. 2018, 20, 356–375. [Google Scholar] [CrossRef] [Green Version]
  39. Dindarloo, S.R.; Siami-Irdemoosa, E. Estimating the unconfined compressive strength of carbonate rocks using gene expression programming. arXiv 2016, arXiv:1602.03854. [Google Scholar]
  40. Amini, E.; Golbaz, D.; Amini, F.; Majidi Nezhad, M.; Neshat, M.; Astiaso Garcia, D. A parametric study of wave energy converter layouts in real wave models. Energies 2020, 13, 6095. [Google Scholar] [CrossRef]
  41. Neshat, M.; Mirjalili, S.; Sergiienko, N.Y.; Esmaeilzadeh, S.; Amini, E.; Heydari, A.; Garcia, D.A. Layout optimisation of offshore wave energy converters using a novel multi-swarm cooperative algorithm with backtracking strategy: A case study from coasts of Australia. Energy 2022, 239, 122463. [Google Scholar] [CrossRef]
  42. Venkata Raju Badanapuri, Design Principles that are involved in the Design of Flow over an Ogee Crest Spillway. Int. J. Sci. Res. 2019, 8, 245–254.
  43. Sarkardeh, H.; Marosi, M.; Roshan, R. Stepped spillway optimization through numerical and physical modeling. Int. J. Energy Environ. 2015, 6, 597. [Google Scholar]
  44. Matozinhos, C.F.; Campagnole dos Santos, A.A. Two-phase CFD simulation of research reactor siphon breakers: A verification, validation and applicability study. Nucl. Eng. Des. 2018, 326, 7–16. [Google Scholar] [CrossRef]
  45. Celik, I.; Ghia, U.; Roache, P.; Freitas, C.; Coloman, H.; Raad, P. Procedure of Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications. J. Fluids Eng. 2008, 130, 078001. [Google Scholar] [CrossRef] [Green Version]
  46. Azimi, H.; Shabanlou, S.; Ebtehaj, I.; Bonakdari, H.; Kardar, S. Combination of computational fluid dynamics, adaptive neuro-fuzzy inference system, and genetic algorithm for predicting discharge coefficient of rectangular side orifices. J. Irrig. Drain. Eng. 2017, 143, 04017015. [Google Scholar] [CrossRef]
  47. Legates, D.R.; McCabe, G.J., Jr. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  48. Morais, V.H.; Gireli, T.Z.; Vatavuk, P. Numerical and experimental models applied to an ogee crest spillway and roller bucket stilling basin. RBRH 2020, 25, 1–12. [Google Scholar] [CrossRef] [Green Version]
  49. Ferreira, C. Gene Expression Programming: Mathematical Modeling by an Artificial Intelligence; Springer: Berlin/Heidelberg, Germany, 2006; Volume 21. [Google Scholar]
  50. Lu, X.; Wang, X.; Zhang, L.; Zhang, T.; Yang, C.; Song, X.; Yang, Q. Improving forecasting accuracy of river flow using gene expression programming based on wavelet decomposition and de-noising. Hydrol. Res. 2018, 49, 711–723. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The ogee−crested spillway.
Figure 1. The ogee−crested spillway.
Water 14 00650 g001
Figure 2. Physical model of Nazloo Spillway.
Figure 2. Physical model of Nazloo Spillway.
Water 14 00650 g002
Figure 3. Isometric View of an ogee-crested.
Figure 3. Isometric View of an ogee-crested.
Water 14 00650 g003
Figure 4. Depiction of the fractional area volume obstacle (FAVOR) approach that is applied to the solid contours.
Figure 4. Depiction of the fractional area volume obstacle (FAVOR) approach that is applied to the solid contours.
Water 14 00650 g004
Figure 5. Boundary conditions and a three-dimensional structural grid model made up of cuboid cells.
Figure 5. Boundary conditions and a three-dimensional structural grid model made up of cuboid cells.
Water 14 00650 g005
Figure 6. Different resolutions used for sensitivity analysis, (a)—cell size is 0.06, (b)—cell size is 0.02, (c)—cell size is 0.01, and (d)—cell size is 0.006.
Figure 6. Different resolutions used for sensitivity analysis, (a)—cell size is 0.06, (b)—cell size is 0.02, (c)—cell size is 0.01, and (d)—cell size is 0.006.
Water 14 00650 g006
Figure 7. Numerical and physical model comparison.
Figure 7. Numerical and physical model comparison.
Water 14 00650 g007
Figure 8. Numerical modeling for simulating flow depth.
Figure 8. Numerical modeling for simulating flow depth.
Water 14 00650 g008
Figure 9. Comparison between estimated and observed values of Q Q d for Train set and Test set.
Figure 9. Comparison between estimated and observed values of Q Q d for Train set and Test set.
Water 14 00650 g009
Figure 10. Expression tree (ET) for presented model (Equation (8)) and the overall discharge is Q = a+b+c. The first line of an equation is represented by (a,b) stands for the equation’s second line, and the third line of the equation is shown by (c).
Figure 10. Expression tree (ET) for presented model (Equation (8)) and the overall discharge is Q = a+b+c. The first line of an equation is represented by (a,b) stands for the equation’s second line, and the third line of the equation is shown by (c).
Water 14 00650 g010
Table 1. Estimation of discretization error for three cell sizes (0.01, 0.018, 0.02).
Table 1. Estimation of discretization error for three cell sizes (0.01, 0.018, 0.02).
ValuesValuesValues
N 1 , N 2 , N 3 1780000, 305212, 2225001780000, 305212, 222500,1780000, 305212, 222500
r 21 1.721.721.72
r 32 1.561.561.56
Q 1 5.4310.2911.805
Q 2 5.5110.211.69
Q 3 5.3710.0811.58
P2.481.791.2
Q e x t 21 5.3810.3411.93
Q a 21 1.47%0.87%0.97%
Q e x t 21 0.5%0.48%1.04%
Q f i n e 21 0.66%0.67%1.31%
Table 2. Parameters of the GEP Model.
Table 2. Parameters of the GEP Model.
Setting ParametersValue
Functions set+, −, /, x 2 , exp, ln, cube root, Atan, Tanh
Chromosomes30
Head size8
Number of genes3
Linking functionAddition
Fitness function error typeRMSE
Mutation rate0.044
Inversion rate0.1
One-point recombination rate0.3
Two-point recombination rate0.3
Gene recombination rate0.1
Gene transposition rate0.1
Table 3. Different GEP operators with errors of prediction for the Training and Testing datasets.
Table 3. Different GEP operators with errors of prediction for the Training and Testing datasets.
TypeOperatorTrain PhasingTest Phasing
R 2 RMSEMAE R 2 RMSEMAE
Option 1+, −, *, /0.91.571.070.861.731.13
Option 2+, −, *, /, x 2 0.911.531.050.861.731.12
Option 3+, −, *, /, x 2 , exp0.91.651.03 0.86 1.59 1.01
Option 4+, −, *, /, x 2 , exp, ln, cube root0.91.591.060.871.611.06
Option 5+, −, *, /, x 2 , exp, ln, cube root, Atan, Tanh, min, max0.9720.850.640.9121.421.12
Table 4. Different run cases of option 5 (equations are provided in Appendix A, Table A1).
Table 4. Different run cases of option 5 (equations are provided in Appendix A, Table A1).
Case NumberTrain PhasingTest Phasing
R 2 RMSEMAE R 2 RMSEMAE
Case 10.91.6110.950.861.490.97
Case 20.921.430.850.891.711.28
Case 30.960.940.740.861.190.87
Case 40.941.210.950.712.111.3
Case 50.911.480.830.871.891.01
Case 60.921.470.940.851.490.93
Case 70.960.970.70.871.471.07
Case 80.931.320.860.821.731.016
Case 90.911.480.950.891.511.02
Case 100.911.541.020.861.691.06
Table 5. The values of the input and constant parameters used in ET.
Table 5. The values of the input and constant parameters used in ET.
ParametersValueCoefficientsValue
d 0 L e H o G1C12.02
d 1 H e H o G3C60.06
d 2 H e p --
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mozaffari, S.; Amini, E.; Mehdipour, H.; Neshat, M. Flow Discharge Prediction Study Using a CFD-Based Numerical Model and Gene Expression Programming. Water 2022, 14, 650. https://doi.org/10.3390/w14040650

AMA Style

Mozaffari S, Amini E, Mehdipour H, Neshat M. Flow Discharge Prediction Study Using a CFD-Based Numerical Model and Gene Expression Programming. Water. 2022; 14(4):650. https://doi.org/10.3390/w14040650

Chicago/Turabian Style

Mozaffari, Sevda, Erfan Amini, Hossein Mehdipour, and Mehdi Neshat. 2022. "Flow Discharge Prediction Study Using a CFD-Based Numerical Model and Gene Expression Programming" Water 14, no. 4: 650. https://doi.org/10.3390/w14040650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop