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;
};