ASPECT
composite_visco_plastic.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2020 - 2022 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_composite_visco_plastic_h
22 #define _aspect_material_model_rheology_composite_visco_plastic_h
23 
24 #include <aspect/global.h>
32 
33 namespace aspect
34 {
35  namespace MaterialModel
36  {
37  using namespace dealii;
38 
39  namespace Rheology
40  {
41 
42  template <int dim>
44  {
45  public:
50 
54  static
55  void
56  declare_parameters (ParameterHandler &prm);
57 
65  void
66  parse_parameters (ParameterHandler &prm,
67  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
68 
76  double
77  compute_viscosity (const double pressure,
78  const double temperature,
79  const std::vector<double> &volume_fractions,
80  const SymmetricTensor<2,dim> &strain_rate,
81  std::vector<double> &partial_strain_rates,
82  const std::vector<double> &phase_function_values = std::vector<double>(),
83  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
84 
93  double
94  compute_composition_viscosity (const double pressure,
95  const double temperature,
96  const unsigned int composition,
97  const SymmetricTensor<2,dim> &strain_rate,
98  std::vector<double> &partial_strain_rates,
99  const std::vector<double> &phase_function_values = std::vector<double>(),
100  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
101 
110  std::pair<double, double>
111  compute_strain_rate_and_derivative (const double creep_stress,
112  const double pressure,
113  const double temperature,
114  const DiffusionCreepParameters diffusion_creep_parameters,
115  const DislocationCreepParameters dislocation_creep_parameters,
116  const PeierlsCreepParameters peierls_creep_parameters,
117  const DruckerPragerParameters drucker_prager_parameters) const;
118 
119  private:
120 
128 
133  std::unique_ptr<Rheology::DiffusionCreep<dim>> diffusion_creep;
134  std::unique_ptr<Rheology::DislocationCreep<dim>> dislocation_creep;
135  std::unique_ptr<Rheology::PeierlsCreep<dim>> peierls_creep;
136  std::unique_ptr<Rheology::DruckerPrager<dim>> drucker_prager;
137 
139 
143 
147  };
148  }
149  }
150 }
151 #endif
std::unique_ptr< Rheology::PeierlsCreep< dim > > peierls_creep
void declare_parameters(ParameterHandler &prm)
std::unique_ptr< Rheology::DislocationCreep< dim > > dislocation_creep
std::unique_ptr< Rheology::DiffusionCreep< dim > > diffusion_creep
std::unique_ptr< Rheology::DruckerPrager< dim > > drucker_prager
Definition: compat.h:88