AqueousReactionsEquilibriumPhysics.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ReactionNetworkPhysicsBase.h"

class AqueousReactionsEquilibriumPhysics : public ReactionNetworkPhysicsBase
{
public:
  static InputParameters validParams();

  AqueousReactionsEquilibriumPhysics(const InputParameters & parameters);

protected:
  virtual void addAuxiliaryKernels() override;
  virtual void addFEKernels() override;

  std::vector<std::vector<Real>> _stos;
  std::vector<std::vector<Real>> _sto_u;
  std::vector<std::vector<std::vector<Real>>> _sto_v;
  std::vector<std::vector<Real>> _weights;
  std::vector<Real> _log_eq_const;
  std::vector<VariableName> _eq_species;
  std::vector<std::vector<bool>> _primary_participation;
  std::vector<std::vector<std::vector<VariableName>>> _coupled_v;
  std::vector<std::vector<VariableName>> _solver_species_involved;

  const std::vector<VariableName> & _pressure_var;
  const RealVectorValue _gravity;
};