[levogl]

Dear All,

I am facing a problem what I unable to solve, and I hope someone of you already met this problem and knows the solution.

Base case: OpenFOAM, multiphase (2 phase) flow (interFoam solver), LES turbulence model, checkMesh says OK.

Problem:
- The case works well under OpenFOAM for Windows 17.06 (created under this)
- Exact the same case (with the necessary modifications what I know) does not work under OpenFOAM for Windows 17.10

I solved, the following issues when I adapted the case to OF for W 17.10:
1. after executing "setFields", alpha.fuel file does not carry the BC-s. But it set well the alpha for the internalField. I solved by manually editing the file, but it is strange I think.

2. The missing keyword "pcorrFinal" issue solved by replacing "pcorr" by "pcorr.*"

After these, in OF 17.10 interFoam start to run, but after the first loop, the solver blows up: residual of "k" is nan meanwhile the initial residual is 1.

I tried to change Co and maxAlphaCo, change tolerances of "k", "p_rgh", set more correctors for PIMPLE but none of them worked.

Did someone face such a problem? Can someone give me some advice what would be the problem?

Log of interFoam:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type LES
Selecting LES turbulence model kEqn
Selecting LES delta type smooth
Selecting LES delta type cubeRootVol
kEqnCoeffs
{
    Ce              1.048;
    Ck              0.094;
}


Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 3.388116e-006, No Iterations 11
GAMGPCG:  Solving for pcorr, Initial residual = 0.04278774, Final residual = 2.661251e-006, No Iterations 4
GAMGPCG:  Solving for pcorr, Initial residual = 0.002614785, Final residual = 5.884809e-006, No Iterations 3
GAMGPCG:  Solving for pcorr, Initial residual = 0.0002990128, Final residual = 3.12335e-006, No Iterations 1
GAMGPCG:  Solving for pcorr, Initial residual = 4.117763e-005, Final residual = 9.902944e-007, No Iterations 1
GAMGPCG:  Solving for pcorr, Initial residual = 7.23902e-006, Final residual = 7.23902e-006, No Iterations 0
time step continuity errors : sum local = 6.399441e-008, global = -1.069007e-009, cumulative = -1.069007e-009
Courant Number mean: 0.004253773 max: 0.8613923

Starting time loop

Courant Number mean: 0.004253773 max: 0.8613923
Interface Courant Number mean: 0 max: 0
deltaT = 0.01
Time = 0.01

PIMPLE: iteration 1
MULES: Solving for alpha.fuel
Phase-1 volume fraction = 0.003936036  Min(alpha.fuel) = 0  Max(alpha.fuel) = 1.000025
MULES: Solving for alpha.fuel
Phase-1 volume fraction = 0.003950298  Min(alpha.fuel) = 0  Max(alpha.fuel) = 1.000049
MULES: Solving for alpha.fuel
Phase-1 volume fraction = 0.003964559  Min(alpha.fuel) = 0  Max(alpha.fuel) = 1.000071
MULES: Solving for alpha.fuel
Phase-1 volume fraction = 0.003978821  Min(alpha.fuel) = 0  Max(alpha.fuel) = 1.000091
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 6.973166e-008, No Iterations 29
GAMG:  Solving for p_rgh, Initial residual = 0.04714752, Final residual = 6.365322e-008, No Iterations 16
GAMG:  Solving for p_rgh, Initial residual = 0.003511602, Final residual = 7.682389e-008, No Iterations 12
GAMG:  Solving for p_rgh, Initial residual = 0.0003714307, Final residual = 9.656306e-008, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 4.977958e-005, Final residual = 7.345298e-008, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 7.690884e-006, Final residual = 7.099504e-008, No Iterations 2
time step continuity errors : sum local = 1.331502e-007, global = 3.895548e-008, cumulative = 3.788647e-008
GAMG:  Solving for p_rgh, Initial residual = 0.4310846, Final residual = 7.909259e-008, No Iterations 20
GAMG:  Solving for p_rgh, Initial residual = 0.02048172, Final residual = 6.497497e-008, No Iterations 13
GAMG:  Solving for p_rgh, Initial residual = 0.001358718, Final residual = 6.007829e-008, No Iterations 9
GAMG:  Solving for p_rgh, Initial residual = 0.000130214, Final residual = 7.332109e-008, No Iterations 5
GAMG:  Solving for p_rgh, Initial residual = 1.497639e-005, Final residual = 5.772462e-008, No Iterations 3
GAMGPCG:  Solving for p_rgh, Initial residual = 2.274067e-006, Final residual = 5.105044e-008, No Iterations 1
time step continuity errors : sum local = 8.880562e-008, global = 1.647652e-008, cumulative = 5.436299e-008
smoothSolver:  Solving for k, Initial residual = 1, Final residual = nan, No Iterations 1000
ExecutionTime = 6.525 s  ClockTime = 6 s