InternalSideIntegralVariablePostprocessor.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 "InternalSideIntegralVariablePostprocessor.h"
#include "MathFVUtils.h"

#include "metaphysicl/raw_type.h"

registerMooseObject("MooseApp", InternalSideIntegralVariablePostprocessor);

InputParameters
InternalSideIntegralVariablePostprocessor::validParams()
{
  InputParameters params = InternalSideIntegralPostprocessor::validParams();
  params.addRequiredCoupledVar("variable",
                               "The name of the variable which this postprocessor integrates");
  params.addClassDescription("Computes an integral on internal sides of the specified variable");
  return params;
}

InternalSideIntegralVariablePostprocessor::InternalSideIntegralVariablePostprocessor(
    const InputParameters & parameters)
  : InternalSideIntegralPostprocessor(parameters),
    MooseVariableInterface<Real>(this,
                                 false,
                                 "variable",
                                 Moose::VarKindType::VAR_ANY,
                                 Moose::VarFieldType::VAR_FIELD_STANDARD),
    _u(coupledValue("variable")),
    _grad_u(coupledGradient("variable")),
    _fv(_fv_variable)
{
  addMooseVariableDependency(&mooseVariableField());

  _qp_integration = !_fv;
}

Real
InternalSideIntegralVariablePostprocessor::computeFaceInfoIntegral(const FaceInfo * const fi)
{
  return MetaPhysicL::raw_value((*_fv_variable)(makeCDFace(*fi), determineState()));
}

Real
InternalSideIntegralVariablePostprocessor::computeQpIntegral()
{
  return _u[_qp];
}

bool
InternalSideIntegralVariablePostprocessor::hasFaceSide(const FaceInfo & fi,
                                                       const bool fi_elem_side) const
{
  mooseAssert(_fv_variable, "Must be non-null");
  return _fv_variable->hasFaceSide(fi, fi_elem_side);
}