[Boku88]

Hallo guys,


i changed in the tutorial of reactingFoam/RAS/Sandia_LTS the blockMesh to study my simulation. Therefore I wanted to run the simulation without the wallthikness and than I put the two inlet pipes together:


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

convertToMeters 0.001;

vertices
(
(0 0 -150) // 0 E-1 jet
(3.6 -0.1571793944724 -150) // 1
(3.6 0.1571793944724 -150) // 2

(3.6 -0.1571793944724 -150) // 3 E-1 pilot ohne Wandung 0,25mm
(3.6 0.1571793944724 -150) // 4
(9.1 -0.3973145804719 -150) // 5
(9.1 0.3973145804719 -150) // 6

(0 0 0) // 7 E0 CH4
(3.6 -0.1571793944724 0) // 8
(3.6 0.1571793944724 0) // 9

(3.6 -0.1571793944724 0) // 10 E0 pilot ohne Wandung
(3.6 0.1571793944724 0) // 11
(9.1 -0.3973145804719 0) // 12
(9.1 0.3973145804719 0) // 13

(150 -6.54914143635 0) // 14 E0 coflow
(150 6.54914143635 0) // 15

(0 0 800) // 16 E1 jet
(3.6 -0.1571793944724 800) // 17
(3.6 0.1571793944724 800) // 18

(3.6 -0.1571793944724 800) // 19 E1 pilot
(3.6 0.1571793944724 800) // 20
(9.1 -0.3973145804719 800) // 21
(9.1 0.3973145804719 800) // 22

(150 -6.54914143635 800) // 23 E1 coflow
(150 6.54914143635 800) // 24
);

blocks
(
hex ( 0 1 2 0 7 8 9 7) (5 1 40) simpleGrading (1 1 1) //Jet
hex ( 3 5 6 4 10 12 13 11) (5 1 40) simpleGrading (1 1 1) //Pilot
hex ( 7 8 9 7 16 17 18 16) (5 1 400) simpleGrading (1 1 2) //Brennraum jet
hex (10 12 13 11 19 21 22 20) (5 1 400) simpleGrading (1 1 2) //Brennraum pilot
hex (12 14 15 13 21 23 24 22) (38 1 400) simpleGrading (3 1 2) //Brennraum coflow
);

boundary
(
inletCH4
{
type patch;
faces
(
(1 0 0 2)
);
}


wallOutside
{
type wall;
faces
(
(14 15 24 23)
);
}

wallTube
{
type wall;
faces
(
(1 2 9 8)
(4 3 10 11)
(5 6 13 12)
);
}

inletPilot
{
type patch;
faces
(
(5 3 4 6)
);
}


inletAir
{
type patch;
faces
(
(14 12 13 15)
);
}

outlet
{
type patch;
faces
(
(16 17 18 16)
(19 21 22 20)
(21 22 24 23)
);
}

axis
{
type empty;
faces
(
(0 7 7 0)
(7 16 16 7)
);
}


frontAndBack_pos
{
type wedge;
faces
(
(2 0 7 9)
(6 4 11 13)

(9 7 16 18)
(13 11 20 22)
(15 13 22 24)
);
}

frontAndBack_neg
{
type wedge;
faces
(
(0 1 8 7)
(3 5 12 10)

(7 8 17 16)

(10 12 21 19)
(12 14 23 21)
);
}


);

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












The blockMesh works and the checkMesh says ok, but if run the ractingFoam I have this error upcoming:











/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-1a0c91b3baa8
Exec : reactingFoam
Date : Nov 23 2018
Time : 18:13:03
Host : "boris-VirtualBox"
PID : 30108
I/O : uncollated
Case : /home/boris/CFD/OpenFOAM/dummy_23
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
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating reaction model

Selecting combustion model EDC<psiChemistryCombustion>
--> FOAM Warning :
From function static Foam::autoPtr<T> Foam::combustionModel::New(typename CombustionModel::reactionThermo&, const Foam::compressibleTurbulenceModel&, const Foam::word&) [with CombustionModel = Foam::CombustionModel<Foam:siReactionThermo>; typename CombustionModel::reactionThermo = Foam:siReactionThermo]
in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/combustionModels/lnInclude/combustionModelTemplates.C at line 69
Template parameters are no longer required when selecting a combustionModel. This information is now obtained directly from the thermodynamics. Actually 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 greyMeanAbsorptionEmission
Selecting scatterModel none
Selecting sootModel none

Starting time loop

--> FOAM Warning :
From function void Foam:article::locate(const vector&, const vector*, Foam::label, bool, Foam::string)
in file particle/particle.C at line 493
"Particle initialised with a location outside of the mesh."
--> FOAM Warning :
From function bool Foam::functionObjectList::read()
in file db/functionObjects/functionObjectList/functionObjectList.C at line 746
Caught FatalError
--> FOAM FATAL ERROR:
Unknown sample set type uniform

The sample set type uniform has been renamed lineUniform

Replace "type uniform;" with "type lineUniform;" for the set x_d_15 in "/home/boris/CFD/OpenFOAM/dummy_23/system/sampleDict.sets"

From function static Foam::autoPtr<Foam::sampledSet> Foam::sampledSet::New(const Foam::word&, const Foam:olyMesh&, const Foam::meshSearch&, const Foam::dictionary&)
in file sampledSet/sampledSet/sampledSet.C at line 161.

Time scales min/max:
Flow = 4.083092965e-05, 0.0001
Temperature = 4.494232837e+307, 4.494232837e+307
Composition = 4.494232837e+307, 4.494232837e+307
Overall = 4.083092965e-05, 0.0001
Time = 1

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 0.9999999754, Final residual = 0.002849317421, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.9999999993, Final residual = 0.0009886161378, No Iterations 1
DILUPBiCGStab: Solving for H2, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for H, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for O, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for O2, Initial residual = 0.99999994, Final residual = 0.003262125792, No Iterations 1
DILUPBiCGStab: Solving for OH, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for H2O, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for HO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for CH4, Initial residual = 1, Final residual = 0.00101207502, No Iterations 1
DILUPBiCGStab: Solving for CO, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DILUPBiCGStab: Solving for CO2, Initial residual = 1, Final residual = 0.00400863106, No Iterations 1
DICPCG: Solving for G, Initial residual = 1, Final residual = 9.182754125e-06, No Iterations 492
DILUPBiCGStab: Solving for h, Initial residual = 1, Final residual = 0.0002459968421, No Iterations 1
min/max(T) = 291, 1880
DICPCG: Solving for p, Initial residual = 0.9999999995, Final residual = 0.00924505048, No Iterations 35
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.987539369e-07, global = -1.213616321e-07, cumulative = -1.213616321e-07
DICPCG: Solving for p, Initial residual = 0.9659863197, Final residual = 9.407598315e-07, No Iterations 109
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.747274333e-08, global = 2.811622849e-10, cumulative = -1.210804698e-07
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.9999999885, Final residual = 0.02094788957, No Iterations 1
bounding epsilon, min: -7.545995368e+13 max: 1.692055105e+13 average: -1.129936214e+12
DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 0.0003996662989, No Iterations 1
bounding k, min: -4.462773291e+10 max: 76401942.73 average: -575399021.1
ExecutionTime = 4.89 s ClockTime = 5 s

Time scales min/max:
Flow = 3.769153522e-07, 0.0001
Temperature = 4.494232837e+307, 4.494232837e+307
Composition = 4.494232837e+307, 4.494232837e+307
Overall = 3.769153522e-07, 0.0001
Time = 2

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 0.6098529666, Final residual = 6.741850122e+11, No Iterations 1000
DILUPBiCGStab: Solving for Uz, Initial residual = 0.6797828315, Final residual = 0.04694567801, No Iterations 1




Can somebody help me out?


Thank you very much