Skip to content

Feature/womersley: Pre-calculate Time invariant portion of Womersley Solution

Kurt Sansom requested to merge kayarre/nektar:feature/Womersley into master

Iterated on the Womersley boundary condition to pre-evaluate the time invariant portions of the analytical solution. Currently this assumes the Radius of the inlet does not change.

Changes include:

  • add NekComplexDouble as a complex type modeled after NekDouble but only in IncNavierStokes.h

  • adding additional arrays in the womersley parameters at initialization.

  • move the coefficient calculation to the SetUpWomersley function

  • removing unnecessary variables

A current deficit is that each field is initialized separately, it would be nice to be able to reduce the duplication that occurs with this implementation. The Womersley test passes.

Edited by Kurt Sansom

Merge request reports