CFD Online Discussion Forums - reactingFoam stops at the beginning without error message

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - reactingFoam stops at the beginning without error message (https://www.cfd-online.com/Forums/openfoam-solving/215056-reactingfoam-stops-beginning-without-error-message.html)

reactingFoam stops at the beginning without error message

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:

Code:

/*---------------------------------------------------------------------------*\

  =========                 |

  \\      /  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.

Hi,

I have the same problem. My version is v2012.

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

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

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

Code:

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.

Code:

/*---------------------------------------------------------------------------*\

  =========                 |

  \\      /  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

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

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,

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).

Hi HPE,

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

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).

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.

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,