<div dir="ltr">Kamal,<div><br></div><div>1) To check for convergence, run your simulation for increasing polynomial orders ( lx1 in the SIZE file) and measure some consistent quantity (drag, TKE, etc). You can plot this quantity against lx1 and this should give you a good picture of how you are converging (without needing to change the physical mesh).</div>
<div><br></div><div>2) If you weren't using temperature in your simulation, you could do put the following code snippet in your userchk</div><div><br></div><div>c-------------------------------------------</div><div>integer i, ntot</div>
<div><br></div><div>ntot = nx1*ny1*nz1*nelv ! total number of points in simulation (on THIS processor)</div><div><br></div><div>do i = 1,ntot</div><div> t(i,1,1,1,1) = 0.5*(vx(i,1,1,1)**2 + vy(i,1,1,1)**2 + vz(i,1,1,1)**2)</div>
<div>enddo</div><div><br></div><div>ifto = .true.</div><div>c-------------------------------------------</div><div><br></div><div>I've assumed here that your density is unity. This computed kinetic energy point-wise and stores it in the temperature field. The part "ifto = .true." tells Nek to output the temperature field whenever it outputs a field file. You could then view this in visit just like you would the velocity field.</div>
<div><br></div><div>Hope this helps!</div><div><br></div><div>Josh</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 15, 2014 at 9:52 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Paul,<br>
<br>
I will try that mapping, I created a straight pipe to do some simulations.<br>
<br>
1) Is there a way to check my spectral convergence, so that I can do some mesh convergence study?<br>
2) Is there a way to plot the kinetic energy in the flow ?<br>
<br>
Thank you,<br>
Kamal<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
On 15/01/2014 12:14, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Hi Kamal,<br>
<br>
For the geometry you describe, I would apply mesh morphing<br>
to a standard straight pipe. You can create such a pipe<br>
using a subset of the mkmesh script provided in the expansion<br>
example.<br>
<br>
If your transition region is conical, then I woud make<br>
certain that the change in slope occurs at element boundaries,<br>
i.e., as in the half cross-section depicted below:<br>
<br>
_______________<br>
/| | |<br>
/ | | |<br>
/| | | |<br>
________________/ | | | |<br>
| | | | | | | |<br>
| | | | | | | |<br>
+-----+-----+---+--+---+------<u></u>+------+----> z<br>
<br>
<br>
Paul<br>
<br>
<br>
<br>
On Wed, 15 Jan 2014, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Paul,<br>
<br>
I looked at the example and it looks very helpful, but for a sudden expansion we have only two parts (i.e) the inlet pipe and the outlet pipe. Whereas in a diverging pipe we have three section (i.e,) The inlet pipe, Diverging section and the outlet pipe.<br>
<br>
<br>
Is it possible for me to create a 2D diverging section using mesh modification and revolve it around the centre axis to get a 3d diverging pipe ?<br>
<br>
<br>
or should I apply the mesh modification directly to the 3d sudden expansion pipe using usrdat2() ?<br>
<br>
Thanks,<br>
Kamal<br>
<br>
<br>
On 14/01/2014 21:09, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Paul,<br>
<br>
Thanks a lot for the update. I will take a look at as soon as possible and keep you updated with it.<br>
<br>
Thank you,<br>
<br>
Kamal<br>
<br>
On Jan 14, 2014, at 9:04 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Kamal,<br>
<br>
I've added an example for a step-type expansion in a pipe. New example directory<br>
is nek5_svn/examples/expansion, with a README.<br>
<br>
If you want a standard diffuser, I would use mesh-morphing, as outlined in the<br>
nek primer, which can be found on this page:<br>
<br>
<a href="http://nek5000.mcs.anl.gov/index.php/UG#Documentation" target="_blank">http://nek5000.mcs.anl.gov/<u></u>index.php/UG#Documentation</a><br>
<br>
I've also automated our turbulent outflow bc treatment -- one simply provides an<br>
expansion factor (typically ~ 1.5) and two persistent arrays to turb_outflow(), which<br>
is now in navier5.f. This can be useful when very strong vortices are leaving the<br>
domain, which is possible for certain diffuser flow conditions (and others).<br>
<br>
Paul<br>
______________________________<u></u>__________<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.<u></u>mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.<u></u>mcs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a>]<br>
Sent: Monday, January 13, 2014 3:35 AM<br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
Subject: [Nek5000-users] Creating a Mesh<br>
<br>
Hi<br>
<br>
Can some one please help me in this because I am new to nek<br>
<br>
1) How can I create a .rea for a specific geometry say ' A sudden<br>
expansion pipe ' or a ' Diverging pipe '<br>
<br>
2) What tool I should use to create that mesh.<br>
<br>
<br>
Thank you<br>
<br>
<br>
kamal<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
<br>
</blockquote>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke
</div>