#include <KKSPhaseConcentrationDerivatives.h>

Inheritance diagram for KKSPhaseConcentrationDerivatives:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	KKSPhaseConcentrationDerivatives (const InputParameters &parameters)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpProperties () override

Protected Attributes
const unsigned int	_num_c
	Number of global concentrations.

const std::vector< VariableName >	_c_names
	Names of global concentrations.

const VariableName	_eta_name
	Phase parameter.

std::vector< std::vector< std::vector< MaterialProperty< Real > * > > >	_dcidb
	Derivative of phase concentrations wrt global concentrations \( \frac d{db} c_i \).

std::vector< std::vector< MaterialProperty< Real > * > >	_dcideta
	Derivative of phase concentrations wrt eta \( \frac d{d{eta}} c_i \).

std::vector< std::vector< std::vector< const MaterialProperty< Real > * > > >	_d2Fidcidbi
	Second derivative of phase concentrations wrt two phase concentrations \( \frac {d^2}{dc_i db_i} F_i \).

const MaterialProperty< Real > &	_prop_h
	Switching function.

const MaterialProperty< Real > &	_prop_dh
	Derivative of switching function.


const std::vector< MaterialPropertyName >	_ci_names
	Phase concentrations.

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


const MaterialName	_Fa_name
	Free energy names.

const MaterialName	_Fb_name

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

Definition at line 15 of file KKSPhaseConcentrationDerivatives.h.

Constructor & Destructor Documentation

◆ KKSPhaseConcentrationDerivatives()

KKSPhaseConcentrationDerivatives::KKSPhaseConcentrationDerivatives ( const InputParameters & parameters )

Second derivative of free energy wrt phase concentrations for use in this material. In _d2Fidcidbi[m][n][l], m is phase index of Fi, n is the species index of ci, l is the species index of bi.

Definition at line 35 of file KKSPhaseConcentrationDerivatives.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _num_c(coupledComponents("global_cs")),
    _c_names(coupledNames("global_cs")),
    _eta_name(getVar("eta", 0)->name()),
    _ci_names(getParam<std::vector<MaterialPropertyName>>("ci_names")),
    _prop_ci(_num_c * 2),
    _dcidb(_num_c),
    _dcideta(_num_c),
    _Fa_name(getParam<MaterialName>("fa_name")),
    _Fb_name(getParam<MaterialName>("fb_name")),
    _d2Fidcidbi(2),
    _prop_h(getMaterialProperty<Real>("h_name")),
    _prop_dh(getMaterialPropertyDerivative<Real>("h_name", _eta_name))
{
  for (const auto m : make_range(_num_c * 2))
    _prop_ci[m] = &getMaterialPropertyByName<Real>(_ci_names[m]);
 
  for (const auto m : make_range(_num_c))
  {
    _dcideta[m].resize(2);
    _dcidb[m].resize(2);
    for (const auto n : make_range(2))
    {
      // Derivative of phase concentration wrt eta. In _dcideta[m][n], m is the species index of
      // ci, n is the phase index of ci
      _dcideta[m][n] = &declarePropertyDerivative<Real>(_ci_names[m * 2 + n], _eta_name);
      _dcidb[m][n].resize(_num_c);
 
      // Derivative of phase concentration wrt global concentration. In _dcidb[m][n][l], m is the
      //  species index of ci, n is the phase index of ci, l is the species index of b
      for (const auto l : make_range(_num_c))
        _dcidb[m][n][l] = &declarePropertyDerivative<Real>(_ci_names[m * 2 + n], _c_names[l]);
    }
  }
 
  for (const auto m : make_range(2))
  {
    _d2Fidcidbi[m].resize(_num_c);
    for (const auto n : make_range(_num_c))
    {
      _d2Fidcidbi[m][n].resize(_num_c);
      for (const auto l : make_range(_num_c))
      {
        if (m == 0)
          _d2Fidcidbi[0][n][l] =
              &getMaterialPropertyDerivative<Real>(_Fa_name, _ci_names[n * 2], _ci_names[l * 2]);
 
        else
          _d2Fidcidbi[1][n][l] = &getMaterialPropertyDerivative<Real>(
              _Fb_name, _ci_names[n * 2 + 1], _ci_names[l * 2 + 1]);
      }
    }
  }
}

Member Function Documentation

◆ computeQpProperties()

void KKSPhaseConcentrationDerivatives::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 97 of file KKSPhaseConcentrationDerivatives.C.

{
  // declare Jacobian matrix A
  Eigen::MatrixXd A(_num_c * 2, _num_c * 2);
 
  A.setZero();
 
  // fill in the non-zero elements in A
  for (const auto m : make_range(_num_c))
  {
    for (const auto n : make_range(_num_c))
    {
      // equal chemical potential derivative equations
      A(m * 2, n * 2) = (*_d2Fidcidbi[0][m][n])[_qp];
      A(m * 2, n * 2 + 1) = -(*_d2Fidcidbi[1][m][n])[_qp];
    }
 
    // concentration conservation derivative equations
    A(m * 2 + 1, m * 2) = 1 - _prop_h[_qp];
    A(m * 2 + 1, m * 2 + 1) = _prop_h[_qp];
  }
 
  A = A.inverse();
 
  // solve linear system of constraint derivatives wrt b for computing dcidb loop through
  // derivatives wrt the ith component; they have the same A, but different k_c
  for (const auto i : make_range(_num_c))
  {
    std::vector<Real> k_c(_num_c * 2);
    std::vector<Real> x_c(_num_c * 2);
 
    // assign the non-zero elements in k_c
    k_c[i * 2 + 1] = 1;
 
    // compute x_c
    for (const auto m : make_range(_num_c * 2))
    {
      for (const auto n : make_range(_num_c * 2))
        x_c[m] += A(m, n) * k_c[n];
    }
 
    // assign the values in x_c to _dcidb
    for (const auto m : make_range(_num_c))
    {
      for (const auto n : make_range(2))
        (*_dcidb[m][n][i])[_qp] = x_c[m * 2 + n];
    }
  }
 
  // solve linear system of constraint derivatives wrt eta for computing dcideta use the same
  // linear matrix as computing dcidb
  std::vector<Real> k_eta(_num_c * 2);
  std::vector<Real> x_eta(_num_c * 2);
 
  // fill in k_eta
  for (const auto m : make_range(_num_c))
  {
    k_eta[m * 2] = 0;
    k_eta[m * 2 + 1] = _prop_dh[_qp] * ((*_prop_ci[m * 2])[_qp] - (*_prop_ci[m * 2 + 1])[_qp]);
  }
 
  // compute x_eta
  for (const auto m : make_range(_num_c * 2))
  {
    for (const auto n : make_range(_num_c * 2))
      x_eta[m] += A(m, n) * k_eta[n];
  }
 
  // assign the values in x_eta to _dcideta
  for (const auto m : make_range(_num_c))
  {
    for (const auto n : make_range(2))
      (*_dcideta[m][n])[_qp] = x_eta[m * 2 + n];
  }
}

◆ validParams()

InputParameters KKSPhaseConcentrationDerivatives::validParams ( )

static

Definition at line 16 of file KKSPhaseConcentrationDerivatives.C.

{
  InputParameters params = DerivativeMaterialInterface<Material>::validParams();
  params.addClassDescription(
      "Computes the KKS phase concentration derivatives wrt global concentrations and order "
      "parameters, which are used in the chain rules in the KKS kernels. This class is intended to "
      "be used with KKSPhaseConcentrationMaterial.");
  params.addRequiredCoupledVar("global_cs", "The interpolated concentrations c, b, etc");
  params.addRequiredCoupledVar("eta", "Order parameter.");
  params.addRequiredParam<std::vector<MaterialPropertyName>>(
      "ci_names",
      "Phase concentrations. The order must match Fa, Fb, and global_cs, for example, c1, "
      "c2, b1, b2, etc");
  params.addRequiredParam<MaterialName>("fa_name", "Fa material object.");
  params.addRequiredParam<MaterialName>("fb_name", "Fb material object.");
  params.addParam<MaterialPropertyName>("h_name", "h", "Switching function h(eta).");
  return params;
}