<div dir="ltr"><div><div><div><div>Hi Barak,<br><br></div>Your err-run.log clearly shows CFL number increases considerably. So, it is indeed a problem of time-stepping stability. So, as you suggest a variable time-stepping might be necessary. you could even write a small routine in .usr file to control your timestep DT. If you decrease your error of divergence, it has nothing to do with simulation time. You can control simulation time with NSTEPS of .rea file. However for explicit method the region of stability of higher accuracy methods are smaller. So, 3rd order accurate method is a little less stable. You could also try changing TORDER = 2 , when your jet develops considerably.<br>
<br></div>I heard lagrangian formuation could be used in nek5000, but unfortunately I do not know that. It will be nice if you post this question of lagrangian part in an open format, so that the best person in the knowledge of this may answer you.<br>
<br></div>Thanks,<br><br></div>Tanmoy<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Dec 9, 2013 at 3:59 PM, <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 Tanmoy,<br>
It is the 4th time that the program stops (probably because of two big dt). See the end of the log of the run. The time step is dt=1e-3. The first test it was 0.02. It seems that the time step has to be reduced when the jet develops considerably.<br>
Perhaps, the best thing to do is to simulate till the time limit when the time step becomes not enough small, to stop, and then restart with the next time step (decreased) till the next stage and so on.<br>
Now I have further questions:<br>
1. I noticed that there is a parameter to control the error of divergence. If I choose smaller value, how it ncrease the simulation time?<br>
2. I am interested in the Lagrangian properties of the flow, i.e. calculating velocity, vorticity or other velocity derivatives at the Lagrangian orbit. Is it possible and which are the available tools?<br>
<br>
Barak<div class="im"><br>
<br>
<br>
On 12/09/2013 02:56 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Hi Barak,<br>
<br>
I am not aware of variable time stepping in nek5000. The time stepping<br>
used is BDFk/EXTk with CFL constraint for stability. You can change<br>
TORDER in .rea file to 2 or 3 to make the time stepping 2nd or third<br>
order accurate. There is a characteristic based time-stepping though...<br>
IFCHAR in .rea file.<br>
<br>
However, if for a value ot DT you are getting convergence, essentially<br>
you should get the same for smaller value of DT. What is the error you<br>
are getting after the program aborts after 1200 iterations? It is still<br>
stored in your log file I suppose.<br>
<br>
Thanks,<br>
Tanmoy<br>
<br>
<br>
On Sun, Dec 8, 2013 at 12:26 PM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a><br></div><div class="im">
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>>> wrote:<br>
<br>
Hi Tanmoy,<br></div><div class="im">
I have divided by 4 and again my run stopped at iteration 1200. My<br>
question is: does nek5000 have variable time-stepping?<br>
Thanks,<br>
Barak<br>
<br>
<br></div>
On 12/07/2013 08:37 PM, nek5000-users@lists.mcs.anl.__<u></u>gov<div class="im"><br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>> wrote:<br>
<br>
Hi Barak,<br>
<br></div><div class="im">
it is because of the fact that the generated file format is<br>
probably not<br>
readable in visit. In the *usrchk* subroutine in your .usr file<br>
write<br>
<br>
these 2 lines<br>
param(66) = 4<br>
param(67) = 4<br>
<br>
before the end of the routine<br>
<br>
Then compile and execute your program. Now you will get output<br>
files in<br>
the form of turbJet.fld001, turbJet.fld002, ...and so on.<br>
type in visnek turbJet in your terminal, and generate<br>
*turbJet.nek5000*.<br>
<br>
This .nek5000 file can be read by visit.<br>
<br>
Thanks,<br>
Tanmoy<br>
<br>
<br>
On Sat, Dec 7, 2013 at 10:33 AM,<br></div><div class="im">
<nek5000-users@lists.mcs.anl._<u></u>_gov<br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>><br>
<mailto:<a href="mailto:nek5000-users@lists." target="_blank">nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>>>> wrote:<br>
<br></div><div class="im">
Hi Tanmoy,<br>
Many thanks. It seems that after dividing by 4 (dt=0.005)<br>
it runs<br>
and I get output files in the form jet0.fxxxxx .<br>
Now I am trying to visualize them using the precompiled<br>
ubuntu11<br>
version of visit. Unfortunately, I get the following error:<br>
<br>
VisIt could not read from the file<br>
<br></div>
"/home/galanti/COMPUTATION/___<u></u>_NEK5000/nek5_svn/TEST-____<u></u>examples/turbJet/jet0.f*<div class="im"><br>
<br>
database".<br>
<br>
The generated error message was:<br>
<br>
There was an error opening<br>
<br></div>
/home/galanti/COMPUTATION/____<u></u>NEK5000/nek5_svn/TEST-____<u></u>examples/turbJet/jet0.f*<div class="im"><br>
<br>
database. It may be an invalid file. VisIt tried using the<br>
following file format readers to open the file: Nek5000<br>
<br>
The following error(s) may be helpful in identifying the<br>
problem:<br>
When attempting to use "Nek", the file matched the general<br>
format<br>
expected by the reader but eventually encountered this<br>
error: "A tag<br>
called filetemplate: must be specified".<br>
<br>
What's wrong?<br>
Many thanks again,<br>
Barak<br>
<br>
<br>
<br></div>
On 12/07/2013 12:42 AM, nek5000-users@lists.mcs.anl.__<u></u>__gov<div class="im"><br>
<br>
<mailto:<a href="mailto:nek5000-users@lists." target="_blank">nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>>> wrote:<br>
<br></div><div class="im">
Hi Barak,<br>
<br>
It is a problem with numerical Stability. Decrease DT<br>
in your<br>
.rea file<br>
until the simulations become stable. You could start<br>
with 1/2 DT to<br>
check if simulation is stable, other wise go for 1/4<br>
DT, and<br>
check with<br>
this simple test.<br>
<br>
Thanks,<br>
Tanmoy<br>
<br>
<br>
On Fri, Dec 6, 2013 at 3:31 PM,<br></div>
<nek5000-users@lists.mcs.anl._<u></u>___gov<br>
<mailto:<a href="mailto:nek5000-users@lists." target="_blank">nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>>><br>
<mailto:<a href="mailto:nek5000-users@lists" target="_blank">nek5000-users@lists</a>.<br>
<mailto:<a href="mailto:nek5000-users@lists" target="_blank">nek5000-users@lists</a>.>_<u></u>___<a href="http://mcs.anl.gov" target="_blank">mcs.anl.gov</a> <<a href="http://mcs.anl.gov" target="_blank">http://mcs.anl.gov</a>><div>
<div class="h5"><br>
<br>
<mailto:<a href="mailto:nek5000-users@lists." target="_blank">nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.<u></u>mcs.anl.gov</a>>>>> wrote:<br>
<br>
Hello,<br>
I am testing the example turbJet, and it seems<br>
that I have some<br>
problem with output. This are the output files:<br>
jet0.f00001<br>
<br>
Moreover, the ran stops with emergency exit at<br>
time 5.64<br>
(CFL too<br>
large?):<br>
<br>
<br>
282 DNORM, DIVEX 6.2789891850172769E-004<br>
8.4967723454134898E-006<br>
282 5.6400E+00 1.1065E+01 Fluid done<br>
282 5.6400000E+00 1.1955953E-03 6.4269380E-04<br>
9.9999991E-01 e2<br>
t max/min: 0.0000000000000000<br>
0.0000000000000000<br>
CFL, Ctarg! 22.504301079496056<br>
2.0000000000000000<br>
call outfld: ifpsco: T<br>
<br>
283 5.6400E+00 Write checkpoint:<br>
0 283 OPEN: jet0.f00001<br>
<br>
283 5.6400E+00 done :: Write checkpoint<br>
file size =<br>
55. MB<br>
avg data-throughput =<br>
572.1MB/s<br>
io-nodes = 1<br>
<br>
<br>
7 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
3 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
5 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
0 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
6 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
4 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
2 Emergency exit: 283 time =<br>
5.6399999999999677<br>
Latest solution and data are dumped for<br>
post-processing.<br>
*** STOP ***<br>
1 Emergency exit: 283 time =<br>
5.6399999999999677<br>
<br>
<br>
<br>
<br>
Thank you,<br>
Barak<br>
<br>
<br></div></div>
______________________________<u></u>_______________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.__<u></u>____gov<br>
<mailto:<a href="mailto:Nek5000-users@lists" target="_blank">Nek5000-users@lists</a>.<br>
<mailto:<a href="mailto:Nek5000-users@lists" target="_blank">Nek5000-users@lists</a>.>_<u></u>___<a href="http://mcs.anl.gov" target="_blank">mcs.anl.gov</a> <<a href="http://mcs.anl.gov" target="_blank">http://mcs.anl.gov</a>><br>
<mailto:<a href="mailto:Nek5000-users@lists." target="_blank">Nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.<u></u>mcs.anl.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-<u></u>users</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-<u></u>users</a>><br>
<br>
<<a href="https://lists.mcs.anl.gov/____mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/___<u></u>_mailman/listinfo/nek5000-<u></u>users</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><u></u>>__><div><div class="h5"><br>
<br>
<<a href="https://lists.mcs.anl.gov/____mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/___<u></u>_mailman/listinfo/nek5000-<u></u>users</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><u></u>><br>
<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>
<<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><u></u>>__>__><br>
<br>
<br>
<br>
<br>
<br>
<br>
______________________________<u></u>_____________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.__<u></u>__gov<br>
<mailto:<a href="mailto:Nek5000-users@lists." target="_blank">Nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.<u></u>mcs.anl.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>
<<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><u></u>><br>
<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>
<<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><u></u>>__><br>
<br>
<br>
______________________________<u></u>_____________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.__<u></u>__gov<br>
<mailto:<a href="mailto:Nek5000-users@lists." target="_blank">Nek5000-users@lists.</a>__<a href="http://mcs.anl.gov" target="_blank"><u></u>mcs.anl.gov</a><br>
<mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.<u></u>mcs.anl.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>
<<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><u></u>><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>
<<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><u></u>>__><br>
<br>
<br>
<br>
<br>
______________________________<u></u>___________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.__<u></u>gov<br>
<mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.<u></u>mcs.anl.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>
<<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><u></u>><br>
<br>
<br>
______________________________<u></u>___________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.__<u></u>gov<br>
<mailto:<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.<u></u>mcs.anl.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>
<<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><u></u>><br>
<br>
<br>
<br>
<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>
</div></div></blockquote>
<br>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a 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></div><br></div>