Re: [Libmesh-users] higher order elements

2013-05-15 Thread John Peterson
On Wed, May 15, 2013 at 2:27 PM, Manav Bhatia wrote: > > It is dNi/dx that I plotted. Actually, for this 1D case, the dNi/dxi > and dNi/dx are simply scaled variants of each other by a constant (the > factor 2/h that you mentioned). The point being that the inverse of this "constant" gets sm

Re: [Libmesh-users] higher order elements

2013-05-15 Thread Manav Bhatia
On Wed, May 15, 2013 at 4:13 PM, John Peterson wrote: > On Wed, May 15, 2013 at 11:58 AM, Manav Bhatia > wrote: > > > > I have attached a pdf which plots the expression for a 2-noded and a > > 3-noded element using Lagrange interpolation functions. The element has > unit > > length. The 2-no

Re: [Libmesh-users] higher order elements

2013-05-15 Thread Manav Bhatia
On Wed, May 15, 2013 at 3:54 PM, lorenzo alessio botti < [email protected]> wrote: > I'm pretty sure that the geometry approximation is responsible of the > convergence degradation. I've experienced this problem with the ringleb > flow (2D Euler) using inviscid wall bcs. If you do p re

Re: [Libmesh-users] higher order elements

2013-05-15 Thread John Peterson
On Wed, May 15, 2013 at 11:58 AM, Manav Bhatia wrote: > > I have attached a pdf which plots the expression for a 2-noded and a > 3-noded element using Lagrange interpolation functions. The element has unit > length. The 2-noded value is constant at 2 for the entire domain of the > element, whi

Re: [Libmesh-users] higher order elements

2013-05-15 Thread lorenzo alessio botti
I'm pretty sure that the geometry approximation is responsible of the convergence degradation. I've experienced this problem with the ringleb flow (2D Euler) using inviscid wall bcs. If you do p refinement instead of h refinement you should see that the entropy error actually stagnates and you get

Re: [Libmesh-users] higher order elements

2013-05-15 Thread Manav Bhatia
On Wed, May 15, 2013 at 1:17 PM, John Peterson wrote: > On Wed, May 15, 2013 at 6:36 AM, Manav Bhatia > wrote: > >> > >> > >> One other thought: do you have some O(h) stabilization terms present? > >> Then if you aren't refining the grid they don't get smaller... > >> > > > > > > There is a GLS

Re: [Libmesh-users] higher order elements

2013-05-15 Thread John Peterson
On Wed, May 15, 2013 at 6:36 AM, Manav Bhatia wrote: >> >> >> One other thought: do you have some O(h) stabilization terms present? >> Then if you aren't refining the grid they don't get smaller... >> > > > There is a GLS stabilization term with the typical tau matrix. I have > experimented with

Re: [Libmesh-users] higher order elements

2013-05-15 Thread Manav Bhatia
> > > > One other thought: do you have some O(h) stabilization terms present? > Then if you aren't refining the grid they don't get smaller... > > There is a GLS stabilization term with the typical tau matrix. I have experimented with different tau definitions, but get the same behavior. Essent

Re: [Libmesh-users] higher order elements

2013-05-14 Thread John Peterson
On May 14, 2013, at 12:37 PM, Manav Bhatia wrote: >> > >> > This is the Gaussian bump problem from the higher-order CFD workshop >> > (http://dept.ku.edu/~cfdku/hiocfd/case_c1.1.html). You are correct that >> > away >> > from the bump the boundary is straight, so linear elements should be fine

Re: [Libmesh-users] higher order elements

2013-05-14 Thread Manav Bhatia
> > > > > This is the Gaussian bump problem from the higher-order CFD workshop > > (http://dept.ku.edu/~cfdku/hiocfd/case_c1.1.html). You are correct that > away > > from the bump the boundary is straight, so linear elements should be > fine. I > > am looking at the entropy error, since the entrop

Re: [Libmesh-users] higher order elements

2013-05-14 Thread John Peterson
On Tue, May 14, 2013 at 11:54 AM, Manav Bhatia wrote: > > > On Tue, May 14, 2013 at 12:59 PM, John Peterson > wrote: >> >> On Tue, May 14, 2013 at 9:25 AM, Manav Bhatia >> wrote: >> > Hi, >> > >> > I am working on higher-order simulation using the hierarchich >> > function >> > on quad8s. >>

Re: [Libmesh-users] higher order elements

2013-05-14 Thread Roy Stogner
On Tue, 14 May 2013, Manav Bhatia wrote: > This is the Gaussian bump problem from the higher-order CFD workshop ( > http://dept.ku.edu/~cfdku/hiocfd/case_c1.1.html). You are correct that away > from the bump the boundary is straight, so linear elements should be fine. > I am looking at the entropy

Re: [Libmesh-users] higher order elements

2013-05-14 Thread Manav Bhatia
On Tue, May 14, 2013 at 12:59 PM, John Peterson wrote: > On Tue, May 14, 2013 at 9:25 AM, Manav Bhatia > wrote: > > Hi, > > > > I am working on higher-order simulation using the hierarchich > function > > on quad8s. > > > > My error-convergence plots give me theoretical convergence for fi

Re: [Libmesh-users] higher order elements

2013-05-14 Thread John Peterson
On Tue, May 14, 2013 at 9:25 AM, Manav Bhatia wrote: > Hi, > > I am working on higher-order simulation using the hierarchich function > on quad8s. > > My error-convergence plots give me theoretical convergence for first > and second order p, but the error stagnates for p > 2. I am speculat

[Libmesh-users] higher order elements

2013-05-14 Thread Manav Bhatia
Hi, I am working on higher-order simulation using the hierarchich function on quad8s. My error-convergence plots give me theoretical convergence for first and second order p, but the error stagnates for p > 2. I am speculating that this might be due to the lower-order geometry, assuming t