On Friday, 28 June 2013 03:10:54 UTC+5:30, Robert Cimrman wrote:On 06/27/2013 10:27 PM, Ankit Mahato wrote:
>
>
> On Friday, 28 June 2013 01:17:19 UTC+5:30, Robert Cimrman wrote:
>>
>> On 06/27/2013 08:17 PM, Ankit Mahato wrote:
>>> Hi R,
>>>
>>> I tried out the Steady Axial convection and diffusion in slug flow with
>>> velocity :math:`U _0` in an insulated pipe. Entry and Exit temperatures
>> are
>>> given.
>>> \\ :math:`T _0 for x \leqslant 0`
>>> \\ :math:`T _1 for x \geqslant L`
>>>
>> https://github.com/animator/sfepy/blob/master/examples/ankit/slug_flow.py
>>> I had the query on how to set the lateral surface of the cylinder to
>>> satisfy the condition \alpha \frac {\partial T}{\partial n} = 0.
>>
>> Short answer: do nothing. Longer answer: [1]
>>
>>
> Oops. I confused \tau D. That was so stupid of me.
I do not understand :)
Sorry Typo I meant \Gamma _D .
>>> Can you please tell what get_pars and get_load_variable is doing.
>>
>> You probably want to use a different right-hand side term (dw_volume_lvf
>> for
>> volume force, the surface force term is missing, but could be easily
>> added, if
>> needed). Otherwise get_pars() computes the position-dependent material
>> parameter (load.f), and get_load_variable() just demonstrates, how to use
>> parameter variables and how to set their values by a function - you do not
>> need
>> that either.
>>
>
> got it!
>
>
>>
>> BTW. look at the patch that adds the dw_convect_v_grad_s term - it
>> captures
>> everything that usually needs to be done to add a term...
>>
>>
> Looking at it R.
>
> So energy eq is working. Should I try and develop few more examples or get
> into the navier strokes part.
IMHO solving the Navier-Stokes will be the "hard" part of the problem, as a new
solver will probably be needed. So the sooner you start with that, the better.
The only problem for now is, that the NS terms work in 3D only...
Yes. Kindly give me a head start to begin it.
Other thing: I fetched your code and the history is rather wild (use gitk --all
to see it) - I suggest you to checkout the current origin master branch, make
your branch (phase_change?) on top of it, and then cherry-pick the commits
adding the example files (maybe rename the directory to examples/phase_change).
There should be no merges... Is it clear? :)
Yes. it is wild as I was trying to fetch the latest code from main sfepy repo and merge it. then things got a bit messy. I will make a new branch for phase_change.
r.
>> Cheers,
>> r.
>>
>>> Regards,
>>> Ankit
>>>
>>
>> [1]
>>
>> http://sfepy.org/doc-devel/tutorial.html#notes-on-solving-pdes-by-the-finite-element-method
>>
>