METODOLOGIA SUGERIDA DE TRABAJO OPENFOAM

Introducción General

Previo a ejecutar la primera simulación en OpenFOAM siempre se recomienda tener un esquema general de trabajo, por tal razón OpenFOAM se sugiere previo a la preparación de casos y ejecución de los mismo cumplir los siguientes ítem.

       ·        Tener en lo posible el Área ò Volumen de control sobre el que se realizara el modelamiento bien definido.

       ·        Seleccionar de forma apropiada las condiciones de control y condiciones de contorno; dando una condición de uniformidad a los casos que se desean trabajar. 

       ·        TENER LAS HIPOTESIS DEL CASO BIEN DEFINIDAS PARA VALIDAR POSTERIOR AL PROCESAMIENTO

       ·        Tener diferentes fundamentos técnicos para la toma de decisiones en cada paso del modelo

·        TECLEAR CON CUIDADO LOS COMANDOS EN LA PLATAFORMA.

Antes de continuar se debe validar que cumplamos con los requisitos mencionados

Posteriormente presentamos la metodología general de trabajo en OpenFOAM, 

esta se encuentra basada en el diagrama que se presenta a continuación

Fuente : Tomado de (Foundation, 2014)

En donde desglosamos a continuación cada fase por separado, siendo el siguiente 

Preprocesamiento

Es aquella fase en la que se procede a elegir el Solver más adecuado para dar solución a nuestro caso  en paralelo se debe construir el mallado de la geometría a partir de los esquemas por elementos (Hexaédricos ò Triangulares)(Jasak, 2009). Permitiendo tener mayor acceso a una esquematización correcta del volumen de control que será simulado.

Para la fase de Preprocesamiento tenemos una oferta amplia de Software para construir mallado como los son

     ·        ANSA

     ·        EnGrid

     ·        GridPro

     ·        PointWise & Gridgen

     ·        Salome-Meca

     ·        FreeCAD

     ·        Gmsh

Siendo los anteriores mencionados los más estables  y Gratuitos en la plataforma UBUNTU 14.06 según (Jasak, 2009) y son ampliamente sugeridos por la red (www. cfd-online.com). En resumen, las siguientes imágenes muestran el proceso completado al que se debe llegar.

Fuente: Tomado de (Junhong, 2004)

En donde se muestra como debe transformarse de un volumen de control como es una elipse a un volumen compuesto por elementos rectangulares que en vista 3d son Hexaedros(Lucchini, 2008). Se debe resaltar que en esta fase dependiendo de la complejidad de la geometría deberá escogerse el tipo de elemento que sea más conveniente para el mallado(Nilsson, 2006).

Procesamiento

Esta fase es la más corta en términos de inversión de tiempo de trabajo y la memoria de trabajo computacional juega un papel importante pues definirá la variable [ Tiempo de Procesamiento (Seg)], siendo uno de los criterios más importantes para muchos diseñadores.

En esta fase se deben ingresar los valores de entrada para el modelo asociado, dependiendo del modelo requerirá sus respectivas variables.

En Windows se puede usar el comando para ingresar los datos con el comando ‘’mcedit’’ para editar desde la misma terminal, en Ubuntu se usa el editor ‘’gedit’’. 

Recomendamos 

*NO EDITAR LOS ARCHIVOS EN WORDPAD O BLOC DE NOTAS PUES CAMBIAN EL FORMATO DEL ARCHIVO; CAUSANDO QUE OPENFOAM NO RECONOZCA EL ARCHIVO Y SE PRODUZCA UN ERROR FATAL. En esta situación sugerimos que use el editor NOTEPAD++ (Editor Libre).

A continuación, se muestra la estructura de carpetas de trabajo de OpenFOAM, esta se debe respetar a toda costa para no tener inconvenientes en un futuro.

Fuente: Captura de Pantalla de OpenFOAM por el Autor

• Donde la carpeta 0, contiene todos los archivos de la información de entrada al modelo

• Donde la carpeta constant, contiene la información asociada a la malla ò volumen de control

• Donde la carpeta system, contiene los archivos asociados al control del Solver & Controles a archivos de post-Procesamiento.

• Posteriormente ingresando al Archivo de Velocidad editamos la velocidad en cada una de las direcciones, recordando que la Velocidad es de carácter Vectorial.

                                         Fuente: Captura de Pantalla de OpenFOAM por el Autor

Antes de continuar se debe calcular 2 campos importantes para controlar la convergencia del caso y el orden de las carpetas como se muestra a continuación.

Calculo del Co ( Numero de Courant)

Para Proceder a calcular el Courant se debe reemplazar los siguientes valores en la siguiente ecuación.

Fuente: Tomado de User-Guide OpenFOAM

• Donde Sigma es deltaT, siendo la variable objetivo para establecer antes de empezar la simulación. La ecuación despejada en términos de deltaT será

Fuente : Tomado de User-Guide OpenFOAM

• Donde Co [ Numero de Courant Deseado] , Se recomienda trabajar con valores menores a  0.5(Junhong, 2004).

• Donde Sigmax [ Longitud Unitaria /# Elementos (Faces)], En caso que sean elementos cuadrados la longitud unitaria se puede asumir como la raíz cuadrada del área mínima; Sugerimos en lo posible usar elementos finitos rectangulares pues se puede tener mayor control que un mallado Triangular; en caso que sean Triangulos se debe calcular usando la regla de (Base*Altura/2).

• Donde Velocidad [ Debe estar en la entrada del sistema y debe expresada en m/Seg],

Posteriormente se procede a guardar la información y se corre el modelo usando el comando desde la carpeta Principal

Fuente: Pantallazo Propio del autor de OpenFOAM 

Posteriormente se da Enter y esperamos a que OpenFOAM termine de procesar la información.

Calculo del Write Interval

En este paso se establece cada cuanto el solver guardara la información en el equipo, se realiza el cálculo a través de la siguiente expresión.

Fuente : Ecuacion Empirica Desarrollada por el autor

• Donde se debe establecer el tiempo de simulación final para este caso lo haremos con 0.5, posteriormente ingresamos el Sigma t. Finalmente sabemos cuántas veces deberá imprimir información, Para lo que podemos ordenar que cada # Veces determinadas que imprima resultados guarde los resultados asociados.

Aplicando la Ecuacion mostrada anteriormente deseamos tener 5 carpetas para este caso aplicando la ecuacion obtuvimos de forma eficiente la cantidad de carpetas que deseábamos

En la imagen se observan las carpetas diferentes a la carpeta 0 , constant y system son las carpetas de trabajo que deseábamos guardara el Solver , validándose de forma correcta el funcionamiento de la ecuacion planteada anteriormente.

Posteriormente ejecutamos el nombre del Solver y obtenemos

Fuente: Pantallazo Propio del autor en OpenFOAM

Como se observa en la imagen anterior se imprimir los resultados hasta que termina todo el procesamiento , una vez se consiguen estos se puede determinar si el solver está trabajando de forma estable con el caso , bajo la condición que el Courant Number mean se mantenga constante a medida que imprime un paso temporal de otro(Rusche, 2011) , también cuando el residual es muy pequeño.

PostProcesamiento

En esta fase de proceder a observar y analizar los detalles de la simulación con paraFoam , para lo que procedemos a ejecutar el siguiente comando terminada la simulación 

Fuente: Pantallazo Propio del autor en OpenFOAM

Posteriormente se abrirá el programa de visualización de datos de OpenFOAM por defecto, segundo a segundo en el tiempo que le hayamos asignado para guardar los resultados.

Fuente: Pantallazo Propio del autor en OpenFOAM

En donde podemos ver los perfiles de la variable que deseemos sobre el volumen de control, debe tenerse en cuenta que el nivel de detalle de la simulación está directamente relacionado con la calidad de refinamiento (Hacer que la malla sea más fina y reporte las variables a un nivel de detalle más alto), En consecuencia, toma un tiempo mayor de procesamiento, pero mejores resultados.

En futuros Entradas se explicará más a detalle cómo se debe realizar el proceso de validación y refinamiento de geometrías. 

Bibliografia

    1.     Foundation, O. (2014). Open∇FOAM - The Open Source CFD Toolbox - User Guide. User Guide, 222. https://doi.org/10.1023/A

     2.     Jasak, H. (2009). Dynamic Mesh Handling in OpenFOAM. 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition, 1–10. https://doi.org/10.2514/6.2009-341

      3.     Junhong, W. (2004). Open Source CFD Solver - OpenFOAM. Lecture, 1–5.

  4.     Lucchini, T. (2008). OpenFOAM programming tutorial. Talk, 1–47. https://doi.org/10.1016/0303-8467(95)00031-E

      5.     Nilsson, H. (2006). Evaluation of OpenFOAM for CFD of Turbulent Flow in Water Turbines. 23rd IAHR Symposium, 1–9.

   6.     Rusche, H. (2011). OpenFOAM Tutorials : Programming Session Tutorials : Programming Session. Training, 1–19.