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

fvDOM model crashes on first iteration

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 19, 2016, 08:25
Unhappy fvDOM model crashes on first iteration
  #1
Member
 
Pedro
Join Date: Nov 2014
Posts: 50
Rep Power: 10
pupo is on a distinguished road
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:

Code:
]
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

Last edited by pupo; May 19, 2016 at 10:08.
pupo is offline   Reply With Quote

Old   May 29, 2016, 09:10
Default
  #2
Member
 
Pedro
Join Date: Nov 2014
Posts: 50
Rep Power: 10
pupo is on a distinguished road
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
pupo is offline   Reply With Quote

Reply

Tags
fvdom, radiation

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Coupling RPI wall boiling model with population balance model in Fluent softice2006 Fluent Multiphase 0 February 1, 2015 16:55
Is it possible to model natural convection in a 2D horizontal model in fluent caitoc FLUENT 1 May 5, 2014 14:32
LRRturbulence Model crashes with symmetryPlanboundaries gschaider OpenFOAM Bugs 2 October 10, 2007 07:53
Advanced Turbulence Modeling in Fluent, Realizable k-epsilon Model Jonas Larsson FLUENT 5 March 13, 2000 04:27


All times are GMT -4. The time now is 20:28.