<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Paul and Stephan<div class=""><br class=""></div><div class="">how does the blow up manifests itself : I didnt dig a lot into the blow up. the flow converges nicely, and once its at steady state I inject the scalar at the inlet through the boundary condition (localized source). it develops a nice turbulent plume, that keeps going for a long time and then suddenly it produces lots of NaN. </div><div class=""><br class=""></div><div class="">this picture is consistent with the idea that sometimes passive scalars develop fronts, where abrupt changes occur across a small scale (bachelor scale; for me ~ kolmogorov scale cause diffusivity = viscosity). the classic Courant condition for advection will miss the presence of fronts, cause these are completely uncorrelated to the velocity field. (see eg PRL 84:2385, 2000). </div><div class=""><br class=""></div><div class="">i am thinking of including a further check to account for the presence of fronts. but i am sure others have faced similar problems, and so i wanted to check whether someone has already implemented something. </div><div class=""><br class=""></div><div class="">in the meantime, i’m using courant 0.25 and see if this is enough, although its clearly suboptimal. </div><div class=""><br class=""></div><div class="">cheers</div><div class="">Agnese</div><div class=""><br class=""></div><div class="">ps. Paul : filtering was already on.</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class=""><br class=""></div><div class=""><div dir="ltr" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div id="x_divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; font-family: Calibri, Helvetica, sans-serif;" class=""><div style="margin-top: 0px; margin-bottom: 0px;" class="">Have you tried reducing dt by, say, 2x ?</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">That should be enough...</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Also, you can try filtering by setting p103=0.05</div><div style="margin-top: 0px; margin-bottom: 0px;" class="">in the .rea file.</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">The velocity and passive scalars have the same CFL restriction, but velocity is regularized by pressure so it tends to be more stable than high Peclet passive scalar simulations.</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">hth,</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Paul</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div></div><hr tabindex="-1" style="display: inline-block; width: 640.90625px;" class=""><div id="x_divRplyFwdMsg" dir="ltr" class=""><font face="Calibri, sans-serif" style="font-size: 11pt;" class=""><b class="">From:</b><span class="Apple-converted-space"> </span><a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a><span class="Apple-converted-space"> </span><<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a>> on behalf of<span class="Apple-converted-space"> </span><a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><span class="Apple-converted-space"> </span><<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>><br class=""><b class="">Sent:</b><span class="Apple-converted-space"> </span>Thursday, June 15, 2017 3:38:03 AM<br class=""><b class="">To:</b><span class="Apple-converted-space"> </span><a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class=""><b class="">Subject:</b><span class="Apple-converted-space"> </span>Re: [Nek5000-users] Courant number for passive scalar</font><div class=""> </div></div></div><font size="2" style="font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-size: 10pt;" class="">How does the blow-up manifest itself?<br class=""><br class="">Try to run it with 2nd order in time and keep the CFL<0.5. This should be stable assuming there is no source term introducing an additional time step constraint.<br class=""><br class="">Does your mesh resolve these steep gradients? <br class=""><br class="">Cheers,<br class="">Stefan<br class=""> <br class="">-----Original message-----<br class="">> From:nek5000-<a href="mailto:users@lists.mcs.anl.gov" class="">users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>><br class="">> Sent: Thursday 15th June 2017 9:29<br class="">> To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="">> Subject: Re: [Nek5000-users] Courant number for passive scalar<br class="">> <br class="">> I havent seen the equivalent of a Courant criterion for passive scalars. <br class="">> what I meant is the scalar is intermittent, so one should have a way to detect the regions where gradients are extremely large (fronts) and make sure to decrease the timestep there.<br class="">> <br class="" />fronts do not correspond to regions where velocities are especially large, so the classic courant criterion will miss them. <br class="">> <br class="" />does this make sense ?<br class="">> any suggestions ?<br class="">> <br class="" />agnese<br class="">> <br class="" /><br class="" /><br class="" />On Jun 14, 2017, at 2:39 PM, <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">mailto:nek5000-users-request@lists.mcs.anl.gov</a>> wrote:<br class="">> <br class="Apple-interchange-newline" />Send Nek5000-users mailing list submissions to<br class="" /> <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">mailto:nek5000-users@lists.mcs.anl.gov</a>><br class="" /><br class="" />To subscribe or unsubscribe via the World Wide Web, visit<br class="" /> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users%3Cbr" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br</a> class="" />or, via email, send a message with subject or body 'help' to<br class="" /> <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a><br class="" /><br class="" />You can reach the person managing the list at<br class="" /> <a href="mailto:nek5000-users-owner@lists.mcs.anl.gov" class="">nek5000-users-owner@lists.mcs.anl.gov</a><br class="" /><br class="" />When replying, please edit your Subject line so it is more specific<br class="" />than "Re: Contents of Nek5000-users digest..."<br class="" />Today's Topics:<br class="" /><br class="" /> 1. Re: Courant for passive scalar ? (<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>)<br class="" /> 2. Re: Courant for passiv<br class=""> e scalar<br class=""> ? (<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>)<br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: Re: [Nek5000-users] Courant for passive scalar ?<br class="" />Date: June 14, 2017 at 10:44:13 AM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />Thanks Jan !<br class="">> <br class="" />so it is confirmed that the tilmestep is based exclusively on the velocity field ?<br class="">> how do people make sure that the passive scalars do not blow up ? <br class="">> anybody ?<br class="">> <br class="" />any help is greatly appreciated !!<br class="">> agnese<br class="">> <br class="" /><br class="" />Agnese Seminara<br class="">> --------------------------------<br class="">> CNRS<br class="">> Institut de Physique de Nice<br class="">> Parc Valrose <br class="">> avenue J Vallot<br class="">> 06108 Nice, France<br class="">> +33 (0) 492 076 775<br class="">> <a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a> <<a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a>><br class="" />On Jun 13, 2017, at 7:00 PM, <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">mailto:nek5000-users-request@lists.mcs.anl.gov</a>> wrote:<br class="">> <br class="Apple-interchange-newline" />Send Nek5000-users mailing list submissions to<br class="" /> <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">mailto:nek5000-users@lists.mcs.anl.gov</a>><br class="" /><br class="" />To subscribe or unsubscribe via the World Wide Web, visit<br class="" /> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users%3Cbr" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br</a> class="" />or, via email, send a message with subject or body 'help' to<br class="" /> <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a><br class="" /><br class="" />You can reach the person managing the list at<br class="" /> <a href="mailto:nek5000-users-owner@lists.mcs.anl.gov" class="">nek5000-users-owner@lists.mcs.anl.gov</a><br class="" /><br class="" />When replying, please edit your Subject line so it is more specific<br class="" />than "Re: Contents of Nek5000-users digest..."<br class="" />Today's Topics:<br class="" /><br class="" /> 1. Re: Courant for passive scalar ? (<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>)<br class="" /> 2. prenek (nek5000-users@<br class=""> lists.mc<br class=""> <a href="http://s.anl.gov/" class="">s.anl.gov</a>)<br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: Re: [Nek5000-users] Courant for passive scalar ?<br class="" />Date: June 12, 2017 at 8:00:26 PM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />Hi,<br class="">> I haven' t tried this myself, but I guess you could write a subroutine that calculates the maximum timestep from the scalar field and sets dt accordingly. Then simple call the routine from userchk.<br class="">> <br class="" />Best,<br class="">> Jan<br class="">> <br class="" />Am 12.06.2017 um 18:15 schrieb <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">mailto:nek5000-users@lists.mcs.anl.gov</a>>:<br class="" /><br class="" />Dear Neks<br class="">> <br class="" />I am running a simulation of a channel, with a passive scalar. Using diffusivity=viscosity, the velocity field converges nicely, whereas the passive scalar blows up at some random point down the simulation. Reynolds number is 10 000 so its in the turbulent regime.<br class="">> <br class="" />I suspect this is due to the fact that the passive scalar is intermittent, so sometimes large fluctuations will occur, totally uncorrelated to velocity fluctuations. to avoid the code to blow up, one would need to setup the timestep based on a courant criterion for the scalar, and not for the velocity. <br class="">> <br class="" />has this already been implemented ?<br class="">> or how do people avoid this ?<br class="">> <br class="" />any suggestion is more than welcome.<br class="">> thanks<br class="">> Agnese<br class="">> <br class="" /><br class="" /><br class="" />Agnese Seminara<br class="">> --------------------------------<br class="">> CNRS<br class="">> Institut de Physique de Nice<br class="">> Parc Valrose <br class="">> avenue J Vallot<br class="">> 06108 Nice, France<br class="">> +33 (0) 492 076 775<br class="">> <a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a> <<a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a>><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">mailto:Nek5000-users@lists.mcs.anl.gov</a>><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>><br class="" /><br class="" /><br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: [Nek5000-users] prenek<br class="" />Date: June 13, 2017 at 12:12:49 PM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />Hi nek,<br class="">> <br class="" />I have a 2D mesh. When I used global refine function, the errors occurred. do you know why?<br class="">> <br class="" /> OCT/Multi-SPLIT<br class="">> quad or multi-split? (q/m)<br class="">> Number of elements after OctSplit operation: 7680<br class="">> Are you sure you want to split (y/n)?<br class="">> inside gencen 7680 F 4<br class="">> done gencen 7680<br class="">> start locglob_lexico: 4 7680 30720 2.00000009E-03<br class="">> locglob: 1 1 30720<br class="">> locglob: 2 3664 30720<br class="">> locglob: 1 7761 30720<br class="">> locglob: 2 7761 30720<br class="">> Performing unique_vertex2 self_chk 30720<br class="">> done locglob_lexico: 7761 7761 30720 4<br class="">> Performing makecell self_chk 30720<br class="">> <br class="" /> RMIN: 0.153545E-03<br class="">> EPS,RMIN: 0.500000E-01 0.153545E-03<br class="">> I/O Error: No String Terminator sent to PRS<br class="">> zeroing y: 1 82 -0.75715434E-10<br class="">> I/O Error: No String Terminator sent to PRS<br class="">> zeroing y: 4 82 -0.14778734E-11<br class="">> I/O Error: No String Terminator sent to PRS<br class="">> zeroing y: 4 83 0.72759583E-10<br class="">> I/O Error: No String Terminator sent to PRS<br class="">> zeroing y: 4 123 0.15279482E-09<br class="">> ......<br class="">> <br class="" />Does it mean my mesh is still not right?<br class="">> <br class="" /><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br class="" /><br class="" /><br class="" /><br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: Re: [Nek5000-users] Courant for passive scalar ?<br class="" />Date: June 14, 2017 at 2:37:35 PM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />How do you define a courant number for a scalar?<br class="" /><br class="" />-----Original message-----<br class="" />From:nek5000-<a href="mailto:users@lists.mcs.anl.gov" class="">users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>><br class="" />Sent: Wednesday 14th June 2017 10:43<br class="" />To: <a href="mailto:nek5000-users@lists.mcs" class="">nek5000-users@lists.mcs</a><br class=""> .<a href="http://anl.gov/" class="">anl.gov</a><br class=""> <br class="" />Subject: Re: [Nek5000-users] Courant for passive scalar ?<br class="" /><br class="" />Thanks Jan !<br class="" /><br class="" />so it is confirmed that the tilmestep is based exclusively on the velocity field ?<br class="" />how do people make sure that the passive scalars do not blow up ? <br class="" />anybody ?<br class="" /><br class="" />any help is greatly appreciated !!<br class="" />agnese<br class="" /><br class="" /><br class="" />Agnese Seminara<br class="" />--------------------------------<br class="" />CNRS<br class="" />Institut de Physique de Nice<br class="" />Parc Valrose <br class="" />avenue J Vallot<br class="" />06108 Nice, France<br class="" />+33 (0) 492 076 775<br class="" /><a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a><<a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a>><br class="" />On Jun 13, 2017, at 7:00 PM, <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a> <<a href="" class=""></a><a href="mailto:nek5000-users-request@lists.mcs.a" class="">mailto:nek5000-users-request@lists.mcs.a</a><br class=""> <a href="http://nl.gov/" class="">nl.gov</a>> <br class=""> wrote:<br class="" /><br class="Apple-interchange-newline" />Send Nek5000-users mailing list submissions to<br class="" /> <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">mailto:nek5000-users@lists.mcs.anl.gov</a>><br class="" /><br class="" />To subscribe or unsubscribe via the World Wide Web, visit<br class="" /> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users%3Cbr" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br</a> class="" />or, via email, send a message with subject or body 'help' to<br class="" /> <a href="mailto:nek5000-users-request@lists.mcs.anl.gov" class="">nek5000-users-request@lists.mcs.anl.gov</a><br class="" /><br class="" />You can reach the person managing the list at<br class="" /> <a href="mailto:nek5000-users-owner@lists.mcs.anl.gov" class="">nek5000-users-owner@lists.mcs.anl.gov</a><br class="" /><br class="" />When replying, please edit your Subject line so it is more specific<br class="" />than "Re: Contents of Nek5000-users digest..."<br class="" />Today's Topics:<br class="" /><br class="" /> 1. Re: Courant for passive scalar ? (<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>)<br class="" /> 2. <br class=""> prenek (<br class=""> nek5000-users@<br class="" /> lists.mc<br class="" /> <a href="http://s.anl.gov/" class="">s.anl.gov</a>)<br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: Re: [Nek5000-users] Courant for passive scalar ?<br class="" />Date: June 12, 2017 at 8:00:26 PM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />Hi,<br class="" />I haven' t tried this myself, but I guess you could write a subroutine that calculates the maximum timestep from the scalar field and sets dt accordingly. Then simple call the routine from userchk.<br class="" /><br class="" />Best,<br class="" />Jan<br class="" /><br class="" />Am 12.06.2017 um 18:15 schrieb <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">mailto:nek5000-users@lists.mcs.anl.gov</a>>:<br class="" /><br class="" />Dear Neks<br class="" /><br class="" />I am running a simulation of a channel, with a passive scalar. Usi<br class=""> ng diffu<br class=""> sivity=viscosity, the velocity field converges nicely, whereas the passive scalar blows up at some random point down the simulation. Reynolds number is 10 000 so its in the turbulent regime.<br class="" /><br class="" />I suspect this is due to the fact that the passive scalar is intermittent, so sometimes large fluctuations will occur, totally uncorrelated to velocity fluctuations. to avoid the code to blow up, one would need to setup the timestep based on a courant criterion for the scalar, and not for the velocity. <br class="" /><br class="" />has this already been implemented ?<br class="" />or how do people avoid this ?<br class="" /><br class="" />any suggestion is more than welcome.<br class="" />thanks<br class="" />Agnese<br class="" /><br class="" /><br class="" /><br class="" />Agnese Seminara<br class="" />--------------------------------<br class="" />CNRS<br class="" />Institut de Physique de Nice<br class="" />Parc<br class=""> Valrose<br class=""> <br class="" />avenue J Vallot<br class="" />06108 Nice, France<br class="" />+33 (0) 492 076 775<br class="" /><a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a> <<a href="http://sites.unice.fr/site/aseminara/" class="">http://sites.unice.fr/site/aseminara/</a>><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">mailto:Nek5000-users@lists.mcs.anl.gov</a>><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a> <<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>><br class="" /><br class="" /><br class="" /><br class="" />From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Subject: [Nek5000-users] prenek<br class="" />Date: June 13, 2017 at 12:12:49 PM GMT+2<br class="" />To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" />Reply-To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="" /><br class="" /><br class="" />Hi nek,<br class="" /><br class="" />I have a 2D mesh. When I used global refine fu<br class=""> nction, <br class=""> the errors occurred. do you know why?<br class="" /><br class="" /> OCT/Multi-SPLIT<br class="" />quad or multi-split? (q/m)<br class="" />Number of elements after OctSplit operation: 7680<br class="" />Are you sure you want to split (y/n)?<br class="" /> inside gencen 7680 F 4<br class="" /> done gencen 7680<br class="" /> start locglob_lexico: 4 7680 30720 2.00000009E-03<br class="" /> locglob: 1 1 30720<br class="" /> locglob: 2 3664 30720<br class="" /> locglob: 1 7761 � A0 30720<br class=""> <br class="" /> locglob: 2 7761 30720<br class="" /> Performing unique_vertex2 self_chk 30720<br class="" /> done locglob_lexico: 7761 7761 30720 4<br class="" /> Performing makecell self_chk 30720<br class="" /><br class="" /> RMIN: 0.153545E-03<br class="" /> EPS,RMIN: 0.500000E-01 0.153545E-03<br class="" /> I/O Error: No String Terminator sent to PRS<br class="" /> zeroing y: 1 82 -0.75715434E-10<br class="" /> I/O Error: No String Terminator sent to PRS<br class="" /> zeroing y: 4 82 -0.14778734E-11<br class="" /> I/O Error: No String Terminator sent to PRS<br class="" /> zeroing y: 4 83 0.72759583E-10<br class=""<br class=""> />�=<br class=""> 0I/O Error: No String Terminator sent to PRS<br class="" /> zeroing y: 4 123 0.15279482E-09<br class="" />......<br class="" /><br class="" />Does it mean my mesh is still not right?<br class="" /><br class="" /><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br class="" /><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br class="" /><br class="" /><br class="" /><br class="" />_______________________________________________<br class="" />Nek5000-users mailing list<br class="" /><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="" /><a href="https://lists.mcs.anl.gov/mailman/listinfo" class="">https://lists.mcs.anl.gov/mailman/listinfo</a><br class=""> /nek5000<br class=""> -users<br class="" /><br class="" />_______________________________________________<br class="">> Nek5000-users mailing list<br class="">> <a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="">> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br class=""></span></font><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class=""></span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">_______________________________________________</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Nek5000-users mailing list</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="mailto:Nek5000-users@lists.mcs.anl.gov" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Nek5000-users@lists.mcs.anl.gov</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a></div></blockquote></div><br class=""></div></body></html>