FVReconstructedPressureGradient.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 "FVGradientMethod.h"
#include "LinearFVGradientReader.h"
#include "MeshChangedInterface.h"
 
#include "libmesh/numeric_vector.h"
#include "libmesh/dense_matrix.h"
#include "libmesh/dense_vector.h"
 
#include <unordered_set>
 
class ElemInfo;
class FaceInfo;
class FEProblemBase;
class RhieChowMassFlux;
 
class FVReconstructedPressureGradient : public FVGradientMethod, public MeshChangedInterface
{
public:
  using GradientContainer = FVGradientMethod::GradientContainer;
  using GradientView = LinearFVGradientReader::GradientContainer;
 
  static InputParameters validParams();
 
  FVReconstructedPressureGradient(const InputParameters & params);
 
  void linkFlowSystem(RhieChowMassFlux & rc, const LinearFVGradientReader & pressure_gradient);
 
  void validateSetup(const RhieChowMassFlux & rc) const;
 
  const GradientMethodName & baseGradientMethodName() const { return _base_gradient_method_name; }
 
  void resetForTimeStep(const RhieChowMassFlux & rc);
 
  void saveLaggedVelocityGradient(RhieChowMassFlux & rc);
 
  void computeCandidateFromCorrectedFlux(const RhieChowMassFlux & rc);
 
  const GradientContainer & reconstructedCandidate(const RhieChowMassFlux & rc) const;
 
  void finalizeCouplingPressureGradient(const RhieChowMassFlux & rc,
                                        const GradientView & base_gradient);
 
  virtual void meshChanged() override;
 
private:
  enum class PressureGradientReconstructionState
  {
    NeedLaggedVelocityGradient,
    NeedCandidate,
    CandidateReady
  };
 
  enum class PressureGradientReconstructionEvent
  {
    Reset,
    SaveLaggedVelocityGradient,
    ReconstructCandidate,
    PublishCandidate
  };
 
  void computeGradientWithoutLimiter(
      SystemBase & system,
      GradientContainer & gradient,
      const std::unordered_set<unsigned int> & variable_numbers) const override;
 
  void resolveGradientMethodDependencies(FEProblemBase & fe_problem) override;
 
  void checkFlowSystem(const RhieChowMassFlux & rc) const;
 
  void copyGradient(const GradientView & source, GradientContainer & destination);
 
  void transition(PressureGradientReconstructionEvent event);
 
  void updateCouplingPressureGradient(const RhieChowMassFlux & rc,
                                      const GradientView & base_gradient,
                                      const GradientContainer & reconstructed_candidate);
 
  RealVectorValue reconstructionVelocityGradient(const RhieChowMassFlux & rc,
                                                 const RealVectorValue & elem_gradient,
                                                 const FaceInfo & fi,
                                                 bool elem_has_info,
                                                 unsigned int velocity_component) const;
 
  void assembleFaceProjection(const RhieChowMassFlux & rc,
                              const ElemInfo & elem_info,
                              const std::vector<RealVectorValue> & elem_velocity_gradients,
                              const FaceInfo * fi,
                              const Point & surface_vector,
                              bool elem_has_info,
                              DenseMatrix<Real> & matrix,
                              DenseVector<Real> & projection_rhs) const;
 
  DenseVector<Real> solveFaceProjection(const DenseMatrix<Real> & matrix,
                                        const DenseVector<Real> & projection_rhs) const;
 
  Real reconstructPressureGradient(const RhieChowMassFlux & rc,
                                   const ElemInfo & elem_info,
                                   unsigned int component,
                                   Real reconstructed_velocity) const;
 
  const GradientMethodName _base_gradient_method_name;
 
  const FVGradientMethod * _base_gradient_method = nullptr;
 
  const Real _gradient_relaxation;
 
  const RhieChowMassFlux * _rhie_chow = nullptr;
 
  const SystemBase * _pressure_system = nullptr;
 
  unsigned int _pressure_variable_number = libMesh::invalid_uint;
 
  std::vector<const SystemBase *> _momentum_systems;
 
  std::vector<const LinearFVGradientReader *> _velocity_gradient_fields;
 
  std::vector<std::vector<std::unique_ptr<NumericVector<Number>>>>
      _lagged_reconstruction_velocity_gradient;
 
  PressureGradientReconstructionState _reconstruction_state =
      PressureGradientReconstructionState::NeedLaggedVelocityGradient;
 
  dof_id_type _last_reconstructed_face_flux_iteration = 0;
 
  GradientContainer _reconstructed_pressure_gradient;
 
  GradientContainer _coupling_pressure_gradient;
};