[Nek5000-users] Calculating very small Lewis numbers in Nekton

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Dec 5 13:31:27 CST 2017


Hello all,

Thanks for the replies, I think my log-file was too big even after
compression and stuck in moderator approval for the list. But here is the
part of the log-file where the problem starts. I have attached a google
drive link which has the log file and .rea file, the Lewis Number is set in
line 133 of the .rea file after parameter 128 (p128). Link -
https://drive.google.com/open?id=1ozsy0IgVKuxDiplGl6CfYDn3GuhDWoWI


Solving for fluid F T T
       2071    Hmholtz TEMP:      1   1.1353E+00   5.5267E+00   7.0000E+00
       2071    Hmholtz PS 1:     13   5.5318E+00   3.8255E+21   7.0000E+00
  2071   100 **ERROR**: Failed in HMHOLTZ: PS 2   6.5359E+15   3.8255E+21
 7.0000E+00
       2071    Hmholtz PS 3:     13   5.5318E+00   3.8255E+21   7.0000E+00
       2071   1.3076E+03  2.0528E+00 Heat done
       2071    Hmholtz VELX:     38   9.8896E-11   4.9469E+02   1.0000E-10
       2071    Hmholtz VELY:     38   9.2808E-11   4.4808E+02   1.0000E-10
       2071    Hmholtz VELZ:     38   8.0264E-11   4.6649E+02   1.0000E-10
        2071 reorthogo:           0
 2071    6 alpha:  9.1687E-01  4.2049E-02  6.7213E-03  2.4412E-03
1.0813E-03  1.7375E-03
  2071   6  1.7386E-01  1.3111E-03  1.3260E+02 alph12
      2071 U-PRES gmres:      5  7.5724E-09  1.0000E-08  4.0064E-08
7.3007E-02  3.7297E-01
        2071  DNORM, DIVEX  5.430060732380546E-007  7.572385434962068E-009
       2071   1.3076E+03  2.4498E+00 Fluid done
  2071 1 flx  1.307563E+03  5.633556E-02  1.121772E+04  1.174017E+04
5.092129E+02
  2071 2 flx  1.307563E+03  5.633556E-02 -1.122518E+04 -1.174681E+04
4.949962E+02
Step   2072, t= 1.3075634E+03, DT= 5.6082644E-05, C=  0.500 6.5792E+03
5.2920E+00
             Solving for heat
 Solving for fluid F T T
       2072    Hmholtz TEMP:      1   1.1351E+00   5.5225E+00   7.0000E+00
  2072   100 **ERROR**: Failed in HMHOLTZ: PS 1   3.2982E+26   3.0368E+40
 7.0000E+00
  2072   100 **ERROR**: Failed in HMHOLTZ: PS 2   4.0048E+36   3.0365E+40
 7.0000E+00
  2072   100 **ERROR**: Failed in HMHOLTZ: PS 3   3.2982E+26   3.0368E+40
 7.0000E+00
       2072   1.3076E+03  4.8085E+00 Heat done
       2072    Hmholtz VELX:     38   9.9602E-11   4.9398E+02   1.0000E-10
       2072    Hmholtz VELY:     38   9.3511E-11   4.4755E+02   1.0000E-10
       2072    Hmholtz VELZ:     38   8.0881E-11   4.6601E+02   1.0000E-10
        2072 reorthogo:           0
 2072    7 alpha:  9.1610E-01  4.8479E-02  7.7574E-03  5.1508E-03
1.8164E-03 -4.1065E-04 -4.3921E-04
  2072   7  1.7378E-01  1.1572E-03  1.5018E+02 alph12
      2072 U-PRES gmres:      4  8.4386E-09  1.0000E-08  3.5390E-08
4.9622E-02  3.2245E-01
        2072  DNORM, DIVEX  5.419045247755792E-007  8.438599514879914E-009
       2072   1.3076E+03  2.5401E+00 Fluid done

Thanks,
Saikat

Saikat Mukherjee,
PhD Student,
Paul Research Group - http://www.me.vt.edu/mpaul/
Engineering Science and Mechanics,
Virginia Tech.

On Tue, Dec 5, 2017 at 1:38 PM, Saikat Mukherjee <msaikat6 at vt.edu> wrote:

> Hello all,
>
> Thanks for the replies, I think my log-file was too big and stuck in
> approval for the list. But here is the part of the log-file where the
> problem starts. I have attached the .rea file, the Lewis Number is set in
> line 133 of the .rea file after parameter 128 (p128).
>
> Thanks,
> Saikat
>
>
>
>
> Solving for fluid F T T
>        2071    Hmholtz TEMP:      1   1.1353E+00   5.5267E+00   7.0000E+00
>        2071    Hmholtz PS 1:     13   5.5318E+00   3.8255E+21   7.0000E+00
>   2071   100 **ERROR**: Failed in HMHOLTZ: PS 2   6.5359E+15   3.8255E+21
>  7.0000E+00
>        2071    Hmholtz PS 3:     13   5.5318E+00   3.8255E+21   7.0000E+00
>        2071   1.3076E+03  2.0528E+00 Heat done
>        2071    Hmholtz VELX:     38   9.8896E-11   4.9469E+02   1.0000E-10
>        2071    Hmholtz VELY:     38   9.2808E-11   4.4808E+02   1.0000E-10
>        2071    Hmholtz VELZ:     38   8.0264E-11   4.6649E+02   1.0000E-10
>         2071 reorthogo:           0
>  2071    6 alpha:  9.1687E-01  4.2049E-02  6.7213E-03  2.4412E-03
> 1.0813E-03  1.7375E-03
>   2071   6  1.7386E-01  1.3111E-03  1.3260E+02 alph12
>       2071 U-PRES gmres:      5  7.5724E-09  1.0000E-08  4.0064E-08
> 7.3007E-02  3.7297E-01
>         2071  DNORM, DIVEX  5.430060732380546E-007  7.572385434962068E-009
>        2071   1.3076E+03  2.4498E+00 Fluid done
>   2071 1 flx  1.307563E+03  5.633556E-02  1.121772E+04  1.174017E+04
> 5.092129E+02
>   2071 2 flx  1.307563E+03  5.633556E-02 -1.122518E+04 -1.174681E+04
> 4.949962E+02
> Step   2072, t= 1.3075634E+03, DT= 5.6082644E-05, C=  0.500 6.5792E+03
> 5.2920E+00
>              Solving for heat
>  Solving for fluid F T T
>        2072    Hmholtz TEMP:      1   1.1351E+00   5.5225E+00   7.0000E+00
>   2072   100 **ERROR**: Failed in HMHOLTZ: PS 1   3.2982E+26   3.0368E+40
>  7.0000E+00
>   2072   100 **ERROR**: Failed in HMHOLTZ: PS 2   4.0048E+36   3.0365E+40
>  7.0000E+00
>   2072   100 **ERROR**: Failed in HMHOLTZ: PS 3   3.2982E+26   3.0368E+40
>  7.0000E+00
>        2072   1.3076E+03  4.8085E+00 Heat done
>        2072    Hmholtz VELX:     38   9.9602E-11   4.9398E+02   1.0000E-10
>        2072    Hmholtz VELY:     38   9.3511E-11   4.4755E+02   1.0000E-10
>        2072    Hmholtz VELZ:     38   8.0881E-11   4.6601E+02   1.0000E-10
>         2072 reorthogo:           0
>  2072    7 alpha:  9.1610E-01  4.8479E-02  7.7574E-03  5.1508E-03
> 1.8164E-03 -4.1065E-04 -4.3921E-04
>   2072   7  1.7378E-01  1.1572E-03  1.5018E+02 alph12
>       2072 U-PRES gmres:      4  8.4386E-09  1.0000E-08  3.5390E-08
> 4.9622E-02  3.2245E-01
>         2072  DNORM, DIVEX  5.419045247755792E-007  8.438599514879914E-009
>        2072   1.3076E+03  2.5401E+00 Fluid done
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> Saikat Mukherjee,
> PhD Student,
> Paul Research Group - http://www.me.vt.edu/mpaul/
> Engineering Science and Mechanics,
> Virginia Tech.
>
> On Mon, Dec 4, 2017 at 3:52 PM, Saikat Mukherjee <msaikat6 at vt.edu> wrote:
>
>> Hey all,
>>
>> Thanks for the replies. Attached are the log file and the .rea file where
>> I set my Lewis Number. The Lewis Number is set after parameter 128 (p128)
>> at line 133 of the .rea file.
>>
>> Saikat
>>
>> Saikat Mukherjee,
>> PhD Student,
>> Paul Research Group - http://www.me.vt.edu/mpaul/
>> Engineering Science and Mechanics,
>> Virginia Tech.
>>
>> On Mon, Dec 4, 2017 at 3:15 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>>
>>> Can you please post your logfile
>>>
>>> On 4 Dec 2017, at 21:08, "nek5000-users at lists.mcs.anl.gov" <
>>> nek5000-users at lists.mcs.anl.gov> wrote:
>>>
>>> Hey Marco,
>>>
>>> The error message is "Failed in HMHOLTZ" for the scalar and NaN value is
>>> being returned.
>>>
>>> Thanks,
>>> Saikat
>>>
>>> Saikat Mukherjee,
>>> PhD Student,
>>> Paul Research Group - http://www.me.vt.edu/mpaul/
>>> Engineering Science and Mechanics,
>>> Virginia Tech.
>>>
>>> On Mon, Dec 4, 2017 at 2:34 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>>>
>>>> Salkat,
>>>>
>>>> could you please be more specific when you write ‘ I have not been
>>>> successful’? Do you get an error message?
>>>>
>>>> Marco
>>>>
>>>> On Dec 4, 2017, at 2:31 PM, nek5000-users at lists.mcs.anl.gov wrote:
>>>>
>>>> Hello Neks,
>>>>
>>>> I am solving a Reaction Advection Diffusion equation, where with my
>>>> non-dimensionalization, the Lewis Number(Le) is the parameter multiplying
>>>> the Laplacian. That is, Le is the  non-dimensional Diffusion co-efficient
>>>> of the species. Currently I can do Le = 0.1, but to reduce it further down
>>>> to 0.01, or 0.001, the spatial resolution has to be increased. I tried
>>>> doing that increasing the GLL points in SIZE file by increasing the lx1,
>>>> lxd, but till now I have been unsuccessful.
>>>>
>>>> I am aware there is a special filtering that can be used( Fischer and
>>>> Mullen, 2001). My question is how to implement that in my .usr file. Is
>>>> there a routine?
>>>>
>>>> Thanks,
>>>> Saikat
>>>>
>>>>
>>>> Saikat Mukherjee,
>>>> PhD Student,
>>>> Paul Research Group - http://www.me.vt.edu/mpaul/
>>>> Engineering Science and Mechanics,
>>>> Virginia Tech.
>>>> _______________________________________________
>>>> 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/20171205/0142d859/attachment-0001.html>


More information about the Nek5000-users mailing list