[n936an]

Hello,

Running a compressible flow over a 737 wing here. As soon as it reaches the first timestep (t=0.1), the final residual for "solving for e" comes out as "nan" and over the next few timesteps all other values also inevitably end up as nan.

I don't have much experience in setting up new cases so I'd appreciate it if you could look at my files and see if I made any mistakes.

Below is my case directory.
https://drive.google.com/file/d/1yR3...ew?usp=sharing

And here is what I see when running the simulation:

Create time

Create mesh for time = 0

Selecting solver fluid
Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
RAS
{
    model           kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

No MRF models present

Selecting patchDistMethod meshWave
Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity

PIMPLE: Convergence criteria found
        p: tolerance 0.0001
        U: tolerance 0.0001
        "(k|omega|epsilon)": tolerance 0.0001


PIMPLE: Operating solver in steady-state mode with 1 outer corrector
PIMPLE: Operating solver in SIMPLE mode



Starting time loop

forcesBase forces:
    Not including porosity effects
forces forces:
    Not including porosity effects
forceCoeffs forces:
    Not including porosity effects
forceCoeffs forces write:
    Cm    = 0.000956058
    Cd    = 0.000574535
    Cl    = 4.41614e-06
    Cl(f) = 0.000958266
    Cl(r) = -0.00095385

Time = 0.1s

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0667967, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0593224, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0677076, No Iterations 2
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
time step continuity errors : sum local = nan, global = nan, cumulative = nan
smoothSolver:  Solving for omega, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = nan, Final residual = nan, No Iterations 1000
ExecutionTime = 66.7093 s  ClockTime = 67 s

forceCoeffs forces write:
    Cm    = -nan
    Cd    = -nan
    Cl    = -nan
    Cl(f) = -nan
    Cl(r) = -nan

I believe a similar topic was posted several years ago but that thread did not receive a full answer it seems.

Many thanks.