[mwaqas]

Hello Everyone

I am using OF2.1.1. to simulate fluidized bed. My geometry is 2D length 100cm width 28cm and initial bed height is 40cm. Initial alphs is 0.6 and restitution coefficient 0.8. I am using algebric KTFG (equilibrium = on), k-apsilon turbulance model with no particle slip boundary condition. Time step is 1e-04 and cell is 0.005*0.005 I am using following settings.
Drag model GidaspowErgunWenYu
viscosityModel Gidaspow; conductivityModel Gidaspow; granularPressureModel Lun; frictionalStressModel JohnsonJackson; radialModel LunSavage;//Gidaspow;

Initially I carried out simulation for 280 micron with 5Umf (0.37m/s) and I got reasonable result.
But when I increase the particle diameter to 728 micron with 5Umf (1.8 m/s). My simulation crashes after 0.5 sec.
I tried many things but could not solve my problem e.g. switched of turbulance changed radial distribution function, coarser and finer mesh, Symala drag model, 3Umf (1.134 m/s),
reduced time step i.e. 1e-5. but could not resolve my problem. The error message is as follows

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9 Foam::fvMatrix<double>::solve() in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
#10
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam

Can you please guide me where I am making mistake and how can i solve this problem. Thank you very much for time.