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

Pressure divergence with interTrackFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 2, 2009, 09:11
Default Hello, I am using interTra
  #1
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Hello,

I am using interTrackFoam to simulate the flow of a fluid which first slides along a slope and then falls freely in a basin. The simulation runs fine until a while. The pressure of the fluid at the point where the slope stops and the fluid falls freely begins to diverge. Could someone give me some advice about how I could fix this problem?

I am using laplaceFaceDecomposition solver for the mesh motion and I did not change the coefficients for the solving of U and p from the hydrofoil tutorial.

If you need some more information, just tell me.

Thank you a lot in advance

Virginie
virginie_e is offline   Reply With Quote

Old   March 2, 2009, 15:04
Default Could you post a few pictures?
  #2
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
Could you post a few pictures? Preferably wire-frame images.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   March 3, 2009, 04:46
Default Hello Sandeep, Here are a
  #3
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Hello Sandeep,

Here are a few pictures. There are the pressure plot, the velocity plot and a zoom on the part where the pressure diverges.
The inlet is the upper left corner, the outlet is the lower right corner. The fluid flows down the "spout" then falls freely over a bath of liquid tin. I simulated the tin with a boundary pressure condition over the free surface.














As you can see the pressure diverges right at the point where the fluid leaves the slope and falls freely. And there does not seem to be any mesh problems… They only come after a while. I wonder if I choose a wrong pressure solver… Could you please give me some advice about it? Thanks a lot.

Virginie
virginie_e is offline   Reply With Quote

Old   March 3, 2009, 05:39
Default Here is my case directory if y
  #4
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Here is my case directory if you want to have a look.

Here is what I do:
fluentMeshToFoam tinpressure.mah
makeFaMesh
interTrackFoam



Thank you for your help!

Virginie
virginie_e is offline   Reply With Quote

Old   March 3, 2009, 05:49
Default Oups, sorry I did not atta
  #5
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Oups,

sorry I did not attach the files correctly

tinpressurecopy.zip

tinpressure.msh.gz
virginie_e is offline   Reply With Quote

Old   March 3, 2009, 05:58
Default Sorry I also forgot my library
  #6
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Sorry I also forgot my library for the pressure condition on the free surface (p=0 above the tin level, otherwise p = rho*g*z basically)

linearTinPressurePerso.zip
virginie_e is offline   Reply With Quote

Old   March 3, 2009, 15:53
Default I only see an entry for SIMPLE
  #7
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
I only see an entry for SIMPLE in fvSolution - is this a steady-state type of simulation? You might want to have a time step restriction using the Courant number too. For bad meshes, you might want to increase the number of nonOrthogonal correctors.

Also, the mesh looks like it was obtained from a previous calculation. Surface tension might behave erratically if the surface normals are ill-behaved.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   March 4, 2009, 04:20
Default Thank you Sandeep for your qui
  #8
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Thank you Sandeep for your quick answer.

I have alreday tried to increase the number of nonOrthogonal correctors but I had the same divergence.
I will try though the time step restriction using the Courant number. That is a good idea.

Yes it is indeed a steady-state type of simulation. Is there a kind of option which I should set? If yes, what kind of option?

Anyway, thanks for your help.
virginie_e is offline   Reply With Quote

Old   March 4, 2009, 06:07
Default Hello, I have tried to set
  #9
Member
 
Virginie Ehrlacher
Join Date: Mar 2009
Posts: 52
Rep Power: 8
virginie_e is on a distinguished road
Hello,

I have tried to set a time step restriction using the Courant number and to increase teh number of nonOrthogonal correctors (30) but I still get the same kind of error. Have you got some other ideas about what I should do?

Otherwise, do you know where I could find some documentation about all the available options that are implemented for SIMPLE? I have tried to search for it but I could not find anything.

Thank you for your time.
virginie_e 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
InterTrackFoam any information rajon OpenFOAM Running, Solving & CFD 29 January 17, 2012 16:50
Wall contact in interTrackFoam virginie_e OpenFOAM Running, Solving & CFD 2 November 8, 2011 12:52
SIMPLE loop in interTrackFoam virginie_e OpenFOAM Running, Solving & CFD 3 March 17, 2009 06:40
Divergence of trailing edge pressure? Bren Main CFD Forum 3 December 6, 2007 18:57
InterTrackFoam error kester OpenFOAM Running, Solving & CFD 10 November 8, 2007 03:55


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