ASPECT
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
steinberger.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_steinberger_h
22 #define _aspect_material_model_steinberger_h
23 
29 
31 #include <deal.II/fe/component_mask.h>
32 
33 namespace aspect
34 {
35  namespace MaterialModel
36  {
37  namespace internal
38  {
48  {
49  public:
53  LateralViscosityLookup(const std::string &filename,
54  const MPI_Comm comm);
55 
59  double lateral_viscosity(double depth) const;
60 
64  int get_nslices() const;
65  private:
69  std::vector<double> values;
70 
74  double min_depth;
75  double delta_depth;
76  double max_depth;
77  };
78 
87  {
88  public:
92  RadialViscosityLookup(const std::string &filename,
93  const MPI_Comm comm);
94 
98  double radial_viscosity(double depth) const;
99 
100  private:
104  std::vector<double> values;
105 
109  double min_depth;
110  double delta_depth;
111  double max_depth;
112  };
113  }
114 
128  template <int dim>
130  {
131  public:
136  void
137  initialize () override;
138 
143  void update() override;
144 
158  bool is_compressible () const override;
167  void
168  evaluate(const MaterialModel::MaterialModelInputs<dim> &in,
169  MaterialModel::MaterialModelOutputs<dim> &out) const override;
170 
178  static
179  void
180  declare_parameters (ParameterHandler &prm);
181 
185  void
186  parse_parameters (ParameterHandler &prm) override;
191  void
192  create_additional_named_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const override;
193 
194 
195  private:
203  double viscosity (const unsigned int q,
204  const std::vector<double> &volume_fractions,
207 
218 
226  std::unique_ptr<ComponentMask> composition_mask;
227 
232 
239 
245  std::unique_ptr<ThermalConductivity::Interface<dim>> thermal_conductivity;
246 
252  std::vector<double> viscosity_prefactors;
254 
258  std::vector<double> average_temperature;
259  unsigned int n_lateral_slices;
260 
265  double min_eta;
266  double max_eta;
268 
272  std::string data_directory;
275 
280  std::unique_ptr<internal::LateralViscosityLookup> lateral_viscosity_lookup;
281 
286  std::unique_ptr<internal::RadialViscosityLookup> radial_viscosity_lookup;
287 
294  void fill_prescribed_outputs (const unsigned int i,
295  const std::vector<double> &volume_fractions,
298 
306 
307  };
308  }
309 }
310 
311 #endif
Rheology::DruckerPrager< dim > drucker_prager_plasticity
Definition: steinberger.h:303
MaterialUtilities::CompositionalAveragingOperation viscosity_averaging_scheme
Definition: steinberger.h:253
std::unique_ptr< internal::LateralViscosityLookup > lateral_viscosity_lookup
Definition: steinberger.h:280
std::vector< double > average_temperature
Definition: steinberger.h:258
Rheology::DruckerPragerParameters drucker_prager_parameters
Definition: steinberger.h:305
std::unique_ptr< ComponentMask > composition_mask
Definition: steinberger.h:226
void declare_parameters(ParameterHandler &prm)
EquationOfState::ThermodynamicTableLookup< dim > equation_of_state
Definition: steinberger.h:231
std::unique_ptr< ThermalConductivity::Interface< dim > > thermal_conductivity
Definition: steinberger.h:245
LateralViscosityLookup(const std::string &filename, const MPI_Comm comm)
std::vector< double > viscosity_prefactors
Definition: steinberger.h:252
std::unique_ptr< internal::RadialViscosityLookup > radial_viscosity_lookup
Definition: steinberger.h:286