Re: [Libmesh-users] Trouble with boundary integration

2017-08-26 Thread Renato Poli
node > if you expect it to be uniform rather than oscillating :) > > > From: Renato Poli > Sent: Friday, August 25, 2017 7:46 PM > To: Mike Marchywka > Cc: John Peterson; libmesh-users > Subject: Re: [Libmesh-users] Trouble with boundary integrat

Re: [Libmesh-users] Trouble with boundary integration

2017-08-26 Thread Mike Marchywka
flux at each node if you expect it to be uniform rather than oscillating :) From: Renato Poli Sent: Friday, August 25, 2017 7:46 PM To: Mike Marchywka Cc: John Peterson; libmesh-users Subject: Re: [Libmesh-users] Trouble with boundary integration On Fri

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
) apply_normal= -3.14973e-05 > > > ____ > From: Renato Poli > Sent: Friday, August 25, 2017 6:12 PM > To: Mike Marchywka > Cc: John Peterson; libmesh-users > Subject: Re: [Libmesh-users] Trouble with boundary integration > &g

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Mike Marchywka
Poli Sent: Friday, August 25, 2017 6:12 PM To: Mike Marchywka Cc: John Peterson; libmesh-users Subject: Re: [Libmesh-users] Trouble with boundary integration Hi Mike, Thanks for the reply. Any hint on how to easily reproduce the shape function calculations by hand (that is, manually producing th

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
On Fri, Aug 25, 2017 at 7:19 PM, John Peterson wrote: > > > On Fri, Aug 25, 2017 at 1:02 PM, Renato Poli wrote: > >> Hi John, >> >> I am still struggling and running out of ideas here. >> > > Are you using FIRST and SECOND, LAGRANGE variables? If so, are you also > changing the mesh when you swi

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread John Peterson
On Fri, Aug 25, 2017 at 1:02 PM, Renato Poli wrote: > Hi John, > > I am still struggling and running out of ideas here. > Are you using FIRST and SECOND, LAGRANGE variables? If so, are you also changing the mesh when you switch back and forth between orders? -- John ---

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
___ > From: Renato Poli > Sent: Friday, August 25, 2017 3:02 PM > To: John Peterson > Cc: libmesh-users > Subject: Re: [Libmesh-users] Trouble with boundary integration > > Hi John, > > I am still struggling and running out of ideas here. > I wrote down some deb

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Mike Marchywka
ith p and h refinements etc. Although maybe someone can clarify a bit. Thanks. From: Renato Poli Sent: Friday, August 25, 2017 3:02 PM To: John Peterson Cc: libmesh-users Subject: Re: [Libmesh-users] Trouble with boundary integration Hi John, I am

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
Hi John, I am still struggling and running out of ideas here. I wrote down some debugging - perhaps you can find something weird. The calculations seem correct, at least to what I could understand so far. But things are zero'ing out. If, in the mesh creation, I declare the variable as "FIRST" or

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
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 wrote: > > > On Fri, Aug 25, 2017 at 10:34 AM, Renato Poli wrote: > >> Hi, >> >> I need some help here ... >> >> I a

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread John Peterson
On Fri, Aug 25, 2017 at 10:34 AM, Renato Poli 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. >

[Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Renato Poli
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 repr