Estoy simulando un flujo incompresible sobre un cilindro en el número de Reynold de 500. Estoy resolviendo la ecuación de Stoke usando el método de corrección de presión. Mi solución se vuelve inestable después de cierto tiempo (aproximadamente 5 segundos).
He intentado refinar mi malla, stepsize (0.05) (asegurándome de que mi CFL <1, aunque estoy usando métodos implícitos)
Mis condiciones de contorno, malla y resultados inestables se muestran en las figuras adjuntas. El dominio es aproximadamente 25 veces más grande que el diámetro del cilindro.
He intentado simular este problema O grid (que se volvió inestable casi de inmediato).
El siguiente enlace contiene las imágenes de las condiciones de contorno y los resultados.
Estaría agradecido si alguien puede compartir sus pensamientos / experiencias sobre este problema. Muchas gracias.
editado:
Disculpas por el error de tipeo:
Estoy usando las siguientes condiciones de contorno: Límite de Neumann
en el límite de Dirichlet
editado:
He aplicado condiciones de límite de velocidad en los nodos alrededor del límite de dirichlet. Además, el nodo de esquina superior derecha e inferior derecha es el límite de dirichlet con velocidad 1.
Después, analicé más profundamente los resultados de la simulación, noté que la inestabilidad comienza a colarse en la unión de entrada / salida.