Skip to content

Remove the BARTH_JESPERSEN slope limiter #2911

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
ManasBagul23 wants to merge 3 commits into from
+1 −42
Closed

Remove the BARTH_JESPERSEN slope limiter #2911

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
d4407e4
Implement the Barth-Jespersen limiter function
ManasBagul23 Sep 16, 2026
5f1ac96
Only skip Barth-Jespersen limiting for a zero projection
ManasBagul23 Sep 17, 2026
97b3a4b
Remove the BARTH_JESPERSEN slope limiter
ManasBagul23 Sep 18, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 0 additions & 2 deletions Common/include/option_structure.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1030,7 +1030,6 @@ enum class LIMITER {
NISHIKAWA_R4 , /*!< \brief Slope limiter using Nishikawa's R4 method (stencil formulation). */
NISHIKAWA_R5 , /*!< \brief Slope limiter using Nishikawa's R5 method (stencil formulation). */
VENKATAKRISHNAN_WANG , /*!< \brief Slope limiter using Venkatakrisnan method, eps based on solution (stencil formulation). */
BARTH_JESPERSEN , /*!< \brief Slope limiter using Barth-Jespersen method (stencil formulation). */
VAN_ALBADA_EDGE , /*!< \brief Slope limiter using Van Albada method (edge formulation). */
SHARP_EDGES , /*!< \brief Slope limiter using sharp edges. */
WALL_DISTANCE /*!< \brief Slope limiter using wall distance. */
Expand All @@ -1042,7 +1041,6 @@ static const MapType<std::string, LIMITER> Limiter_Map = {
MakePair("NISHIKAWA_R4", LIMITER::NISHIKAWA_R4)
MakePair("NISHIKAWA_R5", LIMITER::NISHIKAWA_R5)
MakePair("VENKATAKRISHNAN_WANG", LIMITER::VENKATAKRISHNAN_WANG)
MakePair("BARTH_JESPERSEN", LIMITER::BARTH_JESPERSEN)
MakePair("VAN_ALBADA_EDGE", LIMITER::VAN_ALBADA_EDGE)
MakePair("SHARP_EDGES", LIMITER::SHARP_EDGES)
MakePair("WALL_DISTANCE", LIMITER::WALL_DISTANCE)
Expand Down
3 changes: 0 additions & 3 deletions Common/src/CConfig.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7302,9 +7302,6 @@ void CConfig::SetOutput(SU2_COMPONENT val_software, unsigned short val_izone) {
case LIMITER::VENKATAKRISHNAN_WANG:
cout << "Venkatakrishnan-Wang slope-limiting method, with constant: " << Venkat_LimiterCoeff << "." << endl;
break;
case LIMITER::BARTH_JESPERSEN:
cout << "Barth-Jespersen slope-limiting method." << endl;
break;
case LIMITER::VAN_ALBADA_EDGE:
cout << "Van Albada slope-limiting method implemented by edges." << endl;
break;
Expand Down
31 changes: 0 additions & 31 deletions SU2_CFD/include/limiters/CLimiterDetails.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -137,37 +137,6 @@ struct LimiterHelpers
};


/*!
* \brief Barth-Jespersen specialization.
* \ingroup FvmAlgos
*/
template<>
struct CLimiterDetails<LIMITER::BARTH_JESPERSEN>
{
su2double eps2;

/*!
* \brief Set a small epsilon to avoid divisions by 0.
*/
template<class... Ts>
inline void preprocess(Ts&...) {eps2 = LimiterHelpers<>::epsilon();}

/*!
* \brief No geometric modification for this kind of limiter.
*/
template<class... Ts>
inline su2double geometricFactor(Ts&...) const {return 1.0;}

/*!
* \brief Venkatakrishnan function with a numerical epsilon.
*/
inline su2double limiterFunction(size_t, su2double proj, su2double delta) const
{
return LimiterHelpers<>::venkatFunction(proj, delta, eps2);
}
};


/*!
* \brief Venkatakrishnan specialization.
* \ingroup FvmAlgos
Expand Down
5 changes: 0 additions & 5 deletions SU2_CFD/include/limiters/computeLimiters.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -73,11 +73,6 @@ if (geometry.GetnDim() == 2) {\
END_SU2_OMP_FOR
break;
}
case LIMITER::BARTH_JESPERSEN:
{
INSTANTIATE(LIMITER::BARTH_JESPERSEN);
break;
}
case LIMITER::VENKATAKRISHNAN:
{
INSTANTIATE(LIMITER::VENKATAKRISHNAN);
Expand Down
2 changes: 1 addition & 1 deletion config_template.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1603,7 +1603,7 @@ KIND_MUSCL_RAMP= ITERATION
% Values range from -1 (fully one-sided) to 1 (central difference).
MUSCL_KAPPA_FLOW= 0.0
%
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG, BARTH_JESPERSEN, VAN_ALBADA_EDGE,
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG, VAN_ALBADA_EDGE,
% NISHIKAWA_R3, NISHIKAWA_R4, NISHIKAWA_R5)
SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
%
Expand Down
Loading