|
[Sponsors] |
August 30, 2022, 04:31 |
XiFoam error: Failing on multiple meshes
|
#1 |
New Member
Dilshan Casseer
Join Date: Jan 2018
Posts: 15
Rep Power: 8 |
Hello,
I'm trying to simulate the combustion of a simple case of H2 combustion with air in a square container of 12m x 3m x 3m. I set the ignition at the middle, and tried to simulate a case with o.75 equivalence ratio with inhomogenous method. The mesh is a entirely hexagonal mesh. LES with KEqn SGS model is used. I use pimple with 4 iteration steps. Solver is GAMG. I tried PCG as well, but got same result. I tried different meshes (40mm, 20mm, 30mm) but always I get an error and calculation stops after about 0.02 seconds. At this time I noticed that b goes t negative values, which is impossible. I have used limitedlinear01 as the div scheme for b, so it should essentially limit b within 0 and 1, but it doesn't seem to happen here. I always get the below error. can anyone please shed any light on what the error can be. #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:? #6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? #9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/XiFoam" #10 Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/XiFoam" #11 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/XiFoam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/XiFoam" I can't wrap my head around what is happening and has been stuck here for a couple of weeks, so any help is greatly appreciated. |
|
August 30, 2022, 21:34 |
|
#2 |
New Member
Dilshan Casseer
Join Date: Jan 2018
Posts: 15
Rep Power: 8 |
Anyone can give me some idea on this please?
|
|
March 2, 2023, 00:39 |
|
#3 |
New Member
Sarah
Join Date: Jan 2023
Posts: 12
Rep Power: 3 |
Hi
Did you find any insights on this issue? I am also getting negative values of b, specifically in case of 3D simulations, that doesn't seem to happen with 2D cases. |
|
Tags |
combustion, hydrogen, xifoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Preserve meshes in Fluent parametric run | hwerner2014 | FLUENT | 2 | February 28, 2022 17:29 |
how to set periodic boundary conditions | Ganesh | FLUENT | 15 | November 18, 2020 07:09 |
[mesh manipulation] Importing Multiple Meshes | thomasnwalshiii | OpenFOAM Meshing & Mesh Conversion | 18 | December 19, 2015 19:57 |
OpenFOAM static build on Cray XT5 | asaijo | OpenFOAM Installation | 9 | April 6, 2011 13:21 |
Large 3D tetrahedral meshes | Aldo Bonfiglioli | Main CFD Forum | 4 | August 27, 1999 04:33 |