diff --git a/dumux/freeflow/navierstokes/staggered/localresidual.hh b/dumux/freeflow/navierstokes/staggered/localresidual.hh index 6baa841272ba21006e8d85be608344bf15cfe85d..b797a4b64ac916e63abff9b273b2c4ced633e7f3 100644 --- a/dumux/freeflow/navierstokes/staggered/localresidual.hh +++ b/dumux/freeflow/navierstokes/staggered/localresidual.hh @@ -34,6 +34,14 @@ namespace Dumux { +namespace Impl { +template<class T> +static constexpr bool isRotationalExtrusion = false; + +template<int radialAxis> +static constexpr bool isRotationalExtrusion<RotationalExtrusion<radialAxis>> = true; +} // end namespace Impl + // forward declaration template<class TypeTag, DiscretizationMethod discMethod> class NavierStokesResidualImpl; @@ -172,7 +180,7 @@ public: // Axisymmetric problems in 2D feature an extra source terms arising from the transformation to cylindrical coordinates. // See Ferziger/Peric: Computational methods for fluid dynamics chapter 8. // https://doi.org/10.1007/978-3-540-68228-8 (page 301) - if constexpr (ModelTraits::dim() == 2 && std::is_same_v<Extrusion, RotationalExtrusion<Extrusion::radialAxis>>) + if constexpr (ModelTraits::dim() == 2 && Impl::isRotationalExtrusion<Extrusion>) { if (scvf.directionIndex() == Extrusion::radialAxis) {