[cig-commits] commit 2492 by heister to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Sat Apr 12 06:40:36 PDT 2014
Revision 2492
add test
A trunk/aspect/tests/refine_vel/
A trunk/aspect/tests/refine_vel/screen-output
A trunk/aspect/tests/refine_vel/statistics
A trunk/aspect/tests/refine_vel.prm
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2492&peg=2492
Diff:
Added: trunk/aspect/tests/refine_vel/screen-output
===================================================================
--- trunk/aspect/tests/refine_vel/screen-output (rev 0)
+++ trunk/aspect/tests/refine_vel/screen-output 2014-04-12 13:40:34 UTC (rev 2492)
@@ -0,0 +1,202 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Simulator for Problems in Earth's ConvecTion.
+-- . running in DEBUG 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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+Number of active cells: 837 (on 6 levels)
+Number of degrees of freedom: 11,660 (7,154+929+3,577)
+
+AAA 3
+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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+*** Timestep 2: t=2.32729e+13 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+*** Timestep 4: t=3.81707e+13 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 19 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+*** Timestep 6: t=5.27266e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+*** Timestep 8: t=6.6987e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 18 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+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-refine_vel/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-refine_vel/depth_average.gnuplot
+
+*** Timestep 10: t=8.10056e+13 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 19 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-refine_vel/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-refine_vel/depth_average.gnuplot
+
+AAA 3
+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 | 19s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 11 | 1.16s | 6.1% |
+| Assemble temperature system | 11 | 2.38s | 12% |
+| Build Stokes preconditioner | 6 | 1.27s | 6.7% |
+| Build temperature preconditioner| 11 | 0.0921s | 0.48% |
+| Solve Stokes system | 11 | 0.693s | 3.6% |
+| Solve temperature system | 11 | 0.0691s | 0.36% |
+| Create snapshot | 1 | 0.0481s | 0.25% |
+| Initialization | 2 | 0.0224s | 0.12% |
+| Postprocessing | 11 | 11.4s | 60% |
+| Refine mesh structure, part 1 | 7 | 0.865s | 4.5% |
+| Refine mesh structure, part 2 | 7 | 0.0268s | 0.14% |
+| Setup dof systems | 8 | 0.941s | 4.9% |
++---------------------------------+-----------+------------+------------+
+
Added: trunk/aspect/tests/refine_vel/statistics
===================================================================
--- trunk/aspect/tests/refine_vel/statistics (rev 0)
+++ trunk/aspect/tests/refine_vel/statistics 2014-04-12 13:40:34 UTC (rev 2492)
@@ -0,0 +1,29 @@
+# 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)
+ 0 0.0000e+00 768 7235 3201 0 33 1.5684e+13 output-refine_vel/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
+ 1 1.5684e+13 912 8946 3958 19 24 7.5893e+12 output-refine_vel/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
+ 2 2.3273e+13 912 8946 3958 13 20 7.4936e+12 output-refine_vel/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
+ 3 3.0767e+13 990 9786 4329 12 20 7.4042e+12 output-refine_vel/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
+ 4 3.8171e+13 990 9786 4329 12 19 7.3174e+12 output-refine_vel/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
+ 5 4.5488e+13 1092 10814 4784 14 20 7.2385e+12 output-refine_vel/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
+ 6 5.2727e+13 1092 10814 4784 14 20 7.1624e+12 output-refine_vel/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
+ 7 5.9889e+13 1209 12025 5320 14 19 7.0980e+12 output-refine_vel/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
+ 8 6.6987e+13 1209 12025 5320 14 18 7.0350e+12 output-refine_vel/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
+ 9 7.4022e+13 1356 13502 5975 15 17 6.9837e+12 output-refine_vel/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
+10 8.1006e+13 1356 13502 5975 14 19 6.9340e+12 output-refine_vel/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
Added: trunk/aspect/tests/refine_vel.prm
===================================================================
--- trunk/aspect/tests/refine_vel.prm (rev 0)
+++ trunk/aspect/tests/refine_vel.prm 2014-04-12 13:40:34 UTC (rev 2492)
@@ -0,0 +1,109 @@
+set Dimension = 2
+set CFL number = 1.0
+set End time = 1e300
+set Start time = 0
+set Adiabatic surface temperature = 1613.0
+set Surface pressure = 0
+set Use years in output instead of seconds = false # default: true
+set Nonlinear solver scheme = IMPES
+
+
+subsection Boundary temperature model
+ set Model name = spherical constant
+
+ subsection Spherical constant
+ set Inner temperature = 1613.0
+ set Outer temperature = 1613.0
+ end
+end
+
+subsection Termination criteria
+ set Checkpoint on termination = true
+ set Termination criteria = end step
+ set End step = 10
+
+end
+
+subsection Gravity model
+ set Model name = radial constant
+
+ subsection Radial constant
+ set Magnitude = 9.81
+ end
+end
+
+
+subsection Geometry model
+ set Model name = spherical shell
+
+ subsection Spherical shell
+ set Inner radius = 3481000
+ set Outer radius = 6371000
+ set Opening angle = 90
+ end
+end
+
+
+subsection Initial conditions
+ set Model name = function
+
+ subsection Function
+ set Function expression = 1613.0
+ end
+end
+
+
+subsection Material model
+ set Model name = simple
+
+ subsection Simple model
+ set Reference density = 3300
+ set Reference specific heat = 1250
+ set Reference temperature = 1613 # default: 293
+ set Thermal conductivity = 1e-6 # default: 4.7
+ set Thermal expansion coefficient = 2e-5
+ set Viscosity = 1e22 # default: 5e24
+ end
+end
+
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 0
+ set Initial global refinement = 4
+ set Additional refinement times = 0,0
+
+ set Strategy = velocity
+ set Time steps between mesh refinement = 2
+
+ set Refinement fraction = 0.1
+ set Coarsening fraction = 0.0
+
+
+
+end
+
+subsection Initial conditions
+ set Model name = function
+ subsection Function
+ set Variable names = x,y
+ set Function expression = if((sqrt((x-3e6)^2+(y-2e6)^2)<1e6) , 800.0, 0)
+ end
+end
+
+subsection Model settings
+ set Fixed temperature boundary indicators = 0, 1
+ set Prescribed velocity boundary indicators =
+ set Tangential velocity boundary indicators = 0, 2, 3
+ set Zero velocity boundary indicators = 1
+end
+
+subsection Postprocess
+ set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average
+
+ subsection Visualization
+ set Time between graphical output = 0
+
+ end
+
+end
+
More information about the CIG-COMMITS
mailing list