LevelSetVolume.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "LevelSetVolume.h"

registerMooseObject("LevelSetApp", LevelSetVolume);

InputParameters
LevelSetVolume::validParams()
{
  InputParameters params = ElementVariablePostprocessor::validParams();
  params.addClassDescription(
      "Compute the area or volume of the region inside or outside of a level set contour.");
  params.addParam<Real>(
      "threshold", 0.0, "The level set threshold to consider for computing area/volume.");

  MooseEnum loc("inside=0 outside=1", "inside");
  params.addParam<MooseEnum>("location", loc, "The location of the area/volume to be computed.");
  return params;
}

LevelSetVolume::LevelSetVolume(const InputParameters & parameters)
  : ElementVariablePostprocessor(parameters),
    _threshold(getParam<Real>("threshold")),
    _inside(getParam<MooseEnum>("location") == "inside")
{
}

void
LevelSetVolume::initialize()
{
  _volume = 0;
}

void
LevelSetVolume::execute()
{
  Real cnt = 0;
  Real n = _u.size();

  // Perform the check for inside/outside outside the qp loop for speed
  if (_inside)
  {
    for (_qp = 0; _qp < n; ++_qp)
      if (_u[_qp] <= _threshold)
        cnt++;
  }
  else
  {
    for (_qp = 0; _qp < n; ++_qp)
      if (_u[_qp] > _threshold)
        cnt++;
  }
  _volume += cnt / n * _current_elem_volume;
}

void
LevelSetVolume::finalize()
{
  gatherSum(_volume);
}

Real
LevelSetVolume::getValue() const
{
  return _volume;
}

void
LevelSetVolume::threadJoin(const UserObject & y)
{
  const auto & pps = static_cast<const LevelSetVolume &>(y);
  _volume += pps._volume;
}