https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
BernoulliPressureVariable Class Referenceabstract

A special variable class for pressure which flags faces at which porosity jumps occur as extrapolated or Dirichlet boundary faces. More...

#include <BernoulliPressureVariable.h>

Inheritance diagram for BernoulliPressureVariable:
[legend]

Public Types

typedef typename MooseVariableField< OutputType >::OutputGradient OutputGradient
 
typedef typename MooseVariableField< OutputType >::OutputSecond OutputSecond
 
typedef typename MooseVariableField< OutputType >::OutputDivergence OutputDivergence
 
typedef typename MooseVariableField< OutputType >::FieldVariableValue FieldVariableValue
 
typedef typename MooseVariableField< OutputType >::FieldVariableGradient FieldVariableGradient
 
typedef typename MooseVariableField< OutputType >::FieldVariableSecond FieldVariableSecond
 
typedef typename MooseVariableField< OutputType >::FieldVariableCurl FieldVariableCurl
 
typedef typename MooseVariableField< OutputType >::FieldVariableDivergence FieldVariableDivergence
 
typedef typename MooseVariableField< OutputType >::OutputShape OutputShape
 
typedef typename MooseVariableField< OutputType >::OutputShapeGradient OutputShapeGradient
 
typedef typename MooseVariableField< OutputType >::OutputShapeSecond OutputShapeSecond
 
typedef typename MooseVariableField< OutputType >::OutputShapeDivergence OutputShapeDivergence
 
typedef typename MooseVariableField< OutputType >::FieldVariablePhiValue FieldVariablePhiValue
 
typedef typename MooseVariableField< OutputType >::FieldVariablePhiDivergence FieldVariablePhiDivergence
 
typedef typename MooseVariableField< OutputType >::FieldVariablePhiGradient FieldVariablePhiGradient
 
typedef typename MooseVariableField< OutputType >::FieldVariablePhiSecond FieldVariablePhiSecond
 
typedef Moose::ElemQpArg ElemQpArg
 
typedef Moose::ElemSideQpArg ElemSideQpArg
 
typedef Moose::ElemArg ElemArg
 
typedef Moose::FaceArg FaceArg
 
typedef Moose::StateArg StateArg
 
typedef Moose::NodeArg NodeArg
 
typedef Moose::ElemPointArg ElemPointArg
 
typedef MooseArray< std::vector< OutputShape > > FieldVariablePhiCurl
 
typedef MooseArray< std::vector< OutputShape > > FieldVariableTestValue
 
typedef MooseArray< std::vector< OutputShapeGradient > > FieldVariableTestGradient
 
typedef MooseArray< std::vector< OutputShapeSecond > > FieldVariableTestSecond
 
typedef MooseArray< std::vector< OutputShape > > FieldVariableTestCurl
 
typedef MooseArray< std::vector< OutputShapeDivergence > > FieldVariableTestDivergence
 
typedef typename MooseVariableDataBase< OutputType >::DofValue DofValue
 
typedef typename MooseVariableDataBase< OutputType >::DofValues DofValues
 
typedef typename MooseVariableDataBase< OutputType >::ADDofValue ADDofValue
 
typedef typename MooseVariableDataBase< OutputType >::ADDofValues ADDofValues
 
typedef DataFileName DataFileParameterType
 
typedef FunctorBase< Moose::ADType< OutputType >::type > FunctorType
 
typedef Moose::ADType< OutputType >::type ValueType
 
typedef typename FunctorReturnType< Moose::ADType< OutputType >::type, FunctorEvaluationKind::Gradient >::type GradientType
 
typedef ValueType DotType
 

Public Member Functions

 BernoulliPressureVariable (const InputParameters &params)
 
bool isExtrapolatedBoundaryFace (const FaceInfo &fi, const Elem *elem, const Moose::StateArg &time) const override
 
void initialSetup () override
 
void computeFaceValues (const FaceInfo &) override
 
void computeElemValues () override
 
void computeElemValuesFace () override
 
void computeNeighborValuesFace () override
 
void computeNeighborValues () override
 
void requireQpComputations () const override
 
bool supportsQpBasedLoops () const override
 
bool supportsGeometricInfoBasedLoops () const override
 
virtual void timestepSetup () override
 
virtual void meshChanged () override
 
bool isSeparatorBoundary (const FaceInfo &fi) const
 
virtual bool isExtrapolatedBoundaryFace (const FaceInfo &, const Elem *, const StateArg &) const
 
virtual bool isFV () const override
 
virtual void prepare () override final
 
virtual void prepareNeighbor () override final
 
virtual void prepareAux () override final
 
virtual void reinitNode () override final
 
virtual void reinitNodes (const std::vector< dof_id_type > &) override final
 
virtual void reinitNodesNeighbor (const std::vector< dof_id_type > &) override final
 
virtual void reinitAux () override final
 
virtual void reinitAuxNeighbor () override final
 
virtual void prepareLowerD () override final
 
virtual const dof_id_type & nodalDofIndex () const override final
 
virtual const dof_id_type & nodalDofIndexNeighbor () const override final
 
virtual std::size_t phiSize () const override final
 
virtual std::size_t phiFaceSize () const override final
 
virtual std::size_t phiNeighborSize () const override final
 
virtual std::size_t phiFaceNeighborSize () const override final
 
virtual std::size_t phiLowerSize () const override final
 
virtual void computeLowerDValues () override final
 
virtual void computeNodalNeighborValues () override final
 
virtual void computeNodalValues () override final
 
virtual const std::vector< dof_id_type > & dofIndicesLower () const override final
 
unsigned int numberOfDofs () const override final
 
virtual unsigned int numberOfDofsNeighbor () override final
 
virtual bool isNodal () const override final
 
bool hasDoFsOnNodes () const override final
 
libMesh::FEContinuity getContinuity () const override final
 
virtual bool isNodalDefined () const override final
 
virtual void setNodalValue (const OutputType &value, unsigned int idx=0) override
 
virtual void setDofValue (const DofValue &value, unsigned int index) override
 
void clearDofIndices () override
 
virtual void prepareIC () override
 
virtual const Elem *const & currentElem () const override
 
virtual void getDofIndices (const Elem *elem, std::vector< dof_id_type > &dof_indices) const override
 
virtual const std::vector< dof_id_type > & dofIndices () const final
 
virtual const std::vector< dof_id_type > & dofIndicesNeighbor () const final
 
Moose::FV::InterpMethod faceInterpolationMethod () const
 
void clearAllDofIndices () final
 
const DofValues & nodalVectorTagValue (TagID) const override
 
const DofValues & nodalMatrixTagValue (TagID) const override
 
const FieldVariableValue & vectorTagValue (TagID tag) const override
 
const DofValues & vectorTagDofValue (TagID tag) const override
 
const FieldVariableValue & matrixTagValue (TagID tag) const override
 
const FieldVariableValue & vectorTagValueNeighbor (TagID tag)
 
const FieldVariableValue & matrixTagValueNeighbor (TagID tag)
 
const FieldVariableValue & uDot () const
 
const FieldVariableValue & sln () const override
 
const FieldVariableValue & slnOld () const override
 
const FieldVariableValue & slnOlder () const override
 
const FieldVariableGradient & gradSln () const override
 
const FieldVariableGradient & gradSlnOld () const override
 
const FieldVariableValue & uDotNeighbor () const
 
const FieldVariableValue & slnNeighbor () const override
 
const FieldVariableValue & slnOldNeighbor () const override
 
const FieldVariableGradient & gradSlnNeighbor () const override
 
const FieldVariableGradient & gradSlnOldNeighbor () const override
 
const VariableValue & duDotDu () const
 
const VariableValue & duDotDotDu () const
 
const VariableValue & duDotDuNeighbor () const
 
const VariableValue & duDotDotDuNeighbor () const
 
const ADTemplateVariableValue< OutputType > & adSln () const override
 
const ADTemplateVariableGradient< OutputType > & adGradSln () const override
 
virtual const VectorValue< ADReal > & adGradSln (const Elem *const elem, const StateArg &state, const bool correct_skewness=false) const
 
virtual VectorValue< ADReal > adGradSln (const FaceInfo &fi, const StateArg &state, const bool correct_skewness=false) const
 
virtual VectorValue< ADReal > uncorrectedAdGradSln (const FaceInfo &fi, const StateArg &state, const bool correct_skewness=false) const
 
ADReal getBoundaryFaceValue (const FaceInfo &fi, const StateArg &state, bool correct_skewness=false) const
 
const ADTemplateVariableSecond< OutputType > & adSecondSln () const override
 
const ADTemplateVariableValue< OutputType > & adUDot () const override
 
const ADTemplateVariableValue< OutputType > & adUDotDot () const override
 
const ADTemplateVariableGradient< OutputType > & adGradSlnDot () const override
 
const ADTemplateVariableCurl< OutputType > & adCurlSln () const override
 
const ADTemplateVariableValue< OutputType > & adSlnNeighbor () const override
 
const ADTemplateVariableGradient< OutputType > & adGradSlnNeighbor () const override
 
const ADTemplateVariableSecond< OutputType > & adSecondSlnNeighbor () const override
 
const ADTemplateVariableValue< OutputType > & adUDotNeighbor () const override
 
const ADTemplateVariableValue< OutputType > & adUDotDotNeighbor () const override
 
const ADTemplateVariableGradient< OutputType > & adGradSlnNeighborDot () const override
 
const ADTemplateVariableCurl< OutputType > & adCurlSlnNeighbor () const override
 
virtual void setDofValues (const DenseVector< DofValue > &values) override
 
virtual void setLowerDofValues (const DenseVector< DofValue > &values) override
 
DofValue getElementalValue (const Elem *elem, unsigned int idx=0) const
 
DofValue getElementalValueOld (const Elem *elem, unsigned int idx=0) const
 
DofValue getElementalValueOlder (const Elem *elem, unsigned int idx=0) const
 
virtual void insert (libMesh::NumericVector< libMesh::Number > &vector) override
 
virtual void insertLower (libMesh::NumericVector< libMesh::Number > &vector) override
 
virtual void add (libMesh::NumericVector< libMesh::Number > &vector) override
 
const DofValues & dofValues () const override
 
const DofValues & dofValuesOld () const override
 
const DofValues & dofValuesOlder () const override
 
const DofValues & dofValuesPreviousNL () const override
 
const DofValues & dofValuesNeighbor () const override
 
const DofValues & dofValuesOldNeighbor () const override
 
const DofValues & dofValuesOlderNeighbor () const override
 
const DofValues & dofValuesPreviousNLNeighbor () const override
 
const DofValues & dofValuesDot () const override
 
const DofValues & dofValuesDotNeighbor () const override
 
const DofValues & dofValuesDotOld () const override
 
const DofValues & dofValuesDotOldNeighbor () const override
 
const DofValues & dofValuesDotDot () const override
 
const DofValues & dofValuesDotDotNeighbor () const override
 
const DofValues & dofValuesDotDotOld () const override
 
const DofValues & dofValuesDotDotOldNeighbor () const override
 
const MooseArray< libMesh::Number > & dofValuesDuDotDu () const override
 
const MooseArray< libMesh::Number > & dofValuesDuDotDuNeighbor () const override
 
const MooseArray< libMesh::Number > & dofValuesDuDotDotDu () const override
 
const MooseArray< libMesh::Number > & dofValuesDuDotDotDuNeighbor () const override
 
const ADDofValues & adDofValues () const override
 
const ADDofValues & adDofValuesNeighbor () const override
 
const ADDofValues & adDofValuesDot () const override
 
OutputType getValue (const Elem *elem) const
 
OutputTools< OutputType >::OutputGradient getGradient (const Elem *elem) const
 
bool hasDirichletBC () const
 
std::pair< bool, const FVDirichletBCBase * > getDirichletBC (const FaceInfo &fi) const
 
std::pair< bool, std::vector< const FVFluxBC * > > getFluxBCs (const FaceInfo &fi) const
 
virtual void residualSetup () override
 
virtual void jacobianSetup () override
 
ADReal getElemValue (const Elem *elem, const StateArg &state) const
 
void setActiveTags (const std::set< TagID > &vtags) override
 
bool supportsFaceArg () const override final
 
bool supportsElemSideQpArg () const override final
 
virtual void sizeMatrixTagData () override
 
const MooseArray< OutputType > & nodalValueArray () const override
 
const MooseArray< OutputType > & nodalValueOldArray () const override
 
const MooseArray< OutputType > & nodalValueOlderArray () const override
 
bool computingSecond () const override final
 
bool computingCurl () const override final
 
bool computingDiv () const override final
 
bool usesSecondPhiNeighbor () const override final
 
const FieldVariablePhiValue & phi () const override final
 
const FieldVariablePhiGradient & gradPhi () const override final
 
const FieldVariablePhiSecond & secondPhi () const override final
 
const FieldVariablePhiValue & curlPhi () const override final
 
const FieldVariablePhiDivergence & divPhi () const override final
 
const FieldVariablePhiValue & phiFace () const override final
 
const FieldVariablePhiGradient & gradPhiFace () const override final
 
const FieldVariablePhiSecond & secondPhiFace () const override final
 
const FieldVariablePhiValue & phiFaceNeighbor () const override final
 
const FieldVariablePhiGradient & gradPhiFaceNeighbor () const override final
 
const FieldVariablePhiSecond & secondPhiFaceNeighbor () const override final
 
const FieldVariablePhiValue & phiNeighbor () const override final
 
const FieldVariablePhiGradient & gradPhiNeighbor () const override final
 
const FieldVariablePhiSecond & secondPhiNeighbor () const override final
 
virtual const FieldVariablePhiValue & phiLower () const override
 
unsigned int oldestSolutionStateRequested () const override final
 
virtual ADReal getExtrapolatedBoundaryFaceValue (const FaceInfo &fi, bool two_term_expansion, bool correct_skewness, const Elem *elem_side_to_extrapolate_from, const StateArg &state) const
 
const bool & getTwoTermBoundaryExpansion () const
 
virtual Moose::VarFieldType fieldType () const override
 
virtual bool isArray () const override
 
virtual bool isVector () const override
 
bool usesPhiNeighbor () const
 
bool usesGradPhiNeighbor () const
 
bool hasBlocks (const SubdomainName &name) const
 
bool hasBlocks (const std::vector< SubdomainName > &names) const
 
bool hasBlocks (const std::set< SubdomainName > &names) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool hasBlocks (const SubdomainID id) const override
 
const std::string & componentName (const unsigned int comp) const
 
const std::set< SubdomainID > & activeSubdomains () const
 
bool activeOnSubdomain (SubdomainID subdomain) const
 
bool activeOnSubdomains (const std::set< SubdomainID > &subdomains) const
 
virtual bool needsGradientVectorStorage () const
 
const std::string & arrayVariableComponent (const unsigned int i) const
 
unsigned int number () const
 
const libMesh::FEType & feType () const
 
SystemBase & sys ()
 
const SystemBase & sys () const
 
bool useDual () const
 
const std::vector< dof_id_type > & allDofIndices () const
 
unsigned int totalVarDofs ()
 
Moose::VarKindType kind () const
 
void scalingFactor (const std::vector< Real > &factor)
 
Real scalingFactor () const
 
const std::vector< Real > & arrayScalingFactor () const
 
libMesh::Order order () const
 
unsigned int count () const
 
const libMesh::DofMap & dofMap () const
 
std::vector< dof_id_type > componentDofIndices (const std::vector< dof_id_type > &dof_indices, unsigned int component) const
 
bool eigen () const
 
void eigen (bool eigen)
 
virtual bool enabled () const
 
std::shared_ptr< MooseObject > getSharedPtr ()
 
std::shared_ptr< const MooseObject > getSharedPtr () const
 
bool isKokkosObject () const
 
MooseApp & getMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParameters & parameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
const std::vector< SubdomainName > & blocks () const
 
unsigned int numBlocks () const
 
virtual const std::set< SubdomainID > & blockIDs () const
 
unsigned int blocksMaxDimension () const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
void buildOutputHideVariableList (std::set< std::string > variable_names)
 
const std::set< OutputName > & getOutputs ()
 
virtual void subdomainSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
virtual void customSetup (const ExecFlagType &exec_type) override
 
const ExecFlagEnum & getExecuteOnEnum () const
 
FunctorReturnType< Moose::ADType< OutputType >::type, FET >::type genericEvaluate (const Space &r, const State &state) const
 
const MooseFunctorName & functorName () const
 
void setCacheClearanceSchedule (const std::set< ExecFlagType > &clearance_schedule)
 
bool isInternalFace (const FaceInfo &) const
 
virtual bool isConstant () const
 
virtual bool hasFaceSide (const FaceInfo &fi, const bool fi_elem_side) const override
 
void checkFace (const Moose::FaceArg &face) const
 
ValueType operator() (const ElemArg &elem, const StateArg &state) const
 
ValueType operator() (const FaceArg &face, const StateArg &state) const
 
ValueType operator() (const ElemQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemSideQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemPointArg &elem_point, const StateArg &state) const
 
ValueType operator() (const NodeArg &node, const StateArg &state) const
 
ValueType operator() (const ElemArg &elem, const StateArg &state) const
 
ValueType operator() (const FaceArg &face, const StateArg &state) const
 
ValueType operator() (const ElemQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemSideQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemPointArg &elem_point, const StateArg &state) const
 
ValueType operator() (const NodeArg &node, const StateArg &state) const
 
ValueType operator() (const ElemArg &elem, const StateArg &state) const
 
ValueType operator() (const FaceArg &face, const StateArg &state) const
 
ValueType operator() (const ElemQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemSideQpArg &qp, const StateArg &state) const
 
ValueType operator() (const ElemPointArg &elem_point, const StateArg &state) const
 
ValueType operator() (const NodeArg &node, const StateArg &state) const
 
GradientType gradient (const ElemArg &elem, const StateArg &state) const
 
GradientType gradient (const FaceArg &face, const StateArg &state) const
 
GradientType gradient (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradient (const NodeArg &node, const StateArg &state) const
 
GradientType gradient (const ElemArg &elem, const StateArg &state) const
 
GradientType gradient (const FaceArg &face, const StateArg &state) const
 
GradientType gradient (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradient (const NodeArg &node, const StateArg &state) const
 
GradientType gradient (const ElemArg &elem, const StateArg &state) const
 
GradientType gradient (const FaceArg &face, const StateArg &state) const
 
GradientType gradient (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradient (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradient (const NodeArg &node, const StateArg &state) const
 
DotType dot (const ElemArg &elem, const StateArg &state) const
 
DotType dot (const FaceArg &face, const StateArg &state) const
 
DotType dot (const ElemQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemSideQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemPointArg &elem_point, const StateArg &state) const
 
DotType dot (const NodeArg &node, const StateArg &state) const
 
DotType dot (const ElemArg &elem, const StateArg &state) const
 
DotType dot (const FaceArg &face, const StateArg &state) const
 
DotType dot (const ElemQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemSideQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemPointArg &elem_point, const StateArg &state) const
 
DotType dot (const NodeArg &node, const StateArg &state) const
 
DotType dot (const ElemArg &elem, const StateArg &state) const
 
DotType dot (const FaceArg &face, const StateArg &state) const
 
DotType dot (const ElemQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemSideQpArg &qp, const StateArg &state) const
 
DotType dot (const ElemPointArg &elem_point, const StateArg &state) const
 
DotType dot (const NodeArg &node, const StateArg &state) const
 
GradientType gradDot (const ElemArg &elem, const StateArg &state) const
 
GradientType gradDot (const FaceArg &face, const StateArg &state) const
 
GradientType gradDot (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradDot (const NodeArg &node, const StateArg &state) const
 
GradientType gradDot (const ElemArg &elem, const StateArg &state) const
 
GradientType gradDot (const FaceArg &face, const StateArg &state) const
 
GradientType gradDot (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradDot (const NodeArg &node, const StateArg &state) const
 
GradientType gradDot (const ElemArg &elem, const StateArg &state) const
 
GradientType gradDot (const FaceArg &face, const StateArg &state) const
 
GradientType gradDot (const ElemQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemSideQpArg &qp, const StateArg &state) const
 
GradientType gradDot (const ElemPointArg &elem_point, const StateArg &state) const
 
GradientType gradDot (const NodeArg &node, const StateArg &state) const
 
const MooseArray< GenericReal< is_ad > > & genericDofValues () const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (const std::string &name)
 

Static Public Member Functions

static InputParameters validParams ()
 
static void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
 
static std::string deduceFunctorName (const std::string &name, const InputParameters &params)
 

Public Attributes

 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 

Protected Member Functions

bool isDirichletBoundaryFace (const FaceInfo &fi, const Elem *elem, const Moose::StateArg &time) const override
 
ADReal getDirichletBoundaryFaceValue (const FaceInfo &fi, const Elem *elem, const Moose::StateArg &time) const override
 
std::pair< bool, ADRealVectorValue > elemIsUpwind (const Elem &elem, const FaceInfo &fi, const Moose::StateArg &time) const
 Checks to see whether the provided element is upwind of the provided face.
 
bool isFullyDevelopedFlowFace (const FaceInfo &fi) const
 Returns whether the passed-in FaceInfo corresponds to a fully-developed flow face.
 
void cacheSeparatorBoundaries ()
 Caches the separator boundaries.
 
void clearCaches ()
 
const libMesh::NumericVector< libMesh::Number > & getSolution (const Moose::StateArg &state) const
 
bool doDerivatives () const
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
void initializeKokkosBlockRestrictable ()
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockElements () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockNodes () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockSides () const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockElementID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockNodeID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION auto kokkosBlockElementSideID (Moose::Kokkos::ThreadID tid) const
 
virtual GradientType evaluateGradient (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const NodeArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const NodeArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const FaceArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradient (const NodeArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemSideQpArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemPointArg &, const StateArg &) const
 
virtual DotType evaluateDot (const NodeArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemSideQpArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemPointArg &, const StateArg &) const
 
virtual DotType evaluateDot (const NodeArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemArg &, const StateArg &) const
 
virtual DotType evaluateDot (const FaceArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemQpArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemSideQpArg &, const StateArg &) const
 
virtual DotType evaluateDot (const ElemPointArg &, const StateArg &) const
 
virtual DotType evaluateDot (const NodeArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const FaceArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const NodeArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const FaceArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const NodeArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const FaceArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemSideQpArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const ElemPointArg &, const StateArg &) const
 
virtual GradientType evaluateGradDot (const NodeArg &, const StateArg &) const
 
virtual ValueType evaluate (const ElemArg &elem, const StateArg &state) const =0
 
virtual ValueType evaluate (const FaceArg &face, const StateArg &state) const =0
 
virtual ValueType evaluate (const ElemQpArg &qp, const StateArg &state) const =0
 
virtual ValueType evaluate (const ElemSideQpArg &side_qp, const StateArg &state) const =0
 
virtual ValueType evaluate (const ElemPointArg &elem_point, const StateArg &state) const =0
 
virtual ValueType evaluate (const NodeArg &node, const StateArg &state) const =0
 
std::string deduceFunctorName (const std::string &name) const
 
const Moose::Functor< T > & getFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctor (const std::string &name, THREAD_ID tid)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
bool isFunctor (const std::string &name) const
 
bool isFunctor (const std::string &name, const SubProblem &subproblem) const
 
Moose::ElemArg makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
 
void checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name)
 

Protected Attributes

const Moose::Functor< ADReal > * _u
 The x-component of velocity.
 
const Moose::Functor< ADReal > * _v
 The y-component of velocity.
 
const Moose::Functor< ADReal > * _w
 The z-component of velocity.
 
const Moose::Functor< ADReal > * _eps
 The porosity.
 
const Moose::Functor< ADReal > * _rho
 The density.
 
std::vector< BoundaryName > _pressure_drop_sidesets
 The names of the sidesets which will have associated form loss coefficients.
 
const std::vector< BoundaryID > _pressure_drop_sideset_ids
 The boundary IDs corresponding to the form loss sidesets.
 
std::vector< Real > _pressure_drop_form_factors
 The form loss coefficients corresponding to the sidesets.
 
std::unordered_map< BoundaryID, const INSFVHydraulicSeparatorInterface * > _boundary_id_to_separator
 A container for quick access of hydraulic separator BCs associated with this variable.
 
 usingMooseVariableFieldMembers
 
std::unique_ptr< MooseVariableDataFV< OutputType > > _element_data
 
std::unique_ptr< MooseVariableDataFV< OutputType > > _neighbor_data
 
std::unordered_map< const Elem *, VectorValue< ADReal > > _elem_to_grad
 
bool _two_term_boundary_expansion
 
VectorValue< ADReal > _temp_cell_gradient
 
const bool _cache_cell_gradients
 
Moose::FV::InterpMethod _face_interp_method
 
SystemBase & _sys
 
libMesh::FEType _fe_type
 
unsigned int _var_num
 
unsigned int _index
 
bool _is_eigen
 
Moose::VarKindType _var_kind
 
SubProblem & _subproblem
 
const libMesh::Variable & _variable
 
Assembly & _assembly
 
const libMesh::DofMap & _dof_map
 
std::vector< dof_id_type > _dof_indices
 
MooseMesh & _mesh
 
THREAD_ID _tid
 
const unsigned int _count
 
std::vector< Real > _scaling_factor
 
bool _use_dual
 
bool _is_lower_d
 
std::vector< std::string > _array_var_component_names
 
const bool & _enabled
 
MooseApp & _app
 
Factory & _factory
 
ActionFactory & _action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters & _pars
 
const MaterialData * _blk_material_data
 
const ExecFlagEnum & _execute_enum
 
const ExecFlagType & _current_execute_flag
 
FEProblemBase & _mci_feproblem
 
const TimeIntegrator *const _time_integrator
 
ADReal _ad_real_dummy
 
const Parallel::Communicator & _communicator
 

Private Member Functions

GradientType evaluateGradient (const ElemQpArg &qp_arg, const StateArg &) const override final
 
GradientType evaluateGradient (const ElemArg &elem_arg, const StateArg &) const override final
 
GradientType evaluateGradient (const FaceArg &face, const StateArg &) const override final
 
DotType evaluateDot (const ElemArg &elem, const StateArg &) const override final
 
DotType evaluateDot (const FaceArg &face, const StateArg &) const override final
 
DotType evaluateDot (const ElemQpArg &elem_qp, const StateArg &) const override final
 
ValueType evaluate (const ElemArg &elem, const StateArg &) const override final
 
ValueType evaluate (const FaceArg &face, const StateArg &) const override final
 
ValueType evaluate (const NodeArg &node, const StateArg &) const override final
 
ValueType evaluate (const ElemPointArg &elem_point, const StateArg &state) const override final
 
ValueType evaluate (const ElemQpArg &elem_qp, const StateArg &state) const override final
 
ValueType evaluate (const ElemSideQpArg &elem_side_qp, const StateArg &state) const override final
 
void determineBoundaryToDirichletBCMap ()
 
void determineBoundaryToFluxBCMap ()
 
void lowerDError () const
 
void clearCacheData ()
 
ValueType queryQpCache (unsigned int qp, const libMesh::QBase &qrule, std::vector< std::pair< bool, Moose::ADType< OutputType >::type > > &qp_cache_data, const SpaceArg &space, const StateArg &state) const
 
ValueType queryFVArgCache (std::map< SpaceArg, ValueType > &cache_data, const SpaceArg &space) const
 
const Moose::Functor< T > * defaultFunctor (const std::string &name)
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
virtual bool isADObject () const override
 

Static Private Member Functions

static const hit::Node * getHitNode (const InputParameters &params)
 
static std::string messagePrefix (const InputParameters &params, const bool hit_prefix)
 

Private Attributes

const bool _allow_two_term_expansion_on_bernoulli_faces
 Switch to enable the two-term extrapolation on porosity jump faces.
 
bool _qp_calculations
 Whether to pre-initialize variable data for use in traditional MOOSE quadrature point based objects.
 
bool _dirichlet_map_setup
 
bool _flux_map_setup
 
const libMesh::NumericVector< libMesh::Number > *const & _solution
 
const FieldVariablePhiValue & _phi
 
const FieldVariablePhiGradient & _grad_phi
 
const FieldVariablePhiValue & _phi_face
 
const FieldVariablePhiGradient & _grad_phi_face
 
const FieldVariablePhiValue & _phi_face_neighbor
 
const FieldVariablePhiGradient & _grad_phi_face_neighbor
 
const FieldVariablePhiValue & _phi_neighbor
 
const FieldVariablePhiGradient & _grad_phi_neighbor
 
const Elem * _prev_elem
 
std::unordered_map< BoundaryID, const FVDirichletBCBase * > _boundary_id_to_dirichlet_bc
 
std::unordered_map< BoundaryID, std::vector< const FVFluxBC * > > _boundary_id_to_flux_bc
 
const ParallelParamObject & _parent
 
const MooseBase & _si_moose_base
 
const FEProblemBase * _si_problem
 
std::set< SubdomainID > _blk_ids
 
std::vector< SubdomainID > _vec_ids
 
std::vector< SubdomainName > _blocks
 
const bool _blk_dual_restrictable
 
FEProblemBase * _blk_feproblem
 
const MooseMesh * _blk_mesh
 
const std::set< BoundaryID > _empty_boundary_ids
 
const std::set< BoundaryID > & _boundary_ids
 
THREAD_ID _blk_tid
 
const std::string & _blk_name
 
unsigned int _blk_dim
 
const MooseObject * _moose_object
 
Moose::Kokkos::Array< ContiguousElementID > _kokkos_element_ids
 
Moose::Kokkos::Array< ContiguousNodeID > _kokkos_node_ids
 
Moose::Kokkos::Array< Moose::Kokkos::Pair< ContiguousElementID, unsigned int > > _kokkos_element_side_ids
 
MooseApp & _oi_moose_app
 
OutputWarehouse & _oi_output_warehouse
 
std::set< OutputName > _oi_outputs
 
ExecFlagEnum _empty_execute_enum
 
std::set< ExecFlagType > _clearance_schedule
 
bool _always_evaluate
 
dof_id_type _current_qp_map_key
 
std::vector< std::pair< bool, ValueType > > * _current_qp_map_value
 
std::unordered_map< dof_id_type, std::vector< std::pair< bool, ValueType > > > _qp_to_value
 
dof_id_type _current_side_qp_map_key
 
std::vector< std::vector< std::pair< bool, ValueType > > > * _current_side_qp_map_value
 
std::unordered_map< dof_id_type, std::vector< std::vector< std::pair< bool, ValueType > > > > _side_qp_to_value
 
std::map< ElemArg, ValueType > _elem_arg_to_value
 
std::map< FaceArg, ValueType > _face_arg_to_value
 
std::map< NodeArg, ValueType > _node_arg_to_value
 
MooseFunctorName _functor_name
 
const InputParameters & _fi_params
 
const std::string _fi_name
 
SubProblem *const _fi_subproblem
 
const THREAD_ID _fi_tid
 
std::vector< std::unique_ptr< Moose::Functor< Real > > > _default_real_functors
 
std::vector< std::unique_ptr< Moose::Functor< ADReal > > > _default_ad_real_functors
 

Detailed Description

A special variable class for pressure which flags faces at which porosity jumps occur as extrapolated or Dirichlet boundary faces.

The downwind is flagged as an extrapolated boundary face while the upwind face is flagged as a Dirichlet face. The upwind Dirichlet value is computed using the downwind extrapolated pressure value and the Bernoulli equation

Definition at line 23 of file BernoulliPressureVariable.h.

Constructor & Destructor Documentation

◆ BernoulliPressureVariable()

BernoulliPressureVariable::BernoulliPressureVariable ( const InputParameters &  params)

Definition at line 52 of file BernoulliPressureVariable.C.

53 : INSFVPressureVariable(params),
54 ADFunctorInterface(this),
55 _u(nullptr),
56 _v(nullptr),
57 _w(nullptr),
58 _eps(nullptr),
59 _rho(nullptr),
60 _pressure_drop_sidesets(getParam<std::vector<BoundaryName>>("pressure_drop_sidesets")),
61 _pressure_drop_sideset_ids(this->_mesh.getBoundaryIDs(_pressure_drop_sidesets)),
62 _pressure_drop_form_factors(getParam<std::vector<Real>>("pressure_drop_form_factors")),
63 _allow_two_term_expansion_on_bernoulli_faces(
64 getParam<bool>("allow_two_term_expansion_on_bernoulli_faces"))
65{
66 if (_allow_two_term_expansion_on_bernoulli_faces &&
67 _face_interp_method == Moose::FV::InterpMethod::SkewCorrectedAverage)
68 paramError(
69 "allow_two_term_expansion_on_bernoulli_faces",
70 "Skewness correction with two term extrapolation on Bernoulli faces is not supported!");
71
72 if (_pressure_drop_sidesets.size() != _pressure_drop_form_factors.size())
73 paramError("pressure_drop_sidesets",
74 "The number of sidesets and the number of supplied form losses are not equal!");
75}
BernoulliPressureVariable::_eps
const Moose::Functor< ADReal > * _eps
The porosity.
Definition BernoulliPressureVariable.h:63
BernoulliPressureVariable::_rho
const Moose::Functor< ADReal > * _rho
The density.
Definition BernoulliPressureVariable.h:65
BernoulliPressureVariable::_pressure_drop_form_factors
std::vector< Real > _pressure_drop_form_factors
The form loss coefficients corresponding to the sidesets.
Definition BernoulliPressureVariable.h:74
BernoulliPressureVariable::_pressure_drop_sidesets
std::vector< BoundaryName > _pressure_drop_sidesets
The names of the sidesets which will have associated form loss coefficients.
Definition BernoulliPressureVariable.h:68
BernoulliPressureVariable::_v
const Moose::Functor< ADReal > * _v
The y-component of velocity.
Definition BernoulliPressureVariable.h:59
BernoulliPressureVariable::_w
const Moose::Functor< ADReal > * _w
The z-component of velocity.
Definition BernoulliPressureVariable.h:61
BernoulliPressureVariable::_allow_two_term_expansion_on_bernoulli_faces
const bool _allow_two_term_expansion_on_bernoulli_faces
Switch to enable the two-term extrapolation on porosity jump faces.
Definition BernoulliPressureVariable.h:78
BernoulliPressureVariable::_u
const Moose::Functor< ADReal > * _u
The x-component of velocity.
Definition BernoulliPressureVariable.h:57
BernoulliPressureVariable::_pressure_drop_sideset_ids
const std::vector< BoundaryID > _pressure_drop_sideset_ids
The boundary IDs corresponding to the form loss sidesets.
Definition BernoulliPressureVariable.h:71
MooseVariableFV< Real >::paramError
void paramError(const std::string &param, Args... args) const
MooseVariableFV< Real >::getParam
const T & getParam(const std::string &name) const
MooseMesh::getBoundaryIDs
std::vector< BoundaryID > getBoundaryIDs(const Elem *const elem, const unsigned short int side) const
MooseVariableFV< Real >::_face_interp_method
Moose::FV::InterpMethod _face_interp_method

Member Function Documentation

◆ cacheSeparatorBoundaries()

void INSFVVariable::cacheSeparatorBoundaries ( )
protectedinherited

Caches the separator boundaries.

Definition at line 35 of file INSFVVariable.C.

36{
37 _boundary_id_to_separator.clear();
38 std::vector<FVFluxBC *> bcs;
39
40 const auto base_query = this->_subproblem.getMooseApp()
41 .theWarehouse()
42 .query()
43 .template condition<AttribSystem>("FVFluxBC")
44 .template condition<AttribThread>(_tid);
45
46 for (const auto bnd_id : this->_mesh.getBoundaryIDs())
47 {
48 auto base_query_copy = base_query;
49 base_query_copy.template condition<AttribBoundaries>(std::set<BoundaryID>({bnd_id}))
50 .queryInto(bcs);
51 for (const auto bc : bcs)
52 {
53 const auto separator = dynamic_cast<const INSFVHydraulicSeparatorInterface *>(bc);
54 if (separator)
55 _boundary_id_to_separator.emplace(bnd_id, separator);
56 }
57 }
58}
INSFVHydraulicSeparatorInterface
A base class which serves as a tag for hydraulic separators.
Definition INSFVHydraulicSeparatorInterface.h:19
INSFVVariable::_boundary_id_to_separator
std::unordered_map< BoundaryID, const INSFVHydraulicSeparatorInterface * > _boundary_id_to_separator
A container for quick access of hydraulic separator BCs associated with this variable.
Definition INSFVVariable.h:58
MooseApp::theWarehouse
TheWarehouse & theWarehouse()
SubProblem::getMooseApp
MooseApp & getMooseApp() const
MooseVariableFV< Real >::_subproblem
SubProblem & _subproblem
TheWarehouse::query
Query query()
MooseMeshUtils::getBoundaryIDs
std::vector< BoundaryID > getBoundaryIDs(const libMesh::MeshBase &mesh, const std::vector< BoundaryName > &boundary_name, bool generate_unknown, const std::set< BoundaryID > &mesh_boundary_ids)

Referenced by INSFVVariable::initialSetup(), INSFVVariable::meshChanged(), and INSFVVariable::timestepSetup().

◆ computeElemValues()

void INSFVVariable::computeElemValues ( )
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 75 of file INSFVVariable.h.

76{
77 if (_qp_calculations)
78 MooseVariableFVReal::computeElemValues();
79 else
80 _element_data->setGeometry(Moose::Volume);
81}
INSFVVariable::_qp_calculations
bool _qp_calculations
Whether to pre-initialize variable data for use in traditional MOOSE quadrature point based objects.
Definition INSFVVariable.h:64
MooseVariableFV< Real >::_element_data
std::unique_ptr< MooseVariableDataFV< OutputType > > _element_data
MooseVariableFV< Real >::computeElemValues
virtual void computeElemValues() override

◆ computeElemValuesFace()

void INSFVVariable::computeElemValuesFace ( )
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 84 of file INSFVVariable.h.

85{
86 if (_qp_calculations)
87 MooseVariableFVReal::computeElemValuesFace();
88}
MooseVariableFV< Real >::computeElemValuesFace
virtual void computeElemValuesFace() override

◆ computeFaceValues()

void INSFVVariable::computeFaceValues ( const FaceInfo &  fi)
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 68 of file INSFVVariable.h.

69{
70 if (_qp_calculations)
71 MooseVariableFVReal::computeFaceValues(fi);
72}
MooseVariableFV< Real >::computeFaceValues
virtual void computeFaceValues(const FaceInfo &fi) override

◆ computeNeighborValues()

void INSFVVariable::computeNeighborValues ( )
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 98 of file INSFVVariable.h.

99{
100 if (_qp_calculations)
101 MooseVariableFVReal::computeNeighborValues();
102}
MooseVariableFV< Real >::computeNeighborValues
virtual void computeNeighborValues() override

◆ computeNeighborValuesFace()

void INSFVVariable::computeNeighborValuesFace ( )
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 91 of file INSFVVariable.h.

92{
93 if (_qp_calculations)
94 MooseVariableFVReal::computeNeighborValuesFace();
95}
MooseVariableFV< Real >::computeNeighborValuesFace
virtual void computeNeighborValuesFace() override

◆ elemIsUpwind()

std::pair< bool, ADRealVectorValue > BernoulliPressureVariable::elemIsUpwind ( const Elem &  elem,
const FaceInfo &  fi,
const Moose::StateArg &  time 
) const
protected

Checks to see whether the provided element is upwind of the provided face.

Parameters
elemthe element to check whether it is upwind of the face
fithe face
timeThe time at which to evaluate the velocity
Returns
a pair in which the first member is whether the element is upwind of the face and the second member is the evaluated face superficial velocity

Definition at line 90 of file BernoulliPressureVariable.C.

93{
94 const Moose::FaceArg face{
95 &fi, Moose::FV::LimiterType::CentralDifference, true, false, nullptr, nullptr};
96
97 const VectorValue<ADReal> vel_face{(*_u)(face, time), (*_v)(face, time), (*_w)(face, time)};
98 const bool fi_elem_is_upwind = vel_face * fi.normal() > 0;
99 const bool elem_is_upwind = &elem == &fi.elem() ? fi_elem_is_upwind : !fi_elem_is_upwind;
100
101 return {elem_is_upwind, vel_face};
102}
FaceInfo::normal
const Point & normal() const
FaceInfo::elem
const Elem & elem() const

Referenced by getDirichletBoundaryFaceValue(), isDirichletBoundaryFace(), and isExtrapolatedBoundaryFace().

◆ getDirichletBoundaryFaceValue()

ADReal BernoulliPressureVariable::getDirichletBoundaryFaceValue ( const FaceInfo &  fi,
const Elem *  elem,
const Moose::StateArg &  time 
) const
overrideprotectedvirtual

Reimplemented from MooseVariableFV< Real >.

Definition at line 147 of file BernoulliPressureVariable.C.

150{
151 mooseAssert(isDirichletBoundaryFace(fi, elem, time), "This better be a Dirichlet face");
152
153 if (INSFVPressureVariable::isDirichletBoundaryFace(fi, elem, time))
154 return INSFVPressureVariable::getDirichletBoundaryFaceValue(fi, elem, time);
155
156 const auto [is_jump_face, eps_elem, eps_neighbor] = NS::isPorosityJumpFace(*_eps, fi, time);
157#ifndef NDEBUG
158 mooseAssert(is_jump_face,
159 "If we are not a traditional Dirichlet face, then we must be a jump face");
160#else
161 libmesh_ignore(is_jump_face);
162#endif
163
164 const Moose::FaceArg face_elem{
165 &fi, Moose::FV::LimiterType::CentralDifference, true, false, &fi.elem(), nullptr};
166 const Moose::FaceArg face_neighbor{
167 &fi, Moose::FV::LimiterType::CentralDifference, true, false, fi.neighborPtr(), nullptr};
168
169 const auto [elem_is_upwind, vel_face] = elemIsUpwind(*elem, fi, time);
170
171 const bool fi_elem_is_upwind = elem == &fi.elem() ? elem_is_upwind : !elem_is_upwind;
172 const auto & downwind_face = fi_elem_is_upwind ? face_neighbor : face_elem;
173
174 const auto & upwind_elem = makeElemArg(fi_elem_is_upwind ? fi.elemPtr() : fi.neighborPtr());
175 const auto & downwind_elem = makeElemArg(fi_elem_is_upwind ? fi.neighborPtr() : fi.elemPtr());
176
177 const auto & vel_upwind = vel_face;
178 const auto & vel_downwind = vel_face;
179
180 const auto & eps_upwind = fi_elem_is_upwind ? eps_elem : eps_neighbor;
181 const auto & eps_downwind = fi_elem_is_upwind ? eps_neighbor : eps_elem;
182
183 /*
184 If a weakly compressible material is used where the density slightly
185 depends on pressure. Given that the two-term expansion on Bernoulli faces is
186 enabled, we take use the downwind face assuming the continuity (minimal jump) of the
187 density. This protects against an infinite recursion between pressure and
188 density.
189 */
190 const auto & rho_upwind = (*_rho)(upwind_elem, time);
191 const auto & rho_downwind = (*_rho)(downwind_elem, time);
192
193 const VectorValue<ADReal> interstitial_vel_upwind = vel_upwind * (1 / eps_upwind);
194 const VectorValue<ADReal> interstitial_vel_downwind = vel_downwind * (1 / eps_downwind);
195
196 const auto v_dot_n_upwind = interstitial_vel_upwind * fi.normal();
197 const auto v_dot_n_downwind = interstitial_vel_downwind * fi.normal();
198
199 // Iterate through sidesets to see if we have associated irreversible pressure losses
200 // or not.
201 ADReal factor = 0.0;
202 for (const auto & bd_id : fi.boundaryIDs())
203 for (const auto i : index_range(_pressure_drop_sideset_ids))
204 if (_pressure_drop_sideset_ids[i] == bd_id)
205 factor += _pressure_drop_form_factors[i];
206
207 auto upwind_bernoulli_vel_chunk = 0.5 * rho_upwind * v_dot_n_upwind * v_dot_n_upwind;
208 auto downwind_bernoulli_vel_chunk = 0.5 * rho_downwind * v_dot_n_downwind * v_dot_n_downwind;
209
210 // We add the additional, irreversible pressure loss here.
211 // If it is a contraction we have to use the downwind values, otherwise the upwind values.
212 if (eps_downwind < eps_upwind)
213 downwind_bernoulli_vel_chunk +=
214 0.5 * factor * rho_downwind * v_dot_n_downwind * v_dot_n_downwind;
215 else
216 upwind_bernoulli_vel_chunk += -0.5 * factor * rho_upwind * v_dot_n_upwind * v_dot_n_upwind;
217
218 ADReal p_downwind;
219 if (_allow_two_term_expansion_on_bernoulli_faces)
220 p_downwind = (*this)(downwind_face, time);
221 else
222 p_downwind = (*this)(downwind_elem, time);
223
224 return p_downwind + downwind_bernoulli_vel_chunk - upwind_bernoulli_vel_chunk;
225}
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
for
for(PetscInt i=0;i< nvars;++i)
BernoulliPressureVariable::isDirichletBoundaryFace
bool isDirichletBoundaryFace(const FaceInfo &fi, const Elem *elem, const Moose::StateArg &time) const override
Definition BernoulliPressureVariable.C:126
BernoulliPressureVariable::elemIsUpwind
std::pair< bool, ADRealVectorValue > elemIsUpwind(const Elem &elem, const FaceInfo &fi, const Moose::StateArg &time) const
Checks to see whether the provided element is upwind of the provided face.
Definition BernoulliPressureVariable.C:90
FaceInfo::neighborPtr
const Elem * neighborPtr() const
FaceInfo::elemPtr
const Elem * elemPtr() const
ADFunctorInterface::makeElemArg
Moose::ElemArg makeElemArg(const Elem *elem, bool correct_skewnewss=false) const
MooseVariableFV< Real >::isDirichletBoundaryFace
virtual bool isDirichletBoundaryFace(const FaceInfo &fi, const Elem *elem, const Moose::StateArg &state) const
MooseVariableFV< Real >::getDirichletBoundaryFaceValue
virtual ADReal getDirichletBoundaryFaceValue(const FaceInfo &fi, const Elem *elem, const Moose::StateArg &state) const
NS::isPorosityJumpFace
std::tuple< bool, T, T > isPorosityJumpFace(const Moose::FunctorBase< T > &porosity, const FaceInfo &fi, const Moose::StateArg &time)
Checks to see whether the porosity value jumps from one side to the other of the provided face.
Definition NSFVUtils.C:119
libMesh::index_range
auto index_range(const T &sizable)
libMesh::libmesh_ignore
void libmesh_ignore(const Args &...)
libMesh::if
if(subdm)

◆ initialSetup()

void BernoulliPressureVariable::initialSetup ( )
overridevirtual

Reimplemented from INSFVVariable.

Definition at line 78 of file BernoulliPressureVariable.C.

79{
80 _u = &getFunctor<ADReal>("u", _subproblem);
81 _v = &getFunctor<ADReal>("v", _subproblem);
82 _w = &getFunctor<ADReal>("w", _subproblem);
83 _eps = &getFunctor<ADReal>(NS::porosity, _subproblem);
84 _rho = &getFunctor<ADReal>(NS::density, _subproblem);
85
86 INSFVPressureVariable::initialSetup();
87}
MooseVariableBase::initialSetup
void initialSetup() override
NS::density
static const std::string density
Definition NS.h:34
NS::porosity
static const std::string porosity
Definition NS.h:108

◆ isDirichletBoundaryFace()

bool BernoulliPressureVariable::isDirichletBoundaryFace ( const FaceInfo &  fi,
const Elem *  elem,
const Moose::StateArg &  time 
) const
overrideprotectedvirtual

Reimplemented from MooseVariableFV< Real >.

Definition at line 126 of file BernoulliPressureVariable.C.

129{
130 if (INSFVPressureVariable::isDirichletBoundaryFace(fi, elem, time))
131 return true;
132
133 if (isInternalFace(fi))
134 {
135 if (isSeparatorBoundary(fi))
136 return false;
137
138 if (std::get<0>(NS::isPorosityJumpFace(*_eps, fi, time)))
139 // We choose to apply the Dirichlet condition for the upwind element
140 return std::get<0>(elemIsUpwind(*elem, fi, time));
141 }
142
143 return false;
144}
INSFVVariable::isSeparatorBoundary
bool isSeparatorBoundary(const FaceInfo &fi) const
Definition INSFVVariable.C:82
MooseVariableFV< Real >::isInternalFace
bool isInternalFace(const FaceInfo &) const

Referenced by getDirichletBoundaryFaceValue(), and isExtrapolatedBoundaryFace().

◆ isExtrapolatedBoundaryFace()

bool BernoulliPressureVariable::isExtrapolatedBoundaryFace ( const FaceInfo &  fi,
const Elem *  elem,
const Moose::StateArg &  time 
) const
override

Definition at line 105 of file BernoulliPressureVariable.C.

108{
109 if (isDirichletBoundaryFace(fi, elem, time))
110 return false;
111 if (!isInternalFace(fi))
112 // We are neither a Dirichlet nor an internal face
113 return true;
114
115 if (isSeparatorBoundary(fi))
116 return true;
117
118 if (std::get<0>(NS::isPorosityJumpFace(*_eps, fi, time)))
119 // We choose to extrapolate for the element that is downwind
120 return !std::get<0>(elemIsUpwind(*elem, fi, time));
121
122 return false;
123}

◆ isFullyDevelopedFlowFace()

bool INSFVVariable::isFullyDevelopedFlowFace ( const FaceInfo &  fi) const
protectedinherited

Returns whether the passed-in FaceInfo corresponds to a fully-developed flow face.

Definition at line 110 of file INSFVVariable.C.

111{
112 const auto & face_type = fi.faceType(std::make_pair(this->number(), this->sys().number()));
113
114 mooseAssert(face_type != FaceInfo::VarFaceNeighbors::NEITHER,
115 "I'm concerned that if you're calling this method with a FaceInfo that doesn't have "
116 "this variable defined on either side, that you are doing something dangerous.");
117
118 // If we're defined on both sides of the face, then we're not a boundary face
119 if (face_type == FaceInfo::VarFaceNeighbors::BOTH)
120 return false;
121
122 std::vector<INSFVFullyDevelopedFlowBC *> bcs;
123
124 this->_subproblem.getMooseApp()
125 .theWarehouse()
126 .query()
127 .template condition<AttribINSFVBCs>(INSFVBCs::INSFVFullyDevelopedFlowBC)
128 .template condition<AttribBoundaries>(fi.boundaryIDs())
129 .queryInto(bcs);
130
131 return !bcs.empty();
132}
INSFVBCs::INSFVFullyDevelopedFlowBC
@ INSFVFullyDevelopedFlowBC
FaceInfo::faceType
VarFaceNeighbors faceType(const std::pair< unsigned int, unsigned int > &var_sys) const
FaceInfo::boundaryIDs
const std::set< BoundaryID > & boundaryIDs() const
MooseVariableFV< Real >::number
unsigned int number() const

Referenced by INSFVVelocityVariable::adGradSln(), INSFVVelocityVariable::getExtrapolatedBoundaryFaceValue(), and INSFVVelocityVariable::uncorrectedAdGradSln().

◆ isSeparatorBoundary()

bool INSFVVariable::isSeparatorBoundary ( const FaceInfo &  fi) const
inherited

Definition at line 82 of file INSFVVariable.C.

83{
84 for (const auto bid : fi.boundaryIDs())
85 if (_boundary_id_to_separator.count(bid))
86 return true;
87
88 return false;
89}
MooseVariableFV< Real >::count
unsigned int count() const

Referenced by isDirichletBoundaryFace(), isExtrapolatedBoundaryFace(), and INSFVVariable::isExtrapolatedBoundaryFace().

◆ meshChanged()

void INSFVVariable::meshChanged ( )
overridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 68 of file INSFVVariable.C.

69{
70 MooseVariableFVReal::meshChanged();
71 cacheSeparatorBoundaries();
72}
INSFVVariable::cacheSeparatorBoundaries
void cacheSeparatorBoundaries()
Caches the separator boundaries.
Definition INSFVVariable.C:35
MooseVariableFV< Real >::meshChanged
virtual void meshChanged() override

◆ requireQpComputations()

void INSFVVariable::requireQpComputations ( ) const
inlineoverridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Reimplemented in INSFVVelocityNoQpComputation.

Definition at line 30 of file INSFVVariable.h.

30{ _qp_calculations = true; }

◆ supportsGeometricInfoBasedLoops()

bool INSFVVariable::supportsGeometricInfoBasedLoops ( ) const
inlineoverridevirtualinherited

Reimplemented from MooseVariableField< typename OutputType >.

Definition at line 32 of file INSFVVariable.h.

32{ return true; }

◆ supportsQpBasedLoops()

bool INSFVVariable::supportsQpBasedLoops ( ) const
inlineoverridevirtualinherited

Reimplemented from MooseVariableField< typename OutputType >.

Definition at line 31 of file INSFVVariable.h.

31{ return _qp_calculations; }

◆ timestepSetup()

void INSFVVariable::timestepSetup ( )
overridevirtualinherited

Reimplemented from MooseVariableFV< Real >.

Definition at line 61 of file INSFVVariable.C.

62{
63 MooseVariableFVReal::timestepSetup();
64 cacheSeparatorBoundaries();
65}
MooseVariableFV< Real >::timestepSetup
virtual void timestepSetup() override

◆ validParams()

InputParameters BernoulliPressureVariable::validParams ( )
static

Definition at line 17 of file BernoulliPressureVariable.C.

18{
19 auto params = INSFVPressureVariable::validParams();
20 params += ADFunctorInterface::validParams();
21 params.addRequiredParam<MooseFunctorName>("u", "The x-component of velocity");
22 params.addParam<MooseFunctorName>("v", 0, "The y-component of velocity");
23 params.addParam<MooseFunctorName>("w", 0, "The z-component of velocity");
24 params.addRequiredParam<MooseFunctorName>(NS::porosity, "The porosity");
25 params.addRequiredParam<MooseFunctorName>(NS::density, "The density");
26 params.addParam<std::vector<BoundaryName>>(
27 "pressure_drop_sidesets", {}, "Sidesets over which form loss coefficients are to be applied");
28 params.addParam<std::vector<Real>>(
29 "pressure_drop_form_factors",
30 {},
31 "User-supplied form loss coefficients to be applied over the sidesets listed above");
32 params.addParam<bool>(
33 "allow_two_term_expansion_on_bernoulli_faces",
34 false,
35 "Switch to enable the two-term extrapolation on porosity jump faces. "
36 "WARNING: This might lead to crushes in parallel runs if porosity jump faces are connected "
37 "with one cell (usually corners) due to the insufficient number of ghosted "
38 "layers.");
39
40 // Assuming that this variable is used for advection problems, due to the
41 // utilization of the pressure gradient in the advecting velocity
42 // through the Rhie-Chow interpolation, we have to extend the ghosted layers.
43 params.addRelationshipManager(
44 "ElementSideNeighborLayers",
45 Moose::RelationshipManagerType::GEOMETRIC | Moose::RelationshipManagerType::ALGEBRAIC |
46 Moose::RelationshipManagerType::COUPLING,
47 [](const InputParameters & /*obj_params*/, InputParameters & rm_params)
48 { rm_params.set<unsigned short>("layers") = 3; });
49 return params;
50}
ADFunctorInterface::validParams
static InputParameters validParams()
INSFVPressureVariable::validParams
static InputParameters validParams()
Definition INSFVPressureVariable.C:15
InputParameters::set
T & set(const std::string &name, bool quiet_mode=false)

Member Data Documentation

◆ _allow_two_term_expansion_on_bernoulli_faces

const bool BernoulliPressureVariable::_allow_two_term_expansion_on_bernoulli_faces
private

Switch to enable the two-term extrapolation on porosity jump faces.

Definition at line 78 of file BernoulliPressureVariable.h.

Referenced by BernoulliPressureVariable(), and getDirichletBoundaryFaceValue().

◆ _boundary_id_to_separator

std::unordered_map<BoundaryID, const INSFVHydraulicSeparatorInterface *> INSFVVariable::_boundary_id_to_separator
protectedinherited

A container for quick access of hydraulic separator BCs associated with this variable.

Definition at line 58 of file INSFVVariable.h.

Referenced by INSFVVariable::cacheSeparatorBoundaries(), and INSFVVariable::isSeparatorBoundary().

◆ _eps

const Moose::Functor<ADReal>* BernoulliPressureVariable::_eps
protected

The porosity.

Definition at line 63 of file BernoulliPressureVariable.h.

Referenced by getDirichletBoundaryFaceValue(), initialSetup(), isDirichletBoundaryFace(), and isExtrapolatedBoundaryFace().

◆ _pressure_drop_form_factors

std::vector<Real> BernoulliPressureVariable::_pressure_drop_form_factors
protected

The form loss coefficients corresponding to the sidesets.

Definition at line 74 of file BernoulliPressureVariable.h.

Referenced by BernoulliPressureVariable(), and getDirichletBoundaryFaceValue().

◆ _pressure_drop_sideset_ids

const std::vector<BoundaryID> BernoulliPressureVariable::_pressure_drop_sideset_ids
protected

The boundary IDs corresponding to the form loss sidesets.

Definition at line 71 of file BernoulliPressureVariable.h.

Referenced by getDirichletBoundaryFaceValue().

◆ _pressure_drop_sidesets

std::vector<BoundaryName> BernoulliPressureVariable::_pressure_drop_sidesets
protected

The names of the sidesets which will have associated form loss coefficients.

Definition at line 68 of file BernoulliPressureVariable.h.

Referenced by BernoulliPressureVariable().

◆ _qp_calculations

bool INSFVVariable::_qp_calculations
mutableprivateinherited

Whether to pre-initialize variable data for use in traditional MOOSE quadrature point based objects.

This can be changed to true if objects request variable values in a way that is generated by quadrature point based calculations

Definition at line 64 of file INSFVVariable.h.

Referenced by INSFVVariable::computeElemValues(), INSFVVariable::computeElemValuesFace(), INSFVVariable::computeFaceValues(), INSFVVariable::computeNeighborValues(), INSFVVariable::computeNeighborValuesFace(), INSFVVariable::requireQpComputations(), and INSFVVariable::supportsQpBasedLoops().

◆ _rho

const Moose::Functor<ADReal>* BernoulliPressureVariable::_rho
protected

The density.

Definition at line 65 of file BernoulliPressureVariable.h.

Referenced by initialSetup().

◆ _u

const Moose::Functor<ADReal>* BernoulliPressureVariable::_u
protected

The x-component of velocity.

Definition at line 57 of file BernoulliPressureVariable.h.

Referenced by initialSetup().

◆ _v

const Moose::Functor<ADReal>* BernoulliPressureVariable::_v
protected

The y-component of velocity.

Definition at line 59 of file BernoulliPressureVariable.h.

Referenced by elemIsUpwind(), and initialSetup().

◆ _w

const Moose::Functor<ADReal>* BernoulliPressureVariable::_w
protected

The z-component of velocity.

Definition at line 61 of file BernoulliPressureVariable.h.

Referenced by elemIsUpwind(), and initialSetup().

The documentation for this class was generated from the following files:
Generated on Mon Aug 24 2026 19:11:22 for https://mooseframework.inl.gov by doxygen 1.9.8