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@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@lists.mcs.anl.gov> wrote:

Can you please post your logfile 

On 4 Dec 2017, at 21:08, "nek5000-users@lists.mcs.anl.gov" <nek5000-users@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@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@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@lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

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

_______________________________________________

Nek5000-users mailing list

Nek5000-users@lists.mcs.anl.gov

https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

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