[Nek5000-users] Unrealistic values at element edges

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Dec 10 13:50:01 CST 2014


Hi Paul,

Thank you very much for your reply. I'm doing DNS of turbulent pipe flow and the plots were from a pipe flow with Re_D=5300 (Re_tau=180). Basically, I used the same grid resolution as specified in the paper "Direct numerical simulation of turbulent pipe flow at moderately high Reynolds numbers" which I think should be reasonably good for the given low Reynolds number. And I'm also using PnPn-2, the same as used in the papers which use nek5000.

In this case, shall I use even better grid resolution for this Reynolds number in order to get rid of the wiggles?

For Lambda2, I'll have a look at the filtering subroutines you suggested.

Best regards,
Tony

________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users-request at lists.mcs.anl.gov <nek5000-users-request at lists.mcs.anl.gov>
Sent: 10 December 2014 18:00
To: nek5000-users at lists.mcs.anl.gov
Subject: Nek5000-users Digest, Vol 70, Issue 6

Send Nek5000-users mailing list submissions to
        nek5000-users at lists.mcs.anl.gov

To subscribe or unsubscribe via the World Wide Web, visit
        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
        nek5000-users-request at lists.mcs.anl.gov

You can reach the person managing the list at
        nek5000-users-owner at lists.mcs.anl.gov

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Nek5000-users digest..."


Today's Topics:

   1. Re: Unrealistic values at element edges
      (nek5000-users at lists.mcs.anl.gov)


----------------------------------------------------------------------

Message: 1
Date: Wed, 10 Dec 2014 17:34:21 +0000
From: nek5000-users at lists.mcs.anl.gov
To: "nek5000-users at lists.mcs.anl.gov"
        <nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] Unrealistic values at element edges
Message-ID:
        <mailman.11825.1418232863.3424.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="us-ascii"


Tony,

The wiggles in the quantities your plotting are generally indicative of marginal
resolution for the given Reynolds number.  They will go away with increased resolution.
Lambda2 is not continuous a priori, but converges to a continuous function with increased
resolution (exponentially fast).   Are you using PnPn or PnPn-2 ?   Pressure will generally
look better with PnPn (and be better, in the marginally resolved, i.e., LES, case).

For lambda2 it might be possible to do some filtering to get a smoother plot...(same for
other quantities).   I would suggest going that route --- there are several filtering strategies
but I'm not certain which would be best for your particular case.     You might start with the
g_filter() routine in navier5.f with a diagonal scaling of the wavenumbers like:

  ( 1  1  1   1   1 ... 1  .99 .95  .90 )

say... which reduces the energy in the top 3 modes.  There is also filter_d2()
that you can experiment with.

Paul

________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Wednesday, December 10, 2014 10:28 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Unrealistic values at element edges

Dear Neks,

Sorry that the image wasn't shown. I've attached it here. Hope someone can help.

Thanks and regards,

Tony

________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users-request at lists.mcs.anl.gov <nek5000-users-request at lists.mcs.anl.gov>
Sent: 10 December 2014 16:23
To: nek5000-users at lists.mcs.anl.gov
Subject: Nek5000-users Digest, Vol 70, Issue 4

Send Nek5000-users mailing list submissions to
        nek5000-users at lists.mcs.anl.gov

To subscribe or unsubscribe via the World Wide Web, visit
        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
        nek5000-users-request at lists.mcs.anl.gov

You can reach the person managing the list at
        nek5000-users-owner at lists.mcs.anl.gov

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Nek5000-users digest..."


Today's Topics:

   1. Unrealistic values at element edges
      (nek5000-users at lists.mcs.anl.gov)


----------------------------------------------------------------------

Message: 1
Date: Wed, 10 Dec 2014 16:23:02 +0000
From: nek5000-users at lists.mcs.anl.gov
To: "nek5000-users at lists.mcs.anl.gov"
        <nek5000-users at lists.mcs.anl.gov>
Subject: [Nek5000-users] Unrealistic values at element edges
Message-ID:
        <mailman.11813.1418228589.3424.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="iso-8859-1"

Dear Neks,


I've been running turbulent channel and pipe flows using nek5000 for a while. However, recently when I created movies of lambda2 structures and pressure iso-surfaces in channel and pipe flows I observed some "strange" values at the edge between elements.


[cid:cec5092f-41be-4ae9-a468-a5fa9ab9f7ff]

This is lambda2 structures in turbulent pipe flow.


It is clearly visible from the picture that at some slices in streamwise direction, the values don't seem to be realistic and it feels like the flow is passing through some filters. Interestingly, when I plotted the mesh together with the variables, the vertical noise lines happened to be at the edge between adjacent elements, which is particularly clear in lambda2 plot. Same behaviour was also found in 2D pressure plots and in turbulent channel flow cases. Has anyone observed similar things before in visualisation?


I used both Paraview and VisIt for visualisation and the same thing was observed, so it shouldn't be a problem of software. In this case, I'm just wondering why the result is like this or is it because I set something wrong in my simulations?


Any comments will be appreciated and thank you very much in advance.


Best regards,

Tony
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20141210/8312bf82/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: l2.jpg
Type: image/jpeg
Size: 66013 bytes
Desc: l2.jpg
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20141210/8312bf82/attachment.jpg>

------------------------------

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users


End of Nek5000-users Digest, Vol 70, Issue 4
********************************************


------------------------------

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users


End of Nek5000-users Digest, Vol 70, Issue 6
********************************************


More information about the Nek5000-users mailing list