Upwind numerical flux scheme for the advection equation using a cell-centered finite volume method. More...

#include <AEFVUpwindInternalSideFlux.h>

Inheritance diagram for AEFVUpwindInternalSideFlux:

Public Member Functions
	AEFVUpwindInternalSideFlux (const InputParameters &parameters)

virtual	~AEFVUpwindInternalSideFlux ()

virtual void	calcFlux (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave, std::vector< Real > &flux) const override
	Solve the Riemann problem. More...

virtual void	calcJacobian (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave, DenseMatrix< Real > &jac1, DenseMatrix< Real > &jac2) const override
	Compute the Jacobian matrix. More...

virtual void	execute () override

virtual void	initialize () override

virtual void	finalize () override

virtual void	threadJoin (const UserObject &) override

virtual const std::vector< Real > &	getFlux (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave) const
	Get the flux vector. More...

virtual const DenseMatrix< Real > &	getJacobian (Moose::DGResidualType type, unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave) const
	Get the Jacobian matrix. More...

Protected Attributes
const Real	_velocity
	advective velocity More...

unsigned int	_cached_flux_elem_id
	element ID of the cached flux values More...

unsigned int	_cached_flux_neig_id
	neighbor element ID of the cached flux values More...

unsigned int	_cached_jacobian_elem_id
	element ID of the cached Jacobian values More...

unsigned int	_cached_jacobian_neig_id
	neighbor element ID of the cached Jacobian values More...

std::vector< Real >	_flux
	flux vector of this side More...

DenseMatrix< Real >	_jac1
	Jacobian matrix contribution to the "left" cell. More...

DenseMatrix< Real >	_jac2
	Jacobian matrix contribution to the "right" cell. More...

Detailed Description

Upwind numerical flux scheme for the advection equation using a cell-centered finite volume method.

Definition at line 25 of file AEFVUpwindInternalSideFlux.h.

Constructor & Destructor Documentation

◆ AEFVUpwindInternalSideFlux()

AEFVUpwindInternalSideFlux::AEFVUpwindInternalSideFlux ( const InputParameters & parameters )

Definition at line 25 of file AEFVUpwindInternalSideFlux.C.

   : InternalSideFluxBase(parameters), _velocity(getParam<Real>("velocity"))
 {
 }

◆ ~AEFVUpwindInternalSideFlux()

AEFVUpwindInternalSideFlux::~AEFVUpwindInternalSideFlux ( )

virtual

Definition at line 30 of file AEFVUpwindInternalSideFlux.C.

30 {}

Member Function Documentation

◆ calcFlux()

void AEFVUpwindInternalSideFlux::calcFlux	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave,
		std::vector< Real > &	flux
	)		const

overridevirtual

Solve the Riemann problem.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal
[out]	flux	flux vector across the side

Implements InternalSideFluxBase.

Definition at line 33 of file AEFVUpwindInternalSideFlux.C.

 {
   mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
   mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
  
   // assign the size of flux vector, e.g. = 1 for the advection equation
   flux.resize(1);
  
   // assume a constant velocity on the left
   RealVectorValue uadv1(_velocity, 0.0, 0.0);
  
   // assume a constant velocity on the right
   RealVectorValue uadv2(_velocity, 0.0, 0.0);
  
   // normal velocity on the left and right
   Real vdon1 = uadv1 * dwave;
   Real vdon2 = uadv2 * dwave;
  
   // calculate the so-called a^plus and a^minus
   Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
   Real amins = 0.5 * (vdon2 - std::abs(vdon2));
  
   // finally calculate the flux
   flux[0] = aplus * uvec1[0] + amins * uvec2[0];
 }

◆ calcJacobian()

void AEFVUpwindInternalSideFlux::calcJacobian	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave,
		DenseMatrix< Real > &	jac1,
		DenseMatrix< Real > &	jac2
	)		const

overridevirtual

Compute the Jacobian matrix.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal
[out]	jac1	Jacobian matrix contribution to the "left" cell
[out]	jac2	Jacobian matrix contribution to the "right" cell

Implements InternalSideFluxBase.

Definition at line 66 of file AEFVUpwindInternalSideFlux.C.

 {
   mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
   mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
  
   // assign the size of Jacobian matrix, e.g. = (1, 1) for the advection equation
   jac1.resize(1, 1);
   jac2.resize(1, 1);
  
   // assume a constant velocity on the left
   RealVectorValue uadv1(_velocity, 0.0, 0.0);
  
   // assume a constant velocity on the right
   RealVectorValue uadv2(_velocity, 0.0, 0.0);
  
   // normal velocity on the left and right
   Real vdon1 = uadv1 * dwave;
   Real vdon2 = uadv2 * dwave;
  
   // calculate the so-called a^plus and a^minus
   Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
   Real amins = 0.5 * (vdon2 - std::abs(vdon2));
  
   // finally calculate the Jacobian matrix
   jac1(0, 0) = aplus;
   jac2(0, 0) = amins;
 }

◆ execute()

void InternalSideFluxBase::execute ( )

overridevirtualinherited

Definition at line 40 of file InternalSideFluxBase.C.

41 {

42 }

◆ finalize()

void InternalSideFluxBase::finalize ( )

overridevirtualinherited

Definition at line 45 of file InternalSideFluxBase.C.

46 {

47 }

◆ getFlux()

const std::vector< Real > & InternalSideFluxBase::getFlux	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave
	)		const

virtualinherited

Get the flux vector.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal

Definition at line 55 of file InternalSideFluxBase.C.

 {
   if (_cached_flux_elem_id != ielem || _cached_flux_neig_id != ineig)
   {
     _cached_flux_elem_id = ielem;
     _cached_flux_neig_id = ineig;
  
     calcFlux(iside, ielem, ineig, uvec1, uvec2, dwave, _flux);
   }
   return _flux;
 }

Referenced by AEFVKernel::computeQpResidual().

◆ getJacobian()

const DenseMatrix< Real > & InternalSideFluxBase::getJacobian	(	Moose::DGResidualType	type,
		unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave
	)		const

virtualinherited

Get the Jacobian matrix.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal

Definition at line 73 of file InternalSideFluxBase.C.

 {
   if (_cached_jacobian_elem_id != ielem || _cached_jacobian_neig_id != ineig)
   {
     _cached_jacobian_elem_id = ielem;
     _cached_jacobian_neig_id = ineig;
  
     calcJacobian(iside, ielem, ineig, uvec1, uvec2, dwave, _jac1, _jac2);
   }
  
   if (type == Moose::Element)
     return _jac1;
   else
     return _jac2;
 }

Referenced by AEFVKernel::computeQpJacobian().

◆ initialize()

void InternalSideFluxBase::initialize ( )

overridevirtualinherited

Definition at line 31 of file InternalSideFluxBase.C.

 {
   _cached_flux_elem_id = libMesh::invalid_uint;
   _cached_flux_neig_id = libMesh::invalid_uint;
   _cached_jacobian_elem_id = libMesh::invalid_uint;
   _cached_jacobian_neig_id = libMesh::invalid_uint;
 }

◆ threadJoin()

void InternalSideFluxBase::threadJoin ( const UserObject & )

overridevirtualinherited

Definition at line 50 of file InternalSideFluxBase.C.

51 {

52 }

Member Data Documentation

◆ _cached_flux_elem_id

unsigned int InternalSideFluxBase::_cached_flux_elem_id

mutableprotectedinherited

element ID of the cached flux values

Definition at line 115 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux(), and InternalSideFluxBase::initialize().

◆ _cached_flux_neig_id

unsigned int InternalSideFluxBase::_cached_flux_neig_id

mutableprotectedinherited

neighbor element ID of the cached flux values

Definition at line 117 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux(), and InternalSideFluxBase::initialize().

◆ _cached_jacobian_elem_id

unsigned int InternalSideFluxBase::_cached_jacobian_elem_id

mutableprotectedinherited

element ID of the cached Jacobian values

Definition at line 120 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian(), and InternalSideFluxBase::initialize().

◆ _cached_jacobian_neig_id

unsigned int InternalSideFluxBase::_cached_jacobian_neig_id

mutableprotectedinherited

neighbor element ID of the cached Jacobian values

Definition at line 122 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian(), and InternalSideFluxBase::initialize().

◆ _flux

std::vector<Real> InternalSideFluxBase::_flux

mutableprotectedinherited

flux vector of this side

Definition at line 125 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux().

◆ _jac1

DenseMatrix<Real> InternalSideFluxBase::_jac1

mutableprotectedinherited

Jacobian matrix contribution to the "left" cell.

Definition at line 127 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian().

◆ _jac2

DenseMatrix<Real> InternalSideFluxBase::_jac2

mutableprotectedinherited

Jacobian matrix contribution to the "right" cell.

Definition at line 129 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian().

◆ _velocity

const Real AEFVUpwindInternalSideFlux::_velocity

protected

advective velocity

Definition at line 50 of file AEFVUpwindInternalSideFlux.h.

Referenced by calcFlux(), and calcJacobian().

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

rdg/include/userobjects/AEFVUpwindInternalSideFlux.h
rdg/src/userobjects/AEFVUpwindInternalSideFlux.C

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ AEFVUpwindInternalSideFlux()

◆ ~AEFVUpwindInternalSideFlux()

Member Function Documentation

◆ calcFlux()

◆ calcJacobian()

◆ execute()

◆ finalize()

◆ getFlux()

◆ getJacobian()

◆ initialize()

◆ threadJoin()

Member Data Documentation

◆ _cached_flux_elem_id

◆ _cached_flux_neig_id

◆ _cached_jacobian_elem_id

◆ _cached_jacobian_neig_id

◆ _flux

◆ _jac1

◆ _jac2

◆ _velocity