PorousFlowAqueousPreDisChemistry.C File Reference

Go to the source code of this file.

Functions
	registerMooseObject ("PorousFlowApp", PorousFlowAqueousPreDisChemistry)
template<>
InputParameters	validParams< PorousFlowAqueousPreDisChemistry > ()

Function Documentation

registerMooseObject()

registerMooseObject	(	"PorousFlowApp"	,
		PorousFlowAqueousPreDisChemistry
	)

validParams< PorousFlowAqueousPreDisChemistry >()

template<>

InputParameters validParams< PorousFlowAqueousPreDisChemistry >

(

)

Definition at line 16 of file PorousFlowAqueousPreDisChemistry.C.

{
  InputParameters params = validParams<PorousFlowMaterialVectorBase>();
  params.addRequiredCoupledVar(
      "primary_concentrations",
      "List of MOOSE Variables that represent the concentrations of the primary species");
  params.addRequiredParam<unsigned>("num_reactions",
                                    "Number of equations in the system of chemical reactions");
  params.addParam<bool>("equilibrium_constants_as_log10",
                        false,
                        "If true, the equilibrium constants are written in their log10 form, eg, "
                        "-2.  If false, the equilibrium constants are written in absolute terms, "
                        "eg, 0.01");
  params.addRequiredCoupledVar("equilibrium_constants",
                               "Equilibrium constant for each equation (dimensionless).  If these "
                               "are temperature dependent AuxVariables, the Jacobian will not be "
                               "exact");
  params.addRequiredParam<std::vector<Real>>(
      "primary_activity_coefficients",
      "Activity coefficients for the primary species (dimensionless) (one for each)");
  params.addRequiredParam<std::vector<Real>>(
      "reactions",
      "A matrix defining the aqueous reactions.  The matrix is entered as a long vector: the first "
      "row is "
      "entered first, followed by the second row, etc.  There should be num_reactions rows.  All "
      "primary species should appear only on the LHS of each reaction (and there should be just "
      "one secondary species on the RHS, by definition) so they may have negative coefficients.  "
      "Each row should have number of primary_concentrations entries, which are the stoichiometric "
      "coefficients.  The first coefficient must always correspond to the first primary species, "
      "etc");
  params.addRequiredParam<std::vector<Real>>("specific_reactive_surface_area",
                                             "Specific reactive surface area in m^2/(L solution).");
  params.addRequiredParam<std::vector<Real>>(
      "kinetic_rate_constant",
      "Kinetic rate constant in mol/(m^2 s), at the reference temperature (one for each reaction)");
  params.addRequiredParam<std::vector<Real>>("molar_volume",
                                             "Volume occupied by one mole of the secondary species "
                                             "(L(solution)/mol) (one for each reaction)");
  params.addRequiredParam<std::vector<Real>>("activation_energy",
                                             "Activation energy, J/mol (one for each reaction)");
  params.addParam<Real>("gas_constant", 8.31434, "Gas constant, in J/(mol K)");
  params.addParam<Real>("reference_temperature", 298.15, "Reference temperature, K");
  params.addParam<std::vector<Real>>("theta_exponent",
                                     "Theta exponent.  Defaults to 1.  (one for each reaction)");
  params.addParam<std::vector<Real>>("eta_exponent",
                                     "Eta exponent.  Defaults to 1.  (one for each reaction)");
  params.addPrivateParam<std::string>("pf_material_type", "chemistry");
  params.addClassDescription("This Material forms a std::vector of mineralisation reaction rates "
                             "(L(precipitate)/L(solution)/s) appropriate to the aqueous "
                             "precipitation-dissolution system provided.  Note: the "
                             "PorousFlowTemperature must be measured in Kelvin.");
  return params;
}