# Vortex Mehtod : Internal flow

 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 21, 2005, 23:28 Vortex Mehtod : Internal flow #1 MC Ryu Guest   Posts: n/a Dear , I'm making a code for predicting lid-driven cavity flow using the conventional Lagrangian vortex method. In my results, many particles go outside the box domain, especially near the upper two corners. I think the reason is that there are singular behaviours there. But except the moving wall vorticity and secondary vortex near two lower corners, overall flow patterns and the centerline velocity profiles are nealry same with Ghia's resutls after converging. I think the wrong prediction of the moving wall vorticity profile is due to the outside particles. Am I right ? And please let me know how can I reduce or eliminates the particles going outside ? Thank you, cf) dt = 0.01, Blob size = 1/64 for Re = 100

 June 22, 2005, 00:10 Re: Vortex Mehtod : Internal flow #2 Adrin Gharakhani Guest   Posts: n/a Are you using random walk for diffusion? This would explain so many particles outside the box. How are you dealing with the velocity boundary conditions and diffusion from the wall? A good/correct implementation would not lead to so many particles going out of the domain. And, anyway, particles that move outside should be removed from the computation, as keeping them will lead to incorrect flow field. Adrin Gharakhani

 June 22, 2005, 02:14 Re: Vortex Mehtod : Internal flow #3 MC Ryu Guest   Posts: n/a Thank you for your kindness. I use PSE scheme for diffusion. And for no-slip condition at the wall, diffusion process is iterated. So, strengths of the particels are modified. "Particles that move outside should be removed from the computation" means that outside particels must be eliminated from the computation even though they have large enough strength or they should be forced to have zero strength ? Thank you.

 June 22, 2005, 09:42 2) Vortex Mehtod : Internal flow #4 MC Ryu Guest   Posts: n/a Thank you, I cured some mistakes in my program, so particles don't go outside. But, it is seen some wrong vectors in the region of the right bottom corner, very close to the corner. In the position direction of the velocity vector is north-west, so in the corner region streamline contours are slightly different from Ghia's solution. Please give me some advices.

 June 22, 2005, 19:04 Re: Vortex Mehtod : Internal flow #5 Adrin Gharakhani Guest   Posts: n/a Removing means removing particles from the computations. Setting the strengths to zero has essentially the same effect, but adds to your cost. A zero strength is equivalent to the case where no particles exist - of course when using PSE you will see some nearby particles gaining strength eventually, due to diffusion. There is nothing wrong with removing the particles. What happens is that the loss of vorticity gets compensated for in the next timestep when you apply the no-slip BC. Essentially, you'll have a vortex sheet distribution at the walls with strength that accommodate the loss of particles. Nevertheless, if you are losing particles, then you may still have errors in implementation and/or your timestep size may be too large. A few particles here and there is not a problem, but an accurate solution should not lead to more losses. Adrin Gharakhani

 June 22, 2005, 19:15 Re: 2) Vortex Mehtod : Internal flow #6 Adrin Gharakhani Guest   Posts: n/a Well, since in one day you found problems that significantly improved your errors, I wonder if you will find other errors should you spend some more time!! As for the vector going "in the wrong direction", I'm not so sure, because I don't have the vector field in front of me to make sense. However, barring you have additional errors (quite likely), I don't see a problem with the direction of the vector you described. You're supposed to have a very small secondary vortex there - it just makes physical sense - and the vector direction you mention, sort of, agrees with that. Try higher and lower resolution runs (_after_ you're sure you have no bugs) and see what happens. I won't be surprised one bit that you can capture small scales better than a finite-difference solution! Adrin Gharakhani

 June 22, 2005, 22:07 Thank you !!! Mr. Gharakhani #7 MC Ryu Guest   Posts: n/a I'm very appreciated in your kind answers. I'll try to do following your comments. Thank you again.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post JPBodner Main CFD Forum 1 September 3, 2011 03:07 smschnob Main CFD Forum 1 November 11, 2010 18:12 diaw Main CFD Forum 104 February 16, 2006 06:44 ravi FLUENT 2 March 16, 2005 01:45 Oliver Main CFD Forum 4 September 13, 2000 17:47