[Nek5000-users] turbJet problem
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sun Dec 8 18:56:56 CST 2013
Hi Barak,
I am not aware of variable time stepping in nek5000. The time stepping used
is BDFk/EXTk with CFL constraint for stability. You can change TORDER in
.rea file to 2 or 3 to make the time stepping 2nd or third order accurate.
There is a characteristic based time-stepping though... IFCHAR in .rea file.
However, if for a value ot DT you are getting convergence, essentially you
should get the same for smaller value of DT. What is the error you are
getting after the program aborts after 1200 iterations? It is still stored
in your log file I suppose.
Thanks,
Tanmoy
On Sun, Dec 8, 2013 at 12:26 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
> Hi Tanmoy,
> I have divided by 4 and again my run stopped at iteration 1200. My
> question is: does nek5000 have variable time-stepping?
> Thanks,
> Barak
>
>
> On 12/07/2013 08:37 PM, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hi Barak,
>>
>> it is because of the fact that the generated file format is probably not
>> readable in visit. In the *usrchk* subroutine in your .usr file write
>>
>> these 2 lines
>> param(66) = 4
>> param(67) = 4
>>
>> before the end of the routine
>>
>> Then compile and execute your program. Now you will get output files in
>> the form of turbJet.fld001, turbJet.fld002, ...and so on.
>> type in visnek turbJet in your terminal, and generate *turbJet.nek5000*.
>>
>> This .nek5000 file can be read by visit.
>>
>> Thanks,
>> Tanmoy
>>
>>
>> On Sat, Dec 7, 2013 at 10:33 AM, <nek5000-users at lists.mcs.anl.gov
>> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
>>
>> Hi Tanmoy,
>> Many thanks. It seems that after dividing by 4 (dt=0.005) it runs
>> and I get output files in the form jet0.fxxxxx .
>> Now I am trying to visualize them using the precompiled ubuntu11
>> version of visit. Unfortunately, I get the following error:
>>
>> VisIt could not read from the file
>> "/home/galanti/COMPUTATION/__NEK5000/nek5_svn/TEST-__
>> examples/turbJet/jet0.f*
>>
>> database".
>>
>> The generated error message was:
>>
>> There was an error opening
>> /home/galanti/COMPUTATION/__NEK5000/nek5_svn/TEST-__
>> examples/turbJet/jet0.f*
>>
>> database. It may be an invalid file. VisIt tried using the
>> following file format readers to open the file: Nek5000
>>
>> The following error(s) may be helpful in identifying the problem:
>> When attempting to use "Nek", the file matched the general format
>> expected by the reader but eventually encountered this error: "A tag
>> called filetemplate: must be specified".
>>
>> What's wrong?
>> Many thanks again,
>> Barak
>>
>>
>>
>> On 12/07/2013 12:42 AM, nek5000-users at lists.mcs.anl.__gov
>>
>> <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>>
>> Hi Barak,
>>
>> It is a problem with numerical Stability. Decrease DT in your
>> .rea file
>> until the simulations become stable. You could start with 1/2 DT
>> to
>> check if simulation is stable, other wise go for 1/4 DT, and
>> check with
>> this simple test.
>>
>> Thanks,
>> Tanmoy
>>
>>
>> On Fri, Dec 6, 2013 at 3:31 PM,
>> <nek5000-users at lists.mcs.anl.__gov
>> <mailto:nek5000-users at lists.mcs.anl.gov>
>> <mailto:nek5000-users at lists.__mcs.anl.gov
>>
>> <mailto:nek5000-users at lists.mcs.anl.gov>>> wrote:
>>
>> Hello,
>> I am testing the example turbJet, and it seems that I have
>> some
>> problem with output. This are the output files:
>> jet0.f00001
>>
>> Moreover, the ran stops with emergency exit at time 5.64
>> (CFL too
>> large?):
>>
>>
>> 282 DNORM, DIVEX 6.2789891850172769E-004
>> 8.4967723454134898E-006
>> 282 5.6400E+00 1.1065E+01 Fluid done
>> 282 5.6400000E+00 1.1955953E-03 6.4269380E-04
>> 9.9999991E-01 e2
>> t max/min: 0.0000000000000000 0.0000000000000000
>> CFL, Ctarg! 22.504301079496056 2.0000000000000000
>> call outfld: ifpsco: T
>>
>> 283 5.6400E+00 Write checkpoint:
>> 0 283 OPEN: jet0.f00001
>>
>> 283 5.6400E+00 done :: Write checkpoint
>> file size = 55. MB
>> avg data-throughput =
>> 572.1MB/s
>> io-nodes = 1
>>
>>
>> 7 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 3 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 5 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 0 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 6 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 4 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 2 Emergency exit: 283 time =
>> 5.6399999999999677
>> Latest solution and data are dumped for post-processing.
>> *** STOP ***
>> 1 Emergency exit: 283 time =
>> 5.6399999999999677
>>
>>
>>
>>
>> Thank you,
>> Barak
>>
>>
>> ___________________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.____gov
>> <mailto:Nek5000-users at lists.__mcs.anl.gov
>> <mailto:Nek5000-users at lists.mcs.anl.gov>>
>> https://lists.mcs.anl.gov/____mailman/listinfo/nek5000-users
>> <https://lists.mcs.anl.gov/__mailman/listinfo/nek5000-users>
>> <https://lists.mcs.anl.gov/__mailman/listinfo/nek5000-users
>> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>__>
>>
>>
>>
>>
>>
>>
>> _________________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.__gov
>> <mailto:Nek5000-users at lists.mcs.anl.gov>
>> https://lists.mcs.anl.gov/__mailman/listinfo/nek5000-users
>> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>
>>
>> _________________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.__gov
>> <mailto:Nek5000-users at lists.mcs.anl.gov>
>> https://lists.mcs.anl.gov/__mailman/listinfo/nek5000-users
>> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>
>>
>>
>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131208/05a1218d/attachment-0001.html>
More information about the Nek5000-users
mailing list