ASPECT
diffusion_creep.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2021 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 #ifndef _aspect_material_model_rheology_diffusion_creep_h
22 #define _aspect_material_model_rheology_diffusion_creep_h
23 
24 #include <aspect/global.h>
28 
29 namespace aspect
30 {
31  namespace MaterialModel
32  {
33  using namespace dealii;
34 
35  namespace Rheology
36  {
41  {
46  double prefactor;
51  };
52 
53  template <int dim>
55  {
56  public:
61 
65  static
66  void
67  declare_parameters (ParameterHandler &prm);
68 
76  void
77  parse_parameters (ParameterHandler &prm,
78  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
79 
80 
89  compute_creep_parameters (const unsigned int composition,
90  const std::vector<double> &phase_function_values = std::vector<double>(),
91  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
92 
100  double
101  compute_viscosity (const double pressure,
102  const double temperature,
103  const unsigned int composition,
104  const std::vector<double> &phase_function_values = std::vector<double>(),
105  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
106 
115  std::pair<double, double>
116  compute_strain_rate_and_derivative (const double stress,
117  const double pressure,
118  const double temperature,
119  const DiffusionCreepParameters creep_parameters) const;
120 
121  private:
122 
126  std::vector<double> prefactors_diffusion;
127 
131  std::vector<double> stress_exponents_diffusion;
132 
136  std::vector<double> grain_size_exponents_diffusion;
137 
141  std::vector<double> activation_energies_diffusion;
142 
146  std::vector<double> activation_volumes_diffusion;
147 
151  double grain_size;
152  };
153  }
154  }
155 }
156 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:88