![](https://secure.gravatar.com/avatar/27f6b03a726ccff3f7a5007b1df6bf8c.jpg?s=120&d=mm&r=g)
On 07/02/2013 04:07 PM, Ankit Mahato wrote:
On Tuesday, 2 July 2013 19:23:05 UTC+5:30, Robert Cimrman wrote:
On 07/02/2013 03:46 PM, Ankit Mahato wrote:
On Tuesday, 2 July 2013 19:06:16 UTC+5:30, Ankit Mahato wrote:
On Tuesday, 2 July 2013 13:58:20 UTC+5:30, Robert Cimrman wrote:
Now it remains to implement a robust flow solver. Even this small
example
shows, that the solution is not obtained easily - try decreasing the viscosity, and/or increase the Dirichlet velocity - the solver would not converge.
Yes R,
The solution is not obtained easily. I am looking into it.
PS: Here are blog posts for week 1 & 2 Kindly tell me if this will do before I send it to terri oda:
http://ankitmahato.blogspot.in/2013/07/python-software-foundation-sfepy-gsoc...
http://ankitmahato.blogspot.in/2013/07/python-software-foundation-sfepy-gsoc...
Hi R,
Just wanted your views. Does the problem of Navier-Strokes solver lies with the implementation or the algorithm which is used.
Mostly the algorithm, but it might be also the formulation. I am far from CFD, but people there seem to be preferring a dimensionless form of the incompressible NS equations. It also depends on the discretization/FE spaces used. It's really a broad subject, and there is still no a silver bullet solver. Maybe ask your thesis advisor/colleagues doing CFD? Searching the net is really of no help here, as it returns so many things... Expert advice is needed :)
Yes R dimensionless form of the equations are preferred. I have implemented Patankar's SIMPLE algorithm previously (http://en.wikipedia.org/wiki/SIMPLE_algorithm) [also SIMPLER and SIMPLEC] when I did CFD course but using Finite Difference Method. I will look into it if I find FEM approach. Also I wanted to ask you if we need to stick to FEM for CFD because people use Finite Volume Method [FVM] for CFD.
Yes, we need it in a FEM context, sfepy cannot do FVM. But (big)IMHO those two approaches have things in common, and once you have the matrices and algorithm... It should even be possible to compute fluxes and similar quantities, if needed, though I never tried that and think there is no time to implement FVM in sfepy within your project.
r.