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 :)
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...
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? :)
r.
Cheers, r.
Regards, Ankit
[1]
http://sfepy.org/doc-devel/tutorial.html#notes-on-solving-pdes-by-the-finite...