IsolatedBoundingBoxIC creates several isolated boxes defined by their coordinates in the domain. More...

#include <IsolatedBoundingBoxIC.h>

Inheritance diagram for IsolatedBoundingBoxIC:

Public Member Functions
	IsolatedBoundingBoxIC (const InputParameters &parameters)

Real	value (const Point &p)

Protected Attributes
const Real	_outside
	values outside all the boxes More...

const unsigned int	_nbox
	number of boxes More...

const Real	_int_width
	value of interfacial width More...

const unsigned int	_dim
	dimensionality of the mesh More...

std::vector< Real >	_inside
	values inside the boxes More...


const std::vector< Point >	_c1
	lists of opposite corners More...

const std::vector< Point >	_c2

Detailed Description

IsolatedBoundingBoxIC creates several isolated boxes defined by their coordinates in the domain.

If int_width > zero, the border of the boxes smoothly transitions from the invalue to the outvalue.

Definition at line 25 of file IsolatedBoundingBoxIC.h.

Constructor & Destructor Documentation

◆ IsolatedBoundingBoxIC()

IsolatedBoundingBoxIC::IsolatedBoundingBoxIC ( const InputParameters & parameters )

Definition at line 27 of file IsolatedBoundingBoxIC.C.

   : SmoothMultiBoundingBoxBaseIC(parameters), _outside(getParam<Real>("outside"))
 {
 }

Member Function Documentation

◆ value()

Real IsolatedBoundingBoxIC::value ( const Point & p )

Definition at line 33 of file IsolatedBoundingBoxIC.C.

 {
   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;
 }