Commit faaf60b6 authored by Timo Koch's avatar Timo Koch
Browse files

[test][periodic] Overload neumann as currently necessary

parent ef0c8a5c
......@@ -127,6 +127,37 @@ public:
return source;
* \brief Evaluate the boundary conditions for a neumann
* boundary segment.
* This is the method for the case where the Neumann condition is
* potentially solution dependent
* \param element The finite element
* \param fvGeometry The finite-volume geometry
* \param elemVolVars All volume variables for the element
* \param elemFluxVarsCache Flux variables caches for all faces in stencil
* \param scvf The sub control volume face
* Negative values mean influx.
* E.g. for the mass balance that would be the mass flux in \f$ [ kg / (m^2 \cdot s)] \f$.
* TODO: We actually just want to use this from the base problem. However, the staggered
* problem inherits from FVProblem which uses a different (the default) NumEqVector.
* Therefore not overloading this would return a vector of wrong size.
template<class ElementVolumeVariables, class ElementFluxVariablesCache>
NumEqVector neumann(const Element& element,
const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars,
const ElementFluxVariablesCache& elemFluxVarsCache,
const SubControlVolumeFace& scvf) const
// forward it to the interface with only the global position
return NumEqVector(0.0);
// \}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment