-
Notifications
You must be signed in to change notification settings - Fork 231
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
fix freezing in melt simple material model and improve documentation
- Loading branch information
Showing
7 changed files
with
316 additions
and
22 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,6 @@ | ||
Fixed: The freezing of melt in the melt simple material model is | ||
now more consistent: The freezing rate parameter determines the | ||
freezing of all melt, and melt freezes if both the porosity and | ||
the depletion are larger than the equilibrium melt fraction. | ||
<br> | ||
(Juliane Dannberg, 2018/03/16) |
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
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
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,124 @@ | ||
######################################################### | ||
# This is a test for the melting rate functionality, | ||
# in particular the freezing of melt. | ||
# It uses the melt simple model. | ||
|
||
# The initial porosity has the form of a Gaussian function | ||
# there is a fixed (constant) freezing rate. | ||
|
||
# The amplitude of the Gaussian is 0.2, and the melting time | ||
# scale is chosen in such a way that the fraction of melt | ||
# should be cut in half every year (visible in the maximum | ||
# value of the porosity field). | ||
|
||
############### Global parameters | ||
|
||
set Dimension = 2 | ||
set Start time = 0 | ||
set End time = 2 | ||
set Use years in output instead of seconds = true | ||
set Use direct solver for Stokes system = false | ||
set Adiabatic surface temperature = 1600 | ||
set Nonlinear solver scheme = iterated IMPES | ||
set Max nonlinear iterations = 50 | ||
set Maximum time step = 1 | ||
|
||
set Pressure normalization = surface | ||
set Surface pressure = 0 | ||
|
||
subsection Melt settings | ||
set Melt transport threshold = 0.0 | ||
end | ||
|
||
subsection Material model | ||
|
||
set Model name = melt simple | ||
subsection Melt simple | ||
set Melt extraction depth = 0.0 | ||
set Freezing rate = 1.0 | ||
set Melting time scale for operator splitting = 1.442695041 | ||
end | ||
end | ||
|
||
set Use operator splitting = true | ||
subsection Solver parameters | ||
subsection Operator splitting parameters | ||
set Reaction time step = 0.005 | ||
set Reaction time steps per advection step = 1 | ||
end | ||
end | ||
|
||
############### Parameters describing the model | ||
subsection Geometry model | ||
set Model name = box | ||
|
||
subsection Box | ||
set X extent = 100000 | ||
set Y extent = 100000 | ||
end | ||
end | ||
|
||
subsection Model settings | ||
set Tangential velocity boundary indicators = 0, 1, 2, 3 | ||
set Include melt transport = true | ||
end | ||
|
||
subsection Compositional fields | ||
set Number of fields = 2 | ||
set Names of fields = porosity, peridotite | ||
end | ||
|
||
subsection Gravity model | ||
set Model name = vertical | ||
subsection Vertical | ||
set Magnitude = 10.0 | ||
end | ||
end | ||
|
||
subsection Initial temperature model | ||
set Model name = function | ||
subsection Function | ||
set Function expression = 293 | ||
end | ||
end | ||
|
||
subsection Initial composition model | ||
set Model name = function | ||
|
||
subsection Function | ||
set Variable names = x,y | ||
set Function constants = pi=3.1415926,x0=50000,y0=50000,c=10000 | ||
set Function expression = 0.2 * exp(-((x-x0)*(x-x0)+(y-y0)*(y-y0))/(2*c*c)); 0.0 * exp(-((x-x0)*(x-x0)+(y-y0)*(y-y0))/(2*c*c)) | ||
end | ||
end | ||
|
||
|
||
# The final part of this input file describes how many times the | ||
# mesh is refined and what to do with the solution once computed | ||
subsection Mesh refinement | ||
set Initial adaptive refinement = 0 | ||
set Initial global refinement = 5 | ||
set Time steps between mesh refinement = 0 | ||
end | ||
|
||
subsection Postprocess | ||
|
||
set List of postprocessors = composition statistics, visualization | ||
|
||
subsection Visualization | ||
|
||
set List of output variables = density, viscosity, thermal expansivity, melt material properties | ||
|
||
subsection Melt material properties | ||
set List of properties = permeability, fluid density, compaction viscosity, fluid viscosity | ||
end | ||
|
||
set Number of grouped files = 0 | ||
|
||
set Output format = vtu | ||
|
||
set Time between graphical output = 0 | ||
|
||
end | ||
|
||
end |
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,85 @@ | ||
|
||
Number of active cells: 1,024 (on 6 levels) | ||
Number of degrees of freedom: 32,842 (8,450+2,178+8,450+1,089+4,225+4,225+4,225) | ||
|
||
*** Timestep 0: t=0 years | ||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 0 iterations. | ||
Skipping peridotite composition solve because RHS is zero. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 200+665 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.76497e-16, 1.44358e-16, 0, 0.168633 | ||
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.168633 | ||
|
||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 0 iterations. | ||
Skipping peridotite composition solve because RHS is zero. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 0+0 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.76497e-16, 1.44358e-16, 0, 9.99026e-08 | ||
Relative nonlinear residual (total system) after nonlinear iteration 2: 9.99026e-08 | ||
|
||
|
||
Postprocessing: | ||
Compositions min/max/mass: 2.778e-12/0.2/1.257e+08 // 0/0/0 | ||
Writing graphical output: output-rising_melt_blob_freezing/solution/solution-00000 | ||
|
||
*** Timestep 1: t=1 years | ||
Solving composition reactions in 200 substep(s). | ||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 6 iterations. | ||
Solving peridotite system ... 5 iterations. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 200+726 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.79631e-16, 3.76619e-06, 2.26327e-07, 0.00840578 | ||
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00840578 | ||
|
||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 6 iterations. | ||
Solving peridotite system ... 4 iterations. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 10+0 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.73728e-16, 1.88328e-06, 1.38447e-07, 1.08144e-07 | ||
Relative nonlinear residual (total system) after nonlinear iteration 2: 1.88328e-06 | ||
|
||
|
||
Postprocessing: | ||
Compositions min/max/mass: -7.531e-11/0.09988/6.276e+07 // -0.1252/-1.39e-12/-7.019e+07 | ||
Writing graphical output: output-rising_melt_blob_freezing/solution/solution-00001 | ||
|
||
*** Timestep 2: t=2 years | ||
Solving composition reactions in 200 substep(s). | ||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 6 iterations. | ||
Solving peridotite system ... 5 iterations. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 200+601 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.77111e-16, 3.79756e-06, 7.12047e-08, 0.0313658 | ||
Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0313658 | ||
|
||
Solving temperature system... 0 iterations. | ||
Solving porosity system ... 6 iterations. | ||
Solving peridotite system ... 3 iterations. | ||
Rebuilding Stokes preconditioner... | ||
Solving Stokes system... 8+0 iterations. | ||
Solving for u_f in 12 iterations. | ||
Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.75463e-16, 1.4052e-06, 2.21463e-08, 1.42281e-07 | ||
Relative nonlinear residual (total system) after nonlinear iteration 2: 1.4052e-06 | ||
|
||
|
||
Postprocessing: | ||
Compositions min/max/mass: -8.803e-11/0.04988/3.134e+07 // -0.1877/3.631e-11/-1.052e+08 | ||
Writing graphical output: output-rising_melt_blob_freezing/solution/solution-00002 | ||
|
||
Termination requested by criterion: end time | ||
|
||
|
||
+---------------------------------------------+------------+------------+ | ||
+---------------------------------+-----------+------------+------------+ | ||
+---------------------------------+-----------+------------+------------+ | ||
|
Oops, something went wrong.