21 #ifndef _aspect_material_model_reactive_fluid_transport_h 22 #define _aspect_material_model_reactive_fluid_transport_h 39 namespace MaterialModel
58 bool is_compressible ()
const override;
64 virtual double reference_darcy_coefficient ()
const override;
76 unsigned int q)
const;
88 std::vector<double> &melt_fractions)
const override;
100 void update()
override;
119 parse_parameters (ParameterHandler &prm)
override;
170 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};
171 std::vector<double> csat_peridotite_poly_coeffs {0.00115628, 2.42179};
172 std::vector<double> Td_peridotite_poly_coeffs {-15.4627, 94.9716, 636.603};
174 std::vector<double> LR_gabbro_poly_coeffs {-1.81745, 7.67198, -10.8507, 5.09329, 8.14519};
175 std::vector<double> csat_gabbro_poly_coeffs {-0.0176673, 0.0893044, 1.52732};
176 std::vector<double> Td_gabbro_poly_coeffs {-1.72277, 20.5898, 637.517};
178 std::vector<double> LR_MORB_poly_coeffs {-1.78177, 7.50871, -10.4840, 5.19725, 7.96365};
179 std::vector<double> csat_MORB_poly_coeffs {0.0102725, -0.115390, 0.324452, 1.41588};
180 std::vector<double> Td_MORB_poly_coeffs {-3.81280, 22.7809, 638.049};
182 std::vector<double> LR_sediment_poly_coeffs {-2.03283, 10.8186, -21.2119, 18.3351, -6.48711, 8.32459};
183 std::vector<double> csat_sediment_poly_coeffs {-0.150662, 0.301807, 1.01867};
184 std::vector<double> Td_sediment_poly_coeffs {2.83277, -24.7593, 85.9090, 524.898};
189 const std::vector<double> pressure_cutoffs {10, 26, 16, 50};
191 std::vector<std::vector<double>> devolatilization_enthalpy_changes {LR_peridotite_poly_coeffs, LR_gabbro_poly_coeffs, \
192 LR_MORB_poly_coeffs, LR_sediment_poly_coeffs
195 std::vector<std::vector<double>> water_mass_fractions {csat_peridotite_poly_coeffs, csat_gabbro_poly_coeffs, \
196 csat_MORB_poly_coeffs, csat_sediment_poly_coeffs
199 std::vector<std::vector<double>> devolatilization_onset_temperatures {Td_peridotite_poly_coeffs, Td_gabbro_poly_coeffs, \
200 Td_MORB_poly_coeffs, Td_sediment_poly_coeffs
257 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 min_compaction_visc
double fluid_reaction_time_scale
double tian_max_peridotite_water
double max_compaction_visc
double fluid_compressibility
double tian_max_gabbro_water
double reference_permeability
ReactionModel::Katz2003MantleMelting< dim > katz2003_model
double tian_max_MORB_water