UNIVERSIDAD NACIONAL DE SALTA

FACULTAD DE INGENIERÍA
MATEMÁTICA APLICADA

TRABAJO PRÁCTICO INTEGRADOR

2° Cuatrimestre - Año 2024

Apellido y Nombre
Universidad Nacional de Salta
Facultad de Ingeniería
1er Cuatrimestre 2023

MATEMÁTICA APLICADA
TRABAJO PRÁCTICO INTEGRADOR

MODELOS MATEMATICAS PARA SIMULACION DE LA HIDRODINAMICA

El presente trabajo propone el planteamiento y la resolución de un problema de la
hidrodinámica. Esta rama de la física estudia el movimiento de los fluidos, como los líquidos y los
gases. Se aplica en la construcción de presas, canales, acueductos, entre otros.

Para estudiar la hidrodinámica, se pueden considerar diferentes aproximaciones, tales como:

1) Fluidos incompresibles: se desprecian los cambios de densidad en los fluidos.

2) Régimen estacionario: la velocidad del líquido en cualquier punto es independiente del
tiempo.
3) Flujo laminar: las partículas del fluido se desplazan siguiendo su trayectoria, sin cruzarse
entre sí (es decir no viscoso).

En particular, en este trabajo vamos a aplicarla a la simulación del movimiento de un fluido

dentro de una tubería, considerando el problema en régimen estacionario, pérdida de carga y
condiciones de borde conocidas

Modelo Matemático.
Uno de los principios importantes de la hidrodinámica es el principio de Bernoulli. La ecuación
de Bernoulli relaciona la presión, la energía cinética y la potencial.
El principio de Bernoulli se puede deducir considerando una idealización de la realidad del
fenómeno en estudio. Se toma un modelo que considera conservación de la energía en un fluido
laminar sin pérdidas disipativas (sin pérdida de carga), con flujo constante e incompresible. Así,
para una sección de tubería genérica cuya sección cambia a lo largo del tramo (ver Figura 1) se
tiene:

1 2 1
P1  P2  v2  gh2  v12  gh1 (1)
2 2
Donde:

P1, v1 , h1 : Presión, Velocidad y Altura del fluido en el punto 1

P2, v2 , h2 : Presión, Velocidad y Altura del fluido en el punto 2
 : Densidad del fluido
Figura 1. Ecuación de Bernoulli, [Link]

1 2
La Ec. (1) es equivalente a considerar P  v  gh  cte , lo que permite concluir que no se
2
consideran pérdidas de ningún tipo. Se observa en esta ecuación que al disminuir la sección
aumenta la velocidad del fluido y paralelamente se reduce la presión. No obstante, en tuberías
reales existen pérdidas por fricción y porque el flujo no es perfectamente laminar.
Por otro lado, consideremos el caso de una tubería horizontal, de sección constante por la que
circula un fluido cualquiera. Para un elemento de estudio de longitud dx habrá un gradiente de
presión entre ambas caras. Si no existe variación de altura lo que cambia es la velocidad del
fluido.
Centraremos nuestro estudio en la velocidad y vamos a considerar un caso más apropiado en el
cual las pérdidas por fricción y turbulencia existen y están relacionadas a la velocidad (Figura 2).

A constante
dP

v1 v2

dx

Figura 2. Caso de estudio, modelo con fricción y turbulencia.

En este caso tenemos un gradiente de presión debido a las pérdidas por fricción (proporcionales
a la velocidad) y una resistencia cuadrática (debida al fenómeno de turbulencia), las cuales se
pueden expresar como:

dv
 kv  mv 2 (2)
dx
Donde v(x ) es la velocidad del fluido en función de la posición a lo largo de la tubería, k es el
coeficiente de resistencia lineal debido a la fricción viscosa y m es el coeficiente de resistencia
cuadrática debida a la turbulencia. La Ec. (2) es conocida como la ecuación de Riccati y es un caso
particular de la ecuación diferencial ordinaria no lineal general de Bernoulli
dy
 P (x )y  Q (x )y n con n  2 .
dx

f
Además m  donde:
2D

f es el coeficiente de fricción de Darcy-Weisbach que depende del régimen del flujo y la

rugosidad de la tubería y D es el diámetro de la tubería
 es la densidad del fluido

La fórmula empírica de Colebrook-White dada por

 
1  2.51 
 2 log  D  
f 
 3.7 Re f 


Permite estimar el valor del coeficiente de Darcy-Weisbach. En el caso de un flujo turbulento el

número de Reynolds es Re  4500 y  la rugosidad absoluta del material de la cañería.

CONSIGNA:
Hasta ahora hemos usado el método de Euler explícito (o forward) para resolver estos tipos de
problemas, pero existen otras variantes como es el método de Euler implícito (backward) que,
dependiendo del problema, puede resultar o no en un problema no lineal.

Euler Explícito Euler Implícito

n 1 n 1
y  y n  h (x n , y n ) y  y n  h (x n 1, y n 1 )
Tabla 1: Formulas de los Métodos de Euler

En función de lo anterior se pide que se desarrolle lo siguiente:

a) Identificar de qué tipo de problema se trata (Clasificarlo)
b) Obtener la distribución de velocidad a lo largo de una tubería mediante el
modelo matemático descripto en la Ec. (2) mediante un método analítico.
 c) Resolver el mismo problema a partir de las dos técnicas numéricas descriptas
en la Tabla 1. Implementar dichas técnicas en un programa elaborado en
C++.
d) Realizar un estudio de convergencia de los métodos numéricos.
e) Obtener una tabla de comparación entre los resultados obtenidos con el
método analítico y con ambos métodos numéricos.

DATOS:
GRUPO MATERIAL  k m Longitud Diámetro Velocidad
[1/m] [s/m2] de la de la Inicial
Tubería Tubería
G1 Acero al 0.001 12 m 0.1 m
Carbono
G2 Acero 0.001 10 mv 0.1 m
Inoxidable
G3 Cobre 0.001 13 m 0.12 m
G4 PVC 0.001 15 m 0.12 m
(Policloruro
de Vinilo)
G5 CPVC 0.001 8m 0.1 m
(Cloruro de
Polivinilo
Clorado)
G7 Polietileno de 0.001 9m 0.1 m
Alta Densidad
(PEAD)
G8 Hierro 0.001 11 m 0.12 m
Fundido

Los demás datos que requiera el problema deben ser asumidos por el grupo, ya sea
a partir de consideraciones prácticas o de búsqueda bibliográfica, los cuales deben
ser justificados y las fuentes citadas.

ENTREGABLES:
1) Informe escrito (Monografía) que incluya el desarrollo de cada técnica
utilizada, donde se detallen los criterios y fundamentos de las elecciones
realizadas en los métodos de resolución.
 Además, de los resultados y las simulaciones pedidas se debe incorporar el
código fuente elabroado por el grupo con el algoritmo que le permitió obtener
la respuesta numérica del problema.

2) Una presentación oral que debe responder a la rúbrica de evaluación que se

anexa a este documento

Formato: El establecido por la cátedra (ver plantilla)

Fecha de Consulta y Avance del Informe: Lunes 09/12/24 de 10 a 12 hs
Fecha de Presentación Informe: Miércoles 11/12/24 – Hasta las 12 hs (AV)
Fecha de defensa oral: Viernes 13/12/24 – 10 hs
Defensa: Oral (ver criterios de evaluación)
Forma de Evaluación: Rubrica (ver abajo ANEXOS)

ANEXO I:

Plantilla. Use el siguiente formato para armar su presentación. Si desea incorporar

algún recurso que no está preestablecido en el siguiente formato lo puede hacer
mientras que siga estas pautas generales.

FORMATO
___________________________________________________________________

Títulos (Calibri Light 16, estilo Título MA2023)

Subtítulo (Calibri Light 13, estilo Subtítulo MA2023)
Párrafo (Calibri Cuerpo 13 justificada con 1,5 de espacio, estilo Párrafo MA20223). Utilice
este estilo para los párrafos. Es decir, solamente deberá concentrarse en redactar y luego
en Inicio➡Estilos seleccione el correspondiente a los párrafos.

Imágenes y tablas. Van centradas y con rótulo al pie de la imagen/tabla.

Bibliografía (Calibri Cuerpo 11 espaciado mínimo, estilo Bibliografía MA2023). No olvide citar las fuentes
consultadas para la elaboración de este trabajo. Use las normas APA.

___________________________________________________________________
Criterios de Evaluación. La producción escrita y la defensa oral deben dar cuenta
que usted ha adquirido un lenguaje específico de Matemática Aplicada. Llame a
cada concepto de manera apropiada y clara. Debe demostrar que domina los temas
tratados, es decir se apropió de los saberes disciplinares de Matemática Aplicada.
Usted deberá ser capaz de mostrar su integración de conceptos, durante el cursado
no hubo temas aislados, sino que están relacionados. Exprese las relaciones entre
los conceptos que identifica, esto dará cuenta que no solamente entendió, sino
que, además ganó una herramienta para resolver problemas. Sea criterioso en su
forma de expresión, no necesita inundar de palabras escritas, la capacidad de
síntesis está relacionada con el dominio del tema.

ANEXO II:

Instrumento de Evaluación: Rubrica de Evaluación