InterfaceKernel to enforce a Lagrange-Multiplier based componentwise continuity of a variable gradient. More...

#include <EqualGradientLagrangeInterface.h>

Inheritance diagram for EqualGradientLagrangeInterface:

Public Member Functions
	EqualGradientLagrangeInterface (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual (Moose::DGResidualType type) override

virtual Real	computeQpJacobian (Moose::DGJacobianType type) override

virtual Real	computeQpOffDiagJacobian (Moose::DGJacobianType type, unsigned int jvar) override

Protected Attributes
const unsigned int	_component

const VariableValue &	_lambda
	Lagrange multiplier. More...

const unsigned int	_lambda_jvar

Detailed Description

InterfaceKernel to enforce a Lagrange-Multiplier based componentwise continuity of a variable gradient.

Definition at line 23 of file EqualGradientLagrangeInterface.h.

Constructor & Destructor Documentation

◆ EqualGradientLagrangeInterface()

EqualGradientLagrangeInterface::EqualGradientLagrangeInterface ( const InputParameters & parameters )

Definition at line 30 of file EqualGradientLagrangeInterface.C.

   : InterfaceKernel(parameters),
     _component(getParam<unsigned int>("component")),
     _lambda(getVar("lambda", 0)->sln()),
     _lambda_jvar(getVar("lambda", 0)->number())
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real EqualGradientLagrangeInterface::computeQpJacobian ( Moose::DGJacobianType type )

overrideprotectedvirtual

Definition at line 53 of file EqualGradientLagrangeInterface.C.

 {
   return 0.0;
 }

◆ computeQpOffDiagJacobian()

Real EqualGradientLagrangeInterface::computeQpOffDiagJacobian	(	Moose::DGJacobianType	type,
		unsigned int	jvar
	)

overrideprotectedvirtual

Definition at line 59 of file EqualGradientLagrangeInterface.C.

 {
   if (jvar != _lambda_jvar)
     return 0.0;
  
   // lambda is only solved on the element side
   switch (type)
   {
     case Moose::ElementElement:
       return _phi[_j][_qp] * _grad_test[_i][_qp](_component);
  
     case Moose::NeighborElement:
       return -_phi[_j][_qp] * _grad_test_neighbor[_i][_qp](_component);
  
     default:
       return 0.0;
   }
 }

◆ computeQpResidual()

Real EqualGradientLagrangeInterface::computeQpResidual ( Moose::DGResidualType type )

overrideprotectedvirtual

Definition at line 39 of file EqualGradientLagrangeInterface.C.

 {
   switch (type)
   {
     case Moose::Element:
       return _lambda[_qp] * _grad_test[_i][_qp](_component);
  
     case Moose::Neighbor:
       return -_lambda[_qp] * _grad_test_neighbor[_i][_qp](_component);
   }
  
   mooseError("Internal error.");
 }