Re: [Nek5000-users] Calculating very small Lewis numbers in Nekton
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 <[email protected]> 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 <[email protected]> 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, <[email protected]> wrote:
Can you please post your logfile
On 4 Dec 2017, at 21:08, "[email protected]" < [email protected]> 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, <[email protected]> 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, [email protected] 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 [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________
Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
participants (1)
-
nek5000-users＠lists.mcs.anl.gov