ASPECT
steinberger.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 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_steinberger_h
22 #define _aspect_material_model_steinberger_h
23 
26 
28 #include <deal.II/fe/component_mask.h>
29 
30 namespace aspect
31 {
32  namespace MaterialModel
33  {
34  using namespace dealii;
35 
36  namespace internal
37  {
47  {
48  public:
52  LateralViscosityLookup(const std::string &filename,
53  const MPI_Comm &comm);
54 
58  double lateral_viscosity(double depth) const;
59 
63  int get_nslices() const;
64  private:
68  std::vector<double> values;
69 
73  double min_depth;
74  double delta_depth;
75  double max_depth;
76  };
77 
86  {
87  public:
91  RadialViscosityLookup(const std::string &filename,
92  const MPI_Comm &comm);
93 
97  double radial_viscosity(double depth) const;
98 
99  private:
103  std::vector<double> values;
104 
108  double min_depth;
109  double delta_depth;
110  double max_depth;
111  };
112  }
113 
126  template <int dim>
128  {
129  public:
134  void
135  initialize () override;
136 
141  void update() override;
142 
147  virtual double viscosity (const double temperature,
148  const double pressure,
149  const std::vector<double> &compositional_fields,
150  const SymmetricTensor<2,dim> &strain_rate,
151  const Point<dim> &position) const;
169  bool is_compressible () const override;
178  void
179  evaluate(const MaterialModel::MaterialModelInputs<dim> &in,
180  MaterialModel::MaterialModelOutputs<dim> &out) const override;
181 
189  static
190  void
191  declare_parameters (ParameterHandler &prm);
192 
196  void
197  parse_parameters (ParameterHandler &prm) override;
202  void
203  create_additional_named_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const override;
204 
205 
206  private:
216  double thermal_conductivity (const double temperature,
217  const double pressure,
218  const Point<dim> &position) const;
219 
230 
238  std::unique_ptr<ComponentMask> composition_mask;
239 
244 
251 
257 
262  {
264  p_T_dependent
265  } conductivity_formulation;
266 
271  std::vector<double> conductivity_transition_depths;
272  std::vector<double> reference_thermal_conductivities;
275  std::vector<double> conductivity_exponents;
276  std::vector<double> saturation_scaling;
278 
282  std::vector<double> average_temperature;
283  unsigned int n_lateral_slices;
284 
289  double min_eta;
290  double max_eta;
292 
296  std::string data_directory;
299 
304  std::unique_ptr<internal::LateralViscosityLookup> lateral_viscosity_lookup;
305 
310  std::unique_ptr<internal::RadialViscosityLookup> radial_viscosity_lookup;
311 
318  void fill_prescribed_outputs (const unsigned int i,
319  const std::vector<double> &volume_fractions,
322 
323  };
324  }
325 }
326 
327 #endif
std::vector< double > conductivity_reference_temperatures
Definition: steinberger.h:274
std::vector< double > conductivity_transition_depths
Definition: steinberger.h:271
std::unique_ptr< internal::LateralViscosityLookup > lateral_viscosity_lookup
Definition: steinberger.h:304
std::vector< double > average_temperature
Definition: steinberger.h:282
std::unique_ptr< ComponentMask > composition_mask
Definition: steinberger.h:238
void declare_parameters(ParameterHandler &prm)
EquationOfState::ThermodynamicTableLookup< dim > equation_of_state
Definition: steinberger.h:243
std::vector< double > saturation_scaling
Definition: steinberger.h:276
std::vector< double > conductivity_pressure_dependencies
Definition: steinberger.h:273
std::vector< double > reference_thermal_conductivities
Definition: steinberger.h:272
std::vector< double > conductivity_exponents
Definition: steinberger.h:275
Definition: compat.h:88
std::unique_ptr< internal::RadialViscosityLookup > radial_viscosity_lookup
Definition: steinberger.h:310