[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
popinet at users.sf.net
Fri May 15 02:54:30 UTC 2009
The following commit has been merged in the upstream branch:
commit e04e50f18179ae11a1fd4168ae637b50a78d5d41
Author: Stephane Popinet <popinet at users.sf.net>
Date: Wed Aug 22 13:26:02 2007 +1000
Updated 'shear' VOF advection test case
Including test for mass conservation.
darcs-hash:20070822032602-d4795-48eef891e90ea5998ff7a534bcbcd292a324f140.gz
diff --git a/test/shear/norms.ref b/test/shear/norms.ref
index 81aa016..7f9c6ea 100644
--- a/test/shear/norms.ref
+++ b/test/shear/norms.ref
@@ -1 +1 @@
-T time: 5 first: 2.553e-03 second: 4.056e-02 infty: 1.000e+00 bias: -1.594e-03
+T time: 5 first: 1.662e-04 second: 5.739e-03 infty: 3.802e-01 bias: 7.507e-06
diff --git a/test/shear/shear.gfs b/test/shear/shear.gfs
index 024b3cc..61e4a9e 100644
--- a/test/shear/shear.gfs
+++ b/test/shear/shear.gfs
@@ -60,14 +60,14 @@
# Add tracer T, using a VOF advection scheme.
# The default scheme is a Van-Leer limited, second-order upwind scheme.
- VariableTracerVOF {} T
+ VariableTracerVOF T
- InitFraction {} T (ellipse (0, -.236338, 0.2, 0.2))
+ InitFraction T (ellipse (0, -.236338, 0.2, 0.2))
# Initialize U and V with the vortical shear flow field
Init {} {
- U = sin((x + 0.5)*M_PI)*cos((y + 0.5)*M_PI)
- V = -cos((x + 0.5)*M_PI)*sin((y + 0.5)*M_PI)
+ U = sin((x + 0.5)*M_PI)*cos((y + 0.5)*M_PI)
+ V = -cos((x + 0.5)*M_PI)*sin((y + 0.5)*M_PI)
}
# At t = 2.5 re-initialize U and V with the reversed flow field
@@ -82,7 +82,7 @@
OutputPPM { start = 5 } { convert ppm:- t-5.eps } { v = T }
# Add a new variable
- Variable {} Tref
+ Variable Tref
# Initialize Tref with the initial shape
InitFraction { start = end } Tref (ellipse (0, -.236338, 0.2, 0.2))
@@ -90,7 +90,7 @@
# Output the norms of the difference between T and Tref, stores that into
# new variable DT
OutputErrorNorm { start = end } norms { v = T } {
- s = Tref v = DT
+ s = Tref v = DT
}
OutputPPM { start = end } { convert ppm:- dt-5.eps } { v = DT }
diff --git a/test/shear/shear.sh b/test/shear/shear.sh
index 5ced3fb..6fbac65 100644
--- a/test/shear/shear.sh
+++ b/test/shear/shear.sh
@@ -16,7 +16,8 @@ from check import *
from sys import *
if (Curve('norms',3,5) - Curve('norms.ref',3,5)).max() > 0. or\
(Curve('norms',3,7) - Curve('norms.ref',3,7)).max() > 0. or\
- (Curve('norms',3,9) - Curve('norms.ref',3,9)).max() > 0.:
+ (Curve('norms',3,9) - Curve('norms.ref',3,9)).max() > 0. or\
+ (Curve('norms',3,11) - Curve('norms.ref',3,11)).max() > 0.:
exit(1)
EOF
else
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list