[edwardKGN]

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<double>&, Foam::UList<double> const&, double const&) at ??:?
#6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at ??:?
#7 Foam::kineticTheoryModels::radialModels::LunSavage ::g0prime(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&, Foam::dimensioned<double> const&) const at ??:?
#8 Foam::RASModels::kineticTheoryModel::pPrime() const at ??:?
#9 Foam::MovingPhaseModel<Foam::ThermoPhaseModel<Foam ::phaseModel, Foam::rhoThermo> >::pPrime() const at ??:?
#10 Foam::MomentumTransferPhaseSystem<Foam::phaseSyste m>::DByAfs(Foam::PtrList<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) const at ??:?
#11 Foam::phaseSystem::solve(Foam::PtrList<Foam::Geome tricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::PtrList<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > 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]

Hi there,

Have you solved your problem, because I have same error