Tricrystal2CircleGrainsIC creates a 3 grain structure with 2 circle grains and one matrix grain. More...

#include <Tricrystal2CircleGrainsIC.h>

Inheritance diagram for Tricrystal2CircleGrainsIC:

Public Member Functions
	Tricrystal2CircleGrainsIC (const InputParameters &parameters)

virtual Real	value (const Point &p)

Protected Attributes
const MooseMesh &	_mesh

const unsigned int	_op_num

const unsigned int	_op_index

Point	_bottom_left

Point	_top_right

Point	_range

Detailed Description

Tricrystal2CircleGrainsIC creates a 3 grain structure with 2 circle grains and one matrix grain.

Definition at line 23 of file Tricrystal2CircleGrainsIC.h.

Constructor & Destructor Documentation

◆ Tricrystal2CircleGrainsIC()

Tricrystal2CircleGrainsIC::Tricrystal2CircleGrainsIC ( const InputParameters & parameters )

Definition at line 27 of file Tricrystal2CircleGrainsIC.C.

   : InitialCondition(parameters),
     _mesh(_fe_problem.mesh()),
     _op_num(getParam<unsigned int>("op_num")),
     _op_index(getParam<unsigned int>("op_index"))
 {
   if (_op_num != 3)
     paramError("op_num", "Tricrystal ICs must have op_num = 3");
  
   // Set up domain bounds with mesh tools
   for (unsigned int i = 0; i < LIBMESH_DIM; i++)
   {
     _bottom_left(i) = _mesh.getMinInDimension(i);
     _top_right(i) = _mesh.getMaxInDimension(i);
   }
   _range = _top_right - _bottom_left;
 }

Member Function Documentation

◆ value()

Real Tricrystal2CircleGrainsIC::value ( const Point & p )

virtual

Definition at line 46 of file Tricrystal2CircleGrainsIC.C.

 {
   Point grain_center_left;
   grain_center_left(0) = _bottom_left(0) + _range(0) / 4.0;
   grain_center_left(1) = _bottom_left(1) + _range(1) / 2.0;
   grain_center_left(2) = _bottom_left(2) + _range(2) / 2.0;
  
   Point grain_center_right;
   grain_center_right(0) = _bottom_left(0) + _range(0) * 3.0 / 4.0;
   grain_center_right(1) = _bottom_left(1) + _range(1) / 2.0;
   grain_center_right(2) = _bottom_left(2) + _range(2) / 2.0;
  
   Real radius = _range(0) / 5.0;
   Real dist_left = (p - grain_center_left).norm();
   Real dist_right = (p - grain_center_right).norm();
  
   if ((dist_left <= radius && _op_index == 1) || (dist_right <= radius && _op_index == 2) ||
       (dist_left > radius && dist_right > radius && _op_index == 0))
     return 1.0;
   else
     return 0.0;
 }