PhaseNormalTensor.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 "PhaseNormalTensor.h"
 
registerMooseObject("PhaseFieldApp", PhaseNormalTensor);
 
InputParameters
PhaseNormalTensor::validParams()
{
  InputParameters params = Material::validParams();
  params.addClassDescription("Calculate normal tensor of a phase based on gradient");
  params.addRequiredCoupledVar("phase", "Phase variable");
  params.addRequiredParam<MaterialPropertyName>("normal_tensor_name", "Name of normal tensor");
  return params;
}
 
PhaseNormalTensor::PhaseNormalTensor(const InputParameters & parameters)
  : DerivativeMaterialInterface<Material>(parameters),
    _grad_u(coupledGradient("phase")),
    _normal_tensor(
        declareProperty<RankTwoTensor>(getParam<MaterialPropertyName>("normal_tensor_name")))
{
}
 
void
PhaseNormalTensor::initQpStatefulProperties()
{
  _normal_tensor[_qp].zero();
}
 
void
PhaseNormalTensor::computeQpProperties()
{
  const Real magnitude = _grad_u[_qp].norm();
 
  if (magnitude > 0.0)
  {
    const RealVectorValue vector = _grad_u[_qp] / magnitude;
    _normal_tensor[_qp] = RankTwoTensor::selfOuterProduct(vector);
  }
  else
    _normal_tensor[_qp].zero();
}