Kilian Weishaupt authored Jan 11, 2020 and Ned Coltman committed Jan 23, 2020

* include scvf of velocity dof parallel to boundary
* rename velocityPorousMedium to porousMediumVelocity and return velocity vector
* deprecate old interfaces

* remove previously deprecated interface bjsVelocity
* clean-up includes

changes only affect documentation

Kilian Weishaupt authored Jul 02, 2019 and Kilian Weishaupt committed Jul 19, 2019

* not done for getParallelVelocityFromOtherBoundary -> will raise
  deprecation warning

* account for the tangential velocity gradient along the interface
* refactor some functions, add more docu

* introduce a new function for the slip velocity that also takes into account a later gradient
* introduce beta = alpha/sqrt(K)
* introduce a velocity in the porous medium (0 by default)

Kilian Weishaupt authored Jun 17, 2019 and Melanie Lipp committed Jun 27, 2019

* this should make it easier to properly reconstruct a stress at the boundary

* correctly account for the case when there is a Neumann bc for both
  the normal and tangential part of the stress.

* so far, interior boundaries (e.g., flow around a block)
  were always interpreted as Dirichlet (no-slip) boundaries
* the reduced distance used in the calculation of the gradient
  was already correctly incorporated, we just document it better here
* this concerns the gradient for the shear stress term, e.g., du/dy
  at an edge within the domain

* is defined for all scvfs, not only for those on a boundary
* remove redundant code

Kilian Weishaupt authored Apr 11, 2019 and Ned Coltman committed Apr 25, 2019

* pass only the face normal to our scvf of interest
* calculate distance within function

* when passing the result of partial(sol) to the staggered vtkWriter,
  one has to make sure the type of the sol vec (Dune::MultitypeBlockVector
  storing references) passed to the constructor
  is the same as the one used as template argument
* failing to do so results in a segfault
* this commit ensures the correct types at compile time

* Use boundary velocity v_j in term (v_i*v_j) both for upwinding and
  central average
* do not use upwinding for kovasznay test
* update reference solutions for kovasznay and angeli test

* Remove unused code and simplify

* before, fluxes were called for each eqIdx

Kilian Weishaupt authored Mar 06, 2019 and Ned Coltman committed Mar 11, 2019

* use small indicies for loops
* do not query for dofs using < 0 to check if they exist
* remove superfluous this-> calls
* use efficient types, clean up loops
* move axis and pair data to namespace Detail