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

Inconsistent WSS in rigid fluid pipe mesh? #69

Open
schwarz-e opened this issue Mar 23, 2022 · 4 comments
Open

Inconsistent WSS in rigid fluid pipe mesh? #69

schwarz-e opened this issue Mar 23, 2022 · 4 comments

Comments

@schwarz-e
Copy link

schwarz-e commented Mar 23, 2022

I am running a simple rigid pipe simulation with a mesh generated with a prismatic boundary layer/outer wall.

Screenshot from 2022-03-23 10-56-31

However, the WSS field generated is not uniform. Ignoring effects at the end caps, within one axial slice the WSS can vary by 10-20%.

Screenshot from 2022-03-23 10-58-39

Slice taken at mid-vessel:

Screenshot from 2022-03-23 11-01-40

Here are the input file and mesh for reference:

setup.zip

@vvedula22
Copy link
Contributor

@schwarz-e While the spatial variation is not surprising when using FEM over unstructured meshes, the range is intriguing. Have you tested refining the mesh, particularly in the interior? Can you plot element WSS that is computed over the face (#62)?

@schwarz-e
Copy link
Author

I used Martin's method to calculate the WSS at the faces on a more refined mesh- the mid-slice range is a lot better (only about 5%) but there is a large axial variation when it should be constant.

image

image

image

@vvedula22
Copy link
Contributor

Axial variation is indeed strange. This didn't happen in your previous test case. Is there a chance for a bug in your element-based calculation around the inlet plane?

@schwarz-e
Copy link
Author

schwarz-e commented Dec 14, 2022

It seems like this stems from the velocity solution at the outlet which has non-zero tangential components. This doesn't seem to be a bug, its a consistent steady-state solution when Direchlet at the inlet and Neumman at the outlet is perscribed. However, it's still undesirable.

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

No branches or pull requests

2 participants