PorousFlowActionBase.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "Action.h"
#include "PorousFlowDependencies.h"
#include "MooseEnum.h"
 
#include "libmesh/vector_value.h"
 
class PorousFlowActionBase;
 
template <>
InputParameters validParams<PorousFlowActionBase>();
 
class PorousFlowActionBase : public Action, public PorousFlowDependencies
{
public:
  PorousFlowActionBase(const InputParameters & params);
 
  virtual void act() override;
 
  using Action::addRelationshipManagers;
  virtual void addRelationshipManagers(Moose::RelationshipManagerType when_type) override;
 
protected:
  std::vector<std::string> _included_objects;
 
  const std::string _dictator_name;
 
  const unsigned int _num_aqueous_equilibrium;
 
  const unsigned int _num_aqueous_kinetic;
 
  const RealVectorValue _gravity;
 
  const std::vector<VariableName> & _mass_fraction_vars;
 
  const unsigned _num_mass_fraction_vars;
 
  const std::vector<VariableName> & _temperature_var;
 
  const std::vector<VariableName> & _displacements;
 
  const unsigned _ndisp;
 
  std::vector<VariableName> _coupled_displacements;
 
  const MooseEnum _flux_limiter_type;
 
  const enum class StabilizationEnum { Full, KT } _stabilization;
 
  Moose::CoordinateSystemType _coord_system;
 
  bool _transient;
 
  virtual void addDictator() = 0;
 
  virtual void addUserObjects();
 
  virtual void addAuxObjects();
 
  virtual void addKernels();
 
  virtual void addMaterialDependencies();
 
  virtual void addMaterials();
 
  void addSaturationAux(unsigned phase);
 
  void addDarcyAux(const RealVectorValue & gravity);
 
  void addStressAux();
 
  void addTemperatureMaterial(bool at_nodes);
 
  void addMassFractionMaterial(bool at_nodes);
 
  void addEffectiveFluidPressureMaterial(bool at_nodes);
 
  void addVolumetricStrainMaterial(const std::vector<VariableName> & displacements,
                                   bool consistent_with_displaced_mesh);
 
  void addSingleComponentFluidMaterial(bool at_nodes,
                                       unsigned phase,
                                       bool compute_density_and_viscosity,
                                       bool compute_internal_energy,
                                       bool compute_enthalpy,
                                       const UserObjectName & fp);
 
  void addBrineMaterial(const VariableName xnacl,
                        bool at_nodes,
                        unsigned phase,
                        bool compute_density_and_viscosity,
                        bool compute_internal_energy,
                        bool compute_enthalpy);
 
  void
  addRelativePermeabilityCorey(bool at_nodes, unsigned phase, Real n, Real s_res, Real sum_s_res);
 
  void
  addRelativePermeabilityFLAC(bool at_nodes, unsigned phase, Real m, Real s_res, Real sum_s_res);
 
  void addCapillaryPressureVG(Real m, Real alpha, std::string userobject_name);
 
  void addAdvectiveFluxCalculatorSaturated(unsigned phase,
                                           bool multiply_by_density,
                                           std::string userobject_name);
 
  void addAdvectiveFluxCalculatorUnsaturated(unsigned phase,
                                             bool multiply_by_density,
                                             std::string userobject_name);
 
  void addAdvectiveFluxCalculatorSaturatedHeat(unsigned phase,
                                               bool multiply_by_density,
                                               std::string userobject_name);
 
  void addAdvectiveFluxCalculatorUnsaturatedHeat(unsigned phase,
                                                 bool multiply_by_density,
                                                 std::string userobject_name);
 
  void addAdvectiveFluxCalculatorSaturatedMultiComponent(unsigned phase,
                                                         unsigned fluid_component,
                                                         bool multiply_by_density,
                                                         std::string userobject_name);
 
  void addAdvectiveFluxCalculatorUnsaturatedMultiComponent(unsigned phase,
                                                           unsigned fluid_component,
                                                           bool multiply_by_density,
                                                           std::string userobject_name);
};