ASPECT
diffusion_creep.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2023 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 
56  };
57 
58  template <int dim>
60  {
61  public:
66 
70  static
71  void
72  declare_parameters (ParameterHandler &prm);
73 
81  void
82  parse_parameters (ParameterHandler &prm,
83  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
84 
85 
94  compute_creep_parameters (const unsigned int composition,
95  const std::vector<double> &phase_function_values = std::vector<double>(),
96  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
97 
105  double
106  compute_viscosity (const double pressure,
107  const double temperature,
108  const unsigned int composition,
109  const std::vector<double> &phase_function_values = std::vector<double>(),
110  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
111 
116  std::pair<double, double>
117  compute_strain_rate_and_derivative (const double stress,
118  const double pressure,
119  const double temperature,
120  const DiffusionCreepParameters creep_parameters) const;
121 
126  std::pair<double, double>
127  compute_log_strain_rate_and_derivative (const double log_stress,
128  const double pressure,
129  const double temperature,
130  const DiffusionCreepParameters creep_parameters) const;
131 
132  private:
133 
137  std::vector<double> prefactors_diffusion;
138 
142  std::vector<double> stress_exponents_diffusion;
143 
147  std::vector<double> grain_size_exponents_diffusion;
148 
152  std::vector<double> activation_energies_diffusion;
153 
157  std::vector<double> activation_volumes_diffusion;
158 
162  double grain_size;
163  };
164  }
165  }
166 }
167 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:59
Definition: compat.h:42