21 #ifndef _aspect_material_model_reactive_fluid_transport_h 22 #define _aspect_material_model_reactive_fluid_transport_h 38 namespace MaterialModel
57 bool is_compressible ()
const override;
63 virtual double reference_darcy_coefficient ()
const override;
75 unsigned int q)
const;
87 std::vector<double> &melt_fractions)
const override;
99 void update()
override;
118 parse_parameters (ParameterHandler &prm)
override;
166 std::vector<double> LR_peridotite_poly_coeffs {-19.0609, 168.983, -630.032, 1281.84, -1543.14, 1111.88, -459.142, 95.4143, 1.97246};
167 std::vector<double> csat_peridotite_poly_coeffs {0.00115628, 2.42179};
168 std::vector<double> Td_peridotite_poly_coeffs {-15.4627, 94.9716, 636.603};
170 std::vector<double> LR_gabbro_poly_coeffs {-1.81745, 7.67198, -10.8507, 5.09329, 8.14519};
171 std::vector<double> csat_gabbro_poly_coeffs {-0.0176673, 0.0893044, 1.52732};
172 std::vector<double> Td_gabbro_poly_coeffs {-1.72277, 20.5898, 637.517};
174 std::vector<double> LR_MORB_poly_coeffs {-1.78177, 7.50871, -10.4840, 5.19725, 7.96365};
175 std::vector<double> csat_MORB_poly_coeffs {0.0102725, -0.115390, 0.324452, 1.41588};
176 std::vector<double> Td_MORB_poly_coeffs {-3.81280, 22.7809, 638.049};
178 std::vector<double> LR_sediment_poly_coeffs {-2.03283, 10.8186, -21.2119, 18.3351, -6.48711, 8.32459};
179 std::vector<double> csat_sediment_poly_coeffs {-0.150662, 0.301807, 1.01867};
180 std::vector<double> Td_sediment_poly_coeffs {2.83277, -24.7593, 85.9090, 524.898};
182 std::vector<std::vector<double>> devolatilization_enthalpy_changes {LR_peridotite_poly_coeffs, LR_gabbro_poly_coeffs, \
183 LR_MORB_poly_coeffs, LR_sediment_poly_coeffs
186 std::vector<std::vector<double>> water_mass_fractions {csat_peridotite_poly_coeffs, csat_gabbro_poly_coeffs, \
187 csat_MORB_poly_coeffs, csat_sediment_poly_coeffs
190 std::vector<std::vector<double>> devolatilization_onset_temperatures {Td_peridotite_poly_coeffs, Td_gabbro_poly_coeffs, \
191 Td_MORB_poly_coeffs, Td_sediment_poly_coeffs
238 fluid_solid_reaction_scheme;
double shear_to_bulk_viscosity_ratio
std::unique_ptr< MaterialModel::Interface< dim > > base_model
void declare_parameters(ParameterHandler &prm)
double tian_max_sediment_water
double fluid_reaction_time_scale
double tian_max_peridotite_water
double fluid_compressibility
double tian_max_gabbro_water
double reference_permeability
double tian_max_MORB_water