On Friday, 13 September 2013 16:40:43 UTC+5:30, Robert Cimrman wrote:
Hi Ankit,

thanks for the update. Nice movies!

:)
 

Now we should think how to integrate this into the sfepy package - using
Material for material parameters,

Okie R.
 
Mesh/MeshIO for VTK output, Solver framework
for the solver, numpy(-like) docstring standard for docstrings etc.

On it R.
 
Do you
think you will get to some of those before the deadline?

Yes R. I will get as far as I can before the deadline.
 
As for the finite
volumes used - what are the main obstacles of converting that to finite elements?


A key component which had to be implemented was the enthalpy update scheme. As stated in my blog post, in research papers it was stated and tested keeping FVM base in mind.
Also I am more versed with FVM than FEM so I could implement it in that time frame I had along with the benchmarking without which this code would have had no credibility.
There will definitely be a way using FEM, but wish I had more time in this GSoC timeframe. I have already kept it in my mind as a challenge over which I will ponder as I still have 9 months left to complete my thesis.

Cheers.
  
 
Cheers,
r.

On 09/12/2013 07:32 PM, Ankit Mahato wrote:
> Hi R,
>
> Sorry to keep you waiting.
> Here is the much awaited update -
> http://ankitmahato.blogspot.in/2013/09/python-software-foundation-sfepy-gsoc_12.html
>
> Regards,
> Ankit
>
> On Friday, 6 September 2013 21:54:32 UTC+5:30, Robert Cimrman wrote:
>>
>> On 09/06/2013 03:59 PM, Ankit Mahato wrote:
>>> Hi R,
>>>
>>> In the weak form:
>>> - The navier stokes equation will now have a source term.
>>> - The energy equation will also have a source term (slight mistake in my
>>> previous file. corrected now.) The H field is (=cT).
>>> i need to calculate u,p,H (=cT), gl (liquid fraction) for each element.
>>> Iterating the energy equation to obtain convergence in the nodal liquid
>>> fraction value (source term) is one the key areas where I was having
>> crisis
>>> situation.
>>
>> So is it resolved?
>>
>> If things get too complicated, try some simplifications so that we have
>> some
>> example(s) within SfePy for the GSoC final phase. I will be able to
>> support you
>> fully from 11.9.
>>
>>> The momentum and energy equations are coupled through both the buoyancy
>> and
>>> the source term (Bu). The B depends on the liquid fraction which
>> represents
>>> an extra unknown in the governing equations.
>>> Since the liquid fraction is updated such that the corection term
>> becomes
>>> zero at the convergence of the enthalpy equation.
>>>
>>> I tried searching the internet but could not find any FEM implementation
>> is
>>> this field so I started developing a python solver from scratch in order
>> to
>>> test the algorithm suggested in those papers using finite volume method,
>> as
>>> at least I will have a solver ready to couple it with sfepy.
>>
>> Ok, keep up the good work, and let us know about progress or problems.
>> Once
>> that solver works, you will have something to compare a new FEM solution
>> with.
>>
>> Cheers,
>> r.
>>
>>> regards,
>>> Ankit
>>> On Friday, 6 September 2013 13:21:01 UTC+5:30, Robert Cimrman wrote:
>>>>
>>>> Hi Ankit,
>>>>
>>>> it looks very interesting. Do you have also a weak formulation of the
>>>> complete
>>>> problem? I would like to see the form
>>>>
>>>> Find u, p, ... such that ... holds for all v, q...
>>>>
>>>> to see if all the terms are available. Or do you have already an
>> initial
>>>> implementation? I will be offline from 7. to 11.9. but Vladimir should
>> be
>>>> able
>>>> to answer your questions. I am available today.
>>>>
>>>> r.
>>>>
>>>> On 09/02/2013 10:53 PM, Ankit Mahato wrote:
>>>>> Hi R,
>>>>>
>>>>> Here is the final problem statement file containing all details.
>>>>>
>>>>> regards,
>>>>> Ankit
>>
>>
>