Dear Community
I know that there are richly documented tutorials for the solution of the Stokes equations.
I am curious to know if there any examples for solving the creeping flow problem in 2D axisymmetric geometry using the stream-function formulation. I vaguely remember seeing a demo problem that does this, but I am not able to locate it now.
Any help is greatly appreciated.
Thank You
Warm Regards
The issue is constructing a FE scheme which in some fashion supports a C^1 basis.
This demo serves to show more than the streamfunction formulation of the Stokes problem, but nevertheless it may be a good place to get started. The numerical scheme here weakly enforces C^1 continuity.
@kamensky’s excellent tIGAr may be a good place to look for ideas too. There, one has C^k, k\ge0, continuity for free by virtue of repeated knots in the spline basis knot vector. E.g. the biharmonic demo.
Thank you very much, I shall go through the suggested examples.
Warm Regards
I was wondering if there is any documentation to go with the demo link [streamfunction.py] you have provided above? It would help me understand the code better.
Thank You
Warm Regards