Implementation of constant fields
Hello,
I have to formulate this equation in sfepy:
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 )"""
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 ?
Hi Nicolas,
On 06/30/2016 04:34 PM, Nicolas DELAN wrote:
Hello,
I have to formulate this equation in sfepy:
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
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:
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
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:
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
participants (2)
-
Nicolas DELAN -
Robert Cimrman