A stabilized finite element method for modeling dispersed multiphase flows using orthogonal subgrid scales

We propose a finite-element formulation for simulating multi-component flows occupying the same domain with spatially varying concentrations. Each constituent is assumed to behave as an incompressible Newtonian fluid, and solutions are sought for the velocities and volume fractions of each phase, as...

Descripción completa

Detalles Bibliográficos
Autores: Gravenkamp, Hauke|||0000-0002-2641-6384, Codina, Ramon|||0000-0002-7412-778X, Principe, Ricardo Javier|||0000-0002-1478-2651
Tipo de recurso: artículo
Fecha de publicación:2024
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/407185
Acceso en línea:https://hdl.handle.net/2117/407185
https://dx.doi.org/10.1016/j.jcp.2024.112754
Access Level:acceso abierto
Palabra clave:Computational fluid dynamics
Fluid dynamics
Multiphase flow
Stabilized finite elements
Variational multiscale method
Orthogonal subgrid scales
Dispersed flow
Dinàmica de fluids computacional
Àrees temàtiques de la UPC::Física::Física de fluids::Flux de fluids
Descripción
Sumario:We propose a finite-element formulation for simulating multi-component flows occupying the same domain with spatially varying concentrations. Each constituent is assumed to behave as an incompressible Newtonian fluid, and solutions are sought for the velocities and volume fractions of each phase, as well as the common pressure. Stabilization terms are derived within the framework of the variational multiscale method based on an approximation of the finite-element residual to achieve control of the pressure and volume fractions. We utilize the concept of term-by-term stabilization in conjunction with orthogonal subgrid scales, thus incorporating only those terms of the residual essential to obtain stability and projecting them on a space orthogonal to the finite element space. The resulting system of equations is solved in a monolithic manner, requiring a small number of nonlinear iterations. Several benchmark tests have been performed to confirm the stability and optimal asymptotic convergence rates for linear and higher-order elements using the proposed formulation.