1
Universidad Politécnica Salesiana
Ingeniería Mecatrónica
Proyecto de Métodos Numéricos
Jonathan Quiroz, Pablo Aguirre, Mateo lignia
Resolución de un ejercicio aplicando un método numérico IV. DEFINICIÓN DEL PROBLEMA
a base de programación Lee y Duffy (1976) relacionan el coeficiente de fricción para
el flujo de una suspensión de partículas fibrosas con el numero
Resumen—Con el fin de diseñar un programa para realizar un de Reynolds mediante la siguiente ecuación empírica:
ejercicio utilizamos los diferentes conceptos de métodos numéricos
con apoyo de la programación para crear una script funcional 1 1 5.6
para encontrar la respuesta al problema. = ( ) 𝐿𝑛(𝑅𝑒√𝑓) + (14 − )
√𝑓 𝑘 𝑘
Abstract— In order to design a program to perform an exercise f: coeficiente o factor de fricción.
we use the different concepts of numerical methods with
Re: numero de Reynolds.
programming support to create a functional script to find the
answer to the problem. K: constante que depende de la concentración de la
suspensión.
Palabras clave— método, estructura, script, problema, ¿Cuál es el valor de f, cuando k=0?28 y Re=370?
variables.
V. PARÁMETROS
I. INTRODUCCIÓN PARÁMETROS DE ENTRADA:
E N el presente informe vamos a abordar los principales
aspectos de la creación de un programa para resolver un
ejercicio planteado usando el método de la bisección y un
K: CONSTANTE QUE DEPENDE DE LA CONCENTRACIÓN DE LA
SUSPENSIÓN
programa Matlab para demostrar el funcionamiento, y la
RE: EL NÚMERO DE REYNOLDS
resolución del ejercicio.
PARÁMETROS DE SALIDA:
II. OBJETIVO GENERAL
Crear un programa en base al método de la bisección que F: COEFICIENTE O FACTOR DE FRICCIÓN
facilite la resolución del ejercicio.
VI. CAMPO DEL PROBLEMA
III. OBJETIVOS ESPECÍFICOS
El problema que relaciona el coeficiente de fricción para el
flujo de una suspensión de partículas fibrosas con el número
Demostrar practica y teóricamente cómo funciona el de Reynolds pertenece o hace referencia en el campo de física.
programa
VII. MÉTODO PARA RESOLVER
Diseñar el script por medio de software Matlab. Para la resolución del presente problema matemático, se
escogió el método de la bisección, debido a que este es uno de
Dar a conocer la función del programa. los métodos más sencillos y de fácil intuición para resolver
ecuaciones en una variable, se basa en el teorema del valor
Explicar la razón por la cual se escogió el método medio, el cual establece que toda función continua f en un
numérico. intervalo cerrado [a,b], se elige un intervalo inicial y se va
acortando el intervalo dado hasta llegar a una raíz cercana al
2
punto de corte de dicha función, de esta forma, se asegura la Número de Reynolds:
existencia de al menos una solución de la ecuación f(x)=0. El número de Reynolds (Re) es un número adimensional
utilizado en mecánica de fluidos, diseño de reactores y
fenómenos de transporte para caracterizar el movimiento de
un fluido. Su valor indica si el flujo sigue un modelo laminar o
turbulento.[5]
Bisección
Este método consiste en obtener una mejor aproximación de la
raíz a partir de un intervalo inicial (a,b) en el cual hay un
cambio de signo en la función[2], es decir: f(a)f(b)<0.
Se obtiene el punto medio:
𝑎+𝑏
Ilustración 1:metodo de biseccion(Quintana A.,2005) 𝑥𝑚 =
2
xm es la nueva aproximación a la raíz, y se vuelve a tomar un
VIII. ALCANCE intervalo, pero ahora mas pequeño, considerando que siga
existiendo un cambio de signo en la función[2],, es decir, el
El programa podrá realizar las acciones indicadas: lectura de nuevo intervalo queda determinado por:
variables y aplicación del método de la bisección.
IX. MARCO TEÓRICO
Geoffrey G. Duffy
Geoffrey Duffy es profesor emérito de ingeniería química en
la Universidad de Auckland.[1]
También es propietario y operador de Fiber Dynamics Ltd,
una empresa de consultoría y comercialización especializada
en problemas de secado por infrarrojos y mecánica de fluidos.
[1]
Duffy ha escrito numerosos artículos para Climate Realists ,
un sitio web que promueve las opiniones de los escépticos del Ilustración 2:intervalo de biseccion(Osorio R.,2007)
cambio climático de Nueva Zelanda. [1]
Coeficiente de fricción: El método termina cuando se cumple con alguna condición de
paro, en este programa la condición es la tolerancia:
vincula la oposición al deslizamiento que ofrecen las
superficies de dos cuerpos en contacto según la intensidad del |𝑥(𝑖 + 1) − 𝑥𝑖| ≤ ℇ
apoyo mutuo que experimentan. Es un coeficiente
adimensional. El valor del coeficiente de rozamiento es Este es un método “de encierro”, para aplicarlo se debe contar
característico de cada par de materiales en contacto; no es una con un intervalo inicial, en donde f(a)*f(b) < 0. Este método
propiedad intrínseca de un material. Depende además de requiere de menos pasos en un programa, sin embargo
muchos factores como la temperatura, el acabado de las converge mas lentamente que el de Newton-Raphson[2],.
superficies, la velocidad relativa entre las superficies, etc. La
naturaleza de este tipo de fuerza está ligada a las interacciones Los pasos del método son los siguientes:
de las partículas microscópicas de las dos superficies
implicadas.[3] 1.- Localizar un intervalo que contenga al menos una raíz.
Suspensión de partículas: 2.- Dividir el intervalo en dos partes iguales reteniendo la
mitad en donde f(x) cambia de signo, para conservar al menos
son una serie de diminutos cuerpos sólidos o de gotitas de una raíz.
líquidos dispersos en la atmósfera. Son generadas a partir de
alguna actividad antropogénica o natural.[4] 3.- Repetir el procesó varias veces hasta cumplir con la
tolerancia deseada.
3
si:
f(m) f(b)<0 entonces conservar (m,b) como el sem.
intervalo que contiene al menos una raíz. Diseño de la estructura
A cada paso se le llama “iteración” y reduce el intervalo a la Realizamos el programa utilizando el método de la bisección
mitad. como base para resolver este ejercicio ya que tiene mas
precisión al momento de resolver un ejercicio.
Después de cada iteración el intervalo re reduce a la mitad,
después de n iteraciones, el intervalo original se había Pruebas de funcionamiento
reducido 2n veces, por lo tanto, si el intervalo original es de
tamaño “a” y el criterio de convergencia aplicado al valor En las pruebas realizadas al código del programa se obtuvo
absoluto de la diferencia de dos Xm consecutivas es “ j ”, que este funciona de manera correcta, al cumplir con todos los
entonces se requerían “n” iteraciones donde “n” se calcula con parámetros que son seleccionados y cumple con los
la igualdad de la expresión: requerimientos propuestos para la resolución del ejercicio.
𝑎
≤∈
2𝜋
ln(𝑎)−ln(∈)
de donde: 𝑛 = => iteraciones que se requieren.
ln(2)
X. FUNCIONAMIENTO
El programa disemado para permitir la resolucion del ejercicio
consta del siguiente codigo:
% Metodo de Biseccion
clc
clear all
format short;
xi=input('Introduzca el valor de xi= ');
xu=input('Introduzca el valor de xu= ');
cont=input('Introduzca el numero de
iteraciones= ');
fun=input('Introduzca la funcion
f(x)=','s');
f=inline(fun);
for i=1:cont
xr=(xi+xu)/2;
error = abs(f(xr));
A(i,:)=[i xi xu xr
f(xr) error];
if f(xi)*f(xr)<0
xu=xr;
else
xi=xr; XI. CONCLUSIONES
end
end
Se aplico el metodo de la bisecccion por ser uno de
fprintf('\n \ti \txi \txu \txr \tf(xr) los metodos mas faciles al momento de resolver el
\terror \n') ejercicio.
disp(A)
fprintf('Solucion:\n r=%8.5f\n',xr) Al aplicar el codigo de programacion cumple con los
fprintf('error=%8.5f\n',error) requisitos de la formula de friccion.
%Ecuacion para el ejecicio:
(1/0.28)*(log(3750*sqrt(x)))+(14- Se puede decir que el método de bisección es menos
(5.6/0.28))*(sqrt(x))-1; eficiente que el método de Newton - Raphson, pero
este método es mucho más seguro para garantizar la
convergencia del problema matemático planteado,
4
debido a su número de interacciones, y así se asegura
la existencia de al menos una solución de la ecuación.
XII. BIBLIOGRAFÍA
[1] “Geoffrey Duffy | DeSmog.” [Online]. Available:
[Link]
[Accessed: 20-Nov-2019].
[2] M. C. C. Cunha, Metodos Numericos. 2003.
[3] Carmona, Aníbal Isidoro (2004). «Número de
Reynolds». Aerodinámica y actuaciones del avión.
Thomson Paraninfo.
[4] Chow, Ven Te (1982). Hidráulica de los canales
abiertos.
[5] Crowe, Clayton; Elger, Donald; Williams, Roberson;
Roberson, John (2009). Engineering Fluid Mechanics
[Mecánica de Fluidos Ingeniería] (en inglés) (9.ª
edición). John Wiley & Sons.
