CoupledSwitchingTimeDerivative.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 "CoupledTimeDerivative.h"
#include "ADCoupledTimeDerivative.h"
#include "JvarMapInterface.h"
#include "DerivativeMaterialInterface.h"

// Forward Declaration

template <bool is_ad>
using CoupledSwitchingTimeDerivativeBase =
    typename std::conditional<is_ad, ADCoupledTimeDerivative, CoupledTimeDerivative>::type;

template <bool is_ad>
class CoupledSwitchingTimeDerivativeTempl
  : public DerivativeMaterialInterface<
        JvarMapKernelInterface<CoupledSwitchingTimeDerivativeBase<is_ad>>>
{
public:
  static InputParameters validParams();

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

protected:
  const VariableName _v_name;

  std::vector<MaterialPropertyName> _Fj_names;

  const 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 CoupledSwitchingTimeDerivativeBase<is_ad>::_qp;
};

class CoupledSwitchingTimeDerivative : public CoupledSwitchingTimeDerivativeTempl<false>
{
public:
  CoupledSwitchingTimeDerivative(const InputParameters & parameters);
  virtual void initialSetup() override;

protected:
  virtual Real computeQpResidual() override;
  virtual Real computeQpJacobian() override;
  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
  std::vector<const MaterialProperty<Real> *> _prop_dFjdv;

  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 ADCoupledSwitchingTimeDerivative : public CoupledSwitchingTimeDerivativeTempl<true>
{
public:
  using CoupledSwitchingTimeDerivativeTempl<true>::CoupledSwitchingTimeDerivativeTempl;

protected:
  virtual ADReal precomputeQpResidual() override;
};