| Texto completo | |
| Autor(es): |
Carvalho, Pablo G. S.
;
Devloo, Philippe R. B.
;
Gomes, Sonia M.
Número total de Autores: 3
|
| Tipo de documento: | Artigo Científico |
| Fonte: | FINITE ELEMENTS IN ANALYSIS AND DESIGN; v. 242, p. 20-pg., 2024-09-17. |
| Resumo | |
A two-level version for a recent semi-hybrid-mixed finite element approach for modeling Stokes and Brinkman flows is proposed. In the context of a domain decomposition of the flow region Omega, composite divergence-compatible finite elements pairs in H(div,Omega)xL(2)(Omega) are utilized for discretizing velocity and pressure fields, using the same approach previously adopted for two-level mixed Darcy and stress mixed elasticity models. The two-level finite element pairs of spaces in the subregions may have richer internal resolution than the boundary normal trace. Hybridization occurs by the introduction of an unknown (traction) defined over element boundaries, playing the role of a Lagrange multiplier to weakly enforce tangential velocity continuity and Dirichlet boundary condition. The well-posedness of the method requires a proper choice of the finite element space for the traction multiplier, which can be achieved after a proper velocity FE space enrichment with higher order bubble fields. The method is strongly locally conservative, yielding exact divergence-free velocity fields, demonstrating pressure robustness, and facilitating parallel implementations by limiting the communication of local common data to at most two elements. Easier coupling strategies of finite elements regarding different polynomial degree or mesh widths are permitted, provided that mild mesh and normal trace consistency properties are satisfied. Significant improvement in computational performance is achieved by the application of static condensation, where the global system is solved for coarse primary variables. The coarse primary variables are a piecewise constant pressure variable over the subregions, velocity normal trace and tangential traction over subdomain interfaces, as well as a real number used as a multiplier ensuring global zero-mean pressure. Refined details of the solutions are represented by secondary variables, which are post-processed by local solvers. Numerical results are presented for the verification of convergence histories of the method. (AU) | |
| Processo FAPESP: | 22/04168-2 - CEA/EDF/INRIA Numerical Analysis Summer School |
| Beneficiário: | Sônia Maria Gomes |
| Modalidade de apoio: | Auxílio à Pesquisa - Reunião - Exterior |
| Processo FAPESP: | 17/15736-3 - Centro de Pesquisa em Engenharia em Reservatórios e Gerenciamento de Produção de Petróleo |
| Beneficiário: | Denis José Schiozer |
| Modalidade de apoio: | Auxílio à Pesquisa - Programa Centros de Pesquisa em Engenharia |
| Processo FAPESP: | 22/16636-0 - Computação de alto desempenho com métodos de decomposição de domínio multiescala: escoamentos em meios porosos |
| Beneficiário: | Pablo Giovanni Silva Carvalho |
| Modalidade de apoio: | Bolsas no Brasil - Pós-Doutorado |