Skip to content

Commit b9e8686

Browse files
danieldouglas92danieldouglas92
committed
add test case
1 parent c8caa8a commit b9e8686

File tree

3 files changed

+226
-0
lines changed

3 files changed

+226
-0
lines changed

‎tests/tian_track_fluid_sources.prm‎

Lines changed: 79 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,79 @@
1+
# This model generates a representation of a phase diagram for the equilibrium
2+
# wt% water that can be stored in an upper mantle peridotite (PUM) as shown in
3+
# Figure 11 of Tian et. al., 2019. This model reproduces that phase diagram
4+
# (very low res by default), with an important caveat being that the pressure axis
5+
# is inverted relative to the ASPECT output. A surface pressure of 600 MPa is imposed
6+
# to reproduce the pressure axis,and a lateral temperature gradient is imposed to
7+
# reproduce the temperature axis.
8+
include $ASPECT_SOURCE_DIR/tests/tian_peridotite.prm
9+
set Nonlinear solver scheme = iterated Advection and Stokes
10+
# set Surface pressure = 0.6e9
11+
12+
# The Tian approximation implementation requires that all 4 rock compositions exist
13+
# as compositional fields, though they do not have to all be non zero.
14+
subsection Compositional fields
15+
set Number of fields = 10
16+
set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, peridotite_porosity, gabbro_porosity, MORB_porosity, sediment_porosity
17+
set Compositional field methods = darcy field, field, field, field, field, field, darcy field, darcy field, darcy field, darcy field
18+
set Types of fields = porosity, generic, chemical composition, chemical composition, chemical composition, chemical composition, porosity, porosity, porosity, porosity
19+
end
20+
21+
# Specify that the entire domain is peridotite, with an initial bound water content
22+
# of 11% everywhere, the bound_water will evolve to the parametrized phase diagram
23+
# for peridotite.
24+
subsection Initial composition model
25+
set Model name = function
26+
27+
subsection Function
28+
set Function constants = initial_porosity=0.0, initial_bound_water=0.11
29+
set Function expression = initial_porosity; \
30+
initial_bound_water; \
31+
0.4; \
32+
0.3; \
33+
0.2; \
34+
0.1; \
35+
0; \
36+
0; \
37+
0; \
38+
0
39+
end
40+
end
41+
42+
# No gravity
43+
subsection Gravity model
44+
set Model name = vertical
45+
46+
subsection Vertical
47+
set Magnitude = 0.0
48+
end
49+
end
50+
51+
# Material model, using the simpler base model and selecting the tian approximation
52+
# as the fluid solid reaction scheme for the reactive fluid transport model.
53+
subsection Material model
54+
set Model name = reactive fluid transport
55+
56+
subsection Reactive Fluid Transport Model
57+
set Base model = visco plastic
58+
set Reference fluid density = 2000
59+
subsection Tian 2019 model
60+
set Track source of free water = true
61+
set Maximum weight percent water in peridotite = 11
62+
set Maximum weight percent water in gabbro = 11
63+
set Maximum weight percent water in MORB = 11
64+
set Maximum weight percent water in sediment = 11
65+
end
66+
end
67+
68+
subsection Visco Plastic
69+
set Viscous flow law = diffusion
70+
set Densities = 3000, 3000, 3000, 3300, 3100, 2900, 2300, 3000, 3000, 3000, 3000
71+
set Thermal expansivities = 0
72+
set Minimum viscosity = 1e21
73+
set Maximum viscosity = 1e21
74+
end
75+
end
76+
77+
subsection Melt settings
78+
set Include melt transport = false
79+
end

‎tests/tian_track_fluid_sources/screen-output‎

Lines changed: 93 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,93 @@
1+
-----------------------------------------------------------------------------
2+
-- This is ASPECT --
3+
-- The Advanced Solver for Planetary Evolution, Convection, and Tectonics. --
4+
-----------------------------------------------------------------------------
5+
-- . version 3.1.0-pre (track_fluid_source, 986c9a918)
6+
-- . using deal.II 9.6.2
7+
-- . with 32 bit indices
8+
-- . with vectorization level 1 (SSE2, 2 doubles, 128 bits)
9+
-- . using Trilinos 14.4.0
10+
-- . using p4est 2.3.6
11+
-- . using Geodynamic World Builder 1.0.0
12+
-- . running in DEBUG mode
13+
-- . running with 1 MPI process
14+
-----------------------------------------------------------------------------
15+
16+
-----------------------------------------------------------------------------
17+
-- For information on how to cite ASPECT, see:
18+
-- https://aspect.geodynamics.org/citing.html?ver=3.1.0-pre&mf=1&sha=986c9a918&src=code
19+
-----------------------------------------------------------------------------
20+
Number of active cells: 16 (on 3 levels)
21+
Number of degrees of freedom: 1,078 (162+25+81+81+81+81+81+81+81+81+81+81+81)
22+
23+
*** Timestep 0: t=0 years, dt=0 years
24+
Solving temperature system... 0 iterations.
25+
Skipping porosity composition solve because RHS is zero.
26+
Solving bound_fluid system ... 0 iterations.
27+
Solving peridotite system ... 0 iterations.
28+
Solving gabbro system ... 0 iterations.
29+
Solving MORB system ... 0 iterations.
30+
Solving sediment system ... 0 iterations.
31+
Skipping peridotite_porosity composition solve because RHS is zero.
32+
Skipping gabbro_porosity composition solve because RHS is zero.
33+
Skipping MORB_porosity composition solve because RHS is zero.
34+
Skipping sediment_porosity composition solve because RHS is zero.
35+
Solving Stokes system (GMG)... 0+0 iterations.
36+
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.41953e-16, 0, 1.14893e-16, 1.43204e-16, 8.69635e-17, 1.43204e-16, 1.43204e-16, 0, 0, 0, 0, 0
37+
Relative nonlinear residual (total system) after nonlinear iteration 1: 1.43204e-16
38+
39+
40+
Postprocessing:
41+
Compositions min/max/mass: 0/0/0 // 0.11/0.11/4.4e+09 // 0.4/0.4/1.6e+10 // 0.3/0.3/1.2e+10 // 0.2/0.2/8e+09 // 0.1/0.1/4e+09 // 0/0/0 // 0/0/0 // 0/0/0 // 0/0/0
42+
43+
*** Timestep 1: t=1000 years, dt=1000 years
44+
Solving composition reactions... in 100 substep(s).
45+
Solving temperature system... 6 iterations.
46+
Solving porosity system ... 0 iterations.
47+
Solving bound_fluid system ... 0 iterations.
48+
Solving peridotite system ... 0 iterations.
49+
Solving gabbro system ... 0 iterations.
50+
Solving MORB system ... 0 iterations.
51+
Solving sediment system ... 0 iterations.
52+
Solving peridotite_porosity system ... 0 iterations.
53+
Solving gabbro_porosity system ... 0 iterations.
54+
Solving MORB_porosity system ... 0 iterations.
55+
Solving sediment_porosity system ... 0 iterations.
56+
Solving Stokes system (GMG)... 0+0 iterations.
57+
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.15166e-06, 8.87349e-17, 1.26557e-16, 1.43204e-16, 8.69635e-17, 1.43204e-16, 1.43204e-16, 1.39448e-16, 1.39812e-16, 1.39448e-16, 1.39448e-16, 0
58+
Relative nonlinear residual (total system) after nonlinear iteration 1: 1.15166e-06
59+
60+
61+
Postprocessing:
62+
Compositions min/max/mass: 0.02515/0.1554/5.079e+09 // 0.00301/0.09499/9.701e+08 // 0.4/0.4/1.6e+10 // 0.3/0.3/1.2e+10 // 0.2/0.2/8e+09 // 0.1/0.1/4e+09 // 0.01006/0.06214/2.032e+09 // 0.007546/0.04661/1.524e+09 // 0.00503/0.03107/1.016e+09 // 0.002515/0.01554/5.079e+08
63+
64+
Termination requested by criterion: end time
65+
66+
67+
+----------------------------------------------+------------+------------+
68+
| Total wallclock time elapsed since start | 0.298s | |
69+
| | | |
70+
| Section | no. calls | wall time | % of total |
71+
+----------------------------------+-----------+------------+------------+
72+
| Assemble Stokes system rhs | 2 | 0.00562s | 1.9% |
73+
| Assemble composition system | 20 | 0.071s | 24% |
74+
| Assemble temperature system | 2 | 0.00398s | 1.3% |
75+
| Build Stokes preconditioner | 2 | 0.000314s | 0.11% |
76+
| Build composition preconditioner | 15 | 0.00126s | 0.42% |
77+
| Build temperature preconditioner | 2 | 0.000281s | 0% |
78+
| Initialization | 1 | 0.122s | 41% |
79+
| Postprocessing | 2 | 0.00138s | 0.46% |
80+
| Setup dof systems | 1 | 0.00306s | 1% |
81+
| Setup initial conditions | 1 | 0.00475s | 1.6% |
82+
| Setup matrices | 1 | 0.0031s | 1% |
83+
| Solve Stokes system | 2 | 0.00204s | 0.69% |
84+
| Solve composition reactions | 1 | 0.0676s | 23% |
85+
| Solve composition system | 15 | 0.00337s | 1.1% |
86+
| Solve temperature system | 2 | 0.000535s | 0.18% |
87+
+----------------------------------+-----------+------------+------------+
88+
89+
-- Total wallclock time elapsed including restarts: 0s
90+
-----------------------------------------------------------------------------
91+
-- For information on how to cite ASPECT, see:
92+
-- https://aspect.geodynamics.org/citing.html?ver=3.1.0-pre&mf=1&sha=986c9a918&src=code
93+
-----------------------------------------------------------------------------

‎tests/tian_track_fluid_sources/statistics‎

Lines changed: 54 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,54 @@
1+
# 1: Time step number
2+
# 2: Time (years)
3+
# 3: Time step size (years)
4+
# 4: Number of mesh cells
5+
# 5: Number of Stokes degrees of freedom
6+
# 6: Number of temperature degrees of freedom
7+
# 7: Number of degrees of freedom for all compositions
8+
# 8: Number of nonlinear iterations
9+
# 9: Iterations for temperature solver
10+
# 10: Iterations for composition solver 1
11+
# 11: Iterations for composition solver 2
12+
# 12: Iterations for composition solver 3
13+
# 13: Iterations for composition solver 4
14+
# 14: Iterations for composition solver 5
15+
# 15: Iterations for composition solver 6
16+
# 16: Iterations for composition solver 7
17+
# 17: Iterations for composition solver 8
18+
# 18: Iterations for composition solver 9
19+
# 19: Iterations for composition solver 10
20+
# 20: Iterations for Stokes solver
21+
# 21: Velocity iterations in Stokes preconditioner
22+
# 22: Schur complement iterations in Stokes preconditioner
23+
# 23: Minimal value for composition porosity
24+
# 24: Maximal value for composition porosity
25+
# 25: Global mass for composition porosity
26+
# 26: Minimal value for composition bound_fluid
27+
# 27: Maximal value for composition bound_fluid
28+
# 28: Global mass for composition bound_fluid
29+
# 29: Minimal value for composition peridotite
30+
# 30: Maximal value for composition peridotite
31+
# 31: Global mass for composition peridotite
32+
# 32: Minimal value for composition gabbro
33+
# 33: Maximal value for composition gabbro
34+
# 34: Global mass for composition gabbro
35+
# 35: Minimal value for composition MORB
36+
# 36: Maximal value for composition MORB
37+
# 37: Global mass for composition MORB
38+
# 38: Minimal value for composition sediment
39+
# 39: Maximal value for composition sediment
40+
# 40: Global mass for composition sediment
41+
# 41: Minimal value for composition peridotite_porosity
42+
# 42: Maximal value for composition peridotite_porosity
43+
# 43: Global mass for composition peridotite_porosity
44+
# 44: Minimal value for composition gabbro_porosity
45+
# 45: Maximal value for composition gabbro_porosity
46+
# 46: Global mass for composition gabbro_porosity
47+
# 47: Minimal value for composition MORB_porosity
48+
# 48: Maximal value for composition MORB_porosity
49+
# 49: Global mass for composition MORB_porosity
50+
# 50: Minimal value for composition sediment_porosity
51+
# 51: Maximal value for composition sediment_porosity
52+
# 52: Global mass for composition sediment_porosity
53+
0 0.000000000000e+00 0.000000000000e+00 16 187 81 810 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.00000000e+00 0.00000000e+00 0.00000000e+00 1.10000000e-01 1.10000000e-01 4.40000000e+09 4.00000000e-01 4.00000000e-01 1.60000000e+10 3.00000000e-01 3.00000000e-01 1.20000000e+10 2.00000000e-01 2.00000000e-01 8.00000000e+09 1.00000000e-01 1.00000000e-01 4.00000000e+09 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
54+
1 1.000000000000e+03 1.000000000000e+03 16 187 81 810 1 6 0 0 0 0 0 0 0 0 0 0 0 0 0 2.51518787e-02 1.55352795e-01 5.07898344e+09 3.01021900e-03 9.49916812e-02 9.70052691e+08 4.00000000e-01 4.00000000e-01 1.60000000e+10 3.00000000e-01 3.00000000e-01 1.20000000e+10 2.00000000e-01 2.00000000e-01 8.00000000e+09 1.00000000e-01 1.00000000e-01 4.00000000e+09 1.00607515e-02 6.21411181e-02 2.03159337e+09 7.54556360e-03 4.66058386e-02 1.52369503e+09 5.03037573e-03 3.10705590e-02 1.01579669e+09 2.51518787e-03 1.55352795e-02 5.07898344e+08

0 commit comments

Comments
 (0)