CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   reactingFoam SandiaD_LTS tutorial error (https://www.cfd-online.com/Forums/openfoam-solving/221412-reactingfoam-sandiad_lts-tutorial-error.html)

kit607 October 15, 2019 23:01

reactingFoam SandiaD_LTS tutorial error
 
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.

Code:

/*---------------------------------------------------------------------------*\
  =========                |
  \\      /  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 July 6, 2020 18:13

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

h_kuriakose December 6, 2020 08:34

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 December 6, 2020 18:27

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

kont87 February 9, 2021 08:09

Quote:

Originally Posted by Imed (Post 777019)
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:confused:
Thanks in advance!

kerim December 6, 2021 14:00

1 Attachment(s)
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 December 20, 2021 18:19

"cannot find patchField entry" problems for SandiaD_LTS
 
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::DimensionedField<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 January 11, 2022 10:11

Quote:

Originally Posted by ugata88 (Post 818914)
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::DimensionedField<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 January 13, 2022 07:24

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 January 13, 2022 12:15

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

altair98 January 26, 2022 20:07

Quote:

Originally Posted by kerim (Post 818006)
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 January 27, 2022 04:27

reactingFoam SandiaD_LTS tutorial error
 
Quote:

Originally Posted by altair98 (Post 821049)
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 January 27, 2022 08:40

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 https://www.cfd-online.com/Forums/op...ad-result.html
I wonder if a particular version produces accurate results.
Thank you again

kerim January 27, 2022 09:40

Quote:

Originally Posted by altair98 (Post 821095)
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 https://www.cfd-online.com/Forums/op...ad-result.html
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 January 27, 2022 18:19

Quote:

Originally Posted by kerim (Post 821101)
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 January 28, 2022 01:10

Quote:

Originally Posted by altair98 (Post 821128)
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.


All times are GMT -4. The time now is 18:56.