fractional step method in cylindrical coordinates, help please?
I am trying to solve developing internal flow in a circular duct using a staggered grid and cylindrical navier stokes equations.
At the moment I am using a fractional step method with primitive variables but my code is currently wrong and I think it has to do with the pressure correction step and I think I might be using the wrong equations. Please let me know if you see anything wrong. My process is listed below.
fractional step for advection terms
fractional step for diffusion terms
where Pressure is calculated by
so I'm pretty sure i've narrowed it down to be an error with the poisson pressure solver for cylindrical coordinates.
Is there anyone out there who is familiar marker and cell cfd in cylindrical? or with the mit18086_navierstokes code written by benjamin seibold? I modified his code to cylindrical coordinates.
Also is there anyone out there who might consider looking at my code? I have been staring at it for so long that I might be missing something because I am so used to it.
If someone knows of a simple test case that I can run to debug my poisson solver, please let me know.
|All times are GMT -4. The time now is 23:49.|