Image processing technique for hydraulic
application
1 J. Tellez-Alvarez <Jackson.david.teUez@upc.edu> 1M. Gomez <manuel.gomez@upc.edu> 2 B. Russo <brusso@unizar.es> 1 Institute FLUMEN, Department of Civil and Environmental Engineering Technical University of Catalonia BarcelonaTech, Barcelona, Spain. 2 Group of Hydraulics and Environmental Engineering Technical College of La Almunia (University of Zaragoza) Zaragoza, Spain.
Abstract. In this paper, we present techniques Surface Flow Image Velocimetry (SFIV) that is a practical extension of Particle Image Velocimetry method (PIV). In particular, this technique is to evaluate the behavior on the surface flow for complex flow. SFIV allows to measure complex surface velocity fields for engineering applications. The objective of this work was the application of the SFIV technique to determine the velocity fields and their hydraulic efficiency of grated inlet, making a comparison between the programs able to correlate images like Digiflow or PIVlab program.
Keywords: Digiflow; PIVlab; image processing; grated inlet; MATLAB; Particle Image Processing (PIV); Surface Flow Image Velocimetry (SFIV).
DOI: 10.15514/ISPRAS-2017-29(6)-18
For citation: Tellez-Alvarez J., Gomez M., Russo B. Image processing technique for hydraulic application. ISP RAN/Proc. ISP RAS, vol. 29, issue 6, 2017. pp. 289-298 (in English). DOI: 10.15514/ISPRAS-2017-29(6)-18
1. Introduction
In general, the particle image velocimetry (PIV) method is a classic velocity meter through video cameras, where velocity fields can be obtained by following the definition of velocity u:
u(x, t) = (i)
At
where Ax is the displacement of a marker, located at x at time t, over a short time interval At separating observations of the marker images [1].
Surface Flow Image Velocimetry (SFIV) similar to the Particle Image Velocimetry (PIV) or laser anemometry, that may be used without the need of adding tracer particles to the flow, but following the perturbations of the surface flow as shapes of the vortices or wave flows as defined in [5], for these applications in surface flows local capillary waves, vortices or wakes as in the case of convective or shock induced flows, the vortex cores or the vortex filaments are needed to be able to follow in detail the lagrangian aspects of the flow.
Very high-resolution cameras with high frame speeds are necessary, together with large computer capacity and the implementation of programs coupled to advanced image processing for fluid mechanics, such as Imacal, DigImage and DigiFlow [4, 8].
Technique called Surface Flow Image Velocimetry (SFIV) was developed to obtain the velocity field from the image processing of the surface flow in the vicinity of the grated inlet [6, 7].
The equipment for obtaining these images consists of a camera of high resolution and speed accompanied to a large computer where the image analysis was carried out through the application of the advanced image-processing program for fluid mechanics called Digiflow, developed by the Department of Applied Mathematics and Physical Theory, University of Cambridge and Dalziel Research Partners [2]. In addition, Digiflow is an open software developed for the community of research focus in the area of fluids mechanics and turbulence applications. In order to follow the same line of research and compare the results with Digiflow, we also used a free code for image processing called PIVlab developed by the University of Groningen in the Netherlands. [10, 11], written in the mathematical language MATLAB. Considering that MATLAB language is a private code, but the company Matworks provide a special permission as open software for research community and student, for this reason, now a day is the more common program used in the universities for mathematics application and programming.
2. Experimental Set-Up
The experiment consists of fixed laboratory facilities such as grated inlet platform, 2 reflectors and a grate of study. Mobile installations as the high-resolution camera and speed connected to a large computer.
The platform has dimensions of 5.5 m long and 4 m wide, with a working area of 1.5 m long and 3 meters wide in the vicinity of the grated inlet. The platform is able to modify its longitudinal slope between 0% to 10% and the transverse slope of 0% to 4% [5]. In addition, the maximum approach flow for the experiments is 200 l/s, as shown in Figure1.
In order to obtain high-resolution images of the order of 1280 x 1024 pixels with a capture rate of 150 frames per second, it was decided to place the artificial lighting with two reflectors of power of 500 and 1000 Watt located strategically to provide the best possible illumination in the area of study (Figure 1). In addition, we painted
the gray color the platform, which allowed us to see better the surface flow and avoid reflections. Finally, marks were made around the grate used as reference in the image.
Fig. 1. Area of study
The platform has dimensions of 5.5 m long and 4 m wide, with a working area of 1.5 m long and 3 meters wide in the vicinity of the grated inlet. The platform is able to modify its longitudinal slope between 0% to 10% and the transverse slope of 0% to 4% [5]. In addition, the maximum approach flow for the experiments is 200 l/s, as shown in fig. 2.
In order to obtain high-resolution images of the order of 1280 x 1024 pixels with a capture rate of 150 frames per second, it was decided to place the artificial lighting with two reflectors of power of 500 and 1000 Watt located strategically to provide the best possible illumination in the area of study (fig. 1). In addition, we painted the gray color the platform, which allowed us to see better the surface flow and
291
avoid reflections. Finally, marks were made around the grate used as reference in the image.
3. Surface Flow Image Velocimetry
The Surface Flow Image Velocimetry (SFIV) technique allows the measurement of complex surface flows to obtain the field of velocities in engineering. This methodology has been developed using the physical model of grated inlet located in the Hydraulic Laboratory of the Unicersitat Politécnica of Catalonia BarcelonaTech. The Digiflow algorithm is a method of the synthetic schlieren type, which is a novel technique that gives us a qualitative visualization of the fluctuations of the densities of the frames, and has origin in the classical schlieren method and the moire fringe technique. It also allows us to get good visualization in large study domains [4]. The synthetic schlieren method has different approaches to measure flows within the Digiflow program, such as line refractometry, dot-tracking refractometry, and pattern matching refractometry [3]. The final method is based on PIV techniques and was the method of correlation and image analysis proposed in the SFIV technique with Digiflow.
Fig. 2. Visualization of histogram. Longitudinal slope 0%. Transversal slope 0%. Flowrate 200 l/s
PIVlab is a free access PIV program developed in the MATLAB calculation platform, which allows us to have easier access to this tool. We must take into account that for the image processing using PIVlab, we used the same videos that were used for the analysis with the program Digiflow, only that, in this case, we must take into account that PIVlab does not accept video formats only sequence of fixed images, for this reason, we developed an additional code in MATLAB capable
of decomposing videos into frames, which allowed us to introduce the sequence of images in the PIVlab.
In MATLAB it is possible to generate the histograms of the digital videos, allowing to represent the frequency of intensity of each pixel and the form of the histogram. This is possible through a function within the tools developed in MATLAB, in this case the utility is "imhist".
Fig. 3. Visualization of histogram. Longitudinal slope 10%. Transversal slope 4%. Flowrate 200 l/s.
The intensity scale is in the range of 0 to 256 grayscale. The objective of the visualization of the histograms was to see the difference in form and values for different geometric combinations and a same flow of approximation. In fig. 2 and fig. 3, a comparison of the histograms was made when the platform is completely flat and with a longitudinal slope of 10% and transverse slope of 4%, and these graphs show the significant difference between the intensity values for different videos, which indicates that when we do not have slopes the intensity histogram is more widespread, while in extreme conditions of slope the intensities are accumulating around the mean value of 100. It is evident that this affects the recognition of flow patterns, being one of the reasons why the calculation algorithm must be powerful enough to be able to reproduce the velocity fields for different geometric combinations and approximation flowrates on the platform.
4. Results
After performing several tests with different PIV algorithms, it was finally possible to obtain the velocity fields for 24 different geometries and 5 approximation flows, in order to compare with the Digiflow code. The visualization of the video captured, and the flow approach around the grated inlet is how in the fig. 4. The velocity
fields in the vicinity of the grated inlet are shown in fig. 5 with PIVlab code and fig. 6 using Digiflow program.
Moreover, both algorithm show a good approximation of field velocities around the grated inlet, but the different is on the grated inlet, because Digiflow present better approximation in this area, but also it is not enough, because the holes of the grate represent a behavior more tri-dimensional.
Fig. 4. Visualization of flow approach. Flowrate 200 l/s. Longitudinal slope 10%. Transversal slope 4%.
Fig. 5. Field velocity. Flowrate 200 l/s. Longitudinal slope 10%. Transversal slope 4%. PIVlab code
4. Discussion
Considering that both analysis with different programs have different algorithms, but follow a procedure similar to particle velocimetry (PIV), in both cases the SFIV technique was applied for surface flow, without adding particles as tracers and validating the results with experimental data for the measurement of flowrate collected by the grated inlet. The proposed methodology can become a useful tool to understand the hydraulic behavior of the flow in the surroundings of a catchment of grated inlet where traditional measurement equipment has serious problems and limitations.
Fig. 6. Field velocity. Longitudinal slope 10%. Trasnversal slope 4%. Digiflow program.
As for the difference of the calculation programs, we have noticed that the Digiflow presents us with more advanced options for the analysis and processing of images, since it is a specialized program for visualization of complex flows in fluids mechanics, and it allows us complete stream of videos. In contrast, PIVlab, it has the disadvantage that it is only possible to work with sequence of frames so we need to create a tool capable of decomposing the videos into frames, but at the same time, it is a public access tool written in a common calculation language like is the MATLAB.
4. Conclusion
The Surface Flow Image Velocimetry (SFIV) technique is very useful for determining the fields of surface velocities and the structure of the flow through perturbations of the fluids in surfaces like forms or waves; we must say that works better in a supercritical regime. The proposed methodology of the SFIV technique is
able to obtain velocity fields in the vicinity of the grated inlet from high resolution images.
The results show that the velocity range is between 0.5 m / s and 2.45 m / s, for the 24 geometric combinations studied, with an approximate flow rate of 100 l/s, so we can say that the velocity values calculated through of the Manning and Izzard equation compared to the surface velocities obtained through Digiflow and PIVlab are quite similar.
The SFIV technique, through high resolution video capture and speed, presents certain advantages against traditional instruments of mechanical measures such as Acoustic Doppler Velocimetry (ADV) or EMG Electromagnetic that are able to capture the speed in a single point, while through images it is possible to have a field velocity in the area of focus of the camera, which allows us to better study the behavior of the fluid and thus obtain larger measurement points.
References
[1]. Adrian R.J. (1991) Particle-Imaging Techniques for Experimental Fluid Mechanics. Annual Review of Fluid Mechanics, 23(1), pp.261-304.
[2]. Dalziel R. P. (2012). DigiFlow User Guide. Cambridge, England: Department of Applied Mathematics and Theoretical Physics (DAMTP) - University of Cambridge.
[3]. Dalziel S.B., Hughes G.O. & Sutherland, B.R. (1998) Synthetic Schlieren. Paper 062. Proceedings of the 8th international Symposium on Flow Visualization, In: Carlomagno; Grant (ed), Tennessee, United States.
[4]. Dalziel, S.D. and Redondo J.M. (2007) New visualization and self-similar analysis in experimental turbulence studies. Models, Experiments and Computation in Turbulence. CIMNE, Barcelona, Spain.
[5]. Gómez M., 2008. Curso de Hidrología Urbana. Segunda Ed. Universität Politécnica de Cataluña. Barcelona, España.
[6]. Tellez J., Gómez M., Russo M., Redondo J.M. 2016a. "Metodología para la obtención del campo de velocidad en la proximidad de las rejas de alcantarillado". XXVII Congreso latinoamericano de hidráulica. 28 al 30 de septiembre, 2016. Lima, Perú.
[7]. Tellez J., Gómez M., Russo M., Redondo J.M. 2016b. "Surface Flow Image Velocimetry (SFIV) for Hydraulics Applications" In 18th International Symposium on the Application of Laser and Imaging Techniques to Fluid Mechanics. Julio 4-7. Lisboa, Portugal.
[8]. Tellez J., Russo B., Gómez M. 2015. "Técnica para la obtención del campo de velocidad del flujo superficial en proximidad de rejas de alcantarillado". IV Jornadas de Ingeniería del Agua. Octubre 21-22. Córdoba, España.
[9]. Thielicke W., Stamhuis E. J. (2014). PIVlab - Time-Resolved Digital Particle Image Velocimetry Tool for MATLAB (version: 1.41).
[10]. Thielicke W. & Stamhuis E.J. (2014b). PIVlab - Towards User-friendly, Affordable and Accurate Digital Particle Image Velocimetry in MATLAB. Journal of Open Research Software. 2(1), p.e30. DOI: http://doi.org/10.5334/jors.bl
Технология обработки изображений для гидравлических приложений
1 Дж. Теллез-Альварез <Jackson.david.tellez@upc.edu> 1 М. Гомез <manuel.gomez@upc.edu> 2 Б. Руссо <brusso@unizar.es> 1 Испания, Барселона, Политехнический университет Каталонии, Научно-исследовательский институт FLUMEN, Департамент строительства и охраны окружающей среды 2 Испания, Сарагоса, Университет Сарагосы Технический колледж Ла-Альмунии Группа гидравлики и охраны окружающей среды
Аннотация. В этой статье мы представляем методы Surface Flow Image Velocimetry (SFIV), являющиеся практическим расширением метода Velocimetry Image Particle Image (PIV). В частности, этот метод позволяет оценивать поведение поверхностного потока для комплексного течения. SFIVпозволяет измерять сложные поля поверхностной скорости для инженерных приложений. Целью этой работы было применение метода SFIV для определения полей скоростей и гидравлической эффективности решетчатого водоприемника. Сравнивались результаты программ Digiflow и PIVlab, способных коррелировать изображения.
Ключевые слова: Digiflow; PIVlab; обработка изображений; MATLAB; Particle Image Processing (PIV); Surface Flow Image Velocimetry (SFIV)
DOI: 10.15514/ISPRAS-2017-29(6)-18
For citation: Теллез-Альварез Дж., Гомез М., Руссо Б. Технология обработки изображений для гидравлических приложений. Труды ИСП РАН, том 29, вып. 6, 2017 г., стр. 289-298 (на английском языке). DOI: 10.15514/ISPRAS-2017-29(6)-18
Список литературы
[1]. Adrian R.J. (1991) Particle-Imaging Techniques for Experimental Fluid Mechanics. Annual Review of Fluid Mechanics, 23(1), pp.261-304.
[2]. Dalziel R. P. (2012). DigiFlow User Guide. Cambridge, England: Department of Applied Mathematics and Theoretical Physics (DAMTP) - University of Cambridge.
[3]. Dalziel S.B., Hughes G.O. & Sutherland, B.R. (1998) Synthetic Schlieren. Paper 062. Proceedings of the 8th international Symposium on Flow Visualization, In: Carlomagno; Grant (ed), Tennessee, United States.
[4]. Dalziel, S.D. and Redondo J.M. (2007) New visualization and self-similar analysis in experimental turbulence studies. Models, Experiments and Computation in Turbulence. CIMNE, Barcelona, Spain.
[5]. Gómez M., 2008. Curso de Hidrología Urbana. Segunda Ed. Universitat Politécnica de Cataluña. Barcelona, España.
[6]. Tellez J., Gómez M., Russo M., Redondo J.M. 20^a. "Metodología para la obtención del campo de velocidad en la proximidad de las rejas de alcantarillado". XXVII Congreso latinoamericano de hidráulica. 28 al 30 de septiembre, 201б. Lima, Perú.
[7]. Tellez J., Gómez M., Russo M., Redondo J.M. 20^b. "Surface Flow Image Velocimetry (SFIV) for Hydraulics Applications" In 18th International Symposium on the Application of Laser and Imaging Techniques to Fluid Mechanics. Julio 4-7. Lisboa, Portugal.
[8]. Tellez J., Russo B., Gómez M. 2015. "Técnica para la obtención del campo de velocidad del flujo superficial en proximidad de rejas de alcantarillado". IV Jornadas de Ingeniería del Agua. Octubre 21-22. Córdoba, España.
[9]. Thielicke W., Stamhuis E. J. (2014). PIVlab - Time-Resolved Digital Particle Image Velocimetry Tool for MATLAB (version: 1.41).
[10]. Thielicke W. & Stamhuis E.J. (2014b). PIVlab - Towards User-friendly, Affordable and Accurate Digital Particle Image Velocimetry in MATLAB. Journal of Open Research Software. 2(1), p.e30. DOI: http://doi.org/10.3334/jors.bl