IncompressibleMomentumSPBase.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 "ScalarKernel.h"
#include "ADScalarKernel.h"
#include "FunctorInterface.h"
 
class SinglePhaseFluidProperties;
 
template <bool is_ad>
class IncompressibleMomentumSPBaseTempl
  : public std::conditional<is_ad, ADScalarKernel, ScalarKernel>::type,
    public FunctorInterface
{
  using Base = typename std::conditional<is_ad, ADScalarKernel, ScalarKernel>::type;
 
public:
  IncompressibleMomentumSPBaseTempl(const InputParameters & parameters);
  virtual bool isADObject() const override { return is_ad; };
  static InputParameters validParams();
 
protected:
  virtual GenericReal<is_ad> computeQpResidual() override;
  virtual Real computeQpJacobian() override;
  virtual Real computeQpJacobianMDot();
  virtual Real computeQpJacobianDP();
  virtual GenericReal<is_ad> computeFrictionFactor(const GenericReal<is_ad> & mu,
                                                   const GenericReal<is_ad> & G,
                                                   const GenericReal<is_ad> & Dh,
                                                   const unsigned int j);
  virtual GenericReal<is_ad> massFlowRate() const = 0;
  virtual GenericReal<is_ad> pressureDrop() const = 0;
  const size_t _n_temps;
  std::vector<const VariableValue *> _T;
  const bool _is_implicit;
  const Moose::Functor<GenericReal<is_ad>> & _Pref;
  const SinglePhaseFluidProperties & _fp;
  const size_t _n_segments;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _areas;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _perimeters;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _lengths;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _alphas;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _forms_losses;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _dPps;
  std::vector<const Moose::Functor<GenericReal<is_ad>> *> _roughnesses;
  const Moose::Functor<GenericReal<is_ad>> & _gravity;
};