CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   dieselFoam Error (http://www.cfd-online.com/Forums/openfoam-solving/65488-dieselfoam-error.html)

Rachel June 16, 2009 11:50

dieselFoam Error
 
Hello everyone,

I am trying to solve a problem with following reaction.
Due to help from 2 other threads I have reached a point where the solver starts iterations, however it exits in the first time step.

Any clues why this could be heppening?

Thanks,
Rachel


Code:


Create time
Create mesh for time = 0
 
Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
}
Creating field DpDt
Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 4 and reactions = 1
Reading environmentalProperties
Reading combustion properties
Constructing Spray
--> FOAM Warning :
From function Cloud<ParticleType>::initCloud(const bool checkClass)
in file /home/rachel/OpenFOAM/OpenFOAM-1.5.x/src/lagrangian/basic/lnInclude/CloudIO.C at line 51
Cannot read particle positions file
"/scratch/rachel/coalgas_syn6a1_3spray/0/lagrangian/defaultCloud"
assuming the initial cloud contains 0 particles.
--> FOAM Warning :
From function entry::getKeyword(word& keyword, Istream& is)
in file db/dictionary/entry/entryIO.C at line 72
Reading /scratch/rachel/coalgas_syn6a1_3spray/constant/sprayProperties
found on line 183 the punctuation token '{'
expected either } or EOF
Selecting injectorType unitInjector
Selecting injectorType unitInjector
Selecting injectorType unitInjector
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel standardEvaporationModel
Selecting heatTransferModel RanzMarshall
Selecting wallModel remove
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
Selecting injectorModel hollowConeInjector
Selecting pdfType RosinRammler
Average Velocity for injector 0: 195.235 m/s, injection pressure = 247.667 bar
Average Velocity for injector 1: 195.235 m/s, injection pressure = 247.667 bar
Average Velocity for injector 2: 195.235 m/s, injection pressure = 247.667 bar
Constructing three dimensional spray injection.
Courant Number mean: 0 max: 0.000980191
Starting time loop
Courant Number mean: 0 max: 0.0980191
deltaT = 0.00025
Time = 0.00025
Evolving Spray
Solving chemistry
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 7.0502e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.737e-10, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.01832e-09, No Iterations 2
DILUPBiCG: Solving for CO, Initial residual = 1, Final residual = 6.30454e-07, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
DILUPBiCG: Solving for H2O, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
DILUPBiCG: Solving for H2, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator=(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
#4 main in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
#5 __libc_start_main in "/lib64/libc.so.6"
#6 __gxx_personality_v0 in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
Segmentation fault


mattijs June 16, 2009 14:45

According to your log you have a syntax error in your sprayProperties. Might be a problem.

Rachel June 17, 2009 05:37

Thanks Mattijs,

I was trying to use 3 sprays in the system. Hence I modified sprayProperties.

Now i have switched back to single spray, to test the chemistry. However there is something else thats going wrong and segmentation fault is not that helpful.

Here is the new log


Quote:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5.x |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : /home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam
Date : Jun 17 2009
Time : 10:21:29
Host : euler
PID : 23405
Case : /scratch/rachel/gas6a1
nProcs : 1
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0

Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
}
Creating field DpDt
Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 4 and reactions = 1
Reading environmentalProperties
Reading combustion properties
Constructing Spray
--> FOAM Warning :
From function Cloud<ParticleType>::initCloud(const bool checkClass)
in file /home/rachel/OpenFOAM/OpenFOAM-1.5.x/src/lagrangian/basic/lnInclude/CloudIO.C at line 51
Cannot read particle positions file
"/scratch/rachel/gas6a1/0/lagrangian/defaultCloud"
assuming the initial cloud contains 0 particles.
Selecting injectorType unitInjector
Selecting injectorType unitInjector
Selecting injectorType unitInjector
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel standardEvaporationModel
Selecting heatTransferModel RanzMarshall
Selecting wallModel remove
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
Selecting injectorModel hollowConeInjector
Selecting pdfType RosinRammler
Average Velocity for injector 0: 195.235 m/s, injection pressure = 247.667 bar
Average Velocity for injector 1: 195.235 m/s, injection pressure = 247.667 bar
Average Velocity for injector 2: 195.235 m/s, injection pressure = 247.667 bar
Constructing three dimensional spray injection.
Courant Number mean: 0 max: 0.000980191
Starting time loop
Courant Number mean: 0 max: 0.0980191
deltaT = 0.00025
Time = 0.00025
Evolving Spray
Solving chemistry
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 7.0502e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.737e-10, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.01832e-09, No Iterations 2
DILUPBiCG: Solving for CO, Initial residual = 1, Final residual = 6.30454e-07, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
DILUPBiCG: Solving for H2O, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
DILUPBiCG: Solving for H2, Initial residual = 1, Final residual = 6.30453e-07, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator=(Foam::tmp<Foam::Geometri cField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
#4 main in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
#5 __libc_start_main in "/lib64/libc.so.6"
#6 __gxx_personality_v0 in "/home/rachel/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPOpt/dieselFoam"
Segmentation fault

niklas June 17, 2009 06:59

my guess is that your timestep is too large.

200 m/s injection velocity and 0.25 ms integration step.
that means your parcels will travel approx 5 cm.
The momenum transfer to the gas will be too high.

Rachel June 24, 2009 10:22

Thanks everybody,

It was a problem of setting the right initial/boundary conditions. Since the domain had inlet & outlet, I used BC similar to reactingFoam (from wiki) and then used dieselFoam solver. Now awaiting the results...

Thanks,
Rachel


All times are GMT -4. The time now is 22:41.