ASPECT
grain_boundary_sliding.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2025 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_grain_boundary_sliding_h
22 #define _aspect_material_model_rheology_grain_boundary_sliding_h
23 
24 #include <aspect/global.h>
28 
29 namespace aspect
30 {
31  namespace MaterialModel
32  {
33  namespace Rheology
34  {
39  {
44  double prefactor;
49 
54  };
55 
56  template <int dim>
58  {
59  public:
64 
68  static
69  void
70  declare_parameters (ParameterHandler &prm);
71 
79  void
80  parse_parameters (ParameterHandler &prm,
81  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
82 
83 
92  compute_slide_parameters (const unsigned int composition,
93  const std::vector<double> &phase_function_values = std::vector<double>(),
94  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
95 
104  double
105  compute_viscosity (const double strain_rate,
106  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 
119  double
120  compute_viscosity (const double strain_rate,
121  const double pressure,
122  const double temperature,
123  const double grain_size,
124  const unsigned int composition,
125  const std::vector<double> &phase_function_values = std::vector<double>(),
126  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
127 
128  private:
129 
133  std::vector<double> prefactors;
134 
138  std::vector<double> stress_exponents;
139 
143  std::vector<double> grain_size_exponents;
144 
148  std::vector<double> activation_energies;
149 
153  std::vector<double> activation_volumes;
154 
161  };
162  }
163  }
164 }
165 #endif
void declare_parameters(ParameterHandler &prm)