[Nek5000-users] Smoothness of the derivatives
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jun 25 11:03:01 CDT 2018
Hi,
I guess your only chance is to increase the resolution, which will make
the discontinuities of the derivatives decrease. This is one of the
characteristics of spectral methods; you will see a lack of resolution
directly in the fields. Similar resolution with other methods will not
give you more accurate results, but they appear nicer (smoother).
Indeed, filtering does not help much. But if you want to try, make sure
you use a filter that actually goes beyond one element (e.g. using a
Laplace operator).
Philipp
On 2018-06-25 17:54, nek5000-users at lists.mcs.anl.gov wrote:
> Because I need to analyze the dissipation field
>
> понедельник, 25 июня 2018г., 20:58 +07:00 от
> nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov>:
>
> This is correct, derivatives in SEM are not C0. Your code forces the
> field to be C0 but it does not make it smooth across elements.
> Typically filtering doesn’t help. Why do you want a smooth
> derivative field?
>
>
> > On 25 Jun 2018, at 12:55, "nek5000-users at lists.mcs.anl.gov
> <mailto:nek5000-users at lists.mcs.anl.gov>"
> <nek5000-users at lists.mcs.anl.gov
> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
> >
> > Hi, Neks!
> >
> > 1) I am doing a DNS study of turbulent pipe jets and looking at
> turbulent kinetic energy budgets.
> >
> > I take derivatives via gradm1 subroutine, for example:
> >
> > call gradm1(uxdx,uxdy,uxdz,ux)
> >
> > After all, if I look at my dissipation field (via outpost), I see
> some wrong values on the surface between spectral elements.
> > I want to filter them using:
> >
> > call col2(uxdx,bm1,ntot)
> > call dssum(uxdx,lx1,ly1,lz1)
> > call col2(uxdx,binvm1,ntot)
> >
> > call dsavg(dxux)
> >
> > for every array. But conversely it leads to an increase in gaps.
> What am I doing wrong?
> >
> > 2) How does a filter-function in dssum look and how much points
> does it use for filtering?
> >
> > Best regards,
> > Vlad.
> > _______________________________________________
> >
> > Nek5000-users mailing list
> >
> > Nek5000-users at lists.mcs.anl.gov
> <mailto:Nek5000-users at lists.mcs.anl.gov>
> >
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> >
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov>
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
More information about the Nek5000-users
mailing list