Material class to provide the switching function \( h(\eta) \) for the KKS system. More...

#include <SwitchingFunctionMaterial.h>

Inheritance diagram for SwitchingFunctionMaterial:

Public Member Functions
	SwitchingFunctionMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties ()

Protected Attributes
MooseEnum	_h_order
	Polynomial order of the switching function \( h(\eta) \). More...

const VariableValue &	_eta
	Coupled variable value for the order parameter \( \eta \). More...

unsigned int	_eta_var

VariableName	_eta_name

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

MaterialProperty< Real > &	_prop_f
	Material property to store \( f(\eta) \). More...

MaterialProperty< Real > &	_prop_df
	Material property to store the derivative \( df(\eta)/d\eta \). More...

MaterialProperty< Real > &	_prop_d2f
	Material property to store the second derivative \( d^2f(\eta)/d\eta^2 \). More...

Detailed Description

Material class to provide the switching function \( h(\eta) \) for the KKS system.

See also: KKSPhaseChemicalPotential; KKSCHBulk

Definition at line 27 of file SwitchingFunctionMaterial.h.

Constructor & Destructor Documentation

◆ SwitchingFunctionMaterial()

SwitchingFunctionMaterial::SwitchingFunctionMaterial ( const InputParameters & parameters )

Definition at line 29 of file SwitchingFunctionMaterial.C.

   : OrderParameterFunctionMaterial(parameters), _h_order(getParam<MooseEnum>("h_order"))
 {
 }

Member Function Documentation

◆ computeQpProperties()

void SwitchingFunctionMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 35 of file SwitchingFunctionMaterial.C.

 {
   Real n = _eta[_qp];
   n = n > 1 ? 1 : (n < 0 ? 0 : n);
  
   switch (_h_order)
   {
     case 0: // SIMPLE
       _prop_f[_qp] = 3.0 * n * n - 2.0 * n * n * n;
       _prop_df[_qp] = 6.0 * n - 6.0 * n * n;
       _prop_d2f[_qp] = 6.0 - 12.0 * n;
       break;
  
     case 1: // HIGH
       _prop_f[_qp] = n * n * n * (6.0 * n * n - 15.0 * n + 10.0);
       _prop_df[_qp] = 30.0 * n * n * (n * n - 2.0 * n + 1.0);
       _prop_d2f[_qp] = n * (120.0 * n * n - 180.0 * n + 60.0);
       break;
  
     default:
       mooseError("Internal error");
   }
 }