[cig-commits] [commit] aspect-1.2: remove broken periodic_box for now (a8eccba)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Jan 23 09:52:51 PST 2015
Repository : https://github.com/geodynamics/aspect
On branch : aspect-1.2
Link : https://github.com/geodynamics/aspect/compare/0000000000000000000000000000000000000000...8c1e1e978a72b986edd943b49cf75067db3bca79
>---------------------------------------------------------------
commit a8eccbaba04a5a4bfb076ea3ea5ef640e0cb12e3
Author: Timo Heister <timo.heister at gmail.com>
Date: Fri Jan 23 12:44:08 2015 -0500
remove broken periodic_box for now
>---------------------------------------------------------------
a8eccbaba04a5a4bfb076ea3ea5ef640e0cb12e3
cookbooks/future/periodic_box.prm | 105 -----------------------------
tests/periodic_box.prm | 104 -----------------------------
tests/periodic_box/screen-output | 134 --------------------------------------
tests/periodic_box/statistics | 24 -------
4 files changed, 367 deletions(-)
diff --git a/cookbooks/future/periodic_box.prm b/cookbooks/future/periodic_box.prm
deleted file mode 100644
index 7c97384..0000000
--- a/cookbooks/future/periodic_box.prm
+++ /dev/null
@@ -1,105 +0,0 @@
-set Dimension = 2
-set CFL number = 1.0
-set End time = 1e9
-set Output directory = output
-set Resume computation = false
-set Start time = 0
-set Adiabatic surface temperature = 0
-set Surface pressure = 0
-set Pressure normalization = no
-set Linear solver tolerance = 1.e-7
-set Use years in output instead of seconds = true
-set Nonlinear solver scheme = IMPES
-
-subsection Boundary temperature model
- set Model name = box
- subsection Box
- set Top temperature = 0.0
- set Bottom temperature = 1000.0
- end
-end
-
-
-subsection Discretization
- set Stokes velocity polynomial degree = 2
- set Temperature polynomial degree = 2
- set Use locally conservative discretization = false
- subsection Stabilization parameters
- set alpha = 2
- set beta = 0.078
- set cR = 0.5 # default: 0.11
- end
-end
-
-
-subsection Geometry model
- set Model name = box
- subsection Box
- set X periodic = true
- set X extent = 1.e6
- set Y extent = 5.e5
- set Z extent = 5.e5
- end
-end
-
-
-subsection Gravity model
- set Model name = vertical
- subsection Vertical
- set Magnitude = 10.0
- end
-end
-
-
-subsection Initial conditions
- set Model name = function
- subsection Function
- set Variable names = x,y
- set Function expression = if((sqrt((x-1.e5)^2+(y-4.0e5)^2)<5.0e4) | (sqrt((x-3.e5)^2+(y-2.e5)^2)<1.0e5) , 800.0, 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 = 0.0
- set Thermal conductivity = 4.7
- set Thermal expansion coefficient = 4e-5
- set Viscosity = 1.e20
- end
-end
-
-
-subsection Mesh refinement
- set Additional refinement times =
- set Initial adaptive refinement = 2 # default: 2
- set Initial global refinement = 6 # default: 2
- set Refinement fraction = 0.3
- set Coarsening fraction = 0.03
- set Strategy = thermal energy density
- set Time steps between mesh refinement = 10 # default: 10
-end
-
-
-subsection Model settings
- set Include adiabatic heating = false
- set Include shear heating = false
- set Fixed temperature boundary indicators = bottom, top
- set Prescribed velocity boundary indicators =
- set Tangential velocity boundary indicators = bottom, top
- set Zero velocity boundary indicators =
- set Remove nullspace = net x translation
-end
-
-subsection Postprocess
- set List of postprocessors = visualization
- subsection Visualization
- set List of output variables =
- set Number of grouped files = 1
- set Output format = vtu
- set Time between graphical output = 1.e3
- end
-end
diff --git a/tests/periodic_box.prm b/tests/periodic_box.prm
deleted file mode 100644
index 498f3af..0000000
--- a/tests/periodic_box.prm
+++ /dev/null
@@ -1,104 +0,0 @@
-#Test for running a box in 2D with periodic boundary conditions
-#in the x condition. Also tests the removal of net translations
-#in the periodic direction
-
-set Dimension = 2
-set CFL number = 1.0
-set End time = 1e6
-set Output directory = output
-set Adiabatic surface temperature = 0
-set Surface pressure = 0
-set Pressure normalization = no
-set Linear solver tolerance = 1.e-7
-set Use years in output instead of seconds = true
-set Nonlinear solver scheme = IMPES
-
-subsection Boundary temperature model
- set Model name = constant
- subsection Constant
- set Boundary indicator to temperature mappings = 0:0,1:0,2:0,3:0
- end
-end
-
-
-subsection Discretization
- set Stokes velocity polynomial degree = 2
- set Temperature polynomial degree = 2
- set Use locally conservative discretization = false
- subsection Stabilization parameters
- set alpha = 2
- set beta = 0.078
- set cR = 0.5 # default: 0.11
- end
-end
-
-
-subsection Geometry model
- set Model name = box
- subsection Box
- set X periodic = true
- set X extent = 1.e6
- set Y extent = 5.e5
- set Z extent = 5.e5
- end
-end
-
-
-subsection Gravity model
- set Model name = vertical
- subsection Vertical
- set Magnitude = 10.0
- end
-end
-
-
-subsection Initial conditions
- set Model name = function
- subsection Function
- set Variable names = x,y
- set Function expression = if((sqrt((x-1.e5)^2+(y-4.0e5)^2)<5.0e4) | (sqrt((x-3.e5)^2+(y-2.e5)^2)<1.0e5) , 800.0, 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 = 0.0
- set Thermal conductivity = 4.7
- set Thermal expansion coefficient = 4e-5
- set Viscosity = 1.e20
- end
-end
-
-subsection Mesh refinement
- set Additional refinement times =
- set Initial adaptive refinement = 1 # default: 2
- set Initial global refinement = 5 # default: 2
- set Refinement fraction = 0.3
- set Coarsening fraction = 0.05
- set Strategy = thermal energy density
- set Time steps between mesh refinement = 10 # default: 10
-end
-
-
-subsection Model settings
- set Include adiabatic heating = false
- set Include shear heating = false
- set Fixed temperature boundary indicators = 2,3
- set Prescribed velocity boundary indicators =
- set Tangential velocity boundary indicators = 2,3
- set Zero velocity boundary indicators =
- set Remove nullspace = net x translation
-end
-
-subsection Termination criteria
- set Termination criteria = end step
- set End step = 10
-end
-
-subsection Postprocess
- set List of postprocessors = velocity statistics, basic statistics
-end
diff --git a/tests/periodic_box/screen-output b/tests/periodic_box/screen-output
deleted file mode 100644
index 104c8de..0000000
--- a/tests/periodic_box/screen-output
+++ /dev/null
@@ -1,134 +0,0 @@
------------------------------------------------------------------------------
--- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
--- . version 1.2.pre
--- . running in DEBUG mode
--- . running with 1 MPI process
--- . using Trilinos
------------------------------------------------------------------------------
-
-Number of active cells: 1,024 (on 6 levels)
-Number of degrees of freedom: 13,764 (8,450+1,089+4,225)
-
-*** Timestep 0: t=0 years
- Solving temperature system... 0 iterations.
- Rebuilding Stokes preconditioner...
- Solving Stokes system... 30+4 iterations.
-
-Number of active cells: 493 (on 7 levels)
-Number of degrees of freedom: 7,441 (4,568+589+2,284)
-
-*** Timestep 0: t=0 years
- Solving temperature system... 0 iterations.
- Rebuilding Stokes preconditioner...
- Solving Stokes system... 30+5 iterations.
-
- Postprocessing:
-
- Reference density (kg/m^3): 3300
- Reference gravity (m/s^2): 10
- Reference thermal expansion (1/K): 4e-05
- Temperature contrast across model domain (K): 0
- Model domain depth (m): 500000
- Reference thermal diffusivity (m^2/s): 1.13939e-06
- Reference viscosity (Pas): 1e+20
- Ra number: 0
- k_value: 4.7
- reference_cp: 1250
- reference_thermal_diffusivity: 1.13939e-06
-
- RMS, max velocity: 0.385 m/year, 0.913 m/year
-
-*** Timestep 1: t=4945 years
- Solving temperature system... 13 iterations.
- Solving Stokes system... 30+3 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.386 m/year, 0.916 m/year
-
-*** Timestep 2: t=9825.96 years
- Solving temperature system... 11 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.387 m/year, 0.919 m/year
-
-*** Timestep 3: t=14648.8 years
- Solving temperature system... 12 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.388 m/year, 0.921 m/year
-
-*** Timestep 4: t=19419.1 years
- Solving temperature system... 11 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.388 m/year, 0.923 m/year
-
-*** Timestep 5: t=24142.3 years
- Solving temperature system... 11 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.388 m/year, 0.924 m/year
-
-*** Timestep 6: t=28823.5 years
- Solving temperature system... 11 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.389 m/year, 0.925 m/year
-
-*** Timestep 7: t=33467.4 years
- Solving temperature system... 11 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.389 m/year, 0.925 m/year
-
-*** Timestep 8: t=38078.4 years
- Solving temperature system... 10 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.389 m/year, 0.925 m/year
-
-*** Timestep 9: t=42660.6 years
- Solving temperature system... 10 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.388 m/year, 0.924 m/year
-
-*** Timestep 10: t=47217.6 years
- Solving temperature system... 10 iterations.
- Solving Stokes system... 30+2 iterations.
-
- Postprocessing:
- RMS, max velocity: 0.388 m/year, 0.923 m/year
-
-Number of active cells: 427 (on 7 levels)
-Number of degrees of freedom: 6,638 (4,076+524+2,038)
-
-Termination requested by criterion: end step
-
-
-+---------------------------------------------+------------+------------+
-| Total wallclock time elapsed since start | 25.8s | |
-| | | |
-| Section | no. calls | wall time | % of total |
-+---------------------------------+-----------+------------+------------+
-| Assemble Stokes system | 12 | 4.84s | 19% |
-| Assemble temperature system | 12 | 5.77s | 22% |
-| Build Stokes preconditioner | 2 | 2.53s | 9.8% |
-| Build temperature preconditioner| 12 | 0.125s | 0.49% |
-| Solve Stokes system | 12 | 2.78s | 11% |
-| Solve temperature system | 12 | 0.125s | 0.49% |
-| Initialization | 3 | 0.745s | 2.9% |
-| Postprocessing | 11 | 0.241s | 0.93% |
-| Refine mesh structure, part 1 | 2 | 3.64s | 14% |
-| Refine mesh structure, part 2 | 2 | 0.0247s | 0.096% |
-| Setup dof systems | 3 | 2.07s | 8% |
-+---------------------------------+-----------+------------+------------+
-
diff --git a/tests/periodic_box/statistics b/tests/periodic_box/statistics
deleted file mode 100644
index 12bc39d..0000000
--- a/tests/periodic_box/statistics
+++ /dev/null
@@ -1,24 +0,0 @@
-# 1: Time step number
-# 2: Time (years)
-# 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: Velocity iterations in Stokes preconditioner
-# 9: Schur complement iterations in Stokes preconditioner
-# 10: Time step size (years)
-# 11: RMS velocity (m/year)
-# 12: Max. velocity (m/year)
- 0 0.0000e+00 1024 9539 4225 0 34 54 15 8.5731e+03 0.00000000e+00 0.00000000e+00
- 0 0.0000e+00 493 5157 2284 0 35 64 29 4.9450e+03 3.85023807e-01 9.13275756e-01
- 1 4.9450e+03 493 5157 2284 13 33 43 20 4.8810e+03 3.86043540e-01 9.16294838e-01
- 2 9.8260e+03 493 5157 2284 11 32 33 15 4.8228e+03 3.86882724e-01 9.19028189e-01
- 3 1.4649e+04 493 5157 2284 12 32 35 15 4.7703e+03 3.87564470e-01 9.21097596e-01
- 4 1.9419e+04 493 5157 2284 11 32 35 15 4.7232e+03 3.88090644e-01 9.22830107e-01
- 5 2.4142e+04 493 5157 2284 11 32 34 15 4.6812e+03 3.88460696e-01 9.23894172e-01
- 6 2.8824e+04 493 5157 2284 11 32 34 15 4.6439e+03 3.88676423e-01 9.24695613e-01
- 7 3.3467e+04 493 5157 2284 11 32 36 15 4.6110e+03 3.88738347e-01 9.25031619e-01
- 8 3.8078e+04 493 5157 2284 10 32 39 15 4.5821e+03 3.88653384e-01 9.24749812e-01
- 9 4.2661e+04 493 5157 2284 10 32 36 15 4.5570e+03 3.88423413e-01 9.24251846e-01
-10 4.7218e+04 493 5157 2284 10 32 40 15 4.5354e+03 3.88053832e-01 9.23344025e-01
More information about the CIG-COMMITS
mailing list