ASPECT
visco_plastic.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 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_visco_plastic_h
22 #define _aspect_material_model_visco_plastic_h
23 
28 
29 #include<deal.II/fe/component_mask.h>
30 
31 namespace aspect
32 {
33  namespace MaterialModel
34  {
183  template <int dim>
185  {
186  public:
191  void
192  initialize () override;
193 
195  MaterialModel::MaterialModelOutputs<dim> &out) const override;
196 
207  bool is_compressible () const override;
208 
209  static
210  void
211  declare_parameters (ParameterHandler &prm);
212 
213  void
214  parse_parameters (ParameterHandler &prm) override;
215 
216  void
218 
219  double get_min_strain_rate() const;
220 
226  bool
227  is_yielding (const MaterialModelInputs<dim> &in) const;
228 
229  private:
230 
242  std::unique_ptr<Rheology::ViscoPlastic<dim>> rheology;
243 
244  std::vector<double> thermal_diffusivities;
245 
250 
251  std::vector<double> thermal_conductivities;
252 
257 
261  unsigned int n_phases;
262 
267 
272 
277 
281  std::unique_ptr<MaterialUtilities::PhaseFunctionDiscrete<dim>> phase_function_discrete;
282 
283  };
284 
285  }
286 }
287 
288 #endif
MaterialUtilities::PhaseFunction< dim > phase_function
std::vector< unsigned int > n_phase_transitions_for_each_chemical_composition
bool is_yielding(const MaterialModelInputs< dim > &in) const
std::unique_ptr< Rheology::ViscoPlastic< dim > > rheology
static void declare_parameters(ParameterHandler &prm)
void parse_parameters(ParameterHandler &prm) override
EquationOfState::MulticomponentIncompressible< dim > equation_of_state
bool is_compressible() const override
std::vector< double > thermal_diffusivities
void evaluate(const MaterialModel::MaterialModelInputs< dim > &in, MaterialModel::MaterialModelOutputs< dim > &out) const override
void create_additional_named_outputs(MaterialModel::MaterialModelOutputs< dim > &out) const override
std::unique_ptr< MaterialUtilities::PhaseFunctionDiscrete< dim > > phase_function_discrete
std::vector< double > thermal_conductivities