Estoy trabajando con la biblioteca OpenFOAM C ++ Computational Continuum Mechanics (puede lidiar con la interacción fluido-sólido, flujos MHD ...) que utiliza mallas arbitrarias no estructuradas. Esto fue impulsado por la idea de utilizar la ventaja de la generación rápida (generalmente automática) de mallas no estructuradas para simular problemas en geometrías complejas.
Sin embargo, recientemente me he encontrado con otro enfoque: mallas cartesianas adaptativas de octree con "corte" de celdas, donde el refinamiento agresivo de la malla se utiliza para describir una geometría compleja.
Desde el punto de vista numérico, las mallas carthesianas son mucho más precisas, por lo que mi pregunta es: ¿alguien tiene experiencia en el uso / implementación de uno o ambos enfoques? ¿Cómo se comparan entre sí?
Estoy desarrollando códigos para el flujo de fluido en dos fases y noté, por ejemplo, que la reconstrucción de los gradientes de campo se puede hacer fácilmente más precisa en las mallas carthesianas, mientras que la malla no estructurada requiere una regresión lineal para cambios abruptos en el campo ...