DGConvection.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 "DGConvection.h"
 
registerMooseObject("MooseApp", DGConvection);
 
InputParameters
DGConvection::validParams()
{
  InputParameters params = DGKernel::validParams();
  params.addRequiredParam<RealVectorValue>("velocity", "Velocity vector");
  params.addClassDescription("DG upwinding for the convection");
  return params;
}
 
DGConvection::DGConvection(const InputParameters & parameters)
  : DGKernel(parameters), _velocity(getParam<RealVectorValue>("velocity"))
{
}
 
Real
DGConvection::computeQpResidual(Moose::DGResidualType type)
{
  Real r = 0;
 
  Real vdotn = _velocity * _normals[_qp];
 
  switch (type)
  {
    case Moose::Element:
      if (vdotn >= 0)
        r += vdotn * _u[_qp] * _test[_i][_qp];
      else
        r += vdotn * _u_neighbor[_qp] * _test[_i][_qp];
      break;
 
    case Moose::Neighbor:
      if (vdotn >= 0)
        r -= vdotn * _u[_qp] * _test_neighbor[_i][_qp];
      else
        r -= vdotn * _u_neighbor[_qp] * _test_neighbor[_i][_qp];
      break;
  }
 
  return r;
}
 
Real
DGConvection::computeQpJacobian(Moose::DGJacobianType type)
{
  Real r = 0;
 
  Real vdotn = _velocity * _normals[_qp];
 
  switch (type)
  {
    case Moose::ElementElement:
      if (vdotn >= 0)
        r += vdotn * _phi[_j][_qp] * _test[_i][_qp];
      break;
 
    case Moose::ElementNeighbor:
      if (vdotn < 0)
        r += vdotn * _phi_neighbor[_j][_qp] * _test[_i][_qp];
      break;
 
    case Moose::NeighborElement:
      if (vdotn >= 0)
        r -= vdotn * _phi[_j][_qp] * _test_neighbor[_i][_qp];
      break;
 
    case Moose::NeighborNeighbor:
      if (vdotn < 0)
        r -= vdotn * _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
      break;
  }
 
  return r;
}