# 3.2.18. Non-Conformal DG Boundary Condition

## 3.2.18.1. Conceptual Overview

The non-conformal boundary condition uses the DG approach described by Domino [80] and is currently implemented for turbulence models, continuity, momentum, and heat conduction. The non-conformal boundary condition is applied where you have two domains, and , which share a discontinuous interface with individual sidesets, and . The algorithm is applied in two passes, first iterating over all integration points in and finding the matching face in , then by iterating over the integration points in and finding the matching face in . The generic flux of a scalar, at an integration point on is

(3.531)

where is the diffusive flux, is the mass flux, and is the interior penalty coefficient.

Prior investigations have shown that pressure oscillations can be minimized by using the current integration point normal direction for both diffusive fluxes, so .

The penalty term, is given by

(3.532)

where is the diffusive flux coefficient and is an element length scale.

The advection coefficient, , defines the degree of upwinding to use. A value of results in a fully upwind scheme, while results in a central difference scheme. When using a hybrid approach, this value is calculated locally based on the Peclet number.

### 3.2.18.1.1. Continuity

The mass flow rate at the non-conformal boundary includes the pressure stabilization terms, as

(3.533)

## 3.2.18.2. Performance Considerations

There is a computational cost associated with the use of the non-conformal interface. This is largely due to two tasks: a search to match integration points with opposing faces on both non-conformal boundaries, and the resulting changes to the linear system stencil if the interface moves. Preliminary testing has shown that the cost of reinitializing the linear system with a new stencil is at least an order of magnitude greater than the cost of the search. For this reason, the algorithm implemented in Fuego will do an extra search in order to only reinitialize the linear system when the stencil actually changes. The user can expand the search boxes used in the stencil definition in order to reduce the number of linear system reinitializations by setting the “Search Expansion Factor” in the non-conformal boundary condition specification. This number is the approximate diametrical size increase in the stencil in terms of number of elements.

## 3.2.18.3. Non-Conformal Moving Walls

The non-conformal BC in Fuego supports having a moving surface as one side of the interface. This is implemented by providing a shell block containing nodal displacements on one side of the interface. Then non-conformal search will search for matching faces in opposing blocks and shells. If a shell is found, an appropriate wall boundary condition is applied at the integration point, while if a non-shell element is found the interior algorithms described previously are used.