[infrared]

Hello everyone,

I am a newcomer to using OpenFOAM. I tried to run an impinging jet simulation and modified a previous case. I have the 'constant' folder, 'system' folder, '.foam' file, and a results folder named '533' containing the iteration data. As I understand, we should use the '0' folder instead of copying the previous case data. I tried to modify the necessary files in the '533' folder, but it didn't work well. It works when I run 'blockMesh', but I cannot run the solver (simpleFoam/pisoFoam/icoFoam).

Is there anything wrong? Please give me suggestions. I appreciate your kindness.

Thank you!

[Yann]

What are you modifying between the old and new case?

[infrared]

I want to modify the mesh to a coarser one. The old case is a non-reacting jet, but in the new one, I would like to perform the reacting/combustion process as well.

[Yann]

Alright, you cannot just copy the result folder because its content won't match your new mesh.

you should instead use the mapFields utility to map the results of your first case on your new mesh.

Cheers,
Yann

[infrared]

Originally Posted by Yann

Alright, you cannot just copy the result folder because its content won't match your new mesh.

you should instead use the mapFields utility to map the results of your first case on your new mesh.

Cheers,
Yann

Hello Yann, thank you for your response.

I did it by using the '0' directory instead of the '533' result. Now, I am going to change the previous case from a non-reacting condition to a reacting/combustion condition. I tried to modify it as per the attachments, but it didn't work well. Could you help me fix it?

You can view the attachments here
https://drive.google.com/drive/folde...gU?usp=sharing

Thank you for your kindness.

[Yann]

Hello,

Could you elaborate? What did you do and what is not working?

[infrared]

Originally Posted by Yann

Hello,

Could you elaborate? What did you do and what is not working?

It appears the message when I wrote the command in the terminal as "reactingFoam", the message is

--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Entry 'CHEMKINFile' not found in dictionary "/home/ishrein/OpenFOAM/ishrein-v2312/run/combCoarse/constant/thermophysicalProperties"


file: constant/thermophysicalProperties at line 20 to 118.

From bool Foam::dictionary::readEntry(const Foam::word&, T&, Foam::keyType:ption, Foam::IOobjectOption::readOption) const [with T = Foam::fileName]
in file ./src/OpenFOAM/lnInclude/dictionaryTemplates.C at line 327.

Could you help me to fix it? Thank you for your kindness.

[Yann]

I'm not very familiar with reaction solvers, but I think you get this error because your CHEMKINFile parameter is misplaced in thermophysicalProperties.

Try something like this (basically getting CHEMKINFile and CHEMKINThermoFile out of the mixture dictionary) and see if it solves the error:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.97; // Adjust based on your mixture
    }
    thermodynamics
    {
        Tlow            200;
        Thigh           5000;
        Tcommon         1000;
        Hf              0;
    }
    transport
    {
        As              1.67212e-6;
        Ts              170.672;
    }
    
}

CHEMKINFile         "chem.inp";
CHEMKINThermoFile   "therm.dat";

[...]

[infrared]

Originally Posted by Yann

I'm not very familiar with reaction solvers, but I think you get this error because your CHEMKINFile parameter is misplaced in thermophysicalProperties.

Try something like this (basically getting CHEMKINFile and CHEMKINThermoFile out of the mixture dictionary) and see if it solves the error:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.97; // Adjust based on your mixture
    }
    thermodynamics
    {
        Tlow            200;
        Thigh           5000;
        Tcommon         1000;
        Hf              0;
    }
    transport
    {
        As              1.67212e-6;
        Ts              170.672;
    }
    
}

CHEMKINFile         "chem.inp";
CHEMKINThermoFile   "therm.dat";

[...]

Thank you for your response. The CHEMKINFile and CHEMKINThermoFile parameters are already outside of the mixture dictionary. I don't know why I always encounter error messages like the one listed above.

[Yann]

As long as you get this error:

--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Entry 'CHEMKINFile' not found in dictionary "/home/ishrein/OpenFOAM/ishrein-v2312/run/combCoarse/constant/thermophysicalProperties"

It means something is wrong and the solver cannot find the entry where it is supposed to be. So it's either misplaced, or there is a typo something in the file.

What is your current thermophysicalProperties file?

[infrared]

Originally Posted by Yann

As long as you get this error:

--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Entry 'CHEMKINFile' not found in dictionary "/home/ishrein/OpenFOAM/ishrein-v2312/run/combCoarse/constant/thermophysicalProperties"

It means something is wrong and the solver cannot find the entry where it is supposed to be. So it's either misplaced, or there is a typo something in the file.

What is your current thermophysicalProperties file?

Here is the thermophysicalProperties file:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2112 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

mixture
{
specie
{
nMoles 1;
molWeight 28.97;
}
thermodynamics
{
Tlow 200;
Thigh 5000;
Tcommon 1000;
Hf 0;
}
transport
{
As 1.67212e-6;
Ts 170.672;
}

CHEMKINFile "chem.inp";
CHEMKINThermoFile "therm.dat";
}

species
(
O2
{
molWeight 32.0;
thermodynamics
{
Tlow 200;
Thigh 5000;
Tcommon 1000;
Hf 0;
}
transport
{
As 1.67212e-6;
Ts 170.672;
}
}
H2
{
molWeight 2.016;
thermodynamics
{
Tlow 200;
Thigh 5000;
Tcommon 1000;
Hf 0;
}
transport
{
As 1.67212e-6;
Ts 170.672;
}
}
H2O
{
molWeight 18.01528;
thermodynamics
{
Tlow 200;
Thigh 5000;
Tcommon 1000;
Hf -241.826; // Standard enthalpy of formation in kJ/mol
}
transport
{
As 1.67212e-6;
Ts 170.672;
}
}
N2
{
molWeight 28.0134;
thermodynamics
{
Tlow 200;
Thigh 5000;
Tcommon 1000;
Hf 0;
}
transport
{
As 1.67212e-6;
Ts 170.672;
}
}
);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

What do you think about that?

[Yann]

This is what I was saying in my previous post. Your CHEMKINFile is inside the mixture dictionary:

mixture
{
	specie
	{
		nMoles 1;
		molWeight 28.97;
	}
	thermodynamics
	{
		Tlow 200;
		Thigh 5000;
		Tcommon 1000;
		Hf 0;
	}
	transport
	{
		As 1.67212e-6;
		Ts 170.672;
	}
	
	CHEMKINFile "chem.inp";
	CHEMKINThermoFile "therm.dat";
}

Try to get it out:

mixture
{
	specie
	{
		nMoles 1;
		molWeight 28.97;
	}
	thermodynamics
	{
		Tlow 200;
		Thigh 5000;
		Tcommon 1000;
		Hf 0;
	}
	transport
	{
		As 1.67212e-6;
		Ts 170.672;
	}
}	

CHEMKINFile "chem.inp";
CHEMKINThermoFile "therm.dat";

[infrared]

Originally Posted by Yann

This is what I was saying in my previous post. Your CHEMKINFile is inside the mixture dictionary:

mixture
{
	specie
	{
		nMoles 1;
		molWeight 28.97;
	}
	thermodynamics
	{
		Tlow 200;
		Thigh 5000;
		Tcommon 1000;
		Hf 0;
	}
	transport
	{
		As 1.67212e-6;
		Ts 170.672;
	}
	
	CHEMKINFile "chem.inp";
	CHEMKINThermoFile "therm.dat";
}

Try to get it out:

mixture
{
	specie
	{
		nMoles 1;
		molWeight 28.97;
	}
	thermodynamics
	{
		Tlow 200;
		Thigh 5000;
		Tcommon 1000;
		Hf 0;
	}
	transport
	{
		As 1.67212e-6;
		Ts 170.672;
	}
}	

CHEMKINFile "chem.inp";
CHEMKINThermoFile "therm.dat";

Sorry, the last one is modified. I tried to write the CHEMKINFile both outside and inside the mixture, but I got the same error message.

[Yann]

Alright, I just tried myself.
It is important to carefully read the error messages.

You initial error message states:

--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Entry 'CHEMKINFile' not found in dictionary "/home/ishrein/OpenFOAM/ishrein-v2312/run/combCoarse/constant/thermophysicalProperties"

If you move the CHEMKINFile and CHEMKINThermoFile out of the mixture dictionary, you should get this error:

--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Entry 'CHEMKINTransportFile' not found in dictionary "/home/ishrein/OpenFOAM/ishrein-v2312/run/combCoarse/constant/thermophysicalProperties"

Which is not the same thing. Now the solver cannot find the CHEMKINTransportFile parameter, as it is indeed not defined.

[infrared]

Thank you for your response,
After trying to inspect the Chemkin file problem, I found out that one of the causes is the sum of each chemical compound not being 1; it must be 1, isn't it? However, I also discovered another problem:

Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

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 laminar
Selecting laminar stress model Stokes
Creating reaction model

Selecting combustion model laminar
Selecting chemistry solver
{
solver EulerImplicit;
method standard;
}

StandardChemistryModel: Number of species = 5 and reactions = 1
using integrated reaction rate
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop


BilgerMixtureFraction BilgerMixtureFraction1:
Courant Number mean: 0 max: 0
deltaT = 0.120482
Time = 0.120482

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
[stack trace]
=============
#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#5 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#6 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#7 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#8 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so
#9 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingFoam
#10 ? in /lib/x86_64-linux-gnu/libc.so.6
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingFoam
=============
Floating point exception

Could you give me some suggestions as to why I experienced these problems?
Thank you in advance.

Best regards,
Ir

[Yann]

Hello,

This kind of error usually happens when there is some bad math happening. Typically trying to divide by 0, which might be what's happening here.

Now the question is why is it happening?

The crash is happening when the solver is supposed to solve the velocity and the error message mentions turbulence libraries, so it's most probably related to one of these things.

Since it is crashing at the first iteration, first thing to do is to check your boundary conditions to make sure there is nothing wrong with the initialization and conditions defined.

[infrared]

Hello Yann,

I tried to check the initialization and its boundary condition for days but it didn't work out well. Could you give me some suggestions about that?

Here is the initialization and its boundary condition.

1. Velocity (U)
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (1 0 0);
}

outlet
{
type zeroGradient;
}

top_wall
{
type noSlip;
}

bottom_wall
{
type noSlip;
}

}

2. Pressure
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
bottom_wall
{
type zeroGradient;
}

inlet
{
type zeroGradient;
}

top_wall
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 0;
}
}

3. Temperature
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
inlet
{
type fixedValue;
value uniform 300;
}

outlet
{
type zeroGradient;
}

top_wall
{
type zeroGradient;
}

bottom_wall
{
type zeroGradient;
}
}

4. k
*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.005;

boundaryField
{
bottom_wall
{
type kqRWallFunction;
value uniform 0.1;
}

top_wall
{
type kqRWallFunction;
value uniform 0.1;
}

inlet
{
type fixedValue;
value uniform 0.1;
}

outlet
{
type zeroGradient;
}
}

5. epsilon
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 0.01; // Adjust initial turbulence dissipation rate as required

boundaryField
{
bottom_wall
{
type epsilonWallFunction;
value uniform 0.1;
}

top_wall
{
type epsilonWallFunction;
value uniform 0.1;
}

inlet
{
type fixedValue;
value uniform 0.1;
}

outlet
{
type zeroGradient;
}
}

6. controlDict
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application reactingFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 100;

deltaT 1;

writeControl timeStep;

writeInterval 1;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

maxCo 0.5;

adjustTimeStep yes;

debugSwitches
{
default 1;
general 1;
}

7. fvSchemes
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

8. fvSolution
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}

U
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-05;
relTol 0.1;
}

"(k|epsilon|omega)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-05;
relTol 0.1;
}
}

PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
}

Thank you in advance.

[Yann]

Hello,

Try setting an absolute pressure instead of 0.

Yann

[infrared]

Hello Yann,

Finally, I decided to adopt the tutorial case of 'counterFlowFlame2D' with the reactingFoam solver. The domain is cylindrical with boundary conditions as follows: top wall (noSlip), bottom wall (noSlip), side wall (inletOutlet), and an inlet located at the center of the top wall. The flow is injected through the inlet, and no flow should exit the domain.

I have already performed the simulation, but it appears that the flow was exiting through the bottom wall.

Here are the initial conditions for U, T, and p:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2312 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type turbulentInlet;
fluctuationScale (0.05 0.05 0.05);
referenceField uniform (0 0 1);
alpha 1;
value uniform (0 0 1);
}

outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}

bottom_wall
{
type noSlip;
}

top_wall
{
type fixedValue;
value uniform (0 0 0);
}

spray
{
type fixedValue;
value uniform (0 0 1);
}
}



// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2312 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{
inlet
{
type fixedValue;
value uniform 300;
}

outlet
{
type zeroGradient;
}

bottom_wall
{
type fixedValue;
value uniform 373;
}

top_wall
{
type fixedValue;
value uniform 373;
}
}

// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2312 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
inlet
{
type zeroGradient;
}

outlet
{
type totalPressure;
p0 $internalField;
}

bottom_wall
{
type zeroGradient;
}

top_wall
{
type zeroGradient;
}
}

// ************************************************** *********************** //

Do you have any suggestions for ensuring that there is no flow exiting from the bottom wall, so the flow impinges on the wall?

I really appreciate your help. Thank you.

[Yann]

Hello,

Not sure what's wrong, your BCs look fine to me.

Yann