We show that in the variational multiscale framework, the weak enforcement of essential boundary conditions via Nitsche's method corresponds directly to a particular choice of projection operator. The consistency, symmetry and penalty terms of Nitsche's method all originate from the fine-scale closure dictated by the corresponding scale decomposition. As a result of this formalism, we are able to determine the exact fine-scale contributions in Nitsche-type formulations. In the context of the advection–diffusion equation, we develop a residual-based model that incorporates the non-vanishing fine scales at the Dirichlet boundaries. This results in an additional boundary term with a new model parameter. We then propose a parameter estimation strategy for all parameters involved that is also consistent for higher-order basis functions. We illustrate with numerical experiments that our new augmented model mitigates the overly diffusive behavior that the classical residual-based fine-scale model exhibits in boundary layers at boundaries with weakly enforced essential conditions.@en

We develop the general form of the variational multiscale method in a discontinuous Galerkin framework. Our method is based on the decomposition of the true solution into discontinuous coarse-scale and discontinuous fine-scale parts. The obtained coarse-scale weak formulation includes two types of fine-scale contributions. The first type corresponds to a fine-scale volumetric term, which we formulate in terms of a residual-based model that also takes into account fine-scale effects at element interfaces. The second type consists of independent fine-scale terms at element interfaces, which we formulate in terms of a new fine-scale "interface model." We demonstrate for the one-dimensional Poisson problem that existing discontinuous Galerkin formulations, such as the interior penalty method, can be rederived by choosing particular fine-scale interface models. The multiscale formulation thus opens the door for a new perspective on discontinuous Galerkin methods and their numerical properties. This is demonstrated for the one-dimensional advection-diffusion problem, where we show that upwind numerical fluxes can be interpreted as an ad hoc remedy for missing volumetric fine-scale terms.@en

We initiate the study of the discontinuous Galerkin residual-based variational multiscale (DG-RVMS) method for incorporating subgrid-scale behavior into the finite element solution of hyperbolic problems. We use the one-dimensional viscous Burgers equation as a model problem, as its energy dissipation mechanism is analogous to that of turbulent flows. We first develop the DG-RVMS formulation for a general class of nonlinear hyperbolic problems with a diffusion term, based on the decomposition of the true solution into discontinuous coarse-scale and fine-scale components. In contrast to existing continuous variational multiscale methods, the DG-RVMS formulation leads to additional fine-scale element interface terms. For the Burgers equation, we devise suitable models for all fine-scale terms that do not use ad hoc devices such as eddy viscosities but instead directly follow from the nature of the fine-scale solution. In comparison to single-scale discontinuous Galerkin methods, the resulting DG-RVMS formulation significantly reduces the energy error of the Burgers solution, demonstrating its ability to incorporate subgrid-scale behavior in the discrete coarse-scale system.@en