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

Fix calculation of h/n for 2D cases #158

Open
wants to merge 1 commit into
base: master
Choose a base branch
from
Open

Conversation

Andres-MG
Copy link
Collaborator

We are computing h/n assuming that the cases are always 3D. With this PR, the value of h/n in 2D cases is correct as long as the extruded dimension is 1.

We use this when computing artificial viscosity and LES.

@Andres-MG Andres-MG added the bug Something isn't working label May 16, 2023
@Andres-MG Andres-MG requested a review from loganoz May 16, 2023 13:56
@Andres-MG Andres-MG self-assigned this May 16, 2023
@Andres-MG
Copy link
Collaborator Author

The test NSCylinderSmagorinsky is failing because we were using the 3D definition of delta in a 2D case. How should we handle this?

@Andres-MG
Copy link
Collaborator Author

It turns out that the issue might be in the way we check if a mesh is 2D. I will look into the implementation of CheckIfMeshIs2D:

subroutine HexMesh_CheckIfMeshIs2D(self)

We are apparently considering some meshes as 2D even if not specifically stated in the control file.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

1 participant