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

OpenFOAM solver cannot solve the equation for two of the species

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By dlahaye
  • 1 Post By Amirreza_pro
  • 1 Post By Amirreza_pro

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   January 24, 2022, 12:53
Default OpenFOAM solver cannot solve the equation for two of the species
  #1
Member
 
Amirreza Niazmehr
Join Date: Nov 2018
Posts: 40
Rep Power: 7
Amirreza_pro is on a distinguished road
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:rintStack(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

Last edited by Amirreza_pro; January 25, 2022 at 07:07.
Amirreza_pro is offline   Reply With Quote

 


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ANSYS Meshing] Help with element size sandri_92 ANSYS Meshing & Geometry 14 November 14, 2018 07:54
Species equation in OpenFOAM using VOF method SaufiSangi OpenFOAM Programming & Development 14 June 11, 2018 10:38
OpenFOAM v3.0+ ?? SBusch OpenFOAM 22 December 26, 2016 14:24
Suggestion for a new sub-forum at OpenFOAM's Forum wyldckat Site Help, Feedback & Discussions 20 October 28, 2014 09:04
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01


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