PorousFlowPropertyAux.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 "AuxKernel.h"
#include "PorousFlowDictator.h"
#include "PorousFlowConstants.h"
 
template <bool is_ad>
class PorousFlowPropertyAuxTempl : public AuxKernel
{
public:
  static InputParameters validParams();
 
  PorousFlowPropertyAuxTempl(const InputParameters & parameters);
 
protected:
  virtual Real computeValue() override;
 
private:
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _pressure;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _saturation;
 
  const GenericMaterialProperty<Real, is_ad> * _temperature;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _fluid_density;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _fluid_viscosity;
 
  const GenericMaterialProperty<std::vector<std::vector<Real>>, is_ad> * _mass_fractions;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _relative_permeability;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _enthalpy;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _internal_energy;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _sec_conc;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _mineral_conc;
 
  const GenericMaterialProperty<std::vector<Real>, is_ad> * _mineral_reaction_rate;
 
  const GenericMaterialProperty<Real, is_ad> * _porosity;
 
  const GenericMaterialProperty<RealTensorValue, is_ad> * _permeability;
 
  const MaterialProperty<unsigned int> * _hys_order;
 
  const MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>> *
      _hys_sat_tps;
 
  const MaterialProperty<Real> * _hys_info;
 
  const PorousFlowDictator & _dictator;
 
  const enum class PropertyEnum {
    PRESSURE,
    SATURATION,
    TEMPERATURE,
    DENSITY,
    VISCOSITY,
    MASS_FRACTION,
    RELPERM,
    CAPILLARY_PRESSURE,
    ENTHALPY,
    INTERNAL_ENERGY,
    SECONDARY_CONCENTRATION,
    MINERAL_CONCENTRATION,
    MINERAL_REACTION_RATE,
    POROSITY,
    PERMEABILITY,
    HYSTERESIS_ORDER,
    HYSTERESIS_SATURATION_TURNING_POINT,
    HYSTERETIC_INFO
  } _property_enum;
 
  const unsigned int _phase;
 
  const unsigned int _liquid_phase;
 
  const unsigned int _gas_phase;
 
  const unsigned int _fluid_component;
 
  const unsigned int _secondary_species;
 
  const unsigned int _mineral_species;
 
  const unsigned int _hysteresis_turning_point;
 
  const unsigned int _k_row;
  const unsigned int _k_col;
};
 
typedef PorousFlowPropertyAuxTempl<false> PorousFlowPropertyAux;
typedef PorousFlowPropertyAuxTempl<true> ADPorousFlowPropertyAux;