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

dealii.demon at gmail.com dealii.demon at gmail.com
Tue Aug 13 03:48:22 PDT 2013


Revision 1832

New test.

A   trunk/aspect/tests/composition-passive-tracers/
A   trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu
A   trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu
A   trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu
A   trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu
A   trunk/aspect/tests/composition-passive-tracers/particles.pvd
A   trunk/aspect/tests/composition-passive-tracers/screen-output
A   trunk/aspect/tests/composition-passive-tracers/statistics
A   trunk/aspect/tests/composition-passive-tracers.prm


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

Diff:
Added: trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu
===================================================================
(Binary files differ)

Index: trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu	2013-08-13 10:41:38 UTC (rev 1831)
+++ trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu	2013-08-13 10:47:46 UTC (rev 1832)

Property changes on: trunk/aspect/tests/composition-passive-tracers/particles-00000.0000.vtu
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/xml
\ No newline at end of property
Added: trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu
===================================================================
(Binary files differ)

Index: trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu	2013-08-13 10:41:38 UTC (rev 1831)
+++ trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu	2013-08-13 10:47:46 UTC (rev 1832)

Property changes on: trunk/aspect/tests/composition-passive-tracers/particles-00003.0000.vtu
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/xml
\ No newline at end of property
Added: trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu
===================================================================
(Binary files differ)

Index: trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu	2013-08-13 10:41:38 UTC (rev 1831)
+++ trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu	2013-08-13 10:47:46 UTC (rev 1832)

Property changes on: trunk/aspect/tests/composition-passive-tracers/particles-00006.0000.vtu
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/xml
\ No newline at end of property
Added: trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu
===================================================================
(Binary files differ)

Index: trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu	2013-08-13 10:41:38 UTC (rev 1831)
+++ trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu	2013-08-13 10:47:46 UTC (rev 1832)

Property changes on: trunk/aspect/tests/composition-passive-tracers/particles-00009.0000.vtu
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/xml
\ No newline at end of property
Added: trunk/aspect/tests/composition-passive-tracers/particles.pvd
===================================================================
(Binary files differ)

Index: trunk/aspect/tests/composition-passive-tracers/particles.pvd
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/particles.pvd	2013-08-13 10:41:38 UTC (rev 1831)
+++ trunk/aspect/tests/composition-passive-tracers/particles.pvd	2013-08-13 10:47:46 UTC (rev 1832)

Property changes on: trunk/aspect/tests/composition-passive-tracers/particles.pvd
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/xml
\ No newline at end of property
Added: trunk/aspect/tests/composition-passive-tracers/screen-output
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/screen-output	                        (rev 0)
+++ trunk/aspect/tests/composition-passive-tracers/screen-output	2013-08-13 10:47:46 UTC (rev 1832)
@@ -0,0 +1,195 @@
+Number of active cells: 64 (on 4 levels)
+Number of degrees of freedom: 1,526 (578+81+289+289+289)
+
+*** Timestep 0:  t=0 seconds
+   Solving temperature system... 1 iterations.
+   Solving composition system 1... 0 iterations.
+   Solving composition system 2... 0 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00000
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00000
+
+*** Timestep 1:  t=0.0625 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 7 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 2:  t=0.125 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00001
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00001
+
+*** Timestep 3:  t=0.1875 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 12 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 4:  t=0.25 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00002
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00002
+
+*** Timestep 5:  t=0.3125 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00003
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00003
+
+*** Timestep 6:  t=0.375 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 17 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 7:  t=0.4375 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00004
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00004
+
+*** Timestep 8:  t=0.5 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00005
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00005
+
+*** Timestep 9:  t=0.5625 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 17 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 10:  t=0.625 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 17 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00006
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00006
+
+*** Timestep 11:  t=0.6875 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 12:  t=0.75 seconds
+   Solving temperature system... 11 iterations.
+   Solving composition system 1... 11 iterations.
+   Solving composition system 2... 11 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00007
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00007
+
+*** Timestep 13:  t=0.8125 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 4 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00008
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00008
+
+*** Timestep 14:  t=0.875 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 10 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 4 iterations.
+
+   Postprocessing:
+     Advecting particles: done
+
+*** Timestep 15:  t=0.9375 seconds
+   Solving temperature system... 10 iterations.
+   Solving composition system 1... 10 iterations.
+   Solving composition system 2... 9 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 3 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-composition-passive-tracers/solution-00009
+     Advecting particles:      done. Writing particle graphical output output-composition-passive-tracers/particles-00009
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start    |      36.8s |            |
+|                                             |            |            |
+| Section                         | no. calls |  wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system          |        16 |      3.38s |       9.2% |
+| Assemble composition system     |        32 |      13.4s |        37% |
+| Assemble temperature system     |        16 |      6.69s |        18% |
+| Build Stokes preconditioner     |        16 |      1.87s |       5.1% |
+| Build composition preconditioner|        32 |     0.049s |      0.13% |
+| Build temperature preconditioner|        16 |    0.0235s |     0.064% |
+| Solve Stokes system             |        16 |     0.174s |      0.47% |
+| Solve composition system        |        32 |    0.0463s |      0.13% |
+| Solve temperature system        |        16 |     0.022s |      0.06% |
+| Initialization                  |         2 |     0.046s |      0.12% |
+| Postprocessing                  |        16 |      10.6s |        29% |
+| Setup dof systems               |         1 |     0.106s |      0.29% |
++---------------------------------+-----------+------------+------------+
+

Added: trunk/aspect/tests/composition-passive-tracers/statistics
===================================================================
--- trunk/aspect/tests/composition-passive-tracers/statistics	                        (rev 0)
+++ trunk/aspect/tests/composition-passive-tracers/statistics	2013-08-13 10:47:46 UTC (rev 1832)
@@ -0,0 +1,28 @@
+# 1: Time step number
+# 2: Time (seconds)
+# 3: Number of mesh cells
+# 4: Number of Stokes degrees of freedom
+# 5: Number of temperature degrees of freedom
+# 6: Number of composition degrees of freedom
+# 7: Iterations for temperature solver
+# 8: Iterations for composition solver 1
+# 9: Iterations for composition solver 2
+# 10: Iterations for Stokes solver
+# 11: Time step size (seconds)
+# 12: Visualization file name
+ 0 0.0000e+00 64 659 289 289  1  0  0 18 6.2500e-02 output-composition-passive-tracers/solution-00000 
+ 1 6.2500e-02 64 659 289 289 10 10 11  7 6.2500e-02                                                "" 
+ 2 1.2500e-01 64 659 289 289 10 11 10 18 6.2500e-02 output-composition-passive-tracers/solution-00001 
+ 3 1.8750e-01 64 659 289 289 10 12 11 18 6.2500e-02                                                "" 
+ 4 2.5000e-01 64 659 289 289 10 11 11 18 6.2500e-02 output-composition-passive-tracers/solution-00002 
+ 5 3.1250e-01 64 659 289 289 10 11 11 18 6.2500e-02 output-composition-passive-tracers/solution-00003 
+ 6 3.7500e-01 64 659 289 289 10 11 10 17 6.2500e-02                                                "" 
+ 7 4.3750e-01 64 659 289 289 10 11 10 18 6.2500e-02 output-composition-passive-tracers/solution-00004 
+ 8 5.0000e-01 64 659 289 289 10 10 11 18 6.2500e-02 output-composition-passive-tracers/solution-00005 
+ 9 5.6250e-01 64 659 289 289 10 10 10 17 6.2500e-02                                                "" 
+10 6.2500e-01 64 659 289 289 10 10 10 17 6.2500e-02 output-composition-passive-tracers/solution-00006 
+11 6.8750e-01 64 659 289 289 10 10 10 18 6.2500e-02                                                "" 
+12 7.5000e-01 64 659 289 289 11 11 11 18 6.2500e-02 output-composition-passive-tracers/solution-00007 
+13 8.1250e-01 64 659 289 289 10 10 10  4 6.2500e-02 output-composition-passive-tracers/solution-00008 
+14 8.7500e-01 64 659 289 289 10 10 10  4 6.2500e-02                                                "" 
+15 9.3750e-01 64 659 289 289 10 10  9  3 6.2500e-02 output-composition-passive-tracers/solution-00009 

Added: trunk/aspect/tests/composition-passive-tracers.prm
===================================================================
--- trunk/aspect/tests/composition-passive-tracers.prm	                        (rev 0)
+++ trunk/aspect/tests/composition-passive-tracers.prm	2013-08-13 10:47:46 UTC (rev 1832)
@@ -0,0 +1,117 @@
+#########################################################
+# 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                               = 1
+set Use years in output instead of seconds = false
+set Output directory                       = output-composition-passive-tracers
+
+
+
+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          = 3
+  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 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


More information about the CIG-COMMITS mailing list