[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
s.popinet at niwa.co.nz
Fri May 15 02:52:12 UTC 2009
The following commit has been merged in the upstream branch:
commit 55fdb2b7187b9b007c894afc8d6caf1daf43158f
Author: Stephane Popinet <s.popinet at niwa.co.nz>
Date: Wed Jul 6 07:52:21 2005 +1000
Changed BC for boussinesq example
darcs-hash:20050705215221-fbd8f-1ce2c2378a0e3e87263cf9fd0b8f4ce770f8df70.gz
diff --git a/doc/examples/boussinesq/boussinesq.gfs b/doc/examples/boussinesq/boussinesq.gfs
index a9491c1..867d40c 100644
--- a/doc/examples/boussinesq/boussinesq.gfs
+++ b/doc/examples/boussinesq/boussinesq.gfs
@@ -30,7 +30,7 @@
# Running time: 60 minutes
# Generated files: t.mpg t.eps
#
-2 1 GfsSimulation GfsBox GfsGEdge {} {
+3 2 GfsSimulation GfsBox GfsGEdge {} {
Time { end = 20 }
# Use an initial refinement of 8 levels around the solid boundary
@@ -52,11 +52,11 @@
SurfaceBc T Dirichlet 1
# Adapt the mesh using the vorticity criterion at every timestep
- # down to a maximum level of 8 if y is smaller than 1.25, 3 otherwise.
- # The topmost part of the domain will never have more than
- # 3 levels of refinement and will act as a very efficient "sponge"
- # layer to damp any eddies before they exit the domain.
- AdaptVorticity { istep = 1 } { maxlevel = (y > 1.25 ? 3 : 8) cmax = 1e-2 }
+ # down to a maximum level of 8 if y is smaller than 1.5, 0
+ # otherwise. The topmost part of the domain will not be refined and
+ # will act as a very efficient "sponge" layer to damp any eddies
+ # before they exit the domain.
+ AdaptVorticity { istep = 1 } { maxlevel = (y > 1.5 ? 0 : 8) cmax = 1e-2 }
# Also adapt according to the tracer gradient
AdaptGradient { istep = 1 } { maxlevel = 8 cmax = 5e-2 } T
@@ -90,8 +90,10 @@
EventScript { start = end } { convert -colors 256 t.ppm t.eps ; rm -f t.ppm }
}
GfsBox {}
+GfsBox {}
# The top boundary is a simple outflow condition. This could cause problems
# (eddies getting stuck on the boundary) if the adaptive "sponge" layer was not
# used.
GfsBox { top = BoundaryOutflow }
1 2 top
+2 3 top
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list