On 08/29/2012 02:05 PM, David Libault wrote:
> Hi Robert,
> From the poisson.py I would like to probe the flux of grad(T) over a
> surface (Gamma_Left or Gamma_right), but I am not sure which term to use...
IMHO it's the same term as in my answer to Alec (thread "Question on the
py''") - look at the post_process() function - the
term is d_surface_flux.
> I would use it in electrostatics (temperature is replaced by voltage), to
> compute the resistance of the volume imposing a voltage of 1 across two
> surfaces with dirichlet conditions and computing the electrical current.
> For a ohmic conductor, j = \sigma E = - \sigma grad(V), \sigma being the
> conductivity. j being the current density, the current intensity thru a
> surface is the flux of j thru that surface. For the case of homogenous
> \sigma, the term asked for above would do...
OK, let us know if d_surface_flux works for you - it can take a general
> Le dimanche 26 août 2012 22:55:49 UTC+2, Robert Cimrman a écrit :
>> On 08/26/2012 12:32 PM, Alec Kalinin wrote:
>>> Hi Robert,
>>> Did you mean "linear_elastic_probes.html" instead of
html" example? I found the
>>> "linear_elastic_probes.html" very useful example for my purposes to
>> probe a
>>> solution in the given (x, y, z) points. Also the documentation
>>> "src/sfepy/fem/probes.html" gives all necessary information to help me
>>> implement what I want to do. Thank you!
>> Sorry, I cut&pasted a wrong url, the correct one is . But you found
>> another one that solves the problem.
>>> But, despite this, could you tell me more about low-level way to
>> evaluate a
>>> variable in the given (x, y, z) point?
>> It's exactly how the probes do that: the key function is
>> variable.evaluate_at() , where variable is an unknown or parameter
>> variable. It takes just one compulsory parameter - the coordinates of
>> points in which you wish to evaluate the variable. You can get the
>> variables of a problem by problem.get_variables(), where problem is the
>> second argument of the post_process_hook function.
>> Best regards,
>>  http://sfepy.org/doc-devel/
>>  http://sfepy.org/doc-devel/
>>> On Sunday, August 26, 2012 12:47:50 AM UTC+4, Robert Cimrman wrote:
>>>> Hi Alec,
>>>> On 08/25/2012 05:45 PM, Alec Kalinin wrote:
>>>>> Dear SfePy users,
>>>>> Is it possible to evaluate a solution not only in the FEM mesh node,
>>>>> any arbitrary point in the domain with the given (x, y, z)
>>>> Yes, it is possible. Either, you could use a probe as described in the
>>>> Primer  - the available probes are described in . Or, you could
>>>> directly evaluate a variable in given points - this is a bit low-level
>>>> operation, but I could provide you instructions, if the probes are not
>>>> enough for you.
>>>>> For example, consider Dirichlet problem for Poisson equation. We apply
>>>>> essential boundary conditions on the surface nodes and after the
>>>>> has been solved we have the solution vector, i.e. vector of values in
>>>>> FEM mesh nodes. But I want to know the solution in point v(x, y, z)
>>>>> not FEM mesh node. What is the best way to obtain solution in this
>>>>> Alec Kalinin
>>>>  doc-devel/examples/linear_
>>>>  http://sfepy.org/doc-devel/