Class implementing a flux boundary condition for linear finite volume pressure variables used in the pressure corrector equation which is consistent with the H/A flux and a prescribed boundary velocity. More...

#include <LinearFVPressureFluxBC.h>

Inheritance diagram for LinearFVPressureFluxBC:

[legend]

Public Types
enum	TEST_TYPE

enum	ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }

typedef DataFileName	DataFileParameterType

Public Member Functions
	LinearFVPressureFluxBC (const InputParameters &parameters)
	Class constructor. More...

virtual Real	computeBoundaryValue () const override

virtual Real	computeBoundaryNormalGradient () const override

virtual Real	computeBoundaryValueMatrixContribution () const override

virtual Real	computeBoundaryValueRHSContribution () const override

virtual Real	computeBoundaryGradientMatrixContribution () const override

virtual Real	computeBoundaryGradientRHSContribution () const override

virtual bool	includesMaterialPropertyMultiplier () const override

virtual bool	providesCompleteBoundaryFlux () const override

virtual bool	needsBoundaryNonorthogonalCorrection () const

virtual bool	hasFaceSide (const FaceInfo &fi, bool fi_elem_side) const override

const SubProblem &	subProblem () const

const MooseLinearVariableFV< Real > &	variable () const

void	setupFaceData (const FaceInfo *face_info, const FaceInfo::VarFaceNeighbors face_type)

const FaceInfo *	currentFaceInfo () const

FaceInfo::VarFaceNeighbors	currentFaceType () const

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

virtual const std::set< BoundaryID > &	boundaryIDs () const

const std::vector< BoundaryName > &	boundaryNames () const

unsigned int	numBoundaryIDs () const

bool	hasBoundary (const BoundaryName &name) const

bool	hasBoundary (const std::vector< BoundaryName > &names) const

bool	hasBoundary (const BoundaryID &id) const

bool	hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	isBoundarySubset (const std::set< BoundaryID > &ids) const

bool	isBoundarySubset (const std::vector< BoundaryID > &ids) const

bool	hasBoundaryMaterialProperty (const std::string &prop_name) const

virtual bool	boundaryRestricted () const

const std::set< BoundaryID > &	meshBoundaryIDs () const

virtual bool	checkVariableBoundaryIntegrity () const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (const FunctionName &name) const

bool	hasFunction (const std::string &param_name) const

bool	hasFunctionByName (const FunctionName &name) const

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

const T &	getKokkosFunction (const std::string &name) const

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

bool	hasKokkosFunction (const std::string &param_name) const

bool	hasKokkosFunctionByName (const FunctionName &name) const

UserObjectName	getUserObjectName (const std::string &param_name) const

const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const

const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const

const UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObjectBase &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const

bool	isImplicit ()

Moose::StateArg	determineState () const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const

PenetrationLocator &	getPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

PenetrationLocator &	getQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

NearestNodeLocator &	getNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

NearestNodeLocator &	getQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

bool	requiresGeometricSearch () const

virtual void	meshChanged ()

void	useVectorTag (const TagName &tag_name, VectorTagsKey)

void	useVectorTag (TagID tag_id, VectorTagsKey)

void	useMatrixTag (const TagName &tag_name, MatrixTagsKey)

void	useMatrixTag (TagID tag_id, MatrixTagsKey)

bool	isVectorTagged ()

bool	isMatrixTagged ()

bool	hasVectorTags () const

const std::set< TagID > &	getVectorTags (VectorTagsKey) const

const std::set< TagID > &	getMatrixTags (MatrixTagsKey) const

MooseVariableBase *	mooseVariableBase () const

MooseVariableField< Real > &	mooseVariableField ()

MooseVariableFE< Real > *	mooseVariable () const

MooseVariableFV< Real > *	mooseVariableFV () const

MooseLinearVariableFV< Real > *	mooseLinearVariableFV () const

const std::set< MooseVariableFieldBase *> &	getMooseVariableDependencies () const

std::set< MooseVariableFieldBase *>	checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})

virtual std::set< MooseVariableFieldBase *>	checkVariables (const libMesh::Node &node, const std::set< MooseVariableFieldBase * > &vars_to_check)

std::set< MooseVariableFieldBase *>	checkVariables (const libMesh::Elem &element, const std::set< MooseVariableFieldBase * > &vars_to_check)

void	addMooseVariableDependency (MooseVariableFieldBase *var)

void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)

Moose::FaceArg	makeFace (const FaceInfo &fi, const Moose::FV::LimiterType limiter_type, const bool elem_is_upwind, const bool correct_skewness=false, const Moose::StateArg *state_limiter=nullptr) const

Moose::FaceArg	makeCDFace (const FaceInfo &fi, const bool correct_skewness=false) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

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 bool	restricted (const std::set< BoundaryID > &ids)

static std::string	deduceFunctorName (const std::string &name, const InputParameters &params)

Public Attributes
	ALL

	ANY

	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
Real	computeRequiredPressureFlux () const
	Compute the required boundary pressure flux contribution. More...

Real	computeBoundaryAinv () const
	Compute the scalar A^{-1} coefficient, which is zero before the first momentum assembly. More...

std::string	deduceFunctorName (const std::string &name) const

Real	computeCellToFaceDistance () const

RealVectorValue	computeCellToFaceVector () const

Moose::FaceArg	singleSidedFaceArg (const FaceInfo *fi, Moose::FV::LimiterType limiter_type=Moose::FV::LimiterType::CentralDifference, bool correct_skewness=false) const

Moose::FaceArg	functorFaceArg (const FunctorType &functor, const FaceInfo *fi, Moose::FV::LimiterType limiter_type=Moose::FV::LimiterType::CentralDifference, bool correct_skewness=false) const

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

bool	hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const

KOKKOS_FUNCTION dof_id_type	numKokkosBoundaryNodes () const

KOKKOS_FUNCTION dof_id_type	numKokkosBoundarySides () const

KOKKOS_FUNCTION ContiguousNodeID	kokkosBoundaryNodeID (Moose::Kokkos::ThreadID tid) const

KOKKOS_FUNCTION auto	kokkosBoundaryElementSideID (Moose::Kokkos::ThreadID tid) const

virtual void	addUserObjectDependencyHelper (const UserObjectBase &) const

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &) const

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const

void	prepareVectorTag (Assembly &assembly, unsigned int ivar)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)

void	prepareVectorTagLower (Assembly &assembly, unsigned int ivar)

void	prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)

void	prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const

void	prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)

void	prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)

void	prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const

void	prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)

void	accumulateTaggedLocalResidual ()

void	assignTaggedLocalResidual ()

void	accumulateTaggedLocalMatrix ()

void	accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)

void	accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)

void	accumulateTaggedNonlocalMatrix ()

void	assignTaggedLocalMatrix ()

void	addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, const std::vector< Real > &scaling_factors)

void	addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResiduals (Assembly &assembly, const ADResidualsPacket &packet)

void	addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)

void	addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, const std::vector< Real > &scaling_factors)

void	addJacobian (Assembly &assembly, const ADResidualsPacket &packet)

void	addJacobian (Assembly &assembly, const DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)

void	addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)

void	setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)

void	setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)

void	setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)

virtual const OutputTools< Real >::VariableValue &	value ()

virtual const OutputTools< Real >::VariableValue &	valueOld ()

virtual const OutputTools< Real >::VariableValue &	valueOlder ()

virtual const OutputTools< Real >::VariableValue &	dot ()

virtual const OutputTools< Real >::VariableValue &	dotDot ()

virtual const OutputTools< Real >::VariableValue &	dotOld ()

virtual const OutputTools< Real >::VariableValue &	dotDotOld ()

virtual const VariableValue &	dotDu ()

virtual const VariableValue &	dotDotDu ()

virtual const OutputTools< Real >::VariableGradient &	gradient ()

virtual const OutputTools< Real >::VariableGradient &	gradientOld ()

virtual const OutputTools< Real >::VariableGradient &	gradientOlder ()

virtual const OutputTools< Real >::VariableSecond &	second ()

virtual const OutputTools< Real >::VariableSecond &	secondOld ()

virtual const OutputTools< Real >::VariableSecond &	secondOlder ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTest ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTestFace ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhi ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhiFace ()

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< Real > &	_HbyA_flux
	The H/A flux functor for this BC (can be variable, function, etc) More...

const Moose::Functor< RealVectorValue > &	_Ainv
	The functor for the 1/A tensor serving as a diffusion coefficient. More...

const unsigned short	_dim
	Spatial dimension of the mesh. More...

const bool	_two_term_expansion
	Whether to reconstruct the boundary pressure with the pressure flux and cell gradient. More...

const Moose::Functor< Real > &	_u
	Velocity functors used to prescribe a boundary mass flux. More...

const Moose::Functor< Real > *const	_v

const Moose::Functor< Real > *const	_w

const Moose::Functor< Real > &	_rho
	Density functor used with the prescribed boundary velocity. More...

const THREAD_ID	_tid

SubProblem &	_subproblem

MooseMesh &	_mesh

FEProblemBase &	_fv_problem

MooseLinearVariableFV< Real > &	_var

SystemBase &	_sys

const FaceInfo *	_current_face_info

FaceInfo::VarFaceNeighbors	_current_face_type

const unsigned int	_var_num

const unsigned int	_sys_num

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

GeometricSearchData &	_geometric_search_data

bool	_requires_geometric_search

FEProblemBase &	_mci_feproblem

DenseVector< Number >	_local_re

DenseMatrix< Number >	_local_ke

DenseMatrix< Number >	_nonlocal_ke

bool	_nodal

MooseVariableFE< Real > *	_variable

MooseVariableFV< Real > *	_fv_variable

MooseLinearVariableFV< Real > *	_linear_fv_variable

MooseVariableField< Real > *	_field_variable

Assembly *	_mvi_assembly

const Parallel::Communicator &	_communicator

Detailed Description

Class implementing a flux boundary condition for linear finite volume pressure variables used in the pressure corrector equation which is consistent with the H/A flux and a prescribed boundary velocity.

This allows the pressure boundary flux to be consistent with a non-zero boundary mass flux and the boundary pressure to be reconstructed for a diagonal anisotropic inverse momentum tensor.

Definition at line 21 of file LinearFVPressureFluxBC.h.

Constructor & Destructor Documentation

◆ LinearFVPressureFluxBC()

LinearFVPressureFluxBC::LinearFVPressureFluxBC ( const InputParameters & parameters )

Class constructor.

Parameters

parameters The InputParameters for the object

Definition at line 40 of file LinearFVPressureFluxBC.C.

   : LinearFVAdvectionDiffusionBC(parameters),
     _HbyA_flux(getFunctor<Real>("HbyA_flux")),
     _Ainv(getFunctor<RealVectorValue>("Ainv")),
     _dim(_subproblem.mesh().dimension()),
     _two_term_expansion(getParam<bool>("use_two_term_expansion")),
     _u(getFunctor<Real>("u")),
     _v(parameters.isParamValid("v") ? &getFunctor<Real>("v") : nullptr),
     _w(parameters.isParamValid("w") ? &getFunctor<Real>("w") : nullptr),
     _rho(getFunctor<Real>(NS::density))
 {
   if (_dim >= 2 && !_v)
     paramError("v", "The 'v' boundary velocity functor must be provided for 2D and 3D problems.");
 
   if (_dim >= 3 && !_w)
     paramError("w", "The 'w' boundary velocity functor must be provided for 3D problems.");
 
   if (_dim < 2 && _v)
     paramError("v", "The 'v' boundary velocity functor is only valid in 2D and 3D problems.");
 
   if (_dim < 3 && _w)
     paramError("w", "The 'w' boundary velocity functor is only valid in 3D problems.");
 
   if (_two_term_expansion)
     _var.computeCellGradients();
 }

Member Function Documentation

◆ computeBoundaryAinv()

Real LinearFVPressureFluxBC::computeBoundaryAinv ( ) const

protected

Compute the scalar A^{-1} coefficient, which is zero before the first momentum assembly.

Definition at line 94 of file LinearFVPressureFluxBC.C.

Referenced by computeBoundaryNormalGradient(), and computeBoundaryValue().

 {
   const auto face_arg = singleSidedFaceArg(_current_face_info);
   const auto face_ainv = _Ainv(face_arg, determineState());
   const auto & normal = _current_face_info->normal();
 
   // Match the boundary-normal coefficient used by LinearFVAnisotropicDiffusion:
   // for a diagonal tensor Ainv, the effective normal coefficient is n^T Ainv n.
   Real normal_ainv = 0.0;
   for (const auto i : make_range(_dim))
     normal_ainv += normal(i) * normal(i) * face_ainv(i);
 
   if (normal_ainv < 0.0)
     mooseError("The boundary-normal Ainv coefficient must be nonnegative, but its value is ",
                normal_ainv,
                ".");
 
   return normal_ainv;
 }

◆ computeBoundaryGradientMatrixContribution()

Real LinearFVPressureFluxBC::computeBoundaryGradientMatrixContribution ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 190 of file LinearFVPressureFluxBC.C.

 {
   return 0.0;
 }

◆ computeBoundaryGradientRHSContribution()

Real LinearFVPressureFluxBC::computeBoundaryGradientRHSContribution ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 196 of file LinearFVPressureFluxBC.C.

 {
   return -computeRequiredPressureFlux();
 }

◆ computeBoundaryNormalGradient()

Real LinearFVPressureFluxBC::computeBoundaryNormalGradient ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 144 of file LinearFVPressureFluxBC.C.

Referenced by computeBoundaryValue().

 {
   if (!_two_term_expansion)
     return 0.0;
 
   const Real normal_ainv = computeBoundaryAinv();
   if (normal_ainv == 0.0)
     return 0.0;
 
   const auto state = determineState();
   const auto face_ainv = _Ainv(singleSidedFaceArg(_current_face_info), state);
   const auto elem_info = _current_face_type == FaceInfo::VarFaceNeighbors::ELEM
                              ? _current_face_info->elemInfo()
                              : _current_face_info->neighborInfo();
   const Real boundary_normal_multiplier =
       _current_face_type == FaceInfo::VarFaceNeighbors::ELEM ? 1.0 : -1.0;
   const auto boundary_normal = boundary_normal_multiplier * _current_face_info->normal();
 
   RealVectorValue normal_scaled_ainv;
   for (const auto i : make_range(_dim))
     normal_scaled_ainv(i) = boundary_normal(i) * face_ainv(i);
 
   // The prescribed pressure flux is the complete tensor-weighted flux. Subtract its tangential
   // part only when inverting that flux to reconstruct the boundary-normal pressure gradient.
   const auto tangential_ainv = normal_scaled_ainv - normal_ainv * boundary_normal;
   const Real tangential_pressure_flux = tangential_ainv * _var.gradSln(*elem_info, state);
 
   return (-computeRequiredPressureFlux() - tangential_pressure_flux) / normal_ainv;
 }

◆ computeBoundaryValue()

Real LinearFVPressureFluxBC::computeBoundaryValue ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 115 of file LinearFVPressureFluxBC.C.

Referenced by computeBoundaryValueRHSContribution().

 {
   const auto state = determineState();
   const auto elem_info = _current_face_type == FaceInfo::VarFaceNeighbors::ELEM
                              ? _current_face_info->elemInfo()
                              : _current_face_info->neighborInfo();
 
   if (!_two_term_expansion)
     return _var.getElemValue(*elem_info, state);
 
   const Real normal_ainv = computeBoundaryAinv();
 
   // Ainv is initialized to zero and is first populated by the momentum assembly. Until then, use
   // the cell pressure (1 term expansion) for the initial pressure-gradient
   // calculation.
   if (normal_ainv == 0.0)
     return _var.getElemValue(*elem_info, state);
 
   const Real distance = computeCellToFaceDistance();
   const auto d_cf = computeCellToFaceVector();
   const auto & face_normal = _current_face_info->normal();
   const auto tangential_cell_to_face = d_cf - (d_cf * face_normal) * face_normal;
 
   // Return the 2-term expansion for the boundary value
   return _var.getElemValue(*elem_info, state) + computeBoundaryNormalGradient() * distance +
          _var.gradSln(*elem_info, state) * tangential_cell_to_face;
 }

◆ computeBoundaryValueMatrixContribution()

Real LinearFVPressureFluxBC::computeBoundaryValueMatrixContribution ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 175 of file LinearFVPressureFluxBC.C.

 {
   return 1.0;
 }

◆ computeBoundaryValueRHSContribution()

Real LinearFVPressureFluxBC::computeBoundaryValueRHSContribution ( ) const

overridevirtual

Implements LinearFVAdvectionDiffusionBC.

Definition at line 181 of file LinearFVPressureFluxBC.C.

 {
   const auto elem_info = _current_face_type == FaceInfo::VarFaceNeighbors::ELEM
                              ? _current_face_info->elemInfo()
                              : _current_face_info->neighborInfo();
   return computeBoundaryValue() - _var.getElemValue(*elem_info, determineState());
 }

◆ computeRequiredPressureFlux()

Real LinearFVPressureFluxBC::computeRequiredPressureFlux ( ) const

protected

Compute the required boundary pressure flux contribution.

Definition at line 68 of file LinearFVPressureFluxBC.C.

Referenced by computeBoundaryGradientRHSContribution(), and computeBoundaryNormalGradient().

 {
   const auto face_arg = singleSidedFaceArg(_current_face_info);
   const auto state = determineState();
 
   Real required_pressure_flux = _HbyA_flux(face_arg, state);
   const auto & normal = _current_face_info->normal();
   Real boundary_velocity_dot_normal = _u(face_arg, state) * normal(0);
 
   if (_dim >= 2)
     boundary_velocity_dot_normal += (*_v)(face_arg, state) * normal(1);
 
   if (_dim >= 3)
     boundary_velocity_dot_normal += (*_w)(face_arg, state) * normal(2);
 
   // FaceInfo normals point from element to neighbor, so reverse the prescribed velocity flux when
   // this boundary condition acts on the neighbor side of an internal face.
   const Real boundary_normal_multiplier =
       _current_face_type == FaceInfo::VarFaceNeighbors::ELEM ? 1.0 : -1.0;
   required_pressure_flux +=
       _rho(face_arg, state) * boundary_normal_multiplier * boundary_velocity_dot_normal;
 
   return required_pressure_flux;
 }

◆ includesMaterialPropertyMultiplier()

virtual bool LinearFVPressureFluxBC::includesMaterialPropertyMultiplier ( ) const

inlineoverridevirtual

Reimplemented from LinearFVAdvectionDiffusionBC.

Definition at line 44 of file LinearFVPressureFluxBC.h.

44 { return true; }

◆ providesCompleteBoundaryFlux()

virtual bool LinearFVPressureFluxBC::providesCompleteBoundaryFlux ( ) const

inlineoverridevirtual

Reimplemented from LinearFVAdvectionDiffusionBC.

Definition at line 46 of file LinearFVPressureFluxBC.h.

46 { return true; }

◆ validParams()

InputParameters LinearFVPressureFluxBC::validParams ( )

static

Definition at line 16 of file LinearFVPressureFluxBC.C.

 {
   InputParameters params = LinearFVAdvectionDiffusionBC::validParams();
   params.addClassDescription(
       "Adds a fixed diffusive flux BC which can be used for the assembly of linear "
       "finite volume system and whose normal face gradient values are determined "
       "using the H/A flux and a prescribed boundary velocity. This boundary condition is only "
       "designed to work with advection-diffusion problems.");
   params.addRequiredParam<MooseFunctorName>("HbyA_flux", "The total HbyA face flux value.");
   params.addRequiredParam<MooseFunctorName>(
       "Ainv", "The 1/A where A is the momentum system diagonal vector.");
   params.addParam<bool>(
       "use_two_term_expansion",
       true,
       "Whether to reconstruct the boundary pressure using the pressure flux and cell gradient. If "
       "false, the boundary pressure is approximated by the adjacent cell pressure.");
   params.addRequiredParam<MooseFunctorName>("u", "The x-velocity functor on the boundary.");
   params.addParam<MooseFunctorName>("v", "The y-velocity functor on the boundary.");
   params.addParam<MooseFunctorName>("w", "The z-velocity functor on the boundary.");
   params.addRequiredParam<MooseFunctorName>(
       NS::density, "The density functor used together with the prescribed boundary velocity.");
   return params;
 }

Member Data Documentation

◆ _Ainv

const Moose::Functor<RealVectorValue>& LinearFVPressureFluxBC::_Ainv

protected

The functor for the 1/A tensor serving as a diffusion coefficient.

Definition at line 59 of file LinearFVPressureFluxBC.h.

Referenced by computeBoundaryAinv(), and computeBoundaryNormalGradient().

◆ _dim

const unsigned short LinearFVPressureFluxBC::_dim

protected

Spatial dimension of the mesh.

Definition at line 62 of file LinearFVPressureFluxBC.h.

Referenced by computeBoundaryAinv(), computeBoundaryNormalGradient(), computeRequiredPressureFlux(), and LinearFVPressureFluxBC().

◆ _HbyA_flux

const Moose::Functor<Real>& LinearFVPressureFluxBC::_HbyA_flux

protected

The H/A flux functor for this BC (can be variable, function, etc)

Definition at line 56 of file LinearFVPressureFluxBC.h.

Referenced by computeRequiredPressureFlux().

◆ _rho

const Moose::Functor<Real>& LinearFVPressureFluxBC::_rho

protected

Density functor used with the prescribed boundary velocity.

Definition at line 73 of file LinearFVPressureFluxBC.h.

Referenced by computeRequiredPressureFlux().

◆ _two_term_expansion

const bool LinearFVPressureFluxBC::_two_term_expansion

protected

Whether to reconstruct the boundary pressure with the pressure flux and cell gradient.

Definition at line 65 of file LinearFVPressureFluxBC.h.

Referenced by computeBoundaryNormalGradient(), computeBoundaryValue(), and LinearFVPressureFluxBC().

◆ _u

const Moose::Functor<Real>& LinearFVPressureFluxBC::_u

protected

Velocity functors used to prescribe a boundary mass flux.

Definition at line 68 of file LinearFVPressureFluxBC.h.

Referenced by computeRequiredPressureFlux().

◆ _v

const Moose::Functor<Real>* const LinearFVPressureFluxBC::_v

protected

Definition at line 69 of file LinearFVPressureFluxBC.h.

Referenced by LinearFVPressureFluxBC().

◆ _w

const Moose::Functor<Real>* const LinearFVPressureFluxBC::_w

protected

Definition at line 70 of file LinearFVPressureFluxBC.h.

Referenced by LinearFVPressureFluxBC().

The documentation for this class was generated from the following files:

navier_stokes/include/linearfvbcs/LinearFVPressureFluxBC.h
navier_stokes/src/linearfvbcs/LinearFVPressureFluxBC.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LinearFVPressureFluxBC()

Member Function Documentation

◆ computeBoundaryAinv()

◆ computeBoundaryGradientMatrixContribution()

◆ computeBoundaryGradientRHSContribution()

◆ computeBoundaryNormalGradient()

◆ computeBoundaryValue()

◆ computeBoundaryValueMatrixContribution()

◆ computeBoundaryValueRHSContribution()

◆ computeRequiredPressureFlux()

◆ includesMaterialPropertyMultiplier()

◆ providesCompleteBoundaryFlux()

◆ validParams()

Member Data Documentation

◆ _Ainv

◆ _dim

◆ _HbyA_flux

◆ _rho

◆ _two_term_expansion

◆ _u

◆ _v

◆ _w