SmoothCircleICLinearizedInterface.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 "SmoothCircleICLinearizedInterface.h"

registerMooseObject("PhaseFieldApp", SmoothCircleICLinearizedInterface);

InputParameters
SmoothCircleICLinearizedInterface::validParams()
{
  InputParameters params = SmoothCircleIC::validParams();
  params.addClassDescription(
      "Circle with a smooth interface transformed using the linearized interface function");
  params.addRequiredParam<Real>(
      "bound_value", "Bound value used to keep variable between +/-bound. Must be positive.");
  return params;
}

SmoothCircleICLinearizedInterface::SmoothCircleICLinearizedInterface(
    const InputParameters & parameters)
  : SmoothCircleIC(parameters), _bound(parameters.get<Real>("bound_value"))
{
}

Real
SmoothCircleICLinearizedInterface::value(const Point & p)
{
  // Transform bound into the order parameter space
  const Real trans_bound = (1.0 + std::tanh(-_bound / std::sqrt(2.0))) / 2.0;

  // Get IC value for the order parameter
  const Real trans_value = SmoothCircleIC::value(p);

  // Transform order parameter value into linearized interface variable
  Real value = 0.0;

  if (trans_value < trans_bound)
    value = -_bound;
  else if (trans_value > 1.0 - trans_bound)
    value = _bound;
  else
    value = std::sqrt(2.0) * std::atanh(2.0 * trans_value - 1.0);

  return value;
}