-
Notifications
You must be signed in to change notification settings - Fork 240
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
4 changed files
with
2,141 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,210 @@ | ||
# This checks whether the prescribed dilation option | ||
# for plastic deformation in the visco plastic material | ||
# model successfully runs for 2 nonlinear iterations. | ||
# The test was derived from the kaus_2010 "brick" | ||
# extension benchmark, which was modified to remove | ||
# elasticity, have a lower resolution, and various | ||
# other adapations. | ||
# | ||
# Global parameters | ||
set Dimension = 2 | ||
set Start time = 0 | ||
set End time = 0e3 | ||
set Use years in output instead of seconds = true | ||
set Nonlinear solver scheme = single Advection, iterated Stokes | ||
set Nonlinear solver tolerance = 1e-3 | ||
set Max nonlinear iterations = 2 | ||
set CFL number = 0.5 | ||
set Maximum time step = 1000 | ||
set Output directory = visco_plastic_prescribed_plastic_dilation | ||
set Timing output frequency = 1 | ||
set Pressure normalization = no | ||
|
||
# Solver settings | ||
subsection Solver parameters | ||
subsection Stokes solver parameters | ||
set Stokes solver type = block AMG | ||
set Number of cheap Stokes solver steps = 2000 | ||
set Linear solver tolerance = 1e-7 | ||
set GMRES solver restart length = 50 | ||
set Use full A block as preconditioner = true | ||
end | ||
end | ||
|
||
# Model geometry (40x10 km, 10 km spacing) | ||
subsection Geometry model | ||
set Model name = box | ||
|
||
subsection Box | ||
set X repetitions = 4 | ||
set Y repetitions = 1 | ||
set X extent = 40e3 | ||
set Y extent = 10e3 | ||
end | ||
end | ||
|
||
# Mesh refinement specifications. | ||
# Globally refine to 1.25 km resolution | ||
subsection Mesh refinement | ||
set Initial adaptive refinement = 0 | ||
set Initial global refinement = 3 | ||
set Time steps between mesh refinement = 0 | ||
end | ||
|
||
# Element types | ||
subsection Discretization | ||
set Composition polynomial degree = 2 | ||
set Stokes velocity polynomial degree = 2 | ||
set Temperature polynomial degree = 1 | ||
end | ||
|
||
# Formulation classification | ||
subsection Formulation | ||
set Enable prescribed dilation = true | ||
end | ||
|
||
subsection Mesh deformation | ||
set Mesh deformation boundary indicators = top: free surface | ||
subsection Free surface | ||
set Surface velocity projection = normal | ||
end | ||
set Additional tangential mesh velocity boundary indicators = left, right | ||
end | ||
|
||
# Velocity boundary conditions | ||
# The imposed velocity produces a background strain-rate of 2e-15. | ||
subsection Boundary velocity model | ||
set Tangential velocity boundary indicators = bottom | ||
set Prescribed velocity boundary indicators = left x: function, right x:function | ||
|
||
subsection Function | ||
set Variable names = x,y | ||
set Function constants = cm=0.01, year=1, vel=0.12616 | ||
set Function expression = if (x<20e3 , -vel*cm/year, vel*cm/year); 0; | ||
end | ||
end | ||
|
||
# Number and name of compositional fields | ||
subsection Compositional fields | ||
set Number of fields = 2 | ||
set Names of fields = plastic_strain, inclusion | ||
set Types of fields = strain, chemical composition | ||
set Compositional field methods = particles | ||
set Mapped particle properties = plastic_strain: plastic_strain, \ | ||
inclusion: initial inclusion | ||
end | ||
|
||
# Spatial domain of different compositional fields | ||
subsection Initial composition model | ||
set Model name = function | ||
|
||
subsection Function | ||
set Variable names = x,y | ||
set Function constants = | ||
set Function expression = 0; if (y<=0.625e3 && x>=19.375e3 && x<=20.375e3, 1, 0); | ||
end | ||
end | ||
|
||
# Composition boundary conditions | ||
subsection Boundary composition model | ||
set Fixed composition boundary indicators = | ||
set List of model names = initial composition | ||
end | ||
|
||
# Temperature boundary conditions | ||
subsection Boundary temperature model | ||
set Fixed temperature boundary indicators = bottom, top, left, right | ||
set List of model names = initial temperature | ||
end | ||
|
||
# Temperature initial conditions (isothermal) | ||
subsection Initial temperature model | ||
set Model name = function | ||
|
||
subsection Function | ||
set Function expression = 293 | ||
end | ||
end | ||
|
||
# Gravity model | ||
subsection Gravity model | ||
set Model name = vertical | ||
|
||
subsection Vertical | ||
set Magnitude = 10. | ||
end | ||
end | ||
|
||
# Material model | ||
subsection Material model | ||
|
||
set Model name = visco plastic | ||
|
||
set Material averaging = harmonic average only viscosity | ||
|
||
subsection Visco Plastic | ||
|
||
set Densities = 2700 | ||
set Reference strain rate = 1.e-15 | ||
set Maximum viscosity = 1.e25 | ||
set Minimum viscosity = 1.e20 | ||
|
||
set Prefactors for dislocation creep = 5e-26, 5e-26, 5e-21 | ||
set Stress exponents for dislocation creep = 1.0 | ||
set Activation energies for dislocation creep = 0. | ||
set Activation volumes for dislocation creep = 0. | ||
|
||
set Viscosity averaging scheme = harmonic | ||
|
||
set Angles of internal friction = 30. | ||
set Cohesions = 40.e6, 40.e6, 1.e20 | ||
|
||
set Strain weakening mechanism = plastic weakening with plastic strain only | ||
set Start plasticity strain weakening intervals = 0.0 | ||
set End plasticity strain weakening intervals = 0.1 | ||
set Cohesion strain weakening factors = 0.1 | ||
set Friction strain weakening factors = 1.0 | ||
|
||
set Use plastic damper = false | ||
set Plastic damper viscosity = 1e20 | ||
|
||
end | ||
end | ||
|
||
# Post processing | ||
subsection Postprocess | ||
set List of postprocessors = basic statistics, composition statistics, particles, velocity statistics, visualization | ||
|
||
subsection Visualization | ||
set List of output variables = material properties, strain rate, named additional outputs | ||
set Time between graphical output = 0e3 | ||
set Interpolate output = false | ||
set Output format = gnuplot | ||
subsection Material properties | ||
set List of material properties = density, viscosity | ||
end | ||
end | ||
|
||
subsection Particles | ||
set Minimum particles per cell = 40 | ||
set Maximum particles per cell = 60 | ||
set Load balancing strategy = remove and add particles | ||
set List of particle properties = initial composition, viscoplastic strain invariants | ||
set Interpolation scheme = cell average | ||
set Update ghost particles = true | ||
set Particle generator name = reference cell | ||
subsection Generator | ||
subsection Reference cell | ||
set Number of particles per cell per direction = 7 | ||
end | ||
end | ||
set Time between data output = 0.e3 | ||
set Data output format = none | ||
end | ||
|
||
end | ||
|
||
# Termination criteria | ||
subsection Termination criteria | ||
set Termination criteria = end time | ||
end |
109 changes: 109 additions & 0 deletions
109
tests/visco_plastic_prescribed_plastic_dilation/screen-output
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,109 @@ | ||
----------------------------------------------------------------------------- | ||
-- This is ASPECT -- | ||
-- The Advanced Solver for Planetary Evolution, Convection, and Tectonics. -- | ||
----------------------------------------------------------------------------- | ||
-- . version 2.6.0-pre (add_prescribed_dilation_to_visco_plastic, 8bb195ccb) | ||
-- . using deal.II 9.5.1 | ||
-- . with 32 bit indices and vectorization level 2 (256 bits) | ||
-- . using Trilinos 13.2.0 | ||
-- . using p4est 2.3.2 | ||
-- . using Geodynamic World Builder 0.6.0 | ||
-- . 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.6.0-pre&particles=1&sha=8bb195ccb&src=code | ||
----------------------------------------------------------------------------- | ||
Number of active cells: 256 (on 4 levels) | ||
Number of degrees of freedom: 5,014 (2,210+297+297+1,105+1,105) | ||
|
||
Number of mesh deformation degrees of freedom: 594 | ||
Solving mesh displacement system... 0 iterations. | ||
*** Timestep 0: t=0 years, dt=0 years | ||
Solving mesh displacement system... 0 iterations. | ||
Solving temperature system... 0 iterations. | ||
Advecting particles... done. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 35+0 iterations. | ||
Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1 | ||
|
||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 30+0 iterations. | ||
Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.0197761 | ||
|
||
|
||
Postprocessing: | ||
Compositions min/max/mass: 0/0/0 // 0/0.3265/6.378e+05 | ||
Number of advected particles: 12544 | ||
RMS, max velocity: 0.000769 m/year, 0.00126 m/year | ||
Writing graphical output: output_visco_plastic_prescribed_plastic_dilation/solution/solution-00000 | ||
|
||
|
||
|
||
+----------------------------------------------+------------+------------+ | ||
| Total wallclock time elapsed since start | 3.8s | | | ||
| | | | | ||
| Section | no. calls | wall time | % of total | | ||
+----------------------------------+-----------+------------+------------+ | ||
| Assemble Stokes system | 2 | 0.669s | 18% | | ||
| Assemble temperature system | 1 | 0.258s | 6.8% | | ||
| Build Stokes preconditioner | 2 | 0.256s | 6.7% | | ||
| Build temperature preconditioner | 1 | 0.000402s | 0% | | ||
| Initialization | 1 | 0.631s | 17% | | ||
| Mesh deformation | 1 | 0.0939s | 2.5% | | ||
| Mesh deformation initialize | 1 | 0.0784s | 2.1% | | ||
| Particles: Advect | 2 | 0.0818s | 2.2% | | ||
| Particles: Copy | 1 | 0.00106s | 0% | | ||
| Particles: Generate | 1 | 0.149s | 3.9% | | ||
| Particles: Initialization | 1 | 0.000195s | 0% | | ||
| Particles: Initialize properties | 1 | 0.0535s | 1.4% | | ||
| Particles: Interpolate | 1 | 0.0202s | 0.53% | | ||
| Particles: Sort | 2 | 0.0311s | 0.82% | | ||
| Particles: Update properties | 1 | 0.17s | 4.5% | | ||
| Postprocessing | 1 | 0.194s | 5.1% | | ||
| Setup dof systems | 1 | 0.129s | 3.4% | | ||
| Setup initial conditions | 1 | 0.642s | 17% | | ||
| Setup matrices | 1 | 0.134s | 3.5% | | ||
| Solve Stokes system | 2 | 0.119s | 3.1% | | ||
| Solve temperature system | 1 | 0.000714s | 0% | | ||
+----------------------------------+-----------+------------+------------+ | ||
|
||
-- Total wallclock time elapsed including restarts: 4s | ||
Termination requested by criterion: end time | ||
|
||
|
||
+----------------------------------------------+------------+------------+ | ||
| Total wallclock time elapsed since start | 3.8s | | | ||
| | | | | ||
| Section | no. calls | wall time | % of total | | ||
+----------------------------------+-----------+------------+------------+ | ||
| Assemble Stokes system | 2 | 0.669s | 18% | | ||
| Assemble temperature system | 1 | 0.258s | 6.8% | | ||
| Build Stokes preconditioner | 2 | 0.256s | 6.7% | | ||
| Build temperature preconditioner | 1 | 0.000402s | 0% | | ||
| Initialization | 1 | 0.631s | 17% | | ||
| Mesh deformation | 1 | 0.0939s | 2.5% | | ||
| Mesh deformation initialize | 1 | 0.0784s | 2.1% | | ||
| Particles: Advect | 2 | 0.0818s | 2.2% | | ||
| Particles: Copy | 1 | 0.00106s | 0% | | ||
| Particles: Generate | 1 | 0.149s | 3.9% | | ||
| Particles: Initialization | 1 | 0.000195s | 0% | | ||
| Particles: Initialize properties | 1 | 0.0535s | 1.4% | | ||
| Particles: Interpolate | 1 | 0.0202s | 0.53% | | ||
| Particles: Sort | 2 | 0.0311s | 0.82% | | ||
| Particles: Update properties | 1 | 0.17s | 4.5% | | ||
| Postprocessing | 1 | 0.194s | 5.1% | | ||
| Setup dof systems | 1 | 0.129s | 3.4% | | ||
| Setup initial conditions | 1 | 0.642s | 17% | | ||
| Setup matrices | 1 | 0.134s | 3.5% | | ||
| Solve Stokes system | 2 | 0.119s | 3.1% | | ||
| Solve temperature system | 1 | 0.000714s | 0% | | ||
+----------------------------------+-----------+------------+------------+ | ||
|
||
-- Total wallclock time elapsed including restarts: 4s | ||
----------------------------------------------------------------------------- | ||
-- For information on how to cite ASPECT, see: | ||
-- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&particles=1&sha=8bb195ccb&src=code | ||
----------------------------------------------------------------------------- |
Oops, something went wrong.