DOI:

https://doi.org/10.14483/22487638.14929

Publicado:

2019-04-01

Número:

Vol. 23 Núm. 60 (2019): Abril - Junio

Sección:

Estudio de caso

Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0

Determinación del coeficiente de resistencia de Darcy-Weisbach utilizando el enfoque Newton-Raphson para android 4.0

Autores/as

Palabras clave:

Android, coeficiente de fricción, Colebrook-White, Darcy-Weisbach, tuberías, JavaScript, Newton-Raphson, métodos numéricos (es).

Palabras clave:

Android, coefficient of friction, Colebrook-White, Darcy-Weisbach, head losses in pipes, Newton-Raphson, numerical methods, JavaScript code (en).

Biografía del autor/a

Edgar Orlando Ladino Moreno, Nueva Granada Military University

Civil Engineer, Topographic Engineer, MsC Civil Engineering. Associated Professor Universidad Militar Nueva Granada

César Augusto García Ubaque, University Distrital Francisco José de Caldas

Civil Engineer, MsC Civil Engineering, PhD in Engineering. Associated Professor Universidad Distrital Francisco José de Caldas

María Camila García Vaca, Chemical Engineering

Chemical Engineer, MsC Chemical Engineering

Cómo citar

APA

Ladino Moreno, E. O., García Ubaque, C. A., & García Vaca, M. C. (2019). Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0. Tecnura, 23(60), 52–58. https://doi.org/10.14483/22487638.14929

ACM

[1]
Ladino Moreno, E.O., García Ubaque, C.A. y García Vaca, M.C. 2019. Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0. Tecnura. 23, 60 (abr. 2019), 52–58. DOI:https://doi.org/10.14483/22487638.14929.

ACS

(1)
Ladino Moreno, E. O.; García Ubaque, C. A.; García Vaca, M. C. Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0. Tecnura 2019, 23, 52-58.

ABNT

LADINO MORENO, E. O.; GARCÍA UBAQUE, C. A.; GARCÍA VACA, M. C. Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0. Tecnura, [S. l.], v. 23, n. 60, p. 52–58, 2019. DOI: 10.14483/22487638.14929. Disponível em: https://revistas.udistrital.edu.co/index.php/Tecnura/article/view/14929. Acesso em: 11 abr. 2021.

Chicago

Ladino Moreno, Edgar Orlando, César Augusto García Ubaque, y María Camila García Vaca. 2019. «Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0». Tecnura 23 (60):52-58. https://doi.org/10.14483/22487638.14929.

Harvard

Ladino Moreno, E. O., García Ubaque, C. A. y García Vaca, M. C. (2019) «Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0», Tecnura, 23(60), pp. 52–58. doi: 10.14483/22487638.14929.

IEEE

[1]
E. O. Ladino Moreno, C. A. García Ubaque, y M. C. García Vaca, «Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0», Tecnura, vol. 23, n.º 60, pp. 52–58, abr. 2019.

MLA

Ladino Moreno, E. O., C. A. García Ubaque, y M. C. García Vaca. «Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0». Tecnura, vol. 23, n.º 60, abril de 2019, pp. 52-58, doi:10.14483/22487638.14929.

Turabian

Ladino Moreno, Edgar Orlando, César Augusto García Ubaque, y María Camila García Vaca. «Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0». Tecnura 23, no. 60 (abril 1, 2019): 52–58. Accedido abril 11, 2021. https://revistas.udistrital.edu.co/index.php/Tecnura/article/view/14929.

Vancouver

1.
Ladino Moreno EO, García Ubaque CA, García Vaca MC. Darcy-Weisbach resistance coefficient determination using Newton-Raphson approach for android 4.0. Tecnura [Internet]. 1 de abril de 2019 [citado 11 de abril de 2021];23(60):52-8. Disponible en: https://revistas.udistrital.edu.co/index.php/Tecnura/article/view/14929

Descargar cita

Visitas

1197

Dimensions


PlumX


Descargas

Los datos de descargas todavía no están disponibles.

Recibido: 15 de agosto de 2018; Aceptado: 14 de febrero de 2019

Abstract

Context:

The art¡cle presents the fr¡ct¡on coeff¡c¡ent calculat¡on for pressure p¡pes, from Newton-Raphson numer¡cal method, under JavaScr¡pt source code for Andro¡d 4.0 or h¡gher.

Methodology:

The code development ¡s based on Colebrook-Wh¡te model for turbulent flows. The app also establ¡shes the fr¡ct¡on coeff¡c¡ent for lam¡nar flows from the Po¡seu¡lle equat¡on and determ¡nes the range for the enter flow ¡n a trans¡t¡on zone (unstable flow). An algor¡thm of success¡ve approx¡mat¡ons was ¡mplemented for non-l¡near equat¡on solut¡on based on the relat¡ve roughness and Reynolds number.

Results:

In the Newton-Raphson method, ¡f the ¡n¡t¡al value moves away from zero or presents d¡fferent ¡nflect¡on po¡nts, the probab¡l¡ty of d¡vergence of the method ¡ncreases s¡gn¡f¡cantly. To overcome th¡s anomaly, a fr¡ct¡on coeff¡c¡ent equal to 0.015 was proposed as ¡n¡t¡al (seed) value; thereby we guarantee the convergence for turbulent flows w¡th an approx¡mat¡on of 0.0000001. The values calculated by the app were cal¡brated w¡th the Excel analys¡s tool (Goal Seek).

Conclusion:

The app ¡s ava¡lable for free ¡n Google Play Store, and ¡s a tool to support educat¡onal processes and the profess¡onal f¡eld.

Keywords:

Andro¡d, coeff¡c¡ent of fr¡ct¡on, Colebrook-Wh¡te, Darcy-We¡sbach, head losses ¡n p¡pes, JavaScr¡pt code, Newton-Raphson, numer¡cal methods.

Resumen

Contexto:

El artículo presenta el cálculo de coef¡c¡ente de fr¡cc¡ón para tubos de pres¡ón, ut¡l¡zando el método numér¡co Newton-Raphson, usando cód¡go JavaScr¡pt para Andro¡d 4.0 o super¡or.

Metodología:

El desarrollo del Cód¡go está basado en el modelo de Colebrook-Wh¡te para flujos turbulentos. La apl¡cac¡ón tamb¡én establece el coeficiente de fricción para flujos laminares utiliando la ecuación de Poiseuille, y determina el rango de flujo de entrada en zona de trans¡c¡ón (flujo ¡nestable). Se ¡mplementó un algor¡tmo de aprox¡mac¡ones suces¡vas para soluc¡ón de ecuac¡ones no lineales basado en la rugosidad relativa y en el número de Reynolds.

Resultados:

En el método de Newton-Raphson, la probabilidad de divergencia aumenta considerabemente si el valor inicial es diferente de cero o presenta distintos puntos de inflección. Para superar esta anomalía, se propone un coeficiente de fricción de 0.015 como valor inicial (valor semilla). De esta manera se garantiza la convergencia de flujos turbulentos con aproximación de 0.0000001. Los valores calculados por la aplicación fueron calibrados utilizando la herramienta de análisis de Excel© (Goal Seek).

Conclusión:

La aplicación está disponible en Google Play Store y es gratuita. Es una herramienta para apoyar procesos educativos y trabajos en sector productivo.

Palabras clave:

Android, coeficiente de fricción, Colebrook-White, Darcy-Weisbach, tuberías, JavaScript, Newton-Raphson, métodos numéricos.

INTRODUCTION

In the last decade, the development of apps for mobile devices on Android platform, along with the increase of processing velocity and storage capacity, have originated user- service oriented technologies. This is especially evident in the case of specialized user needs. An increase of software development for mobile devices to solve specific problems in business (m-business), marketplace (m-commerce), and education (m-learning) has been observed (Gasca, Camargo, & Medina, 2013) (Yin, Weng, & Chu, 2012).

Educative Software Engineering (ESE) classifies these computational tools into algorithms and heuristics. Algorithms pose different systematic processes that lead the student to a determined answer. On the other hand, heuristic tools lead to appropriation of knowledge by means of experimentation and autonomous discovering. Synergies between both algorithmic and heuristic tools produce the intelligent tutorial systems development (Galvis A., 1992).

Therefore, the aim of this paper is to present the development of a mobile app in Android platform that calculate Darcy's friction factor, to be used as a tool in both learning and professional processes. JavaScript source code for Android v. 4.0 or higher was used. Poiseuille model was used for laminar flow cases, and Colebrook-White model for turbulent flow cases. Transition zone between the two types of developed flow regimen, where flow present an unstable behavior, was also calculated.

Friction coefficient

According to regulation for potable water and basic sanitation (RAS 2017), hydraulic calculation for pressured systems in Colombia must be effectuated using Darcy-Weisbach 1845 model. It is a widely known empiric equation that determines load loss due to energy dissipation in the form of friction. This friction is the result of the interaction between fluid and its conducing structure, or pipeline (equation (1)).

Where, h f : load loss (m); f: friction coefficient; L: pipe length (m); D: pipe diameter (m); V: flow velocity (m/s); g gravity acceleration (m/s2).

There are numerous models to calculate Darcy's friction coefficient, or friction factor. Most of them are empirical models limited to the range of experimentation in which they were formulated, such as Moody equation. Another example is Swamee-Jain general equation. It presents satisfactory results for most flow conditions, except for smooth turbulent flow, in which a correction factor is necessary (Andrade, 2001).

However, Colebrook-White model has been widely preferred due to its precision within turbulent conditions. Since it is a non-explicit expression, its solution requires the use of numeric methods. In this work, we used Newton-Rapson method for optimizing point integration (equation (2)).

Where, f: Friction coefficient; e: pipe rugosity (m); D: pipe diameter (m); R e : Reynolds Number.

Laminar Flow

Laminar Flow is determined by values of Reynolds number below 2300. It is known that under this flow condition, viscous forces are significant in comparison to inertial forces, and therefore the fluid moves in an overlapping form with independent layers. In this case, the system behaves according to the Newton equation for viscous fluids. (Saldarriaga, 2007).

In figure 1, we see the typical velocity distribution of a laminar flow inside a cylindric pipe. The highest velocity is reached in the center of the tube, and the lowest flow velocities are found in the surface of the pipe, where the fluid is in direct contact with the solid walls. The result of this phenomenon is a parabolic velocity profile. (Saldarriaga, 2007).

Laminar flow scheme

Figure 1: Laminar flow scheme

Source: own work

In the condition of a laminar flow, and therefore a Reynolds Number below 2300, friction coefficient may be calculated from Poiseuille equation (equation (3)).

According to Reynolds' early experiments, for values between 2300 and 4000, approximately, the flow undergoes a transitional behavior. Because of that, this range of values of Reynolds Number is known as the Transition Zone. In this case, neither laminar or turbulent flow are well defined, and thus friction factor cannot be calculated (equation (4)).

Turbulent Flow

Turbulent flow may be defined as a vectoral chaos; in it we find continuous mixture conditions, due to vortex rupture. (Sotelo, 1997). Despite this widely observed phenomenon, recent studies have found vortex behavioral patterns in turbulent systems, which may defy the paradigms of turbulence according to current notions. In fact, tests to model disorganized and complex turbulent flow conditions into series of somehow organized movement patterns, known as "coherent structures", have been made with mixed results. (Dennis, 2015).

Newton-Raphson Method

Newton-Raphson method optimizes the Newton point integration method by defining a more robust convergence equation. Convergence will exist if the difference of absolute value of calculated values in two successive iterations decreases while the number of iterations (n) increases (Quinta & Villalobos, 2005). The algorithm has been generalized in numerous forms to solve non-linear problems, equation systems, non-linear differential, and integral equations (Díaz & Benítez, 1998). However, if the function is not differentiable and presents discontinuity in the calculation interval, the Newton-Raphson method will not converge. In such cases, it will be necessary to apply another numeric method to find the root of the non-linear equation. Newton-Raphson method requires an initial assumed value f 0 (seed value) (equation (5)).

Now, in the problem of calculating the friction factor, Colebrook-White equation becomes the function, as noted below (equation (6)).

And, the first derivative becomes (equation (7)).

Therefore, Newton-Rapson convergence equation will become (equation (8)).

METHODOLOGY

In this work, we followed four main methodological steps (figure 2): (1) to define the numeric method to solve Colebrook-White equation in the Turbulent Flow condition, (2) to stablish the flow diagram of the app's iterative process, (3) to develop the source code in JavaScript, and (4) to calibrate resulting calculations with Goal Seek function in Excel.

methodological steps followed during the development of the app.

Figure 2: methodological steps followed during the development of the app.

Source: own work.

Source code was developed using JavaScript. Once the solution converges, the result appears in a separated window on the screen. The relative rugosity and Reynold Number, which are both dimensionless.

In figure 3, we present the flow diagram of the iterative process the app follows to calculate Darcy's friction factor.

App Flow Diagram

Figure 3: App Flow Diagram

Source: own work.

RESULTS

Resulting calculations for Darcy's friction factor were calibrated by comparing them to solutions calculated through the Goal Seek function in Excel. In table 1, we present the comparison between various calculation with different input data.

Table 1: App Calibration

Source: own work.

Laminar flow conditions

The relative rugosity of a pipe is 0.0000015. If the Reynolds Number is 1550, the app conduces the calculation of friction coefficient using Poiseuille model for Laminar Flow. Therefore, the result will be f = 0,0412903, as shown in figure 4.

Laminar flow (Friction coefficient app)

Figure 4: Laminar flow (Friction coefficient app)

Source: own work.

Transition zone flow conditions

If the relative rugosity of the pipe is 0.0000045 and the Reynolds Number is 3456, the flow is in the transition zone, so the app will indicate exactly that as shown in figure 5.

Unstable flow (Friction coefficient app)

Figure 5: Unstable flow (Friction coefficient app)

Source: own work.

Turbulent flow conditions

If the relative rugosity of a pipe is 0.0000018 and the Reynolds Number is 845203, the app will use Newton-Raphson method to solve Colebrook-Whi-te equation and calculate the friction coefficient. Seed value is already determined to f= 0.015. the result is f = 0,0120250, as shown in figure 6.

Turbulent Flow (Friction coefficient app)

Figure 6: Turbulent Flow (Friction coefficient app)

Source: own work.

The solution for turbulent flow conditions is coherent with the one obtained by the Goal Seek function, shown in the table 2.

Table 2: Friction Coefficient calculation through Excel goal seek

Source: own work

CONCLUSIONS

The app "App cálculo del coeficiente de fricción" accurately calculates Darcy's friction factor for Laminar and Turbulent Flow Conditions. It also indicates when the input data corresponds to the Transition Zone.

The development of mobile apps for educational purposes (m-learning) is a valuable tool for autonomous learning. Therefore, the app "App cálculo del coeficiente de fricción" contributes to Fluid Mechanics and Hydraulics academical environments. It allows students to experiment with varied input data and explore different initial conditions of the problem.

REFERENCES

Andrade, L. &. (2001). Analysis of the Swamee-Jain's equation for the calculation of friction factor. Revista Brasileira de Engenharia Agrícola e Ambiental. https://doi.org/10.1590/s1415-43662001000300030 [Link]

Dennis, D. J. (2015). Coherent structures in wall-bounded turbulence. Anais da Academia Brasileira de Ciências, 1162. https://doi.org/10.1590/0001-3765201520140622 [Link]

Díaz & Benítez, J. M. (1998). Introducción a los métodos numéricos para la resolución de ecuaciones. Cádiz: Ingrasa Artes Gráficas.

Galvis A., P. (1992). Ingeniería de Software Educativo. Bogotá D. C.: Ediciones Uniandes.

Gasca M. C. et al, M. (2013). Metodología para el desarrollo de aplicaciones móviles. Tecnura, 22.

Gómez C. et al, R. (1998). Ingeniería de software educativo con modelaje orientado por objetos: un medio para desarrollar micromundos interactivos. Informática Educativa Uniandes, 37.

Lindley, C. (2013). Java Script Enlightenment. New York: Simon St. Laurent and Meghan Blanchette.

Quinta & Villalobos, P. E. (2005). Métodos Numéricos con aplicaciones en Excel. Barcelona: Reverté S. A.

Saldarriaga, J. (2007). Hidráulica de Tuberías. Bogotá D. C.: Alfaomega.

Sotelo, G. (1997). Hidráulica General. México D. F.: Limusa.

Lad¡no M., E.O, García U., C.A. y García V., M.C. (2019). Darcy-We¡sbach res¡stance coeff¡c¡ent determ¡nat¡on us¡ng Newton-Raphson approach for andro¡d 4.0. Tecnura, 23(60), 52-58. DOI: https://do¡.org/10.14483/22487638.14929

Artículos más leídos del mismo autor/a

1 2 3 4 5 > >>