Staggered scheme local assembler using numeric differentiation and implicit time discretization.
More...
template<std::size_t id, class TypeTag, class Assembler>
class Dumux::SubDomainStaggeredLocalAssembler< id, TypeTag, Assembler, DiffMethod::numeric, true >
The assembly of the cellCenterResidual is done element-wise, the assembly of the face residual is done half-element-wise.
A sketch of what this means can be found in the following image:
Half-element wise assembly means, that integrals are split into contributions from the left and right part of the staggered control volume. For an example term \(\int_{\Omega}\varrho u\text{d}\Omega\) this reads \(\int_{\Omega}\varrho u\text{d}\Omega = \frac{1}{2}\Omega_\text{left}\varrho_\text{left}u+\frac{1}{2}\Omega_\text{right}\varrho_\text{right}u\).
During assembly, \(\frac{1}{2}\Omega_\text{left}\varrho_\text{left}u\) is added to the residual of the staggered control volume \(\Omega\), when the loops reach the scvf within element \(\Omega_\text{left}\). \(\frac{1}{2}\Omega_\text{right}\varrho_\text{right}u\) is added to the residual of \(\Omega\), when the loops reach the scvf within in element \(\Omega_\text{right}\). Other terms are split analogously.
|
| CellCenterResidualValue | assembleCellCenterResidualImpl () |
| FaceResidualValue | assembleFaceResidualImpl (const SubControlVolumeFace &scvf) |
| template<class JacobianMatrixDiagBlock, class GridVariables> |
| CellCenterResidualValue | assembleCellCenterJacobianAndResidualImpl (JacobianMatrixDiagBlock &A, GridVariables &gridVariables) |
| | Computes the derivatives with respect to the given element and adds them to the global matrix.
|
| template<class JacobianMatrixDiagBlock, class GridVariables> |
| auto | assembleFaceJacobianAndResidualImpl (JacobianMatrixDiagBlock &A, GridVariables &gridVariables) |
| | Computes the derivatives with respect to the given element and adds them to the global matrix.
|
| template<class JacobianBlock, class GridVariables> |
| void | assembleJacobianCellCenterCoupling (Dune::index_constant< faceId > domainJ, JacobianBlock &A, const CellCenterResidualValue &origResidual, GridVariables &gridVariables) |
| | Computes the derivatives with respect to the given element and adds them to the global matrix.
|
| template<std::size_t otherId, class JacobianBlock, class GridVariables> |
| void | assembleJacobianCellCenterCoupling (Dune::index_constant< otherId > domainJ, JacobianBlock &A, const CellCenterResidualValue &res, GridVariables &gridVariables) |
| template<class JacobianBlock, class ElementResidualVector, class GridVariables> |
| void | assembleJacobianFaceCoupling (Dune::index_constant< cellCenterId > domainJ, JacobianBlock &A, const ElementResidualVector &origResiduals, GridVariables &gridVariables) |
| | Computes the derivatives with respect to the given element and adds them to the global matrix.
|
| template<std::size_t otherId, class JacobianBlock, class ElementResidualVector, class GridVariables> |
| void | assembleJacobianFaceCoupling (Dune::index_constant< otherId > domainJ, JacobianBlock &A, const ElementResidualVector &res, GridVariables &gridVariables) |
| template<class JacobianMatrixDiagBlock, class GridVariables> |
| void | evalAdditionalDerivatives (const std::vector< std::size_t > &additionalDofDependencies, JacobianMatrixDiagBlock &A, GridVariables &gridVariables) |
| void | bindLocalViews () |
| void | assembleJacobianAndResidual (JacobianMatrixRow &jacRow, SubResidual &res, GridVariablesTuple &gridVariables) |
| | Computes the derivatives with respect to the given element and adds them to the global matrix. The element residual is written into the right hand side.
|
| void | assembleResidual (SubResidual &res) |
| | Assemble the residual only.
|
| CellCenterResidualValue | evalLocalResidualForCellCenter () const |
| | Convenience function to evaluate the complete local residual for the current element. Automatically chooses the the appropriate element volume variables.
|
| CellCenterResidualValue | evalLocalFluxAndSourceResidualForCellCenter () const |
| | Convenience function to evaluate the flux and source terms (i.e, the terms without a time derivative) of the local residual for the current element. Automatically chooses the the appropriate element volume and face variables.
|
| CellCenterResidualValue | evalLocalStorageResidualForCellCenter () const |
| | Convenience function to evaluate storage term (i.e, the term with a time derivative) of the local residual for the current element. Automatically chooses the the appropriate element volume and face variables.
|
| FaceResidualValue | evalLocalResidualForFace (const SubControlVolumeFace &scvf) const |
| | Convenience function to evaluate the local residual for the current face. Automatically chooses the the appropriate element volume and face variables.
|
| FaceResidualValue | evalLocalFluxAndSourceResidualForFace (const SubControlVolumeFace &scvf) const |
| | Convenience function to evaluate the flux and source terms (i.e, the terms without a time derivative) of the local residual for the current element. Automatically chooses the the appropriate element volume and face variables.
|
| FaceResidualValue | evalLocalStorageResidualForFace (const SubControlVolumeFace &scvf) const |
| | Convenience function to evaluate storage term (i.e, the term with a time derivative) of the local residual for the current face. Automatically chooses the the appropriate element volume and face variables.
|
| const Problem & | problem () const |
| ElementFaceVariables & | curElemFaceVars () |
| | The current element volume variables.
|
| ElementFaceVariables & | prevElemFaceVars () |
| | The element volume variables of the provious time step.
|
| CouplingManager & | couplingManager () |
| ElementResidualVector | evalLocalResidual () const |
| | Convenience function to evaluate the complete local residual for the current element. Automatically chooses the the appropriate element volume variables.
|
| ElementResidualVector | evalLocalFluxAndSourceResidual () const |
| | Convenience function to evaluate the flux and source terms (i.e, the terms without a time derivative) of the local residual for the current element. Automatically chooses the the appropriate element volume variables.
|
| ElementResidualVector | evalLocalStorageResidual () const |
| | Convenience function to evaluate storage term (i.e, the term with a time derivative) of the local residual for the current element. Automatically chooses the the appropriate element volume variables.
|
| void | enforceInternalDirichletConstraints (const ApplyFunction &applyDirichlet) |
| | Enforces Dirichlet constraints if enabled in the problem.
|
| const Assembler & | assembler () const |
| | The assembler.
|
| const Element & | element () const |
| | The current element.
|
| bool | elementIsGhost () const |
| | Returns if element is a ghost entity.
|
| const SolutionVector & | curSol () const |
| | The current solution.
|
| FVElementGeometry & | fvGeometry () |
| | The global finite volume geometry.
|
| ElementVolumeVariables & | curElemVolVars () |
| | The current element volume variables.
|
| ElementVolumeVariables & | prevElemVolVars () |
| | The element volume variables of the provious time step.
|
| ElementFluxVariablesCache & | elemFluxVarsCache () |
| | The element flux variables cache.
|
| LocalResidual & | localResidual () |
| | The local residual for the current element.
|
| ElementBoundaryTypes & | elemBcTypes () |
| | The element's boundary types.
|
1.15.0