[kit607]

Hi all. I am using OpenFOAM 6 and I am trying the tutorial SandiaD_LTS. I just use the Allrun file to start the simulation without modifying any parameters. However, there is error for the setFields which shows below.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-fa1285188035
Exec   : setFields
Date   : Oct 16 2019
Time   : 10:47:04
Host   : "mgt-aaibe"
PID    : 30000
I/O    : uncollated
Case   : /home/mgt-aaibe/OpenFOAM/mgt-aaibe-6/run/combustion_test
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 = 1500

Reading setFieldsDict

Setting field default values
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field T not found
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field N2 not found
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field O2 not found
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field CH4 not found

Setting field region values
    Adding cells with center within boxes 1((0 -10 -100) (0.0036 10 0))
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field CH4 not found
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field O2 not found
--> FOAM Warning : 
    From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 117
    Field N2 not found

End

Is my installtion of OF got error or I miss out some step when I run the simulation? Thank you

[Imed]

hi there;
i had the same problem but I could fix it, if you still want fixing the problem I can help

[h_kuriakose]

I know, it's a late reply. But incase you need it, you can find the documentation and simulation files for Sandia D flame tutorial


https://www.cemf.ir/reactingfoam-how...e-by-openfoam/

[Imed]

Hi;
Thank you very much for you replay, it is very useful thanks again.

[kont87]

Originally Posted by Imed

hi there;
i had the same problem but I could fix it, if you still want fixing the problem I can help

Hello! Could you help me? Having the same problem
Thanks in advance!

[kerim]

Dear all. I am using OpenFOAM 7 and I am trying the tutorial SandiaD_LTS without any modification. I just used the Allrun file to start. But I have Floating point exception (core dumped) error after 3719 iteration.

Could you help me?

[ugata88]

Hello everyone, I am new to openfoam and need to simulate nonpremixed combustion in openfoam-9, for this ı select SandiaD_LTS solver

1 .I created axi-symmetric cylinder geometry for 2D case on Salome-Meca 2015, also used seperated small angle (3 degre) wedge planes,

2 .I created patchfiled for inletCH4 (patch), inletair(patch), outlet (patch), frontandBack_pos (wedge) ,frontandBack_neg (wedge)

3. for meshing I used Netgen 2d3d algorithm and 3d parameters hypotheses

4. Exporting .unv mesh file to openfoam

5. I used chemkinToFoam, ideasUnvToFoam, so that constant/polymesh file is appeared, checkMesh is OK ,

6. After that, I arranged system/setFieldsDict for initializing inletCH4 field and running reactinFoam, so I encountered the following message;


--> FOAM FATAL IO ERROR:
Cannot find patchField entry for frontandBack_pos

file: /home/ugur/OpenFOAM/ugur-9/run/SandiaD_LTS/0/T/boundaryField from line 25 to line 52.

From function void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam:imensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /home/ubuntu/OpenFOAM/OpenFOAM-9/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 184.

FOAM exiting

about the message, also I selected wedge patch type for frontandBack_pos and frontandBack_neg in constant / polyMesh /boundary

6
(
inletCH4
{
type patch;
nFaces 7;
startFace 64411;
}
inletair
{
type patch;
nFaces 20;
startFace 64418;
}
outlet
{
type patch;
nFaces 27;
startFace 64438;
}
wallTube
{
type patch;
nFaces 593;
startFace 64465;
}
frontandBack_pos
{
type patch;
nFaces 12980;
startFace 65058;
}
frontandBack_neg
{
type patch;
nFaces 13215;
startFace 78038;
}
)


and, also I selected initial temperature for my boundaryField in 0/T , as like this ;

dimensions [0 0 0 1 0 0 0];

internalField uniform 295;

boundaryField
{
inletCH4
{
type fixedValue;
value uniform 295;
}

wallTube
{
type zeroGradient;
}

inletair
{
type fixedValue;
value uniform 295;
}

outlet
{
type zeroGradient;
}

frontAndBack_pos
{
type wedge;
}

frontAndBack_neg
{
type wedge;
}
}


I feel so helpless and tired. Does anyone have a solution for this issue or have encountered this type of problem?
Thank you for your help in advance.

[JulioPieri]

Originally Posted by ugata88

Hello everyone, I am new to openfoam and need to simulate nonpremixed combustion in openfoam-9, for this ı select SandiaD_LTS solver

1 .I created axi-symmetric cylinder geometry for 2D case on Salome-Meca 2015, also used seperated small angle (3 degre) wedge planes,

2 .I created patchfiled for inletCH4 (patch), inletair(patch), outlet (patch), frontandBack_pos (wedge) ,frontandBack_neg (wedge)

3. for meshing I used Netgen 2d3d algorithm and 3d parameters hypotheses

4. Exporting .unv mesh file to openfoam

5. I used chemkinToFoam, ideasUnvToFoam, so that constant/polymesh file is appeared, checkMesh is OK ,

6. After that, I arranged system/setFieldsDict for initializing inletCH4 field and running reactinFoam, so I encountered the following message;


--> FOAM FATAL IO ERROR:
Cannot find patchField entry for frontandBack_pos

file: /home/ugur/OpenFOAM/ugur-9/run/SandiaD_LTS/0/T/boundaryField from line 25 to line 52.

From function void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam:imensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /home/ubuntu/OpenFOAM/OpenFOAM-9/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 184.

FOAM exiting

about the message, also I selected wedge patch type for frontandBack_pos and frontandBack_neg in constant / polyMesh /boundary

6
(
inletCH4
{
type patch;
nFaces 7;
startFace 64411;
}
inletair
{
type patch;
nFaces 20;
startFace 64418;
}
outlet
{
type patch;
nFaces 27;
startFace 64438;
}
wallTube
{
type patch;
nFaces 593;
startFace 64465;
}
frontandBack_pos
{
type patch;
nFaces 12980;
startFace 65058;
}
frontandBack_neg
{
type patch;
nFaces 13215;
startFace 78038;
}
)


and, also I selected initial temperature for my boundaryField in 0/T , as like this ;

dimensions [0 0 0 1 0 0 0];

internalField uniform 295;

boundaryField
{
inletCH4
{
type fixedValue;
value uniform 295;
}

wallTube
{
type zeroGradient;
}

inletair
{
type fixedValue;
value uniform 295;
}

outlet
{
type zeroGradient;
}

frontAndBack_pos
{
type wedge;
}

frontAndBack_neg
{
type wedge;
}
}


I feel so helpless and tired. Does anyone have a solution for this issue or have encountered this type of problem?
Thank you for your help in advance.

Seems that you have a lower case problem there. Your mesh is defined with frontandBack_pos, while your BCs are frontAndBack_pos

[ugata88]

Julio, thanks for attention first of all, You're right, I realized latter. For meshing I used blockMesh utility instead of Salome, Netgen algorithm. Becasuse of this, I annihilated skewness and ortogonality problems. But, ı have a new problem about initial k-epsilon values

[JulioPieri]

For that you check the OF user guide or other regular resources. There are typical calculations for that!

[altair98]

Originally Posted by kerim

Dear all. I am using OpenFOAM 7 and I am trying the tutorial SandiaD_LTS without any modification. I just used the Allrun file to start. But I have Floating point exception (core dumped) error after 3719 iteration.

Could you help me?

Were you able to solve this issue?
I am facing the same.

[kerim]

Originally Posted by altair98

Were you able to solve this issue?
I am facing the same.

Unfortunately, I haven't solved this problem. The same situation with newer versions of OpenFOAM8 and OpenFOAM9. However, older versions work well!

[altair98]

Thank you for letting me know. Could you mention which version seems to work because OpenFOAM v6 did not produce accurate results as mentioned in Sandia Flame D OpenFoam tutorial case - bad result?
I wonder if a particular version produces accurate results.
Thank you again

[kerim]

Originally Posted by altair98

Thank you for letting me know. Could you mention which version seems to work because OpenFOAM v6 did not produce accurate results as mentioned in Sandia Flame D OpenFoam tutorial case - bad result?
I wonder if a particular version produces accurate results.
Thank you again

There is no ignition in the OF6 standard tutorial. Please try another turbulence model, for example, RNG k-epsilon model.

[altair98]

Originally Posted by kerim

There is no ignition in the OF6 standard tutorial. Please try another turbulence model, for example, RNG k-epsilon model.

Running in OF7 with any other turbulence model produces good results as compared with experimental data from https://tnfworkshop.org/data-archive...edjet/ch4-air/
Thank you for your suggestions!

[kerim]

Originally Posted by altair98

Running in OF7 with any other turbulence model produces good results as compared with experimental data from https://tnfworkshop.org/data-archive...edjet/ch4-air/
Thank you for your suggestions!

I'm glad to hear that. Good luck.