Skip to content

Commit

Permalink
add test
Browse files Browse the repository at this point in the history
  • Loading branch information
@naliboff
naliboff committed Jul 12, 2023
1 parent 1fee6a5 commit abe8007
Show file tree
Hide file tree
Showing 3 changed files with 264 additions and 0 deletions.
124 changes: 124 additions & 0 deletions tests/iterative_dampening_viscosity_no_dampening.prm
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,124 @@
# Global parameters
set Dimension = 2
set Start time = 0
set End time = 0
set Use years in output instead of seconds = true
set Nonlinear solver scheme = iterated Advection and Stokes
set Nonlinear solver tolerance = 1e-5
set Max nonlinear iterations = 4
set Output directory = iterative_dampening_viscosity_no_dampening
set Timing output frequency = 1

# Model geometry (2x2 km, 1 km spacing)
subsection Geometry model
set Model name = box
subsection Box
set X repetitions = 2
set Y repetitions = 2
set X extent = 2e3
set Y extent = 2e3
end
end

# Mesh refinement specifications
subsection Mesh refinement
set Initial adaptive refinement = 0
set Initial global refinement = 0
set Time steps between mesh refinement = 0
end


# Boundary classifications (fixed T boundaries, prescribed velocity)
# The parameters below this comment were created by the update script
# as replacement for the old 'Model settings' subsection. They can be
# safely merged with any existing subsections with the same name.

subsection Boundary temperature model
set Fixed temperature boundary indicators = bottom, top, left, right
end

subsection Boundary velocity model
set Tangential velocity boundary indicators = bottom, top, left, right
end

# Temperature boundary and initial conditions
subsection Boundary temperature model
set List of model names = box
subsection Box
set Bottom temperature = 273
set Left temperature = 273
set Right temperature = 273
set Top temperature = 273
end
end
subsection Initial temperature model
set Model name = function
subsection Function
set Function expression = 273
end
end

# Compositional fields used to track finite strain invariant
subsection Compositional fields
set Number of fields = 1
set Names of fields = viscosity_field
end

# Set the initial composition for the viscosity field
# to the same value as the background viscosity.
# However, his value can be set to any arbitrary
# viscosity and still have the iterative dampening
# scheme work.
subsection Initial composition model
set Model name = function
subsection Function
set Variable names = x,y
set Function expression = 1e26;
end
end

# Material model (values for background material)
subsection Material model
set Model name = visco plastic
subsection Visco Plastic

# Turn on iterative viscosity dampening, but set the
# dampening value to 0 such that no dampening occurs.
set Use iterative viscosity dampening = true
set Iterative viscosity dampening factor = 0

# These strain rate given to the material will
# always be 1e-20, as all boundaries have
# free slip conditions.
set Reference strain rate = 1.e-20
set Minimum strain rate = 1.e-20

# Set the background viscosity to 1e26 P a s
set Viscous flow law = dislocation
set Prefactors for dislocation creep = 5.e-27
set Stress exponents for dislocation creep = 1.0
set Activation energies for dislocation creep = 0.
set Activation volumes for dislocation creep = 0.

# After yielding, the calculated viscosity should be
# 1.e6 [Pa] * ( 1. / (2. * 1.e-20 [1/s]) = 5e25 [Pa s]
set Yield mechanism = drucker
set Angles of internal friction = 0.
set Cohesions = 1.e6

end
end

# Gravity model
subsection Gravity model
set Model name = vertical
subsection Vertical
set Magnitude = 10.0
end
end

# Post processing
subsection Postprocess
set List of postprocessors = basic statistics, composition statistics, material statistics
set Run postprocessors on nonlinear iterations = true
end
117 changes: 117 additions & 0 deletions tests/iterative_dampening_viscosity_no_dampening/screen-output
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,117 @@
-----------------------------------------------------------------------------
-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
-- . version 2.5.0-pre (average_viscosity_nonlinear_itearations, b12e8343d)
-- . using deal.II 9.4.0
-- . with 32 bit indices and vectorization level 1 (128 bits)
-- . using Trilinos 13.2.0
-- . using p4est 2.3.2
-- . running in DEBUG mode
-- . running with 1 MPI process
-----------------------------------------------------------------------------

-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
-- https://aspect.geodynamics.org/citing.html?ver=2.5.0-pre&sha=b12e8343d&src=code
-----------------------------------------------------------------------------
Number of active cells: 4 (on 1 levels)
Number of degrees of freedom: 109 (50+9+25+25)

*** Timestep 0: t=0 years, dt=0 years
Solving temperature system... 0 iterations.
Solving viscosity_field system ... 6 iterations.
Rebuilding Stokes preconditioner...
Solving Stokes system... 2+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.47805e-16, 1, 1
Relative nonlinear residual (total system) after nonlinear iteration 1: 1

Postprocessing:
Compositions min/max/mass: 5e+25/5e+25/2e+32
Average density / Average viscosity / Total mass: 3302 kg/m^3, 5e+25 Pa s, 1.321e+10 kg

Solving temperature system... 0 iterations.
Solving viscosity_field system ... 6 iterations.
Rebuilding Stokes preconditioner...
Solving Stokes system... 12+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.47805e-16, 1, 5.05365e-14
Relative nonlinear residual (total system) after nonlinear iteration 2: 1

Postprocessing:
Compositions min/max/mass: 1e+26/1e+26/4e+32
Average density / Average viscosity / Total mass: 3302 kg/m^3, 5e+25 Pa s, 1.321e+10 kg

Solving temperature system... 0 iterations.
Solving viscosity_field system ... 6 iterations.
Rebuilding Stokes preconditioner...
Solving Stokes system... 13+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.47805e-16, 1, 7.98252e-15
Relative nonlinear residual (total system) after nonlinear iteration 3: 1

Postprocessing:
Compositions min/max/mass: 5e+25/5e+25/2e+32
Average density / Average viscosity / Total mass: 3302 kg/m^3, 5e+25 Pa s, 1.321e+10 kg

Solving temperature system... 0 iterations.
Solving viscosity_field system ... 6 iterations.
Rebuilding Stokes preconditioner...
Solving Stokes system... 12+0 iterations.
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.47805e-16, 1, 8.48347e-15
Relative nonlinear residual (total system) after nonlinear iteration 4: 1

Postprocessing:
Compositions min/max/mass: 1e+26/1e+26/4e+32
Average density / Average viscosity / Total mass: 3302 kg/m^3, 5e+25 Pa s, 1.321e+10 kg




+----------------------------------------------+------------+------------+
| Total wallclock time elapsed since start | 0.719s | |
| | | |
| Section | no. calls | wall time | % of total |
+----------------------------------+-----------+------------+------------+
| Assemble Stokes system | 4 | 0.0226s | 3.1% |
| Assemble composition system | 4 | 0.0333s | 4.6% |
| Assemble temperature system | 4 | 0.0294s | 4.1% |
| Build Stokes preconditioner | 4 | 0.025s | 3.5% |
| Build composition preconditioner | 4 | 0.0003s | 0% |
| Build temperature preconditioner | 4 | 0.000557s | 0% |
| Initialization | 1 | 0.553s | 77% |
| Postprocessing | 4 | 0.00985s | 1.4% |
| Setup dof systems | 1 | 0.0064s | 0.89% |
| Setup initial conditions | 1 | 0.0105s | 1.5% |
| Setup matrices | 1 | 0.00429s | 0.6% |
| Solve Stokes system | 4 | 0.0112s | 1.6% |
| Solve composition system | 4 | 0.00152s | 0.21% |
| Solve temperature system | 4 | 0.00136s | 0.19% |
+----------------------------------+-----------+------------+------------+

-- Total wallclock time elapsed including restarts: 1s
Termination requested by criterion: end time


+----------------------------------------------+------------+------------+
| Total wallclock time elapsed since start | 0.719s | |
| | | |
| Section | no. calls | wall time | % of total |
+----------------------------------+-----------+------------+------------+
| Assemble Stokes system | 4 | 0.0226s | 3.1% |
| Assemble composition system | 4 | 0.0333s | 4.6% |
| Assemble temperature system | 4 | 0.0294s | 4.1% |
| Build Stokes preconditioner | 4 | 0.025s | 3.5% |
| Build composition preconditioner | 4 | 0.0003s | 0% |
| Build temperature preconditioner | 4 | 0.000557s | 0% |
| Initialization | 1 | 0.553s | 77% |
| Postprocessing | 4 | 0.00985s | 1.4% |
| Setup dof systems | 1 | 0.0064s | 0.89% |
| Setup initial conditions | 1 | 0.0105s | 1.5% |
| Setup matrices | 1 | 0.00429s | 0.6% |
| Solve Stokes system | 4 | 0.0112s | 1.6% |
| Solve composition system | 4 | 0.00152s | 0.21% |
| Solve temperature system | 4 | 0.00136s | 0.19% |
+----------------------------------+-----------+------------+------------+

-- Total wallclock time elapsed including restarts: 1s
-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
-- https://aspect.geodynamics.org/citing.html?ver=2.5.0-pre&sha=b12e8343d&src=code
-----------------------------------------------------------------------------
23 changes: 23 additions & 0 deletions tests/iterative_dampening_viscosity_no_dampening/statistics
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,23 @@
# 1: Time step number
# 2: Time (years)
# 3: Time step size (years)
# 4: Number of mesh cells
# 5: Number of Stokes degrees of freedom
# 6: Number of temperature degrees of freedom
# 7: Number of degrees of freedom for all compositions
# 8: Number of nonlinear iterations
# 9: Iterations for temperature solver
# 10: Iterations for composition solver 1
# 11: Iterations for Stokes solver
# 12: Velocity iterations in Stokes preconditioner
# 13: Schur complement iterations in Stokes preconditioner
# 14: Minimal value for composition viscosity_field
# 15: Maximal value for composition viscosity_field
# 16: Global mass for composition viscosity_field
# 17: Average density (kg/m^3)
# 18: Average viscosity (Pa s)
# 19: Total mass (kg)
0 0.000000000000e+00 0.000000000000e+00 4 59 25 25 1 0 6 1 3 2 5.00000000e+25 5.00000000e+25 2.00000000e+32 3.30231000e+03 5.00000000e+25 1.32092400e+10
0 0.000000000000e+00 0.000000000000e+00 4 59 25 25 1 0 6 11 13 13 1.00000000e+26 1.00000000e+26 4.00000000e+32 3.30231000e+03 5.00000000e+25 1.32092400e+10
0 0.000000000000e+00 0.000000000000e+00 4 59 25 25 1 0 6 12 14 14 5.00000000e+25 5.00000000e+25 2.00000000e+32 3.30231000e+03 5.00000000e+25 1.32092400e+10
0 0.000000000000e+00 0.000000000000e+00 4 59 25 25 1 0 6 11 13 13 1.00000000e+26 1.00000000e+26 4.00000000e+32 3.30231000e+03 5.00000000e+25 1.32092400e+10

0 comments on commit abe8007

Please sign in to comment.