# [OpenFOAM v8] Why does Lun Savage radial function fail to converge?

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

 March 8, 2021, 01:23 [OpenFOAM v8] Why does Lun Savage radial function fail to converge? #1 New Member   Edward Ng Kian Guan Join Date: Feb 2021 Posts: 5 Rep Power: 5 I was testing usage of different radial distribution functions by modifying the tutorial case "fluidisedBed". Originally the tutorial case used radial distribution of Sinclair Jackson. The modified tutorial case uses radial distribution of Lun Savage, and had failed to converge at t = 1.41, yielding the error log shown at the bottom of this post. From my understanding, the only difference between Sinclair Jackson's function and Lun Savage function, is the power of the function. I am uncertain why did the simulation failed to converge. The modified case file is at,https://github.com/edwardKGN/fluidis...age_radialDist. Thank you for reading and helping! /start of error log Courant Number mean: 0.198841 max: 0.534878 deltaT = 0.002 Time = 1.408 PIMPLE: Iteration 1 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.105292, Final residual = 3.13208e-10, No Iterations 3 particles fraction, min, max = 0.274977 5.87739e-06 0.578727 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.673921, Final residual = 7.37451e-10, No Iterations 3 particles fraction, min, max = 0.274977 5.88012e-06 0.579236 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.740731, Final residual = 3.65256e-10, No Iterations 3 particles fraction, min, max = 0.274977 5.88293e-06 0.581966 Constructing momentum equations particles min/max T 600 - 600 water min/max T 300 - 600 GAMG: Solving for p_rgh, Initial residual = 0.0370966, Final residual = 0.000304703, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.00313407, Final residual = 4.84552e-09, No Iterations 16 PIMPLE: Iteration 2 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0344012, Final residual = 6.5296e-11, No Iterations 5 particles fraction, min, max = 0.274977 5.8776e-06 0.577119 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0762682, Final residual = 3.26397e-10, No Iterations 5 particles fraction, min, max = 0.274977 5.88029e-06 0.578326 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.141742, Final residual = 1.92945e-10, No Iterations 5 particles fraction, min, max = 0.274977 5.88307e-06 0.580262 Constructing momentum equations particles min/max T 600 - 600 water min/max T 300 - 600 GAMG: Solving for p_rgh, Initial residual = 0.0274588, Final residual = 0.000266116, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.00256743, Final residual = 8.92578e-09, No Iterations 15 PIMPLE: Iteration 3 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.0729752, Final residual = 1.91868e-10, No Iterations 3 particles fraction, min, max = 0.274977 5.87765e-06 0.580288 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.727065, Final residual = 2.14024e-11, No Iterations 3 particles fraction, min, max = 0.274977 5.8804e-06 0.584208 MULES: Solving for alpha.particles DICPBiCGStab: Solving for alpha.particles, Initial residual = 0.772823, Final residual = 4.08709e-10, No Iterations 3 particles fraction, min, max = 0.274977 5.88322e-06 0.586477 Constructing momentum equations particles min/max T 600 - 600 water min/max T 300 - 600 GAMG: Solving for p_rgh, Initial residual = 0.115123, Final residual = 0.00058304, No Iterations 4 GAMG: Solving for p_rgh, Initial residual = 0.00681858, Final residual = 5.59614e-09, No Iterations 17 DILUPBiCGStab: Solving for Theta.particles, Initial residual = 0.855074, Final residual = 5.41886e-07, No Iterations 1 ExecutionTime = 149.94 s Courant Number mean: 0.201101 max: 0.535999 deltaT = 0.002 Time = 1.41 PIMPLE: Iteration 1 MULES: Solving for alpha.particles #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 ? in "/lib/x86_64-linux-gnu/libm.so.6" #4 powf64 in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::pow(Foam::Field&, Foam::UList const&, double const&) at ??:? #6 Foam::tmp > Foam::pow(Foam::tmp > const&, Foam::dimensioned const&) at ??:? #7 Foam::kineticTheoryModels::radialModels::LunSavage ::g0prime(Foam::GeometricField const&, Foam::dimensioned const&, Foam::dimensioned const&) const at ??:? #8 Foam::RASModels::kineticTheoryModel::pPrime() const at ??:? #9 Foam::MovingPhaseModel >::pPrime() const at ??:? #10 Foam::MomentumTransferPhaseSystem::DByAfs(Foam::PtrList > const&, Foam::PtrList > const&) const at ??:? #11 Foam::phaseSystem::solve(Foam::PtrList > const&, Foam::PtrList > const&) at ??:? #12 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" #13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" /end of error log mabrouka likes this.

 July 15, 2022, 09:56 Lun #2 New Member   Abdo Join Date: Jan 2022 Location: Canada Posts: 4 Rep Power: 4 Hi there, Have you solved your problem, because I have same error