IsolatedBoundingBoxIC.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "IsolatedBoundingBoxIC.h"
#include "MooseMesh.h"

registerMooseObject("PhaseFieldApp", IsolatedBoundingBoxIC);

InputParameters
IsolatedBoundingBoxIC::validParams()
{
  InputParameters params = SmoothMultiBoundingBoxBaseIC::validParams();
  params.addClassDescription(
      "Specify variable values inside and outside a list of isolated boxes shaped "
      "axis-aligned regions defined by pairs of opposing corners");
  params.addParam<Real>("outside", 0.0, "The value of the variable outside the boxes");
  return params;
}

IsolatedBoundingBoxIC::IsolatedBoundingBoxIC(const InputParameters & parameters)
  : SmoothMultiBoundingBoxBaseIC(parameters), _outside(getParam<Real>("outside"))
{
}

Real
IsolatedBoundingBoxIC::value(const Point & p)
{
  Real value = SmoothMultiBoundingBoxBaseIC::value(p);

  if (_int_width != 0.0)
  {
    for (unsigned int b = 0; b < _nbox; ++b)
    {
      for (unsigned int i = 0; i < _dim; ++i)
      {
        if (_c1[b](i) < _c2[b](i) && p(i) >= _c1[b](i) - _int_width &&
            p(i) <= _c2[b](i) + _int_width)
        {
          if (i != _dim - 1)
            continue;
          for (unsigned int n = b + 1; n < _nbox; ++n)
          {
            for (unsigned int j = 0; j < _dim; ++j)
            {
              if (p(j) >= _c1[n](j) - _int_width && p(j) <= _c2[n](j) + _int_width)
                mooseError("Partially overlapping boxes are not allowed. Note that this "
                           "includes the overlapping diffused interfaces. For nested boxes, "
                           "use NestedBoundingBoxIC.C.");
            }
          }
          Real f_in = 1.0;
          for (unsigned int j = 0; j < _dim; ++j)
            f_in *= 0.5 * (std::tanh(2.0 * libMesh::pi * (p(j) - _c1[b](j)) / _int_width) -
                           std::tanh(2.0 * libMesh::pi * (p(j) - _c2[b](j)) / _int_width));
          value = _outside + (_inside[b] - _outside) * f_in;
        }
        else if (_c1[b](i) >= _c2[b](i))
          mooseError("The coordinates of the smaller_coordinate_corners are equal to or larger "
                     "than that of "
                     "the larger_coordinate_corners.");
        else
          break;
      }
    }
  }
  return value;
}