Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "MathCTDFreeEnergy.h"
      11             : 
      12             : registerMooseObject("PhaseFieldApp", MathCTDFreeEnergy);
      13             : registerMooseObject("PhaseFieldApp", ADMathCTDFreeEnergy);
      14             : 
      15             : template <bool is_ad>
      16             : InputParameters
      17          94 : MathCTDFreeEnergyTempl<is_ad>::validParams()
      18             : {
      19          94 :   InputParameters params = CompileTimeDerivativesMaterial<1, is_ad, 2>::validParams();
      20          94 :   params.addClassDescription("Material that implements the math free energy using the expression "
      21             :                              "builder and automatic differentiation");
      22         188 :   params.addRequiredCoupledVar("c", "Concentration variable");
      23          94 :   return params;
      24           0 : }
      25             : 
      26             : template <bool is_ad>
      27          72 : MathCTDFreeEnergyTempl<is_ad>::MathCTDFreeEnergyTempl(const InputParameters & parameters)
      28          72 :   : CompileTimeDerivativesMaterial<1, is_ad, 2>(parameters, {"c"})
      29             : {
      30         144 : }
      31             : 
      32             : template <bool is_ad>
      33             : void
      34     1130400 : MathCTDFreeEnergyTempl<is_ad>::computeQpProperties()
      35             : {
      36             :   const auto & [c] = _refs;
      37             :   const auto F = 1.0 / 4.0 * (1.0 + c) * (1.0 + c) * (1.0 - c) * (1.0 - c);
      38             : 
      39     1130400 :   evaluate(F);
      40     1130400 : }