Hi John,
_sys_solution is the system->solution vector.
It is fed by:
_system.solution->localize_to_one( _sys_solution );
On Fri, Aug 25, 2017 at 1:48 PM, John Peterson <jwpeter...@gmail.com> wrote:
>
>
> On Fri, Aug 25, 2017 at 10:34 AM, Renato Poli <rebp...@gmail.com> wrote:
>
>> Hi,
>>
>> I need some help here ...
>>
>> I am having trouble to understand what is going on with a boundary
>> integral.
>> I have solved a Poisson problem for pressure in a 2D mesh, as a model for
>> a
>> incompressible water flow in porous media.
>> The mesh (Tri6 elements) has a near-circular polygon in the middle
>> representing an injection well.
>>
>> RHS of the equations are zero except in boundaries dofs. Dirichlet
>> boundary
>> conditions are imposed using a penalty both on the well and on the
>> external
>> boundaries to represent fixed pressure. (the solution looks ok - seems to
>> work)
>>
>> To validate the solution, I am trying to estimate the water flow both in
>> the well and in the external boundaries (they must be identical and match
>> the analytical solution). I therefore calculate the normal pressure
>> gradient and integrate in the whole boundary path.
>>
>> When I use first order variable, the results are "almost correct". But
>> when
>> I use second order , it is absolutely off.
>>
>> I am using the code below to do the integration. I tried many
>> alternatives,
>> no success.
>>
>> Any idea of what I am doing wrong?
>>
>
> The code seems OK to me, although you should not need to call
> fe_face->reinit() inside the loop over boundary ids.
>
> What is _sys_solution?
>
> --
> John
>
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
