Int. Renew. Energy Dev. 2025, 14. , 1285-1295 | 1285 Contents list available at CBIORE journal website International Journal of Renewable Energy Development Journal homepage: https://ijred. Research Article Renewable energy from irrigation infrastructure: Experimental insights from a Michell-Banki micro-hydropower prototype in a Colombian irrigation district Juan Gonzalo Ardila Maryn . Jean Carlos Acosta Vargas* . Daniela Andrea Narvyez Cortes Programa de ingenierya agrycola. Facultad de ingenierya. Universidad Surcolombiana. Neiva. Huila. Colombia Abstract. This study presents the experimental characterization of a MichellAeBanki micro-hydropower prototype implemented in the Tynel del Ryo Neiva irrigation district (Colombi. , aimed at promoting distributed generation in off-grid rural areas. The system was designed for a nominal flow of 24 mA/s and a net head of 18 meters, capable of delivering up to 20 kW of electrical power. The turbine-generator performance was evaluated under progressive load and variable flow conditions using direct measurements of voltage, current, power, and rotor speed. A custom-built experimental resistor bank (ERB) was implemented to simulate real load scenarios, enabling the analysis of electrical response as the number of luminaires increased. The results revealed a strong linear correlation between rotor speed and generated voltage (RA = 0. , validating the electromechanical design. However, a saturation trend in power output was observed beyond the tenth luminaire, attributed to reduced rotor speed under load. Polynomial regression models were developed to describe the influence of flow rate on speed, voltage, and power. The cubic models for voltage and power achieved coefficients of determination above 97%, with RMSE values of 5. 41 V and 375. 47 W, respectively. Residual plots confirmed the validity of the models and highlighted the importance of operating close to the nominal flow rate to ensure optimal performance. This work demonstrates the feasibility of using MichellAeBanki turbines for rural electrification through irrigation infrastructure. The methodology and findings provide valuable insights for future implementations, emphasizing the need for hydraulic regulation to maintain system efficiency under variable load conditions. Keywords: Distributed generation. Energy saturation. Experimental correlation. Polynomial modeling. @ The author. Published by CBIORE. This is an open access article under the CC BY-SA license . ttp://creativecommons. org/licenses/by-sa/4. 0/). Received: 11th August 2025. Revised: 28th Sept 2025. Accepted: 7th Oct 2025. Available online: 25th Oct 2025 Introduction In off-grid rural areas, limited access to conventional energy sources represents a significant obstacle to social, educational and productive development (Erdiwansyah et al. , 2. This situation especially affects agricultural communities (Bazzana et , 2020. Chel & Kaushik, 2011. Kaygusuz, 2. In this context, harnessing the hydraulic potential of irrigation canals through Michell-Banki type micro-hydropower plants is a viable, low-cost alternative that can be adapted to local conditions (Halder et al. , 2020. Ibayez et al. , 2. These turbines, due to their robust design and capacity to operate efficiently under variable flow rates, allow the transformation of existing infrastructure into sources of distributed electricity generation (Assefa & Tesfay, 2025. Perez-Rodriguez et al. The implementation of this type of technology not only improves the quality of life of irrigation system users, but also strengthens energy self-management processes (LypezGonzylez et al. , 2. , promotes technological innovation in rural areas (Ibayez et al. , 2020. Reyna et al. , 2. and contributes to meeting global targets such as Sustainable Development Goal (SDG) 7, which aims to ensure affordable and sustainable energy for all (Obaideen et al. , 2022. Wendimu et al. , 2. The development of Michell-Banki type micro hydropower plants has been the subject of multiple studies combining numerical, experimental, and applied design approaches. In the Colombian context. Ramyrez Ramyrez & Cerquera Valderrama . identified a hydropower potential of more than 100 kW in the Tynel del Ryo Neiva irrigation district, which led to the design and evaluation of a 20 kW Michell-Banki turbine. Sotto Capera et al. complemented this work and generated characteristic curves of the prototype using CFD simulations. the international level. Derakhshan & Nourbakhsh . proposed methods to characterize pumps operating as turbines, establishing relationships between speed and efficiency that are reflected in the correlation analyses performed in this paper. Sammartano et al. and Elbatran et al. developed hydraulic and electronic regulation schemes for Banki turbines, relating parameters such as power coefficient (C. and tip speed ratio (TSR), which are related to the power and speed curves presented in this work. Verde et al. , although focused on wind turbines, provided methodologies for estimating electromechanical losses in systems with synchronous generators, like the one used in this prototype, which helps interpret the deviations between theoretical and measured Tarimer & Yuzer . and VasiN et al. addressed * Corresponding author Email: acostajeancarlos74@gmail. com (J. Varga. https://doi. org/10. 61435/ijred. ISSN: 2252-4940/A 2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 the design of permanent magnet generators and the selection of turbines based on operating parameters. Finally. Reyna et al. and Niyonzima & Hendrick . highlighted the applicability of Michell-Banki turbines in rural and educational contexts, validating their use in experimental setups similar to the one developed in this project. This background consolidates the technical and scientific framework that supports the experimental characterization presented, and allows the results obtained to be contextualized in terms of the specialized Although MichellAeBanki turbines have been widely studied for rural electrification (Mrope et al. , 2021. Rahman et al. , 2025. Tesfay et al. , 2. , most existing research focuses on numerical simulations (Ardila Cerquera et al. , 2025. Mendonya et al. , 2025. Mereke et al. , 2. and laboratory-scale prototypes under idealized conditions (Macias Rodas et al. , 2022. Romero-Menco et al. , 2024. Sierra-Moreno et al. , 2. There is a lack of experimental studies that assess the real-world performance of these systems when integrated into existing irrigation infrastructure, particularly under progressive load and variable flow scenarios. This gap limits the development of predictive models and practical design guidelines for field implementation. This study provides an experimental characterization of a fullscale MichellAeBanki micro-hydropower system installed in the Tynel del Ryo Neiva irrigation district. Colombia. This study is novel in its integration of hydraulic and electrical measurements under realistic operating conditions, using polynomial modeling to establish functional relationships between flow rate, rotor speed, voltage, and power. These insights contribute to the optimization of rural microgeneration systems and provide a replicable methodology for similar contexts. Specifically, the study analyzes the influence of turbine speed on the quality of the generated electrical energy, the effect of load increase on the hydraulic-mechanical system, and the correlation between valve opening and no-load speed. Characteristic curves are proposed to understand the interaction between flow, power, and voltage, validating the system design and guiding future implementations in non-interconnected rural areas. Method Turbine Selection Rationale The experimental system was installed in the Tynel del Ryo Neiva irrigation district, in the site known as "La Parrilla", located between the municipalities of Campoalegre and Hobo (Huila. Colombi. , at an altitude of approximately 611 meters Fig 2. Hydraulic turbine selection nomogram. Source: MIEM, above sea level, at coordinates 2A36A13. 24" N and 75A22A46. W, based on the WGS 84 international geodetic system, as shown in Figure 1. The selection of the Michell-Banki turbine for this project was based on a combination of hydraulic suitability, economic feasibility, and operational simplicity, particularly in the context of rural electrification under budget constraints. According to the energy potential assessment conducted by Ramyrez and Cerquera . , the site known as AuLa ParrillaAy in the USOIGUA irrigation district offers a flow rate ranging from 0. 24 mA/s and a net head of 18 meters. However, the project had a demonstrative purpose, and therefore the microhydropower plant was designed to use only 0. 24 m3/s of the canal flow. These conditions fall within the operational range of both Kaplan and Michell-Banki turbines, as indicated by standard nomograms used in small hydropower design (Paish. Singhal and Kumar, 2. , as shown in Figure 2, which presents the standard turbine selection nomogram . ead vs. , highlighting the projectAos operating point and the MichellBanki domain. While Kaplan turbines generally offer higher hydraulic efficiency, their design, manufacturing, and maintenance requirements are more complex and costly. contrast, the Michell-Banki turbine presents a more accessible alternative due to its simpler construction, low maintenance requirements, and adaptability to variable flow conditions. These characteristics make it particularly suitable for rural applications where technical resources and financial capacity are limited (Drinkwaard et al. , 2. The decision to implement a Michell-Banki turbine was further supported by the project's funding limitations and the need for a replicable model that could be adopted in similar rural contexts. Microhydro system design and implementation Fig 1. Location of the micro-hydroelectric power plant. Source: Sotto et al. , 2026. The site offers a net head of 18 meters, and the prototype was designed for a flow rate of 0. 24 mA/s to generate approximately 20 kW of electric power. The system includes a lateral intake, a 12-inch-diameter PVC penstock, and a synchronous generator, all integrated into a reinforced structural frame designed to withstand static and dynamic loads (Sotto et al. , 2. The hydraulic and structural components ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 Figure 3. CAD of the turbine rotor, simulated control volume meshes, and CFD results of the performance study. Source: Sotto, et al. , 2023. Figure 5. Fluid inlet box, hydraulic transmission line, and details of the turbine transition and installation. were developed following both empirical guidelines and national construction codes (NSR-. The civil infrastructure includes a intake box, a buried pipeline with concrete anchors spaced every 12 meters, and a turbine housing mounted above the irrigation canal. The structural design was validated using SAP2000 software, ensuring compliance with seismic, wind, and thermal load combinations. The hydraulic and mechanical design of the turbine was carried out following OLADE . optimize turbine performance, a detailed numerical study was conducted using Computational Fluid Dynamics (CFD). The study evaluated the influence of the injectorAos opening angle on turbine efficiency, using a k-O turbulence model and polyhedral mesh to ensure simulation accuracy (Sotto et al. , 2. The CFD simulations highlighted the importance of blade orientation and injector geometry. The turbine was designed with 30 blades and a diameter ratio of 0. However, the use of second-hand oil pipes for blade construction introduced deviations from ideal curvature, which may explain the discrepancy between theoretical and actual performance. Despite this, the design was deemed feasible, and the turbine was built and tested under real conditions as described in the following sections, validating the numerical predictions and confirming the systemAos suitability for rural applications. Figure 3 illustrates the rotor geometry, the control volume meshing strategy employed for the numerical analysis, and representative results from the performance study. These visualizations, adapted from Sotto et al. , offer valuable insights into the interaction between the working fluid and the turbine blades, supporting the design choices and performance After completing the design phase, the construction and installation of the Michell-Banki micro-hydropower system were carried out, adapting the project to the siteAos specific hydraulic and topographic conditions. The Michell-Banki turbine and its casing were fabricated using structural steel, with critical components such as the rotor and blades manufactured from A36 steel and assembled according to detailed CAD The fabrication process involved CNC plasma cutting, precision welding, and dynamic balancing of the rotor to minimize vibrations during operation (Figure . The intake structure was built using reinforced concrete to ensure reliable water capture, and a 12-inch PVC pipeline was installed underground to convey water from the intake to the Special attention was given to the anchoring of the pipeline, with concrete blocks placed at regular intervals to prevent displacement due to hydraulic forces and terrain slope . ee Figure . The transition between the pipeline and the turbine injector required on-site adjustments to account for dimensional variations, ensuring a watertight and efficient The assembled turbine was mounted on a steel support structure designed to withstand both static and dynamic loads. The alignment of the turbine, transmission system, and generator was carefully verified to optimize power transfer and minimize mechanical losses (Figure . Commissioning of the system included the installation of a custom-built test bench equipped with resistive loads and realtime monitoring sensors to evaluate the electrical output under various load conditions, as described in the following section. General description of the experimental procedure Figure 4. Presentation of turbine components: Rotor Ae Casing Ae Structure, in the metalworking manufacturing workshop. The system includes a lateral intake above the sand trap, a 12-inch-diameter, 60-meter-long PVC pipeline, a manually operated 9-position butterfly valve, a V-belt drive system and a three-phase synchronous generator of the STC type. characterize the electrical performance of the turbine-generator ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 ycE = ycOyce ya Figure 6. Experimental resistance bench (ERB) plate. set, an experimental resistor bank (ERB) was designed and implemented to simulate progressive load conditions through the sequential activation of halogen lamps. The ERB was designed as a tool to evaluate the generator under controlled load conditions, following criteria established in standards such as Ie Std 115-2019 e ISO 8528-1. The bank consists of 20 halogen lamps distributed in four columns: the first two columns contain 12 bulbs of 1500 W each . ,000 W tota. , and the remaining two contain 8 bulbs of 500 W each . ,000 W), reaching a maximum load of 22,000 W, equivalent to 110 % of the generatorAos rated power. The lamps are activated using dust-protected rocker switches. All components are mounted on a 1. 1 y 0. 5 m slanted metal plate, mounted on a modular steel frame with elevated rear legs for improved visibility and operation. The system includes a three-phase 40 A thermomagnetic circuit breaker and AWG 10 type THHN wiring, sized according to NTC 2050. Figure 6 shows the plate where the bench elements are embedded. To measure the rotational speed of the turbine rotor and the driving and driven pulleys, a digital photographic tachometer (TFD) was used, applying reflective tape at the measurement The phase-to-phase and line voltages of the generator were measured with a Uni-T UT890D multimeter. For the calculation of the electrical power, a FLUKE 323 True RMS multimeter was used, recording current and voltage at the ERBAos input terminals. The tests were performed under two main conditions: no load . uminaires of. and progressive load . equential lighting of On each experimental day, the butterfly valve was adjusted from positions 1 to 8, making it possible to evaluate the effect of the flow rate on the turbine speed and electrical In the tests with load, the valve opening was kept fixed . osition 5 or . and the lamps were activated sequentially, recording voltage, current and speed at each step. In no-load tests, the rotor speed was measured as a function of valve opening to establish the hydraulic response curve of the system. A stabilization time of 60 seconds was allowed after each setting before taking the measurements. Data processing and analysis The analysis of the data obtained in the field was performed using a spreadsheet to perform statistical regression. The objective was to characterize the behavior of the turbinegenerator system under progressive load conditions, correlating variables such as the number of activated lamps, rotor speed, line voltage, and generated electrical power. Unlike conventional three-phase systems, the ERB operates with two-phase lamps connected between two generator Therefore, the electrical power was not calculated using the conventional three-phase formula but rather using Equation . where: ycOyce is the phase-to-phase voltage supplying each lamp as it was progressively activated, measured either at the generator terminals or at the input of the totalizer, as previously described. ya is the input current to the ERB on either of the two lines supplying the activated lamp. Each luminaire has a power rating . 0 W or 500 W), but the actual voltage measured in the field varied as they were turned on, which affects actual power consumption, as discussed later. To obtain the total accumulated power, the individual power consumptions of each activated lamp were summed, considering the actual voltage at each step. The results were used to construct a 3D graph relating to the number of luminaires lit, the electrical power consumed by each luminaire and the total electrical power consumed, which was considered equivalent to the generated To estimate the percentage of flow entering the turbine in each butterfly valve opening position, the flow coefficient Kv corresponding to valves of 300 mm nominal diameter was used, according to the manufacturer's technical specifications. The value of Kv represents the water flow . n mA/. through the valve at a pressure drop of 1 bar. The values of Kv were taken for each degree of opening . rom 0A to 90A) and normalized with respect to the maximum value (Kv at 90A), thus obtaining the relative flow rate percentage for each position. This ratio was used as a dimensionless variable in the correlation analyses between flow rate, rotor speed and electrical power generated. The procedure allowed linking the hydraulic behavior of the system with the power curves obtained experimentally, without the need to directly measure the flow rate in the field, which is especially useful in rural systems without flow instrumentation installed. The data were organized in matrices that include the number of luminaires on . , the turbine rotor speed . , the phase-to-phase voltage, the input current to the BRE, and the total cumulative electrical power calculated as described above. In addition, the valve positions at which the required nominal voltage (O 220 V) were achieved, and the most representative sections were selected for modeling. As for statistical modeling, polynomial regressions of degree 1, 2 and 3 were applied to model the relationship between the electrical power generated and three independent variables, namely, the number of luminaires lit, the rotor speed, and the voltage between phases. Each regression was performed using the least squares method, using libraries such as PolynomialFeatures to generate model terms and LinearRegression to fit the coefficients. For the interpretation of the models, each model was considered to have the general form of Equation . yc = yca0 yca1 ycu yca2 ycu 2 yca3 ycu 3 where: yc is the dependent variable whose correlation was ycu is the independent variable . , regime, flow percentag. yca0 is the intercept term, included only in some regressions yca1 , yca2 , yca3 are coefficients indicating the sensitivity of power to changes in ycu. The physical meaning of the coefficients, initially that of the linear coefficient a1 represents the change in power per unit increase in the variable, e. , how many watts are gained per activated lamp. The quadratic coefficient yca2 indicates whether the growth is accelerated or decelerated, a negative value may indicate system saturation. And the cubic coefficient yca3 captures more complex effects such as fluctuations or inflection ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 Progressive accumulation of electrical power according to the number of luminaires 9 11 13 15 17 19 GRADUAL ACCUMULATION CONSUMED POWER [W] Figure 7 shows the behavior of the electrical power generated by the Michell-Banki micro power plant as a function of the number of luminaires turned on in the experimental resistor bank. Each bar represents the power consumed by a specific lamp, which decreased as additional lamps were activated, while the cumulative behavior of total power output under increasing load is better illustrated in Figure 8. This analysis showed that, although the system responds almost linearly in the first sections, after a certain point a saturation trend is observed, where the increase in power per additional luminaire begins to decrease. Sinagra et al. reported that the Michell-Banki turbines maintain good efficiency under varying load conditions, but the observed phenomenon is attributed to the reduction of the generator speed under increasing load, which affected the output voltage and. NUMBER OF LUMINAIRE IN THE ARRAY Fig 7. Individual and cumulative distribution of the electrical power consumed by two-phase luminaires. Current [A] Result and DIscussion consequently, the effective power delivered . iscussed in the next sectio. The cumulative curve is key to identifying the point of maximum operating efficiency of the system, and to validate the test bed design as an energy characterization tool. These results form the basis for the regression models presented in later sections. Figure 8 shows the combined evolution of total electrical power output and line current as the number of luminaires turned on in the experimental resistor bank increases. The blue bars represent the cumulative power calculated for each load step, while the orange bars indicate the current measured at the system input. Initially, both variables show proportional growth, reflecting an efficient response of the generating system to increased demand. However, from the tenth luminaire on, a saturation trend is observed in the power, which begins to stabilize and even decrease slightly, while the current continues to increase more steadily. Popescu & Duinea . also presented efficiency and saturation curves for low-scale systems operating outside the design point. This behavior suggests that the system reaches its limiting operating regime, where the voltage begins to drop due to rotor braking, affecting the effective power delivered, which is discussed in the next The divergence between the two curves from this point onwards indicates that, although the system continues to demand current, the energy conversion efficiency is These results are fundamental to understanding the energy conversion phenomenon and to validate the turbinegenerator energy efficiency models. The results presented in Figures 7 and 8 show that, under progressive load conditions, the electrical power supplied by the Michell-Banki micro-station reached a maximum of 0 kW when the first 10 to 12 lamps were activated, stabilizing or even decreasing slightly in the final This value represents just 20% of the nominal design power . kW), which is consistent with the hydraulic conditions under which the measurements were performed. According to Sotto et al. , the system was designed to operate at 400 rpm with a flow rate of 240 L/s, however, for safety reasons, the tests were conducted with the butterfly valve set to position 5, corresponding to 28% of the maximum available flow rate. Niyonzima & Hendrick . also evaluated the performance of a Michell-Banki turbine under real field conditions, including reduced flow rates. Lamesgin & Ali . also discussed the impact of reduced flow on micro-hydropower plant efficiency. This hydraulic limitation resulted in an average speed of 270 rpm, which was insufficient to reach the rated voltage of 220 V under full load conditions. As a result, the power consumed per luminaire was lower than expected, and the system exhibited signs of energy saturation after the tenth Power [W] points in the system. The coefficient of determination ycI2 indicates what percentage of the variability in the data is explained by the model. If ycI2 OO 1. 0, the fit is considered if ycI2 > 0. 9, very good. and if ycI2 < 0. 7, weak, possibly due to noise or a non-polynomial relationship. In this study, most of the models showed values of ycI2 > 0. 94, which indicates a very high correspondence between the experimental parameters evaluated. In addition to the coefficient of determination . cI2 ), which is already discussed as a measure of model fit, additional statistical analyses were performed to reinforce the validity of the correlation models. For each regression . inear, quadratic, or cubi. , the root mean square error (RMSE) and the residualsAo standard deviation were calculated to quantify the average and typical deviations of the experimental data from the model Residual plots were generated to visually assess the distribution and randomness of the errors, to help identify potential systematic deviations or model limitations. For the linear model, a confidence interval for the slope was computed to assess the statistical significance and robustness of the These statistical tools provide a more comprehensive assessment of model performance, allowing for a transparent discussion of the variability, uncertainty, and reliability of correlations, especially considering the limited number of available experimental data points. 9 11 13 15 17 19 NA Lamps on Power [W] Current [A] Fig 8. Ratio between total power supplied and line current as a function of the number of luminaires lit. ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 present study underscores the mechanical limitations of the MichellAeBanki turbine under unregulated flow conditions. This comparison reinforces the need for integrated hydraulic and electrical control systems to sustain performance under variable operating regimes. Influence of rotor speed on the voltage generated under progressive load. Figure 9 presents the joint evolution of the rotational speed of the generator shaft . and the generated inter-phase voltage, as the lamps in the ERB were activated sequentially. This analysis was performed to further investigate the cause of the saturation observed in electrical power, as evidenced in the previous graphs. The orange bars represent rotor speed, while the blue bars indicate the phase-to-phase voltage supplying each activated lamp, measured at the ERB totalizer input. It is observed that, as the load increases, the speed decreases This decrease in speed results in a drop in the voltage generated, especially from the tenth luminaire on. Guo et al. discussed how variable speed affects both efficiency and generated voltage. The joint behavior of these two variables confirms that electrical power is not limited by generator capacity, but by a reduction in speed that directly affects the electromagnetic induction. This finding validates the need to implement hydraulic regulation systems that maintain the speed within optimal ranges, even under high load conditions. Derived from the behavior shown in Figure 9. Figure 10 presents the initial correlation analysis between turbine shaft speed and phase-to-phase voltage generated by the MichellBanki micro power plant. The data show a highly significant direct linear relationship, with Equation . as the fitted linear ycO = 0. And a coefficient of determination ycI2 = 0. 9953, indicating 53% of the voltage variability is explained by rotor This result confirms that, in this system, the voltage generated depends almost exclusively on the rotational speed of the rotor, which is consistent with Faraday's law of electromagnetic induction. Sammartano et al. Sritram & Suntivarakorn . and Niyonzima & Hendrick . had confirmed that the generated voltage is directly related to rotor speed, especially under varying load conditions, and showed similar correlations in low-power systems for rural The positive slope of the model . suggests that, given the direct proportionality, for every 100-rpm increase, the voltage increases approximately 49. 5V, allowing Voltage [V] Regime . luminaire was turned on. The line current, on the other hand, continued to increase to values close to 28 A, indicating that the system was still demanding energy, but without being able to convert it efficiently into useful power. These results demonstrate the need to operate the system at hydraulic conditions close to design to achieve maximum performance and justify the implementation of flow regulation strategies for future tests. The performance of the MichellAeBanki turbine under progressive load conditions revealed a saturation trend in electrical power output beyond the tenth luminaire, attributed to reduced flow and suboptimal rotational speed. While direct efficiency measurements were not conducted, the observed behavior aligns with findings by Niyonzima & Hendrick . , who reported turbine efficiencies between 60% and 75% at 294 rpm under full valve opening in laboratory conditions. Their study also showed significant efficiency drops below 50% when operating at partial valve openings, consistent with the reduced performance observed in the present field tests at 28% flow. Unlike their controlled setup, the current study provides insights into real-world limitations of flow regulation and generator coupling, emphasizing the need for hydraulic control systems to maintain optimal operating conditions. The experimental results showed a saturation trend in power output under reduced flow conditions, with significant efficiency losses when operating at less than 30% of the design flow. This behavior is consistent with the findings of Sinagra et al. who demonstrated through CFD simulations that efficiency drops when the inlet area is reduced without maintaining optimal fluid velocity. Their study proposed a semicircular segment regulation system to preserve constant hydraulic head and inlet velocity, thereby stabilizing turbine performance under variable flow. While our system lacked such regulation, the observed efficiency decline reinforces the importance of hydraulic control strategies in maintaining energy conversion efficiency in micro-hydro installations. The efficiency behavior of the MichellAeBanki turbine under reduced flow conditions exhibited a saturation trend and reduced energy conversion, particularly below 30% of the design flow. This pattern is consistent with the findings of Popescu & Duinea . , who evaluated a centrifugal pump operating as a turbine (PAT) and reported a peak efficiency of 65%, with stable performance above 60% at higher flow rates. Although the machine types differ, both studies underscore the sensitivity of micro-hydropower systems to flow variability and the importance of operating near the design point to achieve optimal efficiency. The comparison reinforces the need for careful hydraulic matching and flow regulation in rural energy This behavior is consistent with the findings of Lamesgin & Ali . , who demonstrated that suboptimal geometries in Archimedes screw turbines resulted in significant power losses, even at higher flow rates. Their study showed that optimization of internal parameters such as pitch and inner radius increased efficiency from 55. 4% to 92. While the turbine types differ, both studies emphasize the critical role of hydraulic design and flow regulation in optimizing energy conversion in micro-hydropower systems. The comparison reinforces the need for tailored optimization strategies to adapt turbine performance to real operating conditions. The saturation trend observed in electrical power output beyond the tenth luminaire aligns with findings by Guo et al. , who demonstrated that variable-speed microhydropower systems exhibit reduced efficiency when rotational speed decreases due to increased load. Their implemented adaptive MPPT control to mitigate this effect, whereas the 1 2 3 4 5 6 7 8 9 1011121314151617181920 NA Lamps on Regime [RPM] Voltage [V] Fig 9. Relationship between rotor speed and voltage generated as a function of the number of luminaires lit. ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 Voltage [V] Regime . Fig 10. Correlation analysis between rotor speed and voltage generated in the Michell-Banki micro power plant. estimation of system behavior under varying hydraulic This finding is key to understanding why, in the tests performed with reduced flow . alve position . , the voltage did not reach the nominal value of 220V, limiting the effective power delivered. It further supports the need to maintain rotor speed within optimal ranges to ensure efficient power The quality of the linear fit was further assessed by calculating the root mean square error (RMSE) and the residualsAo standard deviation. The RMSE, which quantifies the average magnitude of the prediction error, was found to be 8. This value indicates that, on average, the predicted voltages deviate from the measured values by approximately 8 V, which is considered acceptable given the experimental conditions and inherent system variability. The standard deviation of the residuals was 8. 17 V, reflecting the typical dispersion of the errors around their mean. The similarity between the RMSE and the residual standard deviation suggests that the errors are randomly distributed and that the model does not exhibit significant bias. To further evaluate the adequacy of the linear model, a residual plot was generated (Figure . The residuals, defined as the difference between the measured and predicted voltages, are plotted against rotor speed. The distribution of the residuals appears reasonably symmetric around zero, with no clear trend or curvature, indicating that the linear model is appropriate for the observed data. Although a few points exhibit larger deviations, particularly at lower speeds, the overall pattern supports the assumption of linearity and confirms the absence of systematic bias in the model. Residual [V] Regime . Fig 11. Residual plot: Linear model (Voltage vs. The slope of the linear model, constrained to pass through the origin, was estimated using least squares regression and found to be 0. 495 V/rpm. To assess the precision of this estimate, the slopeAos standard error was calculated as 0. V/rpm, based on residual variance and the distribution of rotor speed as the independent variable. Using a two-tailed StudentAos t-distribution with 19 degrees of freedom and a 95% confidence level, the critical t-value was 2. This resulted in a 95% confidence interval for the slope: . 479, 0. V/rpm, indicating that the relationship is statistically significant and that the slope is well constrained by the data. The narrow width of the interval supports the robustness of the model within the tested range. The strong linear correlation between rotor speed and generated voltage (RA = 0. reported in this study aligns with the electromechanical behavior described by Guo et al. , who emphasized the critical role of rotational speed in determining output power and voltage in micro-hydro systems. Although GuoAos setup included advanced electronic control, both studies confirm that maintaining optimal rotor speed is essential for efficient energy delivery, particularly in rural electrification contexts with limited regulation infrastructure. Variable flow studies Figure 12 presents the behavior of the Michell-Banki turbine rotor speed under no-load conditions, as a function of the percentage of available flow . , estimated from the butterfly valve flow coefficient, as reported by Derakhshan & Nourbakhsh . Elbatran et al. , and Sammartano et . who reported regime versus flow rate curves fitted with polynomial or nonlinear models. The experimental data show an increasing trend of the regime as the flow rate increases, with a quadratic polynomial fit curve defined by Equation . ycycyyco = Oe714. The coefficient of determinationycI2 = 0. 9472 indicates that the model explains 94. 72 % of the observed variability, which validates its use to characterize the hydraulic behavior of the The negative quadratic term suggests the existence of a peak efficiency point in rotor speed, beyond which the increase in flow rate does not translate into a proportional increase in velocity, possibly due to effects such as turbulence, cavitation or friction losses. This model allows the identification of the optimal range of hydraulic operation to maximize rotor speed under no-load conditions, which is key to ensure adequate electrical generation when the system is connected to the resistor bank. The largest residuals occur at 16. 8% and 10. 0% flows rates, where the model overestimates rotor speed, while at 28. 5% and 100% flow, the predictions are closer to the measured values. These residuals reflect the nonlinearity of the system and the influence of hydraulic and mechanical factors not captured by the model. Nevertheless, the quadratic fit provides a reasonable approximation of the system behavior under variable flow conditions and supports the use of this model for estimating performance within the tested range. The RMSE, which quantifies the average magnitude of the prediction error, was 51. 81 rpm. This means that, on average, the predicted rotor speed deviates from the measured values by approximately 52 rpm across the tested flow rates. The standard deviation of the residuals was 55. 0 rpm, indicating the typical dispersion of the errors around their mean. These statistical ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 Voltage [V] Regime . 2 Percentage of flow rate. Fig 12. Unloaded rotor speed behavior as a function of available flow rate. Residual [V] Percentage of flow rate . Fig 13. Residuals of the quadratic fit for rotor speed as a function of flow rate. indicators provide a complementary perspective on the modelAos highlighting both the accuracy and consistency of the predictions within the experimental range. Figure 13 shows the residuals of the quadratic fit as a function of the normalized flow rate. The distribution of the residuals does not reveal a clear systematic pattern, although some larger deviations were observed at intermediate flow This suggests that, while the quadratic model captures the overall trend of the data, there are specific conditions where the experimental variability is not fully explained by the model. Nevertheless, the residuals remain within a reasonable range, supporting the adequacy of the quadratic fit for describing the relationship between flow rate and rotor speed within the tested Figure 14 shows the behavior of the inter-phase voltage generated by the Michell-Banki micro-hydropower plant as a function of the percentage of available flow, under full load conditions . lamps activate. Verde et al. reported voltage versus speed curves for synchronous generators with fitted cubic models. The fitted model is a third-degree polynomial function, defined by Equation . ycO = 206. 25yc Oe 457. The coefficient of determination ycI2 = 0. 9798 indicates that the model explains 97. 98 % of the observed variability, which validates its use to represent the behavior of the system under real operating conditions. The positive cubic term suggests that voltage increases at an accelerated rate at higher flow levels, while the negative quadratic term indicates a transition zone Percentage of flow rate . Fig 14. Voltage curve generated as a function of the percentage of flow rate at full load. where the system efficiency declines, likely due to turbulence, cavitation, or hydraulic losses. This nonlinear behavior reflects the complexity of the interaction between the hydraulic flow and the electromechanical response of the generator, and helps identify the optimal flow rate range to achieve the rated voltage of 220 V. The model also shows that, at reduced flow conditions . t less than 50%), the voltage generated is insufficient to supply nominal loads, reinforcing the need to operate the system near its hydraulic design point. The residuals remain within an acceptable range for most flow rates, although larger discrepancies are observed at the lowest flow settings, likely due to increased hydraulic losses or measurement uncertainty. Overall, the results support the suitability of the cubic model for representing the voltage response of the system as a function of hydraulic input. The RMSE of the residuals was 5. 41 V. The standard deviation of the residuals was 5. 56 V, reflecting the typical dispersion of the errors around their mean. These values confirm the reliability of the cubic model within the tested flow range and support its use for estimating voltage behavior under full load conditions. Figure 15 shows the distribution of the residuals from the cubic voltage model as a function of the normalized flow rate. The red dashed line at zero serves as a reference for perfect model fit. The residuals are generally symmetrically distributed around zero, with no clear systematic trend or curvature, confirming that the cubic model adequately represents the relationship between flow rate and voltage within the tested However, slightly larger deviations are observed at the lowest and intermediate flow rates, which may be attributed to increased hydraulic losses or measurement uncertainties under these conditions. Figure 16 shows the behavior of the electrical power generated by the Michell-Banki micro power plant as a function of the percentage of available flow, under full load conditions . lights o. The fitted model is a third-degree polynomial function, defined by Equation . ycO = Oe750. 89yc3 Oe 6769. 8yc2 12280yc The coefficient of determination ycI2 = 0. 9721 indicates that the model explains 97. 21% of the observed variability, which validates its use to represent the behavior of the system under real operating conditions. The positive cubic term suggests that power increases at an accelerated rate towards the upper end of the flow range, while the negative quadratic term indicates a transition zone where the system loses efficiency, possibly due to hydraulic saturation effects, turbulence or generator ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 Residual [W] Residual [V] . 3 Percentage of flow rate . Fig 15. Residuals of the cubic voltage model as a function of normalized flow rate under full load conditions. This model identifies the optimum flow range to maximize the power output, and shows that, at low flow conditions . elow 30%), the power generated is less than 2 kW, well below the nominal design capacity . kW). These results reinforce the need to operate the system close to its hydraulic design point to ensure efficient energy conversion. The cubic regression between flow rate and electrical power output confirms the nonlinear behavior of the turbine-generator system under full load conditions. While the present study is based on field measurements. Sammartano et al. conducted CFD simulations to optimize the geometry of MichellAeBanki turbines, reporting peak efficiencies above 88% for impellers with 35 blades and an attack angle of 22A. Their results showed that efficiency remains above 80% across a wide range of flow rates, which aligns with the saturation trend observed in our experimental data. However, unlike their idealized conditions, our system exhibited reduced performance at low flow rates due to hydraulic limitations and generator coupling losses. This comparison highlights the importance of integrating geometric optimization with realworld constraints to achieve reliable energy conversion in rural micro-hydro systems. The RMSE was 375. 47 W, indicating that the predicted power values deviate from the measured values by approximately 375 W. The standard deviation of the residuals 47 W, reflecting the typical dispersion of the errors around their mean. These values confirm the reliability of the cubic model within the tested flow range and support its use for Power [W] Percentage of flow rate. Fig 16. Curve of electrical power generated as a function of the percentage of flow rate at full load. Percentage of flow rate . Fig 17. Residuals of the cubic power model as a function of normalized flow rate under full load conditions. estimating the power output of the system under full load Figure 17 illustrates the distribution of residuals from the cubic power model (Equation . as a function of the normalized flow rate. Each blue dot represents the deviation between the measured and predicted power at a given flow rate, while the red dashed line at zero serves as a reference for ideal model fit. The residuals are generally distributed around the zero line, with no evident systematic bias, indicating that the model adequately captures the nonlinear relationship between flow rate and power output. However, greater deviations are observed at low and intermediate flow rates, which may be attributed to hydraulic inefficiencies, generator response limitations, or measurement variability. Overall, the residual plot supports the robustness of the cubic model within the tested range and reinforces its suitability for estimating power generation under full load conditions. The polynomial regression models developed in this study to correlate flow rate with rotor speed, voltage, and electrical power exhibit similar behavior to the models proposed by Guo et al. , who employed regression techniques based on Hill Chart measurements to characterize the efficiency of microhydro turbines. Both approaches reveal nonlinear trends and saturation effects when operating outside the design point, underscoring the importance of hydraulic regulation to maintain optimal energy conversion. While GuoAos work focused on simulation and control strategies, the present study validates these patterns through field-based experimentation under realistic load and flow conditions. Polynomial regression models applied in this study to correlate flow rate, rotor speed, and voltage with power output align with the modeling approach used by Sritram and Suntivarakorn . , who employed Response Surface Methodology (RSM) to optimize turbine blade configurations. Both studies confirm that nonlinear models are effective in capturing the complex interactions between hydraulic and mechanical parameters in micro-hydro systems. They are also consistent with the approach proposed by Derakhshan and Nourbakhsh . , who derived second- and third-order polynomial equations to estimate the complete characteristic curves of centrifugal pumps operating as turbines. Both studies emphasize the importance of modeling off-design performance to support field implementation in small hydro sites. Other studies have addressed the modeling of electromechanical systems using regression techniques, such as the works by Wu et al. Siavash et al. , and Balacco . , which apply polynomial models and artificial neural ISSN: 2252-4940/A2025. The Author. Published by CBIORE J. Marin et al Int. Renew. Energy Dev 2025, 14. , 1285-1295 networks to estimate parameters of induction machines, wind turbines with variable ducts, and pumps operating as turbines. While they share a focus on correlating hydraulic and electrical variables, these studies are mostly based on simulated environments or normalized data, which justifies the need for research like the present one, grounded in real field However, to deepen this work, it would be advisable to incorporate artificial intelligence techniques such as neural networks or evolutionary regression, which could enhance the predictive capacity of the model and explore more complex nonlinear relationships among system variables. Conclusion The experimental characterization of the MichellAeBanki micro-hydropower system installed in a Colombian irrigation district was consistent with the principles of hydraulic-electric The system showed a strong linear correlation between rotor speed and generated voltage (RA = 0. validating the electromechanical design and confirming that voltage output is linearly dependent on rotor speed. Under progressive load conditions, the electrical power output showed a saturation trend beyond the tenth luminaire, reaching a maximum of approximately 4. 0 kW . nly 20% of the nominal design capacit. This limitation was attributed to reduced flow conditions during testing . % of maximu. , which prevented the system from reaching its optimal operating Polynomial regression models were effectively applied to describe the influence of flow rate on rotor speed, voltage, and The cubic models achieved coefficients of determination above 97%, with RMSE values of 5. 41 V for voltage and 375. W for power. Residual plots confirmed the robustness of the models and emphasized the need to operate close to the design flow rate to ensure efficient energy conversion. The study addresses a critical gap in the literature by providing field-based experimental data on the performance of MichellAeBanki turbines integrated into existing irrigation The methodology . ased on direct measurements and progressive load testin. offers a reproducible framework for future rural electrification projects. It is recommended that hydraulic regulation systems be implemented to ensure rotor speed remains within optimal ranges, especially under variable load conditions, to maximize system efficiency and reliability. Acknowledgments The authors express their sincere gratitude to Universidad Surcolombiana for the institutional support and resources provided for the development of this research. Likewise, special recognition is given to the USOIGUA Users Association for their valuable collaboration, willingness and contribution of information, which were fundamental for the realization of this Author Contributions: A. : Conceptualization, methodology, formal analysis, resources, project administration, writingAioriginal draft. supervision, writingAioriginal draft, validation, writingAireview and editing N. writingAireview and editing, project administration. All authors have read and agreed to the published version of the manuscript. Conflicts of Interest: The authors declare no conflict of interest. References