www.mooseframework.org
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | List of all members
ACBarrierFunction Class Reference

#include <ACBarrierFunction.h>

Inheritance diagram for ACBarrierFunction:
[legend]

Public Member Functions

 ACBarrierFunction (const InputParameters &parameters)
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual Real computeDFDOP (PFFunctionType type)
 
virtual Real computeQpOffDiagJacobian (unsigned int jvar)
 

Protected Attributes

unsigned int _n_eta
 
const MaterialProperty< Real > & _d2mudvar2
 
const std::vector< VariableName > _vname
 
std::vector< const MaterialProperty< Real > * > _d2mudvardeta
 
const JvarMap & _vmap
 
const NonlinearVariableName _uname
 
const MaterialPropertyName _gamma_name
 
const GenericMaterialProperty< Real, is_ad > & _gamma
 
const GenericMaterialProperty< Real, is_ad > & _dmudvar
 

Private Member Functions

Real calculateF0 ()
 

Detailed Description

Definition at line 43 of file ACBarrierFunction.h.

Constructor & Destructor Documentation

◆ ACBarrierFunction()

ACBarrierFunction::ACBarrierFunction ( const InputParameters parameters)

Definition at line 38 of file ACBarrierFunction.C.

39  : ACBarrierFunctionTempl<false>(parameters),
40  _n_eta(_vals.size()),
41  _d2mudvar2(getMaterialPropertyDerivative<Real>("mu", _uname, _uname)),
42  _vname(getParam<std::vector<VariableName>>("v")),
44  _vmap(getParameterJvarMap("v"))
45 {
46  for (unsigned int i = 0; i < _n_eta; ++i)
47  _d2mudvardeta[i] = &getMaterialPropertyDerivative<Real>("mu", _uname, _vname[i]);
48 }
const std::vector< VariableName > _vname
const NonlinearVariableName _uname
const JvarMap & _vmap
std::vector< const MaterialProperty< Real > * > _d2mudvardeta
const MaterialProperty< Real > & _d2mudvar2

Member Function Documentation

◆ calculateF0()

Real ACBarrierFunction::calculateF0 ( )
private

Definition at line 118 of file ACBarrierFunction.C.

Referenced by computeDFDOP(), and computeQpOffDiagJacobian().

119 {
120  Real var_phase = 0.25 * _u[_qp] * _u[_qp] * _u[_qp] * _u[_qp] - 0.5 * _u[_qp] * _u[_qp];
121  Real eta_phase = 0.0;
122  Real eta_interface = 0.0;
123  Real sum_eta_j2;
124 
125  for (unsigned int i = 0; i < _n_eta; ++i)
126  {
127  sum_eta_j2 = 0.0;
128  for (unsigned int j = i + 1; j < _n_eta; ++j)
129  sum_eta_j2 += (*_vals[j])[_qp] * (*_vals[j])[_qp];
130  // eta_phase += 0.25 * (*_vals[i])[_qp] * (*_vals[i])[_qp] * (*_vals[i])[_qp] * (*_vals[i])[_qp]
131  // -
132  // 0.5 * (*_vals[i])[_qp] * (*_vals[i])[_qp];
133  eta_phase += 0.25 * Utility::pow<4>((*_vals[i])[_qp]) - 0.5 * Utility::pow<2>((*_vals[i])[_qp]);
134  eta_interface +=
135  (_u[_qp] * _u[_qp] + sum_eta_j2) * (*_vals[i])[_qp] * (*_vals[i])[_qp] * _gamma[_qp];
136  }
137  return 0.25 + var_phase + eta_phase + eta_interface;
138 }
const GenericMaterialProperty< Real, is_ad > & _gamma
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")

◆ computeDFDOP()

Real ACBarrierFunction::computeDFDOP ( PFFunctionType  type)
protectedvirtual

Definition at line 51 of file ACBarrierFunction.C.

52 {
53  switch (type)
54  {
55  case Residual:
56  return _dmudvar[_qp] * calculateF0();
57 
58  case Jacobian:
59  {
60  Real df0dvar = 0.0;
61  for (unsigned int i = 0; i < _n_eta; ++i)
62  df0dvar += (*_vals[i])[_qp] * (*_vals[i])[_qp];
63 
64  df0dvar *= 2.0 * _gamma[_qp] * _u[_qp];
65  df0dvar += _u[_qp] * _u[_qp] * _u[_qp] - _u[_qp];
66 
67  return (_d2mudvar2[_qp] * calculateF0() + _dmudvar[_qp] * df0dvar) * _phi[_j][_qp];
68  }
69 
70  default:
71  mooseError("Invalid type passed in");
72  }
73 }
void mooseError(Args &&... args)
const GenericMaterialProperty< Real, is_ad > & _gamma
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< Real > & _d2mudvar2
const GenericMaterialProperty< Real, is_ad > & _dmudvar

◆ computeQpOffDiagJacobian()

Real ACBarrierFunction::computeQpOffDiagJacobian ( unsigned int  jvar)
protectedvirtual

Definition at line 93 of file ACBarrierFunction.C.

94 {
95  const unsigned int j = mapJvarToCvar(jvar);
96  Real sum_etai2 = 0.0;
97  Real df0deta_base = 0.0;
98  Real df0deta = 0.0;
99 
100  for (unsigned int i = 0; i < _n_eta; ++i)
101  if (i != j)
102  sum_etai2 += (*_vals[i])[_qp] * (*_vals[i])[_qp];
103 
104  auto etavar = mapJvarToCvar(jvar, _vmap);
105  if (etavar >= 0)
106  {
107  df0deta_base = (*_vals[etavar])[_qp] * (*_vals[etavar])[_qp] - 1.0 +
108  2.0 * _gamma[_qp] * (_u[_qp] * _u[_qp] + sum_etai2);
109  df0deta = (*_vals[etavar])[_qp] * df0deta_base;
110 
111  return ((*_d2mudvardeta[etavar])[_qp] * calculateF0() + _dmudvar[_qp] * df0deta) *
112  _phi[_j][_qp] * _test[_i][_qp] * _L[_qp];
113  }
114  return 0.0;
115 }
const GenericMaterialProperty< Real, is_ad > & _gamma
const JvarMap & _vmap
std::vector< const MaterialProperty< Real > * > _d2mudvardeta
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
const GenericMaterialProperty< Real, is_ad > & _dmudvar

◆ validParams()

InputParameters ACBarrierFunctionTempl< is_ad >::validParams ( )
staticinherited

Definition at line 18 of file ACBarrierFunction.C.

19 {
21  params.addRequiredParam<MaterialPropertyName>(
22  "gamma", "The interface profile coefficient to use with the kernel");
23  params.addClassDescription("Allen-Cahn kernel used when 'mu' is a function of variables");
24  return params;
25 }
void addRequiredParam(const std::string &name, const std::string &doc_string)
InputParameters validParams()
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _d2mudvar2

const MaterialProperty<Real>& ACBarrierFunction::_d2mudvar2
protected

Definition at line 53 of file ACBarrierFunction.h.

Referenced by computeDFDOP().

◆ _d2mudvardeta

std::vector<const MaterialProperty<Real> *> ACBarrierFunction::_d2mudvardeta
protected

Definition at line 56 of file ACBarrierFunction.h.

Referenced by ACBarrierFunction(), and computeQpOffDiagJacobian().

◆ _dmudvar

const GenericMaterialProperty<Real, is_ad>& ACBarrierFunctionTempl< is_ad >::_dmudvar
protectedinherited

Definition at line 35 of file ACBarrierFunction.h.

Referenced by computeDFDOP(), and computeQpOffDiagJacobian().

◆ _gamma

const GenericMaterialProperty<Real, is_ad>& ACBarrierFunctionTempl< is_ad >::_gamma
protectedinherited

Definition at line 34 of file ACBarrierFunction.h.

Referenced by calculateF0(), computeDFDOP(), and computeQpOffDiagJacobian().

◆ _gamma_name

const MaterialPropertyName ACBarrierFunctionTempl< is_ad >::_gamma_name
protectedinherited

Definition at line 33 of file ACBarrierFunction.h.

◆ _n_eta

unsigned int ACBarrierFunction::_n_eta
protected

◆ _uname

const NonlinearVariableName ACBarrierFunctionTempl< is_ad >::_uname
protectedinherited

Definition at line 32 of file ACBarrierFunction.h.

Referenced by ACBarrierFunction().

◆ _vmap

const JvarMap& ACBarrierFunction::_vmap
protected

Definition at line 57 of file ACBarrierFunction.h.

Referenced by computeQpOffDiagJacobian().

◆ _vname

const std::vector<VariableName> ACBarrierFunction::_vname
protected

Definition at line 55 of file ACBarrierFunction.h.

Referenced by ACBarrierFunction().


The documentation for this class was generated from the following files: