“Año de la recuperación y consolidación de la economía peruana”
UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
(Universidad del Perú, Decana de América)
FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA
Proyecto final de ciclo
Tema: Solución numérica simplificada de las ecuaciones de navier. stokes
para flujo laminar en un canal
Curso
Dinámica de fluidos
DOCENTE:
Felipe Arancel,Llerena
ESTUDIANTE:
Carrión Robles Angel Daniel (23190227)
Lima
2025
Introducción
El estudio del comportamiento de los fluidos es esencial en numerosos campos de la ingeniería y la
ciencia, desde el diseño de sistemas hidráulicos hasta la aerodinámica de vehículos y estructuras. La
mecánica de fluidos proporciona modelos matemáticos fundamentales para describir y predecir el
movimiento de fluidos bajo distintas condiciones. Sin embargo, resolver las ecuaciones completas
de Navier-Stokes puede ser extremadamente complejo, por lo que en muchas situaciones prácticas
se adoptan simplificaciones válidas que permiten obtener soluciones analíticas o computacionales
eficientes.
El presente trabajo tiene como objetivo analizar el flujo laminar entre placas paralelas, mediante el
uso de herramientas computacionales en Python. A través de simulaciones interactivas, se exploran
los efectos de parámetros como la viscosidad dinámica, la densidad del fluido, el gradiente de
presión y la geometría del canal. Estas simulaciones permiten visualizar de forma intuitiva el
comportamiento del flujo, ofreciendo una forma clara de comprender los fundamentos físicos y
matemáticos subyacentes.
2. Marco Teórico
En este proyecto se aborda la solución numérica de las ecuaciones de Navier–Stokes para un flujo
laminar incompresible desarrollado en un canal bidimensional. Pero antes de ir al tema principal se
presentan las ecuaciones de Navier–Stokes, estas describen el movimiento de fluidos viscosos y surgen de
los principios de conservación de masa y cantidad de movimiento. Para fluidos newtonianos
incompresibles
estas ecuaciones se simplifican a:
Y sus simplificaciones se pueden dividir en:
2.1 Flujo Potencial Incompresible (Irrotacional)
Fundamento físico y uso
Según Batchelor (2000), Cuando se trata de flujos incompresibles con baja viscosidad (como aire a alta
velocidad lejos de las paredes), las fuerzas viscosas son insignificantes y el campo de velocidades puede
describirse mediante un potencial escalar ϕ, tal que:
El campo es irrotacional, es decir:
Esto implica que la velocidad cumple con la ecuación de Laplace, una ecuación elíptica bien conocida en
física.
Aplicaciones reales
● Diseño aerodinámico de alas (donde se aplica teoría del perfil delgado).
● Flujos externos sobre estructuras donde la capa límite es delgada.
¿Por qué se usa esta simplificación?
Según Batchelor, esta aproximación es válida cuando el número de Reynolds es alto y la capa límite es
delgada, lo cual permite despreciar la viscosidad fuera de la región cercana a las paredes. La ecuación
resultante permite analizar flujos sin tener que resolver Navier-Stokes completos.
2.2 Flujo Compresible en Toberas (Relación Área–Mach)
Fundamento físico
Según Batchelor (2000), en flujos compresibles unidimensionales, como en una tobera Laval, se aplica la
conservación de masa, momentum y energía. Para un gas ideal con flujo cuasi-unidimensional y sin
fricción, se cumple:
Este resultado relaciona el cambio de área con la velocidad y el número de Mach MMM. El flujo es
subsónico si M<1, sónico si M=1 y supersónico si M>1
Análisis físico
● Si el flujo es subsónico, una disminución del área aumenta la velocidad.
● Si el flujo es supersónico, un aumento del área aumenta la velocidad (efecto de boquilla
divergente).
¿Cómo se determinan los valores?
Usando tablas de propiedades termodinámicas y relaciones de gases ideales para T,P,ρT y considerando
que M=V/a con 𝑎 = γ𝑅𝑇
2.3 Choques en Flujos Supersónicos
Fundamento físico
Según White (2022), es un régimen supersónico, cuando ocurre una perturbación brusca (como un
obstáculo o cambio de área), se puede formar una onda de choque. A través del choque, se conservan
masa, momento y energía, pero hay aumento de entropía (flujo no reversible).
Las ecuaciones de Rankine-Hugoniot gobiernan el salto:
¿Cómo se usa en la práctica?
● Diseño de difusores y toberas de motores de reacción.
● Estimación de pérdidas por choques normales u oblicuos.
¿Por qué se considera un caso particular?
Porque el flujo cambia bruscamente y el análisis debe hacerse en dos zonas: antes y después del choque,
aplicando leyes de conservación.
2.4 Flujo de Poiseuille (entre Placas Paralelas)
Fundamento físico
Según White (2022) cuando el flujo es lento y completamente desarrollado entre dos placas (flujo
laminar), se puede suponer que la velocidad varía sólo en la dirección perpendicular a las placas. Se parte
de Navier–Stokes con estas simplificaciones:
● Estacionario.
● Unidimensional: u=u(y)
● Sin componentes verticales ni aceleración.
Solución:
donde H es la separación entre placas. Es un perfil parabólico con velocidad máxima en el centro.
Análisis físico
● La viscosidad impone una distribución suave de velocidades.
𝑑𝑝
● El gradiente de presión 𝑑𝑥
impulsa el flujo.
● Se usa en microfluidos, canales de laboratorio y validación de software CFD
Fundamentos matemáticos del flujo laminar entre placas paralelas
El caso clásico de flujo laminar incompresible entre dos placas planas paralelas (flujo de Poiseuille
planar) se obtiene al simplificar las ecuaciones de Navier–Stokes bajo ciertas suposiciones. Para un fluido
newtoniano incompresible y estacionario, con densidad constante ρ, las ecuaciones de conservación de
masa y momento en dos dimensiones son:
donde u(x,y) y v(x,y) son las componentes de la velocidad en x y, p(x,y) es la presión y μ la viscosidad
dinámica. Para un canal bidimensional entre placas horizontales separadas una distancia 2h, se considera
adicionalmente:
● Flujo unidireccional totalmente desarrollado: sólo existe componente u(y) en la dirección x
(v=0), y el perfil de velocidad no cambia con x (∂u/∂x=0).
● Ausencia de gradientes verticales de presión: por la geometría, ∂p/∂y=0
● Condición de no deslizamiento (no-slip): la velocidad en las paredes es cero (u(±h)=0).
Bajo estas suposiciones las ecuaciones se simplifican dramáticamente. El término convectivo u ∂u/∂x se
anula por ∂u/∂x=0, y la ecuación de continuidad se satisface trivialmente. Entonces la ecuación de
momento en x queda (ver p. 86–87 de Fox & McDonald (1998)):
Es decir,
con dp/dx constante (flujo impulsado por gradiente de presión fijo). Integrando esta Ecuación de
Poiseuille plana con las condiciones u(h)=u(−h)=0 se obtiene el perfil parabólico clásico (ecuación de
Poiseuille):
Perfil de Poiseuille
Las placas en y=±h cumplen la condición de no-deslizamiento (u=0). El signo de dp/dx negativo indica
que la presión disminuye en la dirección del flujo, de modo que la velocidad resultante es positiva. Este
perfil parabólico es la solución analítica para flujo estacionario y laminar entre placas paralelas fijas, con
velocidad máxima en el centro del canal (en y=0).
En simulaciones numéricas o en formulaciones que parten de la ecuación de Navier–Stokes dividida por
la densidad, también se puede expresar el perfil de Poiseuille en función del gradiente de presión por
unidad de masa. Esto se logra al incluir la densidad ρ\rhoρ y trabajar con la viscosidad cinemática ν=μ/ρ.
Así, el perfil de velocidad entre dos placas paralelas separadas una distancia H se expresa como:
Esta versión es útil cuando el término −1ρd/pdx representa la aceleración volumétrica aplicada al fluido.
La velocidad máxima se alcanza en el centro del canal y se anula en las paredes, respetando la condición
de no deslizamiento. Además podemos calcular la velocidad máxima en el centro y=H/2, y nula en las
paredes. con la siguiente fórmula.
Aplicaciones industriales del flujo laminar en canales
Según Whitesides (2006), el flujo laminar en canales estrechos tiene aplicaciones relevantes en
múltiples industrias debido a su comportamiento predecible y estable. En el campo de la
microfluídica, por ejemplo, dispositivos como los lab‑on‑a‑chip aprovechan el flujo laminar para
manipular líquidos en escalas microscópicas, permitiendo que fluidos distintos fluyan en capas
paralelas sin mezclarse turbulentamente, lo cual es crucial en aplicaciones biomédicas y de análisis
químico de alta precisión .
Por otro lado, según Xie, Zhao y Zhang (2021), en ingeniería térmica los intercambiadores de calor
con microcanales utilizan este tipo de flujo para mejorar la disipación térmica; aunque el régimen
laminar presenta menor mezcla que el turbulento, la gran relación área/volumen y la proximidad de
las paredes permiten una transferencia de calor eficiente, especialmente si se agregan
perturbaciones controladas como cavidades o pulsaciones .
Con lo expresado por Wong et al. (2016), en biotecnología las cámaras de flujo paralelo se emplean
para estudiar el comportamiento celular bajo cizalladura controlada, simulando condiciones
fisiológicas como el flujo sanguíneo, y permiten observar cómo células como las endoteliales
responden y se alinean al flujo, lo que resulta útil para investigaciones biomédicas .
Finalmente, según San Andrés (2010), en el ámbito de la tribología el flujo laminar es fundamental
en sistemas de lubricación como cojinetes y rodamientos; la película de lubricante entre superficies
móviles se analiza usando modelos derivados de las ecuaciones de Navier–Stokes para regímenes
laminares, lo que permite predecir presiones hidrodinámicas y evitar el contacto directo entre
metales .
3. Metodología
3.1 Herramientas empleadas
Para el desarrollo del proyecto se utilizó el lenguaje de programación Python, apoyado en
bibliotecas especializadas como NumPy para cálculos numéricos, Matplotlib y Seaborn para
visualización de datos, y ipywidgets para la creación de controles interactivos. Todo el trabajo fue
realizado en el entorno Jupyter Notebook, lo cual facilitó la integración entre código, texto, gráficos
y simulaciones.
3.2 Procedimiento
Se diseñó un código interactivo que permite modificar dinámicamente los siguientes parámetros del
flujo:
● Viscosidad dinámica (𝜇)
● Gradiente de presión (𝑑𝑝/𝑑𝑥)
● Altura del canal (𝐻)
● Densidad del fluido (𝜌)
A partir de estos datos, el programa calcula la velocidad máxima teórica del flujo laminar y genera
automáticamente el perfil de velocidades. Esta visualización permite estudiar cómo influyen los
distintos parámetros físicos en la forma del flujo.
3.3 Validación del modelo
La herramienta desarrollada fue validada verificando que los resultados obtenidos por el código
coinciden con los cálculos realizados manualmente mediante la fórmula teórica del flujo laminar en
un canal plano. Para ello, se introdujeron valores de prueba y se comprobó que los valores
generados por la simulación concuerdan con las expectativas teóricas, garantizando la confiabilidad
del modelo implementado.
4. Desarrollo del código y simulación
En esta sección se presentan los desarrollos computacionales implementados para simular y analizar el flujo
laminar. Se trabajó con dos lenguajes distintos:
● Un código en MATLAB, orientado a la visualización 3D animada del flujo laminar en una tubería
circular.
● Un código en Python, enfocado en la simulación interactiva y validación cuantitativa del perfil de
velocidades en flujo de Poiseuille en canal plano.
Los códigos completos se incluyen en los anexos del presente trabajo.
4.1 Código en MATLAB: Simulación animada del flujo laminar
Este código tiene como propósito representar visualmente cómo se comportan las partículas dentro
de un canal con flujo laminar. A través de una animación en 3D, se observa el desplazamiento
ordenado de partículas desde la entrada hasta la salida.
Partes principales del código:
● Definición del canal:
Se establecen las
dimensiones físicas del
canal.
● Construcción de la
geometría cilíndrica de la
tubería:
Esta sección genera la superficie 3D de la tubería: theta_cil varía de 0 a 2π, cubriendo todo el contorno
circular (360°). z_cil va desde 0 hasta la longitud total del tubo L. Xc, Yc, Zc convierten de coordenadas
cilíndricas a cartesianas.
● Representación de partículas y su velocidad
Aquí se dibujan las partículas dentro del tubo: (X, Y, Z) son sus coordenadas. 20 es el tamaño de cada
partícula. U es la velocidad axial (en el eje Z), usada
como color (más rápida = más caliente).
'filled' las dibuja sólidas.
● Bucle de animación de flujo
Ciclo infinito de animación: Zt = mod(Z +
U * 0.002 * frame, L) actualiza la posición
axial según la velocidad, simulando
movimiento real. mod(..., L) mantiene las
partículas dentro del dominio (regeneran al
inicio). drawnow actualiza el
gráfico.pause(0.015) da un efecto de
animación fluida.
frame es el contador de tiempo discreto.
● Compilación de Código
Animación de un tubo con partículas en movimiento. El color indica la velocidad, siendo mayor en
el centro. El punto central marca la velocidad máxima umaxen el eje del tubo.
4.2 Código en Python: Cálculo y visualización del perfil de Poiseuille
Este código permite al usuario calcular y visualizar el perfil de velocidad de un flujo laminar plano,
ajustando parámetros físicos clave mediante controles interactivos.
Partes principales del código:
● Importación de librerías:
Se importan librerías
científicas y de
visualización, así como
widgets interactivos.
● Fórmula implementada:
Se calcula la velocidad local en función de la posición vertical yyy, usando la solución teórica de
flujo laminar entre placas paralelas. La densidad ρ se incluye como parámetro adicional para
análisis, aunque no aparece directamente en esta fórmula de velocidad.
● Controles interactivos:
Se definen controles para cambiar la viscosidad, gradiente de presión, altura del canal y densidad.
● Visualización:
El perfil se gráfica dinámicamente
conforme el usuario ajusta los
parámetros.
● Compilación de Código
Cuadro de variación de Parámetros
❖ Altura (H): 0.05 – 0.50 m
Define el tamaño del
dominio transversal del
canal.
❖ Viscosidad cinemática (ν):
1×10⁻⁶ – 1×10⁻² m²/s
Relación entre la viscosidad
dinámica y la densidad del
fluido: ν = μ / ρ. Afecta la
forma del perfil de velocidad,
especialmente en flujos no
estacionarios.
❖ Gradiente de presión (dp/dx): -10 a -1000 Pa/m
La diferencia de presión por unidad de longitud, impulsa el flujo. Puede calcularse como:
dp/dx
❖ Densidad del fluido (ρ): 500 – 2000 kg/m³
Relación masa/volumen del fluido. Aunque no afecta directamente el perfil parabólico de
velocidad en flujo totalmente desarrollado y laminar (como el de Poiseuille), sí es relevante
para analizar otros aspectos como la energía, el número de Reynolds o el caudal másico.
● Perfil de Poiseuille
Este perfil representa la velocidad u(y) del fluido en función de la posición vertical y, mostrando su forma
parabólica típica. La velocidad es máxima en el centro del canal (punto rojo) y nula en las paredes.
● Mapa de Calor: y vs x vs u
Este mapa muestra cómo varía la velocidad uuu a lo largo de la sección transversal del tubo. El color más
cálido indica mayor velocidad (centro del tubo) y el más frío menor velocidad (cerca de las paredes
Resultados y Análisis
Comprobaciòn de funcionamiento de Código Python
● Cambio de viscosidad
Solo se variaron los parámetros de
viscosidad dinámica, según la teoría
estudiada estos deben ser inversamente
proporcionales
Resultado: Al aumentar la Viscosidad la
Velocidad max ,así como la velocidad en
cada punto, disminuye
● Cambio de densidad:
Solo se variaron los parámetros de
gradiente de presión, según la teoría
este parámetro es directamente proporcional a la velocidad
Resultado: se corrobora que al aumentar
la gradiente de presión; la velocidad max
también aumenta
Caso específico de fluido
Caso: Aceite SAE 30 a 100 °C
❖ Viscosidad cinemática (ν): 3.0e-6 ❖ Gradiente de presión (dp/dx): -200
m²/s Pa/m
❖ Densidad (ρ): 875 kg/m³ Altura del canal (H): 0.03 m
Viscosidad dinámica (μ): 0.002625 ❖ Resolución en x (Nx): 100
Pa·s (calculada como ν × ρ) ❖ Resolución en y (Ny): 50
Cálculos Realizados
Conclusiones
En este proyecto se implementó correctamente una solución simplificada de las ecuaciones de
Navier-Stokes para el caso de flujo laminar entre placas paralelas, obteniéndose un perfil parabólico
de velocidades que concuerda con la solución analítica teórica. Las simulaciones realizadas en
Python permitieron verificar el comportamiento físico del flujo bajo la influencia de distintos
parámetros, tales como el gradiente de presión, la viscosidad del fluido y la altura del canal. Se
observó que al aumentar el gradiente de presión (en valor absoluto), la velocidad máxima del flujo
también aumenta, mientras que al incrementar la viscosidad, las velocidades disminuyen
notablemente, evidenciando una mayor resistencia interna del fluido. Asimismo, se comprobó que
la altura del canal tiene un efecto cuadrático sobre la velocidad máxima, lo que resalta la
sensibilidad del sistema ante cambios geométricos. Aunque la densidad no interviene directamente
en la fórmula de velocidad, es necesaria para calcular la viscosidad dinámica cuando se parte de la
cinemática. Finalmente, se concluye que el uso de herramientas computacionales como Python y
MATLAB permitió no solo automatizar los cálculos, sino también visualizar el comportamiento del
flujo de manera clara e interactiva, facilitando la comprensión práctica de la dinámica de fluidos.
Referencias:
Batchelor, G. K. (2000). An Introduction to Fluid Dynamics (reissue ed.). Cambridge University Press.
Recuperado de
https://www.cambridge.org/core/books/an-introduction-to-fluid-dynamics/18AA1576B9C579CE25621E8
0F9266993
White, F. M., & Majdalani, J. (2022). Viscous Fluid Flow (4th ed.). McGraw‑Hill Education. Recuperado
de https://www.mheducation.com/highered/product/Viscous-Fluid-Flow-White.html
Elbadawy, I. (2023). Reliability of different nanofluids and different micro‑heat channel system
configurations. Processes, 11(3), 652. https://doi.org/10.3390/prcs110300652
Fluigent. (2021). Microfluidics for Cell Analysis. Recuperado de
https://www.fluigent.com/research/applications/cell-analysis
San Andres, L. (2009). Fundamentals of fluid film journal bearing [PDF]. OakTrust, Texas A&M
University. Recuperado de https://oaktrust.library.tamu.edu
Wendl, B. (2016?). Fluid Mechanics: Steady laminar flow between two infinite parallel plates [PDF].
Recuperado de https://wendl.weebly.com
Whitesides, G. M. (2006). The origins and the future of microfluidics. Nature, 442(7101), 368–373.
https://www.nature.com/articles/nature05058
Xie, G. F., Zhao, X. M., & Zhang, Y. (2021). Hydraulic and thermal performance of microchannel heat
sinks. Micromachines, 12(3), 245. https://doi.org/10.3390/mi12030245
Wong, A. K., Llanos, P., Boroda, N., Rosenberg, S. R., & Rabbany, S. Y. (2016). A parallel‑plate flow
chamber for mechanical characterization of endothelial cells exposed to laminar shear stress. Cellular
and Molecular Bioengineering, 9(1), 127–138. https://doi.org/10.1007/s12195-015-0424-5
Anexos
https://drive.google.com/drive/folders/1amzAyQNB78eVFaR_C8saGAFuWV-OW3xu
