Skip to content

Commit

Permalink
Update test parameter files
Browse files Browse the repository at this point in the history
  • Loading branch information
@gassmoeller
gassmoeller committed Sep 6, 2023
1 parent 582898e commit 269db98
Show file tree
Hide file tree
Showing 1,002 changed files with 6,417 additions and 19,016 deletions.
78 changes: 43 additions & 35 deletions tests/ISA_rotation_timescale.prm
View file
Original file line number Diff line number Diff line change
Expand Up @@ -8,65 +8,72 @@ set Adiabatic surface temperature = 1570

# geometry
subsection Geometry model
set Model name = box
subsection Box
set X extent = 105000
set Y extent = 70000
set X repetitions = 3
set Y repetitions = 2
end
set Model name = box

subsection Box
set X extent = 105000
set Y extent = 70000
set X repetitions = 3
set Y repetitions = 2
end
end

# gravity
subsection Gravity model
set Model name = vertical

subsection Vertical
set Magnitude = 9.81
end
end

# boundary velocity
subsection Boundary velocity model
set Prescribed velocity boundary indicators = top:function
set Tangential velocity boundary indicators = left
subsection Function
set Function constants = u0=0.05, x0=10000, xc=52500
set Variable names = x,z
set Function expression = if(x<x0,(1-(x/x0-1)*(x/x0-1))*u0,u0); 0
end
set Prescribed velocity boundary indicators = top:function
set Tangential velocity boundary indicators = left

subsection Function
set Function constants = u0=0.05, x0=10000, xc=52500
set Variable names = x,z
set Function expression = if(x<x0,(1-(x/x0-1)*(x/x0-1))*u0,u0); 0
end
end

# boundary traction
subsection Boundary traction model
set Prescribed traction boundary indicators = right:initial lithostatic pressure, bottom:initial lithostatic pressure
subsection Initial lithostatic pressure
set Representative point = 105000, 7000
end
set Prescribed traction boundary indicators = right:initial lithostatic pressure, bottom:initial lithostatic pressure

subsection Initial lithostatic pressure
set Representative point = 105000, 7000
end
end

# initial temperature
subsection Initial temperature model
set Model name = adiabatic
subsection Adiabatic
set Age top boundary layer = 1e8
set Age bottom boundary layer = 1e5
set Amplitude = 0
end
set Model name = adiabatic

subsection Adiabatic
set Age top boundary layer = 1e8
set Age bottom boundary layer = 1e5
set Amplitude = 0
end
end

# boundary temperature
subsection Boundary temperature model
set Fixed temperature boundary indicators = top, bottom
set List of model names = box
subsection Box
set Top temperature = 293
set Bottom temperature = 1570
end
set Fixed temperature boundary indicators = top, bottom
set List of model names = box

subsection Box
set Top temperature = 293
set Bottom temperature = 1570
end
end

# material
subsection Material model
set Model name = simple

subsection Simple model
set Thermal conductivity = 4.7
set Reference density = 3400
Expand All @@ -86,9 +93,10 @@ end

# postprocessing
subsection Postprocess
set List of postprocessors = visualization
subsection Visualization
set List of output variables = ISA rotation timescale
set Output format = gnuplot
end
set List of postprocessors = visualization

subsection Visualization
set List of output variables = ISA rotation timescale
set Output format = gnuplot
end
end
135 changes: 66 additions & 69 deletions tests/Picard_compressible.prm
View file
Original file line number Diff line number Diff line change
Expand Up @@ -41,73 +41,70 @@ subsection Mesh refinement
set Initial global refinement = 4
set Time steps between mesh refinement = 0
set Strategy = minimum refinement function

subsection Minimum refinement function
set Variable names = depth,t
set Coordinate system = depth
set Function expression = if(depth<=200e3, 5,if(depth<=620e3, 5,if(depth<=680e3,5,if(depth>=2586e3,5,5))))
#if(depth<560e3,0,if(depth>760e3,0,6))
# if(depth<560e3,5,if(depth>760e3,5,6))
end
end

##### Set temperature Model
subsection Boundary temperature model
#set Allow fixed temperature on outflow boundaries = true # default
set Fixed temperature boundary indicators = bottom, top

# set Fixed temperature boundary indicators = top, bottom
set List of model names = spherical constant

subsection Spherical constant
set Inner temperature = 3500
set Outer temperature = 273 #273 for surface ############
end
end

subsection Boundary velocity model
set Tangential velocity boundary indicators = inner, outer
# set Zero velocity boundary indicators = outer
end

subsection Nullspace removal
set Remove nullspace = net rotation
set Remove nullspace = net rotation
end


# Number and names of compositional fields
subsection Compositional fields
set Number of fields = 1
set Names of fields = plastic_strain #add strain field
end


##### Set temperature Model
subsection Boundary temperature model
# set Fixed temperature boundary indicators = top, bottom
set List of model names = spherical constant

subsection Spherical constant
set Inner temperature = 3500
set Outer temperature = 273 #273 for surface ############
end
end

#### Temperature model
subsection Initial temperature model
#set Model name = adiabatic
set List of model names = adiabatic, function
set List of model operators = add

subsection Adiabatic
set Age bottom boundary layer = 100.e6
set Age top boundary layer = 100.e6
set Radius = 0
set Amplitude = 0
set Position = center
set Subadiabaticity = 0
end
subsection Function
set Coordinate system = cartesian
set Variable names = x,y,z
set Function constants = R=6371.e3
set Function expression = \
if(sqrt(x*x+y*y)<=R<100e3, rand_seed(1) * 100 -50 , 0)
end
end

subsection Function
set Coordinate system = cartesian
set Variable names = x,y,z
set Function constants = R=6371.e3
set Function expression = \
if(sqrt(x*x+y*y)<=R<100e3, rand_seed(1) * 100 -50 , 0)
end
end

subsection Adiabatic conditions model # need to initialize for initial adiabatic temp.
subsection Adiabatic conditions model # need to initialize for initial adiabatic temp.
#send
set Model name = ascii data

subsection Ascii data model
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
set Data file name = example_isentrope.txt
Expand All @@ -116,28 +113,29 @@ end

subsection Heating model
set List of model names = adiabatic heating, shear heating, radioactive decay
subsection Adiabatic heating

subsection Adiabatic heating
set Use simplified adiabatic heating = true
end

subsection Radioactive decay
set Crust defined by composition = false
set Crust depth = 10e3
set Half decay times = 4.468e9, 0.704e9, 14e9, 1.248e9
set Heating rates = 9.4946e-5, 5.68e-4, 2.6368e-5, 2.8761e-5
set Initial concentrations crust = 0.0, 0.0, 0.0, 0.0
set Initial concentrations mantle =0.0082, 0.000059, 0.024, 0.013
set Initial concentrations mantle = 0.0082, 0.000059, 0.024, 0.013
set Number of elements = 4
end
#set List of model names = compositional heating
end


############################ MATERIAL MODEL ######################################################################

subsection Material model
set Model name = compositing #visco plastic #melt simple
set Material averaging = harmonic average
############################ MATERIAL MODEL ######################################################################

############################ MATERIAL MODEL ######################################################################
subsection Compositing
set Viscosity = visco plastic
set Density = ascii reference profile
Expand All @@ -149,48 +147,49 @@ subsection Material model
set Entropy derivative temperature = ascii reference profile
set Reaction terms = visco plastic
end
subsection Ascii reference profile
set Thermal conductivity = 3.0
subsection Ascii data model
set Data file name = example_isentrope.txt
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
end
end

subsection Visco Plastic
subsection Ascii reference profile
set Thermal conductivity = 3.0

subsection Ascii data model
set Data file name = example_isentrope.txt
set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
end
end

set Define transition by depth instead of pressure = true
set Phase transition Clapeyron slopes = background: 3.0e6|-2.5e6|9.e6
set Phase transition depths = background: 410e3|660e3|2740e3
set Phase transition widths = background: 12e3|23e3|11e3
set Phase transition temperatures = background: 1770|1852|2591
subsection Visco Plastic
set Define transition by depth instead of pressure = true
set Phase transition Clapeyron slopes = background: 3.0e6|-2.5e6|9.e6
set Phase transition depths = background: 410e3|660e3|2740e3
set Phase transition widths = background: 12e3|23e3|11e3
set Phase transition temperatures = background: 1770|1852|2591

# Reference temperature and viscosity
set Reference temperature = 293
# Reference temperature and viscosity
set Reference temperature = 293

# The minimum strain-rate helps limit large viscosities values that arise
# as the strain-rate approaches zero.
# The reference strain-rate is used on the first non-linear iteration
# of the first time step when the velocity has not been determined yet.
set Minimum strain rate = 1.e-20
set Reference strain rate = 1.e-15
# The minimum strain-rate helps limit large viscosities values that arise
# as the strain-rate approaches zero.
# The reference strain-rate is used on the first non-linear iteration
# of the first time step when the velocity has not been determined yet.
set Minimum strain rate = 1.e-20
set Reference strain rate = 1.e-15

# Limit the viscosity with minimum and maximum values
set Minimum viscosity = 1e20
set Maximum viscosity = 1e26
# Limit the viscosity with minimum and maximum values
set Minimum viscosity = 1e20
set Maximum viscosity = 1e26

# Harmonic viscosity averaging
set Viscosity averaging scheme = harmonic
# Harmonic viscosity averaging
set Viscosity averaging scheme = harmonic

# Choose to have the viscosity (pre-yield) follow a dislocation
# diffusion or composite flow law. Here, dislocation is selected
# so no need to specify diffusion creep parameters below, which are
# only used if "diffusion" or "composite" option is selected.
set Viscous flow law = diffusion
# Choose to have the viscosity (pre-yield) follow a dislocation
# diffusion or composite flow law. Here, dislocation is selected
# so no need to specify diffusion creep parameters below, which are
# only used if "diffusion" or "composite" option is selected.
set Viscous flow law = diffusion

# from cizkov
# order: background , mantle_low , mantle_up
# set Thermal diffusivities = 8.384146341463415e-07, 8.384146341463415e-07 , 7.309941520467837e-07
# from cizkov
# order: background , mantle_low , mantle_up
# set Thermal diffusivities = 8.384146341463415e-07, 8.384146341463415e-07 , 7.309941520467837e-07

set Thermal diffusivities = 1.e-6
set Thermal expansivities = background: 3.e-5|2.e-5|2.e-5|2.e-5, plastic_strain: 3.e-5
Expand All @@ -207,18 +206,15 @@ subsection Material model
set Grain size = 1 #1e-3
set Grain size exponents for diffusion creep = 0 #default

### # Plasticity parameters
### # Plasticity parameters
set Angles of internal friction = background: 0|0|0|0, plastic_strain:0 # Angle Mohr
set Cohesions = background: 40e6|1e20|1e20|1e20, plastic_strain:40e6

set Strain weakening mechanism = plastic weakening with plastic strain only
set Start plasticity strain weakening intervals = 1.e30
set Cohesion strain weakening factors = 1
set Friction strain weakening factors = 1
set End plasticity strain weakening intervals = 2.e30

end

end

#######################
Expand All @@ -236,6 +232,7 @@ end
subsection Postprocess
set List of postprocessors = velocity statistics, basic statistics, temperature statistics, \
composition statistics, depth average, heat flux statistics

subsection Depth average
set Number of zones = 3
set Output format = txt
Expand Down

0 comments on commit 269db98

Please sign in to comment.