New Code Saturne tutorial online
- Claus
-
Topic Author
- Offline
- Moderator
-
- Posts: 670
- Thank you received: 34
www.caelinux.org/wiki/index.php/Contrib:...ld_study_CodeSaturne
It is not a complete tutorial, more of a walkthrough of what I did while still sharing the details of the setup.
I would really like some feedback on it - whether the results are correct, what I can do to optimize the work process and the simulation etc.
Otherwise enjoy

/C
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
- Robin Borland
- Offline
- Senior Member
-
- Posts: 48
- Thank you received: 0

Thanks for posting this - it is always good to see some verification models done.
- Robin Borland
- Offline
- Senior Member
-
- Posts: 48
- Thank you received: 0
A couple of specific comments as follows:
1. Using the orifice flow equation from Crane Technical Paper 410, with your orifice size, 100 Pa upstream, and assuming 13mm ID upstream gives a volumetric flowrate of 9.3 *10^-5 M^3/s - the average velocity in the orifice comes out to 32.8 M/s
2. If I model the orifice with an upstream and downstream section and put a velocity in of 0.701 - I get a parabolic type flow profile in the orifice which peaks around 44 M/s - see the attached image
3. The pressure in the outlet section is in the -25 Pa range - however the upstream pressure back calculated to around 750 Pa ( not the 1000 Pa expected) - note these are not exact numbers as I am coming off of the scalar maps
Post edited by: Robin Borland, at: 2009/04/21 20:51<br /><br />Post edited by: Robin Borland, at: 2009/04/21 20:52
- Robin Borland
- Offline
- Senior Member
-
- Posts: 48
- Thank you received: 0
Claus,
A couple of specific comments as follows:
1. Using the orifice flow equation from Crane Technical Paper 410, with your orifice size, 100 Pa upstream, and assuming 13mm ID upstream gives a volumetric flowrate of 9.3 *10^-5 M^3/s - the average velocity in the orifice comes out to 32.8 M/s
2. If I model the orifice with an upstream and downstream section and put a velocity in of 0.701 - I get a parabolic type flow profile in the orifice which peaks around 44 M/s - see the attached image
3. The pressure in the outlet section is in the -25 Pa range - however the upstream pressure back calculated to around 750 Pa ( not the 1000 Pa expected) - note these are not exact numbers as I am coming off of the scalar maps
Post edited by: Robin Borland, at: 2009/04/21 20:51<br /><br />Post edited by: Robin Borland, at: 2009/04/21 20:52
- Claus
-
Topic Author
- Offline
- Moderator
-
- Posts: 670
- Thank you received: 34
Claus,
A couple of specific comments as follows:
1. Using the orifice flow equation from Crane Technical Paper 410, with your orifice size, 100 Pa upstream, and assuming 13mm ID upstream gives a volumetric flowrate of 9.3 *10^-5 M^3/s - the average velocity in the orifice comes out to 32.8 M/s
2. If I model the orifice with an upstream and downstream section and put a velocity in of 0.701 - I get a parabolic type flow profile in the orifice which peaks around 44 M/s - see the attached image
3. The pressure in the outlet section is in the -25 Pa range - however the upstream pressure back calculated to around 750 Pa ( not the 1000 Pa expected) - note these are not exact numbers as I am coming off of the scalar maps
Post edited by: Robin Borland, at: 2009/04/21 20:51<br /><br />Post edited by: Robin Borland, at: 2009/04/21 20:52
Thank you for your comments:
#1 I need to get a better text book for my calculations - the one I have is lacking at best. What does your text book say about the actual case of nozzle size 1.9, upstream pressure of .7kPa downstream dia of 13mm?
#2/#3 Well that's what I expected - if you look at the picture with the manometer, you can see that the gas has to travel through a series of bends and a solenoid valve - hence the pressure should be lower than 1kPa - I just didn't calculate the pressure drop - I lowered the orifice exit speed - but if you say it's 0.7kPa, I should try running the calcs. with that.
#2 (again) inlet velocity of 0.701 ?
thx for taking the time out to test it

Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
- Robin Borland
- Offline
- Senior Member
-
- Posts: 48
- Thank you received: 0
I have attached my Mathcad file for this - I think my flowrate comes out to be 5.6 l/min where as your flowmeter shows 6 - so the velocity that I used is probably low, if I raised it, the inlet pressure would come up a bit also.
The velocity of 0.701 is the hypothetical velocity used on an inlet section 13mm diameter - see V2 in the calculation. Of course this is based on zero compressibility which we know is not true, but an ok assumption for now. I have tried to use the compressibility feature on some other work and struggled with it - I gave up for now, but will circle back around to it.
If you are going to be doing a lot of tests with orifices, valves etc -the Crane 410 paper is a good investment - it is only $30 or $40 if you order it on line. The calculations in this book will get you in the ball-park.
Let me know if you want my HDF file for this.
Best regards
Robin