A Polynomial-Correction Navier-Stokes Characteristic Boundary Condition
Edmond K. Shehadi,Edwin T.A. van der Weide
DOI: https://doi.org/10.1016/j.cpc.2024.109194
IF: 4.717
2024-04-11
Computer Physics Communications
Abstract:In an effort to assess the efficacy of different non-reflecting boundary conditions (NRBCs) in high-order solvers, several NRBCs have been implemented in the context of the discontinuous Galerkin (DG) discretization. The methods considered are based on buffer zone techniques and boundary-imposed conditions. The former mainly utilizes the sponge layer (SL), the perfectly matched layer (PML) and the supersonic layer (SSL), while the latter uses a Riemann extrapolation (RE) technique and a novel polynomial-correction Navier-Stokes characteristic boundary condition (PC-NSCBC). The developed PC-NSCBC method reconstructs an explicit boundary state, instead of modifying the equations on the boundary or specifying ad hoc compatibility conditions on corners and edges. This property pairs well with a Riemann solver, as is typical in a DG method. Moreover, the boundary state can be assembled either from a one-dimensional perspective (NSCBC1D) or using partial knowledge of multi-dimensional effects (NSCBC3D). Four different benchmark problems are considered: a 1D linear acoustic pulse in a quiescent flow, a 2D non-linear pressure pulse in a quiescent flow, an isentropic vortex and a shear-flow vortex shedding instability. Results imply that mainly the PML performs better than most other approaches. This is especially true for acoustically-dominated flows. Alternatively, for flow conditions predominantly less acoustic in nature and strongly non-linear, instabilities in the PML may arise. The SL and SSL both suggest that a larger layer is required to deliver comparable results to the PML. The proposed (PC-)NSCBC methods demonstrate reasonable success with no instabilities reported, which is impressive given the comparatively less computational effort required. Partial inclusion of the transverse terms in the NSCBC3D noticeably enhances the output of the NSCBC1D, albeit further investigation is required to determine an optimal, universal expression for its transverse tuning coefficient.
physics, mathematical,computer science, interdisciplinary applications