LinearWCNSFVMomentumFlux.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"

#include <array>

class RhieChowMassFlux;
class LinearFVBoundaryCondition;
class LinearFVAdvectionDiffusionBC;

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

  LinearWCNSFVMomentumFlux(const InputParameters & params);

  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;

  virtual void setupFaceData(const FaceInfo * face_info) override;

protected:
  Real computeInternalAdvectionElemMatrixContribution();

  Real computeInternalAdvectionNeighborMatrixContribution();

  Real computeInternalStressMatrixContribution();

  Real computeInternalStressRHSContribution();

  Real computeStressBoundaryMatrixContribution(const LinearFVAdvectionDiffusionBC * bc);

  Real computeStressBoundaryRHSContribution(const LinearFVAdvectionDiffusionBC * bc);

  Real computeAdvectionBoundaryMatrixContribution(const LinearFVAdvectionDiffusionBC * bc);

  Real computeAdvectionBoundaryRHSContribution(const LinearFVAdvectionDiffusionBC * bc);

  const unsigned int _dim;

  const RhieChowMassFlux & _mass_flux_provider;

  const Moose::Functor<Real> & _mu;

  const bool _use_nonorthogonal_correction;

  const bool _use_deviatoric_terms;

  const FVAdvectedInterpolationMethod & _adv_interp_method;

  FVAdvectedInterpolationMethod::AdvectedSystemContribution _adv_interp_result;

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

  Real _face_mass_flux;

  Real _boundary_normal_factor;

  Real _stress_matrix_contribution;

  Real _stress_rhs_contribution;

  const unsigned int _index;

  std::array<const MooseLinearVariableFVReal *, 3> _velocity_vars;

  const Moose::CoordinateSystemType _coord_type;

  const unsigned int _rz_radial_coord;

  const MooseLinearVariableFVReal & velocityVar(unsigned int dir) const;
};