ASPECT
dislocation_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_dislocation_creep_h
22 #define _aspect_material_model_rheology_dislocation_creep_h
23 
24 #include <aspect/global.h>
26 
27 namespace aspect
28 {
29  namespace MaterialModel
30  {
31  using namespace dealii;
32 
33  namespace Rheology
34  {
39  {
44  double prefactor;
48  };
49 
50  template <int dim>
52  {
53  public:
58 
62  static
63  void
64  declare_parameters (ParameterHandler &prm);
65 
73  void
74  parse_parameters (ParameterHandler &prm,
75  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
76 
85  compute_creep_parameters (const unsigned int composition,
86  const std::vector<double> &phase_function_values = std::vector<double>(),
87  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
88 
96  double
97  compute_viscosity (const double strain_rate,
98  const double pressure,
99  const double temperature,
100  const unsigned int composition,
101  const std::vector<double> &phase_function_values = std::vector<double>(),
102  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
103 
112  std::pair<double, double>
113  compute_strain_rate_and_derivative (const double stress,
114  const double pressure,
115  const double temperature,
116  const DislocationCreepParameters creep_parameters) const;
117 
118  private:
119 
123  std::vector<double> prefactors_dislocation;
124 
128  std::vector<double> stress_exponents_dislocation;
129 
133  std::vector<double> activation_energies_dislocation;
134 
138  std::vector<double> activation_volumes_dislocation;
139 
140  };
141  }
142  }
143 }
144 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:88