Skip to content
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.

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

Subcell: Add TciOnFdGrid mutator for NewtonianEuler #3199

Merged

Conversation

nilsdeppe
Copy link
Member

Proposed changes

Add the troubled-cell indicator on the FD grid for NewtonianEuler.

Upgrade instructions

Code review checklist

  • The code is documented and the documentation renders correctly. Run
    make doc to generate the documentation locally into BUILD_DIR/docs/html.
    Then open index.html.
  • The code follows the stylistic and code quality guidelines listed in the
    code review guide.
  • The PR lists upgrade instructions and is labeled bugfix or
    major new feature if appropriate.

Further comments

@kidder kidder changed the title Subcell: Add TciOnDgGrid mutator for NewtonianEuler Subcell: Add TciOnFdGrid mutator for NewtonianEuler May 25, 2021
Copy link
Member

@kidder kidder left a comment

Choose a reason for hiding this comment

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

you can squash

min(get(get<Pressure>(dg_grid_prim_vars)))) <
min_pressure_allowed or
evolution::dg::subcell::persson_tci(dg_mass_density, dg_mesh,
persson_exponent, 1.0e-18) or
Copy link
Contributor

Choose a reason for hiding this comment

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

Can you make a constexpr double persson_tci_epsilon = 1.0e-18; as in the other TCI?

* minimum density and pressure on both the DG and subcell mesh are above
* \f$10^{-18}\f$, and runs the Persson TCI on the density and pressure on the
* DG grid. The reason for applying the Persson TCI to both the density and
* pressure is to flag cells at contact discontinuities.
Copy link
Contributor

Choose a reason for hiding this comment

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

Update this to match the newer text in the other TCI

Also, can you explain why you check both grids in a function "TciOnFdGrid" ?

@nilsdeppe
Copy link
Member Author

I've rebased and pushed a fixup. Thanks for the reviews!

Copy link
Member

@kidder kidder left a comment

Choose a reason for hiding this comment

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

you have my permission to squash (including the one additional requested change)

* present in the FD solution.
* - runs the Persson TCI on the density and pressure on the DG grid. The reason
* for applying the Persson TCI to both the density and pressure is to flag
* cells at contact discontinuities. The Persson TCI is only works with
Copy link
Member

Choose a reason for hiding this comment

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

is only -> only

@fmahebert
Copy link
Contributor

LGTM, please squash

@nilsdeppe
Copy link
Member Author

I've squashed. Thanks for reviewing :)

@nilsdeppe nilsdeppe merged commit 9f4dff0 into sxs-collaboration:develop Jun 2, 2021
@nilsdeppe nilsdeppe deleted the subcell_ne_tci_on_fd_grid branch June 2, 2021 21:23
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.

3 participants