On Thu, 4 Jul 2013, Ankit Mahato wrote:

> Hi R,
>
> I did some more digging from the implementation point of view and came
> across some interesting things:
>
> This tutorial demonstrates the solution of Incompressible Navier-Stokes
> Equations using Fenics. it uses Chlorin's method[1] to solve the problem.
> This project was very awesome.�
> http://fenicsproject.org/documentation/dolfin/1.2.0/python/demo/pde/navier-stoke
> s/python/documentation.html

That looks feasible as well, although it is for time-dependent problems. A
stationary solution (if existing) could be obtained by stepping in time
til nothing changes.

(Sidenote: fenics is a cool project with many interesting ideas - good
place for insiration.)

> Then a stackoverflow question where someone says that fenics wasn't fast.
> http://stackoverflow.com/questions/4768045/fluid-flow-heat-transfer-and-python

As fenics is C++, I would say that it will always be faster than sfepy :)
I think that the "slow" in that context meant that the
nonlinearity solution converges slowly.

> Other Implementations:
>  *  Parallel Spectral Numerical Methods/The Two- and Three-Dimensional
>     Navier-Stokes Equations -http://en.wikibooks.org/wiki/Parallel_Spectral_Numerical_Methods/The_Two-_and_Th
>     ree-Dimensional_Navier-Stokes_Equations
>  *  2D Navier-Stokes solver implemented as a Python package with Python
>     modules and C++ extension modules. It uses the finite difference
>     method on a uniform, rectangular grid. It handles single- and
>     two-phase incompressible, Newtonian, laminar flow with obstacles.
>     -https://code.google.com/p/kmkns/

There is a thesis to download - might be interesting.

>  *  Finite Volume Based - http://www.ctcms.nist.gov/fipy/
> Here is a list of Open Source CFD codes. Maybe we can fork a repo and use
> it or learn from it:
> http://www.cfd-online.com/Wiki/Codes

I think it would be easier to follow a paper/thesis, as details in code
often differ. But you can try, yes. Note that sfepy is BSD-licensed, so we
cannot use snippets/functions from GPL-licensed codes.

> According to people iNavier and dolphyn are promising:
> http://www.cfd-online.com/Forums/main/13529-colver-code-c-c.html
>
> Someone was using PyAMG to develop Jacobian-Free Newton-Krylov code to
> solve the Navier Stokes equations :
> https://groups.google.com/forum/#!topic/pyamg-user/HXrXTyvXPpw

This could be really interesting - maybe you could ask the person on how
far that project got?

> This is everything I could harness till now. I hope something useful
> comes out of these. Also I am currently narrowing down and rigorously
> searching a way to implementing SIMPLE in the FE context.

Ok, thanks!

r.