Unconditionally stable operator splitting algorithms for the incompressible magnetohydrodynamics system discretized by a stabilized finite element formulation based on projections
In this article, we propose different splitting procedures for the transient incompressible magnetohydrodynamics (MHD) system that are unconditionally stable. We consider two levels of splitting, on one side we perform the segregation of the fluid pressure and magnetic pseudo-pressure from the vecto...
| Autores: | , , |
|---|---|
| Tipo de recurso: | artículo |
| Estado: | Versión aceptada para publicación |
| Fecha de publicación: | 2013 |
| País: | España |
| Institución: | Universidad de Sevilla (US) |
| Repositorio: | idUS. Depósito de Investigación de la Universidad de Sevilla |
| OAI Identifier: | oai:idus.us.es:11441/168317 |
| Acceso en línea: | https://hdl.handle.net/11441/168317 https://doi.org/10.1002/nme.4392 |
| Access Level: | acceso abierto |
| Palabra clave: | Incompressible magnetohydrodynamics Operator splitting algorithms Fractional step methods Stabilized finite element methods Symmetric projection stabilization |
| Sumario: | In this article, we propose different splitting procedures for the transient incompressible magnetohydrodynamics (MHD) system that are unconditionally stable. We consider two levels of splitting, on one side we perform the segregation of the fluid pressure and magnetic pseudo-pressure from the vectorial fields computation. At the second level, the fluid velocity and induction fields are also decoupled. This way, we transform a fully coupled indefinite multi-physics system into a set of smaller definite ones, clearly reducing the CPU cost. With regard to the finite element approximation, we stick to an unconditionally convergent stabilized finite element formulation because it introduces convection stabilization, allows to cir-cumvent inf-sup conditions (clearly simplifying implementation issues), and is able to capture nonsmooth solutions of the magnetic subproblem. However, residual-based finite element formulations are not suitable for segregation, because they lose the skew-symmetry of the off-diagonal blocks. Therefore, in this work, we have proposed a novel term-by-term stabilization of the MHD system based on projections that is still unconditionally convergent. |
|---|