OpenFOAM: FOAM FATAL IO ERROR - expected Scalar - nan - data.solverPerformance.e

[Foad S Farimani]

hello guys, I have asked this question on cfd-online and reddit too.

I kinda have a same issue as the OP so I post here too.

I'm following this tutorial for compressible fluid. except my geometry is little bit different:

But when I try to solve it with sonicFoam I get the error below:

Time = 0.002

Courant Number mean: 62.7429 max: 120

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

PIMPLE: iteration 1

smoothSolver:  Solving for Ux, Initial residual = 0.738285, Final residual = 0.357953, No Iterations 1000

smoothSolver:  Solving for Uy, Initial residual = 0.0131426, Final residual = 4.93472e-018, No Iterations 1

smoothSolver:  Solving for e, Initial residual = 1, Final residual = nan, No Iterations 1000

smoothSolver:  Solving for p, Initial residual = nan, Final residual = nan, No Iterations 1000

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

time step continuity errors : sum local = nan, global = nan, cumulative = nan

PIMPLE: iteration 2

smoothSolver:  Solving for Ux, Initial residual = nan, Final residual = nan, No Iterations 1000

smoothSolver:  Solving for Uy, Initial residual = nan, Final residual = nan, No Iterations 1000

smoothSolver:  Solving for e, Initial residual = nan, Final residual = nan, No Iterations 1000

--> FOAM FATAL IO ERROR:

wrong token type - expected Scalar, found on line 0 the word 'nan'

file: C:/Users/.../20180124/system/data.solverPerformance.e at line 0.

    From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)

    in file lnInclude/Scalar.C at line 93.

FOAM exiting

I would appreciate if you could help me know what is the problem and how I can solve it

[Mohd Faizal Mohamad]

Hi Foad,

Try to adjust your time step to smaller order.

--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+unsubscribe@googlegroups.com にメールを送信してください。
このグループに投稿するには open...@googlegroups.com にメールを送信してください。
https://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。

[Foad S Farimani]

Dear okamoto,

Thanks for the imediate reply.

I changed the deltaT from 0.002 to 0.0002

and I got a new error:

Time = 0.0016

Courant Number mean: 0.896348 max: 18.7598

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

PIMPLE: iteration 1

smoothSolver:  Solving for Ux, Initial residual = 0.065246, Final residual = 2.3006e-006, No Iterations 2

smoothSolver:  Solving for Uy, Initial residual = 0.0915461, Final residual = 2.71536e-006, No Iterations 3

smoothSolver:  Solving for e, Initial residual = 0.259906, Final residual = 5.07109e-006, No Iterations 14

smoothSolver:  Solving for p, Initial residual = 0.133062, Final residual = 4.84518e-007, No Iterations 9

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

time step continuity errors : sum local = 4.1513e-007, global = -2.97147e-007, cumulative = -1.95105e-006

PIMPLE: iteration 2

smoothSolver:  Solving for Ux, Initial residual = 0.0673966, Final residual = 1.95283e-006, No Iterations 2

smoothSolver:  Solving for Uy, Initial residual = 0.0426551, Final residual = 3.08232e-006, No Iterations 2

smoothSolver:  Solving for e, Initial residual = 0.109499, Final residual = 1.57047e+173, No Iterations 1000

--> FOAM FATAL ERROR:

Maximum number of iterations exceeded: 100

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]

    in file /opt/CFDSupportFOAM4.0/beta/OpenFOAM-dev/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 73.

FOAM aborting

Foad

このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。

[Foad S Farimani]

apparently my initial/boundary value for speed was in meters and too much for such a small geometry. I reduced it and it solves just fine.

F.

[sariew8]

Hi Foad

deltaT 0.0002  ---> 0.00001 or less

good luck!

2018年1月25日木曜日 23時01分20秒 UTC+9 Foad S Farimani: