FunctorPebbleBedDragCoefficients.h

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

#pragma once

#include "FunctorIsotropicDragCoefficients.h"
#include "NS.h"

template <typename Derived>
class FunctorPebbleBedDragCoefficients : public FunctorIsotropicDragCoefficients<Derived>
{
  friend class FunctorIsotropicDragCoefficients<Derived>;

public:
  FunctorPebbleBedDragCoefficients(const InputParameters & parameters);

  static InputParameters validParams();

protected:
  template <typename Space, typename Time>
  ADReal computeDarcyPrefactor(const Space & r, const Time & t);

  template <typename Space, typename Time>
  ADReal computeForchheimerPrefactor(const Space & r, const Time & t);

  template <typename Space, typename Time>
  ADReal computeHydraulicDiameter(const Space & r, const Time & t);

  const Moose::Functor<ADReal> & _eps;

  const Moose::Functor<ADReal> & _rho;

  const Real & _d_pebble;
};

template <typename Derived>
InputParameters
FunctorPebbleBedDragCoefficients<Derived>::validParams()
{
  InputParameters params = FunctorIsotropicDragCoefficients<Derived>::validParams();
  params.addRequiredRangeCheckedParam<Real>(
      NS::pebble_diameter, NS::pebble_diameter + " > 0.0", "Pebble diameter");
  params.addRequiredParam<MooseFunctorName>(NS::porosity, "Porosity");

  params.addParam<MooseFunctorName>(NS::density, NS::density, "Density");
  return params;
}

template <typename Derived>
FunctorPebbleBedDragCoefficients<Derived>::FunctorPebbleBedDragCoefficients(
    const InputParameters & parameters)
  : FunctorIsotropicDragCoefficients<Derived>(parameters),
    _eps(this->template getFunctor<ADReal>(NS::porosity)),
    _rho(this->template getFunctor<ADReal>(NS::density)),
    _d_pebble(this->template getParam<Real>(NS::pebble_diameter))
{
}

template <typename Derived>
template <typename Space, typename Time>
ADReal
FunctorPebbleBedDragCoefficients<Derived>::computeDarcyPrefactor(const Space & r, const Time & t)
{
  const auto Dh = computeHydraulicDiameter(r, t);
  return 1 / (Dh * Dh);
}

template <typename Derived>
template <typename Space, typename Time>
ADReal
FunctorPebbleBedDragCoefficients<Derived>::computeForchheimerPrefactor(const Space & r,
                                                                       const Time & t)
{
  // Why the factor of 2? Because we multiply by a factor of 1/2 in the friction kernel to be
  // consistent with the formulation of Forchheimer in
  // https://holzmann-cfd.com/community/blog-and-tools/darcy-forchheimer and
  // https://www.simscale.com/knowledge-base/predict-darcy-and-forchheimer-coefficients-for-perforated-plates-using-analytical-approach/
  return 2 / computeHydraulicDiameter(r, t);
}

template <typename Derived>
template <typename Space, typename Time>
ADReal
FunctorPebbleBedDragCoefficients<Derived>::computeHydraulicDiameter(const Space & r, const Time & t)
{
  mooseAssert(1.0 - _eps(r, t).value() > 1e-8,
              "Bed hydraulic diameter is ill-defined at porosity of 1.");
  return _eps(r, t) * _d_pebble / (1.0 - _eps(r, t));
}