Scribd
Upload
44 views3 pages

Math. Solution of An Implicite Equation Using The Zero Function Method

The document describes a method called the Zero Function method to solve the Colebrook equation for calculating friction factor implicitly. It defines a function that calculates the left and right hand sides of the Colebrook equation and their difference, called the Zero Function. Using bisection, it iteratively calculates the Zero Function over a range from 0.0000001 to 1 to find the value where the Zero Function equals 0, giving the friction factor. The function takes the relative roughness and Reynolds number as inputs and returns the friction factor.

Uploaded by

hasan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
44 views3 pages

Math. Solution of An Implicite Equation Using The Zero Function Method

The document describes a method called the Zero Function method to solve the Colebrook equation for calculating friction factor implicitly. It defines a function that calculates the left and right hand sides of the Colebrook equation and their difference, called the Zero Function. Using bisection, it iteratively calculates the Zero Function over a range from 0.0000001 to 1 to find the value where the Zero Function equals 0, giving the friction factor. The function takes the relative roughness and Reynolds number as inputs and returns the friction factor.

Uploaded by

hasan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
You are on page 1/ 3

Solution of the Colebrook equation using the Zero Function method An examination of the function indicates t

The range a - b is defined with


a=
Function ZeroFunction(Rrel, Re, x) b=
'The function to be examined, the Colebrook equation, is Y
' x ^ (-0.5) = -0.869 * (Log(ed / 3.71 + 2.51 / (Re * (x ^ 0.5))))
' The Left Hand Side and Right Hand Side of the equation are: Ya
LHS = x ^ (-0.5)
RHS = -0.869 * (Log(Rrel / 3.71 + 2.51 / (Re * (x ^ 0.5))))
'The Zero function is defned as
ZeroFunction = LHS - RHS
End Function
0 a
Function FrictionFactorColebrook(Rrel, Re) Ym
a = 0.0000001 Yb
b=1
eStop = 0.0001 LHS =
line1: RHS =
m = (a + b) / 2 ZF =
Ya = ZeroFunction(Rrel, Re, a)
Yb = ZeroFunction(Rrel, Re, b) Application
Ym = ZeroFunction(Rrel, Re, m) f=
Rrel =
If Abs(Ym) < eStop Then Re =
FrictionFactorColebrook = m f=
Exit Function
ElseIf Ya * Ym < 0 Then
a=a
b=m
Else
a=m
b=b
End If
GoTo line1
End Function

' Function to solve implicite equations '


' '
' The function has to be given in the ZeroFunction function, in two parts: The Right Hand Side (RHS) '
' and the Left Hand Side (LHS) of the equation. ´ '
' '
' Bisection method '
' The method is based upon the definition of a function, the ZeroFunction, a function that shall have '
' a value near zero when evaluated at the solution value. In this case, the friction factor "f" '
' have only one zero. Thai is, the function shall cross the zero axe only once. '
' The calculation will stop when the absolute value of the ZeroFunction is small enough: Abs(Ym) < eStop '
' Here, eStop is a parameter defined according the needs of the equation to be solved. '
' In addition to the function, it is required to input a range of the solution values. '
' The range is defined with its two extreme values "a" and "b". '
' A requirement is that the ZeroFunction has a negative value when evaluated with one of this values, '
' and a positive value when evaluated with the other. '
' That means, the ZeroFunction cross the line ZeroFunction = 0 when going from "a" to "b" '
' '
' In this aplication, the Colebrook equation for the Datcy-Weisbach friction factor is been solved '
' Since the frictionfactor cannot be less than zero, a very small value is asigned to one end of the '
' range (a). Also, since the frictionfactor will never be more than 1, a value of one is asigned to the '
' other end of the range. '
' An evaluation has to be made to be sure that the ZeroFunction takes a positive value with one extreme '
' value and a negative when evaluated with the other value. '
' '
' The friction factor is a function of the relative rugosity (Rrel) and the Reynolds number (Re), '
' with Rrel = Rabs / d '
' Rabs : Absoklute rugosity '
n examination of the function indicates that the desired value is in the range 0 < x < 1
he range a - b is defined with
0.0000001 The average value is
1 m= (a + b) / 2

m b x
1

x ^ (-0.5) (Left hand side)


-0.869 * (LOG(Rrel / 3.71 + 2.51 / (Re * (x ^ 0.5)))) (Right hand side)
LHS - RHS (Zero function)

FrictionFactorColebrook(Rrel,Re)
0.0010 (Relative rugosity)
200,000 (Reynolds number)
#VALUE! (Friction factor)

'

You might also like