Introducción a fluidos computacionales

Discretización, estabilidad y convergencia en simulaciones introductorias de flujo.

Contexto conceptual

La dinámica de fluidos computacional transforma un flujo continuo en un conjunto finito de celdas, variables y ecuaciones algebraicas. En cada celda se aproximan $u$, $v$, presión y otras propiedades. La elección de $h$ decide qué estructuras del flujo pueden resolverse y cuánto se separa $u_num$ de $u_ref$. Una malla gruesa suaviza gradientes y capas límite; una malla fina aumenta el coste y puede exigir reducir $dt$. El problema central no es obtener una imagen vistosa, sino demostrar que $\epsilon _u$ disminuye de forma coherente al refinar la discretización sin cambiar el problema físico.

🟢 Nivel esencial

Una simulación de fluidos divide el dominio en pequeñas celdas y calcula cómo se mueve el fluido en ellas. El $h$ representa la escala espacial de esas celdas. Cuando disminuye, la malla puede describir mejor los cambios de las componentes de velocidad y de la presión, especialmente cerca de paredes o entradas. Por eso, $\epsilon _u$ suele reducirse al refinar. Sin embargo, una malla más fina no corrige condiciones de contorno erróneas, una $u_ref$ poco fiable ni un modelo físico inadecuado. También puede aparecer inestabilidad si el intervalo de avance no se adapta. La práctica correcta consiste en comparar varias mallas manteniendo iguales $Re$, geometría y condiciones de contorno. Una tendencia regular del error apoya la convergencia; una curva irregular indica que otro mecanismo está dominando. El refinamiento es, por tanto, una prueba del modelo numérico, no solo una petición de más detalle.

Además, el refinamiento debe evaluarse con una magnitud definida de forma constante. Comparar una velocidad puntual en una malla con una media en otra produciría una falsa impresión de mejora.

🔵 Nivel formal

Fórmulas básicas

La conservación de masa en un flujo bidimensional incompresible exige que el defecto local medido por $div_u$ sea nulo:

Continuidad para flujo incompresible

El régimen dinámico se caracteriza mediante $Re$, construido con $\rho$, $U$, $L$ y $\mu$:

Número de Reynolds

En simulaciones transitorias, el avance advectivo por celda se controla con $CFL$, relacionado con $dt$ y $h$:

Número CFL advectivo

La discrepancia entre $u_num$ y $u_ref$ se cuantifica mediante:

Error relativo de velocidad

Preguntas reales del alumno

¿Una malla más fina siempre mejora el resultado? Solo si el error espacial domina y el modelo, el solver y el paso temporal son adecuados.

¿Por qué debo mantener fijo $Re$? Porque define el régimen físico. Si cambia, ya no comparas la misma corriente.

¿Qué significa una $div_u$ pequeña? Que la solución satisface aproximadamente la conservación local de masa.

¿Puede una simulación estable ser poco precisa? Sí. $CFL$ adecuado evita ciertas inestabilidades, pero no garantiza resolución espacial suficiente.

¿Qué hago si la referencia vale casi cero? Usa error absoluto o una norma global; el error relativo puntual pierde sentido.

Síntesis final

La fiabilidad de una simulación de flujo se demuestra, no se supone. $h$ debe reducir $\epsilon _u$ con una tendencia coherente, mientras continuidad, régimen y estabilidad permanecen controlados. Una sola malla produce un resultado; varias mallas producen evidencia.