LinearFVAdvection.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 "LinearFVFluxKernel.h"
#include "FVAdvectedInterpolationMethod.h"
#include "FVInterpolationMethodInterface.h"

class LinearFVAdvection : public LinearFVFluxKernel, public FVInterpolationMethodInterface
{
public:
  static InputParameters validParams();

  LinearFVAdvection(const InputParameters & params);

  virtual void setupFaceData(const FaceInfo * face_info) override;
  virtual void initialSetup() override;

  virtual Real computeElemMatrixContribution() override;

  virtual Real computeNeighborMatrixContribution() override;

  virtual Real computeElemRightHandSideContribution() override;

  virtual Real computeNeighborRightHandSideContribution() override;

  virtual Real computeBoundaryMatrixContribution(const LinearFVBoundaryCondition & bc) override;

  virtual Real computeBoundaryRHSContribution(const LinearFVBoundaryCondition & bc) override;

protected:
  const RealVectorValue _velocity;

  const FVAdvectedInterpolationMethod & _adv_interp_method;

  FVAdvectedInterpolationMethod::AdvectedSystemContribution _adv_interp_result;

  VectorValue<Real> _elem_grad_storage;
  VectorValue<Real> _neighbor_grad_storage;

  // Cache for the advected face flux
  Real _adv_face_flux = 0.0;
};