NestedBoundingBoxIC creates several nested boxes defined by their coordinates in the domain. More...

#include <NestedBoundingBoxIC.h>

Inheritance diagram for NestedBoundingBoxIC:

Public Member Functions
	NestedBoundingBoxIC (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

NestedBoundingBoxIC creates several nested 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 NestedBoundingBoxIC.h.

Constructor & Destructor Documentation

◆ NestedBoundingBoxIC()

NestedBoundingBoxIC::NestedBoundingBoxIC ( const InputParameters & parameters )

Definition at line 26 of file NestedBoundingBoxIC.C.

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

Member Function Documentation

◆ value()

Real NestedBoundingBoxIC::value ( const Point & p )

Definition at line 32 of file NestedBoundingBoxIC.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;
           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));
           if (b == _nbox - 1)
             value = _outside + (_inside[b] - _outside) * f_in;
           else
             value = _inside[b + 1] + (_inside[b] - _inside[b + 1]) * f_in;
           goto label;
         }
         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;
       }
     }
   }
 label:
   return value;
 }