
[Sponsors] 
September 30, 2014, 05:18 
Segmentation fault when using reactingFOAM for Fluids

#1 
New Member
Tommy Floessner
Join Date: Sep 2014
Posts: 1
Rep Power: 0 
Dear CFDcommunity,
i tried to simulate a reactive flow of fluids in OF230. As a base i used the tutorialcase for reactingFoam and applied several changes that would describe my problem. The most important change is the mainSolver to rhoReactingFoam. I'll post the logfile for further Informations about the used thermophysical, transport, combustion, etc. models. The Problem in the set up case is a segmentation fault during the initialization of the field. (see logfile) I know that the first thing to check, in case of segmentation faults, is the geometry and boundary conditions. This i tested sucessfully with a methancombustion in the reactingFoam tutorialcase. But i don't want to exclude it as the root of the error by default. I hope somebody can give me one or two hints, that could help solving the problem. If any further Information is needed to get rid of the question tags, just ask. Thanks Tommy Logfile: Code:
/**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.3.0   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ Build : 2.3.0f5222ca19ce6 Exec : rhoReactingFoam Date : Sep 15 2014 Time : 14:06:53 Host : "dreadMS7235" PID : 9818 Case : /home/dread/OpenFOAM/run/rrf_002_poly nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Disallowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading g Creating reaction model Selecting combustion model PaSR<rhoChemistryCombustion> Selecting chemistry type { chemistrySolver ode; chemistryThermo rho; } Selecting thermodynamics package { type heRhoThermo; mixture reactingMixture; transport polynomial; thermo hPolynomial; energy sensibleEnthalpy; equationOfState icoPolynomial; specie specie; } AMI: Creating addressing and weights between 500 source faces and 500 target faces AMI: Patch source sum(weights) min/max/average = 0.999915, 1, 0.999992 AMI: Patch target sum(weights) min/max/average = 1, 1, 1 Selecting chemistryReader foamChemistryReader // end Logfile // additional information #0 Foam::error::printStack(Foam::Ostream&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 at rhoReactionThermos.C:0 #4 Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::correctMassFractions() in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #5 Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::multiComponentMixture(Foam::dictionary const&, Foam::List<Foam::word> const&, Foam::HashPtrTable<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8>, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #6 Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::reactingMixture(Foam::dictionary const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #7 Foam::heThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #8 Foam::rhoReactionThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #9 Foam::autoPtr<Foam::rhoReactionThermo> Foam::basicThermo::New<Foam::rhoReactionThermo>(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #10 Foam::rhoReactionThermo::New(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" #11 Foam::rhoChemistryModel::rhoChemistryModel(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #12 Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::chemistryModel(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #13 Foam::ode<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::ode(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #14 Foam::rhoChemistryModel::addfvMeshConstructorToTable<Foam::ode<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #15 Foam::autoPtr<Foam::rhoChemistryModel> Foam::basicChemistryModel::New<Foam::rhoChemistryModel>(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #16 Foam::rhoChemistryModel::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so" #17 Foam::combustionModels::rhoChemistryCombustion::rhoChemistryCombustion(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so" #18 Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::laminar(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so" #19 Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::PaSR(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so" #20 Foam::combustionModels::rhoCombustionModel::adddictionaryConstructorToTable<Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion> >::New(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so" #21 Foam::combustionModels::rhoCombustionModel::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so" #22 in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/bin/rhoReactingFoam" #23 __libc_start_main in "/lib/i386linuxgnu/libc.so.6" #24 in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/bin/rhoReactingFoam" 

April 24, 2017, 03:11 

#2 
New Member
chang
Join Date: Dec 2016
Posts: 4
Rep Power: 3 
Hello!
I use the rhoReactingFoam to solver a chemistry problem, the thermoType is samilar to yours. I do not know why I can not run my case.I am confused. Can you help you? Or you can upload your case. Thanks you! 

April 24, 2017, 14:40 

#3 
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,246
Blog Entries: 40
Rep Power: 113 
Greetings to all!
@Tommy Floessner: Sorry that no one was able to get to you on this topic. My guess is that there was a division by zero, possibly because the initial temperature for one of the fluids was set to 0 Kelvin. I say this because of these two lines of the output: Code:
#1 Foam::sigFpe::sigHandler(int) in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::correctMassFractions() in "/home/dread/OpenFOAM/OpenFOAM2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so" The second line indicates that the error occurred when updating/calculating/correcting the mass fractions. @C_Cherry: Please provide the complete error message that you have gotten, the same way as Tommy did, because your error message might be similar, but could be due to an entirely different reason. Either way, my answer above might help. Best regards, Bruno
__________________


May 1, 2017, 22:26 

#4 
New Member
chang
Join Date: Dec 2016
Posts: 4
Rep Power: 3 
Thank you for your reply@wyldckat
In my case,the fluid flows in a pipe, the inlet temperature is 473K, the wall temperature is 873K, and the outlet is zeroGradient. The inlet pressure is zeroGradient, and the outlet pressure is 5MPa. The inlet velocity is 0.04244m/s, and the outlet velocity is zeroGradient. I use the rhoReactingFoam to solve a chemistry problem. I hope you can give me some help. I appreciate you very much! Here is the log file. Code:
Create time Create mesh for time = 0 Reading g Creating reaction model Selecting combustion model PaSR<rhoChemistryCombustion> Selecting chemistry type { chemistrySolver EulerImplicit; chemistryThermo rho; } Selecting thermodynamics package { type heRhoThermo; mixture reactingMixture; transport polynomial; thermo hPolynomial; equationOfState icoPolynomial; specie specie; energy sensibleEnthalpy; } Selecting chemistryReader foamChemistryReader chemistryModel: Number of species = 20 and reactions = 1 using integrated reaction rate Reading field U Reading/calculating face flux field phi Creating turbulence model. Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 0.33; sigmak 1; sigmaEps 1.3; Prt 1; } Creating field dpdt Creating field kinetic energy K No finite volume options present Courant Number mean: 4.77042e06 max: 0.000715562 PIMPLE: Operating solver in PISO mode Starting time loop Courant Number mean: 4.77042e06 max: 0.000715562 deltaT = 1.19904e05 Time = 1.19904e05 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCG: Solving for C10H22, Initial residual = 1, Final residual = 7.07715e11, No Iterations 2 DILUPBiCG: Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C2H4, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C2H6, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C3H6, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C3H8, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C4H8, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C4H10, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C5H10, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C5H12, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C6H12, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C6H14, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C7H14, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C7H16, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C8H16, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C8H18, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C9H18, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for C9H20, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 3.83483e08, No Iterations 1 min/max(T) = 473, 873 DICPCG: Solving for p, Initial residual = 1, Final residual = 56.6857, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00423238, global = 7.44309e05, cumulative = 7.44309e05 DICPCG: Solving for p, Initial residual = 0.179419, Final residual = 0.010874, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000257316, global = 4.964e06, cumulative = 7.93949e05 DILUPBiCG: Solving for epsilon, Initial residual = 0.00233739, Final residual = 8.03855e08, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 1.85515e09, No Iterations 2 ExecutionTime = 4.3 s ClockTime = 4 s Courant Number mean: 0.0108272 max: 0.0156534 deltaT = 1.43712e05 Time = 2.63616e05 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCG: Solving for C10H22, Initial residual = 0.227262, Final residual = 1.50724e09, No Iterations 2 DILUPBiCG: Solving for H2, Initial residual = 1.01627e07, Final residual = 1.01627e07, No Iterations 0 DILUPBiCG: Solving for CH4, Initial residual = 7.65904e07, Final residual = 7.65904e07, No Iterations 0 DILUPBiCG: Solving for C2H4, Initial residual = 2.39767e06, Final residual = 1.93074e10, No Iterations 1 DILUPBiCG: Solving for C2H6, Initial residual = 1.2649e06, Final residual = 1.01857e10, No Iterations 1 DILUPBiCG: Solving for C3H6, Initial residual = 3.23123e06, Final residual = 2.60197e10, No Iterations 1 DILUPBiCG: Solving for C3H8, Initial residual = 2.64994e06, Final residual = 2.13387e10, No Iterations 1 DILUPBiCG: Solving for C4H8, Initial residual = 3.67143e06, Final residual = 2.95644e10, No Iterations 1 DILUPBiCG: Solving for C4H10, Initial residual = 1.97929e06, Final residual = 1.59384e10, No Iterations 1 DILUPBiCG: Solving for C5H10, Initial residual = 4.00392e06, Final residual = 3.22418e10, No Iterations 1 DILUPBiCG: Solving for C5H12, Initial residual = 2.98688e06, Final residual = 2.4052e10, No Iterations 1 DILUPBiCG: Solving for C6H12, Initial residual = 5.47903e06, Final residual = 4.41202e10, No Iterations 1 DILUPBiCG: Solving for C6H14, Initial residual = 2.56059e06, Final residual = 2.06193e10, No Iterations 1 DILUPBiCG: Solving for C7H14, Initial residual = 5.53991e06, Final residual = 4.46104e10, No Iterations 1 DILUPBiCG: Solving for C7H16, Initial residual = 2.40867e06, Final residual = 1.9396e10, No Iterations 1 DILUPBiCG: Solving for C8H16, Initial residual = 5.69446e06, Final residual = 4.5855e10, No Iterations 1 DILUPBiCG: Solving for C8H18, Initial residual = 4.57642e07, Final residual = 4.57642e07, No Iterations 0 DILUPBiCG: Solving for C9H18, Initial residual = 2.23377e06, Final residual = 1.79876e10, No Iterations 1 DILUPBiCG: Solving for C9H20, Initial residual = 1.28459e07, Final residual = 1.28459e07, No Iterations 0 DILUPBiCG: Solving for h, Initial residual = 0.70559, Final residual = 1.41858e07, No Iterations 2 min/max(T) = 472.992, 873 DICPCG: Solving for p, Initial residual = 0.0220879, Final residual = 0.813472, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.0237036, global = 9.19515e05, cumulative = 1.25566e05 DICPCG: Solving for p, Initial residual = 0.786122, Final residual = 0.0045141, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000137576, global = 7.97722e06, cumulative = 2.05338e05 DILUPBiCG: Solving for epsilon, Initial residual = 0.00735973, Final residual = 4.25711e07, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.723139, Final residual = 9.64343e07, No Iterations 2 ExecutionTime = 8.52 s ClockTime = 9 s Courant Number mean: 0.0152289 max: 0.0186182 deltaT = 1.72256e05 Time = 4.35873e05 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCG: Solving for C10H22, Initial residual = 0.137911, Final residual = 4.1808e10, No Iterations 2 DILUPBiCG: Solving for H2, Initial residual = 2.2177e07, Final residual = 2.2177e07, No Iterations 0 DILUPBiCG: Solving for CH4, Initial residual = 1.67136e06, Final residual = 7.60635e11, No Iterations 1 DILUPBiCG: Solving for C2H4, Initial residual = 5.19659e06, Final residual = 2.36497e10, No Iterations 1 DILUPBiCG: Solving for C2H6, Initial residual = 2.7415e06, Final residual = 1.24766e10, No Iterations 1 DILUPBiCG: Solving for C3H6, Initial residual = 7.00318e06, Final residual = 3.18716e10, No Iterations 1 DILUPBiCG: Solving for C3H8, Initial residual = 5.74333e06, Final residual = 2.6138e10, No Iterations 1 DILUPBiCG: Solving for C4H8, Initial residual = 7.95722e06, Final residual = 3.62134e10, No Iterations 1 DILUPBiCG: Solving for C4H10, Initial residual = 4.28982e06, Final residual = 1.9523e10, No Iterations 1 DILUPBiCG: Solving for C5H10, Initial residual = 8.67782e06, Final residual = 3.94929e10, No Iterations 1 DILUPBiCG: Solving for C5H12, Initial residual = 6.4736e06, Final residual = 2.94614e10, No Iterations 1 DILUPBiCG: Solving for C6H12, Initial residual = 1.18748e05, Final residual = 5.40424e10, No Iterations 1 DILUPBiCG: Solving for C6H14, Initial residual = 5.54969e06, Final residual = 2.52567e10, No Iterations 1 DILUPBiCG: Solving for C7H14, Initial residual = 1.20067e05, Final residual = 5.46429e10, No Iterations 1 DILUPBiCG: Solving for C7H16, Initial residual = 5.22043e06, Final residual = 2.37583e10, No Iterations 1 DILUPBiCG: Solving for C8H16, Initial residual = 1.23417e05, Final residual = 5.61672e10, No Iterations 1 DILUPBiCG: Solving for C8H18, Initial residual = 9.98669e07, Final residual = 9.98669e07, No Iterations 0 DILUPBiCG: Solving for C9H18, Initial residual = 4.84136e06, Final residual = 2.20331e10, No Iterations 1 DILUPBiCG: Solving for C9H20, Initial residual = 2.80325e07, Final residual = 2.80325e07, No Iterations 0 DILUPBiCG: Solving for h, Initial residual = 0.348893, Final residual = 2.51228e07, No Iterations 2 min/max(T) = 472.999, 873 DICPCG: Solving for p, Initial residual = 0.0961288, Final residual = 0.520595, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00396067, global = 2.93884e05, cumulative = 4.99222e05 DICPCG: Solving for p, Initial residual = 0.356093, Final residual = 0.00296496, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 3.3259e05, global = 7.44155e07, cumulative = 5.06663e05 DILUPBiCG: Solving for epsilon, Initial residual = 0.00521535, Final residual = 4.15496e07, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.470671, Final residual = 1.24424e07, No Iterations 2 ExecutionTime = 13.51 s ClockTime = 14 s Courant Number mean: 0.014602 max: 0.0181641 deltaT = 2.06565e05 Time = 6.42437e05 ...... ...... ...... ...... ...... Courant Number mean: 0.0876726 max: 0.288901 deltaT = 0.002 Time = 0.246 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCG: Solving for C10H22, Initial residual = 0.00637298, Final residual = 5.46549e07, No Iterations 4 DILUPBiCG: Solving for H2, Initial residual = 0.0670608, Final residual = 1.84323e07, No Iterations 6 DILUPBiCG: Solving for CH4, Initial residual = 0.108365, Final residual = 2.97852e07, No Iterations 6 DILUPBiCG: Solving for C2H4, Initial residual = 0.11579, Final residual = 3.18261e07, No Iterations 6 DILUPBiCG: Solving for C2H6, Initial residual = 0.112548, Final residual = 3.09352e07, No Iterations 6 DILUPBiCG: Solving for C3H6, Initial residual = 0.116759, Final residual = 3.20924e07, No Iterations 6 DILUPBiCG: Solving for C3H8, Initial residual = 0.116146, Final residual = 3.19239e07, No Iterations 6 DILUPBiCG: Solving for C4H8, Initial residual = 0.117097, Final residual = 3.21853e07, No Iterations 6 DILUPBiCG: Solving for C4H10, Initial residual = 0.115008, Final residual = 3.16112e07, No Iterations 6 DILUPBiCG: Solving for C5H10, Initial residual = 0.117303, Final residual = 3.22421e07, No Iterations 6 DILUPBiCG: Solving for C5H12, Initial residual = 0.116529, Final residual = 3.20293e07, No Iterations 6 DILUPBiCG: Solving for C6H12, Initial residual = 0.117923, Final residual = 3.24124e07, No Iterations 6 DILUPBiCG: Solving for C6H14, Initial residual = 0.116027, Final residual = 3.18914e07, No Iterations 6 DILUPBiCG: Solving for C7H14, Initial residual = 0.117941, Final residual = 3.24175e07, No Iterations 6 DILUPBiCG: Solving for C7H16, Initial residual = 0.115807, Final residual = 3.18308e07, No Iterations 6 DILUPBiCG: Solving for C8H16, Initial residual = 0.117987, Final residual = 3.243e07, No Iterations 6 DILUPBiCG: Solving for C8H18, Initial residual = 0.101897, Final residual = 2.80074e07, No Iterations 6 DILUPBiCG: Solving for C9H18, Initial residual = 0.115517, Final residual = 3.17512e07, No Iterations 6 DILUPBiCG: Solving for C9H20, Initial residual = 0.0738386, Final residual = 2.02953e07, No Iterations 6 DILUPBiCG: Solving for h, Initial residual = 0.0075258, Final residual = 4.35917e08, No Iterations 4 min/max(T) = 452.174, 873 DICPCG: Solving for p, Initial residual = 0.21668, Final residual = 0.195115, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00845023, global = 0.00029595, cumulative = 0.000458058 DICPCG: Solving for p, Initial residual = 0.125125, Final residual = 0.0586924, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00519283, global = 3.89853e05, cumulative = 0.000419072 DILUPBiCG: Solving for epsilon, Initial residual = 0.0448897, Final residual = 4.74853e08, No Iterations 4 DILUPBiCG: Solving for k, Initial residual = 0.321046, Final residual = 8.4238e07, No Iterations 4 ExecutionTime = 958.48 s ClockTime = 958 s Courant Number mean: 0.152591 max: 0.715227 deltaT = 0.001 Time = 0.247 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCG: Solving for C10H22, Initial residual = 0.00430238, Final residual = 4.93826e07, No Iterations 3 DILUPBiCG: Solving for H2, Initial residual = 0.085043, Final residual = 2.71822e07, No Iterations 5 DILUPBiCG: Solving for CH4, Initial residual = 0.107737, Final residual = 3.44359e07, No Iterations 5 DILUPBiCG: Solving for C2H4, Initial residual = 0.110815, Final residual = 3.54199e07, No Iterations 5 DILUPBiCG: Solving for C2H6, Initial residual = 0.1095, Final residual = 3.49996e07, No Iterations 5 DILUPBiCG: Solving for C3H6, Initial residual = 0.1112, Final residual = 3.55429e07, No Iterations 5 DILUPBiCG: Solving for C3H8, Initial residual = 0.110957, Final residual = 3.54652e07, No Iterations 5 DILUPBiCG: Solving for C4H8, Initial residual = 0.111334, Final residual = 3.55855e07, No Iterations 5 DILUPBiCG: Solving for C4H10, Initial residual = 0.110502, Final residual = 3.53198e07, No Iterations 5 DILUPBiCG: Solving for C5H10, Initial residual = 0.111415, Final residual = 3.56116e07, No Iterations 5 DILUPBiCG: Solving for C5H12, Initial residual = 0.111109, Final residual = 3.55139e07, No Iterations 5 DILUPBiCG: Solving for C6H12, Initial residual = 0.111658, Final residual = 3.56892e07, No Iterations 5 DILUPBiCG: Solving for C6H14, Initial residual = 0.11091, Final residual = 3.54502e07, No Iterations 5 DILUPBiCG: Solving for C7H14, Initial residual = 0.111665, Final residual = 3.56915e07, No Iterations 5 DILUPBiCG: Solving for C7H16, Initial residual = 0.110822, Final residual = 3.54221e07, No Iterations 5 DILUPBiCG: Solving for C8H16, Initial residual = 0.111683, Final residual = 3.56971e07, No Iterations 5 DILUPBiCG: Solving for C8H18, Initial residual = 0.104853, Final residual = 3.35141e07, No Iterations 5 DILUPBiCG: Solving for C9H18, Initial residual = 0.110706, Final residual = 3.53851e07, No Iterations 5 DILUPBiCG: Solving for C9H20, Initial residual = 0.0895874, Final residual = 2.86347e07, No Iterations 5 DILUPBiCG: Solving for h, Initial residual = 0.00516729, Final residual = 9.07852e07, No Iterations 3 min/max(T) = 452.165, 873 DICPCG: Solving for p, Initial residual = 0.140282, Final residual = 2.44093, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.083028, global = 0.000464274, cumulative = 0.000883346 DICPCG: Solving for p, Initial residual = 0.468191, Final residual = 0.0195225, No Iterations 1001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00396294, global = 5.42544e05, cumulative = 0.000829091 DILUPBiCG: Solving for epsilon, Initial residual = 0.0475024, Final residual = 3.84935e07, No Iterations 3 DILUPBiCG: Solving for k, Initial residual = 0.389679, Final residual = 5.24729e07, No Iterations 4 ExecutionTime = 962.36 s ClockTime = 962 s Courant Number mean: 0.172735 max: 0.896032 deltaT = 0.000428571 Time = 0.247429 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/home/wangchang/OpenFOAM/OpenFOAM2.4.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #3 Foam::EulerImplicit<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::solve(Foam::Field<double>&, double&, double&, double&, double&) const at ??:? #4 double Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::solve<Foam::UniformField<double> >(Foam::UniformField<double> const&) at ??:? #5 Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::solve(double) at ??:? #6 Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:? #7 Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:? #8 ? at ??:? #9 __libc_start_main in "/lib/i386linuxgnu/libc.so.6" #10 ? at ??:? Last edited by wyldckat; May 8, 2017 at 19:03. Reason: Added [CODE][/CODE] markers 

April 22, 2018, 12:30 

#5 
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,246
Blog Entries: 40
Rep Power: 113 
Greetings C_Cherry,
My apologies, it's being nearly a year, and I never managed to find the time to look into the problem you were having. Hopefully you managed to solve the issue. If you haven't, the crash seems to have occurred when it was trying to calculate the equation for C10H22, due to not being able to solve one of the polynomial expressions (resulted in an invalid floating point operation, hence the sigFpe)... not sure which polynomial it was exactly, if it was receiving pressure and returning temperature or viceversa. My guess is that the initial pressure inside the domain was set to 0, instead of 1 or 5 MPa, which would somewhat explain the sudden crash due to not being able to calculate things properly, e.g. using a very high pressure as one of the input variables into the polynomial, could not render a valid solution. Best regards, Bruno
__________________


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Segmentation fault when running dieselFoam or dieselEngineFoam in parallel  francesco  OpenFOAM Bugs  4  May 2, 2017 21:59 
Segmentation fault in SU2 V5.0  ygd  SU2  2  March 1, 2017 04:38 
Segmentation fault when running in parallel  Pj.  OpenFOAM Running, Solving & CFD  3  April 8, 2015 08:12 
Segmentation Fault w/ compiled OF 2.2.0  motorBike example  sudo  OpenFOAM Running, Solving & CFD  3  April 2, 2013 17:27 
segmentation fault when installing OF2.1.1 on a cluster  Rebecca513  OpenFOAM Installation  9  July 31, 2012 15:06 