Algoritmo computacional de dinámica de fluidos.
En dinámica de fluidos computacional (CFD), el algoritmo SIMPLE es un procedimiento numérico ampliamente utilizado para resolver las ecuaciones de Navier-Stokes . SIMPLE es un acrónimo de Método semiimplícito para ecuaciones ligadas a presión.
El algoritmo SIMPLE fue desarrollado por el profesor Brian Spalding y su alumno Suhas Patankar en el Imperial College de Londres a principios de los años 1970. Desde entonces, muchos investigadores lo han utilizado ampliamente para resolver diferentes tipos de problemas de flujo de fluidos y transferencia de calor. [1]
Muchos libros populares sobre dinámica de fluidos computacional analizan en detalle el algoritmo SIMPLE. [2] [3]
Una variante modificada es el algoritmo SIMPLER (SIMPLE Revisado), que fue introducido por Patankar en 1979. [4]
Algoritmo
El algoritmo es iterativo . Los pasos básicos en la actualización de la solución son los siguientes:
- Establezca las condiciones de contorno.
- Calcule los gradientes de velocidad y presión.
- Resuelva la ecuación del momento discretizado para calcular el campo de velocidad intermedia.
- Calcule los flujos de masa no corregidos en las caras.
- Resuelva la ecuación de corrección de presión para producir valores de celda de la corrección de presión.
- Actualice el campo de presión: donde urf es el factor de relajación insuficiente para la presión.
![{\displaystyle p^{k+1}=p^{k}+{\text{urf}}\cdot p^{'}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Actualice las correcciones de presión límite .
![{\displaystyle p_{b}^{'}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Corregir los flujos de masa facial:
![{\displaystyle {\dot {m}}_{f}^{k+1}={\dot {m}}_{f}^{*}+{\dot {m}}_{f}^{ '}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Corregir las velocidades de las celdas: ; donde es el gradiente de las correcciones de presión, es el vector de coeficientes centrales para el sistema lineal discretizado que representa la ecuación de velocidad y Vol es el volumen de la celda.
![{\displaystyle {\vec {v}}^{k+1}={\vec {v}}^{*}-{\frac {{\text{Vol}}\ \nabla p^{'}}{ {\vec {a}}_{P}^{v}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\nabla p^{'}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {{\vec {a}}_{P}^{v}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Actualizar la densidad debido a los cambios de presión.
Ver también
Referencias