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<ADReal>(NS::speed)),
    _rho(getFunctor<ADReal>(NS::density)),
    _mu(getFunctor<ADReal>(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() * raw_value(_rho(qp_arg, state)) *
           raw_value(_speed(qp_arg, state)) / raw_value(_mu(qp_arg, state));
  }
  else
  {
    const auto elem_arg = makeElemArg(_current_elem);
    return _current_elem->hmax() * raw_value(_rho(elem_arg, state)) *
           raw_value(_speed(elem_arg, state)) / raw_value(_mu(elem_arg, state));
  }
}