[sosos]

Hello everybody,

I am trying to run a combustion simulation using reactingFoam, however when I start the calculation it loads the solvers and other parameters and that's it. It stops whitout any error message. I am not really what you can call an expert of OpenFoam, so maybe I am missing an obvious thing but I am quite lost. So, thank you for your time and help.

This is what I get:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-d3fd147e6c65
Exec   : reactingFoam
Date   : Feb 21 2019
Time   : 11:13:28
Host   : "DESKTOP-U8LEB3K"
PID    : 619
I/O    : uncollated
Case   : /mnt/d/USN/4th_semester/test/turbulent_combustion
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in PISO mode

Using LTS
Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Prt             0.85;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating reaction model

Selecting combustion model EDC
Selecting chemistry solver
{
    solver          ode;
    method          TDAC;
}

StandardChemistryModel: Number of species = 5 and reactions = 1
Selecting chemistry reduction method DAC
Selecting chemistry tabulation method ISAT
Selecting ODE solver seulex
    using integrated reaction rate
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present

Starting time loop

Time scales min/max:
    Flow        = 0.0001, 0.0001

If you need any other information, please tell me cause I don't really know what can be helpful and what can't.

[evrenykn]

Hi,

I have the same problem. My version is v2012.

[HPE]

Would it be possible for you to reproduce the issue with a minimal working example, so that we can reproduce the issue locally?

[francescomarra]

Dear Sosos,
how did you launch your run? By using scripts (like Allrun) or redirecting the output on file, in my experience it can happen that the last lines of the output can be missed, preventing the possibility to see in which part of the code the error is occurring.
Try to launch reactingFoam without arguments directly (interactively) on the command line to see if further messages arise.
And of course, as indicated by HPE, it is impossible to have an idea of what is happening knowing nothing of what you changed starting from a working tutorial.

Regards,

Francesco

[evrenykn]

Hi Francesco,
Hi HPE,

I think sosos may have stopped following this topic. He created this start of 2019. But I can't pass this step.

I solved the SandiaD_LTS tutorial normally. After this solution, I imported without using Allrun a different model using ideasUnvToFoam and I changed boundary conditions to suit my model in 0 folder.

After the changes, i runned chemkinToFoam, setFields and reactingFoam.

The solver not runned because of this warning

Invalid wall function specification
    Patch type for patch wall must be wall
    Current patch type is patch

After this situation, I changed the wall type from patch to wall in Polymesh>boundary.

Finally, I re-runned the setFields and reactingFoam and I attached the response as below.

Meanwhile, thinking there might be a difference, I changed the version from 2012 to 8.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 8-340defec456f
Exec   : reactingFoam
Date   : Feb 25 2021
Time   : 23:49:04
Host   : "host"
PID    : 4908
I/O    : uncollated
Case   : /mnt/d/OF/x1/x1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Using LTS
Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         multiComponentMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    model           kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model eddyDiffusivity
Creating reaction model

Selecting combustion model EDC
Selecting chemistry solver
{
    solver          ode;
    method          TDAC;
}

StandardChemistryModel: Number of species = 36 and reactions = 219
Selecting chemistry reduction method DAC
element not considered
Selecting chemistry tabulation method ISAT
Selecting ODE solver seulex
    using integrated reaction rate
Creating field dpdt

Creating field kinetic energy K

No MRF models present

Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type radiation
    Source: radiation
Selecting radiationModel P1
Selecting absorptionEmissionModel greyMeanCombustion
Selecting scatterModel none
Selecting soot model none

Starting time loop

Time scales min/max:
    Flow        = 0.0001, 0.0001

[francescomarra]

Dear Evren,
it is very difficult to debug this error without knowing the entire setup of the simulation. From your output, it seems that the error occurs while computing the time scales for the temperature (the next line in the log should be something similar to
Temperature = 4.494232837e+307, 4.494232837e+307
).
I would try to have a uniform initial temperature field without setting a hot zone (the setField command) for the ignition to check if the problem is in the setting of the ignition hot zone. Maybe, by changing the domain and mesh type, some inconsistency occurred.
Best regards,
Francesco

[evrenykn]

Hi again,

I checked Temperature fields, tried to close radiation and solve without allrun script, but I can't changed anything.

In addition I searched "Time scales min/max:" and i found in setRDeltaT.H file. I added a debug code below alphaTemp and alphaY statements.

When I change the codes in that file, shouldn't the text I get on the screen change? I want to see a difference and changed the code from Info<< " Flow = " to Info<< " Flo = " , but it didn't change.

I added the simulation files as below.

Model

If you have a suitable time, Would you be able to examine it, please?

Regards,

[HPE]

Can you run a test by switching off the chemistry/combustion?

If that runs, could you run the case without TDAC (hence just using ODE)?

Then if everything is fine, still let the flow develop a bit without jumping into the chemistry (which is a bit painful to my experience).

[evrenykn]

Hi HPE,

I tried to turn off chemistry. Unfortunately, it didn't solve.

[HPE]

So, one point (chemistry/combustion) has already been isolated out.

Keep switching off the submodules further till the case is OK for the first time, or debug the solver with the simulation you have (this involves coding/compilation etc).

[evrenykn]

Hi again,

I solved the problem but I don't know how I solved it.

I re-created the project in a new folder and I checked from start line by line. The solver started.

[Bk12345]

Hello all,

I know, this is old thread, but my issue is similar. My solver crashed after just starting the 1st iteration and like above mentioned(francescomarra: also temperature something 1e+300). by the way, I have tried most of suggestions like turning off chemsitry, radiation, and also disabling setfields. those didn't work out. HereI am adding screenshot.

Best Regards,