CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

venturi tube help!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 31, 2011, 18:40
Question venturi tube help!
  #1
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Dear Foamers,

I am trying to simulate the 2-D flow of air through a venturi tube (the geometry shown in the image attached). The dimensions of the tube are small; to give you an rough estimate, the outline of the mesh measures 10 mm by 20 mm. The solver I am using is sonicFoam. The starting point of my simulation was 'fowardStep' tutorial case in the sonicFoam subdirectory in tutorials.

My boundary conditions are somewhat similar to the fowardStep case that is already setup, so I changed the mesh geometry from fowardStep to the one I have, and I immediately run into trouble with I run sonicFoam: the error states the max number of iterations exceeded. How do I fix this? I tried making the time step larger but that didn't work.

Now for the real boundary conditions, I want to setup a higher pressure at the inlet, which on the left side of the image, (it's value is 27 % higher than atmospheric pressure) and atmospheric pressure at the outlet; I want to do this so I can flow air at roughly mach 1 through the tube. Does this mean that I'll have to set up the inlet pressure boundary condition to (1.27)*(101325) ? When I setup pressure like this, I get same error. I have no clue how to fix this.

Please help me in any way you can. If you need a reference, my case is just like the foward step case in the sonicFoam directory but with modified geometry and some BC conditions changed. I've read the manual but I still can't get anywhere with this problem. I can provide any information that might be necessary for you to help me.

THANKYOU!

warm regards,
Ammar.
Attached Files
File Type: pdf venturi.pdf (85.0 KB, 67 views)
atareen64 is offline   Reply With Quote

Old   February 9, 2011, 16:16
Default
  #2
New Member
 
Radko Bankras
Join Date: Jul 2010
Location: Almere, The Netherlands
Posts: 5
Rep Power: 6
Radko is on a distinguished road
Ammar,
You should decrease the time step to fix your error. Did you check your modified mesh with checkMesh? You can try to specify your inlet and outlet pressures in 0/p as:

Code:
inletPatch
{
  type fixedValue;
  value uniform (128682.75 0 0); // 27% above atm
}

outletPatch
{
  type waveTransmissive;
  field p; // name of field
  phi phi; // volumetric flux field
  rho rho; // density field
  psi psi; // density field 
  gamma 1.4; // ratio of specific heats
  fieldInf 101325; // far-field value
  lInf 0.1; // distance to far-field
  value uniform 101325;
}
You might want to calculate the critical pressure for choked flows first. I guess your ratio of 1.27 is too small and the flow will not reach Mach 1.

Kind regards, Radko

Last edited by Radko; February 9, 2011 at 16:31. Reason: Forgot code command
Radko is offline   Reply With Quote

Old   February 18, 2011, 12:58
Default
  #3
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Thank you Radko, it works now!

Regards,
Ammar.
atareen64 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
What is best mesh for cylinder tube? olivia FLUENT 7 January 16, 2014 06:55
Tube inside a tube simulation Rashmichem Main CFD Forum 0 December 14, 2010 02:52
HLL Riemann Shock Tube Matlab Problem Luke F Main CFD Forum 1 July 7, 2009 13:44
Venturi tube results make no sense Paul FLUENT 1 June 17, 2008 14:09
Flow of tracer in a tube uma FLUENT 9 October 3, 2003 14:53


All times are GMT -4. The time now is 08:37.