“…As far as the domain tessellation is concerned, automatic adaptive mesh techniques are recognised as an ideal tool (see, e.g., [43]), unless the position of the solution's features is known a priori. In particular, anisotropic grids represent an optimal choice to strike a balance between solution accuracy and numerical efficiency, in the presence of highly directional phenomena (see, e.g., [4,12,18,20,24,37,42]). Anisotropic meshes allow the sharp detection of the solution's directionalities by properly tuning the size, the shape, and the orientation of the mesh elements, in contrast to an isotropic setting, where only the element size can be adjusted.…”