CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

dieselFoam Error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 16, 2009, 11:50
Default dieselFoam Error
  #1
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 7
Rachel is on a distinguished road
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
Rachel is offline   Reply With Quote

Old   June 16, 2009, 14:45
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
According to your log you have a syntax error in your sprayProperties. Might be a problem.
mattijs is offline   Reply With Quote

Old   June 17, 2009, 05:37
Default
  #3
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 7
Rachel is on a distinguished road
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:rintStack(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>:perator=(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
Rachel is offline   Reply With Quote

Old   June 17, 2009, 06:59
Default
  #4
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 691
Rep Power: 18
niklas will become famous soon enough
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.
niklas is offline   Reply With Quote

Old   June 24, 2009, 10:22
Default
  #5
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 7
Rachel is on a distinguished road
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
Rachel is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
Compiling problems with hello worldC fw407 OpenFOAM Installation 21 January 6, 2008 17:38
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
user subroutine error CFDUSER CFX 2 December 9, 2006 06:31


All times are GMT -4. The time now is 01:10.