ReynoldsNumberFunctorAux.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 "ReynoldsNumberFunctorAux.h"
#include "MooseMesh.h"
#include "NS.h"
 
registerMooseObject("NavierStokesApp", ReynoldsNumberFunctorAux);
 
InputParameters
ReynoldsNumberFunctorAux::validParams()
{
  InputParameters params = AuxKernel::validParams();
  params.addClassDescription("Computes rho*u*L/mu.");
  params.addRequiredParam<MooseFunctorName>(NS::speed, "The fluid speed");
  params.addRequiredParam<MooseFunctorName>(NS::density, "The fluid density");
  params.addRequiredParam<MooseFunctorName>(NS::mu, "The fluid dynamic viscosity");
  return params;
}
 
ReynoldsNumberFunctorAux::ReynoldsNumberFunctorAux(const InputParameters & parameters)
  : AuxKernel(parameters),
    _speed(getFunctor<Real>(NS::speed)),
    _rho(getFunctor<Real>(NS::density)),
    _mu(getFunctor<Real>(NS::mu)),
    _use_qp_arg(dynamic_cast<MooseVariableFE<Real> *>(&_var))
{
  if (!_use_qp_arg && !dynamic_cast<MooseVariableFV<Real> *>(&_var))
    paramError(
        "variable",
        "The variable must be a non-vector, non-array finite-volume/finite-element variable.");
 
  if (isNodal())
    mooseError("This AuxKernel only supports Elemental fields");
}
 
Real
ReynoldsNumberFunctorAux::computeValue()
{
  using MetaPhysicL::raw_value;
  const auto state = determineState();
 
  if (_use_qp_arg)
  {
    const Moose::ElemQpArg qp_arg = {_current_elem, _qp, _qrule, _q_point[_qp]};
    return _current_elem->hmax() * _rho(qp_arg, state) * _speed(qp_arg, state) / _mu(qp_arg, state);
  }
  else
  {
    const auto elem_arg = makeElemArg(_current_elem);
    return _current_elem->hmax() * _rho(elem_arg, state) * _speed(elem_arg, state) /
           _mu(elem_arg, state);
  }
}