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

#include <AEFVMaterial.h>

Inheritance diagram for AEFVMaterial:

Public Member Functions
	AEFVMaterial (const InputParameters &parameters)

virtual	~AEFVMaterial ()

Protected Member Functions
virtual void	computeQpProperties ()

Protected Attributes
const VariableValue &	_uc

const SlopeLimitingBase &	_lslope

MaterialProperty< Real > &	_u

Detailed Description

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

Notes:

This material kernel is responsible for invoking the in-cell slope reconstruction and slope limiting based on the cell-average variable values.
The reconstructed/limited linear polynomial variable in each element is then interpolated to the center of each side that surrounds an element.
If a system of governing equations is being solved, the reconstructed/limited vector of variable gradients is only calculated once when the first equation is dealt with, and then cached for use for the rest of the equations.
If reconstruction/limiting scheme is not turned on, the cell-average constant variable is used for flux calculations, which is first-order accurate in space and absolutely stable

Definition at line 43 of file AEFVMaterial.h.

Constructor & Destructor Documentation

◆ AEFVMaterial()

AEFVMaterial::AEFVMaterial ( const InputParameters & parameters )

Definition at line 29 of file AEFVMaterial.C.

   : Material(parameters),
     _uc(coupledValue("u")),
     _lslope(getUserObject<SlopeLimitingBase>("slope_limiting")),
     _u(declareProperty<Real>("u"))
 {
 }

◆ ~AEFVMaterial()

AEFVMaterial::~AEFVMaterial ( )

virtual

Definition at line 37 of file AEFVMaterial.C.

37 {}

Member Function Documentation

◆ computeQpProperties()

void AEFVMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 40 of file AEFVMaterial.C.

 {
   // initialize the variable
   _u[_qp] = _uc[_qp];
  
   // interpolate variable values at face center
   if (_bnd)
   {
     // you should know how many equations you are solving and assign this number
     // e.g. = 1 (for the advection equation)
     unsigned int nvars = 1;
     std::vector<RealGradient> ugrad(nvars, RealGradient(0., 0., 0.));
     ugrad = _lslope.getElementSlope(_current_elem->id());
  
     // get the directional vector from cell center to face center
     RealGradient dvec = _q_point[_qp] - _current_elem->centroid();
  
     // calculate the variable at face center
     _u[_qp] += ugrad[0] * dvec;
  
     // clear the temporary vectors
     ugrad.clear();
   }
   // calculations only for elemental output
   else if (!_bnd)
   {
   }
 }