ASPECT
melt_global.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2015 - 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_melt_global_h
22 #define _aspect_material_model_melt_global_h
23 
27 #include <aspect/melt.h>
28 
29 namespace aspect
30 {
31  namespace MaterialModel
32  {
33  using namespace dealii;
34 
46  template <int dim>
48  {
49  public:
58  bool is_compressible () const override;
59 
60  void evaluate(const typename Interface<dim>::MaterialModelInputs &in,
61  typename Interface<dim>::MaterialModelOutputs &out) const override;
62 
71  void melt_fractions (const MaterialModel::MaterialModelInputs<dim> &in,
72  std::vector<double> &melt_fractions) const override;
73 
78  double reference_darcy_coefficient () const override;
79 
80 
92  static
93  void
94  declare_parameters (ParameterHandler &prm);
95 
99  void
100  parse_parameters (ParameterHandler &prm) override;
105  void
106  create_additional_named_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const override;
107 
108 
109  private:
112  double reference_T;
113  double eta_0;
114  double xi_0;
115  double eta_f;
122  double alpha_phi;
133 
134  // entropy change upon melting
136 
137  virtual
138  double
139  melt_fraction (const double temperature,
140  const double pressure,
141  const double depletion) const;
142  };
143 
144  }
145 }
146 
147 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:88