[en] Incompressible flow simulations are crucial for numerous scientific and engineering applications, from automotive aerodynamics to biomedical device design. Traditional finite volume (FV) methods, such as the cell-centered finite volume (CCFV) approach used by OpenFOAM, face significant challenges related to mesh quality, particularly non-orthogonality and skewness. These issues often result in accuracy loss and lead to a complex and time-consuming mesh generation procedure, especially for complex geometries. This thesis addresses these issues by developing and implementing a face-centered finite volume (FCFV) solver for viscous laminar incompressible flows within the OpenFOAM framework. The FCFV method can be seen as a hybridisable discontinuous Galerkin (HDG) method using the lowest order approximation. It achieves first-order convergence for velocity, pressure, and velocity gradient without requiring reconstruction, thus significantly reducing mesh sensitivity. Additionally, the face-centered finite volume method circumvents the Ladyzhenskaya-Babuška-Brezzi (LBB) condition, enabling monolithic solution strategies that solve for velocity and pressure simultaneously or staggered approaches based on algebraic splitting that avoid the introduction of non-physical boundary conditions, unlike standard OpenFOAM solvers. The primary contribution of this thesis is the seamless integration of an efficient FCFV solver into OpenFOAM, offering a robust alternative for cases dealing with complex geometries and/or distorted meshes. This development includes both steady-state and transient formulations and introduces a novel hybrid pressure FCFV formulation to enhance accuracy for higher Reynolds number simulations. Extensive benchmarking and validation against well-known test problems, always compared with standard OpenFOAM solvers, demonstrate the robustness and accuracy of the FCFV solvers, particularly on meshes with high non-orthogonality and skewness. The results indicate that the FCFV method always maintains optimal first-order convergence for all variables and outperforms the standard CCFV solvers currently available within OpenFOAM in cases with distorted meshes. The developed FCFV solvers offer a valuable alternative for OpenFOAM users dealing with complex geometries, enabling the use of unstructured and stretched meshes without sacrificing accuracy or stability.

[es] Las simulaciones de flujos incompresibles son esenciales para numerosas aplicaciones en ciencia e ingeniería, como por ejemplo en el diseño biomédico o la aerodinámica de vehículos. Los métodos tradicionales de volúmenes finitos (FV), como el utilizado en OpenFOAM centrado en celdas (CCFV), requieren mallas de una cierta calidad, siendo especialmente relevante la ortogonalidad de las caras y la distorsión de las celdas. Estas restricciones generan pérdidas de precisión y, por lo tanto, en geometrías complejas, llevan a procedimientos de generación de mallas tediosos y costosos. Esta tesis aborda estas cuestiones desarrollando una metodología de volúmenes finitos centrada en caras (FCFV) para flujos incompresibles laminares, implementada directamente en OpenFOAM. El método FCFV es una variante de bajo orden del método mixto e híbrido de Galerkin discontinuo, conocido por sus siglas en inglés HDG (Hybridizable Discontinuous Galerkin). Esta metodología FCFV muestra convergencia de primer orden para velocidad, presión y gradiente de velocidad sin necesidad de ninguna reconstrucción reduciendo de esta forma la sensibilidad a la malla. Es importante resaltar que el FCFV satisface la condición de Ladyzhenskaya-Babuška-Brezzi (LBB), evitando, a diferencia de la metodología estándar de OpenFOAM, condiciones de contorno no físicas en la presión y el uso obligatorio de métodos de paso fraccionado para cumplir con LBB. Por consiguiente, el FCFV puede emplear formulaciones monolíticas para velocidad-presión. La principal contribución de esta tesis es la integración eficiente de la metodología FCFV en OpenFOAM, ofreciendo una alternativa robusta para geometrías complejas y/o mallas distorsionadas. Se han implementado formulaciones para problemas estacionarios y transitorios. Además, se propone una novedosa formulación de presión híbrida, también de orden bajo y basada en caras, para mejorar la precisión en simulaciones con altos números de Reynolds. Asimismo, se presenta una extensa serie de pruebas y validaciones que se comparan con la metodología estándar de OpenFOAM, demostrando la robustez y precisión de FCFV, especialmente en mallas con falta de ortogonalidad y alta distorsión. Los resultados indican que el FCFV mantiene la convergencia óptima de primer orden y supera a los CCFV estándar en mallas distorsionadas. La metodología FCFV desarrollada ofrece una alternativa valiosa para usuarios de OpenFOAM con geometrías complejas, permitiendo el uso de mallas no estructuradas y estiradas sin sacrificar precisión ni estabilidad.