1
Modelado de flujo incompresible utilizando
las ecuaciones de Navier-Stokes (febrero de
2022)
Diego Fernando Muñoz Arce 2010032, Santiago Escobar Escobar 2010062
A continuación, se presenta una introducción al modelo de
Resumen discretización utilizado.
El objetivo principal de este artículo es la resolución numérica
de las ecuaciones de Navier-Stokes, para un fluido incomprensible.
Para esto, se presenta el planteamiento numérico con una
discretización por medio del método de diferencias finitas para la II. PARTE TEÓRICA
resolución del acoplamiento entre la velocidad y la presión. Los
resultados obtenidos con la herramienta computacional
desarrollada fueron similares a los esperados. Las ecuaciones de Navier Stokes-Stokes modelan los
movimientos de los fluidos, para representar dichas ecuaciones
Palabras clave: diferencias finitas, Navier-Stokes. Computacionalmente utilizamos el método numérico de
derivación de Newton o Diferencias Finitas que nos lleva a
Abstract utilizar la técnica de la grilla. En el caso de este articulo los
The main goal of this paper is the numerical solution of the resultados están planteados en dos dimensiones y su
Navier-Stokes equations for an incompressible flow. A numerical
approach with a finite differences discretization technique to solve
discretización está dada por las siguientes igualdades
the coupling between velocity and pressure. The results given by
the computational tool are very similar to the expected.
Key words: Finite Differences, Navier-Stokes.
Índice de Términos - Alrededor de cuatro palabras o frases clave
en orden alfabético, separadas por comas. Para obtener una lista
de palabras claves sugeridas, envíe un correo en blanco a
keywords@[Link] o visite el sitio web de IEEE en
[Link]
t
I. INTRODUCCIÓN
E l modelado de flujo turbulento utiliza modelos III. IMPLEMENTACIÓN
matemáticos y computacionales. Los modelos computacionales
y matemáticos unidos se implementan como métodos Se desarrolló un código en Lenguaje Python implementando
numéricos. La discretización de un volumen continuo permite 2 matrices para la velocidad del fluido en X y Y. Se definieron
aproximar soluciones utilizando herramientas los valores de las variables en una configuración inicial para
computacionales. lograr una convergencia rápida, siendo estos vo = 1 (velocidad
Este capítulo presenta conceptos generales detrás de las inicial), h = 1 (tamaño de paso), Ny = 80 y Nx = 30
simulaciones como antecedentes de referencia realizados. (componentes de las matrices U y W, omega = 1.5 (velocidad
Primero se presentan las ecuaciones que gobiernan el flujo angular).
incompresible. Con diferentes configuraciones se pudo observar que el
tiempo de convergencia aumentaba como por ejemplo, si se
2
aumentaba omega = 2, la convergencia se lograría en más de herramienta desarrollada es capaz de resolver las ecuaciones de
200 iteraciones. Con un valor de 1.5 se logra en algo más de 40 Navier-Stokes para flujos laminares e incompresibles. De este
iteraciones y con un valor de 1 se generaba un overflow en el modo, las pruebas numéricas corroboraron las hipótesis
código. planteadas en la solución del modelo matemático y las
Para expresar los resultados gráficamente se implementó la aproximaciones de los métodos en la solución numérica.
librería Matplotlib.
WEBGRAFÍA
IV. RESULTADOS [1] [Link]
100X2013000200003
[2] [Link]
A continuación, veremos 3 de las gráficas de las 50 %20Computational%20Methods%20for%20Fluid%20Dynamics,%203r
iteraciones que se implementaron perteneciendo estas a la d%20Ed%20-%[Link]
[3] [Link]
primera iteración, iteración 25 y última iteración
respectivamente
Se puede apreciar en colores la velocidad del fluido siendo
entre más claro, más rápido.
IX. CONCLUSIÓN
En este documento se estudió el comportamiento de las
ecuaciones de Navier-Stokes en su primera parte haciendo uso
de la discretización por el método de las diferencias finitas. En
primera instancia se presentó el procedimiento para la
obtención del campo de velocidades del fluido incompresible.
El modelo numérico aplicado usa esquemas numéricos para la
discretización espacial.
De las comparaciones de los resultados obtenidos con las
soluciones usadas como referencia, se puede deducir que la