Skip to content

Use the ghost-cell spacing for the z one-sided viscous gradient - #1991

Draft
sbryngelson wants to merge 1 commit into
MFlowCode:masterfrom
sbryngelson:fix/viscous-z-ghost-spacing
Draft

sbryngelson wants to merge 1 commit into
MFlowCode:masterfrom
sbryngelson:fix/viscous-z-ghost-spacing

Conversation

@sbryngelson

Copy link
Copy Markdown
Member

s_compute_fd_gradient's one-sided z derivative at idwbuff(3)%beg
divided by z_cc(beg + 2) - z_cc(is3_viscous%beg). is3_viscous%beg is
idwbuff(3)%beg + 1 there, so the denominator spanned one cell instead
of two and the outermost ghost gradient was doubled on a uniform grid.

Divide by z_cc(beg + 2) - z_cc(beg), matching the x and y stencils.

The 27 3D Viscous tests pass against the existing goldens.


Acknowledgement

  • I confirm this PR meets the above expectations and reflects my own understanding and real-world context.

s_compute_fd_gradient's one-sided z derivative at idwbuff(3)%beg
divided by z_cc(beg + 2) - z_cc(is3_viscous%beg). is3_viscous%beg is
idwbuff(3)%beg + 1 there, so the denominator spanned one cell instead
of two and the outermost ghost gradient was doubled on a uniform grid.

Divide by z_cc(beg + 2) - z_cc(beg), matching the x and y stencils.

The 27 3D Viscous tests pass against the existing goldens.

This branch has not been deployed

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

Labels

None yet

Development

Successfully merging this pull request may close these issues.

1 participant