InternalVolume Class Reference

This class computes the volume of an interior space. More...

#include <InternalVolume.h>

Inheritance diagram for InternalVolume:

Public Member Functions
	InternalVolume (const InputParameters &parameters)
void	initialSetup ()

Protected Member Functions
virtual Real	computeQpIntegral ()
virtual Real	getValue ()

Protected Attributes
const unsigned int	_component
const Real	_scale
const Function &	_addition

Detailed Description

This class computes the volume of an interior space.

This class relies on being handed sidesets that define a closed space.

If the sideset defines an interior surface, the volume reported will be positive. If the sideset defines an exterior surface, the volume reported will be negative. It is therefore possible to compute the net interior volume by including an interior and an exterior surface in the same sideset.

Definition at line 30 of file InternalVolume.h.

Constructor & Destructor Documentation

InternalVolume()

InternalVolume::InternalVolume

(

const InputParameters &

parameters

)

Definition at line 37 of file InternalVolume.C.

  : SideIntegralPostprocessor(parameters),
    _component(getParam<unsigned int>("component")),
    _scale(getParam<Real>("scale_factor")),
    _addition(getFunction("addition"))
{
}

Member Function Documentation

computeQpIntegral()

Real InternalVolume::computeQpIntegral ( )

protectedvirtual

Definition at line 105 of file InternalVolume.C.

{
  // Default scale factor is 1
  Real scale = 1.0;
  if (_coord_sys == Moose::COORD_RSPHERICAL)
  {
    // MOOSE will multiply by 4*pi*r*r
    scale = 1.0 / 3.0;
  }
  else if (_coord_sys == Moose::COORD_RZ && _component == 0)
  {
    // MOOSE will multiply by 2*pi*r
    // Will integrate over z giving 0.5*2*pi*r*r*height
    scale = 0.5;
  }
  else if (_coord_sys == Moose::COORD_RZ && _component == 1)
  {
    // MOOSE will multiply by 2*pi*r
    // Will integrate over r:
    // integral(2*pi*r*height) over r:
    // pi*r*r*height
    scale = 1.0;
  }
  return -scale * _q_point[_qp](_component) * _normals[_qp](_component);
}

getValue()

Real InternalVolume::getValue ( )

protectedvirtual

Definition at line 132 of file InternalVolume.C.

{
  Point p;
  return _scale * SideIntegralPostprocessor::getValue() + _addition.value(_t, p);
}

initialSetup()

void InternalVolume::initialSetup

(

)

Definition at line 72 of file InternalVolume.C.

{
  SideIntegralPostprocessor::initialSetup();
 
  const std::set<SubdomainID> & subdomains = _mesh.meshSubdomains();
  std::set<SubdomainID>::const_iterator iter = subdomains.begin();
  const std::set<SubdomainID>::const_iterator iter_end = subdomains.end();
  for (; iter != iter_end; ++iter)
  {
    const std::set<BoundaryID> & boundaries = _mesh.getSubdomainBoundaryIds(*iter);
    std::set<BoundaryID>::const_iterator bnd = boundaries.begin();
    const std::set<BoundaryID>::const_iterator bnd_end = boundaries.end();
    for (; bnd != bnd_end; ++bnd)
    {
      const std::set<BoundaryID> & b = boundaryIDs();
      std::set<BoundaryID>::const_iterator bit = b.begin();
      const std::set<BoundaryID>::const_iterator bit_end = b.end();
      for (; bit != bit_end; ++bit)
        if (*bit == *bnd)
        {
          Moose::CoordinateSystemType coord_sys = _fe_problem.getCoordSystem(*iter);
          if (_component != 0 && coord_sys == Moose::COORD_RSPHERICAL)
            mooseError("With spherical coordinates, the component must be 0 in InternalVolume.");
 
          if (_component > 1 && coord_sys == Moose::COORD_RZ)
            mooseError(
                "With cylindrical coordinates, the component must be 0 or 1 in InternalVolume.");
        }
    }
  }
}

Member Data Documentation

_addition

const Function& InternalVolume::_addition

protected

Definition at line 43 of file InternalVolume.h.

Referenced by getValue().

_component

const unsigned int InternalVolume::_component

protected

Definition at line 41 of file InternalVolume.h.

Referenced by computeQpIntegral(), and initialSetup().

_scale

const Real InternalVolume::_scale

protected

Definition at line 42 of file InternalVolume.h.

Referenced by getValue().

---

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

• misc/include/postprocessors/InternalVolume.h
• misc/src/postprocessors/InternalVolume.C