I need to prepare or to find, for educational purposes a, so smal as possible, source code (Delphi, Pascal or Fortran I can understand) which enable to estimate the velocity of the AIR (modulus) at the 12 points indicated in the figure I put at internet address http://www.cognitor.com.br/AirVelocity.jpg The input variables are the velocity V at the input window (pressure = 0 at the exit window)and A, B and C dimensions showed in the figure ?

Errors up to 20 % are acceptable . Can you help me indicating where to find it ? Thanks in advance

Sergio, Since your problem is a finite Volume formulation, I would recommend that you develop a code using the SIMPLE method. An Introduction to Computational Fluid Dynamics by Verstegg and Malalasekera is an excellent reference.

Good Luck Scott

Thanks for the reply. That is my problem. I have this book and others but I did not succeeded to arrive to a real code A "small code" which solve my problem will also enable me to see how to mount an algorithm using SIMPLE or SIMPLEC or other. Do you know where I can find this ?

Sergio, Pieter Wesseling has a site with some educational software. I believe there is a Matlab code there that is similar to what you want. The site is http://ta.twi.tudelft.nl/users/wesselin/ Good Luck Scott

Thanks Scott

I will download and give a feedback

Kind Regards

I opened each file and noted they are dealing with the subject but it is impossible to know which part is the relevant one for my original question which is :

I need to prepare or to find, for educational purposes a, so smal as possible, source code (Delphi, Pascal or Fortran I can understand) which enable to estimate the velocity of the AIR (modulus) at the 12 points indicated in the figure I put at internet address http://www.cognitor.com.br/AirVelocity.jpg The input variables are the velocity V at the input window (pressure = 0 at the exit window)and A, B and C dimensions showed in the figure ?

Errors up to 20 % are acceptable . Can you help me indicating where to find it ? Thanks in advance

The impressive thing is that it seems to be easy to solve but I could not find any written solution for this in the last two years.