PenaltyPeriodicSegmentalConstraint.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 "MortarScalarBase.h"
#include "DerivativeMaterialInterface.h"
#include "MooseVariableScalar.h"
#include "Assembly.h"

class PenaltyPeriodicSegmentalConstraint : public DerivativeMaterialInterface<MortarScalarBase>
{
public:
  static InputParameters validParams();
  PenaltyPeriodicSegmentalConstraint(const InputParameters & parameters);

protected:
  virtual void precalculateResidual() override;
  virtual void precalculateJacobian() override;
  virtual void initScalarQpResidual() override;

  virtual Real computeQpResidual(Moose::MortarType mortar_type) override;
  Real computeQpJacobian(Moose::ConstraintJacobianType /*jacobian_type*/,
                         unsigned int /*jvar*/) override
  {
    return 0;
  };

  virtual Real computeScalarQpResidual() override;

  virtual Real computeScalarQpJacobian() override;

  // using MortarScalarBase::computeOffDiagJacobianScalar;

  virtual Real computeQpOffDiagJacobianScalar(Moose::MortarType mortar_type,
                                              const unsigned int jvar) override;

  virtual Real computeScalarQpOffDiagJacobian(Moose::MortarType mortar_type,
                                              const unsigned int jvar) override;

  // Compute concentration jump
  void precalculateMaterial();
  // Compute penalty parameter
  void precalculateStability();

protected:
  Real _temp_jump_global;

  Real _tau_s;

  const MooseVariableScalar * const _kappa_aux_ptr;

  const unsigned int _ka_order;

  const VariableValue & _kappa_aux;

  const Real _pen_scale;
};