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

using mac vel within wall function. #874

Closed
wants to merge 3 commits into from
Closed

Conversation

mbkuhn
Copy link
Contributor

@mbkuhn mbkuhn commented Jul 12, 2023

  • enables use of n+1/2 velocity when it's called for the viscous terms

* enables use of n+1/2 velocity when it's called for the viscous terms
@@ -210,9 +210,6 @@ void ABLVelWallFunc::wall_model(
// const amrex::Real vv = vold_arr(i, j, k, 1);
const amrex::Real uu =
0.5 * (umac_arr(i, j, k) + umac_arr(i + 1, j, k));
// Is this reasonable, to interpolate using the wall
// value? Or should I extrapolate using the two points
// above the wall?
const amrex::Real vv =
0.5 * (vmac_arr(i, j, k) + umac_arr(i, j + 1, k));
Copy link
Contributor

Choose a reason for hiding this comment

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

Please make sure that you don't mix vmac_arr with umac_arr !

Copy link
Contributor Author

Choose a reason for hiding this comment

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

good catch!

Comment on lines +353 to +354
// const amrex::Real uu = vold_arr(i, j, k, 0);
// const amrex::Real vv = vold_arr(i, j, k, 1);

Check notice

Code scanning / CodeQL

Commented-out code Note

This comment appears to contain commented-out code.
Comment on lines +209 to +210
// const amrex::Real uu = vold_arr(i, j, k, 0);
// const amrex::Real vv = vold_arr(i, j, k, 1);

Check notice

Code scanning / CodeQL

Commented-out code Note

This comment appears to contain commented-out code.
@github-actions
Copy link

This PR is stale because it has been open 30 days with no activity. Remove stale label or comment or this will be closed in 7 days.

@github-actions
Copy link

This PR was closed because it has been stalled for 7 days with no activity.

@github-actions github-actions bot closed this Aug 20, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants