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/)
-   -   OpenFOAM solver cannot solve the equation for two of the species (https://www.cfd-online.com/Forums/openfoam-solving/240800-openfoam-solver-cannot-solve-equation-two-species.html)

Amirreza_pro January 24, 2022 12:53
OpenFOAM solver cannot solve the equation for two of the species
 
Hi guys.
I'm working on one of the OpenFOAM's tuturials which its name is SandiaD_LTS (reactingfoam (RAS)). This Model has 36 species but my model has 41 species. I converted my mechanism with ChemkinToFoam to openFoam's format and changed the files in ''constant'' folder accordingly.
The problem is that openFoam's solver can solve the equations for 39 species without any problem, but for two remaining species (A2XXC10H7 and A2XC10H8) it crashes and I get a error which shows max courant number is about 3000 !!!
Does anyone have any idea, how I can resolve this problem?
Here is the error:

Quote:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.2308810194, Final residual = 0.01553386741, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.007643954332, Final residual = 0.0003203595795, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.003917487069, Final residual = 7.599511572e-05, No Iterations 1
DILUPBiCGStab: Solving for SXCH2, Initial residual = 0.09864586648, Final residual = 0.004468863072, No Iterations 2
DILUPBiCGStab: Solving for TXCH2, Initial residual = 0.08400964579, Final residual = 0.004087688696, No Iterations 2
DILUPBiCGStab: Solving for O, Initial residual = 0.002594328112, Final residual = 0.0002288017869, No Iterations 1
DILUPBiCGStab: Solving for H2, Initial residual = 0.0008294417205, Final residual = 6.73505605e-05, No Iterations 1
DILUPBiCGStab: Solving for H, Initial residual = 0.01260778683, Final residual = 0.0006234817519, No Iterations 1
DILUPBiCGStab: Solving for OH, Initial residual = 0.004515186986, Final residual = 0.0004448165445, No Iterations 1
DILUPBiCGStab: Solving for H2O, Initial residual = 0.0003644372557, Final residual = 2.119358504e-05, No Iterations 1
DILUPBiCGStab: Solving for O2, Initial residual = 0.0003114845114, Final residual = 1.876150436e-05, No Iterations 1
DILUPBiCGStab: Solving for HO2, Initial residual = 0.02782765682, Final residual = 0.0005188633784, No Iterations 1
DILUPBiCGStab: Solving for CH, Initial residual = 0.02903526762, Final residual = 0.001284509125, No Iterations 1
DILUPBiCGStab: Solving for CO, Initial residual = 0.0007080496295, Final residual = 6.639013559e-05, No Iterations 1
DILUPBiCGStab: Solving for HCO, Initial residual = 0.1023357629, Final residual = 0.00763330224, No Iterations 1
DILUPBiCGStab: Solving for CH2O, Initial residual = 0.004106367887, Final residual = 0.0001626599125, No Iterations 1
DILUPBiCGStab: Solving for CH3, Initial residual = 0.02556436683, Final residual = 0.001587398739, No Iterations 1
DILUPBiCGStab: Solving for CO2, Initial residual = 0.0003262353295, Final residual = 1.652384498e-05, No Iterations 1
DILUPBiCGStab: Solving for CH4, Initial residual = 0.0003048411883, Final residual = 2.00545891e-05, No Iterations 1
DILUPBiCGStab: Solving for C2H3, Initial residual = 0.07879345427, Final residual = 0.007509705106, No Iterations 1
DILUPBiCGStab: Solving for C2H4, Initial residual = 0.005030547006, Final residual = 0.0002293907421, No Iterations 1
DILUPBiCGStab: Solving for C2H5, Initial residual = 0.03327498135, Final residual = 0.002326852783, No Iterations 1
DILUPBiCGStab: Solving for C2H, Initial residual = 0.02296921496, Final residual = 0.0008335992445, No Iterations 1
DILUPBiCGStab: Solving for HCCO, Initial residual = 0.06650894142, Final residual = 0.006033413946, No Iterations 1
DILUPBiCGStab: Solving for C2H2, Initial residual = 0.002642369275, Final residual = 0.0002598119158, No Iterations 1
DILUPBiCGStab: Solving for C3H3, Initial residual = 0.004317912858, Final residual = 0.0003373478663, No Iterations 1
DILUPBiCGStab: Solving for AXC3H5, Initial residual = 0.02316649507, Final residual = 0.001237440204, No Iterations 1
DILUPBiCGStab: Solving for C2H6, Initial residual = 0.008957664377, Final residual = 0.000620502687, No Iterations 1
DILUPBiCGStab: Solving for PXC3H4, Initial residual = 0.003119606592, Final residual = 0.0001325969569, No Iterations 1
DILUPBiCGStab: Solving for AXC3H4, Initial residual = 0.002883747604, Final residual = 0.0001075566616, No Iterations 1
DILUPBiCGStab: Solving for A1XC6H6, Initial residual = 0.002121768592, Final residual = 0.0001734824224, No Iterations 1
DILUPBiCGStab: Solving for A1XXC6H5, Initial residual = 0.05176928534, Final residual = 0.003402699708, No Iterations 1
DILUPBiCGStab: Solving for C5H5, Initial residual = 0.009168001324, Final residual = 0.0003734758903, No Iterations 1
DILUPBiCGStab: Solving for C3H6, Initial residual = 0.01125870698, Final residual = 0.0004783456533, No Iterations 1
DILUPBiCGStab: Solving for C5H6, Initial residual = 0.005735117599, Final residual = 0.0001768512892, No Iterations 1
DILUPBiCGStab: Solving for A2XC10H8, Initial residual = 0.1584552579, Final residual = 0.009169496086, No Iterations 1
DILUPBiCGStab: Solving for A1C2H2, Initial residual = 0.02461866143, Final residual = 0.0009880261956, No Iterations 2
DILUPBiCGStab: Solving for A1CH2, Initial residual = 0.008684250407, Final residual = 0.0004588292534, No Iterations 1
DILUPBiCGStab: Solving for A1CHO, Initial residual = 0.001207425656, Final residual = 8.051408779e-05, No Iterations 1
DILUPBiCGStab: Solving for A1CH3, Initial residual = 0.003701795297, Final residual = 0.000127854401, No Iterations 1
DILUPBiCGStab: Solving for A1C2HY, Initial residual = 0.001429971375, Final residual = 8.400895735e-05, No Iterations 1
DILUPBiCGStab: Solving for A2XXC10H7, Initial residual = 0.1390251998, Final residual = 0.002329714753, No Iterations 1
DILUPBiCGStab: Solving for A1C2H, Initial residual = 0.03225068257, Final residual = 0.000469650234, No Iterations 1
DICPCG: Solving for G, Initial residual = 0.00348341438, Final residual = 9.732898286e-06, No Iterations 89
DILUPBiCGStab: Solving for h, Initial residual = 0.0005187303513, Final residual = 3.614421769e-05, No Iterations 1
min/max(T) = 290.9900761, 2261.828346
DICPCG: Solving for p, Initial residual = 0.2354451128, Final residual = 0.002122483175, No Iterations 11
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.33393049e-05, global = -9.914027186e-06, cumulative = -3.522183011e-05
DICPCG: Solving for p, Initial residual = 0.1085366126, Final residual = 9.426420545e-07, No Iterations 40
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.1875993e-08, global = -1.78172965e-08, cumulative = -3.52396474e-05
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.03315216171, Final residual = 0.002909196169, No Iterations 1
bounding epsilon, min: -4636388078 max: 8.821650169e+10 average: 166034962.2
DILUPBiCGStab: Solving for k, Initial residual = 0.01703632292, Final residual = 0.00144347403, No Iterations 1
bounding k, min: -12836.36788 max: 679154.7754 average: 1697.827427
ExecutionTime = 179.66 s ClockTime = 179 s

Courant Number mean: 0.03305969975 max: 2961.561288
Time = 0.0487

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.2084757932, Final residual = 0.005058679761, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.02283007525, Final residual = 0.001650918653, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.008530772281, Final residual = 0.0001804308987, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::combustionModels::EDC::correct() at ??:?
#4 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception

dlahaye January 25, 2022 04:07
You can reduce the Courant number in system/controlDict by changing the value of maxCo.

Amirreza_pro January 25, 2022 05:10
Quote:
Originally Posted by dlahaye (Post 820946)
You can reduce the Courant number in system/controlDict by changing the value of maxCo.
Thanks for your reply. I did it but nothing changed.
I would be appreciated if someone could help me by this case.
I stuck from yesterday actually :(

Amirreza_pro January 25, 2022 06:02
Quote:
Originally Posted by dlahaye (Post 820946)
You can reduce the Courant number in system/controlDict by changing the value of maxCo.
Sorry, I made a mistake, I redueced max Co, but I did not activate adjustable Timestep, therefore OpenFoam could not recognize my changes. I activated it and ran the simulation. I'm waiting for the results, I hope it works :)

Amirreza_pro January 25, 2022 07:13
It's a little freaky.
When I choose max Coraunt a number less than 1, It runs succesfully but , I don't have any output of my two remaining species.
When I choose max Courant 1 or more than 1, again I get the same error and the solver does not solve the equations for two remaing species. here is the error for max Courant: 1

Quote:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.1923767326, Final residual = 0.002391625424, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.00391798032, Final residual = 4.889783009e-06, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.01986269349, Final residual = 2.716365484e-05, No Iterations 1
DILUPBiCGStab: Solving for SXCH2, Initial residual = 0.001799432879, Final residual = 0.0001198556525, No Iterations 1
DILUPBiCGStab: Solving for TXCH2, Initial residual = 0.001549822602, Final residual = 0.000101872846, No Iterations 1
DILUPBiCGStab: Solving for O, Initial residual = 0.0002592933617, Final residual = 8.206062757e-06, No Iterations 1
DILUPBiCGStab: Solving for H2, Initial residual = 0.0001317551326, Final residual = 3.643360878e-06, No Iterations 1
DILUPBiCGStab: Solving for H, Initial residual = 0.000303695481, Final residual = 1.479755741e-05, No Iterations 1
DILUPBiCGStab: Solving for OH, Initial residual = 0.0001562788067, Final residual = 5.696827152e-06, No Iterations 1
DILUPBiCGStab: Solving for H2O, Initial residual = 9.925404312e-05, Final residual = 6.535633111e-07, No Iterations 1
DILUPBiCGStab: Solving for O2, Initial residual = 9.595673686e-05, Final residual = 6.321163472e-07, No Iterations 1
DILUPBiCGStab: Solving for HO2, Initial residual = 0.0004574734167, Final residual = 1.332117002e-06, No Iterations 1
DILUPBiCGStab: Solving for CH, Initial residual = 0.003279340848, Final residual = 0.0002206604818, No Iterations 1
DILUPBiCGStab: Solving for CO, Initial residual = 0.0001691771508, Final residual = 3.537159218e-06, No Iterations 1
DILUPBiCGStab: Solving for HCO, Initial residual = 0.002210840106, Final residual = 0.0001045366478, No Iterations 1
DILUPBiCGStab: Solving for CH2O, Initial residual = 0.0003776412299, Final residual = 3.583062651e-06, No Iterations 1
DILUPBiCGStab: Solving for CH3, Initial residual = 0.0003213562774, Final residual = 1.59729028e-05, No Iterations 1
DILUPBiCGStab: Solving for CO2, Initial residual = 0.0001060724704, Final residual = 8.319711586e-07, No Iterations 1
DILUPBiCGStab: Solving for CH4, Initial residual = 9.184257088e-05, Final residual = 5.052811027e-07, No Iterations 1
DILUPBiCGStab: Solving for C2H3, Initial residual = 0.0007250013454, Final residual = 5.122902633e-05, No Iterations 1
DILUPBiCGStab: Solving for C2H4, Initial residual = 8.753288733e-05, Final residual = 2.566754919e-06, No Iterations 1
DILUPBiCGStab: Solving for C2H5, Initial residual = 0.0002128167048, Final residual = 6.544289346e-06, No Iterations 1
DILUPBiCGStab: Solving for C2H, Initial residual = 0.001307048086, Final residual = 4.953380643e-05, No Iterations 1
DILUPBiCGStab: Solving for HCCO, Initial residual = 0.0007889935755, Final residual = 3.047812797e-05, No Iterations 1
DILUPBiCGStab: Solving for C2H2, Initial residual = 0.0001748495501, Final residual = 6.326692874e-06, No Iterations 1
DILUPBiCGStab: Solving for C3H3, Initial residual = 0.0003308128211, Final residual = 1.422627534e-05, No Iterations 1
DILUPBiCGStab: Solving for AXC3H5, Initial residual = 0.0001534220417, Final residual = 5.114761825e-06, No Iterations 1
DILUPBiCGStab: Solving for C2H6, Initial residual = 0.0001505994743, Final residual = 1.784370591e-06, No Iterations 1
DILUPBiCGStab: Solving for PXC3H4, Initial residual = 0.0001273287817, Final residual = 4.208864481e-06, No Iterations 1
DILUPBiCGStab: Solving for AXC3H4, Initial residual = 0.0001282577155, Final residual = 4.126418403e-06, No Iterations 1
DILUPBiCGStab: Solving for A1XC6H6, Initial residual = 0.0002546660534, Final residual = 1.126772027e-05, No Iterations 1
DILUPBiCGStab: Solving for A1XXC6H5, Initial residual = 0.001089981383, Final residual = 4.90641181e-05, No Iterations 1
DILUPBiCGStab: Solving for C5H5, Initial residual = 0.0001420893127, Final residual = 4.549401379e-06, No Iterations 1
DILUPBiCGStab: Solving for C3H6, Initial residual = 7.713533691e-05, Final residual = 2.008487741e-06, No Iterations 1
DILUPBiCGStab: Solving for C5H6, Initial residual = 0.0001759896299, Final residual = 8.6905927e-06, No Iterations 1
DILUPBiCGStab: Solving for A2XC10H8, Initial residual = 0.0002798874123, Final residual = 1.883528242e-05, No Iterations 1
DILUPBiCGStab: Solving for A1C2H2, Initial residual = 0.00129549968, Final residual = 8.608363412e-05, No Iterations 1
DILUPBiCGStab: Solving for A1CH2, Initial residual = 0.0006438713843, Final residual = 2.803087791e-05, No Iterations 1
DILUPBiCGStab: Solving for A1CHO, Initial residual = 0.0008167387307, Final residual = 1.785156524e-05, No Iterations 1
DILUPBiCGStab: Solving for A1CH3, Initial residual = 0.0004428271792, Final residual = 1.423557502e-05, No Iterations 1
DILUPBiCGStab: Solving for A1C2HY, Initial residual = 0.01375418298, Final residual = 0.0002514268187, No Iterations 1
DILUPBiCGStab: Solving for A2XXC10H7, Initial residual = 0.0002159570297, Final residual = 1.119433101e-05, No Iterations 1
DILUPBiCGStab: Solving for A1C2H, Initial residual = 0.008374806879, Final residual = 0.000595453154, No Iterations 1
DICPCG: Solving for G, Initial residual = 0.001704266578, Final residual = 8.63358684e-06, No Iterations 86
DILUPBiCGStab: Solving for h, Initial residual = 0.0001090955023, Final residual = 1.062028985e-06, No Iterations 1
min/max(T) = 289.220187, 2045.724495
DICPCG: Solving for p, Initial residual = 0.04938138081, Final residual = 0.0003789587478, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.321793177e-06, global = 2.612965337e-06, cumulative = 9.82923144e-06
DICPCG: Solving for p, Initial residual = 0.01617178602, Final residual = 8.426708442e-07, No Iterations 28
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.589315925e-08, global = 3.911119184e-09, cumulative = 9.833142559e-06
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.03834642153, Final residual = 0.00223306707, No Iterations 1
bounding epsilon, min: -4.193695319e+10 max: 3.972588841e+10 average: 103253794.7
DILUPBiCGStab: Solving for k, Initial residual = 0.0252667087, Final residual = 0.0002103159889, No Iterations 1
bounding k, min: -338362.8481 max: 346789.2192 average: 925.2429515
ExecutionTime = 189.97 s ClockTime = 190 s

Courant Number mean: 0.01228592804 max: 6.849323873
deltaT = 8.281062879e-06
Time = 0.0491282

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.05794917554, Final residual = 0.0009617219072, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.0005922981551, Final residual = 4.494267625e-07, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.003303595724, Final residual = 6.420447976e-06, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::combustionModels::EDC::correct() at ??:?
#4 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception
I did many tricks, but none of them works :(

Amirreza_pro January 26, 2022 03:39
Update:
I fixed the problem. What I have done:
1- reducing max courant number
2- refining the mesh
3- improving initial conditions
4- changing Temperature limits


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