Enforce gradient continuity between two different variables across a subdomain boundary. More...

#include <InterfaceDiffusionFluxMatch.h>

Inheritance diagram for InterfaceDiffusionFluxMatch:

Public Member Functions
	InterfaceDiffusionFluxMatch (const InputParameters &parameters)

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

virtual Real	computeQpJacobian (Moose::DGJacobianType type) override

Protected Attributes
const Real	_D
	diffusion coefficient More...

const Real	_D_neighbor
	neighbor diffusion coefficient More...

Detailed Description

Enforce gradient continuity between two different variables across a subdomain boundary.

Definition at line 23 of file InterfaceDiffusionFluxMatch.h.

Constructor & Destructor Documentation

◆ InterfaceDiffusionFluxMatch()

InterfaceDiffusionFluxMatch::InterfaceDiffusionFluxMatch ( const InputParameters & parameters )

Definition at line 24 of file InterfaceDiffusionFluxMatch.C.

   : InterfaceDiffusionBase(parameters)
 {
 }

Member Function Documentation

◆ computeQpJacobian()

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

overrideprotectedvirtual

Definition at line 49 of file InterfaceDiffusionFluxMatch.C.

 {
   switch (type)
   {
     case Moose::ElementElement:
       return _D * _grad_phi[_j][_qp] * _normals[_qp] * _test[_i][_qp];
  
     case Moose::NeighborNeighbor:
       return -_D_neighbor * _grad_phi_neighbor[_j][_qp] * _normals[_qp] * _test_neighbor[_i][_qp];
  
     case Moose::ElementNeighbor:
       return -_D_neighbor * _grad_phi_neighbor[_j][_qp] * _normals[_qp] * _test[_i][_qp];
  
     case Moose::NeighborElement:
       return _D * _grad_phi[_j][_qp] * _normals[_qp] * _test_neighbor[_i][_qp];
   }
   mooseError("Internal error.");
 }

◆ computeQpResidual()

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

overrideprotectedvirtual

Definition at line 30 of file InterfaceDiffusionFluxMatch.C.

 {
   // equal gradients means difference is zero
   Real res =
       _D * _grad_u[_qp] * _normals[_qp] - _D_neighbor * _grad_neighbor_value[_qp] * _normals[_qp];
  
   switch (type)
   {
     case Moose::Element:
       return res * _test[_i][_qp];
  
     case Moose::Neighbor:
       return res * _test_neighbor[_i][_qp];
   }
  
   mooseError("Internal error.");
 }

Member Data Documentation

◆ _D

const Real InterfaceDiffusionBase::_D

protectedinherited

diffusion coefficient

Definition at line 30 of file InterfaceDiffusionBase.h.

Referenced by InterfaceDiffusionBoundaryTerm::computeQpJacobian(), computeQpJacobian(), computeQpResidual(), and InterfaceDiffusionBoundaryTerm::computeQpResidual().

◆ _D_neighbor

const Real InterfaceDiffusionBase::_D_neighbor

protectedinherited

neighbor diffusion coefficient

Definition at line 33 of file InterfaceDiffusionBase.h.