CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

PimpleFoam problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 28, 2010, 05:40
Default PimpleFoam problem
  #1
New Member
 
Vadims Geza
Join Date: Apr 2010
Location: Latvia
Posts: 9
Rep Power: 16
Vadims.geza is on a distinguished road
Hi all!

Maybe a dumb question, but I cannot run pimpleFoam. Seems to be something general.
I have done simpleFoam with SST, everything was fine. now I want to run pimpleFoam with same model, but...

[3] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #2 __restore_rt in "/lib64/tls/libc.so.6"
[3] #3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/geza/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"

Have tried different BC's for k and omega. Have also changed schemes. Nothing helps. Every time I start pimple, it goes to large Co and crashes.
My geometry actually is just a domain with walls, inlet and outlet. If I take simple tube with same BC's, everything is fine.
Have also tried different meshes - all of them works fine with simple, but none of them with pimple.
Have also done LES with piso, and it was also fine. When I use mesh and timestep from LES in pimple with SST, crash again.
Please give me some advice!


Thanks in advance!
Vadims.geza is offline   Reply With Quote

Old   July 31, 2010, 06:31
Default
  #2
New Member
 
Robert
Join Date: Apr 2010
Posts: 16
Rep Power: 16
rob3rt is on a distinguished road
Dear Vadims,

I am having a similar problem like you had.

How did you resolve it?

Kind Regards,
Robert
rob3rt is offline   Reply With Quote

Old   December 6, 2011, 05:33
Default
  #3
New Member
 
rene
Join Date: Oct 2011
Posts: 7
Rep Power: 15
krene is on a distinguished road
Dito. I do also have this problems.
krene is offline   Reply With Quote

Old   December 6, 2011, 06:51
Default
  #4
New Member
 
GregorS's Avatar
 
Gregor Seljak
Join Date: Oct 2010
Posts: 21
Rep Power: 16
GregorS is on a distinguished road
Did you change your ddtSchemes from steadyState (used with simpleFoam) to Euler ?
GregorS is offline   Reply With Quote

Old   December 6, 2011, 09:50
Default
  #5
New Member
 
rene
Join Date: Oct 2011
Posts: 7
Rep Power: 15
krene is on a distinguished road
Yes, I tried different ddtSchemes with pimple.

Well I get some results, but this results are completly different to the results out of simpleFoam...
krene is offline   Reply With Quote

Old   December 6, 2011, 12:03
Default
  #6
Senior Member
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0
chegdan will become famous soon enoughchegdan will become famous soon enough
Quote:
Originally Posted by krene View Post
Yes, I tried different ddtSchemes with pimple.

Well I get some results, but this results are completly different to the results out of simpleFoam...
@krene

If you are getting results from a transient solver (pimpleFoam) that differ from your steady state solution from simpleFoam, it might be that you are in a transitional state and not fully turbulent. Does your turbulence dissipate? Is it dissipating from your estimated inlet values of your turbulent parameters? I will add that I'm not an LES expert, but how are you comparing the RANS results from simpleFOAM with the LES of your simulations? Also, RANS simulations are known to be dissipative with respect to turbulence (I have no source on this one, only from my notes at the AIChE conference from the turbulence sessions). Just my two cents. Good luck.
chegdan is offline   Reply With Quote

Old   December 6, 2011, 15:43
Default
  #7
New Member
 
rene
Join Date: Oct 2011
Posts: 7
Rep Power: 15
krene is on a distinguished road
Thanks to Gregor and Dan for your responses.

First I'll introduce you my case, how my case looks like, and afterwards where my problems are.

My case is a cylinder with a diameter of D = 300mm and a length L = 2000mm. There is a inlet 160mm from the top of the tank with a diameter of Di = 3mm. The velocity is assumed to be 1.67m/s which leads to a Reynoldsnumber of around 5000. The fluid is water. At the bottom there is an outlet with a diameter of Do = 25mm. While running the experiment, the pressure at the outlet is constant.

I'm interested to simulate the jet at the inlet. The goal is to estimate the entrainment rate.

I'd like to use a transient solver like PISO or PIMPLE to simulate the transient case. I could simulate the case with SIMPLE, but I'm not so happy with the results.
As BC I do have for U:
Inlet: fixedValue
Outlet: zeroGradient
Wall: fixedValue (0 0 0) ==> with zeroGradient, SIMPLE will work for 400s

BC for p:
Inlet: zeroGradient
Outlet: fixedValue, 0
Wall: fixedValue, 0

With this BC, the results looks OK, but there is also the unlogical behaviour, that there is no backflow near the walls.

If I take this BC to the PIMPLE solver, the jet does not look like a jet, it's more a pencil, which goes through a tart. There is no entrainment and there is also no interaction at the boarder between the turbulent and the laminar layer. The same behaviour can be seen by increasing the velocity to 3.33m/s (Re = 10'000).

Quote:
Originally Posted by chegdan View Post
I will add that I'm not an LES expert, but how are you comparing the RANS results from simpleFOAM with the LES of your simulations?
The results are compared visually, the errors are obviously. Well, I'm not an expert in modelling with LES and RANS. It's the first time I'm working with OpenFOAM, so I do not have to much experiences.

If anyone has time to have a look at my case, I could also upload it here. Thanks a lot.
krene is offline   Reply With Quote

Old   December 6, 2011, 17:19
Default
  #8
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
At the walls, BCs for velocity and pressure should be:

U--->fixedValue (0 0 0)

p--->zeroGradient (not fixedValue!!!)

Regards

V.
vkrastev is offline   Reply With Quote

Old   December 7, 2011, 09:54
Default
  #9
New Member
 
rene
Join Date: Oct 2011
Posts: 7
Rep Power: 15
krene is on a distinguished road
Quote:
Originally Posted by vkrastev View Post
p--->zeroGradient (not fixedValue!!!)
Hi Vesselin,
thanks for your replay.
You're right, p should be zeroGradient.

After setup the first model, I got some bad results, so I adjusted the BC so, that I got the wanted results. And the results were good with fixedValue for p at the wall.

I made a run with SIMPLE, PISO and PIMPLE.

The results out of SIMPLE (RANS) are reasonable, also there is the expected backflow in the tank.

With PISO and PIMPLE (LES), the results are still not useable. Same behaviour as before, there is no jet.

Thank you guys so far for your help!
krene is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 05:43
Incoherent problem table in hollow-fiber spinning Gianni FLUENT 0 April 5, 2008 11:33
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 15:52


All times are GMT -4. The time now is 15:00.