Post-Processing from avg_all files using Nek / Visit
Hi, I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking.. If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ? Regards Shriram
Hi Shiriam, It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation. So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing... Best, Hank On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]> wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Hank, Thanks. I think I can automate (with python) till I get Visit to do a line plot. Do you know a way I can just export the data from line plot to a text file ? I might have about 10 line plots (900 time steps) which I would like to automate. Regards Shriram Jagannathan On 30 December 2010 12:15, <[email protected]> wrote:
Hi Shiriam,
It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation.
So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing...
Best, Hank
On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]> wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Shririam, I would do: OpenDatabase("file.nek5000") # use real filename AddPlot("Pseudocolor", "pressure") # or whatever variable you want DrawPlots() Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2) SetActiveWindow(2) sw = SaveWindowAttributes() sw.format = sw.CURVE sw.family = 0 sw.fileName = "myname_line1" SetSaveWindowAttributes(sw) # saves out myname_line1.curve SaveWindow() DeleteAllPlots() SetActiveWindow(1) Lineout((-1, 1, -1), (2, 4, 5)) SetActiveWindow(2) sw.fileName = "myname_line2" SetSaveWindowAttributes(sw) SaveWindow() DeleteAllPlots() SetActiveWindow(1) .... In VisIt 2.2, this will be much easier ... you won't have to make a Pseudocolor plot and do a Lineout query ... you will be able to add the Lineout directly via a curve plot. Best, Hank On Thu, Dec 30, 2010 at 10:22 AM, <[email protected]> wrote:
Hi Hank,
Thanks. I think I can automate (with python) till I get Visit to do a line plot. Do you know a way I can just export the data from line plot to a text file ? I might have about 10 line plots (900 time steps) which I would like to automate.
Regards Shriram Jagannathan
On 30 December 2010 12:15, <[email protected]> wrote:
Hi Shiriam,
It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation.
So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing...
Best, Hank
On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]>wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Hank, Thank you very much for the script. I will try this now. We have VisIt 2.1 on local machines but the cluster is powered with VisIt 1.12 only. Thanks again ! Regards Shriram Jagannathan On 30 December 2010 12:30, <[email protected]> wrote:
Hi Shririam,
I would do: OpenDatabase("file.nek5000") # use real filename AddPlot("Pseudocolor", "pressure") # or whatever variable you want DrawPlots() Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2) SetActiveWindow(2) sw = SaveWindowAttributes() sw.format = sw.CURVE sw.family = 0 sw.fileName = "myname_line1" SetSaveWindowAttributes(sw) # saves out myname_line1.curve SaveWindow() DeleteAllPlots() SetActiveWindow(1) Lineout((-1, 1, -1), (2, 4, 5)) SetActiveWindow(2) sw.fileName = "myname_line2" SetSaveWindowAttributes(sw) SaveWindow() DeleteAllPlots() SetActiveWindow(1) ....
In VisIt 2.2, this will be much easier ... you won't have to make a Pseudocolor plot and do a Lineout query ... you will be able to add the Lineout directly via a curve plot.
Best, Hank
On Thu, Dec 30, 2010 at 10:22 AM, <[email protected]> wrote:
Hi Hank,
Thanks. I think I can automate (with python) till I get Visit to do a line plot. Do you know a way I can just export the data from line plot to a text file ? I might have about 10 line plots (900 time steps) which I would like to automate.
Regards Shriram Jagannathan
On 30 December 2010 12:15, <[email protected]> wrote:
Hi Shiriam,
It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation.
So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing...
Best, Hank
On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]>wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Hank, It works w/o any issues. Thanks. Regards Shriram Jagannathan On 30 December 2010 12:44, Shriram Jagannathan <[email protected]>wrote:
Hi Hank,
Thank you very much for the script. I will try this now. We have VisIt 2.1 on local machines but the cluster is powered with VisIt 1.12 only. Thanks again !
Regards Shriram Jagannathan
On 30 December 2010 12:30, <[email protected]> wrote:
Hi Shririam,
I would do: OpenDatabase("file.nek5000") # use real filename AddPlot("Pseudocolor", "pressure") # or whatever variable you want DrawPlots() Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2) SetActiveWindow(2) sw = SaveWindowAttributes() sw.format = sw.CURVE sw.family = 0 sw.fileName = "myname_line1" SetSaveWindowAttributes(sw) # saves out myname_line1.curve SaveWindow() DeleteAllPlots() SetActiveWindow(1) Lineout((-1, 1, -1), (2, 4, 5)) SetActiveWindow(2) sw.fileName = "myname_line2" SetSaveWindowAttributes(sw) SaveWindow() DeleteAllPlots() SetActiveWindow(1) ....
In VisIt 2.2, this will be much easier ... you won't have to make a Pseudocolor plot and do a Lineout query ... you will be able to add the Lineout directly via a curve plot.
Best, Hank
On Thu, Dec 30, 2010 at 10:22 AM, <[email protected]>wrote:
Hi Hank,
Thanks. I think I can automate (with python) till I get Visit to do a line plot. Do you know a way I can just export the data from line plot to a text file ? I might have about 10 line plots (900 time steps) which I would like to automate.
Regards Shriram Jagannathan
On 30 December 2010 12:15, <[email protected]> wrote:
Hi Shiriam,
It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation.
So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing...
Best, Hank
On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]>wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Cool! Keep the questions coming and I'll help out with whatever I can. It's very little effort for me to put together little scripts like that... Best, Hank On Thu, Dec 30, 2010 at 11:03 AM, <[email protected]> wrote:
Hi Hank,
It works w/o any issues. Thanks.
Regards Shriram Jagannathan
On 30 December 2010 12:44, Shriram Jagannathan <[email protected]>wrote:
Hi Hank,
Thank you very much for the script. I will try this now. We have VisIt 2.1 on local machines but the cluster is powered with VisIt 1.12 only. Thanks again !
Regards Shriram Jagannathan
On 30 December 2010 12:30, <[email protected]> wrote:
Hi Shririam,
I would do: OpenDatabase("file.nek5000") # use real filename AddPlot("Pseudocolor", "pressure") # or whatever variable you want DrawPlots() Lineout((1,1,1),(2,2,2)) # takes lineout from point (1,1,1) to point (2,2,2) SetActiveWindow(2) sw = SaveWindowAttributes() sw.format = sw.CURVE sw.family = 0 sw.fileName = "myname_line1" SetSaveWindowAttributes(sw) # saves out myname_line1.curve SaveWindow() DeleteAllPlots() SetActiveWindow(1) Lineout((-1, 1, -1), (2, 4, 5)) SetActiveWindow(2) sw.fileName = "myname_line2" SetSaveWindowAttributes(sw) SaveWindow() DeleteAllPlots() SetActiveWindow(1) ....
In VisIt 2.2, this will be much easier ... you won't have to make a Pseudocolor plot and do a Lineout query ... you will be able to add the Lineout directly via a curve plot.
Best, Hank
On Thu, Dec 30, 2010 at 10:22 AM, <[email protected]>wrote:
Hi Hank,
Thanks. I think I can automate (with python) till I get Visit to do a line plot. Do you know a way I can just export the data from line plot to a text file ? I might have about 10 line plots (900 time steps) which I would like to automate.
Regards Shriram Jagannathan
On 30 December 2010 12:15, <[email protected]> wrote:
Hi Shiriam,
It is pretty easy to automate this in VisIt. As we discussed at the user's meeting, VisIt treats each spectral element as a 5x5x5, 7x7x7, etc group of hexahedrons. Within each hexahedron, there is linear interpolation.
So: if you want help automating this in VisIt, that should be easy. As far as accuracy, Nek will clearly do the best job, but hopefully you now understand better what VisIt is doing...
Best, Hank
On Thu, Dec 30, 2010 at 10:12 AM, <[email protected]>wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Shriram, I see that Hank answered some of your questions already. If we define u = <u> + u' then we have <u'>^2 = <u^2> - <u>^2 with <u^2> stored in rms and <u> stored in avg. (Similar forms exist for the cross terms.) postx can load the rms and avg files to give you <u'>^2 -- but not yet with the .f0000 format. Another common approach is to manipulate these by loading all the avg files inside the solver itself. We can set up an example that illustrates this. At that point, you could also interpolate the output onto a finer mesh before passing to visit. Paul On Thu, 30 Dec 2010, [email protected] wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
Z_ kk kk http://cuireadhchunceoil.ie/index.php?option=com_content&view=article&id=58
Hi Paul, Thanks. I did set up a python script (thanks to Hank) in VisIt to do the interpolation and o/p as txt file. Does rmsblah.f0000 file have the other components of Reynolds stress viz. <uv>, <vw>, <wu> ? Yes, we have a routine for getting the fluctuating components from nek and outside of it as well. Regards Shriram On 31 December 2010 09:27, <[email protected]> wrote:
Hi Shriram,
I see that Hank answered some of your questions already.
If we define u = <u> + u' then we have
<u'>^2 = <u^2> - <u>^2
with <u^2> stored in rms and <u> stored in avg. (Similar forms exist for the cross terms.) postx can load the rms and avg files to give you <u'>^2 -- but not yet with the .f0000 format.
Another common approach is to manipulate these by loading all the avg files inside the solver itself. We can set up an example that illustrates this. At that point, you could also interpolate the output onto a finer mesh before passing to visit.
Paul
On Thu, 30 Dec 2010, [email protected] wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
Z_
kk
kk http://cuireadhchunceoil.ie/index.php?option=com_content&view=article&id=58
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Hi Shriram, the rm2 file has uv vw wu Look for "outpost" in navier5.f, which is where the avg_all() routine is found. Note that, because of space limitations, not _all_ of the <uT> <vT> <wT> fields are computed/stored.... but these clearly could be added for someone interested in turbulent Nu computations. Cheers, Paul On Fri, 31 Dec 2010, [email protected] wrote:
Hi Paul,
Thanks. I did set up a python script (thanks to Hank) in VisIt to do the interpolation and o/p as txt file. Does rmsblah.f0000 file have the other components of Reynolds stress viz. <uv>, <vw>, <wu> ? Yes, we have a routine for getting the fluctuating components from nek and outside of it as well.
Regards Shriram
On 31 December 2010 09:27, <[email protected]> wrote:
Hi Shriram,
I see that Hank answered some of your questions already.
If we define u = <u> + u' then we have
<u'>^2 = <u^2> - <u>^2
with <u^2> stored in rms and <u> stored in avg. (Similar forms exist for the cross terms.) postx can load the rms and avg files to give you <u'>^2 -- but not yet with the .f0000 format.
Another common approach is to manipulate these by loading all the avg files inside the solver itself. We can set up an example that illustrates this. At that point, you could also interpolate the output onto a finer mesh before passing to visit.
Paul
On Thu, 30 Dec 2010, [email protected] wrote:
Hi,
I would like to post-process quantities from avg_all files (rms, rm2, avg). Particularly I am looking at line plots (mainly) or contours on a plane. I had a question regarding the approach that I might be taking..
If I load the files, say rmsblah.f0000* into nek using the load subroutine and fill the ux,uy,uz arrays with that of <u^2>, <v^2>, <w^2> and call hpts(), how would the interpolation work ? Would it interpolate <u^2> with the interpolation that is used for velocities ?. The other option would be to do a line plot from Visit and export it to a text file. Also, I am not sure how Visit interpolates. Which is more likely to be accurate ?
Regards Shriram
Z_
kk
kk http://cuireadhchunceoil.ie/index.php?option=com_content&view=article&id=58
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
participants (1)
-
nek5000-users＠lists.mcs.anl.gov