PorousFlowDiffusivityBase.C

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
 
#include "PorousFlowDiffusivityBase.h"
 
template <>
InputParameters
validParams<PorousFlowDiffusivityBase>()
{
  InputParameters params = validParams<PorousFlowMaterialVectorBase>();
  params.addRequiredParam<std::vector<Real>>(
      "diffusion_coeff",
      "List of diffusion coefficients.  Order is i) component 0 in phase 0; ii) "
      "component 1 in phase 0 ...; component 0 in phase 1; ... component k in "
      "phase n (m^2/s");
  params.addClassDescription("Base class for effective diffusivity for each phase");
  params.set<bool>("at_nodes") = false;
  params.addPrivateParam<std::string>("pf_material_type", "diffusivity");
  return params;
}
 
PorousFlowDiffusivityBase::PorousFlowDiffusivityBase(const InputParameters & parameters)
  : PorousFlowMaterialVectorBase(parameters),
 
    _tortuosity(declareProperty<std::vector<Real>>("PorousFlow_tortuosity_qp")),
    _dtortuosity_dvar(
        declareProperty<std::vector<std::vector<Real>>>("dPorousFlow_tortuosity_qp_dvar")),
    _diffusion_coeff(
        declareProperty<std::vector<std::vector<Real>>>("PorousFlow_diffusion_coeff_qp")),
    _ddiffusion_coeff_dvar(declareProperty<std::vector<std::vector<std::vector<Real>>>>(
        "dPorousFlow_diffusion_coeff_qp_dvar")),
    _input_diffusion_coeff(getParam<std::vector<Real>>("diffusion_coeff"))
{
  // Also, the number of diffusion coefficients must be equal to the num_phases * num_components
  if (_input_diffusion_coeff.size() != _num_phases * _num_components)
    paramError("diffusion_coeff",
               "The number of diffusion coefficients entered is not equal to the number of phases "
               "multiplied by the number of fluid components");
 
  if (_nodal_material == true)
    mooseError("PorousFlowRelativeDiffusivity classes are only defined for at_nodes = false");
}
 
void
PorousFlowDiffusivityBase::computeQpProperties()
{
  _diffusion_coeff[_qp].resize(_num_phases);
  _ddiffusion_coeff_dvar[_qp].resize(_num_phases);
  _dtortuosity_dvar[_qp].resize(_num_phases);
 
  for (unsigned int ph = 0; ph < _num_phases; ++ph)
  {
    _diffusion_coeff[_qp][ph].resize(_num_components);
    _ddiffusion_coeff_dvar[_qp][ph].resize(_num_components);
    _dtortuosity_dvar[_qp][ph].assign(_num_var, 0.0);
 
    for (unsigned int comp = 0; comp < _num_components; ++comp)
    {
      _diffusion_coeff[_qp][ph][comp] = _input_diffusion_coeff[ph + comp];
      _ddiffusion_coeff_dvar[_qp][ph][comp].assign(_num_var, 0.0);
    }
  }
}