CrossTermBarrierFunctionMaterial adds free energy contribution on the interfaces between arbitrary pairs of phases in a symmetric way. More...

#include <CrossTermBarrierFunctionMaterial.h>

Inheritance diagram for CrossTermBarrierFunctionMaterial:

[legend]

Public Member Functions
	CrossTermBarrierFunctionMaterial (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	_g_order
	polynomial order of the switching function \( g(\eta) \) More...

std::vector< Real >	_W_ij
	barrier function height matrix More...

unsigned int	_num_eta
	order parameters More...

std::vector< VariableName >	_eta_names

std::vector< const VariableValue * >	_eta


MaterialProperty< Real > &	_prop_g
	Barrier function and its derivatives. More...

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

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


std::vector< const MaterialProperty< Real > * >	_h
	Switching functions and their derivatives. More...

std::vector< std::vector< const MaterialProperty< Real > * > >	_dh

Detailed Description

CrossTermBarrierFunctionMaterial adds free energy contribution on the interfaces between arbitrary pairs of phases in a symmetric way.

Definition at line 24 of file CrossTermBarrierFunctionMaterial.h.

Constructor & Destructor Documentation

◆ CrossTermBarrierFunctionMaterial()

CrossTermBarrierFunctionMaterial::CrossTermBarrierFunctionMaterial ( const InputParameters & parameters )

Definition at line 24 of file CrossTermBarrierFunctionMaterial.C.

   : CrossTermBarrierFunctionBase(parameters)
 {
   // error out if W_ij is not symmetric
   for (unsigned int i = 0; i < _num_eta; ++i)
     for (unsigned int j = 0; j < i; ++j)
       if (_W_ij[_num_eta * i + j] != _W_ij[_num_eta * j + i])
         paramError("W_ij", "Please supply a symmetric W_ij matrix");
 }

Member Function Documentation

◆ computeQpProperties()

void CrossTermBarrierFunctionMaterial::computeQpProperties ( )

protectedvirtual

Reimplemented from CrossTermBarrierFunctionBase.

Definition at line 36 of file CrossTermBarrierFunctionMaterial.C.

 {
   // Initialize properties to zero before accumulating
   CrossTermBarrierFunctionBase::computeQpProperties();
  
   // Sum the components of our W_ij matrix to get constant used in our g function
   for (unsigned int i = 0; i < _num_eta; ++i)
     for (unsigned int j = i + 1; j < _num_eta; ++j)
     {
       const Real ni = (*_eta[i])[_qp];
       const Real nj = (*_eta[j])[_qp];
       const Real Wij = _W_ij[_num_eta * i + j];
  
       switch (_g_order)
       {
         case 0: // SIMPLE
           _prop_g[_qp] += 16.0 * Wij * (ni * ni * nj * nj);
           // first derivatives
           (*_prop_dg[i])[_qp] += 16.0 * Wij * (2 * ni * nj * nj);
           (*_prop_dg[j])[_qp] += 16.0 * Wij * (2 * ni * ni * nj);
           // second derivatives (diagonal)
           (*_prop_d2g[i][i])[_qp] += 16.0 * Wij * (2 * nj * nj);
           (*_prop_d2g[j][j])[_qp] += 16.0 * Wij * (2 * ni * ni);
           // second derivatives (off-diagonal)
           (*_prop_d2g[i][j])[_qp] = 16.0 * Wij * (4 * ni * nj);
           break;
  
         case 1: // LOW
           _prop_g[_qp] += 4.0 * Wij * (ni * nj);
           // first derivatives
           (*_prop_dg[i])[_qp] += 4.0 * Wij * nj;
           (*_prop_dg[j])[_qp] += 4.0 * Wij * ni;
           // second derivatives (diagonal) vanish
           // second derivatives (off-diagonal)
           (*_prop_d2g[i][j])[_qp] = 4.0 * Wij;
           break;
  
         default:
           mooseError("Internal error");
       }
     }
 }