Skip to content

Flather BC - #2024

Draft
mbkuhn wants to merge 64 commits into
kynema:mainfrom
mbkuhn:flather_bc
Draft

Flather BC#2024
mbkuhn wants to merge 64 commits into
kynema:mainfrom
mbkuhn:flather_bc

Conversation

@mbkuhn

@mbkuhn mbkuhn commented Aug 28, 2026

Copy link
Copy Markdown
Contributor

Summary

Pull request type

Please check the type of change introduced:

  • Bugfix
  • Feature
  • Code style update (formatting, renaming)
  • Refactoring (no functional changes, no api changes)
  • Build related changes
  • Documentation content changes
  • Other (please describe):

Checklist

The following is included:

  • new unit-test(s)
  • new regression test(s)
  • documentation for new capability

This PR was tested by running:

  • the unit tests
    • on GPU
    • on CPU
  • the regression tests
    • on GPU
    • on CPU

Additional background

Issue Number:

Comment on lines +446 to +451
const amrex::Real c = std::sqrt(grav_z * interior_h);

const auto Flather_vel =
boundary_val + (ori.isLow() ? -1.0_rt : 1.0_rt) * c /
boundary_h * (interior_h - boundary_h);

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.

I am unsure about the fact I use interior_h to calculate c and then use boundary_h as the multiplier in Flather_vel. The formula in the literature is unclear about where to get h from (interior or boundary) but in some cases, these values are at least assumed to be the same, such that sqrt(g H)/H is written as sqrt(g / H).

, m_v_mac(m_sim.repo().get_field("v_mac"))
, m_vof(m_sim.repo().get_field("vof"))
{
// amrex::ParmParse pp(identifier());
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