Double well phase transformation barrier free energy contribution. More...

#include <MultiBarrierFunctionMaterial.h>

Inheritance diagram for MultiBarrierFunctionMaterial:

Public Member Functions
	MultiBarrierFunctionMaterial (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...

MooseEnum	_g_order
	Polynomial order of the barrier function \( g(\eta) \). More...

bool	_well_only
	zero out g contribution in the eta interval [0:1] 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< MaterialProperty< Real > * >	_prop_d2g

Detailed Description

Double well phase transformation barrier free energy contribution.

\( g(\vec\eta) = \sum_i \eta_i^2(1-\eta_i)^2 \)

Definition at line 25 of file MultiBarrierFunctionMaterial.h.

Constructor & Destructor Documentation

◆ MultiBarrierFunctionMaterial()

MultiBarrierFunctionMaterial::MultiBarrierFunctionMaterial ( const InputParameters & parameters )

Definition at line 33 of file MultiBarrierFunctionMaterial.C.

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

Member Function Documentation

◆ computeQpProperties()

void MultiBarrierFunctionMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 55 of file MultiBarrierFunctionMaterial.C.

 {
   Real g = 0.0;
  
   for (unsigned int i = 0; i < _num_eta; ++i)
   {
     const Real n = (*_eta[i])[_qp];
  
     if (_well_only && n >= 0.0 && n <= 1.0)
     {
       (*_prop_dg[i])[_qp] = 0.0;
       (*_prop_d2g[i])[_qp] = 0.0;
       continue;
     }
  
     switch (_g_order)
     {
       case 0: // SIMPLE
         g += n * n * (1.0 - n) * (1.0 - n);
         (*_prop_dg[i])[_qp] = 2.0 * n * (n - 1.0) * (2.0 * n - 1.0);
         (*_prop_d2g[i])[_qp] = 12.0 * (n * n - n) + 2.0;
         break;
     }
   }
  
   _prop_g[_qp] = g;
 }