Busca avançada
Ano de início
Entree


A semi-hybrid-mixed method for Stokes-Brinkman-Darcy flows with H(div)-velocity fields

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: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING; v. 125, n. 1, p. 36-pg., 2023-10-03.
Resumo

We consider new semi-hybrid-mixed finite element formulations for Stokes-Brinkman problems. Using H(div)-conforming approximate velocity fields, the continuity of normal components over element interfaces is taken for granted, and pressure is searched in discontinuous spaces preserving the divergence compatibility property. Tangential continuity is weakly imposed by a Lagrange multiplier playing the role of tangential traction. The method is strongly mass-conservative, leading to exact divergence-free simulations of incompressible flows. The Lagrange multiplier space requires specific choices according to the velocity approximations implemented in each element geometry. In certain cases, classic divergence-compatible pairs adopted for Darcy flows may require divergence-free bubble enrichment to enforce tangential continuity in some extent, avoiding any extra stabilization technique. Considerable improvement in computational performance is achieved by the application of static condensation: the global system is solved only for a piecewise constant pressure variable, velocity normal trace and tangential traction over inter faces. The remaining solution components are recovered by solving independent local Neumann problems in each element. Numerical results are presented for a set of standard test cases in the field of Stokes-Brinkman-Darcy flows with known analytical solutions. The main convergence properties of the method are verified in the whole range of parameters, from Stokes to Darcy limits, as well as for the combined Stokes-Darcy scenario. The robustness of the method is also demonstrated for more challenging test problems with complex non-conforming meshes with local refinement patterns as well as for curved geometry. (AU)

Processo FAPESP: 16/13637-5 - Advanced numerical methods: recent developments, analysis, and applications
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