On Friday, 5 July 2013 15:39:07 UTC+5:30, Robert Cimrman wrote:
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-s...
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.)
Are you suggesting to use the dolfin module or to use Chlorin's method?
Then a stackoverflow question where someone says that fenics wasn't fast.
http://stackoverflow.com/questions/4768045/fluid-flow-heat-transfer-and-pyth...
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.
Okie :)
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-_an... 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.
Yes, but it is using Difference method.
- 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.
So we can use any code with BSD license. right?
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?
i had already dropped a mail on her email-id. Waiting for her reply.
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.
So, it looks there are lot of things I need to look in depth now.