Monthly Archives: July 2014

Lección 2.3

Lección 2.3 — grid.txt

Los parámetros considerados para la solución de las ecuaciones de gobierno que están en la guía de usuario (Users’ Guide) están especificados en el grid.txt. En el  grid.txt, todas las líneas están constituidas por 3 bloques, exceptuando la línea separadora. Por ejemplo, las tres primeras líneas del grid.txt son algo parecido a lo siguiente:

———- Control File for Flowsquare ver 4.0 (Use SI Unit) ———-
————————- General Control Data ————————
01:cmode 0 // Simulation mode,….

Las dos primeras líneas son separadoras. Así como en la tercera línea, consiste en 3 bloques, menos las separadoras. El primer bloque (01:cmode) denota el nombre de cada parámetro (frecuentemente con nombres similares a la notacion utilizada en Users’ Guide). El segundo bloque (“0″) es el número que se especifica para cada parámetro. Siendo el tercer bloque (// Simulation mode, …) un comentario, el cual no tiene funcionalidad en la simulación. Estos tres bloques deben estar separados por espacio(s).

Algunos de los parámetros en el grid.txt están relacionados con las condiciones de frontera (bc.bmp) y/o las condiciones iniciales (ic.bmp). Aquí tenemos una lista de parámetros especificados en el grid.txt y sus descriptores

  1. ———- Control File for Flowsquare ver 4.0 (Use SI Unit) ———-
  2. ————————- General Control Data ————————
  3. 01:cmode
    Modo de simulación, donde:

    • cmode=0: Simulación de flujo incompresible No-reactivo,
    • cmode=1: Simulación de flujo reactivo premezclado
    • cmode=2: Simulación de flujo reactivo No- premezclado
    • cmode=3: Simulación de flujo Subsónico/supersónico
  4. 02:nx
    Número de puntos del mallado en la direccióndediscretización en x (horizontal).
  5. 03:ny
    Número de puntos del mallado en la direccióndediscretización en y (vertical).
  6. 04:lx
    Dimensión del dominio físico (m) en la dirección x (horizontal).
  7. 05:ly
    Dimensión del dominio físico (m) en la dirección  y (vertical).
  8. 06:sts
    Inicio o reinicio del tiempo de paso de la simulación. Las simulaciones nuevas siempre parten de un sts=0. Se puede reiniciar una simulación desde el punto de salida, desde los archivos de descarga en la carpeta dump.
  9. 07:latts
    Ultimo tiempo de paso de la simulación. (se puede reiniciar la simulación.)
  10. 08:cflfac
    Factor de tiempo de paso (dt: Incremento de tiempo físico durante la simulación). En Flowsquare, dt=lx/(nx-1)/u_max/cflfac, Donde u_max es definido como:

    • 1. Velocidad máxima definida en el grid.txt para cmode=0–2, o
    • 2. Velocidad del sonido para cmode=3.
    • Típicamente se prefiere un cflfac=10~20.
  11. —————- Numerical Scheme, Smoothing & Accuracy —————
  12. 09:iorder
    Técnica numérica utilizada en la simulación. Se puede utilizar un número entre 0 y 5 para el iorder, y cada número significa lo siguiente.

    • iorder=0: Método de Bajo nivel (diferencias de 2º, integral de 1º orden )
    • iorder=1: Método de de alto orden (diferencias de 4 º, integral de 3º)
    • iorder=2: Diferencia de 2º orden y metodo de Lax-Wendroff
    • iorder=3: Diferencia de 4 º orden y metodo de Lax-Wendroff (2° orden)
    • iorder=4: WENO (Lax-Friedrichs flux splitting) y Lax-Wendroff (2° orden)
    • iorder=5: WENO (Steger-Warming flux splitting) y Lax-Wendroff (2° orden)
    • Nota: los iorder=4 y 5 se usan solo para un cmode=3.
  13. 10:nfil
    Intervalos de tiempo de paso de procedimiento de filtrado espacial.
  14. 11:wfil
    Parámetro de relajación de filtrado espacial (Ver sección 4.4 de Users’ Guide).
  15. 12:omega
    Relevante para cmode=0–2. Parámetro de relajación de la ecuación de Poisson. 1.8  es generalmente mejor para una rápida convergencia, pero un omega mas pequeño se puede usar para un flujo con geometría compleja o un flujo reactivo con desprendimiento de calor (ver Seccion 4.3 del Users’ Guide).
  16. 13:peps
    Relevante para cmode=0–2. Tolerancia de convergencia para la ecuación de Poisson (ver Seccion 4.3 del Users’ Guide).
  17. 14:loopmax
    Relevante para  cmode=0–2. Numero máximo de iteraciones. Loopmax es priorizado sobre peps (ver sección 4.3 del Users’ Guide).
  18. 15:wdrho
    Relevante para cmode=1 y 2. wdrho=1 es correcto y generalmente este numero debe ser usado. Pero durante la transicion inicial, valores menores pueden ser usados para la convergencia. 0<=wdrho<=1. ver Eq. (6) del Users’ Guide.
  19. —————— General BC and Global IC (White) ——————-
  20. 16:perikey
    Especifico para fronteras periodicas y su direccion.

    • perikey=0: No es usada una frontera periodica.
    • perikey=1: Es periodica en la direccion x (horizontal).
    • perikey=2: Es periodica en la direccion y (vertical).
    • perikey=3: Eson periodica en las direccionesx e y.
  21. 17:pres0
    En (Pa). Ver Lección 2.2 — ic.bmp
  22. 18:uin0
    En (m/s).
  23. 19:vin0
    En (m/s).
  24. 20:rho0
    En (kg/m^3).
  25. 21:temp0
    En (K).
  26. 22:scalar0
  27. —————— BLUE Local BC and/or IC (optional) —————–
  28. 23:uin1
    En (m/s). Ver Lección 2.1 — bc.bmp
  29. 24:vin1
    En (m/s).
  30. 25:rho1
    En (kg/m^3).
  31. 26:temp1
    En (K).
  32. 27:scalar1
  33. —————— RED Local BC and/or IC (optional) ——————
  34. 28:uin2
    En (m/s). Ver Lección 2.1 — bc.bmp
  35. 29:vin2
    En (m/s).
  36. 30:rho2
    En (kg/m^3).
  37. 31:temp2
    En (K).
  38. 32:scalar2
  39. ———- PINK Local BC and/or IC (pure air flow, optional) ———-
  40. 33:uin3
    En (m/s). Ver Lección 2.1 — bc.bmp
  41. 34:vin3
    En (m/s).
  42. 35:temp3
    En (K).
  43. ————— BLACK Wall Boundary Condition (optional) ————–
  44. 36:tempw
    En (K). Ver Lección 2.1 — bc.bmp
  45. ————– GREEN Moving Boundary Condition (optional) ————-
  46. 37:imb
    Ver Lección 2.1 — bc.bmp
  47. 38:umb
    En (m/s).
  48. 39:vmb
    En (m/s).
  49. 40:tempmb In (K).
  50. ————- YELLOW Scalar Boundary Condition (optional) ————-
  51. 41:scalarT
    En (K). Ver Lección 2.1 — bc.bmp
  52. ———- Transport Properties & Thermochemical Conditions ———–
  53. 42:mu
    Viscosidad dinámica de la mezcla. En (kg/m s). Agua@300K: 8.94E-4, Aire@300K: 18.6E-6. Relevante para cmode=0–2.
  54. 43:R
    Constante especifica de gases en (J/kg K). Air: 286.9. Relevante to cmode=1 and 2.
  55. 44:diff
    Didusividad de la masa de la mezcla en (m^2/s). Si es la unidad se asume el numero de Schmidt, diff~mu/rho y Aire@300K: 20.0E-6. Relevante para cmode=1, 2 y cuando es usada una condición de frontera (B.C.) amarilla.
  56. 45:Tu
    Temperatura de gases sin quemar en (K). Relevante para cmode=1.
  57. 46:Tb
    Temperatura de gases quemados en (K). Relevante para cmode=1 and 2.
  58. ——————- Chemical Reaction (for cmode=1) ——————-
  59. 47:krate
    Constante Pre-exponencial en el mecanismo químico elemental. Ver Eq. (13) de la guía de usuario Users’ Guide.
  60. 48:Trate
    Temperatura de activación en (K).
  61. 49:nrate
    Constante para el mecanismo químico elemental.
  62. 50:cF
    Progreso de la variable donde se localiza la flama (~0.5, solo para propósitos de visualización ).
  63. ————— Non-Premixed Reacting Flow (for cmode=2) ————–
  64. 51:Xst
    Fracción estequiometrica de la mezcla donde la flama se desarrolla (~0.5).
  65. 52:sigma
    Parámetro de relajación para la variación de densidad. 0<=sigma<=1 y 1 es exacto. Ver Eq. (28) del Users’ Guide.
  66. ————————– Display & Output —————————
  67. 53:box
    Tamaño del pixel de cada punto de mallado para cada punto mostrado en la pantalla. Generalmente  nx * nbox es más pequeño que el tamaño de pantalla.
  68. 54:nfig
    Intervalo de tiempo de paso para la generación de imágenes instantáneas (0: ninguna imagen es generada).
  69. 55:nfile
    Intervalo de tiempo de paso para la generación de datos instantáneos de la simulación. (0: ningún dato de simulación es generado).
  70. 56:bcdisp
    • bcdisp=0: Contorno de pared no superpuesto.
    • bcdisp=1: Contorno de pared superpuesto.
  71. 57:idisp
    Visualización de campos.

    • idisp=0: Apagado
    • idisp=1: Densidad.
    • idisp=2: Componente de velocidad en la dirección x (horizontal), u.
    • idisp=3: Componente de velocidad en la dirección y (vertical), v.
    • idisp=4: Rapidez (sqrt(u*u+v*v)).
    • idisp=5: Vorticidad
    • idisp=6: T
    • idisp=7: Gradinte de reaccion (for cmode=1)
    • idisp=8: Progreso de variables para cmode=1, fraccion de mezcla para cmode=2, escalar general para cmode=0, 3.
    • idisp=9: Presión (menos pres0)
    • idisp=10: Fraccion de la mezcla de aire y progreso de variables (cmode=1), E/rho (J/kg) (cmode=3).
  72. 58:cmax
    Valor máximo de la barra de colores para la variable del idisp dado (colocar 0 para la auto escala.)
  73. 59:cmin
    Valor mínimo de la barra de colores para la variable del idisp dado (colocar 0 para la auto escala.)
  74. 60:icolor
    Mapa de colores.

    • icolor=0: Jet
    • Jet
    • icolor=1: Arcoíris
    • Rainbow
    • icolor=2: Nishiki
    • Nishiki
    • icolor=3: Grises
    • Grey
    • icolor=4: Grises (inverso)
    • Grey_inv
    • icolor=5: Calor
    • Hot
    • icolor=6: Marino
    • Sea
    • icolor=7: Hoja
    • Leaf
  75. 61:icont
    Línea de contorno en el frente de reacción (cmode=1 y 2).

    • icont=0: Apagado
    • icont=1: Negro
    • icont=2: Rojo
    • icont=3: Verde
    • icont=4: Azul
    • icont=5: Blanco
  76. 62:linewidth
    Espesor de la línea de contorno en el frente de reacción

    • linewidth=1: 1 pixel
    • linewidth=3: 3 pixeles
    • linewidth=5: 5 pixeles
    • linewidth=7: 7 pixeles
  77. 63:ivec
    Vectores de velocidades y sus colores.

    • ivec=0: Apagado
    • ivec=1: Negro
    • ivec=2: Rojo
    • ivec=3: Verde
    • ivec=4: Azul
    • ivec=5: Blanco
  78. 64:ndiv
    Intervalos de puntos de mallado de los vectores de velocidad (0: ajuste automático).
  79. 65:vecsize
    Tamaño en pixel del vector en flecha (0: ajuste automático).
  80. —————— Lagrangian Trajectory (optional) ——————-
  81. 66:lagkey
    Trayectoria de Lagrange.

    • lagkey=0: Apagado
    • lagkey=1: Partículas iniciadas desde la frontera izquierda.
    • lagkey=2: Partículas iniciadas desde la frontera inferior.
    • lagkey=3: Partículas iniciadas desde la izquierda y derecha de la frontera.
    • lagkey=4: Partículas iniciadas desde la frontera superior e inferior.
  82. 67:lagcolor
    Color de las partículas de Lagrange. 0: Negras y 1: Blancas.
  83. 68:lagsize
    Tamaño en Pixel de las partículas.
  84. 69:nlagra
    Intervalo de tiempo de paso de inicialización de las partículas de Lagrange (>=100).
  85. 70:npart
    Número  partículas de Lagrange (>=1000).
  86. ———————— Body Force (optional) ————————
  87. 71:gfx
    Fuerza del cuerpo en la dirección x (horizontal) debido a la diferencia de densidad en (m/s^2).
  88. 72:gfy
    Fuerza del cuerpo en la dirección y (vertical) debido a la diferencia de densidad en (m/s^2).
  89. 73:dref
    Densidad de referencia.

    • dref=1: Valor máximo de densidades calculadas del grid.txt.
    • dref=2: Valor medio  de densidades calculadas del grid.txt.
    • dref=3: Valor mínimo de densidades calculadas del grid.txt.
  90. ——————- Initial Perturbation (optional) ——————-
  91. 74:pmode
    Interruptor de perturbación inicial. Ver Sección 4.6 de la Users’ Guide.

    • pmode=0: Apagado
    • pmode=1: modo simple
    • pmode=2: modo múltiple
    • pmode=3: multi nodos (amplitud al azar)
  92. 75:umag
    Amplitud de la perturbación de la velocidad en (m/s).
  93. 76:nwave
    Numero de perturbaciones en la dirección x (horizontal).
  94. ——————————– Others ——————————-
  95. 77:nwait
    Tiempo de espera (generalmente 0 para una mayor velocidad computacional).
  96. ———————————————————————–
  97. #End of file

 

Lección 2.2

Lección 2.2 — ic.bmp

Las condiciones iniciales (I.C.) son tan importantes como las condiciones de frontera. En Flowsquare sin embargo, las condiciones iniciales son fijadas deacuerdo a las condiciones de frontera (B.C.) del bc.bmp y los valores iniciales de las variables son especificadas en el grid.txt. El cmode utilizado en el siguiente, es el modo de simulación los cuales son: 0: No-reactivo, 1: reactivo (premezclado), 2: reactivo (No- premezclado), y 3: Subsonico/supersonico dispuesto en el grid.txt.

  • pres0: Presión en (Pa)
  • uin0: Componente inicial de Velocidad en dirección x (horizontal)
  • vin0: Componente inicial de Velocidad en dirección y (vertical)
  • rho0: Densidad inicial (para cmode=0 y 3)
  • temp0: Temperatura inicial (para cmode=1 y 2)
  • scalar0: Fracción de mezcla inicial (para cmode=2)

 

Así que, básicamente, en cualquier lugar en blanco en el bc.bmp se inicializa de acuerdo a los valores antes mencionados a lo largo del dominio. Para muchos de los casos, deberia ser bueno que adicionalmente no se tengan que especificar los campos. Sin embargo, para algunos casos esto no es lo apropiado. Para ciertos casos, las condiciones iniciales (I.C.) pueden ser fijadas por separado usando un mapa de imagen llamado ic.bmp, cuyo tamaño debería ser el mismo que el del dominio (nx x ny píxeles, donde nx y ny son especificadas en el gris.txt), aunque si el tamaño del ic.bmp no esta acorde con las dimensiones del dominio, Flowsquare interpolara la cifra automáticamente para la simulación. Las condiciones iniciales son definidas por diferentes colores así como en el bc.bmp, pero hay algunas diferencias.

  • Color (R,G,B): Descripción del tipo de condiciones iniciales
  • █ Negro (0,0,0): Pared. Velocidad establecida en 0 (cero). Si la temperatura de la pared no esta definida, La temperatura inicial de esta frontera es temp0.
  • Azul (0,0,255): El campo azul esta fijado igual a uin1, vin1, rho1 (cmode=0, 3), temp1 (cmode=1, 2), scalar1 (cmode=2). Estos no están en los extremos del dominio de la simulación a diferencia del bc.bmp.
  • Rojo (255,0,0): El campo rojo esta fijado igual a uin2, vin2, rho2 (cmode=0, 3), temp2 (cmode=1, 2), scalar2 (cmode=2). Estos no están en los extremos del dominio de la simulación a diferencia del bc.bmp.
  • █ Verde (0,255,0): Frontera móvil.
  • Rosa (255,0,255): Corriente de aire puro (Solo para ser usado en el modo (premezclado).
  • Amarillo (255,255,0): Escalar inicial adicional.

El archivo ic.bmp es considerado solo al inicio de cada simulación. Luego, el campo es determinado basado en la solución de las ecuaciones de transporte y el bc.bmp.

 

Lección 2.1

Lección 2.1 — bc.bmp

En esta página, aprenderemos a configurar las condiciones de contorno (siglas en ingles B.C.) para las simulaciones que usted quiera configurar a partir de cero. Sin embargo, se recomienda encontrar un caso similar en los problemas de ejemplo (Sample Problems) y modificar sus archivos de entrada (bc.bmp, grid.txt, etc) para adaptarse a su caso.

Las condiciones de contorno son importantes ya que determina la solución a partir de los límites y esta se propaga hacia todo el dominio. En Flowsquare, usaremos bc.bmp para definir las condiciones de contorno de la simulación, y bc.bmp debe estar preparado para cada cada una de ellas. Con la finalidad de establecer las condiciones de contorno, usaremos una herramienta usual de dibujo como Microsoft Paint y guardaremos la imagen como *.bmp. Por ejemplo, el bc.bmp por defecto para el caso usado en la simulación del flujo de canal 2D es algo como esto:

Sample bc.bmp used in the 2D channel flow simulation.

Ejemplo del bc.bmp usado en la simulación del  flujo de canal 2D.

La figura representa las condiciones de frontera en el dominio de la simulación. El tamaño del bc.bmp debe ser igual a la dimensión del dominio (nx x ny pixeles, donde nx y ny son especificadas en el gris.txt), aunque el tamaño del bc.bmp no coincida con el tamaño del dominio, Flowsquare interpolara la cifra automáticamente en la simulación. En la figura se observa una línea azul a la izquierda del dominio, y dos negras en la parte superior e inferior del dominio. Cada color tiene un significado específico, los siguientes colores se utilizan para especificar varias condiciones de frontera (B.C.) en Flowsquare.

 

  • Color (R,G,B): Descripción del tipo de condiciones de frontera (B.C.).
  •   Negro (0,0,0): Frontera sin deslizamiento, sin flujo y temperatura fija (Paredes).
  •  Azul (0,0,255): Frontera de flujo de entrada.
  •  Rojo (255,0,0): Frontera de flujo de entrada.
  •  Verde (0,255,0): Frontera móvil (sin flujo a temperatura fija).
  •  Rosa (255,0,255): Flujo de Aire puro (solo en modo premezclado).
  •  Amarillo (255,255,0): Frontera escalar adicional.

Específicamente las fronteras azul, roja y rosa, las cuales solo se pueden fijar en los extremos del dominio de la simulación (por ejemplo en (i, j)=(1, j), (nx, j), (i, 1), (i, ny)). Si esos colores de fronteras son usadas dentro del dominio (2<=i<=nx, and 2<=j<=ny), se consideraran como condiciones iniciales (nota: la condición inicial utilizada en bc.bmp tiene prioridad sobre ic.bmp). Cada tipo de fronteras tiene sus parámetros que deben ajustarse en el grid.txt. la mayoría de ellos son opcionales y no todos son para ser especificados por los usuarios. Lo que los usuarios especificarán para cada frontera del grid.txt son: Los nombres de las variables realmente utilizados en grid.txt También se introducen aquí. El cmode utilizado es el modo de simulación, y son:  0: no reactivo, 1: reactivo (premezclado), 2: reactivo (no premezclado) o 3: Sub/Super supersónico, en el grid.txt.

  • █ Negro (0,0,0): Frontera sin deslizamiento, sin flujo y temperatura fija (Paredes).
    tempew: Temperatura en la pared (opcional). Si es distinta de 0 (cero), la temperatura de pared es procesada por la simulación. Si es 0 (cero), la temperatura de pared variara dependiendo de los alrededores del fluido (generalmente a gradiente cero).
  • █ Azul (0,0,255): Frontera de flujo de entrada.
    uin1: Componente de Velocidad en dirección x (horizontal) en la frontera.
    vin1: Componente de Velocidad en dirección y (vertical) en la frontera.
    rho1: Densidad en la frontera (para cmode=0 y 3, si es usada esta frontera).
    temp1: Temperatura en la frontera (para cmode=1 y 2, si es usada esta frontera).
    scalar1: Fracción de Mezcla (para cmode=2, si es usada esta frontera).
  • █ Rojo (255,0,0): Frontera de flujo de entrada.
    uin2: Componente de Velocidad en dirección x (horizontal) en la frontera.
    vin2: Componente de Velocidad en dirección y (vertical) en la frontera.
    rho2: Densidad de la frontera (para cmode=0 y 3, si es usada esta frontera).
    temp2: Temperatura en la frontera (para cmode=1 y 2, si es usada esta frontera).
    scalar2: Fracción de Mezcla (para cmode=2, si es usada esta frontera).
  •  Verde (0,255,0): Frontera móvil (sin flujo a temperatura fija).
    imb: Si este es 1, el movimiento de la frontera es periódico (movimiento oscilante). Hágalo 0 (cero) para el caso contrario.
    umb: Velocidad de desplazamiento en la dirección x (horizontal).
    vmb: Velocidad de desplazamiento en la dirección y (vertical).
    tempmb: Temperatura de la frontera móvil. Si es distinta de 0 (cero), la temperatura de pared es procesada por la simulación. Si es 0 (cero), la temperatura de pared variara dependiendo de los alrededores del fluido (generalmente a gradiente cero).
  •  Rosa (255,0,255): Flujo de aire puro (Solo para ser usado en el modo  premezclado).
    uin3: Componente de Velocidad en dirección x (horizontal) en la frontera.
    vin3: Componente de Velocidad en dirección y (vertical) en la frontera.
    temp3: Temperatura en la frontera (ha de ser fijada si es usada esta frontera es usada).
  • █ Amarillo (255,255,0): Frontera escalar adicional.
    scalarT: Valor escalar de la frontera (corresponde a c en la Eq. (25) sin reaction en la guia de usuario Users’ Guide).

Siempre es más fácil echar un vistazo a varios problemas de ejemplos en Sample Problems para aprender cómo usar cada condición de contorno para conseguir lo que quieres!