OPInterfaceBarrierMaterial is a Free Energy Penalty contribution material that acts on all of the eta_i variables to prevent more than two eta variables going above 0 on an interface. More...

#include <ThirdPhaseSuppressionMaterial.h>

Inheritance diagram for ThirdPhaseSuppressionMaterial:

[legend]

Public Member Functions
	ThirdPhaseSuppressionMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties ()

Protected Attributes
std::string	_function_name
	name of the function of eta (used to generate the material property names) More...

unsigned int	_num_eta
	order parameters More...

std::vector< const VariableValue * >	_eta

MaterialProperty< Real > &	_prop_g
	Barrier functions and their drivatives. More...

std::vector< MaterialProperty< Real > * >	_prop_dg

std::vector< std::vector< MaterialProperty< Real > * > >	_prop_d2g
	Material properties to store the second derivatives. More...

Detailed Description

OPInterfaceBarrierMaterial is a Free Energy Penalty contribution material that acts on all of the eta_i variables to prevent more than two eta variables going above 0 on an interface.

Definition at line 27 of file ThirdPhaseSuppressionMaterial.h.

Constructor & Destructor Documentation

◆ ThirdPhaseSuppressionMaterial()

ThirdPhaseSuppressionMaterial::ThirdPhaseSuppressionMaterial ( const InputParameters & parameters )

Definition at line 26 of file ThirdPhaseSuppressionMaterial.C.

   : DerivativeMaterialInterface<Material>(parameters),
     _function_name(getParam<std::string>("function_name")),
     _num_eta(coupledComponents("etas")),
     _eta(_num_eta),
     _prop_g(declareProperty<Real>(_function_name)),
     _prop_dg(_num_eta),
     _prop_d2g(_num_eta)
 {
   std::vector<std::string> eta_name(_num_eta);
   for (unsigned int i = 0; i < _num_eta; ++i)
   {
     _prop_d2g[i].resize(_num_eta);
     // declare derivative properties, fetch eta values
     eta_name[i] = getVar("etas", i)->name();
   }
   for (unsigned int i = 0; i < _num_eta; ++i)
   {
     _prop_dg[i] = &declarePropertyDerivative<Real>(_function_name, eta_name[i]);
     _eta[i] = &coupledValue("etas", i);
     for (unsigned int j = i; j < _num_eta; ++j)
     {
       _prop_d2g[i][j] = _prop_d2g[j][i] =
           &declarePropertyDerivative<Real>(_function_name, eta_name[i], eta_name[j]);
     }
   }
 }

Member Function Documentation

◆ computeQpProperties()

void ThirdPhaseSuppressionMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 55 of file ThirdPhaseSuppressionMaterial.C.

 {
   // Initialize properties to zero before accumulating
   _prop_g[_qp] = 0.0;
   for (unsigned int i = 0; i < _num_eta; ++i)
   {
     (*_prop_dg[i])[_qp] = 0.0;
     for (unsigned int j = i; j < _num_eta; ++j)
       (*_prop_d2g[i][j])[_qp] = 0.0;
   }
 
   // Create Interface barrier preventing interfaces involving more than two order parameters
   for (unsigned int i = 0; i < _num_eta; ++i)
     for (unsigned int j = 0; j < i; ++j)
       for (unsigned int k = 0; k < j; ++k)
       {
         const Real ni = (*_eta[i])[_qp];
         const Real nj = (*_eta[j])[_qp];
         const Real nk = (*_eta[k])[_qp];
 
         _prop_g[_qp] += ni * ni * nj * nj * nk * nk;
         (*_prop_dg[i])[_qp] += 2 * ni * nj * nj * nk * nk;
         (*_prop_dg[j])[_qp] += 2 * ni * ni * nj * nk * nk;
         (*_prop_dg[k])[_qp] += 2 * ni * ni * nj * nj * nk;
         (*_prop_d2g[i][i])[_qp] += 2 * nj * nj * nk * nk;
         (*_prop_d2g[j][j])[_qp] += 2 * ni * ni * nk * nk;
         (*_prop_d2g[k][k])[_qp] += 2 * ni * ni * nj * nj;
         (*_prop_d2g[i][j])[_qp] += 4 * ni * nj * nk * nk;
         (*_prop_d2g[i][k])[_qp] += 4 * ni * nj * nj * nk;
         (*_prop_d2g[k][j])[_qp] += 4 * ni * ni * nj * nk;
       }
 }