Adaptación de un método pseudo-espectral a condiciones de contorno no-libres

BRIOZZO, C.; RAMOS, I.

Córdoba, Argentina

Taller; VII Taller de Ciencias Planetarias; 2014

Observatorio Astronómico de la Universidad Nacional de Córdoba

Los métodos pseudoespectrales basados en transformada de Fourier compleja presentan rasgos muy atractivos, como una sencilla discretización de los operadores que sin embargo posee resolución espectral compleja, y la simplicidad y velocidad de los códigos resultantes, lo que los vuelve muy convenientes para el estudio de flujos turbulentos. Sin embargo su aplicación a problemas con condiciones de contorno no-libres requiere en general del uso de otras bateas (transformada seno, Chebyshev,etc.) donde muchas de estas ventajas desaparecen. En este trabajo mostraremos cómo estos métodos pueden modificarse para adaptarse a condiciones de contorno de Dirichlet o Neumann, sin perder ninguna de sus ventajas. Como ejemplo aplicaremos el código resultante a la convección a alto número de Rayleigh (R~19^9) en una celda de Rayleigh-Bénard.