[pupo]

I'm running a simulation on OPENFOAM 3.0 with a somewhat complex geometry and fvDOM is crashing on the first iteration. I've followed the error stack but i can't figure out what is wrong.

Mesh:



Files:

https://github.com/pupo162/fvDOMerror

error stack:

]
Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 0.99999998, Final residual = 2.4577608e-05, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.99999998, Final residual = 3.003951e-05, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 8.0562366e-06, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 9.487023e-08, No Iterations 2
Radiation solver iter: 0
#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::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::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 ??:?
#6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#7  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#9  Foam::radiation::radiativeIntensityRay::correct() at ??:?
#10  Foam::radiation::fvDOM::calculate() at ??:?
#11  Foam::radiation::radiationModel::correct() at ??:?
#12  ? at ??:?
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? at ??:?

I've run a a similar project with the exact same parameters for a simpler mesh with no problem, so I'm guessing the problem is in the Mesh. However, the simulation runs and converges fine with radition model: none.

I've tried to improve the mesh quality, reduce the point the number of elements, change the number of direction fvDOM uses (nTheta and nPhi) etc..., the only thing that has shown some effect is seting the "constantAbsorptionCoeficients" in the "RadiationProperties" File to values different 0, which results in the simulation crashing in the second iteration.

If i anyone could have a hint on what the problem could be, I would appreciate any help.

Best Regards,
Pedro

[pupo]

For future reference, It seems the fvDOM model can't handle corners where multiple cells are generated. In the imaga above these are the problematic cells:



An extreme example of what CAN'T happen is given in the image below:



Redesigning the mesh to avoid these cases seems to be the only solution.


Cheers