I have found out where was the problem, i had inverted s and t in the surface_dot expression.


Le vendredi 1 juillet 2016 13:41:35 UTC+2, Nicolas DELAN a écrit :
I have modified my file but it do not work, the app crash sending me an error message :

"ValueError: cannot match arguments! ([('heat_Coeff', 'h'), 't', 's'])"

I don't understand why : i have attached my files

Le vendredi 1 juillet 2016 13:05:58 UTC+2, Robert Cimrman a écrit :
Hi Nicolas,

On 06/30/2016 04:34 PM, Nicolas DELAN wrote:
> Hello,
>
> I have to formulate this equation in sfepy:
>
> <https://lh3.googleusercontent.com/--SRvfhvrVes/V3UtQhHns_I/AAAAAAAAVVI/URP8zAIqx2wcF6Cnmq9oP2-p4VVvax4UQCLcB/s1600/Capture_1.PNG>
>
> So with \partialOmega = gamma:
>
>
>
> """dw_laplace.i.Omega(c,s, t ) + dw_volume_dot.i.Gamma( h , t, s)
>
> =dw_surface_dot.i.Gamma(g, s )"""

 From the attached image it would seem to me more like:

"""dw_laplace.i.Omega(c,s, t ) + dw_surface_dot.i.Gamma( h , t, s)

= dw_surface_integrate.i.Gamma(g, s )"""


>
> My question is: alpha and g are constants of the surface of the airfoil, i
> have two files that contains the values of them for each point of the
> surface (we are in a 2D problem).
>
>
> But i don't know how to define h and alpha in the configuration file as
> variables of the domain gamma with some ebcs conditions ?
>

You need to define the h and g using functions, see [1] - get_pars() in
particular, and functions dict at the bottom (you need just the single
function). In your version of get_pars(), you would have to load the files, and
interpolate the values to the quadrature point coordinates `coors`.

r.

[1] http://sfepy.org/doc-devel/examples/diffusion/poisson_functions.html