[cig-commits] [commit] master: Add a test for Anne Glerum's newly written spherical velocity statistics postprocessor. (22f7bf3)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue May 20 11:54:13 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/8ede49902e2093f35fe062c7fc0f2339a96b9b57...22f7bf397e02ab73095f6873dd2461c456aa8031
>---------------------------------------------------------------
commit 22f7bf397e02ab73095f6873dd2461c456aa8031
Author: Wolfgang Bangerth <bangerth at math.tamu.edu>
Date: Tue May 20 13:51:46 2014 -0500
Add a test for Anne Glerum's newly written spherical velocity statistics postprocessor.
>---------------------------------------------------------------
22f7bf397e02ab73095f6873dd2461c456aa8031
...e_vel.prm => spherical_velocity_statistics.prm} | 5 +-
tests/spherical_velocity_statistics/screen-output | 206 +++++++++++++++++++++
tests/spherical_velocity_statistics/statistics | 32 ++++
3 files changed, 242 insertions(+), 1 deletion(-)
diff --git a/tests/refine_vel.prm b/tests/spherical_velocity_statistics.prm
similarity index 93%
copy from tests/refine_vel.prm
copy to tests/spherical_velocity_statistics.prm
index abb452f..5543b26 100644
--- a/tests/refine_vel.prm
+++ b/tests/spherical_velocity_statistics.prm
@@ -1,3 +1,6 @@
+# like refine_vel.prm, but run the "spherical velocity statistics"
+# postprocessor
+
set Dimension = 2
set CFL number = 1.0
set End time = 1e300
@@ -98,7 +101,7 @@ subsection Model settings
end
subsection Postprocess
- set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average
+ set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average, spherical velocity statistics
subsection Visualization
set Time between graphical output = 0
diff --git a/tests/spherical_velocity_statistics/screen-output b/tests/spherical_velocity_statistics/screen-output
new file mode 100644
index 0000000..62a706b
--- /dev/null
+++ b/tests/spherical_velocity_statistics/screen-output
@@ -0,0 +1,206 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+-- . running in OPTIMIZED mode
+-- . running with 1 MPI process
+-- . using Trilinos
+-----------------------------------------------------------------------------
+
+Number of active cells: 768 (on 5 levels)
+Number of degrees of freedom: 10,436 (6,402+833+3,201)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+3 iterations.
+
+ Postprocessing:
+
+ Reference density (kg/m^3): 3300
+ Reference gravity (m/s^2): 9.81
+ Reference thermal expansion (1/K): 2e-05
+ Temperature contrast across model domain (K): 0
+ Model domain depth (m): 2.89e+06
+ Reference thermal diffusivity (m^2/s): 2.42424e-13
+ Reference viscosity (Pas): 1e+22
+ Ra number: 0
+ k_value: 1e-06
+ reference_cp: 1250
+ reference_thermal_diffusivity: 2.42424e-13
+
+ Writing graphical output: output-spherical_velocity_statistics/solution-00000
+ RMS, max velocity: 1.21e-09 m/s, 3.64e-09 m/s
+ Temperature min/avg/max: 0 K, 101.4 K, 1613 K
+ Heat fluxes through boundary parts: -0.129 W, -0.2469 W, -4.669e-06 W, -4.669e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 8.09e-10 m/s, 8.93e-10 m/s, 1.21e-09 m/s
+
+Number of active cells: 837 (on 6 levels)
+Number of degrees of freedom: 11,660 (7,154+929+3,577)
+
+Number of active cells: 912 (on 6 levels)
+Number of degrees of freedom: 12,904 (7,916+1,030+3,958)
+
+*** Timestep 1: t=1.56836e+13 seconds
+ Solving temperature system... 19 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 24 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00001
+ RMS, max velocity: 1.26e-09 m/s, 3.75e-09 m/s
+ Temperature min/avg/max: -202.2 K, 102.8 K, 1613 K
+ Heat fluxes through boundary parts: -0.1256 W, -0.2457 W, -2.204e-06 W, -6.235e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 8.54e-10 m/s, 9.3e-10 m/s, 1.26e-09 m/s
+
+*** Timestep 2: t=2.32729e+13 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00002
+ RMS, max velocity: 1.29e-09 m/s, 3.8e-09 m/s
+ Temperature min/avg/max: -202.4 K, 103.8 K, 1613 K
+ Heat fluxes through boundary parts: -0.1229 W, -0.2444 W, -1.313e-06 W, -8.754e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 8.76e-10 m/s, 9.49e-10 m/s, 1.29e-09 m/s
+
+Number of active cells: 990 (on 6 levels)
+Number of degrees of freedom: 14,115 (8,658+1,128+4,329)
+
+*** Timestep 3: t=3.07665e+13 seconds
+ Solving temperature system... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00003
+ RMS, max velocity: 1.32e-09 m/s, 3.85e-09 m/s
+ Temperature min/avg/max: -202.6 K, 104.7 K, 1613 K
+ Heat fluxes through boundary parts: -0.1202 W, -0.2435 W, -3.309e-07 W, -1.015e-05 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 8.99e-10 m/s, 9.68e-10 m/s, 1.32e-09 m/s
+
+*** Timestep 4: t=3.81707e+13 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 19 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00004
+ RMS, max velocity: 1.35e-09 m/s, 3.89e-09 m/s
+ Temperature min/avg/max: -202.9 K, 105.6 K, 1613 K
+ Heat fluxes through boundary parts: -0.1178 W, -0.2424 W, 1.713e-07 W, -1.151e-05 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 9.22e-10 m/s, 9.86e-10 m/s, 1.35e-09 m/s
+
+Number of active cells: 1,092 (on 6 levels)
+Number of degrees of freedom: 15,598 (9,568+1,246+4,784)
+
+*** Timestep 5: t=4.54881e+13 seconds
+ Solving temperature system... 14 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00005
+ RMS, max velocity: 1.38e-09 m/s, 3.94e-09 m/s
+ Temperature min/avg/max: -202.9 K, 106.4 K, 1613 K
+ Heat fluxes through boundary parts: -0.116 W, -0.2413 W, 1.648e-06 W, -1.207e-05 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 9.45e-10 m/s, 1e-09 m/s, 1.38e-09 m/s
+
+*** Timestep 6: t=5.27266e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00006
+ RMS, max velocity: 1.41e-09 m/s, 3.98e-09 m/s
+ Temperature min/avg/max: -202.6 K, 107 K, 1613 K
+ Heat fluxes through boundary parts: -0.1147 W, -0.2403 W, 1.338e-06 W, -1.428e-05 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 9.69e-10 m/s, 1.02e-09 m/s, 1.41e-09 m/s
+
+Number of active cells: 1,209 (on 6 levels)
+Number of degrees of freedom: 17,345 (10,640+1,385+5,320)
+
+*** Timestep 7: t=5.9889e+13 seconds
+ Solving temperature system... 14 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 19 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00007
+ RMS, max velocity: 1.44e-09 m/s, 4.01e-09 m/s
+ Temperature min/avg/max: -202.4 K, 107.6 K, 1613 K
+ Heat fluxes through boundary parts: -0.1138 W, -0.2394 W, 6.598e-07 W, -4.841e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 9.92e-10 m/s, 1.04e-09 m/s, 1.44e-09 m/s
+
+*** Timestep 8: t=6.6987e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 18 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00008
+ RMS, max velocity: 1.47e-09 m/s, 4.05e-09 m/s
+ Temperature min/avg/max: -202.1 K, 108.1 K, 1613 K
+ Heat fluxes through boundary parts: -0.1136 W, -0.2386 W, 5.164e-07 W, 5.954e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 1.02e-09 m/s, 1.06e-09 m/s, 1.47e-09 m/s
+
+Number of active cells: 1,356 (on 7 levels)
+Number of degrees of freedom: 19,477 (11,950+1,552+5,975)
+
+*** Timestep 9: t=7.40219e+13 seconds
+ Solving temperature system... 15 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 17 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00009
+ RMS, max velocity: 1.5e-09 m/s, 4.08e-09 m/s
+ Temperature min/avg/max: -201.9 K, 108.5 K, 1613 K
+ Heat fluxes through boundary parts: -0.1138 W, -0.2378 W, 5.627e-07 W, 7.65e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 1.04e-09 m/s, 1.08e-09 m/s, 1.5e-09 m/s
+
+*** Timestep 10: t=8.10056e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 19 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-spherical_velocity_statistics/solution-00010
+ RMS, max velocity: 1.53e-09 m/s, 4.11e-09 m/s
+ Temperature min/avg/max: -201.7 K, 108.9 K, 1613 K
+ Heat fluxes through boundary parts: -0.1144 W, -0.2371 W, 1.18e-06 W, 6.828e-06 W
+ Writing depth average output-spherical_velocity_statistics/depth_average.gnuplot
+ Radial RMS, tangential RMS, total RMS velocity: 1.06e-09 m/s, 1.1e-09 m/s, 1.53e-09 m/s
+
+Number of active cells: 1,512 (on 7 levels)
+Number of degrees of freedom: 21,843 (13,406+1,734+6,703)
+
+Termination requested by criterion: end step
+*** Snapshot created!
+
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 8.91s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 11 | 0.66s | 7.4% |
+| Assemble temperature system | 11 | 1.14s | 13% |
+| Build Stokes preconditioner | 6 | 0.691s | 7.8% |
+| Build temperature preconditioner| 11 | 0.214s | 2.4% |
+| Solve Stokes system | 11 | 1.06s | 12% |
+| Solve temperature system | 11 | 0.0984s | 1.1% |
+| Create snapshot | 1 | 0.126s | 1.4% |
+| Initialization | 2 | 0.0215s | 0.24% |
+| Postprocessing | 11 | 3.9s | 44% |
+| Refine mesh structure, part 1 | 7 | 0.273s | 3.1% |
+| Refine mesh structure, part 2 | 7 | 0.0291s | 0.33% |
+| Setup dof systems | 8 | 0.639s | 7.2% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/spherical_velocity_statistics/statistics b/tests/spherical_velocity_statistics/statistics
new file mode 100644
index 0000000..43b14ff
--- /dev/null
+++ b/tests/spherical_velocity_statistics/statistics
@@ -0,0 +1,32 @@
+# 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: Iterations for temperature solver
+# 7: Iterations for Stokes solver
+# 8: Time step size (seconds)
+# 9: Visualization file name
+# 10: RMS velocity (m/s)
+# 11: Max. velocity (m/s)
+# 12: Minimal temperature (K)
+# 13: Average temperature (K)
+# 14: Maximal temperature (K)
+# 15: Outward heat flux through boundary with indicator 0 (W)
+# 16: Outward heat flux through boundary with indicator 1 (W)
+# 17: Outward heat flux through boundary with indicator 2 (W)
+# 18: Outward heat flux through boundary with indicator 3 (W)
+# 19: Radial RMS velocity (m/s)
+# 20: Tangential RMS velocity (m/s)
+# 21: Total RMS velocity (m/s)
+ 0 0.0000e+00 768 7235 3201 0 33 1.5684e+13 output-spherical_velocity_statistics/solution-00000 1.20522676e-09 3.64484888e-09 0.00000000e+00 1.01429624e+02 1.61300000e+03 -1.28972690e-01 -2.46880441e-01 -4.66932798e-06 -4.66932798e-06 8.09105435e-10 8.93263640e-10 1.20522676e-09
+ 1 1.5684e+13 912 8946 3958 19 24 7.5893e+12 output-spherical_velocity_statistics/solution-00001 1.26276581e-09 3.75246539e-09 -2.02154639e+02 1.02821383e+02 1.61300000e+03 -1.25594741e-01 -2.45686247e-01 -2.20411120e-06 -6.23503138e-06 8.53719710e-10 9.30451578e-10 1.26276581e-09
+ 2 2.3273e+13 912 8946 3958 13 20 7.4936e+12 output-spherical_velocity_statistics/solution-00002 1.29174210e-09 3.80085181e-09 -2.02352134e+02 1.03788568e+02 1.61300000e+03 -1.22866276e-01 -2.44387644e-01 -1.31325397e-06 -8.75368419e-06 8.76307862e-10 9.49042772e-10 1.29174210e-09
+ 3 3.0767e+13 990 9786 4329 12 20 7.4042e+12 output-spherical_velocity_statistics/solution-00003 1.32091007e-09 3.84717699e-09 -2.02556042e+02 1.04688550e+02 1.61300000e+03 -1.20239797e-01 -2.43490988e-01 -3.30878971e-07 -1.01517415e-05 8.99124311e-10 9.67666727e-10 1.32091007e-09
+ 4 3.8171e+13 990 9786 4329 12 19 7.3174e+12 output-spherical_velocity_statistics/solution-00004 1.35017220e-09 3.89186997e-09 -2.02872140e+02 1.05589658e+02 1.61300000e+03 -1.17836300e-01 -2.42360720e-01 1.71333442e-07 -1.15090106e-05 9.22097820e-10 9.86255845e-10 1.35017220e-09
+ 5 4.5488e+13 1092 10814 4784 14 20 7.2385e+12 output-spherical_velocity_statistics/solution-00005 1.37961950e-09 3.93562579e-09 -2.02853244e+02 1.06353285e+02 1.61300000e+03 -1.15998972e-01 -2.41329223e-01 1.64805687e-06 -1.20710820e-05 9.45314170e-10 1.00485376e-09 1.37961950e-09
+ 6 5.2727e+13 1092 10814 4784 14 20 7.1624e+12 output-spherical_velocity_statistics/solution-00006 1.40908707e-09 3.97607791e-09 -2.02621148e+02 1.07035669e+02 1.61300000e+03 -1.14659992e-01 -2.40345976e-01 1.33815066e-06 -1.42771282e-05 9.68640081e-10 1.02335857e-09 1.40908707e-09
+ 7 5.9889e+13 1209 12025 5320 14 19 7.0980e+12 output-spherical_velocity_statistics/solution-00007 1.43861733e-09 4.01221740e-09 -2.02382834e+02 1.07604258e+02 1.61300000e+03 -1.13815411e-01 -2.39423208e-01 6.59831030e-07 -4.84055594e-06 9.92082746e-10 1.04182131e-09 1.43861733e-09
+ 8 6.6987e+13 1209 12025 5320 14 18 7.0350e+12 output-spherical_velocity_statistics/solution-00008 1.46811557e-09 4.04813091e-09 -2.02143532e+02 1.08090112e+02 1.61300000e+03 -1.13580761e-01 -2.38552219e-01 5.16365927e-07 5.95415418e-06 1.01557478e-09 1.06017508e-09 1.46811557e-09
+ 9 7.4022e+13 1356 13502 5975 15 17 6.9837e+12 output-spherical_velocity_statistics/solution-00009 1.49748550e-09 4.07892319e-09 -2.01905106e+02 1.08518814e+02 1.61300000e+03 -1.13779733e-01 -2.37778656e-01 5.62717161e-07 7.64971464e-06 1.03904047e-09 1.07835880e-09 1.49748550e-09
+10 8.1006e+13 1356 13502 5975 14 19 6.9340e+12 output-spherical_velocity_statistics/solution-00010 1.52666686e-09 4.10696634e-09 -2.01667829e+02 1.08911515e+02 1.61300000e+03 -1.14396054e-01 -2.37078396e-01 1.18045064e-06 6.82840632e-06 1.06241895e-09 1.09634743e-09 1.52666686e-09
More information about the CIG-COMMITS
mailing list