[Nek5000-users] Smoothness of the derivatives
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Jun 26 09:03:12 CDT 2018
Vlad,
dssum will just sum up the shared points between elements - so no filtering. Note, currently there is no easy accessable mechanism in place to make your field smooth. As pointed out already modal filtering won't help.
Cheers,
Stefan
-----Original message-----
> From:nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
> Sent: Tuesday 26th June 2018 15:34
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] Smoothness of the derivatives
>
> Hi Philipp!
>
> Thank you for your answer!
>
> In fact I have almost the DNS resolution (the grid spacing is comparable with the Kolmogorov scale from the literature data) and rather small error in the dissipation field. But I want to see a smooth dissipation, so I really want to try to filter it.
>
> As I understood filtering procedure is in dssum subroutine. But looking at it didn't help to me to realize the process of the filtering.
> So could you explain, please, how can I choose the filter to be used? Cause now I think that my filter doesn't go beyond one element and it is the reason why it makes my dissipation field worse.
>
> Best regards,
> Vlad
>
>
>
>
> Понедельник, 25 июня 2018, 23:05 +07:00 от nek5000-users at lists.mcs.anl.gov:
>
> 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 <mailto: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> <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>
> > <mailto:nek5000-users at lists.mcs.anl.gov>"
> > <nek5000-users at lists.mcs.anl.gov <mailto: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>
> > <mailto:Nek5000-users at lists.mcs.anl.gov>
> > >
> > > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <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> <mailto:Nek5000-users at lists.mcs.anl.gov>
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <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 <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 <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