AEFVKernel Class Reference

A dgkernel for the advection equation using a cell-centered finite volume method. More...

#include <AEFVKernel.h>

Inheritance diagram for AEFVKernel:

Public Member Functions
	AEFVKernel (const InputParameters &parameters)
virtual	~AEFVKernel ()

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

Protected Attributes
MooseEnum	_component
	choose an equation More...
const VariableValue &	_uc1
	piecewise constant variable values in cells More...
const VariableValue &	_uc2
const MaterialProperty< Real > &	_u1
	extrapolated variable values at side center More...
const MaterialProperty< Real > &	_u2
const InternalSideFluxBase &	_flux
	flux user object More...

Detailed Description

A dgkernel for the advection equation using a cell-centered finite volume method.

Notes:

1. This dgkernel itself does not do any complex calculation. It gets the internal side flux vector and Jacobian matrix from the internal side flux user object being called.
2. In general, a approximate Riemann solver should be used in the internal flux user object for calculating the flux
3. If a system of governing equations is being solved, the flux vector and Jacobian matrix are calculated only once for the first equation and cached for use for the rest of the equations in the system.
4. On the "left" and "right" states of the internal side, the variable values are interpolated from the reconstructed linear polynomials in the "left" and "right "element, respectively, which are provided from the corresponding material kernel.

Definition at line 43 of file AEFVKernel.h.

Constructor & Destructor Documentation

AEFVKernel()

AEFVKernel::AEFVKernel

(

const InputParameters &

parameters

)

Definition at line 28 of file AEFVKernel.C.

  : DGKernel(parameters),
    _component(getParam<MooseEnum>("component")),
    _uc1(coupledValue("u")),
    _uc2(coupledNeighborValue("u")),
    _u1(getMaterialProperty<Real>("u")),
    _u2(getNeighborMaterialProperty<Real>("u")),
    _flux(getUserObject<InternalSideFluxBase>("flux"))
{
}

~AEFVKernel()

AEFVKernel::~AEFVKernel ( )

virtual

Definition at line 39 of file AEFVKernel.C.

{}

Member Function Documentation

computeQpJacobian()

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

protectedvirtual

Definition at line 68 of file AEFVKernel.C.

{
  // assemble the input vectors, which are
  //   the constant monomial from the current and neighbor elements
  std::vector<Real> uvec1 = {_uc1[_qp]};
  std::vector<Real> uvec2 = {_uc2[_qp]};
 
  // calculate the Jacobian matrices
  const auto & fjac1 = _flux.getJacobian(Moose::Element,
                                         _current_side,
                                         _current_elem->id(),
                                         _neighbor_elem->id(),
                                         uvec1,
                                         uvec2,
                                         _normals[_qp]);
 
  const auto & fjac2 = _flux.getJacobian(Moose::Neighbor,
                                         _current_side,
                                         _current_elem->id(),
                                         _neighbor_elem->id(),
                                         uvec1,
                                         uvec2,
                                         _normals[_qp]);
 
  // distribute the contribution to the current and neighbor elements
  switch (type)
  {
    case Moose::ElementElement:
      return fjac1(_component, _component) * _phi[_j][_qp] * _test[_i][_qp];
 
    case Moose::ElementNeighbor:
      return fjac2(_component, _component) * _phi_neighbor[_j][_qp] * _test[_i][_qp];
 
    case Moose::NeighborElement:
      return -fjac1(_component, _component) * _phi[_j][_qp] * _test_neighbor[_i][_qp];
 
    case Moose::NeighborNeighbor:
      return -fjac2(_component, _component) * _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
  }
 
  return 0.0;
}

computeQpResidual()

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

protectedvirtual

Definition at line 42 of file AEFVKernel.C.

{
  // assemble the input vectors, which are
  //   the reconstructed linear monomial
  //   extrapolated at side center from the current and neighbor elements
  std::vector<Real> uvec1 = {_u1[_qp]};
  std::vector<Real> uvec2 = {_u2[_qp]};
 
  // calculate the flux
  const auto & flux = _flux.getFlux(
      _current_side, _current_elem->id(), _neighbor_elem->id(), uvec1, uvec2, _normals[_qp]);
 
  // distribute the contribution to the current and neighbor elements
  switch (type)
  {
    case Moose::Element:
      return flux[_component] * _test[_i][_qp];
 
    case Moose::Neighbor:
      return -flux[_component] * _test_neighbor[_i][_qp];
  }
 
  return 0.0;
}

Member Data Documentation

_component

MooseEnum AEFVKernel::_component

protected

choose an equation

Definition at line 54 of file AEFVKernel.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_flux

const InternalSideFluxBase& AEFVKernel::_flux

protected

flux user object

Definition at line 68 of file AEFVKernel.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_u1

const MaterialProperty<Real>& AEFVKernel::_u1

protected

extrapolated variable values at side center

Definition at line 64 of file AEFVKernel.h.

Referenced by computeQpResidual().

_u2

const MaterialProperty<Real>& AEFVKernel::_u2

protected

Definition at line 65 of file AEFVKernel.h.

Referenced by computeQpResidual().

_uc1

const VariableValue& AEFVKernel::_uc1

protected

piecewise constant variable values in cells

Definition at line 60 of file AEFVKernel.h.

Referenced by computeQpJacobian().

_uc2

const VariableValue& AEFVKernel::_uc2

protected

Definition at line 61 of file AEFVKernel.h.

Referenced by computeQpJacobian().

---

The documentation for this class was generated from the following files:

• rdg/include/dgkernels/AEFVKernel.h
• rdg/src/dgkernels/AEFVKernel.C