Skip to content

Add optional dissipation for Div_par_mod in Fci - #3467

Open
totork wants to merge 3 commits into
nextfrom
divpar-diss
Open

Add optional dissipation for Div_par_mod in Fci#3467
totork wants to merge 3 commits into
nextfrom
divpar-diss

Conversation

@totork

@totork totork commented Aug 21, 2026

Copy link
Copy Markdown
Contributor

The current Div_par_mod operator does not have any dissipation added in Fci. This PR adds a Rusanov flux. It is optional, so that it can be turned off in physics models for e.g. testing/MMS purposes.

Used to be

#3466

but now on the boutproject repo.

@totork
totork requested a review from dschwoerer August 21, 2026 10:04
@totork

totork commented Aug 21, 2026

Copy link
Copy Markdown
Contributor Author

@dschwoerer Added the exception when it is not Fci

Comment thread include/bout/fv_ops_impl.hxx Outdated
Coordinates* coord = f_in.getCoordinates();
ASSERT1_FIELDS_COMPATIBLE(f_in, v_in);

if (!f_in.isFci() && dissipative) {

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If you move this to the non FCI branch, it will be just

Suggested change
if (!f_in.isFci() && dissipative) {
if (dissipative) {

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done in d5f4814

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants