Colebrook-White phenomenological equation
Publicado por bernardo (15 intervenciones) el 13/01/2021 21:34:00

Valora esta pregunta


0
f = 0.0390
%% Datos
% Constantes
epsilon = 1.5e-4; % [m]
Dh = 0.1; % [m]
Re = 5000; % [-]
% Ecuaciones como simbólicos
syms f;
eq1 = 1 / sqrt(f);
eq2 = -2 * log10(epsilon/(3.71*Dh) + 2.51/(Re*sqrt(f)));
%% Iteración
i = 1; % Inicialización del contador para el bucle
x = 0.1; % Valor de x inicial para empezar el bucle
error = 0.2; % Error absoluto inicial para entrar en el bucle
while error > 0.1
% Calcular f
f = 1 / x^2;
% Evaluar parte izquierda
eq1_sol = double (subs(eq1, f));
% Evaluar parte derecha
eq2_sol = double (subs(eq2, f));
% Calcular error
error = abs(eq1_sol - eq2_sol);
% Aumentar x
x = x + 0.1;
% Aumentar contador
i = i + 1;
end
%% Resultado
fprintf('f = %.3f (%d iteraciones)\n\n',f,i);
syms f positive
e=0.00015; dh=0.10; re=5000;
g= sqrt(1/f)+2*log10(e/(3.71*dh)+2.51/(re*sqrt(f)));
fplot(g,[0.01 0.1])
raiz=fzero(@(f)sqrt(1/f)+2*log10(e/(3.71*dh)+2.51/(re*sqrt(f))) ,[0.000001 100])
>> raizces_de_una_funcion_14_01_21_foro
raiz =
0.0390