ACSwitching.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ACBulk.h"
#include "ADAllenCahnBase.h"

template <bool is_ad>
using ACSwitchingBase = typename std::conditional<is_ad, ADAllenCahnBase<Real>, ACBulk<Real>>::type;

template <bool is_ad>
class ACSwitchingTempl : public ACSwitchingBase<is_ad>
{
public:
  static InputParameters validParams();

  ACSwitchingTempl(const InputParameters & parameters);
  virtual void initialSetup() override;

protected:
  VariableName _etai_name;

  std::vector<MaterialPropertyName> _Fj_names;
  unsigned int _num_j;

  std::vector<const GenericMaterialProperty<Real, is_ad> *> _prop_Fj;

  std::vector<MaterialPropertyName> _hj_names;

  std::vector<const GenericMaterialProperty<Real, is_ad> *> _prop_dhjdetai;

  using ACSwitchingBase<is_ad>::_qp;
  using ACSwitchingBase<is_ad>::_var;
  using ACSwitchingBase<is_ad>::_test;
};

class ACSwitching : public ACSwitchingTempl<false>
{
public:
  ACSwitching(const InputParameters & parameters);

  virtual void initialSetup() override;

protected:
  virtual Real computeDFDOP(PFFunctionType type) override;
  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;

  std::vector<std::vector<const MaterialProperty<Real> *>> _prop_dFjdarg;

  std::vector<const MaterialProperty<Real> *> _prop_d2hjdetai2;

  std::vector<std::vector<const MaterialProperty<Real> *>> _prop_d2hjdetaidarg;
};

class ADACSwitching : public ACSwitchingTempl<true>
{
public:
  using ACSwitchingTempl<true>::ACSwitchingTempl;

protected:
  virtual ADReal computeDFDOP() override;
};