[Nek5000-users] errors

Mon Nov 23 17:13:21 CST 2015

Hi,

The first error was because it did not converge to the level you set (parameter 12) when solving for x velocity. You may check your initial condition or other set up.

As for the second error, you would have a *.sch file every time you run nek, but there is a check to see whether the file exists. You may just delete the *.sch file before you run nek again.

Cheers,
Tony

________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users-request at lists.mcs.anl.gov <nek5000-users-request at lists.mcs.anl.gov>
Sent: 23 November 2015 18:00
To: nek5000-users at lists.mcs.anl.gov
Subject: Nek5000-users Digest, Vol 81, Issue 21

Send Nek5000-users mailing list submissions to
        nek5000-users at lists.mcs.anl.gov

To subscribe or unsubscribe via the World Wide Web, visit
        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
        nek5000-users-request at lists.mcs.anl.gov

You can reach the person managing the list at
        nek5000-users-owner at lists.mcs.anl.gov

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Nek5000-users digest..."

Today's Topics:

   1. Re: low order interpolation (nek5000-users at lists.mcs.anl.gov)
   2. errors (nek5000-users at lists.mcs.anl.gov)

----------------------------------------------------------------------

Message: 1
Date: Sun, 22 Nov 2015 19:44:10 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] low order interpolation
Message-ID:
        <mailman.11018.1448217871.24727.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="utf-8"

Hi Tanmoy

I really don't understand why the interpolation fails.. I hope that someone
can help me because a sample low order interpolation could solve my problem.

Thanks anyway
Alessandro

2015-11-21 23:53 GMT+01:00 <nek5000-users at lists.mcs.anl.gov>:

> Hi Alessandro,
>
> Ok. But, I was not talking about your matlab routines, but rather I was
> talking about your nek5000 postprocessing mode. I guess they are correct,
> since you got good results with Tecplot.
>
> But anyway, the bottom line is, nek5000 uses spectral interpolation
> routine, so even the higher order polynomials are orthogonal, so it is
> unlikely that you will have Runge's phenomenon.
>
> Also, I myself use hpts routines for interpolating a turbulent field
>
> Best Regards,
> Tanmoy
>
> On Fri, Nov 20, 2015 at 10:41 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>> Hi Tanmoy
>>
>> thanks for the quick answer, unfortunately I didn't commit this kinds of
>> errors. I have already used my matlab routine for other cases and it worked
>> well, so I'm pretty sure that there isn't errors. The only differences is
>> that I'm using this routine to interpolate the hpts.in points on more
>> complex grid and on a turbulent field. It seams that I have the Runge's
>> phenomenon. Is this possible with the hpts routine?
>>
>> Sincerely
>> Alessandro
>>
>>
>>
>>
>>
>> 2015-11-20 17:29 GMT+01:00 <nek5000-users at lists.mcs.anl.gov>:
>>
>>> Hi Alessandro,
>>>
>>> I use exactly the same procedure for my case. Just wondering, by any
>>> chance, are there any points in hpts.in which are outside bounds of the
>>> original domain, or you may resize the domain in userdat(), while
>>> postprocessing in nek5000.
>>>
>>> Best Regards,
>>> Tanmoy
>>>
>>> On Fri, Nov 20, 2015 at 6:46 AM, <nek5000-users at lists.mcs.anl.gov>
>>> wrote:
>>>
>>>> Dear neks
>>>>
>>>> I don't know how Nek5000 performs the interpolation with hpts() routine
>>>> but I have some problem when I use this last. Using Matlab I built an
>>>> hpts.in file that contains the points of a structured and not constant
>>>> mesh, then using Nek in postprocessing  mode I evaluate the values of
>>>> velocity an pressure in this points. When I reshape the hpts.out in order
>>>> to visualize my domain, there are many zones that have no-sense values!. I
>>>> performed the same interpolation in tecplot and it works fine.
>>>>
>>>> my question is: there's any way to perform a first order interpolation
>>>> in Nek5000?
>>>>
>>>> Sincerely
>>>> Alessandro
>>>>
>>>> _______________________________________________
>>>> 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
>>>
>>>
>>
>> _______________________________________________
>> 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/20151122/ffe78724/attachment-0001.html>

------------------------------

Message: 2
Date: Mon, 23 Nov 2015 15:32:46 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] errors
Message-ID:
        <mailman.11034.1448289207.24727.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset=utf-8; format=flowed

Hi Neks,
can someone tell me what means these errors?

  1   100 **ERROR**: Failed in HMHOLTZ: VELX   7.5730E+00 8.9997E+04
1.0000E-08
   1.00000000000000002E-008  p22            1           1
   1     1 Helmholtz VELY    F:   1.1292E-11   1.0000E-08 1.0000E-03
1.0000E+02
           1    Hmholtz VELY:      0   1.1292E-11   1.1292E-11 1.0000E-08
   1.00000000000000002E-008  p22            1           1
   1     1 Helmholtz VELZ    F:   6.8078E-12   1.0000E-08 1.0000E-03
1.0000E+02
           1    Hmholtz VELZ:      0   6.8078E-12   6.8078E-12 1.0000E-08
                L1/L2 DIV(V)    :   2.3607E-17   2.1252E-05
                L1/L2 QTL       :   0.0000E+00   0.0000E+00
                L1/L2 DIV(V)-QTL:   2.3607E-17   2.1252E-05
           1   1.0000E-02  1.0853E+01 Fluid done

and:

ERROR: .sch file already exists. Use IFSCHCLOB=F to disable this check
BUT BEWARE!!!!!!  ierr=  1

They are came out at first step of my sim.

Thanks

------------------------------

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

End of Nek5000-users Digest, Vol 81, Issue 21
*********************************************