CrossIC Class Reference

CrossIC creates a C1 continuous initial condition that looks like a cross in the middle of the domain. More...

#include <CrossIC.h>

Inheritance diagram for CrossIC:

Public Member Functions
	CrossIC (const InputParameters &parameters)
virtual Real	value (const Point &p)
virtual RealGradient	gradient (const Point &p)

Public Attributes
const Real	_x1
const Real	_y1
const Real	_x2
const Real	_y2

Protected Member Functions
Number	interfaceValue (Real r)
Number	interfaceDerivative (Real r)

Protected Attributes
Real	_average
Real	_amplitude
Real	_length
Real	_width
Real	_buffer
Real	_interface

Detailed Description

CrossIC creates a C1 continuous initial condition that looks like a cross in the middle of the domain.

Definition at line 43 of file CrossIC.h.

Constructor & Destructor Documentation

CrossIC()

CrossIC::CrossIC

(

const InputParameters &

parameters

)

Definition at line 45 of file CrossIC.C.

  : C1ICBase(parameters),
    _x1(parameters.get<Real>("x1")),
    _y1(parameters.get<Real>("y1")),
    _x2(parameters.get<Real>("x2")),
    _y2(parameters.get<Real>("y2"))
{
}

Member Function Documentation

gradient()

RealGradient CrossIC::gradient ( const Point &  p )

virtual

Definition at line 269 of file CrossIC.C.

{
  Point pxminus = p, pxplus = p, pyminus = p, pyplus = p;
 
  pxminus(0) -= TOLERANCE;
  pyminus(1) -= TOLERANCE;
  pxplus(0) += TOLERANCE;
  pyplus(1) += TOLERANCE;
 
  Number uxminus = value(pxminus), uxplus = value(pxplus), uyminus = value(pyminus),
         uyplus = value(pyplus);
 
  return Gradient((uxplus - uxminus) / 2.0 / TOLERANCE, (uyplus - uyminus) / 2.0 / TOLERANCE);
}

interfaceDerivative()

Number C1ICBase::interfaceDerivative ( Real  r )

protectedinherited

Definition at line 71 of file C1ICBase.C.

{
  Real x = (r - _buffer) / _interface;
 
  if (x < 0.0)
    return 0.0;
  if (x > 1.0)
    return 0.0;
 
  return ((12.0 * x * x - 12.0 * x) * _amplitude);
}

interfaceValue()

Number C1ICBase::interfaceValue ( Real  r )

protectedinherited

Definition at line 58 of file C1ICBase.C.

{
  Real x = (r - _buffer) / _interface;
 
  if (x < 0.0)
    return (_average + _amplitude);
  if (x > 1.0)
    return (_average - _amplitude);
 
  return ((1.0 + 4.0 * x * x * x - 6.0 * x * x) * _amplitude + _average);
}

Referenced by value().

value()

Real CrossIC::value ( const Point &  p )

virtual

Definition at line 55 of file CrossIC.C.

{
  const Real x = (p(0) - _x1) / (_x2 - _x1);
  const Real y = (p(1) - _y1) / (_y2 - _y1);
 
  const Real cmax = _average + _amplitude;
  const Real cmin = _average - _amplitude;
 
  if (x > (0.5 + _length / 2.0 + _buffer + _interface))
  {
    return cmin;
  }
  else if (x > 0.5 + _length / 2.0)
  {
    if (y > 0.5 + _width / 2.0 + _buffer + _interface)
    {
      return cmin;
    }
    else if (y > 0.5 + _width / 2.0)
    {
      Real xd = x - .5 - _length / 2.0;
      Real yd = y - .5 - _width / 2.0;
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else if (y > 0.5 - _width / 2.0)
    {
      return interfaceValue(x - .5 - _length / 2.0);
    }
    else if (y > 0.5 - _width / 2. - _buffer - _interface)
    {
      Real xd = x - .5 - _length / 2.0;
      Real yd = y - .5 + _width / 2.0;
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else
    {
      return cmin;
    }
  }
  else if (x > 0.5 + _width / 2. + 2. * _buffer + _interface)
  {
    if (y > 0.5 + _width / 2 + _buffer + _interface)
      return cmin;
    else if (y > 0.5 + _width / 2.0)
      return interfaceValue(y - .5 - _width / 2.0);
    else if (y > 0.5 - _width / 2.0)
      return cmax;
    else if (y > 0.5 - _width / 2. - _buffer - _interface)
      return interfaceValue(.5 - _width / 2. - y);
    else
      return cmin;
  }
  else if (x > 0.5 + _width / 2.0)
  {
    if (y > 0.5 + _length / 2. + _buffer + _interface)
    {
      return cmin;
    }
    else if (y > 0.5 + _length / 2.0)
    {
      Real xd = x - (.5 + _width / 2.0);
      Real yd = y - (.5 + _length / 2.0);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else if (y > 0.5 + _width / 2. + 2. * _buffer + _interface)
    {
      return interfaceValue(x - .5 - _width / 2.0);
    }
    else if (y > 0.5 + _width / 2.0)
    {
      Real xd = x - (.5 + _width / 2. + 2. * _buffer + _interface);
      Real yd = y - (.5 + _width / 2. + 2. * _buffer + _interface);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(2. * _buffer + _interface - r);
    }
    else if (y > 0.5 - _width / 2.0)
    {
      return cmax;
    }
    else if (y > 0.5 - _width / 2. - 2. * _buffer - _interface)
    {
      Real xd = x - (.5 + _width / 2. + 2. * _buffer + _interface);
      Real yd = y - (.5 - _width / 2. - 2. * _buffer - _interface);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(2. * _buffer + _interface - r);
    }
    else if (y > 0.5 - _length / 2.0)
    {
      return interfaceValue(x - .5 - _width / 2.0);
    }
    else if (y > 0.5 - _length / 2. - _buffer - _interface)
    {
      Real xd = x - (.5 + _width / 2.0);
      Real yd = y - (.5 - _length / 2.0);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else
    {
      return cmin;
    }
  }
  else if (x > 0.5 - _width / 2.0)
  {
    if (y > 0.5 + _length / 2 + _buffer + _interface)
      return cmin;
    else if (y > 0.5 + _length / 2.0)
      return interfaceValue(y - .5 - _length / 2.0);
    else if (y > 0.5 - _length / 2.0)
      return cmax;
    else if (y > 0.5 - _length / 2. - _buffer - _interface)
      return interfaceValue(.5 - _length / 2. - y);
    else
      return cmin;
  }
  else if (x > 0.5 - _width / 2. - 2. * _buffer - _interface)
  {
    if (y > 0.5 + _length / 2. + _buffer + _interface)
    {
      return cmin;
    }
    else if (y > 0.5 + _length / 2.0)
    {
      Real xd = x - (.5 - _width / 2.0);
      Real yd = y - (.5 + _length / 2.0);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else if (y > 0.5 + _width / 2. + 2. * _buffer + _interface)
    {
      return interfaceValue(.5 - _width / 2. - x);
    }
    else if (y > 0.5 + _width / 2.0)
    {
      Real xd = x - (.5 - _width / 2. - 2. * _buffer - _interface);
      Real yd = y - (.5 + _width / 2. + 2. * _buffer + _interface);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(2. * _buffer + _interface - r);
    }
    else if (y > 0.5 - _width / 2.0)
    {
      return cmax;
    }
    else if (y > 0.5 - _width / 2. - 2. * _buffer - _interface)
    {
      Real xd = x - (.5 - _width / 2. - 2. * _buffer - _interface);
      Real yd = y - (.5 - _width / 2. - 2. * _buffer - _interface);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(2. * _buffer + _interface - r);
    }
    else if (y > 0.5 - _length / 2.0)
    {
      return interfaceValue(.5 - _width / 2. - x);
    }
    else if (y > 0.5 - _length / 2. - _buffer - _interface)
    {
      Real xd = x - (.5 - _width / 2.0);
      Real yd = y - (.5 - _length / 2.0);
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else
    {
      return cmin;
    }
  }
  else if (x > 0.5 - _length / 2.0)
  {
    if (y > 0.5 + _width / 2 + _buffer + _interface)
      return cmin;
    else if (y > 0.5 + _width / 2.0)
      return interfaceValue(y - .5 - _width / 2.0);
    else if (y > 0.5 - _width / 2.0)
      return cmax;
    else if (y > 0.5 - _width / 2. - _buffer - _interface)
      return interfaceValue(.5 - _width / 2. - y);
    else
      return cmin;
  }
  else if (x > (.5 - _length / 2. - _buffer - _interface))
  {
    if (y > 0.5 + _width / 2. + _buffer + _interface)
    {
      return cmin;
    }
    else if (y > 0.5 + _width / 2.0)
    {
      Real xd = x - (.5 - _length / 2.0);
      Real yd = y - .5 - _width / 2.0;
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else if (y > 0.5 - _width / 2.0)
    {
      return interfaceValue(.5 - _length / 2. - x);
    }
    else if (y > 0.5 - _width / 2. - _buffer - _interface)
    {
      Real xd = x - (.5 - _length / 2.0);
      Real yd = y - .5 + _width / 2.0;
      Real r = std::sqrt(xd * xd + yd * yd);
      return interfaceValue(r);
    }
    else
      return cmin;
  }
  else
    return cmin;
}

Referenced by gradient().

Member Data Documentation

_amplitude

Real C1ICBase::_amplitude

protectedinherited

Definition at line 52 of file C1ICBase.h.

Referenced by C1ICBase::interfaceDerivative(), C1ICBase::interfaceValue(), and value().

_average

Real C1ICBase::_average

protectedinherited

Definition at line 51 of file C1ICBase.h.

Referenced by C1ICBase::interfaceValue(), and value().

_buffer

Real C1ICBase::_buffer

protectedinherited

Definition at line 55 of file C1ICBase.h.

Referenced by C1ICBase::interfaceDerivative(), C1ICBase::interfaceValue(), and value().

_interface

Real C1ICBase::_interface

protectedinherited

Definition at line 56 of file C1ICBase.h.

Referenced by C1ICBase::interfaceDerivative(), C1ICBase::interfaceValue(), and value().

_length

Real C1ICBase::_length

protectedinherited

Definition at line 53 of file C1ICBase.h.

Referenced by value().

_width

Real C1ICBase::_width

protectedinherited

Definition at line 54 of file C1ICBase.h.

Referenced by value().

_x1

const Real CrossIC::_x1

Definition at line 51 of file CrossIC.h.

Referenced by value().

_x2

const Real CrossIC::_x2

Definition at line 53 of file CrossIC.h.

Referenced by value().

_y1

const Real CrossIC::_y1

Definition at line 52 of file CrossIC.h.

Referenced by value().

_y2

const Real CrossIC::_y2

Definition at line 54 of file CrossIC.h.

Referenced by value().

---

The documentation for this class was generated from the following files:

• phase_field/include/ics/CrossIC.h
• phase_field/src/ics/CrossIC.C