[cig-commits] commit by bangerth to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1835

Add first tracer cookbook.

A   trunk/aspect/cookbooks/composition-passive-tracers.prm


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1835&peg=1835

Diff:
Copied: trunk/aspect/cookbooks/composition-passive-tracers.prm (from rev 1821, trunk/aspect/cookbooks/composition-passive.prm)
===================================================================
--- trunk/aspect/cookbooks/composition-passive-tracers.prm	                        (rev 0)
+++ trunk/aspect/cookbooks/composition-passive-tracers.prm	2013-08-13 11:27:27 UTC (rev 1835)
@@ -0,0 +1,118 @@
+#########################################################
+# This is a variation of the composition-passive.prm
+# parameter file, with the exception that we also
+# advect along a set of tracers and that we do not longer
+# compute temperature statistics and composition statistics
+# as part of postprocessing
+
+set Dimension                              = 2
+set Start time                             = 0
+set End time                               = 20
+set Use years in output instead of seconds = false
+set Output directory                       = output
+
+
+
+subsection Geometry model
+  set Model name = box
+
+  subsection Box
+    set X extent = 2
+    set Y extent = 1
+  end
+end
+
+
+subsection Model settings
+  set Fixed temperature boundary indicators   = 2, 3
+  set Zero velocity boundary indicators       =
+  set Tangential velocity boundary indicators = 0, 1, 2
+  set Prescribed velocity boundary indicators = 3: function
+end
+
+
+subsection Boundary temperature model
+  set Model name = box
+
+  subsection Box
+    set Bottom temperature = 1
+    set Top temperature    = 0
+  end
+end
+
+
+subsection Boundary velocity model
+  subsection Function
+    set Variable names      = x,z,t
+    set Function constants  = pi=3.1415926
+    set Function expression = if(x>1+sin(0.5*pi*t), 1, -1); 0
+  end
+end
+
+
+subsection Gravity model
+  set Model name = vertical
+end
+
+
+subsection Initial conditions
+  set Model name = function
+
+  subsection Function
+    set Variable names      = x,z
+    set Function expression = (1-z)
+  end
+end
+
+
+subsection Material model
+  set Model name = simple
+
+  subsection Simple model
+    set Thermal conductivity          = 1e-6
+    set Thermal expansion coefficient = 1e-4
+    set Viscosity                     = 1
+  end
+end
+
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 0
+  set Initial global refinement          = 5
+  set Time steps between mesh refinement = 0
+end
+
+
+subsection Postprocess
+  set List of postprocessors = visualization, tracers
+
+  subsection Visualization
+    set Time between graphical output = 0.1
+  end
+
+  subsection Tracers
+    set Number of tracers        = 1000
+    set Time between data output = 0.1
+    set Data output format       = vtu
+  end
+end
+
+
+# This is the new part: We declare that there will
+# be two compositional fields that will be
+# advected along. Their initial conditions are given by
+# a function that is one for the lowermost 0.2 height
+# units of the domain and zero otherwise in the first case,
+# and one in the top most 0.2 height units in the latter.
+subsection Compositional fields
+  set Number of fields = 2
+end
+
+subsection Compositional initial conditions
+  set Model name = function
+
+  subsection Function
+    set Variable names      = x,y
+    set Function expression = if(y<0.2, 1, 0) ; if(y>0.8, 1, 0)
+  end
+end