kEpsilonViscosityAux.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 "AuxKernel.h"
#include "INSFVVelocityVariable.h"
#include "NS.h"

class kEpsilonViscosityAux : public AuxKernel
{
public:
  static InputParameters validParams();

  virtual void initialSetup() override;

  kEpsilonViscosityAux(const InputParameters & parameters);

protected:
  virtual Real computeValue() override;

  const unsigned int _dim;

  const Moose::Functor<ADReal> & _u_var;
  const Moose::Functor<ADReal> * _v_var;
  const Moose::Functor<ADReal> * _w_var;

  const Moose::Functor<ADReal> & _k;
  const Moose::Functor<ADReal> & _epsilon;

  const Moose::Functor<ADReal> & _rho;
  const Moose::Functor<ADReal> & _mu;

  const Real _C_mu;

  // Maximum allowable mu_t_ratio : mu/mu_t
  const Real _mu_t_ratio_max;

  const std::vector<BoundaryName> & _wall_boundary_names;

  const bool _bulk_wall_treatment;

  NS::WallTreatmentEnum _wall_treatment;

  const MooseEnum _scale_limiter;

  const bool _newton_solve;

  // -- Parameters of the wall function method

  static constexpr int _MAX_ITERS_U_TAU{50};

  static constexpr Real _REL_TOLERANCE{1e-4};

  std::map<const Elem *, bool> _wall_bounded;
  std::map<const Elem *, std::vector<Real>> _dist;
  std::map<const Elem *, std::vector<const FaceInfo *>> _face_infos;
};