Sunday, July 7, 2013

Python Software Foundation - SfePy : GSoC 2013 : Week 3

Hi Pythonistas,

This week kicked off with the hunt for a robust Navier-Stokes solver for SfePy.
I went through huge amount of resources.

To summarise the journey:

For our Navier-Stokes currently we use the Newton method with backtracking line-search. in OpenFoam and most of the CFD code the linearization approach is based on Patankar's SIMPLE algorithm.[1][2]
I talked to my professor who told me that SIMPLE is used in commercial softwares like FLUENT too.

I found few papers which tells us some other approaches:


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[3] to solve the problem.

Other Implementations:

Here is a list of Open Source CFD codes. Maybe we can fork a repo and use it or learn from it:

According to people iNavier and dolphyn are promising:

Someone was using PyAMG to develop Jacobian-Free Newton-Krylov code to solve the Navier Stokes equations :!topic/pyamg-user/HXrXTyvXPpw


This is everything I could harness this week. There is a lot of things to take care to lock the final solver to be used which I would do the current week. Also I am currently narrowing down and rigorously searching a way to implementing SIMPLE in the FE context.

On a side note I have also been working on a SfePy version for Python 3 and benchmarking the simulation results.

This journey is surely turning out to be awesome!


No comments:

Post a Comment

Comment will be sent for moderation.