Green-Gauss pressure gradient that removes prescribed jumps before interpolating across baffles. More...

#include <FVPressureJumpGreenGaussGradient.h>

Inheritance diagram for FVPressureJumpGreenGaussGradient:

Public Types
typedef std::vector< std::unique_ptr< libMesh::NumericVector< libMesh::Number > > >	GradientContainer

typedef DataFileName	DataFileParameterType

Public Member Functions
	FVPressureJumpGreenGaussGradient (const InputParameters &params)

void	linkFlowSystem (PorousRhieChowMassFlux &rc, const LinearFVGradientReader &pressure_gradient)
	Link this method to the porous Rhie-Chow object that supplies pressure jumps.

void	computeGradient (SystemBase &system, GradientContainer &gradient, const std::unordered_set< unsigned int > &variable_numbers) const

Moose::FV::GradientLimiterType	limiterType () const

virtual void	resolveGradientMethodDependencies (FEProblemBase &fe_problem)

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 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)

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 Types
typedef std::pair< Real, Real >	InternalFaceValues

Protected Member Functions
void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

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

Protected Attributes
const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

const Parallel::Communicator &	_communicator

Private Member Functions
void	computeGradientWithoutLimiter (SystemBase &system, GradientContainer &gradient, const std::unordered_set< unsigned int > &variable_numbers) const override

InternalFaceValues	internalFaceValues (const FaceInfo &fi, const ElemInfo &elem_info, const ElemInfo &neighbor_info, const MooseLinearVariableFV< Real > &variable, Real elem_value, Real neighbor_value) 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 PorousRhieChowMassFlux *	_rhie_chow = nullptr
	Porous Rhie-Chow object supplying baffle locations and signed pressure jumps after linkage.

const LinearFVGradientReader *	_pressure_gradient = nullptr
	Published gradient used to lag the explicit half-cell correction during reconstruction.

const SystemBase *	_pressure_system = nullptr
	Pressure system to which this method is linked.

unsigned int	_pressure_variable_number = libMesh::invalid_uint
	Pressure variable number within the linked system.

const Moose::FV::GradientLimiterType	_limiter_type

const ParallelParamObject &	_parent

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

Detailed Description

Green-Gauss pressure gradient that removes prescribed jumps before interpolating across baffles.

Definition at line 19 of file FVPressureJumpGreenGaussGradient.h.

Constructor & Destructor Documentation

◆ FVPressureJumpGreenGaussGradient()

FVPressureJumpGreenGaussGradient::FVPressureJumpGreenGaussGradient ( const InputParameters & params )

Definition at line 31 of file FVPressureJumpGreenGaussGradient.C.

  : FVGreenGaussGradient(params)
{
}

Member Function Documentation

◆ computeGradientWithoutLimiter()

void FVPressureJumpGreenGaussGradient::computeGradientWithoutLimiter	(	SystemBase &	system,
		GradientContainer &	gradient,
		const std::unordered_set< unsigned int > &	variable_numbers
	)		const

overrideprivatevirtual

Reimplemented from FVGreenGaussGradient.

Definition at line 61 of file FVPressureJumpGreenGaussGradient.C.

{
  if (variable_numbers.size() != 1)
    mooseError("must be used by exactly one pressure variable, but it is used by ",
               variable_numbers.size(),
               ".");
 
  if (_pressure_system)
  {
    mooseAssert(_pressure_system == &system,
                "FVPressureJumpGreenGaussGradient can only compute gradients for its linked "
                "pressure system.");
    mooseAssert(_pressure_variable_number == *variable_numbers.begin(),
                "FVPressureJumpGreenGaussGradient can only compute the linked pressure variable.");
  }
 
  FVGreenGaussGradient::computeGradientWithoutLimiter(system, gradient, variable_numbers);
}

◆ internalFaceValues()

FVGreenGaussGradient::InternalFaceValues FVPressureJumpGreenGaussGradient::internalFaceValues	(	const FaceInfo &	fi,
		const ElemInfo &	elem_info,
		const ElemInfo &	neighbor_info,
		const MooseLinearVariableFV< Real > &	variable,
		Real	elem_value,
		Real	neighbor_value
	)		const

overrideprivatevirtual

Reimplemented from FVGreenGaussGradient.

Definition at line 84 of file FVPressureJumpGreenGaussGradient.C.

{
  if (!_rhie_chow || !_rhie_chow->faceIsBaffle(fi))
    return FVGreenGaussGradient::internalFaceValues(
        fi, elem_info, neighbor_info, variable, elem_value, neighbor_value);
 
  mooseAssert(_pressure_gradient, "A linked pressure gradient is required on pressure-jump faces.");
 
  const Real jump = _rhie_chow->getSignedBaffleJump(fi, /*elem_side=*/true);
  if (!_rhie_chow->pressureDiffusionDataReady())
  {
    // Gradient history is initialized before the pressure-diffusion coefficients exist.
    // Preserve the jump during this startup pass using jump-adjusted interpolation.
    return {Moose::FV::linearInterpolation(elem_value, neighbor_value + jump, fi, true),
            Moose::FV::linearInterpolation(elem_value - jump, neighbor_value, fi, true)};
  }
 
  // Use the published geometric pressure gradient to match the deferred correction in the
  // pressure equation. The reconstructed coupling gradient must not enter this correction.
  const auto face_operator =
      _rhie_chow->pressureJumpInterfaceData(fi,
                                            _pressure_gradient->gradient(elem_info),
                                            _pressure_gradient->gradient(neighbor_info),
                                            fi.faceArea() * fi.faceCoord());
 
  // Eliminate the interface pressures to obtain one conservative mass flux, then recover the
  // one-sided face value used in each cell's Green-Gauss surface sum. Each side's H/A flux enters
  // its face value, so the face pressure carries the body force, e.g. the hydrostatic slope.
  const Real mass_flux = NS::FV::pressureJumpFlux(face_operator, elem_value, neighbor_value) -
                         NS::FV::pressureJumpHbyAFlux(face_operator);
  return {NS::FV::pressureJumpOneSidedFaceValue(
              face_operator, elem_value, mass_flux, /*elem_side=*/true),
          NS::FV::pressureJumpOneSidedFaceValue(
              face_operator, neighbor_value, mass_flux, /*elem_side=*/false)};
}

◆ linkFlowSystem()

void FVPressureJumpGreenGaussGradient::linkFlowSystem	(	PorousRhieChowMassFlux &	rc,
		const LinearFVGradientReader &	pressure_gradient
	)

Link this method to the porous Rhie-Chow object that supplies pressure jumps.

Definition at line 37 of file FVPressureJumpGreenGaussGradient.C.

{
  if (&pressure_gradient.method() != this)
    mooseError("Must be linked using a gradient field produced by that method.");
 
  if (pressure_gradient.system().number() != rc.pressureSystem().number() ||
      pressure_gradient.variableNumber() != rc.pressureVariableNumber())
    mooseError("Must be linked to the pressure variable owned by PorousRhieChowMassFlux '",
               rc.name(),
               "'.");
 
  if (_rhie_chow && _rhie_chow != &rc)
    mooseError("Already linked to PorousRhieChowMassFlux '",
               _rhie_chow->name(),
               "'. Use a separate gradient method for each flow system.");
 
  _rhie_chow = &rc;
  _pressure_gradient = &pressure_gradient;
  _pressure_system = &pressure_gradient.system();
  _pressure_variable_number = pressure_gradient.variableNumber();
}

◆ validParams()

InputParameters FVPressureJumpGreenGaussGradient::validParams ( )

static

Definition at line 20 of file FVPressureJumpGreenGaussGradient.C.

{
  InputParameters params = FVGreenGaussGradient::validParams();
  params.suppressParameter<MooseEnum>("limiter");
  params.addClassDescription(
      "Computes a Green-Gauss pressure gradient that excludes the pressure jump across baffles, "
      "using a one-sided face pressure for each cell, so the zero-thickness jump does not appear "
      "as a mesh-dependent gradient in the adjacent cells.");
  return params;
}