[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit 7403d1c3bd6b9a987891330e668320d1a23eabeb
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Thu Feb 2 08:06:42 2006 +1100

    Stability condition for surface tension was still too restrictive
    
    darcs-hash:20060201210642-d4795-73569ecf32d9032ce9ddc2d126c4bbf0f45d9a4d.gz

diff --git a/src/tension.c b/src/tension.c
index 953f1cb..fade65a 100644
--- a/src/tension.c
+++ b/src/tension.c
@@ -233,14 +233,14 @@ static gdouble gfs_source_tension_stability (GfsSourceGeneric * s,
   p.c = t->c;
   gfs_domain_cell_traverse (GFS_DOMAIN (sim), FTT_PRE_ORDER, FTT_TRAVERSE_LEAFS, -1,
 			    (FttCellTraverseFunc) min_max_alpha, &p);
-  h = 2.*ftt_level_size (p.depth);
+  h = ftt_level_size (p.depth);
   if (p.alpha) {
     gdouble rhom = (1./p.amin + 1./p.amax)/2.;
     gdouble rhod = p.amin < p.amax ? p.amin/p.amax : p.amax/p.amin;
-    return sqrt (rhom*rhod*h*h*h/(2.*M_PI*sigma));
+    return 3.*sqrt (rhom*rhod*h*h*h/sigma);
   }
   else
-    return sqrt (h*h*h/(2.*M_PI*sigma));
+    return 3.*sqrt (h*h*h/sigma);
 }
 
 static void gfs_source_tension_class_init (GfsSourceGenericClass * klass)

-- 
Gerris Flow Solver