#include <KKSPhaseConcentrationMultiPhaseDerivatives.h>

Inheritance diagram for KKSPhaseConcentrationMultiPhaseDerivatives:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	KKSPhaseConcentrationMultiPhaseDerivatives (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 std::vector< VariableName >	_eta_names
	Phase parameters.

const unsigned int	_num_j
	Number of phase parameters.

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

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

std::vector< MaterialName >	_Fj_names
	Free energy names.

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 \).

std::vector< std::vector< const MaterialProperty< Real > * > >	_dhjdetai
	Derivatives of switching functions.


std::vector< const MaterialProperty< Real > * >	_prop_ci
	Phase concentrations.

std::vector< MaterialPropertyName >	_ci_names


std::vector< MaterialPropertyName >	_hj_names
	Switching functions.

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

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 14 of file KKSPhaseConcentrationMultiPhaseDerivatives.h.

Constructor & Destructor Documentation

◆ KKSPhaseConcentrationMultiPhaseDerivatives()

KKSPhaseConcentrationMultiPhaseDerivatives::KKSPhaseConcentrationMultiPhaseDerivatives ( const InputParameters & parameters )

Definition at line 36 of file KKSPhaseConcentrationMultiPhaseDerivatives.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _num_c(coupledComponents("global_cs")),
    _c_names(coupledNames("global_cs")),
    _eta_names(coupledNames("all_etas")),
    _num_j(coupledComponents("all_etas")),
    _prop_ci(_num_c * _num_j),
    _ci_names(getParam<std::vector<MaterialPropertyName>>("ci_names")),
    _dcidetaj(_num_c),
    _dcidb(_num_c),
    _Fj_names(getParam<std::vector<MaterialName>>("Fj_names")),
    _d2Fidcidbi(_num_j),
    _hj_names(getParam<std::vector<MaterialPropertyName>>("hj_names")),
    _prop_hj(_num_j),
    _dhjdetai(_num_j)
 
{
  for (const auto m : make_range(_num_c * _num_j))
    _prop_ci[m] = &getMaterialPropertyByName<Real>(_ci_names[m]);
 
  for (const auto m : make_range(_num_c))
  {
    _dcidb[m].resize(_num_j);
    _dcidetaj[m].resize(_num_j);
 
    for (const auto n : make_range(_num_j))
    {
      _dcidb[m][n].resize(_num_c);
      _dcidetaj[m][n].resize(_num_j);
 
      // 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, and l is the species index of b
      for (const auto l : make_range(_num_c))
        _dcidb[m][n][l] = &declarePropertyDerivative<Real>(_ci_names[n + m * _num_j], _c_names[l]);
 
      // Derivative of phase concentration wrt eta. In _dcidetaj[m][n][l], m is the species index
      // of ci, n is the phase index of ci, and l is the phase of etaj
      for (const auto l : make_range(_num_j))
        _dcidetaj[m][n][l] =
            &declarePropertyDerivative<Real>(_ci_names[n + m * _num_j], _eta_names[l]);
    }
  }
 
  // 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.
  for (const auto m : make_range(_num_j))
  {
    _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))
        _d2Fidcidbi[m][n][l] = &getMaterialPropertyDerivative<Real>(
            _Fj_names[m], _ci_names[m + n * _num_j], _ci_names[m + l * _num_j]);
    }
  }
 
  for (const auto m : make_range(_num_j))
  {
    _prop_hj[m] = &getMaterialPropertyByName<Real>(_hj_names[m]);
 
    _dhjdetai[m].resize(_num_j);
 
    for (const auto n : make_range(_num_j))
      _dhjdetai[m][n] = &getMaterialPropertyDerivative<Real>(_hj_names[m], _eta_names[n]);
  }
}

Member Function Documentation

◆ computeQpProperties()

void KKSPhaseConcentrationMultiPhaseDerivatives::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 109 of file KKSPhaseConcentrationMultiPhaseDerivatives.C.

{
  // declare Jacobian matrix A
  Eigen::MatrixXd A(_num_c * _num_j, _num_c * _num_j);
 
  // initialize all elements in A to be zero
  A.setZero();
 
  // fill in the non-zero elements in A
  for (const auto m : make_range(_num_c))
  {
    // equal chemical potential derivative equations
    for (const auto n : make_range(_num_j - 1))
    {
      for (const auto l : make_range(_num_c))
      {
        A(m * _num_j + n, n + l * _num_j) = (*_d2Fidcidbi[n][m][l])[_qp];
        A(m * _num_j + n, n + l * _num_j + 1) = -(*_d2Fidcidbi[n + 1][m][l])[_qp];
      }
    }
 
    // concentration conservation derivative equations
    for (const auto n : make_range(_num_j))
      A((m + 1) * _num_j - 1, m * _num_j + n) = (*_prop_hj[n])[_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_j * _num_c);
    std::vector<Real> x_c(_num_j * _num_c);
 
    // assign non-zero elements in k_c
    k_c[i * _num_j + _num_j - 1] = 1;
 
    // compute x_c
    for (const auto m : make_range(_num_j * _num_c))
    {
      for (const auto n : make_range(_num_j * _num_c))
        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(_num_j))
        (*_dcidb[m][n][i])[_qp] = x_c[m * _num_j + n];
    }
  }
 
  // solve linear system of constraint derivatives wrt eta for computing dcidetaj use the same
  // linear matrix as computing dcidb
  for (const auto i : make_range(_num_j))
  {
    std::vector<Real> k_eta(_num_j * _num_c);
    std::vector<Real> x_eta(_num_j * _num_c);
 
    // assign non-zero elements in k_eta
    for (const auto m : make_range(_num_c))
    {
      Real sum = 0.0;
 
      for (const auto n : make_range(_num_j))
        sum += (*_dhjdetai[n][i])[_qp] * (*_prop_ci[m * _num_j + n])[_qp];
 
      k_eta[m * _num_j + _num_j - 1] = -sum;
    }
 
    // compute x_eta
    for (const auto m : make_range(_num_j * _num_c))
    {
      for (const auto n : make_range(_num_j * _num_c))
        x_eta[m] += A(m, n) * k_eta[n];
    }
 
    // assign the values in x_eta to _dcidetaj
    for (const auto m : make_range(_num_c))
    {
      for (const auto n : make_range(_num_j))
        (*_dcidetaj[m][n][i])[_qp] = x_eta[m * _num_j + n];
    }
  }
}

◆ validParams()

InputParameters KKSPhaseConcentrationMultiPhaseDerivatives::validParams ( )

static

Definition at line 16 of file KKSPhaseConcentrationMultiPhaseDerivatives.C.

{
  InputParameters params = DerivativeMaterialInterface<Material>::validParams();
  params.addClassDescription(
      "Computes the KKS phase concentration derivatives wrt global concentrations and order "
      "parameters, which are used for the chain rule in the KKS kernels. This class is intended to "
      "be used with KKSPhaseConcentrationMultiPhaseMaterial.");
  params.addRequiredCoupledVar("global_cs", "The interpolated concentrations c, b, etc");
  params.addRequiredCoupledVar("all_etas", "Order parameters.");
  params.addRequiredParam<std::vector<MaterialPropertyName>>(
      "ci_names",
      "Phase concentrations. They must have the same order as Fj_names and global_cs, for "
      "example, c1, c2, b1, b2.");
  params.addRequiredParam<std::vector<MaterialName>>(
      "Fj_names", "Free energy material objects in the same order as all_etas.");
  params.addRequiredParam<std::vector<MaterialPropertyName>>(
      "hj_names", "witching functions in the same order as all_etas.");
  return params;
}