Resumen

En este trabajo se implementa un método Lagrange-Galerkin estabilizado por presión en un sistema paralelo de descomposición de dominio, y se demuestra que la nueva estrategia de estabilización es efectiva para simulaciones con números de Reynolds y números de Rayleigh grandes. La simetría de la matriz de rigidez permite resolver los problemas de interfaz del sistema lineal mediante el método conjugado precondicionado, y se aplica un precondicionador de dominio equilibrado incompleto a los problemas acoplados de flujo térmico. La metodología muestra una buena eficiencia paralela y una alta escalabilidad numérica, y el nuevo solucionador se valida mediante la comparación con soluciones exactas y resultados de referencia disponibles. Ocupa menos memoria que los solucionadores clásicos de tipo producto; además, es capaz de resolver problemas de más de 30 millones de grados de libertad en un día en un clúster de PC de 80 núcleos.

• Materias:Ecuaciones diferenciales Dinámica de fragmentación Medios porosos Circuito abierto
• Subjects:Differential equations Fragmentation dynamics Porous media Open circuit
• Palabras claves:Lagrange-galerkin; Descomposición de dominio paralelo; Estrategia de estabilización; Número de reynolds; Número de rayleigh; Escalabilidad numérica
• Keywords:Lagrange-galerkin; Parallel domain decomposition; Stabilization strategy; Reynolds number; Rayleigh number; Numerical scalability