<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Paul, <br>
<br>
Thanks paul. I tried working with the diverging mesh modification.<br>
<br>
I tried to use this form of equation to modify the mesh <br>
<br>
z = mx + b<br>
z = ny + c <br>
<br>
x = ( z - b ) / m <br>
y = (z - c) / n<br>
<br>
Do you think would work ? or is there any other way to input that ?<br>
<br>
Thanks, <br>
<br>
Kamal<br>
<br>
<br>
<div class="moz-cite-prefix">On 17/01/2014 14:14,
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
</div>
<blockquote
cite="mid:mailman.10711.1389964496.3278.nek5000-users@lists.mcs.anl.gov"
type="cite">
<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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.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>
+-----+-----+---+--+---+------+------+---->
z<br>
<br>
<br>
Paul<br>
<br>
<br>
<br>
On Wed, 15 Jan 2014, <a moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.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 moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.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
moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.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 moz-do-not-send="true"
href="http://nek5000.mcs.anl.gov/index.php/UG#Documentation"
target="_blank">http://nek5000.mcs.anl.gov/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>
________________________________________<br>
From: <a moz-do-not-send="true"
href="mailto:nek5000-users-bounces@lists.mcs.anl.gov"
target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>
[<a moz-do-not-send="true"
href="mailto:nek5000-users-bounces@lists.mcs.anl.gov"
target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>]
on behalf of <a moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.gov</a>
[<a moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.gov</a>]<br>
Sent: Monday, January 13, 2014 3:35 AM<br>
To: <a moz-do-not-send="true"
href="mailto:nek5000-users@lists.mcs.anl.gov"
target="_blank">nek5000-users@lists.mcs.anl.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>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
</blockquote>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Nek5000-users@lists.mcs.anl.gov"
target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a moz-do-not-send="true"
href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users"
target="_blank">https://lists.mcs.anl.gov/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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Nek5000-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
</pre>
</blockquote>
<br>
</body>
</html>