ASPECT
melt_global.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2015 - 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_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  {
44  template <int dim>
47  public ::aspect::SimulatorAccess<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 
76  std::vector<double> &melt_fractions,
77  const MaterialModel::MaterialModelOutputs<dim> *out = nullptr) const override;
78 
83  double reference_darcy_coefficient () const override;
84 
85 
97  static
98  void
99  declare_parameters (ParameterHandler &prm);
100 
104  void
105  parse_parameters (ParameterHandler &prm) override;
110  void
112 
113 
114  private:
117  double reference_T;
118  double eta_0;
119  double xi_0;
120  double eta_f;
127  double alpha_phi;
138 
139  // entropy change upon melting
141 
142  virtual
143  double
144  melt_fraction (const double temperature,
145  const double pressure,
146  const double depletion) const;
147  };
148 
149  }
150 }
151 
152 #endif
void parse_parameters(ParameterHandler &prm) override
void create_additional_named_outputs(MaterialModel::MaterialModelOutputs< dim > &out) const override
double reference_darcy_coefficient() const override
bool is_compressible() const override
void melt_fractions(const MaterialModel::MaterialModelInputs< dim > &in, std::vector< double > &melt_fractions, const MaterialModel::MaterialModelOutputs< dim > *out=nullptr) const override
virtual double melt_fraction(const double temperature, const double pressure, const double depletion) const
static void declare_parameters(ParameterHandler &prm)
void evaluate(const typename Interface< dim >::MaterialModelInputs &in, typename Interface< dim >::MaterialModelOutputs &out) const override