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

2D Poiseuille flow #61

Open
lucafadigati opened this issue Oct 10, 2022 · 3 comments
Open

2D Poiseuille flow #61

lucafadigati opened this issue Oct 10, 2022 · 3 comments

Comments

@lucafadigati
Copy link

Good morning,

blastFoam has been tested on simple case: 2D Poiseuille flow, where the pressure drop has been set between the inlet and the outlet of the 2D pipe. The maximum velocity expected at the channel centerline should be 1.24997 mm/s (such value has been computed analitically using 2D Poiseuille formula and then verified with rhoCentralFoam solver); actually blastFoam computed a velocity of 1.3773 mm/s. The simulation files are provided in attachment.

Thanks,
Luca

2DPoiseuilleCase.zip

@jheylmun
Copy link
Collaborator

Hi Luca,
Thank you for bringing this to my attention (and sorry for the long delay). I am still trying to figure out why the results differ, since every thing should be the same. I will let you know I do figure out what is causing the discrepancies though. Out of curiosity, have you tried a higher Mach number flow as well?

Thanks,
Jeff

@carlo-fiorina
Copy link

Hi. Was this issue resolved?

@lucafadigati
Copy link
Author

Hi,
I am sorry but I have been busy and forgot to check github.
I have run a simple simulation inside a channel, the flow is subsonic at the inlet and it becomes supersonic after the nozzle throat.
I tried to use blastReactingFoam because in the futures I would run some simulations with chemical reactions. The issue that I have found is related to the viscosity. I am imposing a noSlip condition on the walls patch but the velocity does not go to zero, it remain constant across the nozzle section.
Maybe I am using your solver in the wrong way. I have imposed a static pressure at the inlet and "pressureInletOutletVelocity" for U field.
I know your solver was designed for explosion but I don't understand what mistake I made in the simulation settings.

Best regards,
Luca

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

3 participants