PorousRhieChowMassFlux.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 "RhieChowMassFlux.h"
#include "FaceCenteredMapFunctor.h"
#include <unordered_set>
 
class PressureJumpModel;
namespace NS::FV
{
struct PressureJumpInterfaceData;
}
 
class PorousRhieChowMassFlux : public RhieChowMassFlux
{
public:
  static InputParameters validParams();
 
  PorousRhieChowMassFlux(const InputParameters & params);
 
  void
  linkMomentumPressureSystems(const std::vector<LinearSystem *> & momentum_systems,
                              LinearSystem & pressure_system,
                              const std::vector<unsigned int> & momentum_system_numbers) override;
 
  void computeHbyA(bool verbose) override;
 
  Real getFaceSidePorosity(const FaceInfo & fi, bool elem_side, const Moose::StateArg & time) const;
 
  Real getSignedBaffleJump(const FaceInfo & fi, bool elem_side) const;
 
  bool faceIsBaffle(const FaceInfo & fi) const;
 
  bool hasPressureBaffles() const { return !_pressure_jump_models.empty(); }
 
  RealVectorValue baffleFaceVelocity(const FaceInfo & fi, bool elem_side) const;
 
  Real baffleNormalVelocity(const FaceInfo & fi) const;
 
  Real cellPressureDiffusionCoefficient(const ElemInfo & elem_info, unsigned int component) const;
 
  bool pressureDiffusionDataReady() const;
 
  bool pressureDiffusionUsesNonorthogonalCorrection() const;
 
  NS::FV::PressureJumpInterfaceData
  pressureJumpInterfaceData(const FaceInfo & fi,
                            const RealVectorValue & elem_gradient,
                            const RealVectorValue & neighbor_gradient,
                            Real face_area) const;
 
  void initCouplingField() override;
 
  void meshChanged() override;
 
private:
  void setupPorousMeshInformation();
 
  void updateBaffleJumps();
 
  const PressureJumpModel * getPressureJumpModel(const FaceInfo & fi) const;
 
  const Moose::Functor<Real> & _eps;
 
  std::vector<const PressureJumpModel *> _pressure_jump_models;
 
  const Real _pressure_jump_relaxation;
 
  std::vector<const LinearFVGradientReader *> _velocity_gradients;
 
  FaceCenteredMapFunctor<Real, std::unordered_map<dof_id_type, Real>> & _baffle_jump;
};