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 LinearFVGradientReader;
 
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;
 
  const LinearFVGradientReader * const _advected_gradient_field;
 
  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 LinearFVGradientReader * _gradient_field;
 
  std::array<const LinearFVGradientReader *, 3> _velocity_gradient_fields;
 
  const Moose::CoordinateSystemType _coord_type;
 
  const unsigned int _rz_radial_coord;
 
  const MooseLinearVariableFVReal & velocityVar(unsigned int dir) const;
 
  const LinearFVGradientReader & velocityGradientField(unsigned int dir) const;
};