Investigation of the Different Models of Elliptical-Lopac Gate Performance under Submerged Flow Conditions

The effect of the elliptic ratio on the hydraulic performance of the LOPAC gate has been investigated using experimental and Computational Fluid Dynamics CFD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$CFD$$\end{document} simulations (Flow-3D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3D$$\end{document} Hydro). Experiments and simulations were carried out at three different flow discharges, three different submerged ratios, and five different elliptic ratios. In this context, the CFD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$CFD$$\end{document} model was first calibrated and verified using the measured data, and then CFD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$CFD$$\end{document} simulation was performed. The ratio of upstream/downstream flow depth, the flow discharge coefficient, and the energy dissipation through the gate have been calculated and analyzed, and the three-dimensional features of the flow have been described. Based on the results, the elliptical LOPAC gates with circular weirs determine extreme values of the controlled flow parameters. Asymmetric recirculation downstream of the gate is sometimes observed in the model predictions. The use of Coriolis coefficients relevant to the entire cross-section, here deployed just for the circular and the traditional rectangular Lopac gates, has allowed a concise way to report and compare complex flows for any experimental condition.


Introduction
It is becoming more necessary for societies to manage water resources more accurately and equitably stirring researchers to seek structures with better and more accurate performance. This can be achieved by optimizing existing structures or by introducing new ones. Several solutions have been introduced and tested to distribute and regulate flow in rivers, such as weirs, sills, embankments, abutments, diversions, sluice or side gates, etc. It would be difficult, and beyond the purposes of this paper, to summarize the various specific pros and cons of each of them. They interact with morphology, with bioecological issues (e.g., fish transit), they involve costs, maintenance, the possibility of maneuver, and changes in their placement, and they behave differently at high and low discharges. Almost all of these aspects still deserve scientific and technical deepening, which may be nowadays faced also through the use of Computational Fluid Dynamics ( CFD ). In fact, recent advances in computer models, besides the boosting up of computer hardware resources, are gaining more and more acceptance as complementary methods for laboratory and field studies.
Gate is one of the hydraulic structures that are used to distribute and regulate the flow discharge and flow depth. There are multiple types of gates, each one having its advantages and disadvantages, and the specific gate type can be selected according to the application (Weyer 2001). Some extensive laboratory studies have been done to show the pros and cons of them (Negm 1998;Negm and Al-Brahim 2002;Akbari et al. 2019;Cox et al. 2015;Salmasi and Abraham 2021) and numerous experimental and numerical studies have been conducted regarding these issues as well (Long et al. 1991;Ferro 2000;Ma et al. 2001;Raiford andKhan 2013, andGumus et al. 2015;Abhash and Pandey 2021). Long et al. (1991) carried out experiments of submerged jets under sluice gates and compared them in satisfactory way with the numerical results. Ma et al. (2001) also investigated the experimental data of Long et al. (1991) by using a turbulent model in the 2D vertical plane. The deviations they observed were attributed to the three-dimensional nature of the flow. The same conclusions, for similar flow conditions, were inferred by Raiford and Khan (2013), who used two models of flow turbulence, the K − E and the Renormalization group ( RNG ), showing that the RNG model gives better vertical profiles of velocity, though both models estimated a water level above the experimental value.
The LOPAC gate, exploited as a flow distribution and regulation structure, gave acceptable performances in executive projects, matching project expectations, such as in the Middle Rio Grande project in New Mexico (Oad and Kinzli 2006). In the 1980s, Peter Langemann designed the LOPAC gate to assist tail end irrigators with balancing fluctuating water supplies. By combining a LOPAC with a hydraulic actuator, Aqua Systems 2000 Inc. (2013) has provided a flexible and economical solution for controlling water in small-and medium-sized canals. After the introduction of the sharp-crested elliptical weir by Cox et al. (2014), the effects of different elliptic ratios on the discharge rate were investigated and relations for the discharge rate were proposed. Afterward, just a few studies have been done on this type of gate. In 2021, Pilbala et al. experimentally studied the elliptical-LOPAC gate ( EL ) with a given shape factor, finding that the discharge coefficient of the flow through the gate in the elliptical case is higher than in the rectangular case. Shafai Bejestan et al. (2020) measured the discharge coefficient associated to the opening of an elliptical Lopac gate gradually converging to the center of channel. Kheybar et al. (2021) used in a laboratory flume a compound structure obtained by installing the Lopac gate at the end of a sudden contraction, made by a short wall orthogonal to the bank. The effects of the presence of the sudden contraction on the discharge coefficient and on energy dissipation have been measured in a set of different opening angles of the Lopac gate, confirming trends already shown by Pilbala et al. (2021) without the sudden contraction.
In all the previous experiments on the elliptical LOPAC gate, only a fixed elliptic ratio (elliptic ratio equal to 1) has been considered. As it happens for the other hydraulic works causing planimetric deviations of riverine flows, such as abutments (Coleman et al. 2005) or bridge piles (Akhlaghi et al. 2019), geometrical changes of the aspect ratio of the Lopac gates have effects on all flow features (on energy dissipation, in particular).
In this study, the primary objective is to investigate gate performance with different elliptical coefficients, submergence ratios, and discharges on the upstream flow depth, discharge coefficient, and energy dissipation.

Experimental Setup Runs, and Measurements
The experiments were performed in a flume that is 0.8 m wide, 0.8 m high, and 10 m long with a very mild (slightly varying around 0.005) constant longitudinal slope, placed in the hydraulic laboratory of Shahid Chamran University of Ahvaz (Fig. 1). The elliptical ratio of a Lopac gate is defined as the ratio n = H∕W , as shown in Fig. 1b. The LOPAC gate ( Fig. 1) is made of galvanized material with a thickness of 2 mm and is installed at 5 m Fig. 1 a) Plan view of the experimental flume, b) the dimension of experimental EL gate ( n=1), c) the EL gate installed in the flume (downstream view), d) the EL gate installed in the flume (side view), e) mesh generating and boundary conditions of the model distance from the upstream end. The inlet discharge was measured by an ultrasonic flow meter (Farasanj Abzar CO.LTDUF6001-S / N: 940,088) with an error of 0.1 l∕s . At the downstream end of the flume, the flow was controlled by a sluice gate mounting a fullwidth rectangular sharp-edged weir. It permitted to regulate the submerged condition of the gate. A mesh basket was installed at the upstream part of the flume to calm and uniformly distribute the inlet flow. In total, 9 experiments were carried out at three discharges (25, 35, and 45 l∕s ) and three submergence ratios (S = y 2 y 1 ) , that are 0.7, 0.8, 0.9 . Important enough is where y 1 and y 2 were measured. The former is taken just upstream of the gate, along the M 1 profile taking place there. As it is a slowly varying profile, it does not matter to be precise in the position where to take the measure; we chose to stay about 1 m upstream of the gate. As far as y 2 is concerned, it is measured just upstream of the downstream uplifted sluice gate by a point gauge, downstream of the 3D recirculation (shown in Fig. 6) caused by the LOPAC gate. In such a way it is controlled only by the sluice gate opening at a given discharge.
The gate opening angle Ɵ is set constant at 35 degrees ( Fig. 1). In each experiment, after that the discharge is adjusted to the desired value, and the selected submergences were imposed by the sluice gate at the end of the flume. After steady conditions were reached, the measurements were performed using a depth-measuring device at specific points.
Herein we resume the flow variables necessary to describe each run: Q is the discharge, y 1 is the water depth upstream of the gate; y 2 is the water depth downstream of the gate, b g is the opening of the gate (Fig. 1), B is the channel width, Ɵ is the angle of the wings of the gate, g is the acceleration of gravity, is the water-specific mass, is the fluid viscosity, e s is the roughness of the bed, s is the longitudinal slope. In our experimental and numerical runs, some variables are constant, e.g., b g ,B Ɵ s,g,ρ,μ,e s . Some of them are set independently of any others, e.g., Q, y 2 , n . Finally, y 1 is just measured, as a result of the configuration the flow reached at steady state. As far as the targets of the paper are concerned, we are interested to evaluate the discharge coefficient C d of a LOPAC gate and the relative energy dissipation ΔE E 1 of the flow through the gate, in submerged conditions. These two quantities will be defined later.
We end up with the following expected dependence of the two targeted quantities of this paper: where, the Froude number is set as F r = Q∕ � By 2 √ gy 2 � , the Reynolds number as (4) and ΔE = E 1 − E 2 , with E = y + Q 2 2gB 2 y 2 . The value of Re is varying in a small interval, between 31,250 and 56,250, then the roughness does not play an important role. F r varies from 0.077 to 0.184, well below unity. This is why we are focusing more on n , the effects of which will be shown hereafter.
The governing equation for this type of gate is the energy balance (Eq. 3) between the cross sections where y 1 and y 2 are measured. Developing it, and using the continuity equation ( Q is constant in any cross-section), we obtain Eq. 4, expressing the flow discharge, under submerged and steady flow conditions: (1) In Eq. 3, C c is a contraction coefficient. The higher the value of C c (which is expected to be smaller than unity), the smaller the difference between the area of the vena contracta and the cross-section. C d in Eq. 4 is directly linked to the coefficient C c of Eq. 3 through the following Eq. 5.
The experiments are designed to indirectly obtain C d from Eq. 4, and C c from Eq. 5. The dependence of C d on the elliptic shape factor, n is therefore worked out through the dependence of C c on it.
As a final remark, we emphasize that the size of the laboratory flume implies values of the Reynolds number of the flows we worked with are smaller than those relevant to real cases. On one side this represents a limitation to export and exploit directly our results in the design of real structures. Nevertheless, our aim is to describe and measure the effects on flow induced by different gates under different conditions, so we deem that our results can keep a general validity also in the perspective of real applications.

Numerical Approach and Runs
A part of the obtained experimental data was first used for calibration and validation of the Flow-3D Hydro model Flow Science Inc et al. (2018). The validated model then was executed for different gate geometries and flow conditions to get the final data we needed to complete our analysis. For the calibration stage, the elliptical LOPAC gate with an elliptic ratio equal to 1 was considered, and the three discharge values were run in the experiments. By refining the mesh size and by selecting the turbulence model, the numerical model was calibrated so that the results of the simulation were mesh independent and compared well with the measured values of water level within the required accuracy. Subsequently, further four LOPAC gates with elliptic ratios ( n = H∕W ) equal to 0.5 ( n = 0.

Turbulence Closures
In this study, four turbulence models of RNG , K − E , K − , and Large Eddy Simulation ( LES ) have been used to measure their accuracy in estimating the characteristics of flow through the EL gate. The equations of RNG and K − E models are similar; however, equation constants are found empirically in the standard K − E model and explicitly in the RNG model. The RNG model is known to describe low-intensity turbulence and flows having strong shear regions more accurately than the standard K − E model, so the former model has wider applicability. The K − model is deemed to be accurate in particular to describe flow domains near wall boundaries and flows with streamwise pressure gradients, like jets and wakes. In the LES model, the effects of small scale turbulence are represented by an eddy viscosity. In this study, in Section 4.1, all these turbulence models have been exploited and compared.

Mesh Generation and Boundary Conditions
Three mesh blocks (the first mesh block upstream of the gate, the middle mesh block around the LOPAC gate, and the final one downstream of it) (Fig. 1e) have been employed. The mesh dimension is constant throughout any block, and it is the same in the first and third blocks, whereas it is constant but finer in the second block. The mesh sizes, the large and the small ones have been determined based on three criteria: simulation time, error rate, and capabilities of the applied software/hardware. The extension of the first, second, and third blocks, in the x , y , and z directions, are 4 × 0.8 × 0.4, 0.237 × 0.8 × 0.4, and 4.763 × 0.8 × 0.4 m 3 , respectively (Fig. 1e).
The boundary conditions used in the model consisted of imposing the flow discharge at the upstream cross-section and the flow depth at the downstream cross-section. In both boundaries, the pressure has been distributed according to the hydrostatic assumption. Moreover, the common boundaries of the mesh blocks were of the (or transmissive) type; floor boundaries and all the flume walls were selected to be of the wall type; and the free surface was of the symmetry type with a null relative pressure (Fig. 1e).

The Sensitivity Analysis of the Mesh Resolution
A mandatory choice in the numerical analysis concerns the mesh accuracy, or, in other words, the independency of the results on the mesh size. We did some preliminary work devoted to selecting the optimum size of the cells that allowed us to get results that are mesh-independent within an assigned error. In this study, we evaluated 3 different mesh sizes ( M0 , M1 , and M2 ) to assess the mesh size. In all cases, the adopted turbulent model was RNG , which, as described hereafter, was selected among the other ones. In each mesh, the cells in blocks 1 and 3 are the same, whereas, in mesh block number 2, which contains the gate, the cells are finer. The cell size of mesh blocks 1 and 2 for M0 were Δx = Δy = Δz = 0.02695 m and Δx = Δy = Δz = 0.00667 m, for M1 Δx = Δy = Δz = 0.0245 m and Δx = Δy = Δz = 0.0061 m, and finally, for M2, Δx = Δy = Δz = 0.02205 m and Δx = Δy = Δz = 0.00554 m. The results have shown that M1 mesh, compared with M0 , gives results with a relative error greater than 0.1% (0.424), while, comparing M2 and M1 , the relative error, is smaller than 0.1% (0.083). M1 , therefore, has been chosen for the simulation with the total number of cells being less than 500,000 cells.

Choice of the Turbulent Model and Validation
The results of the numerical simulations of the Elliptical-LOPAC gate with an ellipticity ratio n = 1 , in comparison with the experimental data, is generally satisfying for all four turbulence closures. Table 1 shows the error rate for various criteria, RMSE (Root Mean Square Error); MAPE (Mean absolute percentage error); RE (relative error); MAD (Median absolute deviation); and MSE (Mean squared error), applied to flow depths at three specified points (upstream, on the gate, and downstream of the LOPAC-elliptical gate), against the relevant measured data, for the various turbulent models in runs at a discharge of 45 l∕s and three submergence ratios. Results show that among the different turbulence models, the best performance, in comparison with the measured flume data, was obtained by the RNG and K − E models; RNG was chosen for the rest of the runs. To validate the RNG model, we then exploited it also in the simulation of the flow through the Elliptical-LOPAC gate for another discharge rate, which is Q = 25l∕s , and the results were in an acceptable range of errors. Hence, in the end, the RNG model was selected to simulate all the cases, also those for which relevant experimental data are not available.

Free Surface Profiles Along the Flume for the Various Elliptic Ratios ( n ) of the Lopac Gates
One of the criteria for assessing the effects induced on the flow by different gate types is by checking the value of the upstream depth. The higher variations of the upstream depths associated with given variations of the discharge, keeping constant y 2 , the more accurate is the use of the characteristic law (Eq. 4) to calculate the discharge (i.e., the depth measurement error rate is less). We have exploited the CFD results to observe how the elliptic ratios n are affecting the submergence ratios, keeping constant the downstream depth y 2 . Figure 2 illustrates the calculated longitudinal profiles of the flow for different elliptic ratios ( n ). In these Figures, we represent the profile just in a reach from 1 m upstream till 4 m downstream of the gate (placed at x = 0 in the Figures), along the longitudinal symmetric axis; out of this reach the profile is about uniform. The profiles in Fig. 2 indicate that elliptical Lopac gates with elliptic ratios n less than 1 ( n < 1 ) have a better performance because they present higher flow depths ( y 1 ) upstream of the gate (e.g., with Q = 45l∕s , y 1 = 0.37 m at n = 0.5 , y 1 = 0.35 m at n = 1 and y 1 = 0.357 m at n = 1.5).
The dependence of y 1 on n , at given flow conditions (same Q and S ) is monotonously decreasing as n increases in the range n < 1 whereas the opposite takes place for n > 1 . Furthermore, the rate of change of y 1 with n, is greater for n less than 1. Hence, gates with elliptic ratios with n = 1 have the least effects on the increase of the flow head upstream of the gate, whereas gates with elliptic ratios smaller than 1 have shown the best performance in enhancing the flow depth y 2 , upstream of the gate.

Coriolis' Coefficient ( ˛,ˇ)
Equations 6 and 7 report the well-known Coriolis compensations coefficients, , and , that are used to express the flux of momentum and mechanical energy, respectively, in terms of cross-section, and averaged velocity. Hence, their behaviors are a function of the longitudinal position, only, and reflect the variations of velocities and depths across the flow section, lumping in the strong three-dimensional features induced by the gate. In particular, they take into account the asymmetries of the flow observed in the recirculating zones, downstream of the gate, and the flux overtopping the curved, sharp-crested weir of the Lopac gate. It is also reasonable and confirmed by the results, that the higher the Coriolis coefficients are around the gate, the higher the energy dissipation will be. All the results we present are obtained by running the Flow-3D code.
In Fig. 3a, b it is the role of the elliptic coefficient n on and to be focused on. We see that n = 1 presents the lowest values of and throughout the flow domain. To highlight the role of the flow overtopping the weir we compare the EL gate with a traditional rectangular Lopac ( RL ) gate (Pilbala et al. 2021) which conveys the flow between its vertical sides. In this case, the contraction of the vena downstream of the gate is more pronounced than for a Lopac gate. In the comparison, we changed only the shape of the gate, while all the other parameters (i.e., the flow discharge, the width of the gate and the channel, the depth y 2 , etc.) are the same. The results, reported in Figs. 3c, d, indicate that the y 1 depth is higher for the rectangular gate, and so is therefore also the energy dissipation of the flow through the gate. Accordingly, the Coriolis coefficients are higher for the rectangular gate, being the difference between the RL and EL gates bigger when the submergence is higher.

Elliptic Ratio Effect on Energy Dissipation ( E)
The relative flow energy dissipation, ΔE = , is important to characterize the flow behavior around Lopac gates and to direct the gate selection. E 1 is the Energy before the gate ( E 1 = y 1 + 1 V 2 1 2g ) and E 2 is Energy after the gate ( E 2 = y 2 + 2 V 2 2 2g ). As shown in Fig. 4, for a certain type of gate (given n ), the relative energy dissipation decreases with S . This can be explained by observing that by increasing the submergence ratios S , the Froude number of the flow in front of the gate decreases; hence, a smaller amount of kinetic energy to be dissipated through the gate is available. The reference gate ( n = 1 ) allows a minimum energy dissipation in all flow conditions, as displayed in Fig. 4. An explanation of this effect is not straightforward, but quite important as a criterion to select the gate.
We already know (Pilbala et al. 2021) that the circular gate profile ( n = 1 ) allows more friendly management than the rectangular (n → ∞) one. In this study, we see (Fig. 4) that n = 1 corresponds to the minimum relative energy dissipation independently of Q and S . In particular, the minimum relative energy dissipation, obtained with n = 1 , submergence ratio 0.9, and Q = 25l∕s the (minimum) ( ΔE ), was equal to 7.6%. Vice versa, the highest energy loss was around 31%, reached with n = 0.5 , submergence ratio 0.7, and Q = 45l∕s (Fig. 4). Furthermore, we observe that the gates with n = 0.5 and n = 0.7 have always dissipated more relative energy than the cases n = 1.3 and n = 1.5 . The role of Q in the value of relative dissipation is generally less important than that caused by S and n.

Elliptic Ratios Effect on Energy Dissipation Coefficient ( C c ) and Discharge Coefficient ( C d )
This section is aimed at investigating the effect of different elliptic ratios n on contraction coefficient and discharge coefficient, C c and C d , respectively. C c , defined by Eqs. 4 and 5, as above highlighted, is also connected to the amount of energy dissipation. This coefficient is close to 1, with small variations as Q varies (Fig. 5), meaning that the model choice for the local (Borda type) energy dissipation is quite appropriate. According to Fig. 5, the large variation of C c , from the minimum and maximum discharge, is relevant to the gate with n = 0.5 and S = 0.9 . By the three panels of Fig. 5, we can also appreciate that S = 0.9 determines the largest deviations from unity, as Q varies. The discharge coefficient, C d , defined in Eqs. 4 and 5, represents a specific coefficient of the discharge passing through the gate. Equation 5 actually consists of a characteristic formula to evaluate the discharge and has a prominent role in classifying the hydraulic performance of the gate. As mentioned in the previous section, the gate model n = 1 has created the lowest increment of the upstream head, which means that it allows the passage of a flow with the highest discharge compared to EL gates with any other value of n . Figure 5 shows all this and presents C d values that, though depending on C c (Eqs. 4 and 5), have variations with Q , generally smaller than the corresponding variations of the C c values. This is quite good as far as the easiness and accuracy of this evaluation formula are concerned. In particular, with n = 1 , the values of C d are about 0.5, close to 1, and slightly higher than 2.1 for submergence ratios S equal to 0.7, 0.8, and 0.9, respectively (Fig. 5). Also, this figure illustrates that the discharge coefficients obtained with n < 1 are smaller than those obtained with n > 1 , at all S ( S = 0.7, 0.8, 0.9 ). EL with n = 1 , therefore, confirms to be a reference among the other EL , as it was associated with the maximum or minimum value of any parameter involved in the phenomenon.

Three-dimensional Flow Pattern in the Gate Range
The simulations of the elliptical LOPAC gate show that, at all discharges and submergence ratios, in the contraction area after the gate, the lowest velocities are displayed for the case n = 1 , the highest velocities for n < 1 , as reported in Fig. 6a. It also shows the vortices forming downstream of the Lopac gates, at n = 0.5, 1, and 1.5, on both sides of the flume. These vortices are quite similar, but not identical, for the three gates. The stagnant water 0 0 Our results indicate that the stagnant water space area behind the gate with an elliptic ratio 0.5 has the minimum extension. An important aspect, which somehow unexpectedly came out from most of the simulation, is the asymmetric feature. This is not triggered by the geometry, boundary, or initial conditions, which all would favor a symmetric solution. The simulation starts from an arbitrary symmetric initial condition, which is not important for the analysis. Later on, the flow autonomously evolves toward its final configuration, which could be steady or not, and, if not, periodic (Varchanis et al. 2020, Williamson and Brown 1998, Tavelli and Dumbser 2014 or chaotic. From the observation of the movies that present the predicted evolution in time of the flow, we observed that the flow reaches, in all conditions, a steady state, independently of the initial conditions. Furthermore, we see that the flow may deviate, downstream of the gate, toward either side of the channel, whereas it keeps the symmetry upstream of the gate and also inside it. Figure 6 clearly illustrates this issue. Boundary and initial conditions are not playing a role in this outcome, as is proven by the fact that in some runs the symmetric configuration of the flow is maintained. We are deeming this feature is interesting and important to be reported, and we are aware that it could be responsible for some consequences; for instance, the stress field around one bank is increased compared to the symmetric solution. We also wonder whether some sediment transport would create positive or negative feedback, contributing to enhancing or damping the flow deviations. From our runs, we saw that there are conditions that permit to achieve symmetric numerical solutions. One of these is the case with a wide enough opening angle, and another one is the case with a high enough submergence ratio. An increase of the angle (and so of the width at the end of the gate) or an increase of the submergence ratio, favor the maintenance of a symmetric solution, from the start of the simulation till the reaching of its final steady state. Figure 6 illustrate this point. In particular, Fig. 6b, c compare the steady solution for the case of two gate angles, given all the other quantities (discharge, submergence, type of gate, bed slope, channel width, roughness, etc.), while in Fig. 6d, e are the submergence rate to vary.
Given the computational effort required for a single run and the amount of quantities to tune (including here also the turbulent model), it has not been possible to go further in the analysis of the setting allowing or not a symmetric solution. All this acknowledged, it was not even the focus of this work to go in-depth and enlighten the mathematical/numerical origins of the asymmetry.

Conclusion
The Elliptical LOPAC gate, introduced in 2021, can be considered as a new generation of weir-gates, with specific performance features, to be used for the management of irrigation canals. In this study, using laboratory information and CFD code, the elliptic ratio effect of elliptical LOPAC gates on the upstream flow control, discharge coefficient, and energy dissipation of the flow through the gate were evaluated. The following items can be mentioned as key research results: 1. Elliptic ratios n less than 1 ( n = 0.5, n = 0.7 ) yield higher values of the flow depth in the upstream flow than ratios n larger than 1 ( n = 1.3, n = 1.5 ); the lowest upstream flow depth was obtained with the elliptic ratio n = 1. 2. The highest discharge coefficient ( C d =2.1) was obtained for the elliptic ratio 1 and the minimum was C d =0.3 for the elliptic ratio 0.5. The results show that elliptic ratios greater than 1 ( n = 1.3, n = 1.5 ) determine a higher discharge coefficient than ratios n smaller than 1 ( n = 0.5, n = 0.7). 3. The lowest relative energy dissipation rate (i.e., 7.9%) was obtained for the gate with the elliptic ratio n = 1 . The energy dissipation of the EL gates with n less than 1 ( n = 0.5, n = 0.7 ) is higher than the energy dissipation of EL gates with n larger than 1(n = 1.3, n = 1.5); 4. Steady symmetric or non-symmetric recirculating flows, downstream of the Lopac gate, have been predicted by the numerical solutions. An increase of the b g ∕B ratio, or an increase of the submergence ratios, favors the maintenance of a symmetric steady solution. 5. By exploiting a formulation of the Coriolis compensation coefficients for the momentum and the mechanical energy fluxes ( and , respectively) involving the entire crosssection, we found a way to lump into single parameters the complex three-dimensional features of the flow, especially around and immediately downstream of the gate. Compared to rectangular gates, the ellipticity of the shape for the gate with n = 1 allows a reduction of the values of the Coriolis coefficients throughout the flow domain. Also, among elliptical the Lopac gates, n = 1 yields the lowest values of and , which is consistent with the energy dissipation rates above presented and discussed.
Funding Open access funding provided by Università degli Studi di Trento within the CRUI-CARE Agreement.
-We are grateful to the Research Council of Shahid Chamran University of Ahvaz for financial support (GN:SCU.WH98.343).
-This work has been supported by the Italian PRIN 2017 project Enterprising (2017SEB7Z8).

Data Availability
The data and code that support the study are available from the corresponding author upon reasonable request.

Declarations
Competing Interests All authors declare they have no financial interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.