Modeling of steep layers in singularly perturbed diffusion–reaction equation via flexible fine-scale basis

We present a new stabilized method for the diffusion–reaction equation which develops sharp boundary and/or internal layers for the reaction-dominated case (i.e. singularly perturbed case). The method relies on an improved expression for the stabilization parameter that is derived via the fine-scale...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Computer methods in applied mechanics and engineering 2020-12, Vol.372 (C), p.113343, Article 113343
Hauptverfasser: Masud, Arif, Anguiano, Marcelino, Harari, Isaac
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We present a new stabilized method for the diffusion–reaction equation which develops sharp boundary and/or internal layers for the reaction-dominated case (i.e. singularly perturbed case). The method relies on an improved expression for the stabilization parameter that is derived via the fine-scale variational formulation facilitated by the variational multiscale (VMS) framework. The proposed fine-scale basis consists of enrichment functions which may be nonzero at element edges. The derived stabilization parameter enjoys spatial variation over element interiors and along inter-element boundaries that helps model the rapid variation of the residual of the Euler–Lagrange equations over the domain. This feature also facilitates consistent stabilization across boundary and internal layers as well as capturing anisotropic refinement effects. The method is able to better satisfy the maximum principle as compared to other existing methods. We present a priori mathematical analysis of the stability and convergence of the method for the diffusion–reaction equation. Optimal convergence on meshes comprised of linear triangles and bilinear quadrilateral elements are presented for smooth problems, as well as for problems with steep boundary layers. Stability and accuracy features of the method for problems with discontinuous forcing function, internal layers, and boundary layers are shown and its performance on unstructured and distorted meshes comprised of quadrilateral and triangular elements is highlighted. •VMS stabilized method for singularly perturbed diffusion–reaction equation.•Fine-scale basis with multiple local functions that are nonzero at element edges.•Method eliminates oscillations, captures steep layers sharply, and adapts to directionality.•A priori mathematical analysis of the stability and convergence of the method.•Optimal convergence and stability shown on structured, stretched, and unstructured meshes.
ISSN:0045-7825
1879-2138
DOI:10.1016/j.cma.2020.113343