INSFVTKEDSourceSink.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 "FVElementalKernel.h"
#include "MathFVUtils.h"
#include "INSFVMomentumResidualObject.h"
#include "INSFVVelocityVariable.h"
#include "NS.h"

class INSFVTKEDSourceSink : public FVElementalKernel
{
public:
  static InputParameters validParams();

  virtual void initialSetup() override;

  INSFVTKEDSourceSink(const InputParameters & parameters);

protected:
  ADReal computeQpResidual() override;

protected:
  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> & _rho;

  const Moose::Functor<ADReal> & _mu;

  const Moose::Functor<ADReal> & _mu_t;

  const std::vector<BoundaryName> & _wall_boundary_names;

  const bool _linearized_model;

  NS::WallTreatmentEnum _wall_treatment;

  const Real _C1_eps;

  const Real _C2_eps;

  const Real _C_mu;

  // Production Limiter Constant
  const Real _C_pl;

  std::map<const Elem *, Real> _symmetric_strain_tensor_norm_old;
  std::map<const Elem *, Real> _old_destruction;

  std::map<const Elem *, Real> _pevious_nl_sol;


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

  const bool _newton_solve;
};