Lagrange multiplier "FaceKernel" that is used in conjunction with EqualGradientLagrangeInterface. More...

#include <EqualGradientLagrangeMultiplier.h>

Inheritance diagram for EqualGradientLagrangeMultiplier:

Public Member Functions
	EqualGradientLagrangeMultiplier (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
	x,y,z component of the gradient to constrain More...

unsigned int	_neighbor_jvar
	variable to control gradient on the slave side of the interface More...

const Real	_jacobian_fill
	compensate Jacobian fill term from NullKernel More...


const VariableGradient &	_grad_element_value
	variable to control gradient on the master side of the interface More...

unsigned int	_element_jvar

Detailed Description

Lagrange multiplier "FaceKernel" that is used in conjunction with EqualGradientLagrangeInterface.

Definition at line 23 of file EqualGradientLagrangeMultiplier.h.

Constructor & Destructor Documentation

◆ EqualGradientLagrangeMultiplier()

EqualGradientLagrangeMultiplier::EqualGradientLagrangeMultiplier ( const InputParameters & parameters )

Definition at line 34 of file EqualGradientLagrangeMultiplier.C.

   : InterfaceKernel(parameters),
     _component(getParam<unsigned int>("component")),
     _grad_element_value(getVar("element_var", 0)->gradSln()),
     _element_jvar(getVar("element_var", 0)->number()),
     _neighbor_jvar(_neighbor_var.number()),
     _jacobian_fill(getParam<Real>("jacobian_fill"))
 {
 }

Member Function Documentation

◆ computeQpJacobian()

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

overrideprotectedvirtual

Definition at line 55 of file EqualGradientLagrangeMultiplier.C.

 {
   if (type == Moose::ElementNeighbor)
     return -_grad_phi_neighbor[_j][_qp](_component) * _test[_i][_qp];
  
   if (type == Moose::ElementElement)
     return -_jacobian_fill;
  
   return 0.0;
 }

◆ computeQpOffDiagJacobian()

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

overrideprotectedvirtual

Definition at line 67 of file EqualGradientLagrangeMultiplier.C.

 {
   if (type == Moose::ElementElement && jvar == _element_jvar)
     return _grad_phi[_j][_qp](_component) * _test[_i][_qp];
  
   return 0.0;
 }

◆ computeQpResidual()

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

overrideprotectedvirtual

Definition at line 45 of file EqualGradientLagrangeMultiplier.C.

 {
   if (type == Moose::Element)
     return (_grad_element_value[_qp](_component) - _grad_neighbor_value[_qp](_component)) *
            _test[_i][_qp];
  
   return 0.0;
 }