# Interfoam (OF 1.7) : pressure evolution, impact, 2D computation

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

December 16, 2010, 06:07
Interfoam (OF 1.7) : pressure evolution, impact, 2D computation
#1
New Member

Christophe Kassiotis
Join Date: Mar 2009
Location: Paris
Posts: 17
Rep Power: 10
Dear All,

I'm currently trying to perform a 2D computation for validation purpose with interFoam. The problem is quite simple : it is a modified dambreak problem with a wedge obstacle.

I'm using different meshes with around 7e3, 3e4, 1e6 and 8e6 cells. The two fluids are considered to be laminar, and the high density flow has a large viscosity (1e-2). I'm giving you the evolution of :

$\int_S p dS$

computed with libforces at the right wall (force-r.png) of the bassin and the left wall of the wedge (force-o.png). As you can see, before the impact on the right wall, the results are in accordance. The problem is at the impact ( t~0.6s ) : the pressure pic does not seems to converge. I'm even more worried by the fact that the impact on the right wall modify the pressure evolution on the left of the wedge.

This pressure problem does not seem to influence too much the following of the computation for coarse meshes, but leads to divergence for fine meshes. Moreover, I never had this problem with 3D computation.

My guess : the air captured under the tongue (t=0.4 to t = 0.7s for instance) can't escape in 2D, and has a two large pressure (see pressure at t=0.6, given in attachement). The same problem occurs with the dambreak example given as a tutorial in OpenFOAM

What do you think about it ? Do you now a way to avoid this pressure pic ?

I give you my sources as well in attachement (without the mesh that is too heavy).
Attached Images
 forces-o.png (3.6 KB, 23 views) forces-r.png (4.7 KB, 23 views) p-0400ms.png (91.2 KB, 34 views) p-0600ms.jpg (10.1 KB, 34 views)
Attached Files
 wedge-src.tar.gz (2.6 KB, 9 views)

 December 21, 2010, 10:55 #2 New Member   Christophe Kassiotis Join Date: Mar 2009 Location: Paris Posts: 17 Rep Power: 10 I tried : - change the boundary condition at the bottom right - reduce the maxCo number (to 0.05) - change the pressure solvers (pcorr, p_rgh, p_rghFinal) and add relaxation (but this does not seem to be a good idea, as I can't find p_rgh.relax() anywhere in the interFoam directory with a short grep). Nothing give for the moment satisying results. Does someone have the same problem ?

 December 21, 2010, 16:09 #3 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 Hi, sorry if I ask. You partially answered my question, but I have a doubt from your fvSolution. Did you try to use a zero relTol for p_rgh, so that the actual tolerance is achieved? If you want a correct transient solution, under-relaxing is a big no with PISO algorithm, since interFoam does not perform outer iterations ;-) Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post kwardle OpenFOAM 2 March 4, 2015 05:29 Dan Moskal Main CFD Forum 0 October 24, 2002 22:02 Abhi Main CFD Forum 12 July 8, 2002 09:11 Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 11:52 DS & HB Main CFD Forum 0 January 8, 2000 16:00