LinearFVFluxKernel.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 "LinearFVKernel.h"
#include "FaceArgInterface.h"

class LinearFVFluxKernel : public LinearFVKernel, public FaceArgProducerInterface
{
public:
  static InputParameters validParams();

  LinearFVFluxKernel(const InputParameters & params);

  virtual void addMatrixContribution() override;

  virtual void addRightHandSideContribution() override;

  virtual bool hasFaceSide(const FaceInfo & fi, bool fi_elem_side) const override;

  virtual void setupFaceData(const FaceInfo * face_info);

  void setCurrentFaceArea(const Real area) { _current_face_area = area; };

  virtual Real computeElemMatrixContribution() = 0;

  virtual Real computeNeighborMatrixContribution() = 0;

  virtual Real computeElemRightHandSideContribution() = 0;

  virtual Real computeNeighborRightHandSideContribution() = 0;

  virtual Real computeBoundaryMatrixContribution(const LinearFVBoundaryCondition & bc) = 0;

  virtual Real computeBoundaryRHSContribution(const LinearFVBoundaryCondition & bc) = 0;

protected:
  Moose::FaceArg singleSidedFaceArg(
      const FaceInfo * fi,
      Moose::FV::LimiterType limiter_type = Moose::FV::LimiterType::CentralDifference,
      bool correct_skewness = false) const;

  const FaceInfo * _current_face_info;

  Real _current_face_area;

  FaceInfo::VarFaceNeighbors _current_face_type;

  bool _cached_matrix_contribution;

  bool _cached_rhs_contribution;

  const bool _force_boundary_execution;

  DenseVector<dof_id_type> _dof_indices;

  DenseMatrix<Real> _matrix_contribution;

  DenseVector<Real> _rhs_contribution;
};