ParallelElectricFieldInterface.C

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
 
#include "ParallelElectricFieldInterface.h"
 
registerMooseObject("ElectromagneticsApp", ParallelElectricFieldInterface);
 
InputParameters
ParallelElectricFieldInterface::validParams()
{
  InputParameters params = VectorInterfaceKernel::validParams();
  params.addClassDescription("VectorInterfaceKernel that implements the condition "
                             "$\\vec{E}_{1}^{\\parallel} - \\vec{E}_{2}^{\\parallel} = 0$");
  return params;
}
 
ParallelElectricFieldInterface::ParallelElectricFieldInterface(const InputParameters & parameters)
  : VectorInterfaceKernel(parameters)
{
}
 
Real
ParallelElectricFieldInterface::computeQpResidual(Moose::DGResidualType type)
{
  _u_parallel = -_normals[_qp].cross(_normals[_qp].cross(_u[_qp]));
  _secondary_parallel = -_normals[_qp].cross(_normals[_qp].cross(_neighbor_value[_qp]));
 
  Real res = 0;
 
  switch (type)
  {
    case Moose::Element:
      res = _test[_i][_qp] * (_u_parallel - _secondary_parallel);
      break;
 
    case Moose::Neighbor:
      res = _test_neighbor[_i][_qp] * -(_u_parallel - _secondary_parallel);
      break;
  }
 
  return res;
}
 
Real
ParallelElectricFieldInterface::computeQpJacobian(Moose::DGJacobianType type)
{
  _phi_u_parallel = -_normals[_qp].cross(_normals[_qp].cross(_phi[_j][_qp]));
  _phi_secondary_parallel = -_normals[_qp].cross(_normals[_qp].cross(_phi_neighbor[_j][_qp]));
 
  Real jac = 0;
 
  switch (type)
  {
    case Moose::ElementElement:
      jac = _test[_i][_qp] * _phi_u_parallel;
      break;
 
    case Moose::NeighborNeighbor:
      jac = _test_neighbor[_i][_qp] * _phi_secondary_parallel;
      break;
 
    case Moose::NeighborElement:
      jac = _test_neighbor[_i][_qp] * -_phi_u_parallel;
      break;
 
    case Moose::ElementNeighbor:
      jac = _test[_i][_qp] * -_phi_secondary_parallel;
      break;
  }
 
  return jac;
}