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

Convergence issues in extremely low Re problems

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 28, 2014, 11:53
Default Convergence issues in extremely low Re problems
  #1
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 452
Rep Power: 23
santiagomarquezd will become famous soon enough
Hi all! I've been working in a extremely low Reynolds problems, namely Re = 1.16x10^−6. The problem is very simple a 2D water drop in a box filled with petroleum, the parameters are as follows:

rho_w=1000
nu_w=1e-6
rho_p=800
nu_p=3.75e-4
sigma=0.02
g=10
b=600e-6 (box' sides)
dw=100e-6 (water droplets)
dx=10E-6 (structured mesh)

the idea is to predict the sedimentation velocity which is v=3.7x10^-6. (More details and comparisons with other codes here: http://cimec.org.ar/ojs/index.php/mc...view/4648/4580). I've found the following limitations:

1. The curvature calculation procedure is basically useless, the errors in curvature values are greater than 50% (while in Gerris, f.e. they are below 1%). Due to this, I changed the surface tension to sigma=0.00002, since this change doesn't modify the sedimentation velocity but dumps the effect of the surface tension term. Without doing this the parasitic current are so high that the problem cannot be solved.

2. With these new case (low sigma) the solver still has problems, now with convergence, it is extremely difficult to go down in residual values (hundreds of PISO steps due to many outer PIMPLE loops).

3. Even converging to residuals O(1E-13) the final predicted velocity is high respect to the theoretical value, approximately two times.

4. I programmed a interStokesFoam, which doesn't have the convection term div(rho U U) and the velocity prediction was improved, but not at all.

So that, some questions:

a. Why these kind of problem are so hard to solve? I've found that Gerris (Fractional Step p-V coupling) or other codes with monolitic p-V coupling perform really much better (see paper ref)
b. What kind of settings and/or solver changes could improve the quality and speed of the solution (actual setting are similar to damBreak problem but with a lot of PIMPLE, PISO loops and exiting by converge criteria).

Thanks in advance!
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd 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
wet steam model convergence problems adeln FLUENT 1 July 8, 2013 18:49
No convergence method available in ANYS FSI problems? acasas Structural Mechanics 11 May 12, 2013 15:39
force convergence problems in CFX 6DOF rigid body solver ajay_ks CFX 8 March 25, 2013 05:02
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 06:59
convergence problems david aquilina Siemens 4 October 27, 2000 10:18


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