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]

> 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.

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...

Cheers,
r.

> Regards,
> Ankit
>

[1]
http://sfepy.org/doc-devel/tutorial.html#notes-on-solving-pdes-by-the-finite-element-method