We carry out a stability and convergence analysis for the fully discrete scheme obtained by combining a finite or virtual element spatial discretization with the upwind-discontinuous Galerkin time-stepping applied to the time-dependent advection–diffusion equation. A space–time streamline-upwind Petrov–Galerkin term is used to stabilize the method. More precisely, we show that the method is inf–sup stable with constant independent of the diffusion coefficient, which ensures the robustness of the method in the convection- and diffusion-dominated regimes. Moreover, we prove optimal convergence rates in both regimes for the error in the energy norm. An important feature of the presented analysis is the control in the full L2(0,T;L2(Ω)) norm without the need of introducing an artificial reaction term in the model. We finally present some numerical experiments in (3+1)-dimensions that validate our theoretical results.
BEIRAO DA VEIGA, L., Dassi, F., Gomez, S. (2025). SUPG-stabilized time-DG finite and virtual elements for the time-dependent advection–diffusion equation. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 436(1 March 2025) [10.1016/j.cma.2024.117722].
SUPG-stabilized time-DG finite and virtual elements for the time-dependent advection–diffusion equation
Lourenco Beirao da Veiga;Franco Dassi;Sergio Gomez
2025
Abstract
We carry out a stability and convergence analysis for the fully discrete scheme obtained by combining a finite or virtual element spatial discretization with the upwind-discontinuous Galerkin time-stepping applied to the time-dependent advection–diffusion equation. A space–time streamline-upwind Petrov–Galerkin term is used to stabilize the method. More precisely, we show that the method is inf–sup stable with constant independent of the diffusion coefficient, which ensures the robustness of the method in the convection- and diffusion-dominated regimes. Moreover, we prove optimal convergence rates in both regimes for the error in the energy norm. An important feature of the presented analysis is the control in the full L2(0,T;L2(Ω)) norm without the need of introducing an artificial reaction term in the model. We finally present some numerical experiments in (3+1)-dimensions that validate our theoretical results.File | Dimensione | Formato | |
---|---|---|---|
Beirao da Veiga-2025-Comp Meth Appl Mech Engineer-VoR.pdf
accesso aperto
Tipologia di allegato:
Publisher’s Version (Version of Record, VoR)
Licenza:
Creative Commons
Dimensione
2.55 MB
Formato
Adobe PDF
|
2.55 MB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.