On 07/09/2013 12:29 PM, Ankit Mahato wrote:
I went through the codes I mentioned above: The following solvers were good but were licensed under GPL DUNS 2d/3d navier Stokes - GNU GPL license 2.0 Channelflow - GNU GPL2 OpenFlower - GPLv2 PETSc-FEM - GPLv2 Tochnog - GPL
Ankit, I think Fenics project's application cbcpdesys might be the best choice as it is in python and solvers are mainly called from python. Has a few turbulence models already coded. Might be a good idea to look how they did it (it is FE based). None of the codes you listed is python based. Another FE based solver that may be scripted with python is Fluidity from http://amcg.ese.ic.ac.uk/index.php?title=Fluidity again just to investigate how they do it.
You might want to think about using petsc solver suite (from Argonne National Labs) to be called from python as in fenics/dolfin. It has quite a few choices for solving nolinear systems. Its license seems very liberal (like bsd).
You'll need to solve not NS but Reynolds averaged NS. This ends up adding additional equations related to turbulence. Just NS by itself will diverge as the Reynolds number goes above a certain value. All of this already done in cbcpdesys. It is GNU LGPL so should be ok to use with bsd. But you won't be copying verbatim anyways.
just a few thoughts. I use fluent at work and many open cfd codes at home. Regards, Osman