RZSymmetry Class Reference

Interface class for enabling objects to be RZ symmetric about arbitrary axis. More...

#include <RZSymmetry.h>

Inheritance diagram for RZSymmetry:

Public Member Functions
	RZSymmetry (const MooseObject *moose_object, const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeCircumference (const RealVectorValue &pt)

Protected Attributes
RealVectorValue	_axis_point
	A point on the axis of symmetry.
const RealVectorValue &	_axis_dir
	The direction of the axis of symmetry.
const Real &	_offset
	Radial offset of the axis of symmetry.

Detailed Description

Interface class for enabling objects to be RZ symmetric about arbitrary axis.

Definition at line 19 of file RZSymmetry.h.

Constructor & Destructor Documentation

RZSymmetry()

RZSymmetry::RZSymmetry	(	const MooseObject *	moose_object,
		const InputParameters &	parameters
	)

Definition at line 26 of file RZSymmetry.C.

  : _axis_point(0., 0., 0.),
    _axis_dir(parameters.get<RealVectorValue>("axis_dir")),
    _offset(parameters.get<Real>("offset"))
{
  const FEProblemBase * fe_problem =
      moose_object->isParamValid("_fe_problem_base")
          ? moose_object->getParam<FEProblemBase *>("_fe_problem_base")
          : nullptr;
  if (fe_problem != nullptr)
  {
    const MooseMesh & mesh = fe_problem->mesh();
 
    const BlockRestrictable * blk_restr = dynamic_cast<const BlockRestrictable *>(moose_object);
    const BoundaryRestrictable * bnd_restr =
        dynamic_cast<const BoundaryRestrictable *>(moose_object);
 
    if (blk_restr != nullptr)
    {
      const std::set<SubdomainID> & blks = blk_restr->meshBlockIDs();
 
      for (auto & b : blks)
      {
        if (fe_problem->getCoordSystem(b) != Moose::COORD_XYZ)
          mooseError(
              moose_object->name(),
              ": This is a THM-specific object can be applied only on subdomains with Cartesian "
              "coordinate system. If your domain has RZ cooridnate system associated with it, you "
              "need to use the object without the RZ suffix to obtain the desired result.");
      }
    }
    if (bnd_restr != nullptr)
    {
      std::set<SubdomainID> blks;
      const std::set<BoundaryID> & bnds = bnd_restr->boundaryIDs();
      for (auto & bid : bnds)
      {
        std::set<SubdomainID> conn_blks = mesh.getBoundaryConnectedBlocks(bid);
        for (auto & cb : conn_blks)
          blks.insert(cb);
      }
 
      for (auto & b : blks)
      {
        if (fe_problem->getCoordSystem(b) != Moose::COORD_XYZ)
          mooseError(
              moose_object->name(),
              ": This is a THM-specific object can be applied only on subdomains with Cartesian "
              "coordinate system. If your domain has RZ cooridnate system associated with it, you "
              "need to use the object without the RZ suffix to obtain the desired result.");
      }
    }
  }
 
  const Point pt = parameters.get<Point>("axis_point");
  _axis_point = Point(pt(0), pt(1), pt(2));
}

Member Function Documentation

computeCircumference()

Real RZSymmetry::computeCircumference ( const RealVectorValue &  pt )

protectedvirtual

Definition at line 85 of file RZSymmetry.C.

{
  RealVectorValue v = (pt - _axis_point);
  const Real r = v.cross(_axis_dir).norm() / _axis_dir.norm();
  return 2 * libMesh::pi * (r + _offset);
}

Referenced by LayeredAverageRZ::computeIntegral(), ADHeatStructureEnergyRZ::computeQpIntegral(), HeatStructureEnergyRZ::computeQpIntegral(), ADElementIntegralMaterialPropertyRZ::computeQpIntegral(), ADSideFluxIntegralRZ::computeQpIntegral(), FunctionElementIntegralRZ::computeQpIntegral(), FunctionSideIntegralRZ::computeQpIntegral(), HeatRateConductionRZ::computeQpIntegral(), HeatRateConvectionRZ::computeQpIntegral(), HeatRateExternalAppConvectionRZ::computeQpIntegral(), HeatRateHeatFluxRZ::computeQpIntegral(), HeatRateRadiationRZ::computeQpIntegral(), SideFluxIntegralRZ::computeQpIntegral(), ConvectionHeatTransferRZBC::computeQpJacobian(), ExternalAppConvectionHeatTransferRZBC::computeQpJacobian(), RadiativeHeatFluxRZBC::computeQpJacobian(), CoupledForceRZ::computeQpJacobian(), CoupledForceRZ::computeQpOffDiagJacobian(), ADHeatStructureHeatSourceRZ::computeQpResidual(), ADConvectionHeatTransferRZBC::computeQpResidual(), ADExternalAppConvectionHeatTransferRZBC::computeQpResidual(), ADHSHeatFluxRZBC::computeQpResidual(), ADRadiativeHeatFluxRZBC::computeQpResidual(), ConvectionHeatTransferRZBC::computeQpResidual(), ExternalAppConvectionHeatTransferRZBC::computeQpResidual(), HeatStructure2DCouplerRZBC::computeQpResidual(), HeatStructure2DRadiationCouplerRZBC::computeQpResidual(), HSCoupler2D2DRadiationRZBC::computeQpResidual(), RadiativeHeatFluxRZBC::computeQpResidual(), CoupledForceRZ::computeQpResidual(), LayeredAverageRZ::execute(), ADHeatConductionRZ::precomputeQpResidual(), and ADHeatConductionTimeDerivativeRZ::precomputeQpResidual().

validParams()

InputParameters RZSymmetry::validParams ( )

static

Definition at line 17 of file RZSymmetry.C.

{
  InputParameters params = emptyInputParameters();
  params.addRequiredParam<Point>("axis_point", "A point on the axis of RZ symmetry.");
  params.addRequiredParam<RealVectorValue>("axis_dir", "The direction of the axis of RZ symmetry.");
  params.addParam<Real>("offset", 0., "Radial offset of the axis of RZ symmetry.");
  return params;
}

Referenced by ADConvectionHeatTransferRZBC::validParams(), ADExternalAppConvectionHeatTransferRZBC::validParams(), ADHSHeatFluxRZBC::validParams(), ADRadiativeHeatFluxRZBC::validParams(), ConvectionHeatTransferRZBC::validParams(), ExternalAppConvectionHeatTransferRZBC::validParams(), HeatStructure2DCouplerRZBC::validParams(), HeatStructure2DRadiationCouplerRZBC::validParams(), HSCoupler2D2DRadiationRZBC::validParams(), RadiativeHeatFluxRZBC::validParams(), ADHeatConductionRZ::validParams(), ADHeatConductionTimeDerivativeRZ::validParams(), ADHeatStructureHeatSourceRZ::validParams(), CoupledForceRZ::validParams(), ADElementIntegralMaterialPropertyRZ::validParams(), ADHeatStructureEnergyRZ::validParams(), ADSideFluxIntegralRZ::validParams(), FunctionElementIntegralRZ::validParams(), FunctionSideIntegralRZ::validParams(), HeatRateConductionRZ::validParams(), HeatRateConvectionRZ::validParams(), HeatRateExternalAppConvectionRZ::validParams(), HeatRateHeatFluxRZ::validParams(), HeatRateRadiationRZ::validParams(), HeatStructureEnergyRZ::validParams(), SideFluxIntegralRZ::validParams(), and LayeredAverageRZ::validParams().

Member Data Documentation

_axis_dir

const RealVectorValue& RZSymmetry::_axis_dir

protected

The direction of the axis of symmetry.

Definition at line 30 of file RZSymmetry.h.

Referenced by computeCircumference().

_axis_point

RealVectorValue RZSymmetry::_axis_point

protected

A point on the axis of symmetry.

Definition at line 28 of file RZSymmetry.h.

Referenced by computeCircumference(), and RZSymmetry().

_offset

const Real& RZSymmetry::_offset

protected

Radial offset of the axis of symmetry.

Definition at line 32 of file RZSymmetry.h.

Referenced by computeCircumference().

---

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

• thermal_hydraulics/include/base/RZSymmetry.h
• thermal_hydraulics/src/base/RZSymmetry.C