Andrzej Jaeschke

Matthias Möller, Andrzej Jaeschke

This work extends the high-resolution isogeometric analysis approach established for scalar transport equations to the equations of gas dynamics. The group finite element formulation is adopted to obtain an efficient assembly procedure for the standard Galerkin approximation, which is stabilized by adding artificial viscosities proportional to the spectral radius of the Roe-averaged flux-Jacobian matrix. Excess stabilization is removed in regions with smooth flow profiles with the aid of algebraic flux correction \cite{KBNII}. The underlying principles are reviewed and it is shown that linearized FCT-type flux limiting \cite{Kuzmin2009} originally derived for nodal low-order finite elements ensures positivity-preservation for high-order B-Spline discretizations.

Andrzej Jaeschke, Matthias Möller

Isogeometric analysis was applied very successfully to many problem classes like linear elasticity, heat transfer and incompressible flow problems but its application to compressible flows is very rare. However, its ability to accurately represent complex geometries used in industrial applications makes IGA a suitable tool for the analysis of compressible flow problems that require the accurate resolution of boundary layers. The convection-diffusion solver presented in this chapter, is an indispensable step on the way to developing a compressible flow solver for complex viscous industrial flows. It is well known that the standard Galerkin finite element method and its isogeometric counterpart suffer from spurious oscillatory behaviour in the presence of shocks and steep solution gradients. As a remedy, the algebraic flux correction paradigm is generalized to B-Spline basis functions to suppress the creation of oscillations and occurrence of non-physical values in the solution. This work provides early results for scalar conservation laws and lays the foundation for extending this approach to the compressible Euler equations.