ASPECT
frank_kamenetskii.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2020 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_frank_kamenetskii_h
22 #define _aspect_material_model_rheology_frank_kamenetskii_h
23 
24 #include <aspect/global.h>
27 
28 namespace aspect
29 {
30  namespace MaterialModel
31  {
32  using namespace dealii;
33 
34  namespace Rheology
35  {
47  template <int dim>
49  {
50  public:
55 
59  static
60  void
61  declare_parameters (ParameterHandler &prm);
62 
66  void
67  parse_parameters (ParameterHandler &prm);
68 
72  double
73  compute_viscosity (const double temperature,
74  const unsigned int composition,
75  const double pressure = std::numeric_limits<double>::infinity(),
76  const double density = std::numeric_limits<double>::infinity(),
77  const double gravity = std::numeric_limits<double>::infinity()) const;
78 
79  private:
84 
88  std::vector<double> prefactors_frank_kamenetskii;
89 
94 
95  std::vector<double> reference_temperatures;
96  std::vector<double> reference_pressures;
97  };
98  }
99  }
100 }
101 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:59
Definition: compat.h:42