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

#include <MixedSwitchingFunctionMaterial.h>

Inheritance diagram for MixedSwitchingFunctionMaterial:

Public Member Functions
	MixedSwitchingFunctionMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties ()

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

Real	_weight
	Weight parameter of mixed-type 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 MixedSwitchingFunctionMaterial.h.

Constructor & Destructor Documentation

◆ MixedSwitchingFunctionMaterial()

MixedSwitchingFunctionMaterial::MixedSwitchingFunctionMaterial ( const InputParameters & parameters )

Definition at line 32 of file MixedSwitchingFunctionMaterial.C.

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

Member Function Documentation

◆ computeQpProperties()

void MixedSwitchingFunctionMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 40 of file MixedSwitchingFunctionMaterial.C.

 {
   Real n = _eta[_qp];
   n = n > 1 ? 1 : (n < 0 ? 0 : n);
  
   switch (_h_order)
   {
     case 0: // MIX234
       _prop_f[_qp] =
           n * n * (3.0 * _weight + (4.0 - 6.0 * _weight) * n + (3.0 * _weight - 3.0) * n * n);
       _prop_df[_qp] =
           n * (6.0 * _weight + (12.0 - 18.0 * _weight) * n + (12.0 * _weight - 12.0) * n * n);
       _prop_d2f[_qp] = 6.0 * _weight + ((24.0 - 36.0 * _weight) + (36.0 * _weight - 36.0) * n) * n;
       break;
  
     case 1: // MIX246
       _prop_f[_qp] =
           n * n * (2.0 * _weight + n * n * ((3.0 - 4.0 * _weight) + (2.0 * _weight - 2.0) * n * n));
       _prop_df[_qp] =
           n * (4.0 * _weight + n * n * ((12.0 - 16.0 * _weight) + (12.0 * _weight - 12.0) * n * n));
       _prop_d2f[_qp] =
           4.0 * _weight + n * n * ((36.0 - 48.0 * _weight) + (60.0 * _weight - 60.0) * n * n);
       break;
  
     default:
       mooseError("Internal error");
   }
 }