PFCRFFEnergyDensity Class Reference

#include <PFCRFFEnergyDensity.h>

Inheritance diagram for PFCRFFEnergyDensity:

Public Member Functions
	PFCRFFEnergyDensity (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue ()

Protected Attributes
unsigned int	_order
std::vector< const VariableValue * >	_vals
Real	_a
Real	_b
Real	_c
unsigned int	_num_exp_terms
MooseEnum	_log_approach
Real	_tol

Detailed Description

Definition at line 20 of file PFCRFFEnergyDensity.h.

Constructor & Destructor Documentation

PFCRFFEnergyDensity()

PFCRFFEnergyDensity::PFCRFFEnergyDensity

(

const InputParameters &

parameters

)

Definition at line 33 of file PFCRFFEnergyDensity.C.

  : AuxKernel(parameters),
    _order(coupledComponents("v")),
    _a(getParam<Real>("a")),
    _b(getParam<Real>("b")),
    _c(getParam<Real>("c")),
    _num_exp_terms(getParam<unsigned int>("num_exp_terms")),
    _log_approach(getParam<MooseEnum>("log_approach")),
    _tol(getParam<Real>("tol"))
{
  _vals.resize(_order);
  for (unsigned int i = 0; i < _order; ++i)
    _vals[i] = &coupledValue("v", i);
}

Member Function Documentation

computeValue()

Real PFCRFFEnergyDensity::computeValue ( )

protectedvirtual

Definition at line 49 of file PFCRFFEnergyDensity.C.

{
  Real val = 0.0;
  switch (_log_approach)
  {
    case 0: // approach using tolerence
      if (1.0 + (*_vals[0])[_qp] < _tol)
        val += ((1.0 + _tol) * std::log(1 + _tol)) - _tol;
      else
        val += ((1.0 + (*_vals[0])[_qp]) * std::log(1 + (*_vals[0])[_qp])) - (*_vals[0])[_qp];
      break;
 
    case 1: // approach using cancellation
      val += ((1.0 + (*_vals[0])[_qp]) * std::log(1.0 + (*_vals[0])[_qp])) - (*_vals[0])[_qp];
      break;
 
    case 2: // approach using Taylor Series Expansion
      Real coef = 1.0;
 
      for (unsigned int i = 2; i < (2 + _num_exp_terms); i++)
      {
        if (i == 2)
          coef = _c;
        else if (i == 3)
          coef = _a;
        else if (i == 4)
          coef = _b;
        else
          coef = 1.0;
 
        val +=
            coef * (std::pow(-1.0, Real(i)) / (i * (i - 1))) * std::pow((*_vals[0])[_qp], Real(i));
      }
      break;
  }
 
  // Loop Through Variables
  Real sumL = 0.0;
  for (unsigned int i = 1; i < _order; ++i)
    sumL += (*_vals[i])[_qp] * 0.5;
 
  val -= ((*_vals[0])[_qp] * sumL);
 
  return val;
}

Member Data Documentation

_a

Real PFCRFFEnergyDensity::_a

protected

Definition at line 31 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_b

Real PFCRFFEnergyDensity::_b

protected

Definition at line 32 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_c

Real PFCRFFEnergyDensity::_c

protected

Definition at line 33 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_log_approach

MooseEnum PFCRFFEnergyDensity::_log_approach

protected

Definition at line 35 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_num_exp_terms

unsigned int PFCRFFEnergyDensity::_num_exp_terms

protected

Definition at line 34 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_order

unsigned int PFCRFFEnergyDensity::_order

protected

Definition at line 28 of file PFCRFFEnergyDensity.h.

Referenced by computeValue(), and PFCRFFEnergyDensity().

_tol

Real PFCRFFEnergyDensity::_tol

protected

Definition at line 36 of file PFCRFFEnergyDensity.h.

Referenced by computeValue().

_vals

std::vector<const VariableValue *> PFCRFFEnergyDensity::_vals

protected

Definition at line 29 of file PFCRFFEnergyDensity.h.

Referenced by computeValue(), and PFCRFFEnergyDensity().

---

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

• phase_field/include/auxkernels/PFCRFFEnergyDensity.h
• phase_field/src/auxkernels/PFCRFFEnergyDensity.C