[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