AEFVMaterial.C

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
 
#include "AEFVMaterial.h"
#include "MooseMesh.h"
 
#include "libmesh/quadrature.h"
 
using namespace libMesh;
 
registerMooseObject("RdgApp", AEFVMaterial);
 
InputParameters
AEFVMaterial::validParams()
{
  InputParameters params = Material::validParams();
  params.addClassDescription(
      "A material kernel for the advection equation using a cell-centered finite volume method.");
  params.addRequiredCoupledVar("u", "Cell-averge variable");
  params.addRequiredParam<UserObjectName>("slope_limiting", "Name for slope limiting user object");
  return params;
}
 
AEFVMaterial::AEFVMaterial(const InputParameters & parameters)
  : Material(parameters),
    _uc(coupledValue("u")),
    _lslope(getUserObject<SlopeLimitingBase>("slope_limiting")),
    _u(declareProperty<Real>("u"))
{
}
 
AEFVMaterial::~AEFVMaterial() {}
 
void
AEFVMaterial::computeQpProperties()
{
  // initialize the variable
  _u[_qp] = _uc[_qp];
 
  // interpolate variable values at face center
  if (_bnd)
  {
    // you should know how many equations you are solving and assign this number
    // e.g. = 1 (for the advection equation)
    unsigned int nvars = 1;
    std::vector<RealGradient> ugrad(nvars, RealGradient(0., 0., 0.));
    ugrad = _lslope.getElementSlope(_current_elem->id());
 
    // get the directional vector from cell center to face center
    RealGradient dvec = _q_point[_qp] - _current_elem->vertex_average();
 
    // calculate the variable at face center
    _u[_qp] += ugrad[0] * dvec;
 
    // clear the temporary vectors
    ugrad.clear();
  }
  // calculations only for elemental output
  else if (!_bnd)
  {
  }
}