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

rhoSimpleFoam - test case

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 29, 2011, 12:20
Angry rhoSimpleFoam - test case
  #1
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Hello Everyone!

I posted this question earlier but got no replies, so I am posting this again in hopes of getting some help or insights. I have a rectangular pipe, with an inlet on one side and an outlet on the other side, higher pressure on the inlet and lower pressure on the outlet, so the fluid should flow from the inlet to the outlet.

The solver I wanna use is rhoSimpleFoam. I got some really good results with sonicFoam but I want to use a steady-state solver. Right now I am getting some reasonable profiles but the solution diverges soon, chiefly because I think I am not setting up the relaxation factors in fvSolution correctly.

Kindly look at my case and tell me what I could be possibly doing wrong. This is a test case and should be very easy to figure out (ironically). Please please help.

Case file
https://rapidshare.com/files/4597932...SimpleFoam.zip

Regards,
~Ammar.
atareen64 is offline   Reply With Quote

Old   April 29, 2011, 13:22
Default
  #2
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
Rename your inlet to inlet and outlet to outlet and I hope you will notice one mistake yourself.
This will be noticable once you turn on turbulence though.

The main problem is the outlet bc for temperature.
set it to zeroGradient.
you cant have fixedValue temperature on both inlet and outlet.
niklas is offline   Reply With Quote

Old   April 29, 2011, 13:27
Default
  #3
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Thanks Niklas, doing it right now.
atareen64 is offline   Reply With Quote

Old   April 29, 2011, 14:15
Default
  #4
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Thank you Niklas, the results look a lot a lot better now.

Quick question: why can't I specify T at both the inlet and outlet? does that over-constrain the problem?

Also I changed the inlet and outlet names to inlet and outlet. The only thing I noticed was that my 'k' outlet condition seemed suspicious. The simulation still ran but I started getting 'nan' after a while. Could you give me a hint as to what I might be doing wrong in the turbulent case?

If not I'll look for other boundary conditions. Thanks A LOT for you're help!

regards,
Ammar.
atareen64 is offline   Reply With Quote

Old   April 29, 2011, 16:16
Default
  #5
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
Quote:
Originally Posted by atareen64 View Post
Thank you Niklas, the results look a lot a lot better now.

Quick question: why can't I specify T at both the inlet and outlet? does that over-constrain the problem?
yes.

Quote:
Originally Posted by atareen64 View Post
Also I changed the inlet and outlet names to inlet and outlet. The only thing I noticed was that my 'k' outlet condition seemed suspicious. The simulation still ran but I started getting 'nan' after a while. Could you give me a hint as to what I might be doing wrong in the turbulent case?
you have the inlet bc on outlet and vice versa. just switch them around (on k and epsilon)
niklas is offline   Reply With Quote

Old   May 1, 2011, 17:17
Smile
  #6
Member
 
Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 6
atareen64 is on a distinguished road
Niklas,

Thank you so much for your help. I see the mistakes I was making and I have corrected them. I really appreciate your help! One last question:

In my laminar case, I set the end time to be 3000. However at 3000, the velocity and pressure profiles look really bizarre, however at around t = 1000, the profiles look very good and comparable to a simulation carried in solidWorks. I find this very strange and can't seem to think why this is happening. Should I just let it run for much longer? Do you have any ideas why this might me happening?

Thanks again for your help!
warm 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
RhoSimpleFoam tutorial test case system dictionaries flavio_pergolesi OpenFOAM Running, Solving & CFD 2 April 17, 2014 08:41
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 14 March 8, 2011 09:19
Help Required With Simple Test Case steph79 OpenFOAM Pre-Processing 4 August 3, 2010 07:45
3D TRANSITION TEST CASE venkatesh4386@gmail.com FLUENT 0 March 9, 2009 13:04
Porous Media test case Alex FLUENT 0 April 9, 2006 08:23


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